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/0205421 · Physical Review A

Modeling Bose-Einstein condensed gases at finite temperatures withN-body simulations

B. Jackson, E. Zaremba

Published Sep 17, 2002• 79 citations
arXivHTMLPDF

Abstract

We consider a model of a dilute Bose-Einstein condensed gas at finite temperatures, where the condensate coexists in a trap with a cloud of thermal excitations. Within the Zaremba, Nikuni, and Griffin formalism, the dynamics of the condensate is described by a generalized Gross-Pitaevskii equation, while the thermal cloud is represented by a semiclassical kinetic equation. Our numerical approach simulates the kinetic equation using a cloud of representative test particles, while collisions are treated by means of a Monte Carlo sampling technique. A full three-dimensional split-operator fast Fourier transform method is used to evolve the condensate wave function. We give details regarding the numerical methods used and discuss simulations carried out to test the accuracy of the numerics. We use this scheme to simulate the monopole mode in a spherical trap. The dynamical coupling between the condensate and thermal cloud is responsible for frequency shifts and damping of the condensate collective mode. We compare our results to previous theoretical approaches, not only to confirm the reliability of our numerical scheme, but also to check the validity of approximations which have been used in the past.

DOI
10.1103/physreva.66.033606

Related papers

Dynamical Simulations of Trapped Bose Gases at Finite Temperatures

Jackson, B., E. Zaremba

In this paper, we develop a numerical procedure for investigating the dynamics of trapped Bose gases based on the ZGN theory. The dynamical equations used consist of a generalized Gross-Pitaevskii equation for the condensate order parameter and a semiclassical kinetic equation for the thermal cloud. The former is solved using a fast Fourier transform split-operator technique while the Boltzmann equation is treated by means of N-body simulations. The two components are coupled by mean fields as well as collisional processes that transfer atoms between the two. This scheme has been applied to a model equilibration problem, dipole oscillations in isotropic traps and scissors modes in anisotropic traps. In the case of the latter, the frequencies and damping rates of the condensate mode have been extracted from the simulations for a wide range of temperatures. Good agreement with recent experiments has been found.

arXiv:cond-mat/0106652•Jun 29, 2001•1 citation

Simulation Method for Evaporative Cooling of Trapped Bose Gases at Finite Temperatures

Emiko Arahata, Tetsuro Nikuni

Abstract We have devised a simulation method for studying the evaporative cooling process of trapped Bose-Einstein condensates at finite temperatures, leveraging the Zaremba-Nikuni-Griffin (ZNG) formalism. This framework incorporates the generalized Gross-Pitaevskii (GP) equation and a semiclassical kinetic equation governing the behavior of the thermal cloud. Our approach builds upon the numerical techniques pioneered by Jackson and Zaremba, which involve simulating the kinetic equation using test particles. A pivotal aspect of our method lies in emulating the evaporative cooling process by selectively removing test particles with high energy. Through extensive simulations, we demonstrate the efficacy of our approach in accurately describing condensate growth during evaporative cooling. Furthermore, we extend our simulations to investigate the formation of vortex lattices during evaporative cooling in the presence of a rotating thermal cloud. This comprehensive analysis provides valuable insights into the intricate dynamics of trapped Bose-Einstein condensates undergoing evaporative cooling processes.

arXiv:2310.03250•Springer proceedings in physics•Jan 1, 2026•0 citations

Finite-Temperature Simulations of the Scissors Mode in Bose-Einstein Condensed Gases

B. Jackson, E. Zaremba

The dynamics of a trapped Bose-condensed gas at finite temperatures is described by a generalized Gross-Pitaevskii equation for the condensate order parameter and a semiclassical kinetic equation for the thermal cloud, solved using N-body simulations. The two components are coupled by mean fields as well as collisional processes that transfer atoms between the two. We use this scheme to investigate scissors modes in anisotropic traps as a function of temperature. Frequency shifts and damping rates of the condensate mode are extracted, and are found to be in good agreement with recent experiments.

arXiv:cond-mat/0105465•Physical Review Letters•Aug 20, 2001•61 citations

Damping and revivals of collective oscillations in a finite-temperature model of trapped Bose-Einstein condensation

B. Jackson, Charles S. Adams

We utilize a two-gas model to simulate collective oscillations of a Bose-Einstein condensate at finite temperatures. The condensate is described using a generalized Gross-Pitaevskii equation, which is coupled to a thermal cloud modeled by a Monte Carlo algorithm. This allows us to include the collective dynamics of both the condensed and noncondensed components self-consistently. We simulate quadrupolar excitations, and measure the damping rate and frequency as a function of temperature. We also observe revivals in condensate oscillations at high temperatures, and in the thermal cloud at low temperature. Extensions of the model to include nonequilibrium effects and describe more complex phenomena are discussed.

arXiv:cond-mat/0009468•Physical Review A•Apr 17, 2001•16 citations

Finite-temperature theory of the scissors mode in a Bose gas using the moment method

Tetsuro Nikuni

We use a generalized Gross-Pitaevskii equation for the condensate and a semiclassical kinetic equation for the noncondensate atoms to discuss the scissors mode in a trapped Bose-condensed gas at finite temperatures. Both equations include the effect of collisions between the condensate and noncondensate atoms. We solve the coupled moment equations describing oscillations of the quadrupole moments of the condensate and noncondensate components to find the collective-mode frequencies and collisional damping rates as a function of temperature. Our calculations extend those of Gu\'ery-Odelin and Stringari at $T=0$ and in the normal phase. They complement the numerical results of Jackson and Zaremba, although Landau damping is left out of our approach. Our results are also used to calculate the quadrupole response function, which is related to the moment of inertia. It is shown explicitly that the moment of inertia of a trapped Bose gas at finite temperatures involves a sum of an irrotational component from the condensate and a rotational component from the thermal-cloud atoms.

arXiv:cond-mat/0109149•Physical Review A•Feb 11, 2002•15 citations

Coupled equations for the condensate and thermal cloud

Allan Griffin, Tetsuro Nikuni, E. Zaremba

In contrast with Chapter 2, in this chapter we include the dynamics of the thermal cloud. As noted in Chapter 1, we treat the noncondensate atoms using the simplest microscopic model approximation that captures the important physics. In particular, we consider only temperatures high enough ( T ≥ 0.4 T BEC ) that the noncondensate atoms can be described by a particle-like Hartree–Fock (HF) spectrum. To extend the analysis to very low temperatures is in principle straightforward (see Chapter 7). However, the details are more complicated since the excitations of the thermal cloud take on a collective aspect (i.e. a Bogoliubov-type quasiparticle spectrum must be used). In trapped Bose gases, the HF single-particle spectrum gives a good approximation down to much lower temperatures than in the case of uniform Bose gases, as first emphasized by Giorgini et al . (1997). In Section 3.1, we derive a generalized form of the Gross–Pitaevskii equation for the Bose order parameter Φ( r , t ) that is valid at finite temperatures. It involves terms that are coupled to the noncondensate component (the thermal cloud) and thus its solution in general requires one to know the equations of motion for the dynamics of the noncondensate atoms. In Section 3.2 we restrict ourselves to finite temperatures high enough that the noncondensate atoms can be described by a quantum kinetic equation for the single-particle distribution function f ( p, r , t ). A detailed microscopic derivation of this kind of kinetic equation is given in Chapters 6 and 7 using the Kadanoff–Baym Green's function formalism. A characteristic feature of a Bose-condensed gas is that the kinetic equation governing f ( p, r , t ) involves a collision integral C 12 [ f , Φ] describing collisions between condensate and noncondensate atoms.

Cambridge University Press eBooks•Feb 19, 2009•0 citations
We consider a model of a dilute Bose-Einstein condensed gas at finite temperatures, where the condensate coexists in a trap with a cloud of thermal excitations. Within the Zaremba, Nikuni, and Griffin formalism, the dynamics of the condensate is described by a generalized Gross-Pitaevskii equation, while the thermal cloud is represented by a semiclassical kinetic equation. Our numerical approach simulates the kinetic equation using a cloud of representative test particles, while collisions are treated by means of a Monte Carlo sampling technique. A full three-dimensional split-operator fast Fourier transform method is used to evolve the condensate wave function. We give details regarding the numerical methods used and discuss simulations carried out to test the accuracy of the numerics. We use this scheme to simulate the monopole mode in a spherical trap. The dynamical coupling between the condensate and thermal cloud is responsible for frequency shifts and damping of the condensate collective mode. We compare our results to previous theoretical approaches, not only to confirm the reliability of our numerical scheme, but also to check the validity of approximations which have been used in the past.
In this paper, we develop a numerical procedure for investigating the dynamics of trapped Bose gases based on the ZGN theory. The dynamical equations used consist of a generalized Gross-Pitaevskii equation for the condensate order parameter and a semiclassical kinetic equation for the thermal cloud. The former is solved using a fast Fourier transform split-operator technique while the Boltzmann equation is treated by means of N-body simulations. The two components are coupled by mean fields as well as collisional processes that transfer atoms between the two. This scheme has been applied to a model equilibration problem, dipole oscillations in isotropic traps and scissors modes in anisotropic traps. In the case of the latter, the frequencies and damping rates of the condensate mode have been extracted from the simulations for a wide range of temperatures. Good agreement with recent experiments has been found.
Abstract We have devised a simulation method for studying the evaporative cooling process of trapped Bose-Einstein condensates at finite temperatures, leveraging the Zaremba-Nikuni-Griffin (ZNG) formalism. This framework incorporates the generalized Gross-Pitaevskii (GP) equation and a semiclassical kinetic equation governing the behavior of the thermal cloud. Our approach builds upon the numerical techniques pioneered by Jackson and Zaremba, which involve simulating the kinetic equation using test particles. A pivotal aspect of our method lies in emulating the evaporative cooling process by selectively removing test particles with high energy. Through extensive simulations, we demonstrate the efficacy of our approach in accurately describing condensate growth during evaporative cooling. Furthermore, we extend our simulations to investigate the formation of vortex lattices during evaporative cooling in the presence of a rotating thermal cloud. This comprehensive analysis provides valuable insights into the intricate dynamics of trapped Bose-Einstein condensates undergoing evaporative cooling processes.
The dynamics of a trapped Bose-condensed gas at finite temperatures is described by a generalized Gross-Pitaevskii equation for the condensate order parameter and a semiclassical kinetic equation for the thermal cloud, solved using N-body simulations. The two components are coupled by mean fields as well as collisional processes that transfer atoms between the two. We use this scheme to investigate scissors modes in anisotropic traps as a function of temperature. Frequency shifts and damping rates of the condensate mode are extracted, and are found to be in good agreement with recent experiments.
We utilize a two-gas model to simulate collective oscillations of a Bose-Einstein condensate at finite temperatures. The condensate is described using a generalized Gross-Pitaevskii equation, which is coupled to a thermal cloud modeled by a Monte Carlo algorithm. This allows us to include the collective dynamics of both the condensed and noncondensed components self-consistently. We simulate quadrupolar excitations, and measure the damping rate and frequency as a function of temperature. We also observe revivals in condensate oscillations at high temperatures, and in the thermal cloud at low temperature. Extensions of the model to include nonequilibrium effects and describe more complex phenomena are discussed.
We use a generalized Gross-Pitaevskii equation for the condensate and a semiclassical kinetic equation for the noncondensate atoms to discuss the scissors mode in a trapped Bose-condensed gas at finite temperatures. Both equations include the effect of collisions between the condensate and noncondensate atoms. We solve the coupled moment equations describing oscillations of the quadrupole moments of the condensate and noncondensate components to find the collective-mode frequencies and collisional damping rates as a function of temperature. Our calculations extend those of Gu\'ery-Odelin and Stringari at T=0T=0T=0 and in the normal phase. They complement the numerical results of Jackson and Zaremba, although Landau damping is left out of our approach. Our results are also used to calculate the quadrupole response function, which is related to the moment of inertia. It is shown explicitly that the moment of inertia of a trapped Bose gas at finite temperatures involves a sum of an irrotational component from the condensate and a rotational component from the thermal-cloud atoms.
In contrast with Chapter 2, in this chapter we include the dynamics of the thermal cloud. As noted in Chapter 1, we treat the noncondensate atoms using the simplest microscopic model approximation that captures the important physics. In particular, we consider only temperatures high enough ( T ≥ 0.4 T BEC ) that the noncondensate atoms can be described by a particle-like Hartree–Fock (HF) spectrum. To extend the analysis to very low temperatures is in principle straightforward (see Chapter 7). However, the details are more complicated since the excitations of the thermal cloud take on a collective aspect (i.e. a Bogoliubov-type quasiparticle spectrum must be used). In trapped Bose gases, the HF single-particle spectrum gives a good approximation down to much lower temperatures than in the case of uniform Bose gases, as first emphasized by Giorgini et al . (1997). In Section 3.1, we derive a generalized form of the Gross–Pitaevskii equation for the Bose order parameter Φ( r , t ) that is valid at finite temperatures. It involves terms that are coupled to the noncondensate component (the thermal cloud) and thus its solution in general requires one to know the equations of motion for the dynamics of the noncondensate atoms. In Section 3.2 we restrict ourselves to finite temperatures high enough that the noncondensate atoms can be described by a quantum kinetic equation for the single-particle distribution function f ( p, r , t ). A detailed microscopic derivation of this kind of kinetic equation is given in Chapters 6 and 7 using the Kadanoff–Baym Green's function formalism. A characteristic feature of a Bose-condensed gas is that the kinetic equation governing f ( p, r , t ) involves a collision integral C 12 [ f , Φ] describing collisions between condensate and noncondensate atoms.