Topological aspects in spinor Bose-Einstein condensates

Integrator for general spin-s Gross-Pitaevskii systems

We provide an algorithm, i-SPin 2, for evolving general spin-s Gross-Pitaevskii or nonlinear Schrödinger systems carrying a variety of interactions, where the 2s+1 components of the "spinor" field represent the different spin-multiplicity states. We consider many nonrelativistic interactions up to quartic order in the Schrödinger field (both short and long range, and spin-dependent and spin-independent interactions), including explicit spin-orbit couplings. The algorithm allows for spatially varying external and/or self-generated vector potentials that couple to the spin density of the field. Our work can be used for scenarios ranging from laboratory systems such as spinor Bose-Einstein condensates (BECs), to cosmological or astrophysical systems such as self-interacting bosonic dark matter. As examples, we provide results for two different setups of spin-1 BECs that employ a varying magnetic field and spin-orbit coupling, respectively, and also collisions of spin-1 solitons in dark matter. Our symplectic algorithm is second-order accurate in time, and is extensible to the known higher-order-accurate methods.

A numerical scheme for the ground state of rotating spin-1 Bose-Einstein condensates

We study the existence of nontrivial solution branches of three-coupled Gross-Pitaevskii equations (CGPEs), which are used as the mathematical model for rotating spin-1 Bose-Einstein condensates (BEC). The Lyapunov-Schmidt reduction is exploited to test the branching of nontrivial solution curves from the trivial one in some neighborhoods of bifurcation points. A multilevel continuation method is proposed for computing the ground state solution of rotating spin-1 BEC. By properly choosing the constraint conditions associated with the components of the parameter variable, the proposed algorithm can effectively compute the ground states of spin-1 [Formula: see text] and [Formula: see text] under rapid rotation. Extensive numerical results demonstrate the efficiency of the proposed algorithm. In particular, the affect of the magnetization on the CGPEs is investigated.

Vortex-lattice formation in a spin-orbit coupled rotating spin-1 condensate

We study the vortex-lattice formation in a rotating Rashba spin-orbit (SO) coupled quasi-two-dimensional (quasi-2D) hyper-fine spin-1 spinor Bose-Einstein condensate (BEC) in the x-y plane using a numerical solution of the underlying mean-field Gross-Pitaevskii equation. In this case, the non-rotating Rashba SO-coupled spinor BEC can have topological excitation in the form of vortices of different angular momenta in the three components, e.g. the (0, +1, +2)- and (-1, 0, +1)-type states in ferromagnetic and anti-ferromagnetic spinor BEC: the numbers in the parenthesis denote the intrinsic angular momentum of the vortex states of the three components with the negative sign denoting an anti-vortex. The presence of these states with intrinsic vorticity breaks the symmetry between rotation with vorticity along the z and -z axes and thus generates a rich variety of vortex-lattice and anti-vortex-lattice states in a rotating quasi-2D spin-1 spinor ferromagnetic and anti-ferromagnetic BEC, not possible in a scalar BEC. For weak SO coupling, we find two types of symmetries of these states - hexagonal and 'square'. The hexagonal (square) symmetry state has vortices arranged in closed concentric orbits with a maximum of 6, 12, 18… (8, 12, 16…) vortices in successive orbits. Of these two symmetries, the square vortex-lattice state is found to have the smaller energy.

Spin-Nematic Vortex States in Cold Atoms

The (pseudo)spin degrees of freedom greatly enriches the physics of cold atoms. This is particularly so for systems with high spins (i.e., spin quantum number larger than 1/2). For example, one can construct not only the rank-1 spin vector, but also the rank-2 spin tensor in high spin systems. Here we propose a simple scheme to couple the spin tensor and the center-of-mass orbital angular momentum in a spin-1 cold atom system and show that this leads to a new quantum phase of the matter: the spin-nematic vortex state that features vorticity in an SU(2) spin-nematic tensor subspace. Under proper conditions, such states are characterized by quantized topological numbers. Our work opens up new avenues of research in topological quantum matter with high spins.

Double-layer Bose-Einstein condensates: A quantum phase transition in the transverse direction, and reduction to two dimensions

We revisit the problem of the reduction of the three-dimensional (3D) dynamics of Bose-Einstein condensates, under the action of strong confinement in one direction (z), to a 2D mean-field equation. We address this problem for the confining potential with a singular term, viz., V_{z}(z)=2z^{2}+ζ^{2}/z^{2}, with constant ζ. A quantum phase transition is induced by the latter term, between the ground state (GS) of the harmonic oscillator and the 3D condensate split in two parallel noninteracting layers, which is a manifestation of the "superselection" effect. A realization of the respective physical setting is proposed, making use of resonant coupling to an optical field, with the resonance detuning modulated along z. The reduction of the full 3D Gross-Pitaevskii equation (GPE) to the 2D nonpolynomial Schrödinger equation (NPSE) is based on the factorized ansatz, with the z -dependent multiplier represented by an exact GS solution of the 1D Schrödinger equation with potential V_{z}(z). For both repulsive and attractive signs of the nonlinearity, the 2D NPSE produces GS and vortex states, that are virtually indistinguishable from the respective numerical solutions provided by full 3D GPE. In the case of the self-attraction, the threshold for the onset of the collapse, predicted by the 2D NPSE, is also virtually identical to its counterpart obtained from the 3D equation. In the same case, stability and instability of vortices with topological charge S=1, 2, and 3 are considered in detail. Thus, the procedure of the spatial-dimension reduction, 3D → 2D, produces very accurate results, and it may be used in other settings.

Stripe and supersolid phases of spin-orbit coupled spin-2 Bose-Einstein condensates in an optical lattice

We study the ground-state phases of two-dimensional spin-orbit coupled spin-2 Bose-Einstein condensates in a one-dimensional spin-dependent optical lattice. Due to the competition among optical lattice, spin-orbit coupling and spin-exchange interaction, the exotic ground-state phases are found, i.e. three types of the stripe phases and three types of the supersolid phases. The spin-exchange interaction can adjust the direction of the stripe in the stripe phase and generate various vortex lattice structures in the supersolid phase, which shows that the spin-exchange interaction plays an important role in the formation of the stripe and supersolid phases of spin-orbit coupled spin-2 Bose-Einstein condensates in an optical lattice.

Rotating Atomic Quantum Gases with Light-Induced Azimuthal Gauge Potentials and the Observation of the Hess-Fairbank Effect

We demonstrate synthetic azimuthal gauge potentials for Bose-Einstein condensates from engineering atom-light couplings. The gauge potential is created by adiabatically loading the condensate into the lowest energy Raman-dressed state, achieving a coreless vortex state. The azimuthal gauge potentials act as effective rotations and are tunable by the Raman coupling and detuning. We characterize the spin textures of the dressed states, in agreements with the theory. The lowest energy dressed state is stable with a 4.5-s half-atom-number-fraction lifetime. In addition, we exploit the azimuthal gauge potential to demonstrate the Hess-Fairbank effect, the analogue of Meissner effect in superconductors. The atoms in the absolute ground state has a zero quasiangular momentum and transits into a polar-core vortex when the synthetic magnetic flux is tuned to exceed a critical value. Our demonstration serves as a paradigm to create topological excitations by tailoring atom-light interactions where both types of SO(3) vortices in the |⟨F[over →]⟩|=1 manifold, coreless vortices and polar-core vortices, are created in our experiment. The gauge field in the stationary Hamiltonian opens a path to investigating rotation properties of atomic superfluids under thermal equilibrium.

Spin-Orbital-Angular-Momentum Coupled Bose-Einstein Condensates

We demonstrate coupling between the atomic spin- and orbital-angular momentum (OAM) of the atom's center-of-mass motion in a Bose-Einstein condensate (BEC). The coupling is induced by Raman-dressing lasers with a Laguerre-Gaussian beam and creates coreless vortices in an F=1 ^{87}Rb spinor BEC. We observe correlations between spin and OAM in the dressed state and characterize the spin texture; the result is in good agreement with the theory. In the presence of the Raman field, our dressed state is stable for 0.1 s or longer, and it decays due to collision-induced relaxation. As we turn off the Raman beams, the vortex cores in the bare spin |m_{F}=1⟩ and |-1⟩ split. These spin-OAM coupled systems with the Raman-dressing approach have great potential for exploring new topological textures and quantum states.