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/0011431 · Physical Review Letters

Simulations of Bose Fields at Finite Temperature

Matthew J. Davis, S. A. Morgan, Keith R. Burnett

Published Sep 27, 2001• 228 citations
arXivHTMLPDF

Abstract

We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.

DOI
10.1103/physrevlett.87.160402

Related papers

Simulations of thermal Bose fields in the classical limit

Matthew J. Davis, S. A. Morgan, Keith R. Burnett

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

arXiv:cond-mat/0201571•Physical Review A•Nov 25, 2002•103 citations

Microcanonical temperature for a classical field: Application to Bose-Einstein condensation

Matthew J. Davis, S. A. Morgan

We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a uv cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behavior of the specific heat.

arXiv:cond-mat/0307155•Physical Review A•Nov 24, 2003•62 citations

Gross–Pitaevskii description of superfluid dynamics at finite temperature: A short review of recent results

Marc Brächet

The Gross–Pitaevskii equation (GPE) describes the dynamics of superflows and Bose–Einstein Condensates (BEC) at very low temperature. When a truncation of Fourier modes is performed, the resulting truncated GPE (TGPE) can also describe the correct thermal behavior of a Bose gas, as long as all relevant modes are highly occupied [M.J. Davis, S.A. Morgan, K. Burnett, Simulations of Bose fields at finite temperature, Phys. Rev. Lett. 87 (16) (2001) 160402]. We review some of our groupʼs recent GPE- and TGPE-based numerical studies of superfluid dynamics and BEC stability. The relations with experiments are discussed.

Comptes Rendus Physique•Nov 1, 2012•12 citations

Projected Gross-Pitaevskii equation for harmonically confined Bose gases at finite temperature

P. B. Blakie, Matthew J. Davis

We extend the projected Gross-Pitaevskii equation formalism of Davis et al. [Phys. Rev. Lett. 87, 160402 (2001)] to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.

arXiv:cond-mat/0410496•Physical Review A•Dec 7, 2005•106 citations

Dynamics of thermal Bose fields in the classical limit

Matthew J. Davis, R. J. Ballagh, Keith R. Burnett

We develop an approximate formalism suitable for performing simulations of the thermal dynamics of interacting Bose gases. The method is based on the observation that when the lowest-energy modes of the Bose field operator are highly occupied, they may be treated classically to a good approximation. We derive a finite-temperature Gross-Pitaevskii equation (FTGPE) for these modes which is coupled to an effective reservoir described by quantum kinetic theory. We discuss each of the terms that arise in this GPE, and their relevance to experimental systems. We then describe a simpler projected GPE that may be useful in simulating thermal Bose condensates. This classical method could be applied to other Bose fields.

arXiv:cond-mat/0107515•Journal of Physics B Atomic Molecular and Optical Physics•Nov 13, 2001•76 citations

Mean-field theory for excitations of trapped Bose condensates at finite temperatures

M. Rusch, Keith R. Burnett

We describe a method for studying the dynamics of a condensate in a weakly interacting, dilute Bose gas which is based on a finite-temperature extension of the Gross-Pitaevskii equation. We present a general finite-temperature linear response formalism for nonuniform systems and give explicit results for homogeneous ones, where we can link our treatment to well-known results. This treatment leads to a consistent ${V}^{2}$ calculation of the excitations for a weakly interacting Bose gas, including Landau and Beliaev damping as well as collisional lifetimes. We do not make use of the Popov approximation, but include pair and triplet anomalous averages to second order in the interaction potential and thus obtain a gapless theory for the excitations of a weakly interacting homogeneous gas. We also show how the formalism can be used to obtain self-consistent theories. We explicitly derive gapless equations for a weakly coupled finite-temperature condensate with a static thermal component and give a clear criterion of their validity in trapped systems.

Physical Review A•May 1, 1999•25 citations
We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.
We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.
We show that the projected Gross-Pitaevskii equation (PGPE) can be mapped exactly onto Hamilton's equations of motion for classical position and momentum variables. Making use of this mapping, we adapt techniques developed in statistical mechanics to calculate the temperature and chemical potential of a classical Bose field in the microcanonical ensemble. We apply the method to simulations of the PGPE, which can be used to represent the highly occupied modes of Bose condensed gases at finite temperature. The method is rigorous, valid beyond the realms of perturbation theory, and agrees with an earlier method of temperature measurement for the same system. Using this method we show that the critical temperature for condensation in a homogeneous Bose gas on a lattice with a uv cutoff increases with the interaction strength. We discuss how to determine the temperature shift for the Bose gas in the continuum limit using this type of calculation, and obtain a result in agreement with more sophisticated Monte Carlo simulations. We also consider the behavior of the specific heat.
The Gross–Pitaevskii equation (GPE) describes the dynamics of superflows and Bose–Einstein Condensates (BEC) at very low temperature. When a truncation of Fourier modes is performed, the resulting truncated GPE (TGPE) can also describe the correct thermal behavior of a Bose gas, as long as all relevant modes are highly occupied [M.J. Davis, S.A. Morgan, K. Burnett, Simulations of Bose fields at finite temperature, Phys. Rev. Lett. 87 (16) (2001) 160402]. We review some of our groupʼs recent GPE- and TGPE-based numerical studies of superfluid dynamics and BEC stability. The relations with experiments are discussed.
We extend the projected Gross-Pitaevskii equation formalism of Davis et al. [Phys. Rev. Lett. 87, 160402 (2001)] to the experimentally relevant case of thermal Bose gases in harmonic potentials and outline a robust and accurate numerical scheme that can efficiently simulate this system. We apply this method to investigate the equilibrium properties of the harmonically trapped three-dimensional projected Gross-Pitaevskii equation at finite temperature and consider the dependence of condensate fraction, position, and momentum distributions and density fluctuations on temperature. We apply the scheme to simulate an evaporative cooling process in which the preferential removal of high-energy particles leads to the growth of a Bose-Einstein condensate. We show that a condensate fraction can be inferred during the dynamics even in this nonequilibrium situation.
We develop an approximate formalism suitable for performing simulations of the thermal dynamics of interacting Bose gases. The method is based on the observation that when the lowest-energy modes of the Bose field operator are highly occupied, they may be treated classically to a good approximation. We derive a finite-temperature Gross-Pitaevskii equation (FTGPE) for these modes which is coupled to an effective reservoir described by quantum kinetic theory. We discuss each of the terms that arise in this GPE, and their relevance to experimental systems. We then describe a simpler projected GPE that may be useful in simulating thermal Bose condensates. This classical method could be applied to other Bose fields.
We describe a method for studying the dynamics of a condensate in a weakly interacting, dilute Bose gas which is based on a finite-temperature extension of the Gross-Pitaevskii equation. We present a general finite-temperature linear response formalism for nonuniform systems and give explicit results for homogeneous ones, where we can link our treatment to well-known results. This treatment leads to a consistent V2{V}^{2}V2 calculation of the excitations for a weakly interacting Bose gas, including Landau and Beliaev damping as well as collisional lifetimes. We do not make use of the Popov approximation, but include pair and triplet anomalous averages to second order in the interaction potential and thus obtain a gapless theory for the excitations of a weakly interacting homogeneous gas. We also show how the formalism can be used to obtain self-consistent theories. We explicitly derive gapless equations for a weakly coupled finite-temperature condensate with a static thermal component and give a clear criterion of their validity in trapped systems.