We explore the time evolution of two-component Bose-Einstein condensates (BECs), quasi one dimensional with respect to their spinor dynamics, following a quench from one component BECs with a $\mathrm{U}(1)$ order parameter into two-component condensates with a $\mathrm{U}(1)\ifmmode\times\else\texttimes\fi{}{Z}_{2}$ order parameter. In our case, these two spin components have a propensity to phase separate, i.e., they are immiscible. Remarkably, these spin degrees of freedom can equivalently be described as a single-component attractive BEC. A spatially uniform mixture of these spins is dynamically unstable, rapidly amplifying any quantum or pre-existing classical spin fluctuations. This coherent growth process drives the formation of numerous spin-polarized domains, which are far from the system's ground state. At much longer times these domains grow in size, coarsening, as the system approaches equilibrium. The experimentally observed time evolution is consistent with our fully three-dimensional stochastic-projected Gross--Pitaevskii calculation.
We explore the time evolution of two-component Bose-Einstein condensates (BECs), quasi one dimensional with respect to their spinor dynamics, following a quench from one component BECs with a U(1) order parameter into two-component condensates with a U(1)\ifmmode×\else\texttimes\fiZ2 order parameter. In our case, these two spin components have a propensity to phase separate, i.e., they are immiscible. Remarkably, these spin degrees of freedom can equivalently be described as a single-component attractive BEC. A spatially uniform mixture of these spins is dynamically unstable, rapidly amplifying any quantum or pre-existing classical spin fluctuations. This coherent growth process drives the formation of numerous spin-polarized domains, which are far from the system's ground state. At much longer times these domains grow in size, coarsening, as the system approaches equilibrium. The experimentally observed time evolution is consistent with our fully three-dimensional stochastic-projected Gross--Pitaevskii calculation.
The study of quantum degenerate gases has many applications in topics such as condensed matter dynamics, precision measurements and quantum phase transitions. We built an apparatus to create 87Rb Bose-Einstein condensates (BECs) and generated, via optical and magnetic interactions, novel quantum systems in which we studied the contained phase transitions. For our first experiment we quenched multi-spin component BECs from a miscible to dynamically unstable immiscible state. The transition rapidly drives any spin fluctuations with a coherent growth process driving the formation of numerous spin polarized domains. At much longer times these domains coarsen as the system approaches equilibrium. For our second experiment we explored the magnetic phases present in a spin-1 spin-orbit coupled BEC and the contained quantum phase transitions. We observed ferromagnetic and unpolarized phases which are stabilized by the spin-orbit coupling’s explicit locking between spin and motion. These two phases are separated by a critical curve containing both first-order and second-order transitions joined at a critical point. The narrow first-order transition gives rise to long-lived metastable states. For our third experiment we prepared independent BECs in a double-well potential, with an artificial magnetic field between the BECs. We transitioned to a single BEC by lowering the barrier while expanding the region of artificial field to cover the resulting single BEC. We compared the vortex distribution nucleated via conventional dynamics to those produced by our procedure, showing our dynamical process populates vortices much more rapidly and in larger number than conventional nucleation.
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 study the dynamics of domain formation and coarsening in a binary Bose-Einstein condensate that is quenched across a miscible-immiscible phase transition. The late-time evolution of the system is universal and governed by scaling laws for the correlation functions. We numerically determine the scaling forms and extract the critical exponents that describe the growth rate of domain size and autocorrelations. Our data are consistent with inviscid hydrodynamic domain growth, which is governed by a universal dynamical critical exponent of 1/z=0.68(2). In addition, we analyze the effect of domain wall configurations which introduce a nonanalytic term in the short-distance structure of the pair correlation function, leading to a high-momentum "Porod" tail in the static structure factor, which can be measured experimentally.
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 explore the defect formation and universally critical dynamics in two-dimensional (2D) two-component Bose-Einstein condensates(BECs) subjected to two types of potential traps: a homogeneous trap and a harmonic trap.We focus on the non-equilibrium universal dynamics of the miscible-immiscible phase transition with both linear and nonlinear quenching types.Although there exists spatial independence of the critical point, we find that the inhomogeneity of trap doesn't affect the phase transition of system and the critical exponents can still be explained by the homogeneous Kibble-Zurek mechanism. By analyzing the Bogoliubov excitations, we establish a power-law relationship between the domain correlation length, the phase transition delay, and the quench time.Furthermore, through real-time simulations of phase transition dynamics, the formation of domain defects and the delay of phase transition in non-equilibrium dynamics are demonstrated, along with the corresponding universal scaling of correlation length and phase transition delay for various quench time and quench coefficients, which align well with our analytical predictions.Our study confirms that the universality class of two-component BECs remains unaffected by dimensionality, while the larger nonlinear coefficients effectively suppress non-adiabatic excitations, offering a novel perspective for addressing adiabatic evolution.
The spatial structure of the spinor Bose-Einstein condensates with the spin degrees of freedom is analyzed based on the generalized Gross-Pitaevskii (GP) equation in light of the present spin-domain experiment on ${m}_{F}=\ifmmode\pm\else\textpm\fi{}1,$ and 0 of the hyperfine state F=1 of 23Na atom gases. The GP solutions in three- and one-spatial dimensional cases reproduce the observed spin-domain structures, revealing the length scale associated with the existence of the weak interaction of the spin-spin channel, other than the ordinary coherence length related to the density-density channel. The obtained domain structure in GP is compared with the result in the Thomas-Fermi approximation. The former solution is found to better describe the observed features than the latter.