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.
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.
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 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.
In this Thesis is presented a study on dynamical properties of mixtures of ultraold Bose gases. The behaviour of this system in different regimes is analysed: with and without coherent coupling between the two components, in homogeneous and harmonic shaped trapping potentials and in different dimensions and geometries. Most of the results presented here have been obtained by means of numerical solutions of coupled Gross-Pitaevskii equations and have been compared with theoretical predictions (and sometimes experiments), describing the same phenomena. In particualr the stability of persistent currents in a two-component Bose-Einstein condensate in a toroidal trap is studied in both the miscible and the immiscible regime. In the miscible regime we show that superflow decay is related to linear instabilities of the spin-density Bogoliubov mode. We find a region of partial stability, where the flow is stable in the majority component while it decays in the minority component. We also characterize the dynamical instability appearing for a large relative velocity between the two components. In the immiscible regime the stability criterion is modified and depends on the specific density distribution of the two components. The effect of a coherent coupling between the two components is also discussed. A study on the collective modes of the minority component of a highly unbalanced Bose-Bose mixture is also presented. In the immiscible case we find that the ground state can be a two-domain walls soliton. Although the mode frequencies are continuous at the transition, their behaviour is very different with respect to the miscible case. The dynamical behaviour of the solitonic structure and the frequency dependence on the inter- and intra-species interaction is numerically studied using coupled Gross-Pitaevskii equations. The results of the study on the static and the dynamic response of coherently coupled two component Bose-Einstein condensates due to a spin-dipole perturbation is also sown. The static dipole susceptibility is determined and is shown to be a key quantity to identify the second order ferromagnetic transition occurring at large inter-species interactions. The dynamics, which is obtained by quenching the spin-dipole perturbation, is very much affected by the system being paramagnetic or ferromagnetic and by the correlation between the motional and the internal degrees of freedom. In the paramagnetic phase the gas exhibits well defined out-of-phase dipole oscillations, whose frequency can be related to the susceptibility of the system using a sum rule approach. In particular in the interaction SU (2) symmetric case, when all the two-body interactions are the same, the external dipole oscillation coincides with the internal Rabi flipping frequency. In the ferromagnetic case, where linear response theory is not applicable, the system shows highly non-linear dynamics. In particular we observe phenomena related to ground state selection: the gas, initially trapped in a domain wall configuration, reaches a final state corresponding to the magnetic ground state plus small density ripples. Interestingly, the time during which the gas is unable to escape from its initial configuration is found to be proportional to the square root of the wall surface \ntension.
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 study the interactions between two atomic species in a binary Bose-Einstein condensate to revisit the conditions for miscibility, oscillatory dynamics between the species, steady-state solutions, and their stability. By employing a variational approach for a quasi-one-dimensional, two-atomic species condensate, we obtain equations of motion for the parameters of each species: amplitude, width, position, and phase. A further simplification leads to a reduction of the dynamics into a simple classical Newtonian system where components oscillate in an effective potential with a frequency that depends on the harmonic trap strength and the interspecies coupling parameter. We develop explicit conditions for miscibility that can be interpreted as a phase diagram that depends on the harmonic trap's strength and the interspecies coupling parameter. We numerically illustrate the bifurcation scenario whereby nontopological, phase-separated states of increasing complexity emerge out of a symmetric state as the interspecies coupling is increased. The symmetry-breaking dynamical evolution of some of these states is numerically monitored and the associated asymmetric states are also explored.
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.