We simulate the growth of a dual species Bose-Einstein condensate using a Gross-Pitaevskii equation with an additional gain term giving rise to the growth. Such growth occurs during simultaneous evaporative cooling of a mixture of two gases. The ground state of a dual condensate is normally either a miscible mixture, or an immiscible phase with two spatially separated components. In a cigar trap the ground state typically consists of one component in the center, and the other component flanking it. Our simulations show that when the condensates are formed in a cigar trap and the mixture is phase separated, then the final state upon the end of the growth is generally far from the true ground state of the system. Instead it consists of multiple, interleaved bubbles of the two species. Such a pattern was observed recently in an experiment by Wieman's group at JILA [Papp, Pino, and Wieman, Phys. Rev. Lett. 101, 040402 (2008)], and our simulations are in good qualitative agreement with the experiment. We explain the pattern formation as due to the onset of modulation instability during growth, and study the dependence of the final state pattern on various parameters of the system.
We simulate the growth of a dual species Bose-Einstein condensate using a Gross-Pitaevskii equation with an additional gain term giving rise to the growth. Such growth occurs during simultaneous evaporative cooling of a mixture of two gases. The ground state of a dual condensate is normally either a miscible mixture, or an immiscible phase with two spatially separated components. In a cigar trap the ground state typically consists of one component in the center, and the other component flanking it. Our simulations show that when the condensates are formed in a cigar trap and the mixture is phase separated, then the final state upon the end of the growth is generally far from the true ground state of the system. Instead it consists of multiple, interleaved bubbles of the two species. Such a pattern was observed recently in an experiment by Wieman's group at JILA [Papp, Pino, and Wieman, Phys. Rev. Lett. 101, 040402 (2008)], and our simulations are in good qualitative agreement with the experiment. We explain the pattern formation as due to the onset of modulation instability during growth, and study the dependence of the final state pattern on various parameters of the system.
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 introduce a phenomenological mean-field model to describe the growth of immiscible two-species atomic Bose-Einstein condensates towards some equilibrium. Our model is based on the coupled Gross-Pitaevskii equations with the addition of dissipative terms to account for growth. While our model may be applied generally, we take a recent Rb-Cs experiment [McCarron et al., Phys. Rev. A 84 011603(R) (2011)] as a case study. As the condensates grow, they can pass through ranging transient density structures which can be distinct from the equilibrium states, although such a model always predicts the predominance of one condensate species over longer evolution times.
The dynamics of multiple domain formation caused by the modulation instability of two-component Bose-Einstein condensates in an axially symmetric trap are studied by numerically integrating the coupled Gross-Pitaevskii equations. The modulation instability induced by the intercomponent mean-field coupling occurs in the out-of-phase fluctuation of the wave function and leads to the formation of multiple domains that alternate from one domain to another, where the phase of one component jumps across the density dips where the domains of the other exist. This behavior is analogous to a soliton train, which explains the origin of the long lifetime of the spin domains observed by Miesner et al. [Phys. Rev. Lett. 82, 2228 (1999)]].
We investigate nonlinear dynamics induced by the modulation instability of a two-component mixture in an atomic Bose-Einstein condensate. The nonlinear dynamics is examined using numerical simulations of the time-dependent coupled Gross-Pitaevskii equations. The unstable modulation grows from initially miscible condensates into various types of vector solitary waves, depending on the combinations of the sign of the coupling constants (intracomponent and intercomponent). We discuss the detailed features of the modulation instability, dynamics of solitary wave formation, and an analogy with the collapsing dynamics in a single-component condensate with attractive interactions.
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 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.