Simple one-body term in variational quantum Monte Carlo calculations for solids and atoms

Paola Gori‐Giorgi, ALESSANDRO FEDERICO, Giovanni B. Bachelet

R. Gaudoin, Maziar Nekovee, W. M. C. Foulkes +2 more

Quantum Monte Carlo simulations of interacting electrons in solids often use Slater-Jastrow trial wave functions with Jastrow factors containing one- and two-body terms. In uniform systems the long-range behavior of the two-body term may be deduced from the random-phase approximation (RPA) of Bohm and Pines. Here we generalize the RPA to nonuniform systems. This gives the long-range behavior of the inhomogeneous two-body correlation term and provides an accurate analytic expression for the one-body term. It also explains why Slater-Jastrow trial wave functions incorporating determinants of Hartree-Fock or density-functional orbitals are close to optimal even in the presence of an RPA Jastrow factor. After adjusting the inhomogeneous RPA Jastrow factor to incorporate the known short-range behavior, we test it using variational Monte Carlo calculations. We find that the most important aspect of the two-body term is the short-range behavior due to electron-electron scattering, although the long-range behavior described by the RPA should become more important at high densities.

Dyutiman Das, Richard M. Martin, Malvin H. Kalos

The quantum Monte Carlo (QMC) technique is an extremely powerful method to treat many-body systems. Usually the quantum Monte Carlo method has been applied in cases where the interaction potential has a simple analytic form, like the 1/r Coulomb potential. However, in a complicated environment as in a semiconductor heterostructure, the evaluation of the interaction itself becomes a nontrivial problem. Obtaining the potential from any grid-based finite-difference method for every walker and every step is infeasible. We demonstrate an alternative approach of solving the Poisson equation by a classical Monte Carlo calculation within the overall quantum Monte Carlo scheme. We have developed a modified "walk on spheres" algorithm using Green's function techniques, which can efficiently account for the interaction energy of walker configurations, typical of quantum Monte Carlo algorithms. This stochastically obtained potential can be easily incorporated with variational, diffusion, and other Monte Carlo techniques. We demonstrate the validity of this method by studying a simple problem, the polarization of a helium atom in the electric field of an infinite capacitor.

Håvard Sandsdalen

Markus Holzmann

Clio Johnson

Variational and diffusion Monte Carlo (VMC and DMC) are widely used techniques in the Monte Carlo study of many-body quantum systems. Through efficient use of ran- dom sampling to solve many-dimensional integrals, they can be used to reliably solve the Schrödinger equation. Within the wider ecosystems of electronic structure calculation, quantum chemistry, and quantum simulation, quantum Monte Carlo (QMC) methods are used in tandem with density functional theory (DFT) and other computational methods to deliver insight into a vast array of condensed matter systems. Of particular interest to this work are systems in the intermediate correlation regime, where methods such as simple band theory and the Hubbard model break down. This thesis proposes and tests an explicit three-body backflow term to improve the accuracy of VMC and DMC calculations by incorporating higher-order correlation in order to reduce fixed-node er- ror in DMC, finding statistically significant improvements in ground state energies and improved variance in a water molecule. Then, QMC is used to calculate accurate charge densities for weak and intermediately correlated materials, to test commonly used dens- ity functional approximations in a controlled environment, providing insight into where these approximations perform well, and where they fail. Finally, DMC is used to generate ground state energy data for a dipolar Fermi gas of dressed Rydberg atoms interacting in the intermediate correlation regime. This data is parametrised in terms of interaction strength, for use in generating exchange-correlation functionals, to enable the DFT study of inhomogeneous systems of dressed Rydberg atoms.