We review the topic of quantized vortices in multicomponent Bose–Einstein condensates of dilute atomic gases, with an emphasis on the two-component condensates. First, we review the fundamental structure, stability and dynamics of a single vortex state in a slowly rotating two-component condensates. To understand recent experimental results, we use the coupled Gross–Pitaevskii equations and the generalized nonlinear sigma model. An axisymmetric vortex state, which was observed by the JILA group, can be regarded as a topologically trivial skyrmion in the pseudospin representation. The internal, coherent coupling between the two components breaks the axisymmetry of the vortex state, resulting in a stable vortex molecule (a meron pair). We also mention unconventional vortex states and monopole excitations in a spin-1 Bose–Einstein condensate. Next, we discuss a rich variety of vortex states realized in rapidly rotating two-component Bose–Einstein condensates. We introduce a phase diagram with axes of rotation frequency and the intercomponent coupling strength. This phase diagram reveals unconventional vortex states such as a square lattice, a double-core lattice, vortex stripes and vortex sheets, all of which are in an experimentally accessible parameter regime. The coherent coupling leads to an effective attractive interaction between two components, providing not only a promising candidate to tune the intercomponent interaction to study the rich vortex phases but also a new regime to explore vortex states consisting of vortex molecules characterized by anisotropic vorticity. A recent experiment by the JILA group vindicated the formation of a square vortex lattice in this system.
We investigate vortex states of immiscible two-component Bose-Einstein condensates under rotation through numerical simulations of the coupled Gross-Pitaevskii equations. For strong intercomponent repulsion, the two components undergo phase separation to form several density domains of the same component. In the presence of the rotation, the nucleated vortices are aligned between the domains to make up winding chains of singly quantized vortices, a vortex sheet, instead of periodic vortex lattices. The vortices of one component are located at the region of the density domains of the other component, which results in the serpentine domain structure. The sheet configuration is stable as long as the imbalance of the intracomponent parameter is small. We employ a planar sheet model to estimate the distance between neighboring sheets, determined by the competition between the surface tension of the domain wall and the kinetic energy of the superflow via quantized vortices. By comparing the several length scales in this system, the phase diagram of the vortex state is obtained.
We investigate the structure of vortex states in rotating two-component Bose-Einstein condensates with equal intracomponent but varying intercomponent-coupling constants. A phase diagram in the intercomponent-coupling versus rotation-frequency plane reveals rich equilibrium structures of vortex states. As the ratio of intercomponent to intracomponent couplings increases, the interlocked vortex lattices undergo phase transitions from triangular to square, to double-core lattices, and eventually develop interwoven "serpentine" vortex sheets with each component made up of chains of singly quantized vortices.
When a two-component Bose-Einstein condensate is placed into rotation, a lattice of vortices and cores appear. The geometry of this lattice (triangular or square) varies according to the rotational value and the intercomponent coupling strengths. In this paper, assuming a Thomas-Fermi regime, we derive a point energy that allows us to determine for which values of the parameters the lattice goes from triangular to square. It turns out that the separating curve in the phase diagram agrees fully with the complete numerical simulations of the Gross-Pitaevskii equations. We also derive a formula for the critical velocity of appearance of the first vortex and prove that the first vortex always appears first in the component with largest support in the case of two disks, and we give a criterion in the case of disk and annulus.
We investigate two kinds of coreless vortices with axisymmetric and nonaxisymmetric configurations in rotating two-component Bose-Einstein condensates. Starting from the Gross-Pitaevskii energy functional in a rotating frame, we derive a nonlinear sigma model generalized to the two-component condensates. In terms of a pseudospin representation, an axisymmetric vortex and a nonaxisymmetric one correspond to spin textures referred to as a ``skyrmion'' and a ``meron-pair,'' respectively. A variational method is used to investigate the dependence of the sizes of the stable spin textures on system parameters, and the optimized variational function is found to reproduce well the numerical solution. In the SU(2) symmetric case, the optimal skyrmion and meron-pair are degenerate and transform to each other by a rotation of the pseudospin. An external rf field that couples coherently the hyperfine states of two components breaks the degeneracy in favor of the meron-pair texture due to an effective transverse pseudomagnetic field. The difference between the intracomponent and intercomponent interactions yields a longitudinal pseudomagnetic field and a ferromagnetic or an antiferromagnetic pseudospin interaction, leading to a meron-pair texture with an anisotropic distribution of vorticity.
We study numerically the dynamical instabilities and splitting of singly and doubly quantized composite vortices in two-component Bose-Einstein condensates harmonically confined to quasi two dimensions. In this system, the vortices become pointlike composite defects that can be classified in terms of an integer pair $({\ensuremath{\kappa}}_{1},\phantom{\rule{0.16em}{0ex}}{\ensuremath{\kappa}}_{2})$ of phase winding numbers. Our simulations based on zero-temperature mean-field theory reveal several vortex splitting behaviors that stem from the multicomponent nature of the system and do not have direct counterparts in single-component condensates. By calculating the Bogoliubov excitations of stationary axisymmetric composite vortices, we find nonreal excitation frequencies (i.e., dynamical instabilities) for the singly quantized
(1,1) and $(1,\phantom{\rule{0.16em}{0ex}}\ensuremath{-}1)$ vortices and for all variants of doubly quantized vortices, which we define by the condition ${max}_{j=1,2}|{\ensuremath{\kappa}}_{j}|=2$. While the short-time predictions of the linear Bogoliubov analysis are confirmed by direct time integration of the Gross-Pitaevskii equations of motion, the time integration also reveals intricate long-time decay behavior not captured by the linearized dynamics. First, the $(1,\phantom{\rule{0.16em}{0ex}}\ifmmode\pm\else\textpm\fi{}1)$ vortex is found to be unstable against splitting into a
(1,0) vortex and a $(0,\phantom{\rule{0.16em}{0ex}}\ifmmode\pm\else\textpm\fi{}1)$ vortex. Second, the
(2,1) vortex exhibits a two-step decay process in which its initial splitting into a
(2,0) vortex and a
(0,1) vortex is followed by the off-axis splitting of the
(2,0) vortex into two
(1,0) vortices. Third, the $(2,\phantom{\rule{0.16em}{0ex}}\ensuremath{-}2)$ vortex is observed to split into a $(\ensuremath{-}1,\phantom{\rule{0.16em}{0ex}}1)$ vortex, three
(1,0) vortices, and three $(0,\phantom{\rule{0.16em}{0ex}}\ensuremath{-}1)$ vortices. Each of these splitting processes is the dominant decay mechanism of the respective stationary composite vortex for a wide range of intercomponent interaction strengths and relative populations of the two condensate components and should be amenable to experimental detection. Our results contribute to a better understanding of vortex physics, hydrodynamic instabilities, and two-dimensional quantum turbulence in multicomponent superfluids.