--- title: Gross–Pitaevskii Vortex Dynamics url: https://www.emergentmind.com/topics/gross-pitaevskii-vortex-dynamics type: topic --- # Gross–Pitaevskii Vortex Dynamics Gross–Pitaevskii vortex dynamics concerns the motion, interaction, deformation, and reconnection of quantized vortices in solutions of the Gross–Pitaevskii (GP) equation. In the GP description of a dilute, zero-temperature Bose–Einstein condensate, the condensate is represented by a single complex order parameter $\psi$, and vortices appear as isolated zeros of $|\psi|$ around which the phase winds by $2\pi n$. In two spatial dimensions these defects are point-like in the Jacobian or vorticity measure; in three dimensions they form vortex filaments that may close into rings or more complicated topological structures such as torus knots. Across asymptotic analysis, reduced-order models, direct simulations, and tracking algorithms, the subject links Hamiltonian point-vortex motion, filament dynamics, reconnection physics, Kelvin-wave cascades, and the effects of trapping, inhomogeneity, dissipation, and multicomponent structure [1307.7250]. ## 1. Governing equations, conserved quantities, and vortex observables In dimensional form, the zero-temperature GP equation for a condensate wavefunction $\psi$ is $$ i\hbar\,\partial_t\psi = -\frac{\hbar^2}{2m}\nabla^2\psi + \frac{4\pi\hbar^2 a_s}{m}|\psi|^2\psi + V_{\rm ext}\psi, $$ with boson mass $m$, $s$-wave scattering length $a_s$, and external potential $V_{\rm ext}$. In the homogeneous setting $V_{\rm ext}\equiv 0$, one convenient nondimensionalization uses the healing length $\xi=\sqrt{1/(8\pi a_s\rho_\infty)}$, the time scale $m/(8\pi\hbar a_s\rho_\infty)$, and density scale $\rho_\infty$, yielding $$ i\,\partial_t\psi + \tfrac12\nabla^2\psi - \tfrac12|\psi|^2\psi = 0. $$ The corresponding GP energy functional is $$ E[\psi]=\int\Bigl(\tfrac12|\nabla\psi|^2+\tfrac14|\psi|^4\Bigr)\,dV, $$ and excitation energies are measured relative to the uniform ground state [1307.7250]. Trapped condensates introduce a confining potential and a different natural scaling. For an axisymmetric harmonic trap, $$ V_{\rm trap}(\mathbf r)=\tfrac12 m(\omega_r^2 r_\perp^2+\omega_z^2 z^2), $$ one standard dimensionless form is $$ i\frac{\partial\psi}{\partial t} = \Bigl[ -\frac12\nabla^2 +\frac12(\lambda^2 r_\perp^2+z^2) +g|\psi|^2 \Bigr]\psi, \qquad \lambda=\omega_r/\omega_z, $$ with $\int |\psi|^2\,d^3r=1$ preserved [1309.7012]. Several asymptotic and rigorous studies instead use the small-core form $$ i\,\partial_t\psi = -\Delta\psi +\frac{1}{\varepsilon^2}\bigl(|\psi|^2-1\bigr)\psi \quad\text{or}\quad i\,\varepsilon^2\partial_t u+\Delta u+(1-|u|^2)u=0, $$ where $\varepsilon$ is the vortex-core scale. In this regime the Ginzburg–Landau energy and the supercurrent $$ j(\psi)=\frac12\bigl(\overline\psi\,\nabla\psi-\psi\,\nabla\overline\psi\bigr) $$ play a central role, while the Jacobian $$ J\psi=\tfrac12\nabla\times j(\psi)\approx \sum_j \pi d_j\,\delta_{a_j} $$ encodes vortex positions $a_j$ and degrees $d_j=\pm1$ [2404.02133]. The hydrodynamic representation $\psi=\sqrt{\rho}\,e^{i\phi}$ identifies $\rho=|\psi|^2$ as the condensate density and $\mathbf v=\nabla\phi$ as the superfluid velocity away from cores. For computational geometry, the pseudo-vorticity field $$ \Omega=\nabla\psi_R\times\nabla\psi_I $$ is especially useful: it is tangent to the vortex filament and remains smooth even though the physical vorticity is singular on the nodal set. This quantity underlies spectrally accurate filament-tracking algorithms in periodic domains [1604.03595]. ## 2. Reduced vortex laws: point vortices, inhomogeneity, and filaments In the small-core limit on the plane, the leading-order dynamics of well-separated vortices are governed by the classical Helmholtz–Kirchhoff or Kirchhoff–Onsager point-vortex system. For degree-$\pm1$ vortices at positions $a_j^0(t)$, one asymptotic form is $$ \dot a_j^0 = 2\sum_{k\neq j} d_k\, \frac{(a_j^0-a_k^0)^\perp}{|a_j^0-a_k^0|^2}, $$ with Hamiltonian $$ K(a)=2\sum_{j25.17$. Increasing the optical-lattice strength expands the survival window to larger norms [1210.4040]. Annular geometries make nucleation and boundary injection particularly explicit. In a two-dimensional ring potential with Gaussian radial profile and a rotating sinusoidal perturbation, $$ V_{\rm rot}(r,\theta,t)=\varepsilon\,V_{\rm pot}(r)\sin(n_\theta\theta-\Omega t), $$ a finite-volume Strang-splitting computation with benchmark parameters $m=10$, $V_0=100$, $\gamma=100$, $r_{\min}=0.6$, $r_{\max}=1.4$, $\varepsilon=0.05$, and $n_\theta=6$ observed no vortices for $\Omega=\pi$ over $t\in[0,3]$, whereas for $\Omega=10\pi/3$ vortices nucleated from both the inner and outer edges and, by $t=3$, produced 12 vortices: six of index $+1$ near the outer edge and six of index $-1$ near the inner edge [2404.10360]. These results delimit several distinct notions of stability in GP vortex dynamics. In rotating traps, stability refers to persistent co-rotating equilibria or crystals; in quenched attractive systems it refers to survival against collapse and against topological demotion; in annular stirring problems it refers to nucleation thresholds, edge-mediated injection, and the long-time arrangement of defects. A plausible implication is that GP vortices are best viewed not as a single universal defect class but as topological excitations whose effective dynamics depend sensitively on the competition among core energy, background density, trapping geometry, and forcing protocol. ## 5. Tracking vortex filaments, turbulent tangles, and reduced-order numerics Direct numerical study of GP vortex dynamics depends on the ability to identify vortex cores accurately in the complex field. In periodic domains, a robust method is to locate points where $\Re\psi=0$ and $\Im\psi=0$ by Newton–Raphson iteration in the plane normal to the pseudo-vorticity $\Omega=\nabla\psi_R\times\nabla\psi_I$. One first finds low-density seed points, constructs the normal plane using the local tangent $\hat t=\Omega/|\Omega|$, iterates to machine precision, and then marches along the filament with step size $\zeta\lesssim O(\xi)$ until closure or periodic re-entry. From the reconstructed curves one computes curvature $$ \kappa(s)=\Bigl|\frac{d^2r}{ds^2}\Bigr|, $$ torsion $$ \tau(s)= \frac{\bigl(\frac{dr}{ds}\times\frac{d^2r}{ds^2}\bigr)\cdot\frac{d^3r}{ds^3}} {\bigl|\frac{d^2r}{ds^2}\bigr|^2}, $$ and Kelvin-wave occupations from Fourier transforms of transverse displacements [1604.03595]. Applied to a GP-driven turbulent tangle generated from a Taylor–Green configuration, this tracking framework resolves the total vortex-line density $$ \mathcal L(t)=\frac{1}{V}\sum_i\oint_{\mathcal C_i}ds $$ and shows Vinen decay, $$ \frac{d\mathcal L}{dt} = -\chi_2\frac{\Gamma}{2\pi}\mathcal L^2, $$ with best-fit coefficient $\chi_2\approx0.65$, in quantitative agreement with low-temperature He II experiments. The same data reveal rare linked rings during decay, curvature PDFs with exponential tails at large $\kappa$, torsion PDFs with symmetric $|\tau|^{-3}$ tails, and a Kelvin-wave spectrum $$ n_k\propto k^{-11/3} $$ in the inertial range, consistent with the L’vov–Nazarenko weak-wave-turbulence prediction [1605.00567]. For two-dimensional small-core problems, vortex tracking also enables reduced-order simulation. A numerical strategy based on the reduced Hamiltonian point-vortex system integrates the ODEs for vortex centers, reconstructs a smoothed GP field from the canonical harmonic map and precomputed radial core profile, and proves rigorous supercurrent-error estimates of the form $$ \|J\psi_\varepsilon(t)-\pi\sum_j d_j\delta_{a_j(t)}\|_{\dot W^{-1,1}} \lesssim \varepsilon^\beta, \qquad \|j(\psi_\varepsilon(t))-j(u^*(a(t),d))\|_{L^{4/3}} \lesssim \varepsilon^\gamma, $$ up to the first vortex-collision time. In practice, the reconstructed field reproduces the density dips and far-field phase at a tiny fraction of the cost of a full GP solve when $\varepsilon$ is small [2404.02133]. Boundary conditions can themselves be a dynamical tool. Quasi-periodic boundary conditions generalize ordinary periodicity by keeping $|\psi|^2$ strictly periodic while allowing prescribed phase windings compatible with net vorticity. In this framework the phase jump over a lattice translation $\mathbf L$ is fixed by the total vortex count $\mathbf N_v$, and Kelvin’s circulation theorem together with the boundary matching implies that the net vortex number is conserved. Real-time simulations use a third-order Adams–Bashforth integrator and nearest-neighbour vortex matching, and they reveal bulk phenomena unavailable in periodic boxes with zero net circulation, including perfectly periodic vortex arrays, vortex depinning, pair nucleation, and Kármán vortex streets. The same formulation is proposed as a toy model for neutron-star bulk vortex dynamics [2509.15298]. ## 6. Multicomponent, stochastic, and nonequilibrium generalizations In multicomponent condensates beyond the Thomas–Fermi regime, variational reductions show that vortex cores may acquire inertia. For a two-dimensional $N$-component GP system with one vortex per component, elimination of phase-gradient collective coordinates yields an effective Hamiltonian $$ H = \sum_{j=1}^N \Bigl[ \frac12(\mathbf p_j-n_j\mathbf A_j)^2 +\frac{k}{2}|\mathbf R_j|^2 \Bigr] + \sum_{i