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/0108042 · Journal of Physics B Atomic Molecular and Optical Physics

An exact stochastic field method for the interacting Bose gas at thermal equilibrium

Iacopo Carusotto, Yvan Castin

Published Nov 23, 2001• 37 citations
arXivHTMLPDF

Abstract

We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.

DOI
10.1088/0953-4075/34/23/305

Related papers

Canonical ensemble of an interacting Bose gas: Stochastic matter fields and their coherence

Sigmund Heller, Walter T. Strunz

We present a novel quantum stochastic evolution equation for a matter field describing the canonical state of a weakly interacting ultracold Bose gas. In the ideal gas limit our approach is exact. This numerically very stable equation suppresses high-energy fluctuations exponentially, which enables us to describe condensed and thermal atoms within the same formalism. We present applications to ground-state occupation and fluctuations, density profile of ground-state and thermal cloud, and ground-state number statistics. Our main aim are spatial coherence properties which we investigate through the determination of interference contrast and spatial density correlations. Parameters are taken from actual experiments.

arXiv:1101.2617•Europhysics Letters (EPL)•Mar 1, 2013•3 citations

Stochastic field equation for the canonical ensemble of an ideal Bose gas

Sigmund Heller, Walter T. Strunz

We present a novel norm preserving stochastic evolution equation for a Bose field at finite temperature. Ensemble averages are quantum expectation values in the canonical ensemble. This numerically very stable equation suppresses high-energy fluctuations exponentially, which prevent cutoff problems from occurring. We present 3D simulations for an ideal gas in various trapping potentials and focus on ground-state occupation numbers and spatial correlation functions for a wide range of temperatures above and below the critical temperature. Although rigorously valid for non-interacting bosons only, we argue that weakly interacting Bose gases may also be amenable to this approach, in the usual mean-field approximation.

Journal of Physics B Atomic Molecular and Optical Physics•Mar 30, 2009•8 citations

Stochastic field equation for the canonical ensemble of a Bose gas

Sigmund Heller, Walter T. Strunz

We present a novel norm preserving stochastic evolution equation for a Bose field. Ensemble averages are quantum expectation values in the canonical ensemble. This numerically very stable equation suppresses high-energy fluctuations exponentially, preventing cutoff problems to occur. We present 3D simulations for an ideal gas in various trapping potentials focussing on ground state occupation numbers and spatial correlation functions for a wide range of temperatures above and below the critical temperature. Although rigorously valid for non-interacting Bosons only, we argue that weakly interacting Bose gases may also be amenable to this approach, in the usual mean-field approximation.

arXiv:0809.1002•Sep 5, 2008•0 citations

Comparison between microscopic methods for finite-temperature Bose gases

S. P. Cockburn, Antonio Negretti, Nick P. Proukakis +1 more

We analyze the equilibrium properties of a weakly interacting, trapped quasi-one-dimensional Bose gas at finite temperatures and compare different theoretical approaches. We focus in particular on two stochastic theories: a number-conserving Bogoliubov (NCB) approach and a stochastic Gross-Pitaevskii equation (SGPE) that have been extensively used in numerical simulations. Equilibrium properties like density profiles, correlation functions, and the condensate statistics are compared to predictions based upon a number of alternative theories. We find that due to thermal phase fluctuations, and the corresponding condensate depletion, the NCB approach loses its validity at relatively low temperatures. This can be attributed to the change in the Bogoliubov spectrum, as the condensate gets thermally depleted, and to large fluctuations beyond perturbation theory. Although the two stochastic theories are built on different thermodynamic ensembles (NCB, canonical; SGPE, grand-canonical), they yield the correct condensate statistics in a large Bose-Einstein condensate (BEC) (strong enough particle interactions). For smaller systems, the SGPE results are prone to anomalously large number fluctuations, well known for the grand-canonical, ideal Bose gas. Based on the comparison of the above theories to the modified Popov approach, we propose a simple procedure for approximately extracting the Penrose-Onsager condensate from first- and second-order correlation functions that is both computationally convenient and of potential use to experimentalists. This also clarifies the link between condensate and quasicondensate in the Popov theory of low-dimensional systems.

arXiv:1012.1512•Physical Review A•Apr 21, 2011•55 citations

Stochastic simulation of a finite-temperature one-dimensional Bose gas: From the Bogoliubov to the Tonks-Girardeau regime

Bernd G Schmidt, Lev I. Plimak, Michael Fleischhauer

We present an ab initio stochastic method for calculating thermal properties of a trapped, one-dimensional Bose gas covering the whole range from weak to strong interactions. Discretization of the problem results in a Bose-Hubbard--type Hamiltonian, whose imaginary time evolution is made computationally accessible by stochastic factorization of the kinetic energy. To achieve convergence for low enough temperatures such that quantum fluctuations are essential, the stochastic factorization is generalized to blocks, and ideas from density-matrix renormalization are employed. We compare our numerical results for density and first-order correlations with analytic predictions.

arXiv:cond-mat/0411332•Physical Review A•Apr 1, 2005•13 citations

Condensation ofNbosons. III. Analytical results for all higher moments of condensate fluctuations in interacting and ideal dilute Bose gases via the canonical ensemble quasiparticle formulation

Vitaly V. Kocharovsky, Vl. V. Kocharovsky, Marlan O. Scully

We consider equilibrium canonical-ensemble fluctuations of the Bose-Einstein condensate (BEC) in a trap using the particle-number conserving operator formalism of Girardeau and Arnowitt [Phys. Rev. 113, 755 (1959)] and the concept of the canonical-ensemble quasiparticles. We find explicit expressions for the characteristic function and all cumulants of the ground-state occupation statistics both for the dilute weakly interacting and ideal Bose gases. We present analytical formulas for the moments of the ground-state occupation fluctuations in the weakly interacting and ideal Bose gases in a box (``homogeneous gas'') as well as in the ideal Bose gas in an arbitrary power-law trap, in particular, in an arbitrary harmonic trap. We calculate the effect of Bogoliubov coupling between quasiparticles on suppression of the ground-state occupation fluctuations at moderate temperatures and their enhancement at very low temperatures. We clarify a crossover between ideal-gas and weakly interacting gas statistics, which is governed by a pair-correlation, squeezing mechanism. The important conclusion is that in most cases the ground-state occupation fluctuations are anomalously large and are not Gaussian even in the thermodynamic limit.

Physical Review A•Apr 13, 2000•68 citations
We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.
We present a novel quantum stochastic evolution equation for a matter field describing the canonical state of a weakly interacting ultracold Bose gas. In the ideal gas limit our approach is exact. This numerically very stable equation suppresses high-energy fluctuations exponentially, which enables us to describe condensed and thermal atoms within the same formalism. We present applications to ground-state occupation and fluctuations, density profile of ground-state and thermal cloud, and ground-state number statistics. Our main aim are spatial coherence properties which we investigate through the determination of interference contrast and spatial density correlations. Parameters are taken from actual experiments.
We present a novel norm preserving stochastic evolution equation for a Bose field at finite temperature. Ensemble averages are quantum expectation values in the canonical ensemble. This numerically very stable equation suppresses high-energy fluctuations exponentially, which prevent cutoff problems from occurring. We present 3D simulations for an ideal gas in various trapping potentials and focus on ground-state occupation numbers and spatial correlation functions for a wide range of temperatures above and below the critical temperature. Although rigorously valid for non-interacting bosons only, we argue that weakly interacting Bose gases may also be amenable to this approach, in the usual mean-field approximation.
We present a novel norm preserving stochastic evolution equation for a Bose field. Ensemble averages are quantum expectation values in the canonical ensemble. This numerically very stable equation suppresses high-energy fluctuations exponentially, preventing cutoff problems to occur. We present 3D simulations for an ideal gas in various trapping potentials focussing on ground state occupation numbers and spatial correlation functions for a wide range of temperatures above and below the critical temperature. Although rigorously valid for non-interacting Bosons only, we argue that weakly interacting Bose gases may also be amenable to this approach, in the usual mean-field approximation.
We analyze the equilibrium properties of a weakly interacting, trapped quasi-one-dimensional Bose gas at finite temperatures and compare different theoretical approaches. We focus in particular on two stochastic theories: a number-conserving Bogoliubov (NCB) approach and a stochastic Gross-Pitaevskii equation (SGPE) that have been extensively used in numerical simulations. Equilibrium properties like density profiles, correlation functions, and the condensate statistics are compared to predictions based upon a number of alternative theories. We find that due to thermal phase fluctuations, and the corresponding condensate depletion, the NCB approach loses its validity at relatively low temperatures. This can be attributed to the change in the Bogoliubov spectrum, as the condensate gets thermally depleted, and to large fluctuations beyond perturbation theory. Although the two stochastic theories are built on different thermodynamic ensembles (NCB, canonical; SGPE, grand-canonical), they yield the correct condensate statistics in a large Bose-Einstein condensate (BEC) (strong enough particle interactions). For smaller systems, the SGPE results are prone to anomalously large number fluctuations, well known for the grand-canonical, ideal Bose gas. Based on the comparison of the above theories to the modified Popov approach, we propose a simple procedure for approximately extracting the Penrose-Onsager condensate from first- and second-order correlation functions that is both computationally convenient and of potential use to experimentalists. This also clarifies the link between condensate and quasicondensate in the Popov theory of low-dimensional systems.
We present an ab initio stochastic method for calculating thermal properties of a trapped, one-dimensional Bose gas covering the whole range from weak to strong interactions. Discretization of the problem results in a Bose-Hubbard--type Hamiltonian, whose imaginary time evolution is made computationally accessible by stochastic factorization of the kinetic energy. To achieve convergence for low enough temperatures such that quantum fluctuations are essential, the stochastic factorization is generalized to blocks, and ideas from density-matrix renormalization are employed. We compare our numerical results for density and first-order correlations with analytic predictions.
We consider equilibrium canonical-ensemble fluctuations of the Bose-Einstein condensate (BEC) in a trap using the particle-number conserving operator formalism of Girardeau and Arnowitt [Phys. Rev. 113, 755 (1959)] and the concept of the canonical-ensemble quasiparticles. We find explicit expressions for the characteristic function and all cumulants of the ground-state occupation statistics both for the dilute weakly interacting and ideal Bose gases. We present analytical formulas for the moments of the ground-state occupation fluctuations in the weakly interacting and ideal Bose gases in a box (``homogeneous gas'') as well as in the ideal Bose gas in an arbitrary power-law trap, in particular, in an arbitrary harmonic trap. We calculate the effect of Bogoliubov coupling between quasiparticles on suppression of the ground-state occupation fluctuations at moderate temperatures and their enhancement at very low temperatures. We clarify a crossover between ideal-gas and weakly interacting gas statistics, which is governed by a pair-correlation, squeezing mechanism. The important conclusion is that in most cases the ground-state occupation fluctuations are anomalously large and are not Gaussian even in the thermodynamic limit.