The theoretical description of non-equilibrium Bose--Einstein condensate (BEC) is one of the main challenges in modern statistical physics and kinetics. The non-equilibrium nature of BEC makes it impossible to employ the well-established formalism of statistical mechanics. We develop a framework for the analytical description of a non-equilibrium phase transition to BEC that, in contrast to previously developed approaches, takes into account the infinite number of the states. We consider the limit of fast thermalization and obtain an analytical expression for the full density matrix of a non-equilibrium BEC which also covers the equilibrium case. For the particular cases of 2D and 3D, we investigate the non-equilibrium formation of BEC by finding the temperature dependence of the ground state occupation and second-order coherence function. We show that for a given pumping rate, the macroscopic occupation of the ground state and buildup of coherence occur at different temperatures. Moreover, the buildup of coherence strongly depends on the pumping scheme. We also investigate the condensate linewidth and show that the Schawlow--Townes law holds for BEC in 3D and does not hold for BEC in 2D.
The theoretical description of non-equilibrium Bose--Einstein condensate (BEC) is one of the main challenges in modern statistical physics and kinetics. The non-equilibrium nature of BEC makes it impossible to employ the well-established formalism of statistical mechanics. We develop a framework for the analytical description of a non-equilibrium phase transition to BEC that, in contrast to previously developed approaches, takes into account the infinite number of the states. We consider the limit of fast thermalization and obtain an analytical expression for the full density matrix of a non-equilibrium BEC which also covers the equilibrium case. For the particular cases of 2D and 3D, we investigate the non-equilibrium formation of BEC by finding the temperature dependence of the ground state occupation and second-order coherence function. We show that for a given pumping rate, the macroscopic occupation of the ground state and buildup of coherence occur at different temperatures. Moreover, the buildup of coherence strongly depends on the pumping scheme. We also investigate the condensate linewidth and show that the Schawlow--Townes law holds for BEC in 3D and does not hold for BEC in 2D.
The theoretical description of non-equilibrium Bose–Einstein condensate (BEC) is one of the main challenges in modern statistical physics and kinetics. The non-equilibrium nature of BEC makes it impossible to employ the well-established formalism of statistical mechanics. We develop a framework for the analytical description of a non-equilibrium phase transition to BEC that, in contrast to previously developed approaches, takes into account the infinite number of continuously distributed states. We consider the limit of fast thermalization and obtain an analytical expression for the full density matrix of a non-equilibrium ideal BEC which also covers the equilibrium case. For the particular cases of 2D and 3D, we investigate the non-equilibrium formation of BEC by finding the temperature dependence of the ground state occupation and second-order coherence function. We show that for a given pumping rate, the macroscopic occupation of the ground state and buildup of coherence may occur at different temperatures. Moreover, the buildup of coherence strongly depends on the pumping scheme. We also investigate the condensate linewidth and show that the Schawlow–Townes law holds for BEC in 3D and does not hold for BEC in 2D.
When a gas of bosonic atoms is cooled sufficiently near to absolute zero, it undergoes a startling transition from a collection of randomly-behaved individual atoms to a phase in which thousands, even millions, of atoms act in collective unison. In this new phase, known as a Bose-Einstein condensate (BEC), the atoms all occupy the exact same quantum mechanical ground state and can be collectively described by a single phase-coherent wavefunction of macroscopic size. But how exactly does this harmonious behaviour arise within a collection of incoherent atoms as their temperature decreases? Does the condensate wavefunction grow from a single central nucleus, or does it emerge simultaneously throughout the system? And how does the situation differ when the gas is cooled suddenly, away from thermodynamic equilibrium? Can non-equilibrium behaviour be inherited by the condensate that forms? In this thesis we aim to answer these and other questions by studying the formation dynamics of Bose-Einstein condensates and the development of phase coherence. The research described herein makes use of several theoretical tools, which we lay out on the workbench in Chapter 2: Semi-classical Hartree-Fock theory provides a mean-field description of a condensate and thermal cloud at equilibrium; Quantum kinetic theory is used to calculate the dynamics of condensate growth into an equilibrium state, from a non-equilibrium thermal cloud; The stochastic projected Gross-Pitaevskii equation (SPGPE) employs a classical-field description to encapsulate the non-equilibrium dynamics of the condensate and low-energy region of the atom cloud, coupled to an equilibrium thermal cloud. This thesis is composed of four principal studies, each constituting a separate chapter.In Chapter 3 we study the formation of a Bose-Einstein condensate in a cigar-shaped 3D harmonic trap, induced by the controlled addition of an attractive ``dimple'' potential along the weak axis in a recent experiment [Garrett et al., Phys. Rev. A 83, 013630 (2011)]. In this manner condensation is induced without cooling due to a localized increase in the phase space density. We perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self-consistent semi-classical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time dependence of the condensate fraction with calculations based on quantum kinetic theory. In Chapter 4 we study the onset of phase coherence during the formation of a BEC following a sudden quench to below the transition temperature. Using the SPGPE, we simulate a recent experiment performed at ANU in which a sudden RF quench is applied to a dilute gas of metastable Helium atoms confined in an elongated 3D harmonic magnetic trap. At regularly-spaced intervals following the quench, broadband RF pulses are applied to homogeneously outcouple small (<5%) pulses of atoms, which expand and fall to a microchannel plate detector. The second-order correlation function at the detector, postulated to provide a measure of phase coherence of the in-situ cloud from which they were outcoupled, is measured versus time and compared with the condensate growth curve. A delay was observed in the relaxation of the second-order correlation function, relative to the peak in the number of condensate atoms. In our simulations we incorporate the effects of periodic outcoupling both the condensate and thermal cloud, and monitor the time-dependence of first- and second-order correlation functions in both position and momentum space. We observe long-range phase coherence in the final equilibrium state, though mild persistent phase fluctuations are present. We also observe the curious phenomenon of ``superbunching'' - an enhancement of density fluctuations - at specific values of the momentum, which we relate to anomalous correlations of fluctuations. Despite obtaining a good quantitative fit to the experimentally-measured condensate growth curve in our SPGPE simulations, we do not observe any evidence of the delay in coherence observed in the experiment, which remains unexplained. In Chapter 5 we use the SPGPE to study the crossover from a true BEC with long-range phase coherence to a quasicondensate dominated by phase fluctuations, for the case of a Bose gas in an highly-elongated 3D harmonic trap. We use parameters corresponding to an experiment in which a sudden quench resulted in the formation of a non-equilibrium condensate, the dynamics of which we study in the chapter that follows. Sweeping through a range of temperatures and chemical potentials within the crossover region, we reveal the behaviour of phase fluctuations through analysis of correlation functions. We compare the ideal-gas approximation for the transition temperature with the more accurate value obtained via the Binder cumulant. By comparing the phase coherence length with the spatial extent across which density fluctuations are suppressed, we identify the BEC/quasicondensate crossover point. We compare directly with predictions of mean-field theory [Petrov et al., Phys. Rev. Lett. 87, 050404 (2001)] and find that phase fluctuations persist to lower temperatures in our SPGPE simulations. We calculate anomalous correlations of the fluctuation field and show that the resulting correction to the mean-field condensate self-interaction energy is significant within the quasicondensate regime. In Chapter 6 we study the formation of a BEC into a non-equilibrium state by simulating an experiment in which a dilute gas of Rubidium-87 atoms in a highly-elongated 3D harmonic magnetic trap is rapidly quenched to below the transition temperature. Following the quench, a condensate is observed to form with an excessively elongated shape, and subsequently undergo damped quadrupole shape oscillations [Shvarchuck et al., Phys. Rev. Lett. 89, 270404 (2002)]. We simulate this experiment by modifying the SPGPE to incorporate a hydrodynamic thermal cloud, parametrized by axially and temporally-dependent temperature and chemical potential. Comparison against simulations without this axial dependence suggests the hydrodynamic behaviour of the thermal cloud plays a crucial role in the formation of a non-equilibrium condensate. Analyses of the first-order correlation function and the condensate wavefunction reveal the spontaneous formation of quantum vortices, as well as strong phase fluctuations at equilibrium indicating quasicondensate behaviour.
Bose–Einstein condensation (BEC) is discussed at the level of an advanced course of statistical thermodynamics, clarifying some formal and physical aspects that are usually not covered by the standard pedagogical literature. The non-conventional approach adopted starts by showing that the continuum limit, in certain cases, cancels out the crucial role of the bosonic ground level. If so, a correct treatment of the problem, including the ground level population N 0 by construction, leads to BEC in a straightforward way. For a density of states of the form G ( )∝ γ , the chemical potential µ is explicitly calculated as a function of the temperature T and of the number N of bosons, for various significant values of the positive exponent γ. In the thermodynamic limit, in which the boson number N diverges and BEC is a sharp process, the chemical potential µ is a singular function of T at the critical temperature T B , determined by an appropriate critical exponent. The condensate population N 0 is studied analytically and numerically as a function of the temperature, for various values of N and for different γ. This provides an accurate description of the way BEC approaches the character of a sharp phase transition. Some aspects of the real experiments on BEC, involving a finite number of bosons, are also illustrated.
In this Letter, we give an analytical quantum description of a nonequilibrium polariton Bose-Einstein condensate (BEC) based on the solution of the master equation for the full polariton density matrix in the limit of fast thermalization. We find the density matrix of a nonequilibrium BEC, that takes into account quantum correlations between all polariton states. We show that the formation of BEC is accompanied by the build-up of cross-correlations between the ground state and the excited states reaching their highest values at the condensation threshold. Despite the nonequilibrium nature of polariton systems, we show the average population of polariton states exhibits the Bose-Einstein distribution with an almost zero effective chemical potential above the condensation threshold similar to an equilibrium BEC. We demonstrate that above threshold the effective temperature of polaritons drops below the reservoir temperature.
A nonequilibrium approach to the dynamics and statistics of the condensate of an ideal N-atom Bose gas cooling via interaction with a thermal reservoir using the canonical ensemble is developed. We derive simple analytical expressions for the canonical partition function and equilibrium distribution of the number of atoms in the ground state of a trap under different approximations, and compare them with exact numerical results. The N-particle constraint associated with the canonical ensemble is usually a burden. In the words of Kittel, ``in the investigation of the Bose-Einstein\dots{}laws it is very inconvenient to impose the restriction that the number of particles in the subsystem shall be held constant.'' But in the present approach, based on the analogy between a second-order phase transition and laser threshold behavior, the N-particle constraint makes the problem easier. We emphasize that the present work provides another example of a case in which equilibrium (detailed balance) solutions to nonequilibrium equations of motion provide a useful supplementary approach to conventional statistical mechanics. We also discuss some dynamical and mesoscopic aspects of Bose-Einstein condensation. The conclusion is that the present analytical (but approximate) results, based on a nonequilibrium approach, are in excellent agreement with exact (but numerical) results. The present analysis has much in common with the quantum theory of the laser.
We consider the ordering kinetics in a strongly non-equilibrium state of a (weakly) interacting Bose gas, characterized, on one hand, by large occupation numbers, and, on the other hand, by the absence of long-range order. Up to higher-order corrections in inverse occupation numbers, the evolution is described by non-linear Schrodinger equation with a turbulent initial state. The ordering process is rather rich and involves a number of qualitatively different regimes that take place in different regions of energy space. Specially addressed is the case of evolution in an external potential.