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.
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.
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.
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.
A recent paper of Gardiner [Phys. Rev. A 56, 1414 (1997)] introduces a particle-number-conserving Bogoliubov method for the excitation spectrum of a Bose-condensed gas, for use in theories of recently experimentally produced trapped atomic Bose condensates. Gardiner's approach is compared and contrasted with the 1959 Girardeau-Arnowitt theory [Phys. Rev. 113, 755 (1959)], to which it is closely related and which is also fully number conserving. The number-conserving Bogoliubov quasiparticle operators of the Girardeau-Arnowitt theory satisfy Bose commutation relations exactly so long as states with the condensate totally depleted are neglected, whereas those of Gardiner satisfy Bose commutation relations only in an approximation that deteriorates progressively as the condensate is depleted.
This paper extends an earlier quantum kinetics treatment for dilute, weakly interacting, partially Bose-Einstein condensed gases, presented by the author elsewhere (Proukakis N P and Burnett K 1996 J. Res. Natl Inst. Stand. Technol. 101 457), by consistently treating the dynamics of the uncondensed atoms to the same level of approximation as the condensed atoms. Our method is based on a hierarchy of coupled equations of motion for the condensate mean field and fluctuations around this mean field, truncated to second order in the (effective two-body) interatomic potential, and with suitable decoupling approximations for higher-order correlations. By applying perturbation theory in the Hartree-Fock-Bogoliubov basis, we re-derive the quantum kinetic theory of Walser et al (Walser R, Williams J, Cooper J and Holland M 1999 Phys. Rev. A 59 3878), which further indicates the consistency of our treatment with the Kadanoff-Baym non-equilibrium Green function formalism for trapped gases.
This paper extends an earlier quantum kinetics treatment for dilute, weakly interacting, partially Bose-Einstein condensed gases, presented by the author elsewhere (Proukakis N P and Burnett K 1996 J. Res. Natl Inst. Stand. Technol. 101 457), by consistently treating the dynamics of the uncondensed atoms to the same level of approximation as the condensed atoms. Our method is based on a hierarchy of coupled equations of motion for the condensate mean field and fluctuations around this mean field, truncated to second order in the (effective two-body) interatomic potential, and with suitable decoupling approximations for higher-order correlations. By applying perturbation theory in the Hartree-Fock-Bogoliubov basis, we re-derive the quantum kinetic theory of Walser et al (Walser R, Williams J, Cooper J and Holland M 1999 Phys. Rev. A 59 3878), which further indicates the consistency of our treatment with the Kadanoff-Baym non-equilibrium Green function formalism for trapped gases.
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.