The study of nonequilibrium dynamics in Bose-Einstein condensed gases has been dominated by the zero-temperature mean field Gross-Pitaevskii formalism. Motivated by recent experiments on the reflection of condensates from silicon surfaces, we revisit the so called classical field description of condensate dynamics, which incorporates the effects of quantum noise and can also be generalized to include thermal effects. The noise is included in a stochastic manner through the initial conditions. We show that the inclusion of such noise is important in the quantitative description of the recent reflection experiments.
The study of nonequilibrium dynamics in Bose-Einstein condensed gases has been dominated by the zero-temperature mean field Gross-Pitaevskii formalism. Motivated by recent experiments on the reflection of condensates from silicon surfaces, we revisit the so called classical field description of condensate dynamics, which incorporates the effects of quantum noise and can also be generalized to include thermal effects. The noise is included in a stochastic manner through the initial conditions. We show that the inclusion of such noise is important in the quantitative description of the recent reflection experiments.
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.
We present a theory for the description of energy relaxation in a nonequilibrium condensate of bosonic particles. The approach is based on coupling to a thermal bath of other particles (e.g., phonons in a crystal, or noncondensed atoms in a cold atom system), which are treated with a Monte Carlo type approach. Together with a full account of particle-particle interactions, dynamic driving, and particle loss, this offers a complete description of recent experiments in which Bose-Einstein condensates are seen to relax their energy as they propagate in real space and time. As an example, we apply the theory to the solid-state system of microcavity exciton polaritons, in which nonequilibrium effects are particularly prominent.
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 study the low-energy collective oscillations of a dilute Bose gas at finite temperature in the collisionless regime. By using a time-dependent mean-field scheme we derive for the dynamics of the condensate and noncondensate components a set of coupled equations, which we solve perturbatively to second order in the interaction coupling constant. This approach is equivalent to the finite-temperature extension of the Beliaev approximation and includes corrections to the Gross-Pitaevskii theory due to both quantum and thermal fluctuations. For a homogeneous system we explicitly calculate the temperature dependence of the velocity of propagation and damping rate of zero sound. In the case of harmonically trapped systems in the thermodynamic limit, we calculate, as a function of temperature, the frequency shift of the low-energy compressional and surface modes.
We review the stochastic Gross-Pitaevskii approach for non-equilibrium finite temperature Bose gases, focussing on the formulation of Stoof; this method provides a unified description of condensed and thermal atoms, and can thus describe the physics of the critical fluctuation regime. We discuss simplifications of the full theory, which facilitate straightforward numerical implementation, and how the results of such stochastic simulations can be interpreted, including the procedure for extracting phase-coherent (`condensate') and density-coherent (`quasi-condensate') fractions. The power of this methodology is demonstrated by successful ab initio modelling of several recent atom chip experiments, with the important information contained in each individual realisation highlighted by analysing dark soliton decay within a phase-fluctuating condensate.
We study the time-evolution of initially trapped Bose-Einstein condensates in the Gross-Pitaevskii regime. Under a physically motivated assumption on the energy of the initial data, we show that condensation is preserved by the many-body evolution and that the dynamics of the condensate wave function can be described by the time-dependent Gross-Pitaevskii equation. With respect to previous works, we provide optimal bounds on the rate of condensation (i.e. on the number of excitations of the Bose-Einstein condensate). To reach this goal, we combine the method of \cite{LNS}, where fluctuations around the Hartree dynamics for N-particle initial data in the mean-field regime have been analyzed, with ideas from \cite{BDS}, where the evolution of Fock-space initial data in the Gross-Pitaevskii regime has been considered.