Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0009220 · Physical Review A

Quantum corrections to the ground-state energy of a trapped Bose-Einstein condensate: A diffusion Monte Carlo calculation

D. Blume, Chris H. Greene

Published May 4, 2001• 73 citations
arXivHTMLPDF

Abstract

The diffusion Monte Carlo method is applied to describe a trapped atomic Bose-Einstein condensate at zero temperature, fully quantum mechanically and nonperturbatively. For low densities, ${n(0)a}^{3}<~2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}3}$ $[n(0),$ peak density; a, s-wave scattering length], our calculations confirm that the exact ground-state energy for a sum of two-body interactions depends to a good approximation on only the atomic-physics parameter a, and no other details of the two-body model potential. Corrections to the mean-field Gross-Pitaevskii energy range from being essentially negligible to about 20% for $N=2--50$ particles in the trap with positive s-wave scattering length $a=100--10000\mathrm{a}.\mathrm{u}.$ Our numerical calculations confirm that inclusion of an additional effective potential term in the mean-field equation, which accounts for quantum fluctuations [see, e.g., E. Braaten and A. Nieto, Phys. Rev. B 56, 14 745 (1997)], leads to a greatly improved description of trapped Bose gases.

DOI
10.1103/physreva.63.063601

Related papers

Correlation effects in the ground state of trapped atomic Bose gases

Wirawan Purwanto, Shiwei Zhang

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 (LDA).

arXiv:cond-mat/0506219•Physical Review A•Nov 10, 2005•30 citations

Ground state of a homogeneous Bose gas: A diffusion Monte Carlo calculation

S. Giorgini, Jordi Boronat, J. Casulleras

We use a diffusion Monte Carlo method to calculate the lowest-energy state of a uniform gas of bosons interacting through different model potentials, both strictly repulsive and with an attractive well. We explicitly verify that at low density the energy per particle follows a universal behavior fixed by the gas parameter ${\mathrm{na}}^{3}.$ In the regime of densities typical for experiments in trapped Bose-condensed gases, the corrections to the mean-field energy greatly exceed the differences due to the details of the potential.

arXiv:cond-mat/9902185•Physical Review A•Dec 1, 1999•149 citations

BOSE-EINSTEIN CONDENSATION IN TRAPS: A QUANTUM MONTE CARLO STUDY

Jonathan L. DuBois

We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.

•Jan 1, 2002•1 citation

Quantum corrections to the ground state of a trapped Bose-Einstein condensate

Eric Braaten, Agustín Nieto

In the mean-field approximation, the number density $\ensuremath{\rho}(\mathbf{r})$ for the ground state of a Bose-Einstein condensate trapped by an external potential $V(\mathbf{r})$ satisfies a classical field equation called the Gross-Pitaevskii equation. We show that quantum corrections to $\ensuremath{\rho}$ are dominated by quantum fluctuations with wavelengths of order $1/\sqrt{\ensuremath{\rho}a}$, where $a$ is the $S$-wave scattering length. By expanding the equations for the Hartree-Fock approximation to second order in the gradient expansion, we derive local correction terms to the Gross-Pitaevskii equation for $\ensuremath{\rho}$ that take into account the dominant effects of quantum fluctuations. We also show that the gradient expansions for the density and for the condensate break down at fourth order and at second order, respectively.

arXiv:cond-mat/9707199•Physical review. B, Condensed matter•Dec 1, 1997•74 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, 2005•0 citations

Beyond the Fermi pseudopotential: A modified Gross-Pitaevskii equation

Haixiang Fu, Yuzhu Wang, Bo Gao

We present an effective potential and the corresponding modified Gross-Pitaevskii equation, which account for the energy dependence of the two-body scattering amplitude through an effective-range expansion. For the ground-state energy of a trapped condensate, the theory leads to what we call a shape-dependent confinement correction that improves agreements with diffusion Monte Carlo calculations. The theory illustrates, for relatively strong confinement and/or high density, how the shape dependence on atom-atom interaction can come into play in a many-atom quantum system.

arXiv:cond-mat/0302331•Physical Review A•May 27, 2003•73 citations
The diffusion Monte Carlo method is applied to describe a trapped atomic Bose-Einstein condensate at zero temperature, fully quantum mechanically and nonperturbatively. For low densities, n(0)a3< 2\ifmmode×\else\texttimes\fi10\ensuremath−3{n(0)a}^{3}<~2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}3}n(0)a3< 2\ifmmode×\else\texttimes\fi10\ensuremath−3 [n(0),[n(0),[n(0), peak density; a, s-wave scattering length], our calculations confirm that the exact ground-state energy for a sum of two-body interactions depends to a good approximation on only the atomic-physics parameter a, and no other details of the two-body model potential. Corrections to the mean-field Gross-Pitaevskii energy range from being essentially negligible to about 20% for N=2−−50N=2--50N=2−−50 particles in the trap with positive s-wave scattering length a=100−−10000a.u.a=100--10000\mathrm{a}.\mathrm{u}.a=100−−10000a.u. Our numerical calculations confirm that inclusion of an additional effective potential term in the mean-field equation, which accounts for quantum fluctuations [see, e.g., E. Braaten and A. Nieto, Phys. Rev. B 56, 14 745 (1997)], leads to a greatly improved description of trapped 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 (LDA).
We use a diffusion Monte Carlo method to calculate the lowest-energy state of a uniform gas of bosons interacting through different model potentials, both strictly repulsive and with an attractive well. We explicitly verify that at low density the energy per particle follows a universal behavior fixed by the gas parameter na3.{\mathrm{na}}^{3}.na3. In the regime of densities typical for experiments in trapped Bose-condensed gases, the corrections to the mean-field energy greatly exceed the differences due to the details of the potential.
We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.
In the mean-field approximation, the number density $\ensuremath{\rho}(\mathbf{r})$ for the ground state of a Bose-Einstein condensate trapped by an external potential V(r)V(\mathbf{r})V(r) satisfies a classical field equation called the Gross-Pitaevskii equation. We show that quantum corrections to $\ensuremath{\rho}$ are dominated by quantum fluctuations with wavelengths of order $1/\sqrt{\ensuremath{\rho}a}$, where aaa is the SSS-wave scattering length. By expanding the equations for the Hartree-Fock approximation to second order in the gradient expansion, we derive local correction terms to the Gross-Pitaevskii equation for $\ensuremath{\rho}$ that take into account the dominant effects of quantum fluctuations. We also show that the gradient expansions for the density and for the condensate break down at fourth order and at second order, respectively.
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.
We present an effective potential and the corresponding modified Gross-Pitaevskii equation, which account for the energy dependence of the two-body scattering amplitude through an effective-range expansion. For the ground-state energy of a trapped condensate, the theory leads to what we call a shape-dependent confinement correction that improves agreements with diffusion Monte Carlo calculations. The theory illustrates, for relatively strong confinement and/or high density, how the shape dependence on atom-atom interaction can come into play in a many-atom quantum system.