The formation of an equilibrium quantum state from an uncorrelated thermal one through the dynamical crossing of a phase transition is a central question of non-equilibrium many-body physics. During such crossing, the system breaks its symmetry by establishing numerous uncorrelated regions separated by spontaneously-generated defects, whose emergence obeys a universal scaling law with the quench duration. Much less is known about the ensuing re-equilibrating or "coarse-graining" stage, which is governed by the evolution and interactions of such defects under system-specific and external constraints. In this work we perform a detailed numerical characterization of the entire non-equilibrium process, addressing subtle issues in condensate growth dynamics and demonstrating the quench-induced decoupling of number and coherence growth during the re-equilibration process. Our unique visualizations not only reproduce experimental measurements in the relevant regimes, but also provide valuable information in currently experimentally-inaccessible regimes.

We explore a way of producing the Rashba spin-orbit coupling (SOC) for ultracold atoms by using a two-component (spinor) atomic Bose-Einstein condensate (BEC) confined in a bilayer geometry. The SOC of the Rashba type is created if the atoms pick up a pi phase after completing a cyclic transition between four combined spin-layer states composed of two spin and two layer states. The cyclic coupling of the spin-layer states is carried out by combining an intralayer Raman coupling and an interlayer laser assisted tunneling. We theoretically determine the ground-state phases of the spin-orbit-coupled BEC for various strengths of the atom-atom interaction and the laser-assisted coupling. It is shown that the bilayer scheme provides a diverse ground-state phase diagram. In an intermediate range of the atom-light coupling two interlacing lattices of half-skyrmions and half-antiskyrmions are spontaneously created. In the strong-coupling regime, where the SOC of the Rashba type is formed, the ground state represents plane-wave or standing-wave phases depending on the interaction between the atoms. A variational analysis is shown to be in good agreement with the numerical results.

We study the sensitivity of coupled condensate formation dynamics on the history of initial stochastic domain formation in the context of instantaneously quenched elongated harmonically trapped immiscible two-component atomic Bose gases. The spontaneous generation of defects in the fastest condensing component, and subsequent coarse-graining dynamics, can lead to a deep oscillating microtrap into which the other component condenses, thereby establishing a long-lived composite defect in the form of a dark-bright solitary wave. We numerically map out diverse key aspects of these competing growth dynamics, focusing on the role of shot-to-shot fluctuations and global parameter changes (initial state choices, quench parameters, and condensate growth rates), with our findings also qualitatively confirmed by realistic finite-duration quenches. We conclude that phase-separated structures observable on experimental time scales are likely to be metastable states whose form is influenced by the stability and dynamics of the spontaneously emerging dark-bright solitary wave.

We theoretically explore atomic Bose–Einstein condensates (BECs) subject to position-dependent spin–orbit coupling (SOC). This SOC can be produced by cyclically laser coupling four internal atomic ground (or metastable) states in an environment where the detuning from resonance depends on position. The resulting spin–orbit coupled BEC (SOBEC) phase separates into domains, each of which contain density modulations—stripes—aligned either along the x or y direction. In each domain, the stripe orientation is determined by the sign of the local detuning. When these stripes have mismatched spatial periods along domain boundaries, non-trivial topological spin textures form at the interface, including skyrmions-like spin vortices and anti-vortices. In contrast to vortices present in conventional rotating BECs, these spin-vortices are stable topological defects that are not present in the corresponding homogenous stripe-phase SOBECs.

Abstract Transcritical flow of a stratified fluid past a broad localised topographic obstacle is studied analytically in the framework of the forced extended Korteweg–de Vries, or Gardner, equation. We consider both possible signs for the cubic nonlinear term in the Gardner equation corresponding to different fluid density stratification profiles. We identify the range of the input parameters: the oncoming flow speed (the Froude number) and the topographic amplitude, for which the obstacle supports a stationary localised hydraulic transition from the subcritical flow upstream to the supercritical flow downstream. Such a localised transcritical flow is resolved back into the equilibrium flow state away from the obstacle with the aid of unsteady coherent nonlinear wave structures propagating upstream and downstream. Along with the regular, cnoidal undular bores occurring in the analogous problem for the single-layer flow modelled by the forced Korteweg–de Vries equation, the transcritical internal wave flows support a diverse family of upstream and downstream wave structures, including kinks, rarefaction waves, classical undular bores, reversed and trigonometric undular bores, which we describe using the recent development of the nonlinear modulation theory for the (unforced) Gardner equation. The predictions of the developed analytic construction are confirmed by direct numerical simulations of the forced Gardner equation for a broad range of input parameters.

We present an overview of our recent theoretical studies on the quantum phenomena of the spin-1 Bose-Einstein condensates, including the phase diagram, soliton solutions and the formation of the topological spin textures. A brief exploration of the effects of spin-orbit coupling on the ground-state properties is given. We put forward proposals by using the transmission spectra of an optical cavity to probe the quantum ground states: the ferromagnetic and polar phases. Quasi-one-dimension solitons and ring dark solitons are studied. It is predicted that characteristics of the magnetic solitons in optical lattice can be tuned by controlling the long-range light-induced and static magnetic dipoledipole interactions; solutions of single-component magnetic and single-, two-, three-components polar solitons are found; ring dark solitons in spin-1 condensates are predicted to live longer lifetimes than that in their scalar counterparts. In the formation of spin textures, we have considered the theoretical model of a rapidly quenched and fast rotating trapped spin-1 Bose-Einstein condensate, whose dynamics can be studied by solving the stochastic projected Gross-Pitaevskii equations. Spontaneous generation of nontrivial topological defects, such as the hexagonal lattice skyrmions and square lattice of half-quantized vortices was predicted. In particular, crystallization of merons (half skyrmions) can be generated in the presence of spin-orbit coupling.

We develop modulation theory for undular bores (dispersive shock waves) in the framework of the Gardner, or extended Korteweg-de Vries (KdV), equation, which is a generic mathematical model for weakly nonlinear and weakly dispersive wave propagation, when effects of higher order nonlinearity become important. Using a reduced version of the finite-gap integration method we derive the Gardner-Whitham modulation system in a Riemann invariant form and show that it can be mapped onto the well-known modulation system for the Korteweg-de Vries equation. The transformation between the two counterpart modulation systems is, however, not invertible. As a result, the study of the resolution of an initial discontinuity for the Gardner equation reveals a rich phenomenology of solutions which, along with the KdV-type simple undular bores, include nonlinear trigonometric bores, solibores, rarefaction waves, and composite solutions representing various combinations of the above structures. We construct full parametric maps of such solutions for both signs of the cubic nonlinear term in the Gardner equation. Our classification is supported by numerical simulations.

The non-equilibrium dynamics of a rapidly quenched spin-1 Bose gas with spin-orbit coupling is studied. By solving the stochastic projected Gross-Pitaevskii equation, we show that crystallization of merons can occur in a spinor condensate of ^{87}Rb. Analytic form and stability of the crystal structure are given. Likewise, inverted merons can be created in a spin-polarized spinor condensate of ^{23}Na. Our studies provide a chance to explore the fundamental properties of meron-like matter.

Equilibrium vortex formation in rotating binary Bose gases with a rotating frequency higher than the harmonic trapping frequency is investigated theoretically. We consider the system being evaporatively cooled to form condensates and a combined numerical scheme is applied to ensure the binary system being in an authentic equilibrium state. To keep the system stable against the large centrifugal force of ultrafast rotation, a quartic trapping potential is added to the existing harmonic part. Using the Thomas-Fermi approximation, a critical rotating frequency Omega(c) is derived, which characterizes the structure with or without a central density hole. Vortex structures are studied in detail with rotation frequency both above and below Omega(c) and with respect to the miscible, symmetrically separated, and asymmetrically separated phases in their nonrotating ground-state counterparts.

We investigate the spontaneous generation of crystallized topological defects via the combining effects of fast rotation and rapid thermal quench on spin-1 Bose-Einstein condensates (BECs). By solving the stochastic projected Gross-Pitaevskii equation, we show that, when the system reaches equilibrium, a hexagonal lattice of skyrmions and a square lattice of half-quantized vortices can be formed in a ferromagnetic and antiferromagnetic spinor BEC, respectively, which can be imaged by using the polarization-dependent phase-contrast method.

We investigate the dynamics of turbulent flow in a two-dimensional trapped Bose-Einstein condensate by solving the Gross-Pitaevskii equation numerically. The development of the quantum turbulence is activated by the disruption of an initially embedded vortex quadrupole. By calculating the incompressible kinetic-energy spectrum of the superflow, we conclude that this quantum turbulent state is characterized by the Kolmogorov-Saffman scaling law in the wave-number space. Our study predicts the coexistence of two subinertial ranges responsible for the energy cascade and enstrophy cascade in this prototype of two-dimensional quantum turbulence.

Formation of stationary three-dimensional (3D) wave patterns generated by a small pointlike impurity moving through a Bose-Einstein condensate with supersonic velocity is studied. Asymptotic formulae for a stationary far-field density distribution are obtained. Comparison with the three-dimensional numerical simulations demonstrates that these formulae are accurate enough already at the distances from the obstacle equal to a few wavelengths.

We investigate the dynamics of an unstable vortex ring in a pancake-shaped Bose-Einsten condensate by solving the Gross-Pitaevskii equation numerically. It is found that a quasisteady turbulent state with long relaxation time can be achieved through the disruption of a perturbed vortex ring in the condensate owing to the bending-wave instability. We verify that this quantum turbulent state is characterized by Kolmogorov energy spcetrum.

We propose a numerical scheme for obtaining the stationary vortex-ring solutions of the Gross–Pitaevskii (GP) equation for an axisymmetrically trapped Bose–Einstein condensate (BEC). The effective energy functional and the associated GP equation are derived by assuming a trial phase profile for the wavefunction that is subject to the condition of circulation quantization on the rz plane. The wavefunction of the vortex ring is determined by solving the ground state of the effective GP equation numerically. Application of our method to the formation of a three-dimensional Skyrmion in a trapped two-component BEC is demonstrated.

Based on a velocity formula derived by matched asymptotic expansion, we investigate the dynamics of a circular vortex ring in an axisymmetric Bose-Einstein condensate in the Thomas-Fermi limit. The trajectory for an axisymmetrically placed and oriented vortex ring is entirely determined, revealing that the vortex ring generally precesses in condensate. The linear instability due to bending waves is investigated both numerically and analytically. General stability boundaries for various perturbed wavenumbers are computed. In particular, the excitation spectrum and the absolutely stable region for the static ring are analytically determined.

We investigate theoretically the transition temperature of a weakly interacting F=1 spinor Bose-Einstein condensate (BEC) formed in an optical trap. The Hartree-Fock-Bogoliubov formalism for the spinor BEC is developed based on the mean-field theory for inhomogeneous systems at finite temperatures. The shift of the transition temperature due to finite-size effects and the two-body interaction is obtained simultaneously by analytically solving the coupled eigenvalue equations for quasiparticle excitations and amplitudes under Hartee-Fock-Bogoliubov and other relevant approximations. When applying this to the recent experiment of the all-optical formation of a spinor BEC with a condensate of 3.5x10(4) Rb-87 atoms [Phys. Rev. Lett. 87, 010404 (2001)], our calculations show that the condensation occurs at a temperature about 177 nK. An analytic result for the shift of the transition temperature in the presence of a weak magnetic field is also derived.

A hydrodynamic description is used to study the zero-temperature properties of a trapped spinor Bose-Einstein condensate in the presence of a uniform magnetic field. We show that, in the case of antiferromagnetic spin-spin interaction, the polar and ferromagnetic configurations of the ground state can coexist in the trap. These two phases are spatially segregated in such a way that the polar state occupies the inner part while the ferromagnetic state occupies the outer part of the atomic cloud. We also derive a set of coupled hydrodynamic equations for the number density and spin density excitations of the system. It is shown that these equations can be analytically solved for the system in an isotropic harmonic trap and a constant magnetic field. Remarkably, the related low lying excitation spectra are completely determined by the solutions in the region occupied by the polar state. We find that, within the Thomas-Fermi approximation, the presence of a constant magnetic field does not change the excitation spectra which still possess the similar form of that obtained by Stringari.

We present a method for the direct measurement of the Wigner-function matrix for complex vibronic states of a trapped atom, that is suited to analyse the entanglement between two motional degrees of freedom and the internal electronic dynamics. It is a generalisation of the method for the determination of vibronic quantum states [S. Wallentowitz, R.L. de Matos Filho, W. Vogel, Phys. Rev. A 56, 1205 (1997)] in conjunction with the scheme for the direct observation of the Wigner function of a single motional degree of freedom [L.G. Lutterbach, L. Davidovich, Phys. Rev. Lett. 78, 2547 (1997)]. The major advantage of the present method is that it reduces the experimental efforts substantially. On the other hand, it is demonstrated that the nonlinear vibronic coupling necessary for this method turns out to be its main limitation.