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:1401.5661 · Physical Review A

Number-conserving approaches to n -component Bose-Einstein condensates

Peter Mason, S. A. Gardiner

Published Apr 23, 2014• 10 citations
arXivHTMLPDF

Abstract

We develop the number-conserving approach, which has previously been used in a single-component Bose-Einstein condensed dilute atomic gas, to describe consistent coupled condensate and noncondensate number dynamics, to an $n$-component condensate. The resulting system of equations is comprised, for each component, of a generalized Gross-Pitaevskii equation coupled to modified Bogoliubov--de Gennes equations. Lower order approximations yield general formulations for multicomponent Gross-Pitaevskii equations, and systems of multicomponent Gross-Pitaevskii equations coupled to multicomponent modified number-conserving Bogoliubov--de Gennes equations. The analysis is left general, such that, in the $n$-component condensate, there may or may not be mutually coherent components. An expansion in powers of the ratio of noncondensate-to-condensate particle numbers for each coherent set is used to derive the self-consistent, second-order, dynamical equations of motion. The advantage of the analysis developed in this article is in its applications to dynamical instabilities that appear when two (or more) components are in conflict and where a significant noncondensed fraction of atoms is expected to appear.

DOI
10.1103/physreva.89.043617

Related papers

Number-conserving approach to a minimal self-consistent treatment of condensate and noncondensate dynamics in a degenerate Bose gas

S. A. Gardiner, S. A. Morgan

We describe a number-conserving approach to the dynamics of Bose-Einstein condensed dilute atomic gases. This builds upon the works of Gardiner [Phys. Rev. A 56, 1414 (1997)] and Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. We consider what is effectively an expansion in powers of the ratio of noncondensate to condensate particle numbers, rather than inverse powers of the total number of particles. This requires the number of condensate particles to be a majority, but not necessarily almost equal to the total number of particles in the system. We argue that a second-order treatment of the relevant dynamical equations of motion is the minimum order necessary to provide consistent coupled condensate and noncondensate number dynamics for a finite total number of particles, and show that such a second-order treatment is provided by a suitably generalized Gross-Pitaevskii equation, coupled to the Castin-Dum number-conserving formulation of the Bogoliubov--de Gennes equations. The necessary equations of motion can be generated from an approximate third-order Hamiltonian, which effectively reduces to second order in the steady state. Such a treatment as described here is suitable for dynamics occurring at finite temperature, where there is a significant noncondensate fraction from the outset, or dynamics leading to dynamical instabilities, where depletion of the condensate can also lead to a significant noncondensate fraction, even if the noncondensate fraction is initially negligible.

arXiv:cond-mat/0610623•Physical Review A•Apr 23, 2007•70 citations

Second-order number-conserving description of nonequilibrium dynamics in finite-temperature Bose-Einstein condensates

T. P. Billam, Peter Mason, S. A. Gardiner

While the Gross-Pitaevskii equation is well established as the canonical dynamical description of atomic Bose-Einstein condensates (BECs) at zero temperature, describing the dynamics of BECs at finite temperatures remains a difficult theoretical problem, particularly when considering low-temperature, nonequilibrium systems in which depletion of the condensate occurs dynamically as a result of external driving. In this paper, we describe a fully time-dependent numerical implementation of a second-order, number-conserving description of finite-temperature BEC dynamics. This description consists of equations of motion describing the coupled dynamics of the condensate and noncondensate fractions in a self-consistent manner, and is ideally suited for the study of low-temperature, nonequilibrium, driven systems. The $\ensuremath{\delta}$-kicked-rotor BEC provides a prototypical example of such a system, and we demonstrate the efficacy of our numerical implementation by investigating its dynamics at finite temperature. We demonstrate that the qualitative features of the system dynamics at zero temperature are generally preserved at finite temperatures, and predict a quantitative finite-temperature shift of resonance frequencies which would be relevant for, and could be verified by, future experiments.

arXiv:1207.2821•Physical Review A•Mar 27, 2013•25 citations

Particle-number-conserving Bogoliubov approximation for Bose-Einstein condensates using extended catalytic states

Jiang Zhang, Carlton M. Caves

We encode the many-body wave function of a Bose-Einstein condensate (BEC) in the $N$-particle sector of an extended catalytic state. This catalytic state is a coherent state for the condensate mode and an arbitrary state for the modes orthogonal to the condensate mode. Going to a time-dependent interaction picture where the state of the condensate mode is displaced to the vacuum, we can organize the effective Hamiltonian by powers of ${N}^{\ensuremath{-}1/2}$. Requiring the terms of order ${N}^{1/2}$ to vanish gives the Gross-Pitaevskii equation. Going to the next order, ${N}^{0}$, we derive equations for the number-conserving Bogoliubov approximation, first given by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. In contrast to other approaches, ours is well suited to calculating the state evolution in the Schr\"odinger picture; moreover, it is straightforward to generalize our method to multicomponent BECs and to higher-order corrections.

arXiv:1503.02132•Physical Review A•Mar 14, 2016•6 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

Quantitative derivation of the two-component Gross–Pitaevskii equation in the hard-core limit with uniform-in-time convergence rate

Jacky J. Chong, Jinyeop Lee, Zhiwei Sun

Abstract We derive, in the Gross–Pitaevskii (GP) limit, the time-dependent two-component GP equation as an effective description of the dynamics of a dilute two-component Bose gas near its ground state, which exhibits a two-component Bose–Einstein condensate. Our main result establishes a uniform-in-time bound on the convergence rate between the many-body dynamics and the effective description, explicitly quantified in terms of the particle number N . This improves upon the works of Michelangeli and Olgiati 2014 J. Nonlinear Math. Phys. 24 466–4; Olgiati 2017 Effective Non-Linear Dynamics of Binary Condensates and Open Problems pp 239–56) by providing a sharper, N -dependent, time-independent convergence rate, and also implies a uniform-in-time bound for the one-component case. Furthermore, our approach extends the framework of Benedikter et al (2015 Commun. Pure Appl. Math., 68 1399–482) to the multi-component Bose gas in the hard-core limit setting. More specifically, we develop the necessary Bogoliubov theory to analyze the dynamics of multi-component Bose gases in the GP regime.

arXiv:2501.18787•Nonlinearity•Aug 25, 2026•1 citation

Nonequilibrium kinetic theory for trapped binary condensates

Matthew Edmonds, Kean Loon Lee, Nick P. Proukakis

We derive a nonequilibrium finite-temperature kinetic theory for a binary mixture of two interacting atomic Bose-Einstein condensates and use it to explore the degree of hydrodynamicity attainable in realistic experimental geometries. Based on the standard separation-of-time-scales argument of kinetic theory, the dynamics of the condensates of the multicomponent system are shown to be described by dissipative Gross-Pitaevskii equations self-consistently coupled to corresponding quantum Boltzmann equations for the noncondensate atoms: On top of the usual mean-field contributions, our scheme identifies a total of eight distinct collisional processes, whose dynamical interplay is expected to be responsible for the system's equilibration. In order to provide their first characterization, we perform a detailed numerical analysis of the role of trap frequency and geometry on collisional rates for experimentally accessible mixtures of $^{87}\mathrm{Rb}\ensuremath{-}^{41}\mathrm{K}$ and $^{87}\mathrm{Rb}\ensuremath{-}^{85}\mathrm{Rb}$, discussing the extent to which the system may approach the hydrodynamic regime with regard to some of those processes as a guide for future experimental investigations of ultracold Bose gas mixtures.

arXiv:1507.05078•Physical Review A•Dec 3, 2015•6 citations
We develop the number-conserving approach, which has previously been used in a single-component Bose-Einstein condensed dilute atomic gas, to describe consistent coupled condensate and noncondensate number dynamics, to an nnn-component condensate. The resulting system of equations is comprised, for each component, of a generalized Gross-Pitaevskii equation coupled to modified Bogoliubov--de Gennes equations. Lower order approximations yield general formulations for multicomponent Gross-Pitaevskii equations, and systems of multicomponent Gross-Pitaevskii equations coupled to multicomponent modified number-conserving Bogoliubov--de Gennes equations. The analysis is left general, such that, in the nnn-component condensate, there may or may not be mutually coherent components. An expansion in powers of the ratio of noncondensate-to-condensate particle numbers for each coherent set is used to derive the self-consistent, second-order, dynamical equations of motion. The advantage of the analysis developed in this article is in its applications to dynamical instabilities that appear when two (or more) components are in conflict and where a significant noncondensed fraction of atoms is expected to appear.
We describe a number-conserving approach to the dynamics of Bose-Einstein condensed dilute atomic gases. This builds upon the works of Gardiner [Phys. Rev. A 56, 1414 (1997)] and Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. We consider what is effectively an expansion in powers of the ratio of noncondensate to condensate particle numbers, rather than inverse powers of the total number of particles. This requires the number of condensate particles to be a majority, but not necessarily almost equal to the total number of particles in the system. We argue that a second-order treatment of the relevant dynamical equations of motion is the minimum order necessary to provide consistent coupled condensate and noncondensate number dynamics for a finite total number of particles, and show that such a second-order treatment is provided by a suitably generalized Gross-Pitaevskii equation, coupled to the Castin-Dum number-conserving formulation of the Bogoliubov--de Gennes equations. The necessary equations of motion can be generated from an approximate third-order Hamiltonian, which effectively reduces to second order in the steady state. Such a treatment as described here is suitable for dynamics occurring at finite temperature, where there is a significant noncondensate fraction from the outset, or dynamics leading to dynamical instabilities, where depletion of the condensate can also lead to a significant noncondensate fraction, even if the noncondensate fraction is initially negligible.
While the Gross-Pitaevskii equation is well established as the canonical dynamical description of atomic Bose-Einstein condensates (BECs) at zero temperature, describing the dynamics of BECs at finite temperatures remains a difficult theoretical problem, particularly when considering low-temperature, nonequilibrium systems in which depletion of the condensate occurs dynamically as a result of external driving. In this paper, we describe a fully time-dependent numerical implementation of a second-order, number-conserving description of finite-temperature BEC dynamics. This description consists of equations of motion describing the coupled dynamics of the condensate and noncondensate fractions in a self-consistent manner, and is ideally suited for the study of low-temperature, nonequilibrium, driven systems. The $\ensuremath{\delta}$-kicked-rotor BEC provides a prototypical example of such a system, and we demonstrate the efficacy of our numerical implementation by investigating its dynamics at finite temperature. We demonstrate that the qualitative features of the system dynamics at zero temperature are generally preserved at finite temperatures, and predict a quantitative finite-temperature shift of resonance frequencies which would be relevant for, and could be verified by, future experiments.
We encode the many-body wave function of a Bose-Einstein condensate (BEC) in the NNN-particle sector of an extended catalytic state. This catalytic state is a coherent state for the condensate mode and an arbitrary state for the modes orthogonal to the condensate mode. Going to a time-dependent interaction picture where the state of the condensate mode is displaced to the vacuum, we can organize the effective Hamiltonian by powers of ${N}^{\ensuremath{-}1/2}$. Requiring the terms of order N1/2{N}^{1/2}N1/2 to vanish gives the Gross-Pitaevskii equation. Going to the next order, N0{N}^{0}N0, we derive equations for the number-conserving Bogoliubov approximation, first given by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. In contrast to other approaches, ours is well suited to calculating the state evolution in the Schr\"odinger picture; moreover, it is straightforward to generalize our method to multicomponent BECs and to higher-order corrections.
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.
Abstract We derive, in the Gross–Pitaevskii (GP) limit, the time-dependent two-component GP equation as an effective description of the dynamics of a dilute two-component Bose gas near its ground state, which exhibits a two-component Bose–Einstein condensate. Our main result establishes a uniform-in-time bound on the convergence rate between the many-body dynamics and the effective description, explicitly quantified in terms of the particle number N . This improves upon the works of Michelangeli and Olgiati 2014 J. Nonlinear Math. Phys. 24 466–4; Olgiati 2017 Effective Non-Linear Dynamics of Binary Condensates and Open Problems pp 239–56) by providing a sharper, N -dependent, time-independent convergence rate, and also implies a uniform-in-time bound for the one-component case. Furthermore, our approach extends the framework of Benedikter et al (2015 Commun. Pure Appl. Math., 68 1399–482) to the multi-component Bose gas in the hard-core limit setting. More specifically, we develop the necessary Bogoliubov theory to analyze the dynamics of multi-component Bose gases in the GP regime.
We derive a nonequilibrium finite-temperature kinetic theory for a binary mixture of two interacting atomic Bose-Einstein condensates and use it to explore the degree of hydrodynamicity attainable in realistic experimental geometries. Based on the standard separation-of-time-scales argument of kinetic theory, the dynamics of the condensates of the multicomponent system are shown to be described by dissipative Gross-Pitaevskii equations self-consistently coupled to corresponding quantum Boltzmann equations for the noncondensate atoms: On top of the usual mean-field contributions, our scheme identifies a total of eight distinct collisional processes, whose dynamical interplay is expected to be responsible for the system's equilibration. In order to provide their first characterization, we perform a detailed numerical analysis of the role of trap frequency and geometry on collisional rates for experimentally accessible mixtures of $^{87}\mathrm{Rb}\ensuremath{-}^{41}\mathrm{K}$ and $^{87}\mathrm{Rb}\ensuremath{-}^{85}\mathrm{Rb}$, discussing the extent to which the system may approach the hydrodynamic regime with regard to some of those processes as a guide for future experimental investigations of ultracold Bose gas mixtures.