Auxiliary-field Quantum Monte Carlo Method Research Articles

A study of H+H2 and several H-bonded molecules by phaseless auxiliary-field quantum Monte Carlo with plane wave and Gaussian basis sets

The authors present phaseless auxiliary-field (AF) quantum Monte Carlo (QMC) calculations of the ground states of some hydrogen-bonded systems. These systems were selected to test and benchmark different aspects of the new phaseless AF QMC method. They include the transition state of H+H(2) near the equilibrium geometry and in the van der Walls limit, as well as the H(2)O, OH, and H(2)O(2) molecules. Most of these systems present significant challenges for traditional independent-particle electronic structure approaches, and many also have exact results available. The phaseless AF QMC method is used either with a plane wave basis with pseudopotentials or with all-electron Gaussian basis sets. For some systems, calculations are done with both to compare and characterize the performance of AF QMC under different basis sets and different Hubbard-Stratonovich decompositions. Excellent results are obtained using as input single Slater determinant wave functions taken from independent-particle calculations. Comparisons of the Gaussian based AF QMC results with exact full configuration interaction show that the errors from controlling the phase problem with the phaseless approximation are small. At the large basis-size limit, the AF QMC results using both types of basis sets are in good agreement with each other and with experimental values.

Auxiliary-field quantum Monte Carlo study of first- and second-row post-d elements

A series of calculations for the first- and second-row post-d elements (Ga-Br and In-I) are presented using the phaseless auxiliary-field quantum Monte Carlo (AF QMC) method. This method is formulated in a Hilbert space defined by any chosen one-particle basis and maps the many-body problem into a linear combination of independent-particle solutions with external auxiliary fields. The phase/sign problem is handled approximately by the phaseless formalism using a trial wave function, which in our calculations was chosen to be the Hartree-Fock solution. We used the consistent correlated basis sets of Peterson et al. [J. Chem. Phys. 119, 11099 (2003); 119, 11113 (2003)], which employ a small-core relativistic pseudopotential. The AF QMC results are compared with experiment and with those from density functional (generalized gradient approximation and B3LYP) and CCSD(T) calculations. The AF QMC total energies agree with CCSD(T) to within a few millihartrees across the systems and over several basis sets. The calculated atomic electron affinities, ionization energies, and spectroscopic properties of dimers are, at large basis sets, in excellent agreement with experiment.

Correlation effects in the ground state of trapped atomic Bose gases

We study the effects of many-body correlations in trapped ultracold atomic Bose gases. We calculate the ground state of the gas using a ground-state auxiliary-field quantum Monte Carlo (QMC) method [Phys. Rev. E 70, 056702 (2004)]. We examine the properties of the gas, such as the energetics, condensate fraction, real-space density, and momentum distribution, as a function of the number of particles and the scattering length. We find that the mean-field Gross-Pitaevskii (GP) approach gives qualitatively incorrect result of the kinetic energy as a function of the scattering length. We present detailed QMC data for the various quantities, and discuss the behavior of GP, modified GP, and the Bogoliubov method under a local density approximation.

Adaptive sampling approach to the negative-sign problem in the auxiliary-field quantum Monte Carlo method

We propose a sampling method to calculate the ground state of interacting quantum systems. This approximate method, which we call the adaptive sampling quantum Monte Carlo (ASQMC) method, utilizes information from the short projection time density matrix derived from the Monte Carlo steps. The ASQMC method circumvents the negative-sign problem by introducing a conditional probability distribution function and guiding walkers to avoid a set of canceling fields. With the ASQMC method, the negative-sign ratio is greatly reduced. The ratio becomes zero in the limit $\ensuremath{\Delta}\ensuremath{\tau}$ goes to zero without prior knowledge of the exact ground state wave function. Comparisons with numerical results obtained by using other methods are made, and we observe that the ASQMC method gives accurate results with less than $0.5%$ error for the energy over wide regions of physical parameter values.

A discrete Hubbard-Stratonovich decomposition for general, fermionic two-body interactions

A scheme is presented to decompose the exponential of a two-body operator in a discrete sum over exponentials of one-body operators. This discrete decomposition can be used instead of the Hubbard-Stratonovich transformation in auxiliary-field quantum Monte-Carlo methods. As an illustration, the decomposition is applied to the Hubbard model, where it is equivalent to the discrete Hubbard-Stratonovich transformation introduced by Hirsch, and to the nuclear pairing Hamiltonian.

AUXILIARY FIELDS AND THE SIGN PROBLEM

The auxiliary-field quantum Monte Carlo method is reviewed. The Hubbard-Stratonovich transformation converts an interacting Hamiltonian into a non-interacting Hamiltonian in a time-dependent stochastic field, allowing calculation of the resulting functional integral by Monte Carlo methods. The method is presented in a sufficiently general form to be applicable to any Hamiltonian with one- and two-body terms, with special reference to the Heisenberg model and one- and many-band Hubbard models. Many physical correlation functions can be related to correlation functions of the auxiliary field; general results are given here. Issues relating to the choice of auxiliary fields are addressed; operator product identities change the relative dimensionalities of the attractive and repulsive parts of the interaction. Frequently the integrand is not positive-definite, rendering numerical evaluation unstable. If the auxiliary field violates time-reversal invariance, the integrand is complex and this sign problem becomes a phase problem. The origin of this sign or phase is examined from a number of geometrical and other viewpoints and illustrated by simple examples: the phase problem by the spin (1/2) Heisenberg model, and the sign problem by the attractive SU(N) Hubbard model on a triangular molecule with negative hopping integrals. In the latter case, widely studied in the Jahn-Teller literature, the sign is due neither to fermions nor spin, but to frustration. This system is used to illustrate a number of suggested interpretations of the sign problem.