In this thesis we study the nonlinear dynamics of spin-1 and spin-2 Bose-Einstein condensates, with particular application to antiferromagnetic systems exhibiting nematic (beyond magnetic) order. Firstly, we give a derivation of the spinor energy functionals with a focus on the connections between the nonlinear terms. We derive a hierarchy of nonlinear irreducible multipole observables sensitive to different levels of nematic order, and explore the various nematic states in terms of their multipolar order, representations of their symmetries, and topological defects.\n\nWe then develop an exact solution to the nonlinear dynamics of spinor Bose-Einstein condensates. We use this solution to construct efficient and accurate numerical algorithms to evolve the spinor Gross-Pitaevskii equation in time. We demonstrate the advantages of our algorithms with several 1D numerical test problems, comparing with existing methods in the literature.\n\nWe apply our numerical methods to simulating quenches of the condensate between various antiferromagnetic phases for spin-1 and spin-2. For spin-1, we carry out quenches for a theoretical uniform system in 2D, and then specialize to the parameters used in a recent harmonically trapped experiment in 3D. We connect the long-time coarsening growth law of the relevant order parameter to the decay of half-quantum vortices, which are the relevant topological defects of the ground state.\n\nFor the spin-2 system, we investigate a novel quench from two different quadrupolar-nematic phases to an octupolar-nematic “cyclic” phase which supports 1/3 fractional vortices. We develop appropriate order parameter observables which couple to the spin and superfluid currents generated by these defects, and show that a new growth law appears with exponent 1/3.
In this thesis we study the nonlinear dynamics of spin-1 and spin-2 Bose-Einstein condensates, with particular application to antiferromagnetic systems exhibiting nematic (beyond magnetic) order. Firstly, we give a derivation of the spinor energy functionals with a focus on the connections between the nonlinear terms. We derive a hierarchy of nonlinear irreducible multipole observables sensitive to different levels of nematic order, and explore the various nematic states in terms of their multipolar order, representations of their symmetries, and topological defects.\n\nWe then develop an exact solution to the nonlinear dynamics of spinor Bose-Einstein condensates. We use this solution to construct efficient and accurate numerical algorithms to evolve the spinor Gross-Pitaevskii equation in time. We demonstrate the advantages of our algorithms with several 1D numerical test problems, comparing with existing methods in the literature.\n\nWe apply our numerical methods to simulating quenches of the condensate between various antiferromagnetic phases for spin-1 and spin-2. For spin-1, we carry out quenches for a theoretical uniform system in 2D, and then specialize to the parameters used in a recent harmonically trapped experiment in 3D. We connect the long-time coarsening growth law of the relevant order parameter to the decay of half-quantum vortices, which are the relevant topological defects of the ground state.\n\nFor the spin-2 system, we investigate a novel quench from two different quadrupolar-nematic phases to an octupolar-nematic “cyclic” phase which supports 1/3 fractional vortices. We develop appropriate order parameter observables which couple to the spin and superfluid currents generated by these defects, and show that a new growth law appears with exponent 1/3.
Spinor Bose-Einstein condensates exhibit both superfluid and magnetic order, and accommodate phases with rich symmetry properties and topological defects. Transitions between these phases can be induced by tuning external fields. In this thesis we explore the dynamics of order formation in a quasi-2D spin-1 ferromagnetic condensate following a quench from an unmagnetised phase to one of three ferromagnetic phases. The ferromagnetic phases exhibit distinct symmetry properties (easy-plane, easy-axis or isotropic) and support distinct topological defects. In each phase we observe scale invariant ordering and identify the relevant topological defect affecting the order parameter growth. We find that each phase is characterised by a distinct dynamic critical exponent. In the easy-plane phase we identify a persistent turbulent cascade that affects spin ordering long after all topological defects have annihilated. In addition to our exploration of phase ordering dynamics, we study a microscopic model of spin vortex dynamics in the easy-plane phase. Our work provides a comprehensive theoretical study of phase ordering in a conservative system, provides a thorough foundation for studies of phase ordering in antiferromagnetic and higher spin condensates, and offers prospects for further research into fundamental questions regarding the ordering properties of spin systems. Our work is pertinent to current experiments, which have explored the initial stages of phase ordering in both ferromagnetic and antiferromagnetic spin-1 condensates.
We consider the formation of order in a quasi-two-dimensional antiferromagnetic spin-1 condensate quenched from an easy-axis to an easy-plane nematic phase. We define the relevant order parameter to quantify the spin-nematic degrees of freedom and study the evolution of the spin-nematic and superfluid order during the coarsening dynamics using numerical simulations. We observe dynamical scaling in the late-time dynamics, with both types of order extending across the system with a diffusive growth law. We identify half-quantum vortices as the relevant topological defects of the ordering dynamics and demonstrate that the growth of both types of order is determined by the mutual annihilation of these vortices.
This thesis discusses non-equilibrium phenomena in antiferromagnetic spinor quantum gases. We investigate, both experimentally and theoretically, the quench dynamics of antiferromagnetic spinor Bose-Einstein condensates in the vicinity of a zero temperature quantum phase transition. A microwave dressing field was used as a novel tool to access the polar to antiferrmagnetic quantum phase transition at the quadratic Zeeman shift, q=0. By performing instantaneous quenches across the phase transition,a dynamical instability was induced in the atomic cloud. This lead to pair formation and rapid amplification of spin +/-1 atoms through spin-mixing collisions. The spatial ordering kinetics in the vicinity of the phase transition was quantified by performing the local spin measurement. The coarsening dynamics was observed for the instability for q < 0, as it nucleated small domains that grew to the axial size of the cloud. Rich nonequilibrium behavior was observed in the form of growth and decay of both nematic and magnetic spin waves, following the quench. The spatiotemporal evolution was characterized through two particle correlations between atoms in each pair of spin states, which revealed dramatic differences between the dynamics of the spin correlations and those of the spin populations. \nFinally, we demonstrate the remarkable influence of linear Zeeman term on non-equilibrium dynamics of spinor Bose-Einstein condensates. We show that contrary to prior understanding magnetic field gradient actually suppresses the tendency of spin 0 state to separate into spin +/-1 state. This reconciles the Bogoliubov discrepancies and leads to a dramatic sharpening of the transition point, resulting in a resolution at the 1 Hz level. Our results point to the use of dynamics, rather than equilibrium quantities for high precision measurements of phase transitions in quantum gases.
Due to the coexistence of phase coherence and magnetic order, Bose-Einstein condensates with internal degrees of freedom, i.e., spinor BECs, are considered to serve as an ideal table-top playground for the studies of various phenomena in different fields of physics.They have been widely used in studying spin textures and topological defects, coherent spin dynamics, quantum symmetry breaking, quantum phase transitions, etc.However, the theoretical framework used in these works is based on the mean-field theory in which the effects of quantum fluctuations are ignored.Despite the success of the mean-field theory in the description of spinor BECs, there are particular features whose account must rely on a beyond-mean-field theory as the effects of quantum fluctuations become significant.In this thesis, we develop the Beliaev theory, which is a Green's function approach beyond the Bogoliubov theory, for spinor BECs to study the effects of quantum fluctuations on the condensates' phase structures, elementary excitations, stabilities, and phase transitions.First, we point out that the Bogoliubov theory fails to capture the metastable states associated with first-order quantum phase transitions in spin-2 BECs.In fact, we show that these metastable states are induced by quantum fluctuations, and only by going to the next-order approximation, i.e., the Beliaev theory, can we show that they indeed appear around the phase boundaries.Besides the fluctuation-induced metastability, we find another important class of first-order quantum phase transitions in both spin-1 and spin-2 BECs.In this class, there is no metastable state to all orders of approximation since it is prohibited by a high symmetry of the Hamiltonian at the phase boundary, resulting in a characteristic flat energy landscape.Despite being first-order phase transitions, the flat energy landscape leads to the criticality in the dynamics of the condensate through these transitions in a manner similar to second-order phase transitions.Second, we succeed in deriving the analytic expression for the energy gap of the so-called quasi-Nambu-Goldstone (quasi-NG) modes, which are excitations not generated by spontaneous symmetry breaking, in the nematic phase of spin-2 BECs.Although quasi-NG modes are gapless at the mean-field level, we prove that with quantum corrections they acquire a nonzero energy gap.From the obtained magnitude of the energy gap, we can evaluate the critical temperature above which a topological defect such as a vortex of spin nematicity would decay by emitting thermally excited quasi-NG modes.We also study how the propagation of quasi-NG modes in space is affected by the particle-number density fluctuations of the condensate.Third, we calculate the damping rates of various types of quasiparticles in a spin-2 BEC including phonons, magnons, and quasi-NG modes.They actually have finite lifetime and decay via numerous channels of collision with the condensate atoms.Using either the Fermi's golden rule or the spinor Beliaev theory, we obtain the analytic expressions for the damping rates of phonons and magnons, from which the power-law dependence of the damping rates on the momentum is deduced.In contrast, the damping of quasi-NG modes is suppressed due to the energy conservation.
Since their experimental discovery in 1995, dilute Bose–Einstein condensates have proven to be excellent platforms for studying interacting many-body quantum systems. The gaseous condensates can be conveniently controlled in space and time with lasers and magnetic fields. Theoretically, these low-temperature systems can be described starting from first principles, and their dynamics is accurately captured by the mean-field approach. Dilute condensates with internal spin degrees of freedom may host various topological defects including vortices, monopoles, and skyrmions. Precise external control fields can be used to engineer and study these topological states. \n \nIn this thesis, we investigate the creation, stability, and dynamics of various topological defects in spinor Bose–Einstein condensates. A combination of numerical and analytical methods is used and a part of the research is carried out in collaboration with an experimental research group. The majority of the research involves numerically solving the mean-field Gross–Pitaevskii equation for a spin-1 condensate. The computations are accelerated with graphics processing units. \n \nWe find that the ground state of a ferromagnetic spin-1 condensate supports a Dirac monopole in its synthetic magnetic field in the presence of a quadrupole magnetic field. The energetics of various stationary states of spin-orbit coupled condensates are analyzed and a robust method to observe these exotic states in experiments is proposed. The creation of topological skyrmions is simulated and a very good quantitative agreement with a recent experiment is obtained. Importantly, the first experimental observations of both, Dirac monopoles and topological point defects in a system governed by a quantum field are presented. Subsequently, the decay dynamics of the point defect is studied, revealing the decay into a ferromagnetic state supporting the Dirac monopole. \n \nThe impact of the research conducted in this thesis is not limited to gaseous Bose–Einstein condensates, but due to the universality of quantum mechanics and topological defects, our work provides knowledge, ideas, and inspiration across research fields.
In this paper, we systematically review mathematical models, theories and numerical methods for ground states and dynamics of spinor Bose-Einstein condensates (BECs) based on the coupled Gross-Pitaevskii equations (GPEs). We start with a pseudo spin-1/2 BEC system with/without an internal atomic Josephson junction and spin-orbit coupling including (i) existence and uniqueness as well as non-existence of ground states under different parameter regimes, (ii) ground state structures under different limiting parameter regimes, (iii) dynamical properties, and (iv) efficient and accurate numerical methods for computing ground states and dynamics. Then we extend these results to spin-1 BEC and spin-2 BEC. Finally, extensions to dipolar spinor systems and/or general spin-F (F>=3) BEC are discussed.