Auxiliary field Monte Carlo for quantum many-body ground states

G Sugiya, S. E. Koonin

Published Jun 18, 2007• 0 citations

Abstract

Abstract This paper describes a different Q11C method for solving the Schrödinger equation, one unlike variational, diffusion, and Green’s function methods. It can be applied to many-body systems to produce solutions for ground states. For systems of fermions, antisymmetrization can be imposed to overcome the problem of node location, but the instability of the method is a serious shortcoming, so that the method as described is limited to treating relatively simple systems. Nevertheless, this auxiliary field quantum 1,fonte Carlo (AFMC) method was demonstrated to be successful for test cases reported in this paper, and in time has been improved to allow treatment of realistic systemsa. The method makes use of the imaginary-time propagation of an initial wavefunction toward the ground-state function. It requires use of a trial wavefunction which is the product of single-particle orbitals for bosons or an antisymmetrized product for fermions. A Hubbard-Stratonovich transformation of the propagator changes the many-body problem to a problem of a collection of single particles interacting with a fluctuating (auxiliary) field which replaces the many-body interactions. The authors developed an integral expression for the energy which, after discretization, could be evaluated by a Metropolis Monte Carlo integration forward in time. The algorithm was tested for one-dimensional systems of particles: 6 to 20 bosons with varied interactions, 2 to 6 neutron-proton pairs, others with repulsive potentials. Importance sampling was found to be critical and computational requirements to be large. For spin systems on a lattice the AFMC method was considered feasible, but for general multidimensional systems the method appeared to be limited by its computation requirements.

DOI: 10.1093/oso/9780195310108.003.0043

Related papers

Auxiliary Field Monte-Carlo for Quantum Many-Body Systems

Gayle A. Sugiyama

An algorithm is developed for determining the exact ground state properties of quantum many-body systems which is equally applicable to bosons and fermions. The Schroedinger eigenvalue equation for the ground state energy is recast into the form of a many-dimensional integral through the use of the Hubbard-Stratonovitch representation of the imaginary time many- body evolution operator. The resulting functional integral is then evaluated stochastically. The algorithm is tested for an exactly soluble boson system and is then extended to include fermions and repulsive potentials. Importance sampling is crucial to the success of the method, particularly for more complex systems. Improved computational efficiency is attained by performing the calculations in momentum space.

Jan 7, 20190 citations

Quantum Monte Carlo method for the ground state of many-boson systems

Wirawan Purwanto, Shiwei Zhang

We formulate a quantum Monte Carlo (QMC) method for calculating the ground state of many-boson systems. The method is based on a field-theoretical approach, and is closely related to existing fermion auxiliary-field QMC methods which are applied in several fields of physics. The ground-state projection is implemented as a branching random walk in the space of permanents consisting of identical single-particle orbitals. Any single-particle basis can be used, and the method is in principle exact. We illustrate this method with a trapped atomic boson gas, where the atoms interact via an attractive or repulsive contact two-body potential. We choose as the single-particle basis a real-space grid. We compare with exact results in small systems and arbitrarily sized systems of untrapped bosons with attractive interactions in one dimension, where analytical solutions exist. We also compare with the corresponding Gross-Pitaevskii (GP) mean-field calculations for trapped atoms, and discuss the close formal relation between our method and the GP approach. Our method provides a way to systematically improve upon GP while using the same framework, capturing interaction and correlation effects with a stochastic, coherent ensemble of noninteracting solutions. We discuss various algorithmic issues, including importance sampling and the back-propagation technique for computing observables, and illustrate them with numerical studies. We show results for systems with up to N approximately 400 bosons.

arXiv:physics/0403146Physical Review ENov 9, 200475 citations

Green’s function Monte Carlo for few fermion problems

David M. Arnow, Malvin H. Kalos, Michael A. Lee +1 more

The Green’s function Monte Carlo method used for obtaining exact solutions to the Schrödinger equation of boson systems is generalized to treat systems of several fermions. We show that when it is possible to select eigenfunctions of the Hamiltonian based on physical symmetries, the GFMC method can be used to yield the lowest energy state of that symmetry. In particular, the lowest totally antisymmetric eigenfunction, the fermion ground state, can be obtained. Calculations on several two- and three-body model problems show the method to be computationally feasible for few-body systems.

The Journal of Chemical PhysicsDec 1, 198299 citations

AUXILIARY FIELDS AND THE SIGN PROBLEM

John H. Samson

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.

International Journal of Modern Physics CJun 1, 19954 citations

Quantum Monte Carlo Method For Boson Ground States: Application To Trapped Bosons With Attractive And Repulsive Interactions

Wirawan Purwanto

We formulate a quantum Monte Carlo (QMC) method for calculating the ground state of many-boson systems. The method is based on a field-theoretical approach, and is closely related to existing fermion auxiliary-field QMC methods which are applied in several fields of physics. The ground-state projection is implemented as a branching random walk in the space of permanents consisting of identical single-particle orbitals. Any single-particle basis can be used, and the method is in principle exact. We apply this method to an atomic Bose gas, where the atoms interact via an attractive or repulsive contact two-body potential parametrized by the s-wave scattering length. We choose as the single-particle basis a real-space grid. We compare with exact results in small systems, and arbitrarily-sized systems of untrapped bosons with attractive interactions in one dimension, where analytical solutions exist. Our method provides a way to systematically improve upon the mean-field Gross-Pitaevskii (GP) method while using the same framework, capturing interaction and correlation effects with a stochastic, coherent ensemble of non-interacting solutions. to study the role of many-body correlations in the ground state, we examine the properties of the gas, such as the energetics, condensate fraction, and the density and momentum distributions as a function of the number of particles and the scattering length, both in the homogenous and trapped gases. Results are presented for systems with up to 1000 bosons. Comparing our results to the mean-field GP results, we find significant departure from mean field at large positive scattering lengths. The many-body correlations tend to increase the kinetic energy and reduce the interaction energy compared to GP. In the trapped gases, this results in a qualitatively different behavior as a function of the scattering length. Possible experimental observation is discussed.

Jan 1, 20050 citations

Computational Studies Of Strongly Correlated Quantum Matter

Hao Shi

The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.

Jan 1, 20170 citations