In this paper, we present a method for simulating the kinetic evolution of a dilute gas of atoms that are cooled below the critical temperature for Bose-Einstein condensation. Our method gives insight into the formulation of physical kinetics by illustrating directly the decomposition of the distribution function into an infinite sum of single-particle trajectories. This approach is valid for the entire range of phase-space densities, although we limit the discussion here to exclude the region where the condensate fraction is close to unity and the effect of the mean field is significant. We present explicit calculations of finite number effects on equilibrium, the dynamic build-up of the ground state, and simulations of evaporative cooling.
In this paper, we present a method for simulating the kinetic evolution of a dilute gas of atoms that are cooled below the critical temperature for Bose-Einstein condensation. Our method gives insight into the formulation of physical kinetics by illustrating directly the decomposition of the distribution function into an infinite sum of single-particle trajectories. This approach is valid for the entire range of phase-space densities, although we limit the discussion here to exclude the region where the condensate fraction is close to unity and the effect of the mean field is significant. We present explicit calculations of finite number effects on equilibrium, the dynamic build-up of the ground state, and simulations of evaporative cooling.
We have simulated the evaporative cooling of a dilute gas of Bose particles including quantum statistics using a Monte Carlo method. This approach can model situations which are far away from quasiequilibrium such as occur during forced evaporative cooling. We have also simulated the dynamical formation of Bose-Einstein condensate for homogeneous and inhomogeneous Bose gases under the random-phase approximation. It was found that the rate of accumulation of Bose particles into low-energy states through collisions is increased by forced evaporation; and a macroscopic population at ground state can be reached at a time scale characterized by classical collision time for an inhomogeneous gas in a harmonic potential. We present the results of simulations for the evaporative cooling and formation of a Bose-Einstein condensate in one-, two-, and three-dimensional position cuts and energy cuts.
The problem of understanding how a coherent, macroscopic Bose- Einstein condensate (BEC) emerges from the cooling of a thermal Bose gas has attracted significant theoretical and experimental interest over several decades. The pioneering achievement of BEC in weakly interacting dilute atomic gases in 1995 was followed by a number of experimental studies examining the growth of the BEC number, as well as the development of its coherence. More recently, there has been interest in connecting such experiments to universal aspects of nonequilibrium phase transitions, in terms of both static and dynamical critical exponents. Here, the spontaneous formation of topological structures such as vortices and solitons in quenched cold-atom experiments has enabled the verification of the Kibble-Zurek mechanism predicting the density of topological defects in continuous phase transitions, first proposed in the context of the evolution of the early universe. This chapter reviews progress in the understanding of BEC formation and discusses open questions and future research directions in the dynamics of phase transitions in quantum gases. Introduction The equilibrium phase diagram of the dilute Bose gas exhibits a continuous phase transition between condensed and noncondensed phases. The order parameter characteristic of the condensed phase vanishes above some critical temperature T c and grows continuously with decreasing temperature below this critical point. However, the dynamical process of condensate formation has proved to be a challenging phenomenon to address both theoretically and experimentally. This formation process is a crucial aspect of Bose systems and of direct relevance to all condensates discussed in this book, despite their evident system-specific properties. Important questions leading to intense discussions in the early literature include the time scale for condensate formation and the role of inhomogeneities and finite-size effects in “closed” systems. These issues are related to the concept of spontaneous symmetry breaking, its causes, and implications for physical systems (see, for example, Chapter 5 by Snoke and Daley). In this chapter, we give an overview of the dynamics of condensate formation and describe the present understanding provided by increasingly well-controlled cold-atom experiments and corresponding theoretical advances over the past twenty years. We focus on the growth of BECs in cooled Bose gases, which, from a theoretical standpoint, requires a suitable nonequilibrium formalism.
Evaporative cooling of trapped atoms is studied based on the quantum kinetic theory. We extend the classical analysis of Luiten et al. [Phys. Rev. A 53, 381 (1996)] to treat Bose statistics of alkali-metal atoms. Dynamics in evaporative cooling process is described by a kinetic equation for the truncated Bose-Einstein distribution function under the assumption that the system is close to the thermal equilibrium. This approach is applicable to the ``slow'' evaporative cooling normally adopted in most experiments for the efficient production of Bose-Einstein condensates. Time-evolution calculations explain the current experiment well with sodium atoms using a slow evaporative cooling and demonstrate a rapid growth of condensate in the cooling process.
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.
Using a functional formulation of the Keldysh theory we investigate the time evolution of a dilute Bose gas through the critical point. We determine the critical region of the gas and show that if the system is quenched inside this region, Bose-Einstein condensation is nucleated on a time scale of O(\ensuremath{\Elzxh}/kBTc) by means of a coherent population of the one-particle ground state. We also examine the subsequent buildup of the condensate density, taking place on a much longer time scale. However, in the experimentally interesting case of spin-polarized atomic hydrogen or cesium it is still short compared to the lifetime of the gas.
Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in L2-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. Full-text embargoed until: 2026-01-01