Quench dynamics and emergence of phase separation in two-component atomic Bose gases at zero temperature and above the Bose-Einstein-condensation critical temperature
We study the dynamics of two-component atomic Bose gases initially in a mixture encountering a sudden quench of the interspecies interactions. The dynamics above the critical temperature ${T}_{c}$ is studied using a leading-order large-$N$ approximation that predicted a phase transition from mixing to phase separation as a function of the interspecies coupling. Here, we explore the dynamics of this phase transition following a quench and compare our results to those found at zero temperature using the time-dependent Gross-Pitaevskii equations which ignore quantum and thermal fluctuations. In the regime above ${T}_{c}$ where no condensate is present, however, the time evolution of the densities following the quench exhibits features similar to those found at zero temperature where only the condensates contribute to the densities. When the interspecies interaction jumps above the critical value, we observe dynamical transitions from a homogeneous mixture to a phase-separated structure for both cases. Our simulations suggest that at temperatures above ${T}_{c}$ where no condensate is present this dynamical transition should still be observable in experiments.
We study the dynamics of two-component atomic Bose gases initially in a mixture encountering a sudden quench of the interspecies interactions. The dynamics above the critical temperature Tc is studied using a leading-order large-N approximation that predicted a phase transition from mixing to phase separation as a function of the interspecies coupling. Here, we explore the dynamics of this phase transition following a quench and compare our results to those found at zero temperature using the time-dependent Gross-Pitaevskii equations which ignore quantum and thermal fluctuations. In the regime above Tc where no condensate is present, however, the time evolution of the densities following the quench exhibits features similar to those found at zero temperature where only the condensates contribute to the densities. When the interspecies interaction jumps above the critical value, we observe dynamical transitions from a homogeneous mixture to a phase-separated structure for both cases. Our simulations suggest that at temperatures above Tc where no condensate is present this dynamical transition should still be observable in 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 discuss the mean-field theories obtained from the leading order in a large-N approximation for one- and two-component dilute Bose gases. For a one-component Bose gas, this approximation has the following properties: the Bose-Einstein condensation (BEC) phase transition is second order but the critical temperature Tc is not shifted from the noninteracting gas value T0. The spectrum of excitations in the BEC phase resembles the Bogoliubov dispersion with the usual coupling constant replaced by the running coupling constant which depends on both temperature and momentum. We then study two-component Bose gases with both inter- and intraspecies interactions and focus on the stability of the mixture state above Tc. Our mean-field approximation predicts an instability from the mixture state to a phase-separated state when the ratio of the interspecies interaction strength to the intraspecies interaction strength (assuming equal strength for both species) exceeds a critical value. At high temperature this is a structural transition and the global translational symmetry is broken. Our work complements previous studies on the instability of the mixture phase in the presence of BEC.
The miscibility of two interacting quantum systems is an important testing ground for the understanding of complex quantum systems. Two-component Bose-Einstein condensates enable the investigation of this scenario in a particularly well controlled setting. In a homogeneous system, the transition between mixed and separated phases is fully characterized by a miscibility parameter based on the ratio of intra- to interspecies interaction strengths. Here we show, however, that this parameter is no longer the optimal one for trapped gases, for which the location of the phase boundary depends critically on atom numbers. We demonstrate how monitoring of damping rates and frequencies of dipole oscillations enables the experimental mapping of the phase diagram by numerical implementation of a fully self-consistent finite-temperature kinetic theory for binary condensates. The change in damping rate is explained in terms of surface oscillation in the immiscible regime, and counterflow instability in the miscible regime, with collisions becoming only important in the long time evolution.
We examine the dynamics associated with the miscibility-immiscibility transition of trapped two-species Bose-Einstein condensates (TBECs) of dilute atomic gases in the presence of vortices. In particular, we consider TBECs of Rb hyperfine states and Rb-Cs mixture. In the case of a singly charged vortex in only one of the condensates, there is enhancement when the vortex is present in the species which occupies the edges at phase separation. But suppression occurs when the vortex is in the species which occupies the core region. To examine the role of the vortex, we quench the interspecies interactions to drive the TBECs from miscible to immiscible phase, and use the time-dependent Gross-Pitaevskii equation to probe the phenomenon of phase separation. We also examine the effects of higher charged vortex to phase separation.
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.
The miscibility condition for a binary mixture of two interacting Bose-Einstein condensates is shown to be deeply affected by interaction driven thermal fluctuations. These give rise to a first order phase transition to a demixed phase with full spatial separation of the two condensates, even if the mixture is miscible at zero temperature. Explicit predictions for the isothermal compressibility, the spin susceptibility, and the phase transition temperature T_{M} are obtained in the framework of Popov theory, which properly includes beyond mean-field quantum and thermal fluctuations in both the spin and density channels. For a mixture of two sodium condensates occupying the hyperfine states |F=1⟩,|m_{F}=1⟩ and |F=1,m_{F}=-1⟩, respectively, T_{M} is predicted to occur at about 0.7 times the usual BEC critical temperature.