The microscopic theory of the time evolution of an interacting Bose gas being cooled to the temperature T<T c is presented. The considerations are based on the exact (in classical limit) equations of motion derived by Langer. It is shown that at the initial stages of relaxation the sharp peak appears at the low energy tail of the distribution function. The following time evolution of this peak leads to a delta function formation. When this peak is too narrow the Boltzmann equation is no longer applicable and one must apply the exact equations of motion or the quasi-Boltzmann equation derived in this work which takes into account the coherent effects. The evolution of the distribution function and of the energy spectrum in the system is given by the solution of the equation of motion. The development of a new phase is of explosive character and has a lot of features in common with instabilities in non-linear plasmas, optics, etc., and thus is of quite a general nature. In addition it is shown that the phenomenological models widely used for the description of the dynamics of the lambda -transition in a critical region can be derived from microscopic equations of motion.
The microscopic theory of the time evolution of an interacting Bose gas being cooled to the temperature T<T c is presented. The considerations are based on the exact (in classical limit) equations of motion derived by Langer. It is shown that at the initial stages of relaxation the sharp peak appears at the low energy tail of the distribution function. The following time evolution of this peak leads to a delta function formation. When this peak is too narrow the Boltzmann equation is no longer applicable and one must apply the exact equations of motion or the quasi-Boltzmann equation derived in this work which takes into account the coherent effects. The evolution of the distribution function and of the energy spectrum in the system is given by the solution of the equation of motion. The development of a new phase is of explosive character and has a lot of features in common with instabilities in non-linear plasmas, optics, etc., and thus is of quite a general nature. In addition it is shown that the phenomenological models widely used for the description of the dynamics of the lambda -transition in a critical region can be derived from microscopic equations of motion.
The kinetics of the formation of coherent correlation properties associated with Bose condensation is studied in detail. The evolution of a nonequilibrium state with no “condensate-seed” is related to a hierarchy of relaxation times. At the first stage, a particle flux in energy space toward low energies sets in. The evolution in this case is described by a nonlinear Boltzmann equation, with a characteristic time given by interparticle collisions. When the particles which will later form the condensate have a kinetic energy which is less than the potential energy, a quasicondensate starts to form. In this stage, fluctuations of the density (but not of the phase) are suppressed and short-range coherent correlation properties are governed by the equation of motion for a quasiclassical complex field. The next stage is connected with the formation of the long-range order. The time for forming topological order and therefore genuine superfluidity proves to be dependent on the system size. The off-diagonal long-range order, arising after the attenuation of long-wave phase fluctuations, has a size-dependent relaxation time as well. Introduction The problem of Bose–Einstein Condensation (BEC) kinetics, being interesting in itself, has acquired a special significance in connection with experimental efforts to observe this condensation in a number of systems with particles with a finite lifetime. Such systems include spin-polarized atomic hydrogen [1], excitons [2] and biexcitons [3] in semiconductors and, more recently, laser-cooled atomic systems [4].
Kinetics of Bose condensation of an ideal Bose gas and of the classical Bose gas of hard spheres is considered. Deviations from the equilibrium state are not assumed to be small. An analytic solution for the problem of an ideal Bose gas is found and the role of collisions between the Bose particles themselves is discussed qualitatively. It is shown that Bose condensation is a slowing down process which can be completed at $t\ensuremath{\rightarrow}\ensuremath{\infty}$ while the formation of the high-energy tail of the distribution function is relatively fast. The transition time is finite if the nuclei of new phase are present in the beginning of 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.
Summary form only given. The dilute gas experiments on the Bose-Einstein condensation of magnetically confined alkali vapors have already generated widespread interest across a broad spectrum of traditional fields of physics. They have provided an exciting stage on which to build bridges linking the theory of complex and collective phenomena in superconducting and superfluid systems with the single particle microscopic physics described in quantum optics and laser physics. The approach I discuss is based on a derivation of the quantum kinetic theory from first principles. Based on a microscopic picture of the weakly interacting bosonic gas, a generalized kinetic theory for a Markovian many-body density operator is derived. The density operator is functionally dependent on a few key variables, which must be carefully chosen. These quantities serve as master variables and determine the system's evolution on a coarse-grained time scale. The weak interactions allow a perturbative expansion of the evolution, from which kinetic equations that describe the dynamics of the expectation value of any single-time observable are obtained. The final structure of the theory has many analogies with the quantum theory of the laser, such as amplitude and phase noise for the macroscopic field, squeezing dynamics, and a nonlinear master equation describing the underlying evolution.
This thesis manuscript is devoted to different aspects of the non-equilibrium dynamics of a moving two-dimensional ultra-cold Bose gas, under the simultaneous action of a spatially disordered potential and interactions between particles. In the absence of interactions, the momentum distribution of the gas is the sum of two well-known components: a classical, diffusive background made of particles multiply scattered by the disorder, and a coherent backscattering interference peak (CBS peak), consequence of the weak localization phenomenon. In the presence of weak interactions, mean field simulations show that these two contributions evolve in time. This evolution is described by two coupled kinetic equations, which we derive using a microscopic quantum transport theory. These predict the destruction of the CBS peak by an effective dephasing mechanism, and the thermalization of the diffusive background. As the interaction strength increases, the effects of disorder become negligible, and the dynamics of the gas resembles that of a superfluid. In this regime, the gas quickly enters a so-called "pre-thermal" state: the spatial correlations g_1 decrease algebraically and propagate along a light cone, in accordance with theoretical predictions. Finally, at long times, our simulations describe the global evolution of the system towards a thermal equilibrium state, characterized by a Kosterlitz-Thouless transition separating a normal fluid state from a superfluid state. Some properties of this phase transition, resulting from the non-equilibrium evolution in the presence of disorder, are discussed.
The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.