Stochastic method of optimizing the basis for variational calculations of many-particle systems

V. I. Kukulin

Paolo Recchia, Debabrota Basu, M. Gattobigio +2 more

Quantum few-body systems are deceptively simple. Indeed, with the notable exception of a few special cases, their associated Schrodinger equation cannot be solved analytically for more than two particles. One has to resort to approximation methods to tackle quantum few-body problems. In particular, variational methods have been proposed to ease numerical calculations and obtain precise solutions. One such method is the Stochastic Variational Method, which employs a stochastic search to determine the number and parameters of correlated Gaussian basis functions used to construct an ansatz of the wave function. Stochastic methods, however, face numerical and optimization challenges as the number of particles increases. We introduce a family of gradient variational methods that replace stochastic search with gradient optimization. We comparatively and empirically evaluate the performance of the baseline Stochastic Variational Method, several instances of the gradient variational method family, and some hybrid methods for selected few-body problems. We show that gradient and hybrid methods can be more efficient and effective than the Stochastic Variational Method. We discuss the role of singularities, oscillations, and gradient optimization strategies in the performance of the respective methods.

Vladimir M. Krasnopolsky, V. I. Kukulin

Paolo Recchia, Debabrota Basu, M. Gattobigio +2 more

Abstract Quantum few-body systems are deceptively simple. Indeed, with the notable exception of a few special cases, their associated Schrödinger equation cannot be solved analytically for more than two particles. One has to resort to approximation methods to tackle quantum few-body problems. In particular, variational methods have been proposed to ease numerical calculations and obtain precise solutions. One such method is the Stochastic Variational Method, which employs a stochastic search to determine the number and parameters of correlated Gaussian basis functions used to construct an ansatz of the wave function. Stochastic methods, however, face numerical and optimization challenges as the number of particles increases. We introduce a family of gradient variational methods that replace stochastic search with gradient optimization. We comparatively and empirically evaluate the performance of the baseline Stochastic Variational Method, several instances of the gradient variational method family, and some hybrid methods for selected few-body problems. We show that gradient and hybrid methods can be more efficient and effective than the Stochastic Variational Method. We discuss the role of singularities, oscillations, and gradient optimization strategies in the performance of the respective methods.

Jay D. Mancini, Vassilios Fessatidis, Samuel P. Bowen +1 more

Kálmán Varga, Y. Suzuki

Precise variational solutions are given for problems involving diverse fermionic and bosonic (N=2--7)-body systems. The trial wave functions are chosen to be combinations of correlated Gaussians, which are constructed from products of the single-particle Gaussian wave packets through an integral transformation, thereby facilitating fully analytical calculations of the matrix elements. The nonlinear parameters of the trial function are chosen by a stochastic technique. The method has proved very efficient, virtually exact, and it seems feasible for any few-body bound-state problems emerging in nuclear or atomic physics.