Stable and unstable vortex knots in a trapped Bose-Einstein condensate

TL;DR

This paper investigates the dynamics and stability of vortex knots in a trapped Bose-Einstein condensate, revealing stable solutions, potential singularities, and parameter-dependent stability regions through analytical and numerical methods.

Contribution

It introduces a simplified model for vortex knot dynamics in BECs, finds stable solutions for certain parameters, and explores the conditions leading to instability and singularity formation.

Findings

Stable vortex knot solutions for P=3 found at epsilon=0.

Dynamics indicating possible finite-time singularities.

Wide stability regions identified for P=2 configurations.

Abstract

The dynamics of a quantum vortex torus knot ${\cal T}_{P,Q}$ and similar knots in an atomic Bose-Einstein condensate at zero temperature in the Thomas-Fermi regime has been considered in the hydrodynamic approximation. The condensate has a spatially nonuniform equilibrium density profile $\rho(z,r)$ due to an external axisymmetric potential. It is assumed that $z_*=0$, $r_*=1$ is a maximum point for function $r\rho(z,r)$, with $\delta (r\rho)\approx-(\alpha-\epsilon) z^2/2 -(\alpha+\epsilon) (\delta r)^2/2$ at small $z$ and $\delta r$. Configuration of knot in the cylindrical coordinates is specified by a complex $2\pi P$-periodic function $A(\varphi,t)=Z(\varphi,t)+i [R(\varphi,t)-1]$. In the case $|A|\ll 1$ the system is described by relatively simple approximate equations for re-scaled functions $W_n(\varphi)\propto A(2\pi n+\varphi)$, where $n=0,\dots,P-1$, and $iW_{n,t}=-(W_{n,\varphi\varphi}+\alpha W_n -\epsilon W_n^*)/2-\sum_{j\neq n}1/(W_n^*-W_j^*)$. At $\epsilon=0$, numerical examples of stable solutions as $W_n=\theta_n(\varphi-\gamma t)\exp(-i\omega t)$ with non-trivial topology have been found for $P=3$. Besides that, dynamics of various non-stationary knots with $P=3$ was simulated, and in some cases a tendency towards a finite-time singularity has been detected. For $P=2$ at small $\epsilon\neq 0$, rotating around $z$ axis configurations of the form $(W_0-W_1)\approx B_0\exp(i\zeta)+\epsilon C(B_0,\alpha)\exp(-i\zeta) + \epsilon D(B_0,\alpha)\exp(3i\zeta)$ have been investigated, where $B_0>0$ is an arbitrary constant, $\zeta=k_0\varphi -\Omega_0 t+\zeta_0$, $k_0=Q/2$, $\Omega_0=(k_0^2-\alpha)/2-2/B_0^2$. In the parameter space $(\alpha, B_0)$, wide stability regions for such solutions have been found. In unstable bands, a recurrence of the vortex knot to a weakly excited state has been noted to be possible.