Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1708.09202 · Physical Review A

Motion of vortices in inhomogeneous Bose-Einstein condensates

Andrew J. Groszek, David M. Paganin, Kristian Helmerson, Tapio Simula

Published Feb 12, 2018• 55 citations
arXivHTMLPDF

Abstract

We derive a general and exact equation of motion for a quantized vortex in an inhomogeneous two-dimensional Bose-Einstein condensate. This equation expresses the velocity of a vortex as a sum of local ambient density and phase gradients in the vicinity of the vortex. We perform Gross-Pitaevskii simulations of single-vortex dynamics in both harmonic and hard-walled disk-shaped traps, and find excellent agreement in both cases with our analytical prediction. The simulations reveal that, in a harmonic trap, the main contribution to the vortex velocity is an induced ambient phase gradient, a finding that contradicts the commonly quoted result that the local density gradient is the only relevant effect in this scenario. We use our analytical vortex velocity formula to derive a point-vortex model that accounts for both density and phase contributions to the vortex velocity, suitable for use in inhomogeneous condensates. Although good agreement is obtained between Gross-Pitaevskii and point-vortex simulations for specific few-vortex configurations, the effects of nonuniform condensate density are in general highly nontrivial, and are thus difficult to efficiently and accurately model using a simplified point-vortex description.

DOI
10.1103/physreva.97.023617

Related papers

Velocity of vortices in inhomogeneous Bose–Einstein condensates

Halvor Møll Nilsen, Gordon Baym, C. J. Pethick

We derive, from the Gross-Pitaevskii equation, an exact expression for the velocity of any vortex in a Bose-Einstein condensate, in equilibrium or not, in terms of the condensate wave function at the center of the vortex. In general, the vortex velocity is a sum of the local superfluid velocity, plus a correction related to the density gradient near the vortex. A consequence is that in rapidly rotating, harmonically trapped Bose-Einstein condensates, unlike in the usual situation in slowly rotating condensates and in hydrodynamics, vortices do not move with the local fluid velocity. We indicate how Kelvin's conservation of circulation theorem is compatible with the velocity of the vortex center being different from the local fluid velocity. Finally, we derive an exact wave function for a single vortex near the rotation axis in a weakly interacting system, from which we derive the vortex precession rate.

arXiv:cond-mat/0508195•Proceedings of the National Academy of Sciences•May 15, 2006•32 citations

Dissipative motion of vortices in spatially inhomogeneous Bose-Einstein condensates

T. I. Zueva

A dissipative model of vortex motion within a two-dimensional Gross-Pitaevskii equation is studied. With the asymptotic coordination method of solutions, a system of ordinary first-order differential equations describing vortex motion in rotating Bose-Einstein condensate was derived. The model takes into account the effect of various external factors: magnetic trap frequencies, number of particles in the condensate, angular rotational velocity of the condensate, dissipation parameter, etc. In special cases of dissipation-free motion, the results are coordinated with the known results of other authors. The addition of dissipation generalizes the known equations and makes it possible to see the vortex motion to the equilibrium points and to define the equilibrium configurations of any number of vortices. The model is illustrated by a large number of examples.

Low Temperature Physics•Jan 1, 2019•3 citations

Dynamics of straight vortex filaments in a Bose–Einstein condensate with the Gaussian density profile

V. P. Ruban

The dynamics of interacting quantized vortex filaments in a rotating Bose–Einstein condensate existing in the Thomas–Fermi regime at zero temperature and obeying the Gross–Pitaevskii equation has been considered in the hydrodynamic “nonelastic” approximation. A noncanonical Hamilton equation of motion for the macroscopically averaged vorticity has been derived for a smoothly inhomogeneous array of filaments (vortex lattice) taking into account spatial nonuniformity of the equilibrium density of the condensate, which is determined by the trap potential. The minimum of the corresponding Hamiltonian describes the static configuration of the deformed vortex lattice against the preset density background. The condition of minimum can be reduced to a nonlinear second-order partial differential vector equation for which some exact and approximate solutions are obtained. It has been shown that if the condensate density has an anisotropic Gaussian profile, the equation of motion for the averaged vorticity has solutions in the form of a vector exhibiting a nontrivial time dependence, but homogeneous in space. An integral representation has also been obtained for the matrix Green function that determines the nonlocal Hamiltonian of a system of several quantized vortices of an arbitrary shape in a Bose–Einstein condensate with the Gaussian density. In particular, if all filaments are straight and oriented along one of the principal axes of the ellipsoid, we have a finitedimensional reduction that can describe the dynamics of the system of pointlike vortices against an inhomogeneous background. A simple approximate expression is proposed for the 2D Green function with an arbitrary density profile and is compared numerically with the exact result in the Gaussian case. The corresponding approximate equations of motion, describing the long-wavelength dynamics of interacting vortex filaments in condensates with a density depending only on transverse coordinates, have been derived.

arXiv:1612.00165•Journal of Experimental and Theoretical Physics•Jun 1, 2017•11 citations

Motion of quantum vortices on inhomogeneous backgrounds

Peter Mason, Natalia G. Berloff

The motion of quantized vortices on two-dimensional inhomogeneous density backgrounds with boundaries is considered numerically and asymptotically. We show that a Hamiltonian group relation together with the method of images and an approximation for the density background is useful to understand vortex motion. We analyze the vortex motion on a variety of background configurations motivated by experiments on trapped Bose-Einstein condensates including linear and quadratic traps. We show that close to the center of the condensate cloud in the Thomas-Fermi regime the vortex is moving predominantly due to the image vortex that is effectively shifted because of the large density depletion of the condensate at the boundary. Close to the boundary the vortex moves mostly due to a large background density gradient across the vortex core. The vortex velocity is a nonlinear combination of these effects and of the global shape of the background density. We find the complete families of the traveling coherent structures in an infinitely long two-dimensional channel and derive an approximation of the velocity of a single vortex moving close to the center of the channel in the Thomas-Fermi regime.

Physical Review A•Mar 20, 2008•28 citations

Properties of vortices in Bose–Einstein condensates

Amandine Aftalion

One of the key issues related to superfluidity is the existence of quantized vortices. Following recent experiments on Bose–Einstein condensates exhibiting vortices, we investigate the behavior of the wave function which solves the Gross–Pitaevskii equation. For a rotating Bose–Einstein condensate in a harmonic trap, we give a simplified expression of the Gross–Pitaevskii energy in the Thomas–Fermi regime, which only depends on the number and shape of the vortex lines. This allows us to study in detail the structure of the line of a single quantized vortex, which has a U or S shape. S type vortices exist for all values of the angular velocity Ω but are not minimizers of the energy while U vortices are minimizers and exist only for Ω sufficiently large. Finally, we study the drag created by the movement of a detuned laser beam in a condensate and the nucleation of vortices in the low density region close to the boundary layer of the cloud.

Comptes Rendus Physique•Jan 1, 2004•4 citations

Free expansion of Bose-Einstein condensates with quantized vortices

F. Dalfovo, Michele Modugno

The expansion of Bose-Einstein condensates with quantized vortices is studied by solving numerically the time-dependent Gross-Pitaevskii equation at zero temperature. For a condensate initially trapped in a spherical harmonic potential, we confirm previous results obtained by means of variational methods showing that, after releasing the trap, the vortex core expands faster than the radius of the atomic cloud. This could make the detection of vortices feasible, by observing the depletion of the density along the axis of rotation. We find that this effect is significantly enhanced in the case of anisotropic disk-shaped traps. The results obtained as a function of the anisotropy of the initial configuration are compared with the analytic solution for a noninteracting gas in three dimensions as well as with the scaling law predicted for an interacting gas in two dimensions.

arXiv:cond-mat/9907102•Physical Review A•Jan 14, 2000•62 citations
We derive a general and exact equation of motion for a quantized vortex in an inhomogeneous two-dimensional Bose-Einstein condensate. This equation expresses the velocity of a vortex as a sum of local ambient density and phase gradients in the vicinity of the vortex. We perform Gross-Pitaevskii simulations of single-vortex dynamics in both harmonic and hard-walled disk-shaped traps, and find excellent agreement in both cases with our analytical prediction. The simulations reveal that, in a harmonic trap, the main contribution to the vortex velocity is an induced ambient phase gradient, a finding that contradicts the commonly quoted result that the local density gradient is the only relevant effect in this scenario. We use our analytical vortex velocity formula to derive a point-vortex model that accounts for both density and phase contributions to the vortex velocity, suitable for use in inhomogeneous condensates. Although good agreement is obtained between Gross-Pitaevskii and point-vortex simulations for specific few-vortex configurations, the effects of nonuniform condensate density are in general highly nontrivial, and are thus difficult to efficiently and accurately model using a simplified point-vortex description.
We derive, from the Gross-Pitaevskii equation, an exact expression for the velocity of any vortex in a Bose-Einstein condensate, in equilibrium or not, in terms of the condensate wave function at the center of the vortex. In general, the vortex velocity is a sum of the local superfluid velocity, plus a correction related to the density gradient near the vortex. A consequence is that in rapidly rotating, harmonically trapped Bose-Einstein condensates, unlike in the usual situation in slowly rotating condensates and in hydrodynamics, vortices do not move with the local fluid velocity. We indicate how Kelvin's conservation of circulation theorem is compatible with the velocity of the vortex center being different from the local fluid velocity. Finally, we derive an exact wave function for a single vortex near the rotation axis in a weakly interacting system, from which we derive the vortex precession rate.
A dissipative model of vortex motion within a two-dimensional Gross-Pitaevskii equation is studied. With the asymptotic coordination method of solutions, a system of ordinary first-order differential equations describing vortex motion in rotating Bose-Einstein condensate was derived. The model takes into account the effect of various external factors: magnetic trap frequencies, number of particles in the condensate, angular rotational velocity of the condensate, dissipation parameter, etc. In special cases of dissipation-free motion, the results are coordinated with the known results of other authors. The addition of dissipation generalizes the known equations and makes it possible to see the vortex motion to the equilibrium points and to define the equilibrium configurations of any number of vortices. The model is illustrated by a large number of examples.
The dynamics of interacting quantized vortex filaments in a rotating Bose–Einstein condensate existing in the Thomas–Fermi regime at zero temperature and obeying the Gross–Pitaevskii equation has been considered in the hydrodynamic “nonelastic” approximation. A noncanonical Hamilton equation of motion for the macroscopically averaged vorticity has been derived for a smoothly inhomogeneous array of filaments (vortex lattice) taking into account spatial nonuniformity of the equilibrium density of the condensate, which is determined by the trap potential. The minimum of the corresponding Hamiltonian describes the static configuration of the deformed vortex lattice against the preset density background. The condition of minimum can be reduced to a nonlinear second-order partial differential vector equation for which some exact and approximate solutions are obtained. It has been shown that if the condensate density has an anisotropic Gaussian profile, the equation of motion for the averaged vorticity has solutions in the form of a vector exhibiting a nontrivial time dependence, but homogeneous in space. An integral representation has also been obtained for the matrix Green function that determines the nonlocal Hamiltonian of a system of several quantized vortices of an arbitrary shape in a Bose–Einstein condensate with the Gaussian density. In particular, if all filaments are straight and oriented along one of the principal axes of the ellipsoid, we have a finitedimensional reduction that can describe the dynamics of the system of pointlike vortices against an inhomogeneous background. A simple approximate expression is proposed for the 2D Green function with an arbitrary density profile and is compared numerically with the exact result in the Gaussian case. The corresponding approximate equations of motion, describing the long-wavelength dynamics of interacting vortex filaments in condensates with a density depending only on transverse coordinates, have been derived.
The motion of quantized vortices on two-dimensional inhomogeneous density backgrounds with boundaries is considered numerically and asymptotically. We show that a Hamiltonian group relation together with the method of images and an approximation for the density background is useful to understand vortex motion. We analyze the vortex motion on a variety of background configurations motivated by experiments on trapped Bose-Einstein condensates including linear and quadratic traps. We show that close to the center of the condensate cloud in the Thomas-Fermi regime the vortex is moving predominantly due to the image vortex that is effectively shifted because of the large density depletion of the condensate at the boundary. Close to the boundary the vortex moves mostly due to a large background density gradient across the vortex core. The vortex velocity is a nonlinear combination of these effects and of the global shape of the background density. We find the complete families of the traveling coherent structures in an infinitely long two-dimensional channel and derive an approximation of the velocity of a single vortex moving close to the center of the channel in the Thomas-Fermi regime.
One of the key issues related to superfluidity is the existence of quantized vortices. Following recent experiments on Bose–Einstein condensates exhibiting vortices, we investigate the behavior of the wave function which solves the Gross–Pitaevskii equation. For a rotating Bose–Einstein condensate in a harmonic trap, we give a simplified expression of the Gross–Pitaevskii energy in the Thomas–Fermi regime, which only depends on the number and shape of the vortex lines. This allows us to study in detail the structure of the line of a single quantized vortex, which has a U or S shape. S type vortices exist for all values of the angular velocity Ω but are not minimizers of the energy while U vortices are minimizers and exist only for Ω sufficiently large. Finally, we study the drag created by the movement of a detuned laser beam in a condensate and the nucleation of vortices in the low density region close to the boundary layer of the cloud.
The expansion of Bose-Einstein condensates with quantized vortices is studied by solving numerically the time-dependent Gross-Pitaevskii equation at zero temperature. For a condensate initially trapped in a spherical harmonic potential, we confirm previous results obtained by means of variational methods showing that, after releasing the trap, the vortex core expands faster than the radius of the atomic cloud. This could make the detection of vortices feasible, by observing the depletion of the density along the axis of rotation. We find that this effect is significantly enhanced in the case of anisotropic disk-shaped traps. The results obtained as a function of the anisotropy of the initial configuration are compared with the analytic solution for a noninteracting gas in three dimensions as well as with the scaling law predicted for an interacting gas in two dimensions.