Vortex Lattice Locking in Rotating Two-Component Bose-Einstein Condensates

The vortex density of a rotating superfluid, divided by its particle mass, dictates the superfluid's angular velocity through the Feynman relation. To find how the Feynman relation applies to superfluid mixtures, we investigate a rotating two-component Bose-Einstein condensate, composed of bosons with different masses. We find that in the case of sufficiently strong interspecies attraction, the vortex lattices of the two condensates lock and rotate at the drive frequency, while the superfluids themselves rotate at two different velocities, whose ratio is the ratio between the particle mass of the two species. In this paper, we characterize the vortex-locked state, establish its regime of stability, and find that it surives within a disk smaller than a critical radius, beyond which vortices become unbound, and the two Bose-gas rings rotate together at the frequency of the external drive.

After the first experimental realization of Bose-Einstein condensates (BECs) of alkali atoms, their study has experienced enormous advancements [1]. Among the major threads of investigation in BECs has been the study of vortices both experimentally [2,3,4] and theoretically [5]. It is known from the classic works of Onsager and Feynman [6,7] that superfluids rotate by nucleating vortices. When there are several vortices present, they form a triangular Abrikosov vortex lattice [8], with density given by where m is the mass of a constituent boson and Ω is the rate at which the superfluid -which rotates with the vortex lattice -is being rotated (see, for instance, Ref. [9]). The so-called "Feynman relation" (1) states that, on average, a uniform superfluid rotates like a rigid body. It has been shown that corrections to Eq. (1) due to the typical experimental situation of nonuniform superfluid density resulting from a harmonic trap are small [10] but experimentally observable [11]. Vortex physics becomes much more intriguing in multi-component BEC's. So far, the investigation of vortex lattices in multicomponent BECs utilized different hyperfine levels of the constituent atoms to obtain multi-component condensates (e.g., [12,13]). Thus the mass of all condensate components was identical, and the generalization of Eq. (1) straightforward.
In this paper we ask what are the consequences of the Feynman condition, Eq. (1), in a system of interacting two-component rotating condensates with different masses. We find that for sufficiently large attractive interactions, the Feynman condition leads to a novel state. The two components, rather than rotate together at the drive frequency, rotate at angular velocities Ω 1,2 inversely proportional to their masses, m 1,2 , such that: while the vortex lattices of the two components lock at the drive angular velocity Ω v , lying between Ω 1 and Ω 2 (see Fig. 1). Qualitatively, the attractive interspecies interaction leads to an attraction between vortices of the two flavors. If it is sufficiently strong, vortices pair, and the lowest-energy vortex configuration then occurs when the vortex lattices of the two flavors are "locked" together, rotate at the same rate, and have essentially the same density. Eq. (1) then reads: where ρ 1,2 v are the vortex densities of the two flavors. This state is strongly related to experiments in Ref. 14, where a vortex lattice was locked to an optical lattice with a similar periodicity.
As we show below, this state survives within a finite disk about the center of the rotating condensate. The relative motion between the vortices and the condensate gives rise to a Magnus force that opposes the interspecies vortex attraction. Beyond a critical distance the Magnus force ( Fig. 1) becomes larger then the maximal pairing force, and the vortices become unbound. In this region, the two condensates and their vortices rotate together at the drive frequency Ω v ; the vortex densities in the two flavors are no longer equal, but reflect themselves the mass ratio: Below we derive the characteristics and conditions for the vortex locking state.
Energetics. The energy of weakly interacting Bose-Einstein condensates is well-described by the Gross-Pitaevskii functional for the condensate wave function ψ α (α = 1, 2) for the two atomic species [1]. We consider the situation in which the two condensates are stirred at the same rate Ω v . Transforming to the rotating frame, our problem becomes time independent. The energy of the two-component system in the rotating frame is given by describes the energy for each atomic species and g α is the intraspecies coupling for bosons of flavor α. The interspecies interaction is given by where n α = |ψ α | 2 is the density of flavor α, V trap is the external trapping potential, and the z-component of the angular momentum operator is L z = −i ∂ ∂ϕ (where ϕ is the azimuthal coordinate).
The energy E α of one BEC component is minimized via the nucleation of a vortex lattice rotating with the external drive Ω v . In the following, we assume that the coherence length ξ α = h 2 2mαgαn α 0 is much shorter than the characteristic distance between vortices. Each vortex is then well-described by a small core region of size ξ α , at which the superfluid density drops to zero and its phase field accounts for the flow around the vortex. The vortex core region gives rise to a small constant energy, and we can account for the phase field by writing ψ α = √ n α e iθα , where θ α determines a lattice of vortices with unit winding number at the positions {r α i }. We assume that we can write θ α as a sum of the different vortex contributions θ α = θ α 1 + θ α 2 + . . . . With these assumptions, E α can be written in terms of the positions of the vortices as where we have dropped terms that do not depend on the positions of the vortices and have also neglected effects due to the nonuniform superfluid density [10]. The first term in Eq. (6) is the usual logarithmic interaction between vortices, and the second is the centripetal energy, reflecting the fact that vortices toward the edge of the cloud carry less angular momentum relative to the center of the cloud. In a single-component rotating BEC, the balance of the two terms gives the Feynman condition (1). The equations describe charged particles interacting in two dimensions with a uniform background charge of opposite sign.
The energy E 12 arising from the interspecies interaction energy is less straightforward to evaluate. Unlike the intraspecies logarithmic interaction, this nonuniversal interaction depends on the details of the short-distance density variations around the vortex cores. For instance, in [15] to study the interaction of a vortex with an optical lattice, a step function having the width of the BEC coherence length was taken. In this work we take a Gaussian depletion around the vortex core: so that the system will be amenable to analytic treatment. For a single vortex this depletion gives the correct behavior at short distances, but not the long distance behavior, in which the density due to a single isolated vortex heals as ξ 2 /r 2 . As we show in the Appendix, the combined density variations due to the vortex lattice on scales larger than the inter-vortex separation combines only to change the chemical potential (which is proportional to the density correction), byh 2 2m π 2 ρ 2 v r 2 , which just reflects the kinetic energy associated with uniform rotation of the condensate. This can be shown to have a negligible effect on the vortex pairs, which are the focus of this work. Neglecting this piece of the density fluctuation is also consistent with neglecting the intraspecies core-core interactions which is a standard approximation [5]. Indeed, the short distance region is the relevant one for vortex locking; two-species vortex pairs become unbound once their separation is comparable to the coherence lengths, where the approximate form we take for the density profiles is still valid. Evaluating the interaction integral in Eq. (5), and keeping only the contributions due to the interactions between pairs of vortices between different species, gives We note that to prevent phase separation, the criterion |g 1 g 2 | > |g 12 | 2 must be satisfied. Equations (6) and (8) now give the energetics of the system purely in terms of the positions of the vortices.
Forces. An understanding of the locked phases can be obtained by considering the forces acting on the vortices. Equation (6) leads to the well-known Magnus force acting on a vortex of species α [9]: whereκ is a unit vector pointing out of the plane, n α 0 is the equilibrium superfluid density for species α (evaluated away from the vortex core), v α SF =h mα ∇θ α is the superfluid velocity for species α (with the vortex on which the force operates excluded from θ α ), and v v is the velocity of the vortex. On the other hand, the force arising from the energy in Eq. (8) provides an attractive force between two vortices of different species. It has the form where d is the displacement vector between vortices. Let us first briefly consider the unlocked case where the vortex interspecies interaction force is small. Since the force counteracting the Magnus force is too small, to bind vortex pairs, we must have F mag = 0 for any isolated vortex, which implies that the superfluid velocity must be the same as the vortex velocity. That is, the vortex lattice rotates with the superfluid. Thus, because the two vortex lattices rotate at the same frequency, the two superfluids rotate together at that frequency. Accordingly, for this case, the vortex densities are not equal: We next consider the other extreme, in which the two vortex lattices are locked. Our approach is to consider the forces acting on a bound pair of vortices at distance r from the center of the trap (see Fig. 1). As stated before, the Magnus force for the locked state is nonzero because the superfluids are rotating at different rates. Balancing the forces gives F α mag = F α rstr . Because the restoring force acting on either species has the same magnitude, we obtain F 1 mag = F 2 mag . Noting that v α = Ω α r and v v = Ω v r (we are assuming that r is much larger than the distance between the two vortices) gives n 1 0 (Ω v − Ω 1 ) = n 2 0 (Ω 2 − Ω v ). This, along with the condition m 1 Ω 1 = m 2 Ω 2 (from ρ 1 v = ρ 2 v ) gives the following relation between the angular velocities: (11) Note that unlike the restoring force, the Magnus force grows linearly with distance from the center of the condensate. Thus, at some critical distance r c , the pairs of vortices invariably become unbound from each other. for radii r > r c the unlocked phase applies, and after a short healing region, the two condensates rotate at the same frequency. An expression for r c can be obtained by equating the Magnus force with the maximum possible value for the restoring force: Note that (12) diverges when the masses are equal. In addition, the bound pairs of vortices are pulled further apart at increasing distances from the center of the condensate. The interspecies vortex separation x v satisfies which is valid for r < r c . This introduces a small correction to the vortex density and creates a small shear in the motion of the two condensates.
For the vortex-locked state to be stable up to the critical radius r c , the superfluid velocities in the rotating frame of the vortex lattice, |v α SF − v v |, must not exceed the critical velocity of the superfluid. Otherwise, it would be possible to create elementary excitations from the flow of the superfluid around the vortices. For the two coupled superfluids, the Bogoliubov elementary excitations are given by where E α k = (ε α k ) 2 + 2g α n α 0 ε α k and ε α k =h 2 k 2 2mα . It is then straightforward to compute the critical velocity v c = min k ( Ω k hk ); one obtains v c = min α g α n α The superfluid velocity of species 1 or 2 in the vortex lattice frame evaluated at r c (where it is maximal) is given by The condition |v α SF − v v | < v c must be checked so that the vortex-locked state is stable against the creation of elementary excitations. For instance, it can be shown that the system is stable against creating such elementary excitations if the conditions 1 10 ≤ n 0 1 n 0 2 ≤ 10 and 1 10 ≤ m1 m2 ≤ 10 are satisfied. Numerical Simulations. Now that the expected types of phases have been discussed, we minimize the total energy E = E 1 + E 2 + E 12 as a function of the vortex positions given by Eqs. (6) and (8). The ability to compute analytical expressions for the gradients of the energy γ=28 γ=17 γ=5 and the above procedure is repeated until it converges.
To simplify the analysis, we restrict our attention to the case in which the equilibrium densities and healing lengths of the two condensate components are equal: n 0 ≡ n 1 0 = n 2 0 , ξ ≡ ξ 1 = ξ 2 . Motivated by the example of a 133 Cs-87 Rb condensate [16], we fix the mass ratio to be m 1 /m 2 = 1.5. The vortex interaction strength is parametrized by γ = |g12|n0 hΩv , which we vary while keeping the quantitiesh π Ωv mα 1 πξ 2 for α = 1, 2 fixed. (The total energy has been scaled byhΩ v πn 0 ξ 2 .) We set the ratio of the "average" vortex lattice constant to the coherence length a lat /ξ=10 (which is consistent with typical experiments). We define a lat by m1+m2 , and consider a system with 43 vortices of each species. The results for such a calculation are shown in Fig. 2. We also plot our prediction for the critical radius at which the vortex pairs become unbound [see Eq. (12)], which for equal densities is This prediction agrees quite well with our numerical results. Experimental realization. A very promising candidate for the realization of these locked states is a 133 Cs-87 Rb condensate mixture [16], which has a mass ratio of about 1.5. One has exquisite control over the self-scattering length of Cesium [17], and a Cs-Rb mixture is also expected to exhibit interspecies Feshbach resonances. This allows one to control the interaction g 12 over a wide range; such interspecies resonances have recently been identified for Li-Na [18] and Rb-K [19] mixtures.
Conclusions. In this paper, we described a novel state of rotating interacting condensates with unequal masses in the Thomas-Fermi regime. The possibility of locking the two individual vortex lattices yields a remarkable demonstration of the nonintuitive behavior of superfluids: the two gasses, rather than equilibrating to the same speed, prefer to move at different angular velocities that are inversely proportional to their masses. The vortexlocking of the different-mass condensates is also an example of synchronization: a phenomenon that is ubiquitous in physics, biology, and other fields [20]. Already in single-mass mixtures, a rich variety of vortex dynamics arises from the extra degrees-of-freedom, resulting in such effects as the formation of square vortex lattice, as well as topologically nontrivial defects such as skyrmions or hedgehogs [13,21]. To investigate these effects in the different-mass mixtures, as well as to better establish the locked state we proposed in this manuscript, this system must be numerically investigated by solving the appropriate dynamical Gross-Pitaevskii equations. Such numerical investigation will also allow one to find the preferred lattice geometry of the locked vortex-lattice. Other directions for future study involve dynamical aspects such as Tkachenko modes [22]) of the locked state, as well as whether a similar state survives in the Landau regime of large vortex density.
ence Foundation under Grant No. PHY05-51164 (RB, GR, HPB). We also acknowledge useful discussions with Simon Cornish, Michael Cross, Peter Engels, and Erich Mueller.

APPENDIX
In this appendix we discuss the change in the superfluid density as a result of a vortex lattice, and its effect on the validity of approximation (7) and the resulting expression for the interspecies vortex attraction, Eq. (8) and (5).
First consider a single species which has the energy functional We write ψ = f e iθ where f is real and θ contains information about the positions of the vortices as θ = i θ i . By varying f we obtain the equation determining the superfluid density which minimizes E for the particular vortex configuration. First let us consider a single vortex taken to be at the origin ∇θ =ẑ ×r r . The long-distance behavior r ≫ ξ is obtained from the Thomas-Fermi approximation (neglecting the ∇ 2 f in Eq. 21) and we obtain where we have neglected the contribution from the trapping potential. This implies that the suppression of the density is due to the kinetic energy in the supercurrent, which counters the condensation energy of the BEC. Eq. (22) seems to imply that at large distances from a vortex core, the interspecies vortex-vortex interaction will include a persistent power-law component, and die off only as 1/r 2 , rather than as an exponential. But the observation that the power law decay reflects the currentinduced superfluid suppression, allows as to ignore the power-law decay in a many-vortex situation, with the argument as follows. Let us consider a vortex lattice, and evaluate f at a position which is several coherence lengths away from any vortex. This allows us to invoke the continuum approximation and write where ρ v is the density of the vortices. Each individual contribution alone in this sum would yield a ∝ 1 r 2 dependence in the density corrections, but the vector sum of the velocities of all vortices squared is not a simple sum of the 1/r 2 corrections. Inserting this into the equation for the density profile one finds: Thus the combined vortex effect re-normalizes the trapping potential and does not need to be explicitly taken into account. To get an estimate for the magnitude of such a renormalization, one can compare this term with the chemical potential where a lat is the vortex lattice constant which is small for typical experiments.
For a two-component BEC, the situation is similar for vortices which are many coherence lengths away from each other. On the other hand, when the cores of the different types of overlap, their interaction needs to be explicitly calculated, and the continuum approximation cannot be used. This is the case for paired-vortex configurations. Since the combined effect of far-away vortices on a locked pair is small, the locking depends only on the short distance density profile taken. Had we used the step-function potential interaction between two vortices of Ref. 15 the results would only differ from the Gaussian depletion Eq. 7 by small quantitative amounts.