Generating vortex rings in Bose-Einstein condensates in the line-source approximation
M. Guilleumas1 , D.M. Jezek2 , R. Mayol1 , M. Pi1 , M. Barranco1
arXiv:cond-mat/0202312v1 19 Feb 2002
d’Estructura i Constituents de la Mat` eria, Facultat de F? ?sica, Universitat de Barcelona, E-08028 Barcelona, Spain 2 Departamento de F? ?sica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, RA-1428 Buenos Aires, and Consejo Nacional de Investigaciones Cient? ??cas y T? ecnicas, Argentina (February 1, 2008)
We present a numerical method for generating vortex rings in BoseEinstein condensates con?ned in axially symmetric traps. The vortex ring is generated using the line-source approximation for the vorticity, i.e., the rotational of the super?uid velocity ?eld is di?erent from zero only on a circumference of given radius located on a plane perpendicular to the symmetry axis and coaxial with it. The particle density is obtained by solving a modi?ed Gross-Pitaevskii equation that incorporates the e?ect of the velocity ?eld. We discuss the appearance of density pro?les, the vortex core structure and the vortex nucleation energy, i.e., the energy di?erence between vortical and ground-state con?gurations. This is used to present a qualitative description of the vortex dynamics. 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.40.Vs, 67.40.Db
Typeset using REVTEX 1
Since 1999, when vortex lines in a trapped Bose-Einstein condensate (BEC) were ?rst experimentally obtained , their study has received great experimental and theoretical interest as it constitutes a clear signature of super?uidity e?ects in these con?ned systems. Some remarkable experimental achievements are, among others, the study of the dynamics of single vortex lines  and the formation of small  and large  vortex arrays. A review of the research done in this ?eld is presented in Ref. . It is worth to note that the experimental situation in BEC regarding the formation and detection of vortices is at variance with that in helium II droplets. For these drops, a paradigm of super?uid ?nite systems, the problem of nucleating vortices and their stability has only been addressed very recently from the theoretical point of view (see Refs. [6–8]), and their experimental detection is still an open question. Vortex rings are vortices whose core is a closed loop with quantized circulation around it . They are complex topological structures that have attracted and will continue to attract some experimental and theoretical interest. Based on numerical simulations, di?erent methods have been proposed to generate vortex rings in BEC. As in bulk liquid helium, vortex rings may be produced introducing an impurity in the condensate with a de?nite velocity, whose displacement causes a vortex ring . Another mechanism  consists in using dynamical instabilities in the condensate to cause dark solitons to decay into vortex rings. A well controlled method to produce vortex rings by electromagnetically induced atomic transitions in two-component condensates has also been put forward . The method proposed in Ref.  has been successfully applied in Ref.  to generate vortex rings experimentally. Rather than proposing a method that could be implemented to create a vortex ring experimentally, our aim here is to set up a numerical method simple yet accurate enough to generate quantized vortex rings with a de?nite radius R in one-component condensates. We restrict our study to vortex states such that the divergence of their velocity is vanishingly small. This assumption yields analytical expressions for the velocity ?eld, and provides a fair approximation to the super?uid ?ow around the vortex core for rings in the bulk of large condensates (Thomas-Fermi limit). We have considered large condensates at zero temperature and therefore, dissipation has not been taken into account. They are axially symmetric about the z axis, and have the z = 0 plane as symmetry plane, and may host a vortex ring of radius R, coaxial with the symmetry axis of the trap, and placed on a plane at a distance Z ≥ 0 from the symmetry plane. Although generalization to situations with more than one such vortex rings is straightforward, we have not considered this possibility. This paper is organized as follows. In Sect. II we present the method used to generate vortices. The way to obtain the velocity ?eld of a vortex ring is described in Sect. III, and the explicit expressions of the velocity components are given in an Appendix. Sect. IV is devoted to the analysis of the particle density pro?les and vortex nucleation energies, which allows a qualitative description of their dynamics. We also present results obtained from a completely di?erent, more involved method we have set up to generate ring vortices that permits to test the approximation of zero divergence of the velocity. Finally, a brief summary of the results is presented in Sect. V.
II. ENERGY FUNCTIONAL
We consider a weakly interacting Bose-condensed gas con?ned in a harmonic trap Vext (r) at zero temperature. In the Gross-Pitaevskii (GP) theory, the ground state (g.s.) energy of the condensate is given by the functional  E[ Ψ ] = h ?2 g dr |?Ψ|2 + Vext (r)|Ψ|2 + |Ψ|4 2m 2 , (1)
where Ψ(r) is the condensate wave function. The ?rst term in Eq. (1) is the kinetic energy of the condensate, the second term is the harmonic oscillator energy arising from the trapping potential, and the third term is the mean-?eld interaction energy. The coupling constant is g = 4πh ? 2 a/m, where a is the s-wave scattering length, and m is the atomic mass. The number of atoms in the condensate is dr|Ψ|2 = N . The g.s. wave function is determined by solving the GP equation obtained minimizing the energy functional. The wave function Ψ can be written in terms of the particle density ρ(r) = |Ψ(r)|2 and phase S (r) as Ψ(r) = ρ(r) exp[iS (r)] , (2)
and the super?uid velocity is given by v = (? h/m)?S . In this work we will use the equivalent quantum hydrodynamic description of the condensate in terms of the density and the super?uid velocity , since it allows a straightforward generalization of the energy functional Eq. (1) to include vortex states. Using Eq. (2) it follows that
v E [ Ψ ] = E0 [ ρ ] + Ekin [ ρ, v ] .
The ?rst term is only density-dependent: E0 [ ρ ] = dr h ? g √ |? ρ|2 + Vext (r)ρ + ρ2 2m 2 . (4)
The ?rst term in E0 [ρ] is the quantum kinetic energy. The second term of the energy functional Eq. (3) corresponds to the kinetic energy associated to the ?ow velocity v, and is given by 1 v Ekin [ ρ, v ] = m 2 dr ρ(r) v2(r) . (5)
The g.s. wave function in the absence of vortices has a spatially constant phase and therefore zero velocity. It can be obtained by minimizing E0 [ ρ ]. We consider condensates with 2 2 2 positive scattering length and axially symmetric traps Vext (r) = m[ω⊥ (x2 + y 2 ) + ωz z ]/2 2 2 2 2 [= m(ω⊥ r + ωz z )/2 in cylindrical (r, z ) coordinates], with di?erent values of the asymmetry parameter λ = ωz /ω⊥ . The trap harmonic frequency ω⊥ provides a length scale for the system, a⊥ = (? h/mω⊥ )1/2 . We will use a⊥ , h ? ω⊥ , and N/a3 ⊥ as units of length, energy, and density, respectively. If we consider large condensates in which the Thomas-Fermi (TF) approximation holds , the quantum kinetic energy can be neglected compared to the interaction energy in 3
Eq. (4), and the g.s. density of the condensate in the absence of vortices is given by 2 2 ρ0 (r, z ) = ?(1 ? r 2 /RTF ? z 2 /ZTF )/g in the region where this expression is positive, and zero elsewhere. The TF extents of the condensate in the radial and axial directions are 2 1/2 RTF = (2?/mω⊥ ) and ZTF = RTF /λ. The chemical potential ? is ?xed by normalization ?=h ? ω⊥ (15λaN/a⊥ )2/5 /2. We recall that the validity of the Thomas-Fermi approximation is guaranteed if Na/a⊥ ? 1. We turn now our attention to the case of condensates with vortex states characterized by a given irrotational velocity ?eld associated with a non-vanishing quantized circulation. The total energy of the system is given by the energy functional Eq. (3). If the velocity ?eld is known, for a given number of particles the density pro?le can be obtained minimizing Eq. (3). This yields the equation ? h ? 2 ?2 1 + Vext (r) + g | ψ (r) |2 + mv2 ψ (r) = ? ψ (r) , 2m 2 (6)
√ where ψ = ρ is the modulus of the complex wave function Eq. (2). This is the GP equation expressed in terms of the hydrodynamic variables. In the presence of a quantized vortex, the density of the system drops to zero at its core, whose size is characterized by a healing length ξ . For very large condensates the healing length can be approximated by ξ = (8πρ0 a)?1/2 , where ρ0 is the density of the condensate before creating the vortex. For a centered vortex line in the TF approximation ρ0 = ?/g , and the corresponding healing length ξ0 is a⊥ ξ0 = RTF RTF
with ξ0 ? a⊥ ? RT F . In this approximation a local healing length can be de?ned  as ξ (r, z ) = ξ0 / 1 ? (r/RTF )2 ? (z/ZTF )2 . Note that the size of the core is larger for a vortex in the low-density region.
III. SUPERFLUID VELOCITY FIELD
The velocity ?eld around a straight vortex line has an analytical expression when the vortex is along the symmetry axis [5,16]. Approximate analytical expressions can be found in the case of vortex lines o? the symmetry axis in large and very elongated, quasi-two dimensional condensates (see Ref.  and references therein). However, vortex lines generally bend in three-dimensional condensates [5,18] which renders impracticable an analytical treatment of the velocity ?eld. In the case of quantized vortex rings we proceed as Schwarz and Jang have done in the case of helium II  (see also Ref. ). For a vortex ring characterized by the values of (R, Z ) already de?ned, circulation number n = 1, 2, . . . and quantum circulation k0 = n h/m, we write the vorticity in the line-source approximation: ? , ω = k 0 δ (r ? R ) δ (z ? Z ) φ (8)
? is the unit vector in the azimuthal where (r, z, φ) are the cylindrical coordinates, and φ ? = (? sin φ, cos φ, 0). The super?uid velocity ?eld that arises from this distributed direction φ 4
vorticity ful?lls ω = ? × v. Hence, the velocity ?eld around the vortex is irrotational except on the vorticity line, where the density of the condensate is zero. If ? · v = 0 to a good approximation , a velocity vector potential A(r) can be introduced such that v = ?× A. If the vorticity is speci?ed, A(r) is determined by the equation ω = ? × (? × A ) , ?, with whose integral solution is  A(r) = A0 (r, z )φ A0 (r, z ) = k0 R 4π
cos φ ′ dφ ′ r 2 + R2 ? 2 rR cos φ ′ + (z ? Z )2
The radial and z -components of the velocity are obtained as vr = ? vz = ?A0 ?z (11) 1 ? (rA0 ) , r ?r
whereas the azimuthal component of the velocity ?eld around the vortex ring Eq. (8) is zero. We give in the Appendix the general expressions of A0 , vr and vz written in terms of hypergeometric functions . For a con?ned system, the existence of a boundary and the fact that the density is inhomogeneous may have an e?ect on the actual velocity ?eld . In our case, since the condensate extends up to ‘in?nite distances’, there is no need to introduce any image vortex to ensure that there is no particle ?ow across the boundary due to the super?uid motion. Moreover, it has been shown  that in the TF limit the corrections to the velocity due to density inhomogeneities can be safely neglected. Thus, we approximate the velocity ?eld by Eqs. (10) and (11) (see Appendix), but have restricted ourselves to study vortex ring con?gurations (R, Z ) whose vorticity line is inside the domain where the density is positive in the TF approximation, i.e., only vortex rings (R, Z ) satisfying the condition (R/RTF )2 + (Z/ZTF )2 < 1 are considered. We call the (R, Z ) line de?ned by the condition (R/RTF )2 + (Z/ZTF)2 = 1 the TF boundary. Once the the super?uid velocity ?eld has been ?xed, the density pro?le of the condensate is obtained solving Eq. (6).
We present numerical results for large condensates in the Thomas-Fermi limit hosting a quantized vortex ring with circulation number equal to one (nucleating vortex rings with n > 1 is energetically less favorable [9,24]). For a singly quantized vortex ring con?guration (R, Z ) with vorticity given by Eq. (8) and circulation number n = 1, we have computed the velocity ?eld, Eq. (11), and have obtained the density pro?le of the con?ned condensate by solving Eq. (6) in cilyndrical coordinates using the imaginary time method . Note that although we are in the TF limit, we do solve the complete GP equation to obtain the density pro?les. 5
Figure 1 shows density pro?les in the z = 0 plane as function of r for the experimental parameters of Ref. , that is, N = 3 × 105 atoms of 87 Rb (scattering length a = 5.82 × 10?9 cm) in a spherical trap with ω⊥ /2π = 7.8 Hz. We have plotted two con?gurations with a vortex ring located in Z = 0 having a radius R = 2 a⊥ (dotted line), and R = 4 a⊥ (dot-dashed line), respectively. The density pro?le of the condensate without vortex is also shown (solid line). The density is zero on the vorticity line, as seen in the density pro?les. As expected, the size of the core is of the order of the healing length. Indeed, for this condensate one has RTF ? 5.8 a⊥ and ξ0 ? 0.17 a⊥ . Assuming that the vortex diameter is twice the local healing length, and using the TF expressions we get as core diameters the values 2 ×ξ (2, 0) = 0.37 a⊥ , and 2 × ξ (4, 0) = 0.47 a⊥ . One can see from Fig. 1 that the core sizes are in agreement with these estimates. When the radius of the vortex ring increases, the core lies at a lower density region and therefore its size increases. To study the vortex ring energetics we have chosen a larger condensate made of N = 106 atoms of 87 Rb con?ned in an axially symmetric trap with axial frequency ωz /2π = 220 Hz and √three di?erent geometries, namely spherical, disk-shaped and cigar-shaped, with λ = 1, 8 and 0.2, respectively. For this condensate, we show in Fig. 2 several equidensity lines (arbitrary values) in the y = 0 plane. The top panel corresponds to the cigar-shaped trap (λ = 0.2, RTF = 8.7 a⊥ , ZTF = 43.3 a⊥ ), the middle panel to the spherically symmetric √ trap (λ = 1, RTF = ZTF = 10.2 a⊥ ), and the bottom panel to the disk-shaped trap (λ = 8, RTF = 45.5 a⊥ , ZTF = 11.3 a⊥ ). The x and z axes are in units of a⊥ . For all geometries, the vortex ring con?guration is (R = 3.1 a⊥ , Z = 0). The intersection between the vorticity line and the y = 0 plane appears as two dark dots, indicating the steep density depression around the vortex core. The presence of the vorticity causes a drastic distortion of the density with respect to the g.s. pro?le. Yet, for large condensates it is a rather local e?ect, as can be seen from Fig. 2, and also from Fig. 1 for the smaller condensate. The nucleation of a vortex has an energy cost, since the energy of a condensate with a vortex is always larger than the energy of the condensate without it, EGS . The vortex nucleation energy, E ? EGS , corresponding to the condensate of the spherical trap in Fig. 2 is plotted in Fig. 3 as a function of Z for di?erent values of R. Fixed R, the nucleation energy is maximum when the ring is in the z = 0 plane, and it decreases as Z increases, since displacing the ring outwards implies that the super?uid ?ow a?ects less atoms in the condensate. Note that to have the vorticity line within the TF boundary, the R-curves end at di?erent values of Z/ZTF given by 1 ? (R/RTF )2 . The nucleation energy does not vanish when the ring is located on the TF boundary as there still exists a super?uid velocity ?eld which produces some e?ect. Since we treat the super?uid velocity as an external ?eld, the nucleation energy of very super?cial vortex con?gurations may be somewhat overestimated. This drawback is apparent beyond the TF boundary, and shows up as a long-tailed nucleation energy as a function of either Z/ZT F or R/RT F , although it goes eventually to zero. The nucleation energy is plotted in Fig. 4 as a function of R for di?erent Z values. The smaller ring we have considered has R/RTF ? 0.1. The middle panel corresponds to the same spherical condensate as in Fig. 3. Besides, the top panel corresponds to the cigar-shaped trap, and the bottom panel to the disk-shaped trap. The top curve in all panels corresponds to vortex rings with Z = 0, and for a given trap geometry (?xed panel), the higher the curve, the lower the Z . As in Fig. 3, we have only described vortex ring con?gurations inside the 6
TF boundary. For this reason, ?xed a Z value the corresponding nucleation energy curve stops at the value R/RTF = 1 ? (Z/ZTF)2 . Figure 4 shows that for each asymmetry parameter λ, given a Z there exists a vortex con?guration corresponding to a certain radius R that maximizes the nucleation energy. This radius increases as Z decreases. In the (R, Z ) two-dimensional con?guration space, the nucleation energy E (R, Z ) ? EGS has only one maximum at (Req /RTF ? 0.6, Z = 0). It is in agreement with the results obtained in Ref.  for a spherically symmetric trap using a simpli?ed model of a vortex ring. This maximum corresponds to the only stationary, although unstable, vortex ring con?guration (equivalent to a vortex line nucleated along the z -axis), whose location is almost independent of λ but its value depends on it, decreasing as λ does. This means that the energy to nucleate a vortex ring characterized by the values of (R, Z ) is lower the more elongated the trap is. For each asymmetry parameter, the ‘iso-nucleation energy’ lines in the (R, Z ) con?guration space corresponding to a value not too far from the maximum value are closed curves around the maximum, and the allowed ring radii are in the range de?ned by the intersection of a straight line representing the given nucleation energy and the Z = 0 curve in Fig. 4. If we relax the condition that the vorticity line must be inside the TF boundary, this should be always like this . Otherwise, the isoenergy lines eventually do not close. Yet, the allowed ring radii are determined in a similar way. We have considered the vortex ring as if it were a static object. However, in a trapped condensate a vortex ring will move with a velocity that results from the interplay between the e?ect caused by the inhomogeneity of the condensate and self-induced e?ects arising from its own local curvature [5,9,26]. A qualitative analysis of the dynamics of coaxial vortex rings can be carried out from Figs. 3 and 4, as proposed in Ref. , and tested in a real time dynamics for two-component condensates [12,26]. It can be seen that a coaxial vortex ring oscillates inside the condensate, moving towards the surface and instead of being annihilated at the boundary, the ring generates a background ?ow in the rest of the cloud that draws it back towards the center. The ring dynamics in a trapped condensate results in an oscillation along the symmetry axis of the trap. The radius of the vortex ring grows or shrinks depending on its Z -position, in such a way that the ring motion almost follows a trajectory with constant energy in the (R, Z ) con?guration space. Finally, we have used a completely di?erent method, numerically more involved, to generate quantized vortex rings. It does not make any assumption on the super?uid velocity ?eld, nor on the radius or position of the vortex ring either, and can be used to test the results we have obtained based on the validity of the approximation ? · v = 0. We have proceeded by analogy with the case of a vortex line, and have added a constraint to the energy functional that forces to form a toroidal-like hole in the condensate with a ?ux of atoms around its core with zero azimuthal velocity component. We have used as ? axis, i.e., φ ?·L ≡ Lφ = ih constraining operator the angular momentum about the φ ? (r?z ?z?r ) . Introducing the functional E[ Ψ ] = dr h ?2 g |?Ψ|2 + Vext (r)|Ψ|2 + |Ψ|4 + ?φ 2m 2 7 d r Ψ?Lφ Ψ , (12)
where the Lagrange multiplier ?φ can be understood as the local angular velocity about the azimuthal direction, we have obtained the associated GP equation. For a given ?φ , we have solved the partial di?erential equations obeyed by the real and imaginary part of the condensate wave function, which are coupled by the presence of the constraint. These equations have been discretized on a (r, z ) mesh using seven-point formulas to represent the di?erential operators, and have been solved using again the imaginary time method . Taking as an example ?φ = 0.2 ω⊥ , we have obtained a coaxial vortex ring con?guration with its vorticity located on the z = 0 plane whose equidensity lines cannot be distinguished from these represented in the middle panel of Fig. 2 at the scale of this ?gure (whose parameters were chosen before to allow us to make this statement). Figure 5 presents a comparison between the two methods. It shows the particle density pro?le as a function of r at z = 0 for the spherical condensate of Fig. 3. The solid line is the con?guration that minimizes Eq. (12). The density goes to zero in R = 3.1 a⊥ , and the vortex ring has an energy equal to 37.164 h ? ω⊥ . The dashed line is the con?guration obtained by the method of Sect. 2 with the ring vorticity Eq. (8) placed in (R = 3.1 a⊥ , Z = 0) which has an energy of 37.167 h ? ω⊥ . Thus, the energy and density pro?le of both con?gurations are almost identical, and the velocity ?elds are also in good agreement. We have also checked that minimizing the functional Eq. (12) yields quantized ring vortices with n = 1. To this end, we have proceeded to a direct numerical integration of the circulation v · dl using several arbitrary closed paths in the (r, z ) ‘plane’. We have found that the circulation for paths enclosing the vortex core is h/m with a good accuracy, and is zero otherwise. The velocity ?eld has been calculated from the wave function recalling that the current ?eld is j(r) = ρ v = (? h/2mi)[Ψ? ?Ψ ? (?Ψ? )Ψ]. The current ?eld j(r) in the y = 0 plane is plotted in arbitrary units in Fig. 6. The axes are in units of a⊥ . It can be clearly distinguished in this ?gure the position of the vortex ring core at Z = 0 and R ? 3.1 a⊥ .
We have developed a method to study vortex rings hosted in con?ned condensates in the TF limit. It provides analytical expressions for the super?uid velocity ?eld that satisfy the irrotational condition and the quantization of the circulation of the super?uid ?ow. The velocity around the vortex core is introduced as an external ?eld in the Gross-Pitaevskii energy functional, which is minimized to obtain the wave function of the condensate. Using this method, we have computed the density pro?le and nucleation energy of vortex rings in one-component condensates. We have found that the size of the vortex ring core is of the order of the healing length as in the case of vortex lines. The presence of the vortex causes a sizeable distortion of the density, as atoms are pushed o? the core, but the e?ect is rather local. The analysis of the dependence of the nucleation energy on the position and radius of the ring allows one to conclude that a Hamiltonian dynamics will lead to oscillations of the vortex ring along the symmetry axis of the condensate, increasing and decreasing its radius in accordance with previous predictions. A dissipative dynamics would cause the vortex to decay, for example disappearing across the border of the condensate, as occurs for vortex lines . 8
The method provides a handleable way to generate vortex rings in con?ned condensates that can be used as a starting point to study the dynamics. A detailed calculation of the vortex ring dynamics in a con?ned condensate, either Hamiltonian or dissipative, is beyond the scope of the present work and will be developed elsewhere.
This work has been performed under Grants No. PB98-1247 from DGESIC, Spain, and No. 2000SGR-00024 from Generalitat de Catalunya. D.M.J. acknowledges the CONICET (Argentina) and the Generalitat de Catalunya ACI program for ?nancial support.
In this Appendix we give the expressions and method we have used to obtain the velocity ?eld. From Eqs. (10-11) one has: vr (r, z ) = π [(r + R)2 + (z ? Z )2 ]3/2 k 0 R (z ? Z )
2 cos2 φ ? 1 dφ , (1 ? s cos2 φ)3/2
where we have de?ned s ≡ 4rR/ [(r + R)2 + (z ? Z )2 ]. The integral Eq. (A1) is written in terms of hypergeometric functions  yielding: vr (r, z ) = Similarly, vz (r, z ) = [F (3/2, 1/2; 2; s) ? F (1/2, 1/2; 1; s)] 2 r [(r + R)2 + (z ? Z )2 ]1/2 k0 R [(r + R) F (1/2, 3/2; 1; s) ? (r + 2R) F (3/2, 3/2; 2; s) + 2 [(r + R)2 + (z ? Z )2 ]3/2 3R + F (5/2, 3/2; 3; s)] (A3) 2 k0 R 2 [(r + R )2 + (z ? Z )2 ]1/2 k0 R 2 [(r + R)2 + (z ? Z )2 ]3/2 k 0 R (z ? Z ) [F (3/2, 3/2; 2; s) ? F (1/2, 3/2; 1; s)] . (A2)
and A0 (r, z ) = [F (3/2, 1/2; 2; s) ? F (1/2, 1/2; 1; s)] . (A4)
The hypergeometric functions entering Eqs. (A2-A4) can be written in terms of the complete elliptic integrals E and K  as follows: F (1/2, 1/2; 1; s) = 2 √ K ( s) π (A5)
√ 2 E ( s) F (1/2, 3/2; 1; s) = F (3/2, 1/2; 1; s) = π 1?s 9
F (3/2, 1/2; 2; s) = F (1/2, 3/2; 2; s) =
√ √ 4 [K( s) ? E( s)] πs
√ √ 4 E ( s) ? K ( s) F (3/2, 3/2; 2; s) = πs 1 ? s F (5/2, 3/2; 3; s) = √ 16 2 ? s √ E ( s ) ? 2 K ( s) . 3πs2 1 ? s
To evaluate E and K we have used polynomial approximations .
 M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, C.E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 83, 2498 (1999).  B. P. Anderson, P. C. Haljan, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 85, 2857 (2000).  K. W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000).  J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ketterle. Science 292, 476 (2001).  A. L. Fetter and A. A. Svidzinsky, J. Phys.: Condens. Matter 13, R135 (2001).  G. H. Bauer, R. J. Donnelly, and W. F. Vinen, J. Low Temp. Phys. 98, 47 (1995).  F. Dalfovo, R. Mayol, M. Pi, and M. Barranco, Phys. Rev. Lett. 85, 1028 (2000).  R. Mayol, M. Pi, M. Barranco, and F. Dalfovo, Phys. Rev. Lett. 87, 145301 (2001).  R. J. Donnelly, Quantized Vortices in Helium II, (Cambridge University Press, Cambridge, 1991).  B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. A 60, 4882 (1999).  D. L. Feder, M. S. Pindzola, L. A. Collins, B. I. Schneider, and C. W. Clark, Phys. Rev. A 62, 053606 (2000).  J. Ruostekoski and J. R. Anglin, Phys. Rev. Lett. 86, 3934 (2001).  B. P. Anderson, P. C. Haljan, C. A. Regal, D. L. Feder, L. A. Collins, C. W. Clark, and E. A. Cornell, Phys. Rev. Lett. 86, 2926 (2001).  F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999).  E. Lundh and P. Ao, Phys. Rev. A 61, 063612 (2000).  F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996).  M. Guilleumas and R. Graham, Phys. Rev. A 64, 033607 (2001).  J. J. Garc? ?a-Ripoll and V. M. P? erez-Garc? ?a, Phys. Rev. A 63, 041603(R) (2001).  K. K. Schwarz and P. S. Jang, Phys. Rev. A 8, 3199 (1973).  D. M. Jezek, M. Pi, M. Barranco, R. J. Lombard, and M. Guilleumas, J. Low Temp. Phys. 112, 303 (1998).  Note that since the density is not uniform, the assumption ? · v ≈ 0 does not imply that the system is incompressible.  I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, New York, 1980).  J. R. Anglin, e-print cond-mat/0110389.  T. Winiecki, J. F. McCann, and C. S. Adams, Europhys. Lett., 48, 476 (1999).  M. Pi, A. Emperador, M. Barranco, and F. Garcias, Phys. Rev. B 63, 115316 (2001).  B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. A 61, 013604 (1999). ?·L+L·φ ?)/2.  Actually, we have used the hermitian form Lφ = (φ  P. O. Fedichev and G. V. Shlyapnikov Phys. Rev. A 60, R1779 (1999).  M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Dover Pu., New York, 1970).
ρ (r,z=0) (10 N/a⊥)
10 5 0 0 1 2 3 4 5 6 7
FIG. 1. Density pro?le as a function of r at z = 0 for a condensate in a spherical trap with the experimental parameters of Ref.  (aN/a⊥ = 440) hosting a singly quantized vortex ring at Z = 0. The dotted line corresponds to a vortex ring with R = 2 a⊥ , and the dot-dashed line to a vortex ring with R = 4 a⊥ . The solid line is the g.s. density pro?le in the absence of vortices. The densities have been normalized to one.
9 6 3 0 -3 -6
6 3 0 -3 -6 9 6 3 0 -3 -6 -9 -9 -6 -3 0 3 6 9
FIG. 2. Vortex ring state equidensity lines, in arbitrary units, plotted in the y = 0 plane (axes are in units of a⊥ ). The ring is at (R = 3.1 a⊥ , Z = 0) in a conden6 87 sate with N = 10 atoms of Rb, con?ned in an axially symmetric trap with axial frequency ωz /2π = 220 Hz and di?erent trap geometries. The top panel corresponds to a cigar-shaped trap (λ = 0.2, RTF = 8.7 a⊥ , ZTF = 43.3 a⊥ ), the middle panel to a spherically √ symmetric trap (λ = 1, RTF = 10.2 a⊥ ), and the bottom panel to a disk-shaped trap (λ = 8, RTF = 45.5 a⊥ , ZTF = 11.3 a⊥ ).
_ E - EGS ( hω⊥)
R=2 R=10 R=1 R=9
Z / ZTF
FIG. 3. Nucleation energy of a vortex ring (R, Z ) as a function of Z for di?erent values of the radius. The labels of the curves are in units of a⊥ . The condensate has N = 106 atoms of 87 Rb and it is con?ned in an spherically symmetric trap with frequency ωz /2π = 220 Hz. The TF radius is RTF = ZTF = 10.2 a⊥ .
0.06 0.05 0.04 0.03 0.02 0.01
_ ( hω⊥)
0.15 0.10 0.05
E - EGS
0.30 0.20 0.10 0.00 0.0
R / RTF
FIG. 4. Nucleation energy of a vortex ring (R, Z ) as a function of R for di?erent values of Z . The condensate and trap geometries are these of Fig. 2. The top panel corresponds to the cigar-shaped trap, the middle panel to the spherically symmetric trap, and the bottom panel to the disk-shaped trap. From top to bottom, and in units of a⊥ , the curves correspond to Z = 0 to 40 in Z -steps of 5 (top panel); Z = 0 to 10 in Z -steps of 1 (middle panel); Z = 0, 1, 2, and 3 (bottom panel).
6 5 4 3
ρ (r,z=0) (10 N/a⊥)
2 1 0 0 1 2 3 4 5 6 7 8 9 10 11
FIG. 5. Density pro?le as function of r at z = 0 for the N = 106 condensate in a spherically symmetric trap with a singly quantized vortex ring. The solid line corresponds to the density pro?le that minimizes the functional Eq. (12) with constraint ?φ = 0.2 ω⊥ . The dashed line is the density obtained imposing a ring vorticity at (R = 3.1 a⊥ , Z = 0). The densities have been normalized to one.
2 1 0 -1 -2 0 1 2 3 4 5 6
FIG. 6. Current ?eld around the vortex core corresponding to the vortex ring state calculated taking ?φ = 0.2 ω⊥ in Eq. (12). The current is in arbitrary units and the axes are in units of a⊥ .