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:1711.00280 · Journal of Experimental and Theoretical Physics

Stable and Unstable Vortex Knots in a Trapped Bose Condensate

V. P. Ruban

Published Mar 1, 2018• 5 citations
arXivHTMLPDF

Abstract

The dynamics of a quantum vortex toric knot TP,Q and other analogous knots in an atomic Bose condensate at zero temperature in the Thomas–Fermi regime is considered in the hydrodynamic approximation. The condensate has a spatially inhomogeneous equilibrium density profile ρ(z, r) due to the action of an external axisymmetric potential. It is assumed that z*= 0, r*= 1 is the point of maximum of function rρ(z, r), so that δ(rρ) ≈ –(α–)z2/2–(α + )(δr)2/2 for small z and δr. The geometrical configuration of a knot in the cylindrical coordinates is determined by a complex 2πP-periodic function A(ϕ, t) = Z(ϕ, t) + i[R(ϕ, t))–1]. When |A| ≪ 1, the system can be described by relatively simple approximate equations for P rescaled functions $${W_n}(\varphi ) \propto A(2\pi n + \varphi ):i{W_{n,t}} = - ({W_{n,\varphi \varphi }} + \alpha {W_n} - \in W_n^*)/2 - \sum\nolimits_{j \ne n} {1/(W_n^* - W_j^*)} $$ . For = 0, examples of stable solutions of type W n = θ n (ϕ–γt)exp(–iωt) with a nontrivial topology are found numerically for P = 3. In addition, the dynamics of various unsteady knots with P = 3 is modeled, and the tendency to the formation of a singularity over a finite time interval is observed in some cases. For P = 2 and small ≠ 0, configurations of type W0–W1 ≈ B0exp(iζ) + C(B0, α)exp(–iζ) + D(B0, α)exp(3iζ), where B0 > 0 is an arbitrary constant, ζ = k0ϕ–Ω0t + ζ0, k0 = Q/2, and Ω0 = (–α)/2–2/B02, which rotate about the z axis, are investigated. Wide stability regions for such solutions are detected in the space of parameters (α, B0). In unstable zones, a vortex knot may return to a weakly excited state.

DOI
10.1134/s1063776118030196

Related papers

Quasistable quantum vortex knots and links in anisotropic harmonically trapped Bose-Einstein condensates

Christopher Ticknor, V. P. Ruban, P. G. Kevrekidis

The long-time existence of topologically nontrivial configurations of quantum vortices in the form of torus knots and links in trapped Bose-Einstein condensates is demonstrated numerically within the three-dimensional Gross-Pitaevskii equation with an external anisotropic parabolic potential. We identify parametric domains of trap anisotropy, characterized by the axial over planar frequency ratio $\ensuremath{\lambda}\ensuremath{\approx}1.5--1.6$, where the lifetime of such quasistationary rotating vortex structures is many hundreds of typical rotation times. This suggests the potential experimental observability of the structures. We quantify the relevant lifetimes as a function of the model parameters (e.g., $\ensuremath{\lambda}$) and initial condition parameters of the knot profile.

arXiv:1903.02042•Physical Review A•Jun 10, 2019•13 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

Oscillations of rotating trapped Bose-Einstein condensates

Armen Sedrakian, Ira Wasserman

The tensor-virial method is applied for a study of oscillation modes of uniformly rotating Bose-Einstein condensed gases, whose rigid-body rotation is supported by an vortex array. The second-order virial equations are derived in the hydrodynamic regime for an arbitrary external harmonic trapping potential assuming that the condensate is a superfluid at zero temperature. The axisymmetric equilibrium shape of the condensate is determined as a function of the deformation of the trap; its domain of stability is bounded by the constraint \ensuremath{\Omega}1 on the rotation rate (measured in units of the trap frequency ${\ensuremath{\omega}}_{0}).$ The oscillations of the axisymmetric condensate are stable with respect to the transverse-shear and toroidal modes of oscillations, corresponding to the $l=2,|m|=1,2$ surface deformations. The eigenfrequencies of the modes are real and represent undamped oscillations. The condensate is also stable against quasiradial pulsation modes $(l=2,m=0),$ and its oscillations are undamped, if the superflow is assumed incompressible. In the compressible case we find that for a polytropic equation of state, the quasiradial oscillations are unstable when $\ensuremath{\gamma}(3\ensuremath{-}{\ensuremath{\Omega}}^{2})1\ensuremath{-}3{\ensuremath{\Omega}}^{2},$ and are stable otherwise. Thus, a dilute Bose gas, whose equation of state is polytropic with \ensuremath{\gamma}=2 to leading order in the diluteness parameter, is stable irrespective of the rotation rate. In nonaxisymmetric traps, the equilibrium constrains the (dimensionless) deformation in the plane orthogonal to the rotation to the domain ${A}_{2}{\ensuremath{\Omega}}^{2}$ with \ensuremath{\Omega}1. The second-harmonic-oscillation modes in nonaxisymmetric traps separate into two classes that have even or odd parity with respect to the direction of the rotation axis. Numerical solutions show that these modes are stable in the parameter domain where equilibrium figures exist.

arXiv:cond-mat/0007394•Physical Review A•May 10, 2001•17 citations

Stability of a Vortex in a Trapped Bose-Einstein Condensate

Anatoly A. Svidzinsky, Alexander L. Fetter

Based on the method of matched asymptotic expansions and on a time-dependent variational analysis, we study the dynamics of a vortex in a three-dimensional disk-shaped nonaxisymmetric condensate in the Thomas-Fermi limit. Both methods show that a vortex in a trapped Bose-Einstein condensate has formally unstable normal mode(s) with positive normalization and negative frequency, corresponding to a precession of the vortex line around the center of the trap. In a rotating trap, the solution becomes stable above an angular velocity Omega(m), characterizing the onset of metastability with respect to small transverse displacements of the vortex from the central axis.

arXiv:cond-mat/9811348•Physical Review Letters•Jun 26, 2000•154 citations

Three-Dimensional Numerical Simulation of Long-Lived Quantum Vortex Knots and Links in a Trapped Bose Condensate

V. P. Ruban

The dynamics of the simplest vortex knots, “unknots,” and torus links in an atomic Bose condensate at zero temperature in an anisotropic harmonic trap has been simulated numerically within the three-dimensional Gross–Pitaevskii equation. It has been found that such quasistationary rotating vortex structures exist for a very long time in wide ranges of the parameters of the system. This new result is qualitatively consistent with a previous prediction based on a simplified one-dimensional model approximately describing the motion of knotted vortex filaments.

arXiv:1809.07527•Journal of Experimental and Theoretical Physics Letters•Nov 1, 2018•15 citations

Parametric instability of oscillations of a vortex ring in a z-periodic Bose condensate and return to the initial state

V. P. Ruban

The dynamics of deformations of a quantum vortex ring in a Bose condensate with the periodic equilibrium density ρ(z) = 1 − ϵ cos z has been considered in the local induction approximation. Parametric instabilities of normal modes with the azimuthal numbers ±m at the energy integral E near the values $$E_m^{\left( p \right)} = 2m\sqrt {{m^2} - 1} /p$$ , where p is the order of resonance, have been revealed. Numerical experiments have shown that the amplitude of unstable modes with m = 2 and p = 1 can sharply increase already at ϵ ~ 0.03 to values about unity. Then, after several fast oscillations, fast return to a weakly perturbed state occurs. Such a behavior corresponds to the integrable Hamiltonian H ∝ σ(E 2 (1) − E)(|b +|2 + |b -|2)-ϵ(b + b - + b +*b -*)+u(|b +|4 + |b -|4)+w|b +|2|b -|2 for two complex envelopes b ±(t). The results have been compared to parametric instabilities of the vortex ring in the condensate with the density ρ(z, r) = 1 − r 2 − αz 2, which occur at α ≈ 8/5 and 16/7.

arXiv:1706.04348•Journal of Experimental and Theoretical Physics Letters•Aug 1, 2017•8 citations
The dynamics of a quantum vortex toric knot TP,Q and other analogous knots in an atomic Bose condensate at zero temperature in the Thomas–Fermi regime is considered in the hydrodynamic approximation. The condensate has a spatially inhomogeneous equilibrium density profile ρ(z, r) due to the action of an external axisymmetric potential. It is assumed that z*= 0, r*= 1 is the point of maximum of function rρ(z, r), so that δ(rρ) ≈ –(α–)z2/2–(α + )(δr)2/2 for small z and δr. The geometrical configuration of a knot in the cylindrical coordinates is determined by a complex 2πP-periodic function A(ϕ, t) = Z(ϕ, t) + i[R(ϕ, t))–1]. When |A| ≪ 1, the system can be described by relatively simple approximate equations for P rescaled functions Wn(φ)∝A(2πn+φ):iWn,t=−(Wn,φφ+αWn−∈Wn∗)/2−∑j≠n1/(Wn∗−Wj∗){W_n}(\varphi ) \propto A(2\pi n + \varphi ):i{W_{n,t}} = - ({W_{n,\varphi \varphi }} + \alpha {W_n} - \in W_n^*)/2 - \sum\nolimits_{j \ne n} {1/(W_n^* - W_j^*)} Wn​(φ)∝A(2πn+φ):iWn,t​=−(Wn,φφ​+αWn​−∈Wn∗​)/2−∑j=n​1/(Wn∗​−Wj∗​) . For = 0, examples of stable solutions of type W n = θ n (ϕ–γt)exp(–iωt) with a nontrivial topology are found numerically for P = 3. In addition, the dynamics of various unsteady knots with P = 3 is modeled, and the tendency to the formation of a singularity over a finite time interval is observed in some cases. For P = 2 and small ≠ 0, configurations of type W0–W1 ≈ B0exp(iζ) + C(B0, α)exp(–iζ) + D(B0, α)exp(3iζ), where B0 > 0 is an arbitrary constant, ζ = k0ϕ–Ω0t + ζ0, k0 = Q/2, and Ω0 = (–α)/2–2/B02, which rotate about the z axis, are investigated. Wide stability regions for such solutions are detected in the space of parameters (α, B0). In unstable zones, a vortex knot may return to a weakly excited state.
The long-time existence of topologically nontrivial configurations of quantum vortices in the form of torus knots and links in trapped Bose-Einstein condensates is demonstrated numerically within the three-dimensional Gross-Pitaevskii equation with an external anisotropic parabolic potential. We identify parametric domains of trap anisotropy, characterized by the axial over planar frequency ratio $\ensuremath{\lambda}\ensuremath{\approx}1.5--1.6$, where the lifetime of such quasistationary rotating vortex structures is many hundreds of typical rotation times. This suggests the potential experimental observability of the structures. We quantify the relevant lifetimes as a function of the model parameters (e.g., $\ensuremath{\lambda}$) and initial condition parameters of the knot profile.
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 tensor-virial method is applied for a study of oscillation modes of uniformly rotating Bose-Einstein condensed gases, whose rigid-body rotation is supported by an vortex array. The second-order virial equations are derived in the hydrodynamic regime for an arbitrary external harmonic trapping potential assuming that the condensate is a superfluid at zero temperature. The axisymmetric equilibrium shape of the condensate is determined as a function of the deformation of the trap; its domain of stability is bounded by the constraint \ensuremath{\Omega}1 on the rotation rate (measured in units of the trap frequency ${\ensuremath{\omega}}_{0}).$ The oscillations of the axisymmetric condensate are stable with respect to the transverse-shear and toroidal modes of oscillations, corresponding to the l=2,∣m∣=1,2l=2,|m|=1,2l=2,∣m∣=1,2 surface deformations. The eigenfrequencies of the modes are real and represent undamped oscillations. The condensate is also stable against quasiradial pulsation modes (l=2,m=0),(l=2,m=0),(l=2,m=0), and its oscillations are undamped, if the superflow is assumed incompressible. In the compressible case we find that for a polytropic equation of state, the quasiradial oscillations are unstable when $\ensuremath{\gamma}(3\ensuremath{-}{\ensuremath{\Omega}}^{2})1\ensuremath{-}3{\ensuremath{\Omega}}^{2},$ and are stable otherwise. Thus, a dilute Bose gas, whose equation of state is polytropic with \ensuremath{\gamma}=2 to leading order in the diluteness parameter, is stable irrespective of the rotation rate. In nonaxisymmetric traps, the equilibrium constrains the (dimensionless) deformation in the plane orthogonal to the rotation to the domain ${A}_{2}{\ensuremath{\Omega}}^{2}$ with \ensuremath{\Omega}1. The second-harmonic-oscillation modes in nonaxisymmetric traps separate into two classes that have even or odd parity with respect to the direction of the rotation axis. Numerical solutions show that these modes are stable in the parameter domain where equilibrium figures exist.
Based on the method of matched asymptotic expansions and on a time-dependent variational analysis, we study the dynamics of a vortex in a three-dimensional disk-shaped nonaxisymmetric condensate in the Thomas-Fermi limit. Both methods show that a vortex in a trapped Bose-Einstein condensate has formally unstable normal mode(s) with positive normalization and negative frequency, corresponding to a precession of the vortex line around the center of the trap. In a rotating trap, the solution becomes stable above an angular velocity Omega(m), characterizing the onset of metastability with respect to small transverse displacements of the vortex from the central axis.
The dynamics of the simplest vortex knots, “unknots,” and torus links in an atomic Bose condensate at zero temperature in an anisotropic harmonic trap has been simulated numerically within the three-dimensional Gross–Pitaevskii equation. It has been found that such quasistationary rotating vortex structures exist for a very long time in wide ranges of the parameters of the system. This new result is qualitatively consistent with a previous prediction based on a simplified one-dimensional model approximately describing the motion of knotted vortex filaments.
The dynamics of deformations of a quantum vortex ring in a Bose condensate with the periodic equilibrium density ρ(z) = 1 − ϵ cos z has been considered in the local induction approximation. Parametric instabilities of normal modes with the azimuthal numbers ±m at the energy integral E near the values Em(p)=2mm2−1/pE_m^{\left( p \right)} = 2m\sqrt {{m^2} - 1} /pEm(p)​=2mm2−1​/p , where p is the order of resonance, have been revealed. Numerical experiments have shown that the amplitude of unstable modes with m = 2 and p = 1 can sharply increase already at ϵ ~ 0.03 to values about unity. Then, after several fast oscillations, fast return to a weakly perturbed state occurs. Such a behavior corresponds to the integrable Hamiltonian H ∝ σ(E 2 (1) − E)(|b +|2 + |b -|2)-ϵ(b + b - + b +*b -*)+u(|b +|4 + |b -|4)+w|b +|2|b -|2 for two complex envelopes b ±(t). The results have been compared to parametric instabilities of the vortex ring in the condensate with the density ρ(z, r) = 1 − r 2 − αz 2, which occur at α ≈ 8/5 and 16/7.