Multiquantum vortices in dilute atomic Bose-Einstein condensates confined in long cigar-shaped traps are known to be both energetically and dynamically unstable. They tend to split into single-quantum vortices even in the ultralow temperature limit with vanishingly weak dissipation, which has also been confirmed in the recent experiments [Y. Shin et al., Phys. Rev. Lett. 93, 160406 (2004)] utilizing the so-called topological phase engineering method to create multiquantum vortices. We study the stability properties of multiquantum vortices in different trap geometries by solving the Bogoliubov excitation spectra for such states. We find that there are regions in the trap asymmetry and condensate interaction strength plane in which the splitting instability of multiquantum vortices is suppressed, and hence they are dynamically stable. For example, the doubly quantized vortex can be made dynamically stable even in spherical traps within a wide range of interaction strength values. We expect that this suppression of vortex-splitting instability can be experimentally verified.
Recently, the splitting of a topologically created doubly quantized vortex into two singly quantized vortices was experimentally investigated in dilute atomic cigar-shaped Bose-Einstein condensates [Y. Shin, Phys. Rev. Lett. 93, 160406 (2004)10.1103/PhysRevLett.93.160406]. In particular, the dependency of the splitting time on the peak particle density was studied. We present results of theoretical simulations which closely mimic the experimental setup. We show that the combination of gravitational sag and time dependency of the trapping potential alone suffices to split the doubly quantized vortex in time scales which are in good agreement with the experiments.
We study the instability of a doubly quantized vortex topologically imprinted on 23 Na condensate, as reported in recent experiment (Shin et al 2004 Phys. Rev. Lett. 93 160406). We have performed numerical simulations using three-dimensional Gross–Pitaevskii equation with classical thermal noise. Splitting of a doubly quantized vortex turns out to be a process that is very sensitive to the presence of thermal atoms. We observe that even very small thermal fluctuations, corresponding to 10–15% of thermal atoms, cause the decay of doubly quantized vortex into two singly quantized vortices in tens of milliseconds. As in the experiment, the lifetime of doubly quantized vortex is a monotonic function of the interaction strength.
We study the nonlinear dynamics of the splitting of a doubly quantized vortex in a trapped Bose-Einstein condensate. The dynamics is studied in detail by solving the Gross-Pitaevskii equation. The main dynamical features are explained in terms of a nonlinear three-level system, to which we find an analytical solution. This result is compared to linear Bogoliubov analysis and to the full numerical time evolution. It is concluded that the time scale for the splitting is mainly determined by the instability of the linearized system, and nonlinear effects contribute logarithmically.
The dynamical instabilities and ensuing dynamics of singly and doubly quantized vortex states of Bose-Einstein condensates with attractive interactions are investigated using full three-dimensional numerical simulations of the Gross-Pitaevskii equation. With increasing the strength of attractive interactions, a series of dynamical instabilities such as quadrupole, dipole, octupole, and monopole instabilities emerge. The most prominent instability depends on the strength of interactions, the geometry of the trapping potential, and deviations from the axisymmetry due to external perturbations. Singly quantized vortices split into two clusters and subsequently undergo split-merge cycles in a pancake-shaped trap, whereas the split fragments immediately collapse in a spherical trap. Doubly quantized vortices are always unstable to disintegration of the vortex core. If we suddenly change the strength of interaction to within a certain range, the vortex splits into three clusters, and one of the clusters collapses after a few split-merge cycles. The vortex split can be observed using a current experimental setup of Leanhardt et al. [Phys. Rev. Lett. 89, 190403 (2002)].
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.