The non-equilibrium dynamics of trapped ultracold atomic gases, or mixtures thereof, is an extremely rich subject. Despite 20 years of studies, and remarkable progress mainly on the experimental front, numerous open question remain, related to the growth, relaxation and thermalisation of such systems, and there is still no universally accepted theory for their theoretical description. In this paper we discuss one of the state-of-the-art kinetic approaches, which gives an intuitive picture of the physical processes happening at the microscopic scale, being broadly applicable both below and above the critical region (but not within the critical region itself, where fluctuations become dominant and symmetry breaking takes place). Specifically, the 'Zaremba–Nikuni–Griffin' (ZNG) scheme provides a self-consistent description of the coupling between the condensate and the thermal atoms, including the collisions between these two subsystems. It has been successfully tested against experiments in various settings, including investigation of collective modes (e.g. monopole, dipole and quadrupole modes), dissipation of topological excitations (solitons and vortices) as well as surface evaporative cooling. Here, we show that the ZNG model can capture two important aspects of non-equilibrium dynamics for both single-component and two-component BECs: the Kohn mode (the undamped dipole oscillation independent of interactions and temperature) and (re)thermalisation leading to condensate growth following sudden evaporation. Our simulations, performed in a spherically symmetric trap reveal (i) an interesting two-stage dynamics and the emergence of a prominent monopole mode in the evaporative cooling of a single-component Bose gas, and (ii) the long thermalisation time associated with the sympathetic cooling of a realistic two-component mixture. Related open questions arise about the mechanisms and the nature of thermalisation in such systems, where further controlled experiments are needed for benchmarking.
The non-equilibrium dynamics of trapped ultracold atomic gases, or mixtures thereof, is an extremely rich subject. Despite 20 years of studies, and remarkable progress mainly on the experimental front, numerous open question remain, related to the growth, relaxation and thermalisation of such systems, and there is still no universally accepted theory for their theoretical description. In this paper we discuss one of the state-of-the-art kinetic approaches, which gives an intuitive picture of the physical processes happening at the microscopic scale, being broadly applicable both below and above the critical region (but not within the critical region itself, where fluctuations become dominant and symmetry breaking takes place). Specifically, the 'Zaremba–Nikuni–Griffin' (ZNG) scheme provides a self-consistent description of the coupling between the condensate and the thermal atoms, including the collisions between these two subsystems. It has been successfully tested against experiments in various settings, including investigation of collective modes (e.g. monopole, dipole and quadrupole modes), dissipation of topological excitations (solitons and vortices) as well as surface evaporative cooling. Here, we show that the ZNG model can capture two important aspects of non-equilibrium dynamics for both single-component and two-component BECs: the Kohn mode (the undamped dipole oscillation independent of interactions and temperature) and (re)thermalisation leading to condensate growth following sudden evaporation. Our simulations, performed in a spherically symmetric trap reveal (i) an interesting two-stage dynamics and the emergence of a prominent monopole mode in the evaporative cooling of a single-component Bose gas, and (ii) the long thermalisation time associated with the sympathetic cooling of a realistic two-component mixture. Related open questions arise about the mechanisms and the nature of thermalisation in such systems, where further controlled experiments are needed for benchmarking.
In this thesis we study the dynamics of dark solitons, vortices and Josephson effects, in trapped atomic Bose-Einstein condensates. We firstly demonstrate a novel, sound mediated, long-range interaction mechanism between dark solitons using the Gross-Pitaevskii equation. Then we present the effect of finite temperatures on the dynamics of vortices and Josephson effects using a generalised form of this equation, which additionally includes a source term modelling the coupling of the condensate to a fully dynamical thermal cloud, described by a quantum Boltzmann equation. The latter formalism is known as the Zaremba-Nikuni-Griffin scheme (ZNG). The sound-mediated interaction between dark solitons is such that the speed and trajectory of one soliton in a condensate confined by a harmonic trap can be significantly modified by the presence of a second soliton in the same trap. By confining the two solitons to spatially separated subregions of a double well trap, we show how this effect can be magnified. In particular we find it can be large, and therefore detectable, in experimentally relevant geometries, including high-periodicity optical lattices in which the ability of solitons to act as both absorbers and emitters of energy becomes apparent. At finite temperatures, we model the dynamics via the ZNG scheme. Firstly, we give a detailed description of the scheme for numerically solving these coupled equations including all collisional terms. Then we study the effect of a dynamical thermal cloud on vortex dynamics, in particular, focussing on the experimentally-relevant quantities of precession frequencies, decay rates and vortex core brightness. The changes in these are found to increase with increasing temperature for a pancake-shaped geometry, in particular for trapping parameters of a recent experiment of Freilich et al. (Science 329, 1182 (2010)). Particle exchanging collisions between the condensate and thermal cloud are found to be crucial in determining the rate of decay of a precessing vortex. We further show that, rotation of the thermal cloud can be a mechanism for radially translating the position of a vortex, a scenario under current investigation. Finally, we analyse the effect of finite temperatures on population dynamics in a double well atom chip experiment by LeBlanc et al. (Phys. Rev. Lett. 106, 025302 (2011)). Focussing on the coupled evolution in the absence of collisions, we find that it generates significant damping which is, however, less than that observed in the experiment. A detailed comparison in the presence of collisions is beyond the scope of this thesis.
We review a self-consistent scheme for modelling trapped weakly-interacting quantum gases at temperatures where the condensate coexists with a significant thermal cloud. This method has been applied to atomic gases by Zaremba, Nikuni, and Griffin, and is often referred to as ZNG. It describes both mean-field-dominated and hydrodynamic regimes, except at very low temperatures or in the regime of large fluctuations. Condensate dynamics are described by a dissipative Gross-Pitaevskii equation (or the corresponding quantum hydrodynamic equation with a source term), while the non-condensate evolution is represented by a quantum Boltzmann equation, which additionally includes collisional processes which transfer atoms between these two subsystems. In the mean-field-dominated regime collisions are treated perturbatively and the full distribution function is needed to describe the thermal cloud, while in the hydrodynamic regime the system is parametrised in terms of a set of local variables. Applications to finite temperature induced damping of collective modes and vortices in the mean-field-dominated regime are presented.
In this Thesis we study steady state solutions and dynamical evolutions of two– component atomic Bose–Einstein Condensates. We initially investigate the equilibrium properties of condensate mixtures in harmonic trapping potentials at zero temperature. Subsequently we simulate the coupled growth of these condensates by inclusion of damping terms. Finally, we investigate the evolution of coupled Bose gases via the so-called classical–field method. A recent experiment [D. J. McCarron et al., Phys. Rev. A, 84, 011603(R) (2011)] achieved Bose–Einstein Condensation of a two–species 87Rb–133Cs phase segregated mixture in harmonic trapping potentials. Depending on relative atom numbers of the two species, three distinct regimes of density distributions were observed. For these experimental parameters, we investigate the corresponding time–independent ground state solutions through numerical simulations of the coupled Gross–Pitaevskii equations. By including experimentally relevant shifts between the traps, we observe a range of structures including ‘ball and shell’ formations and axially/radially separated states. These are found to be very sensitive to the trap shifts. For all three experimental regimes, our numerical simulations reveal good qualitative agreement. The observed experimental profiles cannot be guaranteed to be fully equilibrated. This, coupled with the rapid sympathetic cooling of the experimental system, leads to a situation where growth may play a determining factor in the density structures formed. To investigate this further, we introduce phenomenological damping to describe the associated condensate growth/decay, revealing a range of transient structures. However, such a model always predicts the predominance of one condensate species over longer evolution times. Work undertaken by collaborators with the more elaborate Stochastic Projected Gross–Pitaevskii equations, which can describe condensate formation by coupling to a heat bath, predicts the spontaneous formation of dark–bright solitons. Motivated by this, we show how the presence of solitons can affect the condensate distribution, thus highlighting the overall dynamical role in the emerging patterns. Finally, we use classical field methods to analyse the evolution of non trapped Bose gases from strongly nonequilibrium initial distributions. The contrast between miscible (overlapping) and immiscible (phase segregated) components gives rise to important distinctions for condensate fractions and the formation of domains and vortices. In addition, splitting the particles of a single component thermalised state into two components is investigated. We then study the effects of suddenly quenching the strength of the interspecies interactions. Under suitable conditions, this quench generates isotropic vortex tangles. While this tangle subsequently decays over time, we propose how a repeat sequence of quenches at regular intervals could be employed to drive the tangle, thereby potentially providing a novel route to the generation of quantum turbulence.
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.
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.
We construct a nonequilibrium theory for the dynamics of two interacting finite-temperature atomic Bose-Einstein condensates and use it to numerically estimate the relative rates of the arising collisional processes near equilbrium. The condensates are described by dissipative Gross-Pitaevskii equations, coupled to quantum Boltzmann equations for the thermal atoms. The density-density interactions between atoms in different components facilitate a number of transport processes of relevance to sympathetic cooling: in particular, considering realistic miscible and immiscible trapped atomic $^{87}\mathrm{Rb}\text{\ensuremath{-}}^{41}\mathrm{K}$ and $^{87}\mathrm{Rb}\text{\ensuremath{-}}^{85}\mathrm{Rb}$ condensate mixtures, we highlight the dominance of an intercomponent scattering process associated with collisional ``exchange'' of condensed and thermal atoms between the components close to equilibrium.