Vortices in the twodimensional dipolar Bose gas
Abstract
We present vortex solutions for the homogeneous twodimensional BoseEinstein condensate featuring dipolar atomic interactions, mapped out as a function of the dipolar interaction strength (relative to the contact interactions) and polarization direction. Stable vortex solutions arise in the regimes where the fully homogeneous system is stable to the phonon or roton instabilities. Close to these instabilities, the vortex profile differs significantly from that of a vortex in a nondipolar quantum gas, developing, for example, density ripples and an anisotropic core. Meanwhile, the vortex itself generates a mesoscopic dipolar potential which, at distance, scales as and has an angular dependence which mimics the microscopic dipolar interaction.
1 Introduction
A landmark experiment in quantum gas research was the achievement in 2005 of a dipolar BoseEinstein condensate (BEC) made of Cr atoms [1]. Other groups have subsequently used novel methods to make Cr condensates [2], as well as dipolar condensates composed of Dy [3] and Er [4]. The dipolar nature of these atoms arises from their sizeable magnetic dipole moment which for Cr is , for Er is , and for Dy is , where denotes the Bohr magneton. These numbers should be compared to the value found in the alkali atoms which are the standard workhorses used in most ultracold atomic experiments. While the dominant interatomic interactions in alkali atom condensates are the isotropic van der Waals interactions which fall off as , the effects of magnetic dipoledipole interactions, which are anisotropic and fall off as , have been seen in the BECs composed of the more strongly magnetic atoms listed above (note that the strength of the magnetic dipoledipole interaction goes as ). In particular, magnetostriction of a condensate has been observed [5] as well as dwave collapse and explosion [6]. Of particular relevance to this paper are theoretical studies of dipolar condensates which predict roton excitations [7, 8, 9], akin to those observed in Helium II. Although the experimental observation of rotons in dipolar BECs with static dipoledipole interactions remains illusive, their effects have been seen in BECs with laserinduced dipoledipole interactions [10]. For a comprehensive review of the field of dipolar BECs up to 2009 we refer the reader to the article [11].
Despite the enhanced dipoledipole interactions in Cr relative to alkali atoms, the dipolar interactions in Cr still play a secondary role compared to the van der Waals interactions unless something is done to suppress the latter. Fortuitously, this has proved possible using the technique of Feshbach resonance [12] which has allowed the creation of a purely dipolar condensate [13]. This method has also been applied to Er [4]. Polar molecules are another promising system for observing (electric) dipolar interactions due to their potentially very large electric dipole moments. These systems are now close to quantum degeneracy [14].
The quantumcoherent nature of BECs means that they support vortices characterised by a core of vanishing density and a 2quantized singularity of the quantum mechanical phase. Diverse vortex structures have been generated in BECs, from single vortices to vortex rings and vortex pairs, to vortex lattices and turbulent states [15, 16]. While vortices have not yet been observed in dipolar BECs, simulations indicate that vortex rings were formed in the dwave collapse experiment [6]. Theoretical work has also revealed that dipolar interactions modify the vortex properties, introducing density ripples about the core and an elliptical core [17, 18, 19, 20], as well as affecting the critical rotation frequency for vortex generation [21] and the associated dynamical instability [22, 23, 24]. The critical obstacle velocity for vortex nucleation is likewise modified [25] as well as the structure of vortex lattices [26].
Here we examine the single vortex solutions in an infinite twodimensional dipolar Bose gas where the dipoles are aligned by an external field, expanding on some aspects of our recent work [20]. We begin by establishing the stability regimes of the homogeneous 2D system. We then explore the vortex solutions as a function of dipolar strength and polarization angle, and end with an analysis of the dipolar potential generated by the vortex.
2 MeanField Model
We consider a dilute BEC composed of atoms of mass with dipole moment . These atomic dipoles are polarized in a common direction by an external magnetic field, with the polarization axis taken to lie at angle to the axis and in the plane. For the purposes of GrossPitaevskii meanfield theory the total interaction between two atoms at and can be represented by the potential [11, 27],
(1) 
The first term is a pseudopotential which models the van der Waals interactions through a contact potential with coupling strength , where is the swave scattering length. The second term is the bare dipoledipole interaction, where is the angle between the polarization direction and the interparticle vector . For magnetic dipoles the dipolar coupling strength is , where is the permeability of free space.
Working in the limit of zero temperature, we describe the condensate by the meanfield wavefunction , which specifies the atomic density . The wavefunction obeys the dipolar GrossPitaevskii equation,
(2) 
where is the external potential acting on the condensate and is the chemical potential. The van der Waals interactions lead to the purely local term in the above Hamiltonian, while the dipolar interactions are incorporated via the nonlocal potential given by
(3) 
where is the dipolar interaction potential, the second term in Eq. (1).
The reduction to an effectively 2D BEC, homogeneous in one plane, is performed by considering the external potential to be harmonic trapping in the axial direction only, , where is the corresponding trap frequency. The axial trapping is assumed to be so strong that becomes frozen into the axial ground harmonic state, a Gaussian of width . Formally, this “2D meanfield regime” is reached if , where is the average 2D density [28]. This scenario of an effectively 2D homogeneous dipolar BEC is illustrated in Fig. 1.
Within the 2D meanfield regime one can decompose the wavefunction according to , where is the 2D wavefunction in the plane. Integrating over then leads to the effective 2D dipolar GPE [25], equivalent to Eq. (2) with the following 3D quantities replaced by their 2D analogs: , , and . The reduced dipolar potential can be expressed as,
(4) 
where is the effective 2D dipolar interaction potential. Rather than evaluating directly through Eq. (4), it is more convenient to evaluate it using the convolution theorem as . Here the Fourier transform of has been shown elsewhere to be [25, 30],
(5) 
with , , and . The axial width of the condensate plays a more important role for 2D dipolar condensates than their nondipolar counterparts since it specifies how much of the outofplane dipolar interaction is experienced by the condensate.
The strength of the dipolar interactions can be parameterized via the ratio [11]
(6) 
We note that can be tuned between and via a Feshbach resonance while can be varied below its natural value, including to negative values, via magnetic field rotation [11]. Hence it is feasible to consider and both negative and positive .
Numerical solutions are obtained by propagating the 2D dipolar GPE in imaginary time using a RungaKutta algorithm on a grid extending over a region , where is the healing length (defined in the following section). At the boundary of the numerical box, the value of the wavefunction is fixed to the uniform value to mimic an infinite system (Dirichlet boundary conditions).
3 The Homogeneous 2D Dipolar Gas
In the absence of any potential in the plane the stable ground state solutions have uniform density, which we denote . This generates a uniform dipolar potential given by . The dipolar GPE then reduces to [25],
(7) 
where is the chemical potential of the homogeneous system. This relation is similar to the result for a homogeneous 2D nondipolar BEC , but where the dipolar interactions effectively introduce an additional contribution to the contact interaction . The angular dependence accounts for the “tilt” of the dipoles. Taking the dipoles to be conventional in the sense that , then for the dipoles lie sidebyside and experience the repulsive part of the dipoledipole interaction, while for , the dipoles lie headtotail in the plane and experience the attractive part of the interaction. For , the sign of the dipoledipole interaction becomes reversed, such that the arrangement of dipoles becomes attractive and the arrangement becomes repulsive. Note that, at the “magic angle” , the inplane dipolar interaction is zero (for either sign of ).
A natural length scale for the homogeneous system is the dipolar healing length, . In our analysis, we set the axial length scale of the condensate to be , which satisfies the criteria for the system to be twodimensional .
3.1 The phonon instability
Perturbations in the density and velocity of the condensate, with frequency and wavevector , satisfy the Bogoluibov excitation spectrum [11],
When is less than zero (i.e. the net local contract interactions are attractive), then in the phonon limit , the perturbations possess imaginary energies, signifying the unstable exponential growth of the perturbations over time. This is the phonon instability, familiar from the conventional condensate with attractive contact interactions [29].
From Eq. (), the instability condition means that [8, 30, 11, 33],
(8) 
Let us consider the implications for the cases of and separately, with the results summarized in Fig. 2(a) and (b) respectively. For the system undergoes the phonon instability (PI) when the inplane dipolar contact interactions are attractive (when either and , or and ) and of sufficient magnitude to overcome the repulsive vdW contact interactions (characterized by ). Dividing Eq. (8) by () and introducing we obtain the condition for the phonon instability to be . We define the threshold value of as (red lines in Fig. 2). Note the importance of the magic angle , at which the instability threshold diverges. It then follows that (for ) the phonon instability arises for when , and for when , as depicted in Fig. 2(a) (pink shaded regions).
Meanwhile, for the system will be stable to the phonon instability when inplane dipolar interactions are repulsive and of greater magnitude than the attractive vdW interactions. The instability criterion is then for and for , as depicted by in Fig. 2(b) (pink shaded regions).
3.2 The roton instability
In addition to the phonon instability, a dipolar BEC also experiences the roton instability. Dipolar interactions in a BEC can induce a roton dip in the excitation spectrum at finite momenta [7, 8, 18], reminiscent of rotons in superfluid Helium [35]. For certain parameters, the roton “softens” to zero energy, inducing a collapse instability at finite momentum, associated with the runaway alignment of dipoles in the headtotail configuration (for or sidebyside (for ). Close to the RI, ripples in the BEC density arise when the roton mode mixes into the ground state [17, 18].
We assess the regimes of roton stability/instability by numerically seeking solutions of the 2D dipolar GPE where the initial condition is a homogeneous wavefunction modulated with sinusoidal perturbations in and , . This perturbation facilitates the development of the roton instability, where it exists. Where stationary solutions are not obtained (and the system is in a regime which does not suffer the phonon instability) it is deemed to possess the roton instability. The roton instability regimes are depicted in the space of in Fig. 2(a) and (b) (blue shaded regions).
Again, we interpret the and cases separately. For the roton instability is induced by the attractive part of the dipolar interaction. For (dipoles polarized parallel to the axis) the condensate, confined to a narrow region about , cannot probe this attractive part of and no roton instability exists [30]. For increased polarization angle (while still in the regime ), the attractive part of can be experienced by the condensate and so a regime of roton instability does occur. This falls in the parameter regime and (i.e. ). Meanwhile, for , a region of roton instability occurs for negative , i.e. . Although the dipoles are approaching an headtotail configuration, since the dipoles repel in the plane and attract axially, with the latter attraction inducing the roton instability.
Now we consider the roton instability for . For the presence of attractive contact interactions can support a roton and an associated roton instability [34]. In contrast to the case, for the roton instability exists for all (barring the magic angle). For the roton instability occurs in a band of in the negative halfplane, i.e. . As for , the roton instability here is induced by the attractive contact interactions. As is made increasingly negative (corresponding to increasingly positive ) the roton instability is removed due to the stabilizing effect of the large, repulsive inplane dipolar interactions. For , a region of roton instability occurs for positive , i.e. negative . As for , case discussed above, the roton instability here is driven by the attractive axial component of the dipolar interactions.
Looking at the overall stability/instability of the system depicted in Fig. 2, it is noteworthy that, for , regimes of stability exist for all , while for the only stable regime lies for small and negative .
3.3 Role of
It follows from the analysis in Section 3.1 that the regimes of phonon stability/instability, and the threshold between them, is independent of the axial size of the condensate, . However, the roton instability is affected by . For a narrower (wider) condensate, the regimes of roton instability in the phase diagram shrink (expand). This is because the atoms experience less (more) of the outofplane component of . To illustrate the quantitative dependence, take for example . Denoting the threshold for the roton instability as the solutions are unstable for , where . For , and then , and , respectively.
4 Vortex Solutions
We generate singlycharged vortices at the origin of our otherwise homogeneous system by setting the phase of the wavefunction to be (corresponding to a singlycharged vortex) during imaginary time propagation of the dipolar GPE.
In the nondipolar system, a vortex has a size characterized by the conventional healing length , where . In the dipolar system, the vortex core size is of the order of the dipolar healing length . This is the motivation for expressing length in these units. However, it should be noted that this length scale is a function of the system parameters, and even diverges as .
We examine the structure of the vortex solutions for three representative polarization angles: (a) , (b) and (c) (each of which fall into distinct regimes of the stability diagram). Figures 3(a), (b) and (c) plot the vortex solutions for the respective polarization angle, as a function of . The lefthand side of the main plots shows the density profile along (with ) and the righthand side shows the density profile along (with ). Note that we find that the introduction of a vortex has no observable effect on the stability of the system, which maintains the stability diagram plotted in Fig. 2, and hence regions with no vortex solutions occur in the plot.

(a) : Since the dipoles are perpendicular to the plane, the dipolar potential, and hence density profile, are circularlysymmetric [2(a)]. For the vortex has the standard appearance familiar from nondipolar condensates [see left inset of Fig. 2(a)], consisting of a circularlysymmetric core of vanishing density of width [29].
For and for most of the parameter space the vortex solution is essentially identical to that for (due to our choice of as the length scale). Here the system behaves like a nondipolar system but with modified contact interactions, and the nonlocal effects of the dipolar interaction are insignificant. However, there exist two limiting regimes where the vortex solutions deviate significantly. Firstly, as one approaches from above, the vortex core becomes increasingly narrow (in units of ) and takes on a modified profile [see middle inset of Fig. 2(a)]. This is associated with the cancellation of explicit contact interactions in the system—the vdW interactions cancel the contact contribution from the dipolar interactions. Indeed, the total contact interaction disappears altogether at the PI threshold, ). Secondly, as the RI is approached from below, circularlysymmetric density ripples appear around the vortex [see middle inset of Fig. 2(c)]. The ripples decay with distance from the core and have an amplitude of up to . These ripples are associated with the roton mode mixing into the ground state [17, 18], and are physically associated with the energetic favourability of dipoles aligning headtotail, but where the van der Waals interactions and axial trapping prevent collapse.

(b) : The dipoles are now polarized offaxis, i.e. tipped in the direction of the axis. As a result, the circularsymmetry of the dipolar potential and density is lost, as evident from the asymmetry of the family of solutions in Fig. 3(b). For () the dipoles attract (repel) along and repel (attract) along . The local dipolar potential is still net positive throughout (since ). The vortex core is elongated along for , and reverses for . Close to the RI, density ripples form about the vortex (of large amplitude up to ) with wavefronts parallel to the axis. This alignment is due to the preferred headtotail alignment of dipoles which are tipped along , and related to the anisotropic mixing of the roton into the ground state [25].

(c) : For () the dipoles are net attractive (repulsive) along and repulsive (attractive) along . The anisotropy of the vortex core and ripples is reversed compared to .
5 Dipolar Potential of a Vortex
5.1 Dipoles Perpendicular to the Plane (=0)
A region of depleted density due to the presence of a vortex in an otherwise uniform system may be alternatively viewed as a lump of “antidipoles” (whose directions have been reversed) sitting in empty space [32, 21]. From this point of view a vortex treated as an object in its own right will generate its own mesoscopic dipolar potential. We now seek to estimate this dipolar potential analytically and compare our estimate against exact numerical calculations.
Consider first the simplified case of dipoles polarized perpendicular to the plane, . In Fig. 4(a) we show the dipolar potential in the vicinity of the vortex core for various values of . These values range from large negative (, black line) to large positive (, red line). In between we include a case close to the roton instability where the vortex possesses large ripples (, blue line) and close to the phonon instability, where the vortex core becomes narrow (, green line).
Far away from the core approaches , the homogeneous result. For , is dominated by the core structure and ripples (where present). Figure 4(b) plots, on a loglog scale, the decay of towards the homogeneous value. Far away from the core, decays as (grey line).
The longrange behaviour of the dipolar potential can be obtained analytically from the vortex ansatz [38, 21],
(9) 
where . This ansatz gives both the correct scaling of the density as and the correct decay to the homogeneous density as [29]. We evaluate the dipolar potential via the convolution result and by taking advantage of the cylindrical symmetry afforded by the scenario to employ Hankel transformations. The Hankel transform of Eq. (9) is,
where is a modified Bessel function of the second kind. We expand the dipolar interaction in terms of the condensate width parameter up to firstorder giving,
Then, up to first order in and third order in , the dipolar potential generated by the vortex ansatz is,
(10) 
with constants and . Recall is the dipolar potential at infinity.
To analyse this result it is useful to consider as the sum of a local term , proportional to density and given by,
(11) 
and nonlocal term . The nonlocal term can be viewed as arising from a fictitious electrostatic potential associated with the dipoles [37, 28] and has a complicated form given elsewhere [31]. It depends on variations of density, becoming zero for a homogeneous system.
The first bracketed term in Eq. (10) corresponds to the local contribution to the vortex dipolar potential. This is evident since, at longrange (), the vortex ansatz becomes , i.e. it decays as to the background density. Note that the swave interactions generate a local potential with identical form (proportional to density) [36]. The second bracketed term in Eq. (10) describes the nonlocal contribution, i.e. the contribution to the potential arising from the longrange contributions from all antidipoles in the vortex. This vanishes in the true 2D limit since the volume of antidipoles in the vortex core vanishes. This illustrates the point that the dipolar potential generated by a vortex is not a topological quantity like the potential associated with the hydrodynamic flow around the vortex (which leads to the logarithmic interaction between two vortices), but instead depends on the number of dipoles excluded from the region by the presence of the vortex. Also, the fact that the nonlocal vortex potential scales dominantly as , and not , informs us that the vortex does not strictly behave as a pointlike collection of dipoles at longrange. This is due to the slow powerlaw recovery of the vortex density to . This point is reinforced by noting that for sharplytruncated density profiles, e.g. a circular “hole”, and exponentiallydecaying density profiles, e.g. or , the dipolar potential does decay as at longrange. Due to the fast decay of the density, the inhomogeneity then really does resemble a giant dipole at longrange. This further confirms the subtle but critical role of the slow powerlaw decay of the vortex density. Nevertheless, whatever the subtleties of the nonlocal part of , its long range behaviour is dominated by the local scaling. This is confirmed by our numerical solutions, whose decays to the background value as at longrange, as shown in Fig. 4(b).
5.2 Dipoles with a Projection in the Plane ()
When the polarization direction is tilted away from the vertical, the dipolar interaction becomes anisotropic in the plane, resulting in the anisotropic vortex cores discussed in Section 4 above. The dipolar potential will then be anisotropic and have three contributions: an isotropic contribution from the homogeneous background density, a local contribution from the anisotropic density profile of the vortex, and a nonlocal anisotropic contribution arising from the antidipoles associated with the rest of the vortex.
We analyse the ensuing dipolar potential through the exemplar of and . The density profile of this vortex [Fig. 5(a)] shows a core elongated in with adjacent density ripples aligned in . The dipolar potential of the vortex [Fig. 5(b)] is indeed anisotropic, with a raised potential along and a reduced potential along (relative to the background dipolar potential ). Indeed, the angular dependence of , relative to , resembles that of the dipoledipole interaction itself . Thus, at least in its angular dependence, the vortex shares qualitatively the characteristics of a mesoscopic dipole.
Away from the core, decays towards the homogeneous value . At longrange , we again see the decay of developing a behaviour [Fig. 5(c)], consistent with the local contribution to from the vortex core. The same behaviour is achieved along the and directions; up to midscale distances , is dominated by the core structure, while at longer ranges, a common behaviour is obtained.
6 Conclusions
In this paper we have explored the stability and vortex solutions of the infinite twodimensional dipolar BoseEinstein condensate as a function of the strength and polarization angle of the dipoles. We mapped out the parameter regimes where a homogeneous density distribution in the plane is unstable towards either long wavelength phonon fluctuations or finite wavelength roton fluctuations. The instability is always towards configurations where the dipoles are arranged headtotail in order to benefit from the attractive portion of the dipolar interaction. The longrange and anisotropic nature of the dipolar interaction can lead to significant modifications to the structure of a single vortex, particularly in and around the core, including the formation of ripples in the core periphery and, for dipoles polarized offaxis, an anisotropic density profile. Furthermore, the inhomogeneous density of dipoles associated with the vortex generates a potential which decays slowly from the core and, for dipoles polarized offaxis, varies anisotropically in space. At longrange, this potential has and contributions, and an angular dependence . The longrange dipolar potential emanating from a single vortex in turn modifies the vortexvortex interaction [20] and is expected to support new regimes of vortex dynamics and collective structures in manyvortex systems.
DHJO acknowledges support from NSERC (Canada).
References
 [1] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler and T. Pfau, Phys. Rev. Lett. 94, 160401 (2005).
 [2] Q. Beaufils, R. Chicireanu, T. Zanon, B. LaburtheTolra, E. Maréchal, L. Vernac, J.C. Keller, and O. Gorceix, Phys. Rev. A 77, 061601(R) (2008).
 [3] M. Lu, N. Q. Burdick, S. H. Youn and B. L. Lev, Phys. Rev. Lett. 107, 190401 (2011).
 [4] K. Aikawa et al., Phys. Rev. Lett. 108, 210401 (2012).
 [5] T. Lahaye et al., Nature 448, 672 (2007).
 [6] T. Lahaye et al., Phys. Rev. Lett. 101, 080401 (2008).
 [7] D. H. J. O’Dell, S. Giovanazzi and G. Kurizki, Phys. Rev. Lett. 90, 110402 (2003).
 [8] L. Santos, G. V. Shlyapnikov and M. Lewenstein, Phys. Rev. Lett. 90, 250403 (2003).
 [9] S. Ronen, D. C. E. Bortolotti, and J. L. Bohn, Phys. Rev. Lett. 98, 030406 (2007).
 [10] R. Mottl, F. Brennecke, K. Baumann, R. Landig, T. Donner, and T. Esslinger, Science 336, 1570 (2012).
 [11] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein and T. Pfau, Rep. Prog. Phys. 72, 126401 (2009).
 [12] J. Werner, A. Griesmaier, S. Hensler, J. Stuhler, T. Pfau, A. Simoni and E. Tiesinga, Phys. Rev. Lett. 94, 183201 (2005).
 [13] T. Koch, T. Lahaye, J. Metz, B. Frölich, A. Griesmaier and T. Pfau, Nature Physics 4, 218 (2008).
 [14] K.K. Ni et al., Nature 464, 1324 (2010).
 [15] “Emergent Nonlinear Phenomena in BoseEinstein condensates”, P. G. Kevrekidis, D. J. Frantzeskakis and R. CarreteroGonzález (Eds.) (Springer, Berlin, 2008).
 [16] A. L. Fetter, J. Low Temp. Phys. 161, 445 (2010).
 [17] S. Yi and H. Pu, Phys. Rev. A 73, 061602(R) (2006).
 [18] R. M. Wilson, S. Ronen, J. L. Bohn, and H. Pu, Phys. Rev. Lett. 100, 245302 (2008); R. M. Wilson, S. Ronen, and J. L. Bohn, Phys. Rev. A 79, 013621 (2009).
 [19] M. Abad, M. Guilleumas, R. Mayol, M. Pi and D. M. Jezek, Phys. Rev. A 79, 063622 (2009).
 [20] B. C. Mulkerin, R. M. W. van Bijnen, A. M. Martin and N. G. Parker, Phys. Rev. Lett. 111, 170402 (2013).
 [21] D. H. J. O’Dell and C. Eberlein, Phys. Rev. A 75, 013604 (2007).
 [22] R. M. W. van Bijnen, D. H. J. O’Dell, N. G. Parker, and A. M. Martin, Phys. Rev. Lett. 98, 150401 (2007).
 [23] R. M. W. van Bijnen, A. J. Dow, D. H. J. O’Dell, N. G. Parker and A. M. Martin, Phys. Rev. A 80, 033617 (2009).
 [24] R. Kishor Kumar and P. Muruganandam, J. Phys. B 45, 215301 (2012).
 [25] C. Ticknor, R. M. Wilson and J. L. Bohn, Phys. Rev. Lett. 106, 065301 (2011).
 [26] N. R. Cooper, E. H. Rezayi and S. H. Simon, Phys. Rev. Lett. 95 200402 (2005); J. Zhang and H. Zhai, Phys. Rev. Lett. 95, 200403 (2005).
 [27] S. Yi and L. You, Phys. Rev. A 63, 053607 (2001).
 [28] N. G. Parker and D. H. J. O’Dell, Phys. Rev. A 78, 041601(R) (2008).
 [29] “BoseEinstein Condensation in Dilute Gases”, C. J. Pethick and H. Smith (Cambrige University Press, Cambridge, 2002)
 [30] U. R. Fischer, Phys. Rev. A 73, 031602 (2006)
 [31] Y. Cai, M. Rosenkranz, Z. Lei, W. Bao, Phys. Rev. A 82, 043623 (2010).
 [32] M. Klawunn, R. Nath, P. Pedri, and L. Santos, Phys. Rev. Lett. 100, 240403 (2008).
 [33] R. M. W. van Bijnen, N. G. Parker, S. J. J. M. F. Kokkelmans, A. M. Martin and D. H. J. O’Dell, Phys. Rev. A 82, 033612 (2010).
 [34] M. Klawunn and L. Santos, Phys. Rev. A 80, 013611 (2009).
 [35] “Quantized Vortices in Helium II”, R. J. Donnelly (Cambridge University Press, Cambridge, 1991)
 [36] H. Wu and D. W. L. Sprung, Phys. Rev. A 49, 4305 (1994).
 [37] D. H. J. O’Dell, S. Giovanazzi and C. Eberlein, Phys. Rev. Lett. 92, 250401 (2004).
 [38] A. L. Fetter, in Lectures in Theoretical Physics, K. T. Mahanthappa and W. E. Brittin (eds.) (Gordon and Breach, New York, 1969).