---
title: Multi-Vortex States and Interactions
url: https://www.emergentmind.com/topics/multi-vortex
type: topic
---

# Multi-Vortex States and Interactions

Searching arXiv for recent and foundational papers relevant to the different technical meanings of “multi-vortex.”
“Multi-vortex” denotes a class of states, solutions, or configurations containing more than one vortex degree of freedom, but its precise meaning depends strongly on context. In nonlinear optics it refers to dissipative soliton complexes composed of several phase singularities sustained by structured gain landscapes [1109.0139]; in quantum mechanics it can mean multi-valued wavefunctions built from knotted vortex filaments [1911.00499]; in fluid dynamics it denotes interacting assemblies of vortices, including deformable viscous vortex particles and polygonal co-rotating configurations [1010.2475], [1506.05227]; in Bose–Einstein condensates it encompasses bifurcating polygonal vortex states, vortex crystals, multiply charged vortices that split into many singly charged lines, and bound vortex clusters in multicomponent gases [1701.01494], [1708.04336], [1505.00616], [1504.03203], [1505.06227], [1609.05489]; in superconductivity it includes controllable multi-flux configurations and vortex matter with multi-scale interactions [1905.06303], [1605.00524], [1404.4305]; and in photonics and electromagnetics it also describes multiple vector-vortex beams and multiplexed orbital-angular-momentum channels [1303.4487], [2502.17479]. The common structure is the coexistence and interaction of several vortical singularities, circulations, or topological windings within a single composite system, but the governing equations, stabilization mechanisms, and observable consequences vary substantially across subfields.

## 1. Optical and photonic multi-vortex states

In nonlinear optics, one of the most explicit uses of the term is the multi-vortex soliton complex supported by structured gain in a focusing Kerr medium with strong two-photon absorption [1109.0139]. The governing equation is a dissipative nonlinear Schrödinger equation,
\[
i\,\frac{\partial q}{\partial \xi} = -\frac{1}{2}\left(\frac{\partial^2 q}{\partial \eta^2} + \frac{\partial^2 q}{\partial \zeta^2}\right) - |q|^2\,q + i\,p\,R(\eta,\zeta)\,q - i\,\alpha\,|q|^2\,q,
\]
with annular gain channels
\[
R(\eta,\zeta) = \sum_{k=1}^{N} \exp\!\left[-\frac{\big(r_k(\eta,\zeta) - r_c\big)^2}{w^2}\right].
\]
Here the multi-vortex state is a dissipative attractor produced by a double balance: diffraction against Kerr self-focusing, and localized gain against two-photon absorption [1109.0139]. The topological charge is
\[
m = \frac{1}{2\pi} \oint \nabla \phi \cdot d\boldsymbol{\ell},
\]
and the work distinguishes “vortex twins” with identical charges and “vortex anti-twins” with opposite charges. Stable complexes with two, three, or four vortices exist above a minimal gain and minimal ring separation. In two-ring geometries, vortex twins are always asymmetric in intensity, whereas anti-twins can be symmetric near threshold and become asymmetric as gain increases or ring separation decreases [1109.0139].

The same broad optical theme reappears in integrated photonics, where “multi-vortex” refers not to interacting singularities in one nonlinear beam, but to the simultaneous generation of several vector-vortex beams by concentric ring resonators [1303.4487]. In that setting the emitted orbital angular momentum order is fixed by the angular phase-matching rule
\[
\ell = m - q,
\]
with \(m\) the ring azimuthal mode number and \(q\) the number of grating elements. Several concentric rings, each with its own vertically displaced bus waveguide, generate multiple OAM channels and arbitrary qudit superpositions \(|\psi\rangle = \sum_\ell c_\ell |\ell\rangle\) [1303.4487]. The same device can also operate as an OAM sorter with an average cross-talk of \(-10\,\mathrm{dB}\) between different OAM channels [1303.4487].

A related, but distinct, optical development is the Bessel-lattice necklace-vortex soliton. In azimuthally modulated Bessel potentials, the field evolves according to
\[
i\,\partial_{z}\psi + \tfrac{1}{2}\nabla_{\perp}^{2}\psi + \gamma\,|\psi|^{2}\psi + p\,V(r,\theta)\,\psi = 0,
\]
with
\[
V_{0}(r,\theta) = J_{1}^{2}\!\big(2\,b_{\mathrm{lin}}\,r\big)\,\cos(k\,\theta).
\]
Here “multi-vortex” means a compound state with \(M=k\) azimuthal beads per ring and total topological charge \(S\), so the adjacent-bead phase step is \(\Delta\phi = 2\pi S/M\) [2602.18703]. The lattice discretizes the azimuthal degree of freedom and can stabilize high-winding-number necklace vortices that would be azimuthally unstable in uniform Kerr media. Stable sextupoles, breathing octupoles, and 12-pole necklaces with \(S=3\) in a finite stability window are reported [2602.18703]. This suggests that, within nonlinear optics, “multi-vortex” spans both multi-core dissipative solitons and discretized many-bead vortex necklaces, unified by the role of structured media in suppressing azimuthal breakup.

## 2. Quantum and hydrodynamic formulations

In single-particle quantum mechanics, “multi-vortex” can mean multi-valued vortex solutions of the Schrödinger equation [1911.00499]. The wavefunction is written in Madelung form,
\[
\psi(x,t) = \sqrt{\rho(x,t)}\,e^{i S(x,t)/\hbar},
\]
with Bohmian velocity \(v = \nabla S/m\) and quantum potential
\[
Q(x,t) = - \frac{\hbar^2}{2m}\frac{\Delta\sqrt{\rho}}{\sqrt{\rho}}.
\]
The paper constructs \(S\) from vortex-filament potentials familiar from incompressible Euler flow, using a Seifert surface as the branch-cut structure, and shows that analytic continuation around a filament multiplies the wavefunction by \(e^{i m\Gamma N_w/\hbar}\) [1911.00499]. The nodal set is enforced by
\[
I_n^f(x) = \int_0^l \frac{d\sigma}{|x - R_f(\sigma)|^n}, \qquad
\psi(x,0) = \Phi(x)\,\frac{e^{i m \phi_f(x)/\hbar}}{I_n^f(x)},
\]
with \(n\ge 2\), so knotted or linked vortex filaments become nodal lines of \(\psi\) [1911.00499]. The analysis argues that multi-valued free-particle states radiate unless the single-valuedness condition \(m\Gamma/\hbar = 2\pi K\) holds, and it further proposes radiation signatures in Aharonov–Bohm configurations when the induced phase is not an integer multiple of \(2\pi\) [1911.00499].

In classical viscous fluid dynamics, multi-vortex states appear as interacting localized vorticity concentrations. The multi-moment vortex method represents the vorticity by Hermite–Gaussian moments,
\[
\omega(x,t) = \sum_{k_1,k_2\ge 0} M[k_1,k_2;t] \,\phi_{k_1,k_2}(x,t;\lambda),
\]
with the base Gaussian
\[
\phi_{00}(x,t;\lambda) = \frac{1}{\pi \lambda^2} \exp\!\left(-\frac{|x|^2}{\lambda^2}\right), \qquad \lambda^2 = \lambda_0^2 + 4\nu t.
\]
For a multi-vortex flow,
\[
\omega(x,t) = \sum_{j=1}^n \omega^j(x - x^j(t),t),
\]
and each particle carries a hierarchy of Hermite moments encoding deformation and anisotropy [1010.2475]. The central advance is a closed ODE system for self- and cross-interactions that avoids direct Biot–Savart evaluation for deformed particles. In merger problems this framework captures deformation, spiral-arm formation, and centroid drift far better than non-deforming Gaussian blobs [1010.2475].

A complementary hydrodynamic problem is the merger of \(N\ge 3\) co-rotating vortices initially arranged on a regular polygon [1506.05227]. In that case the multi-vortex interaction is qualitatively distinct from the two-vortex problem because an annular vortical structure forms as an intermediate state. The axisymmetric annulus is described by the generalized Lamb–Oseen solution
\[
\omega(r,t)
=\frac{\Gamma_{o}}{4\pi \nu t}\exp\!\left[-\frac{a^2+r^2}{4\nu t}\right]\,I_{0}\!\left(\frac{a\,r}{2\nu t}\right),
\]
which diffuses inward until it becomes a standard Lamb–Oseen vortex [1506.05227]. At moderate Reynolds number, this annular stage slows merger significantly. At higher Reynolds number, odd and even vortex numbers behave differently: even-\(N\) configurations evolve via pairing, whereas odd-\(N\) cases exhibit rapid and chaotic collapse [1506.05227]. A plausible implication is that multi-vortex interactions alter inverse-cascade timing in two-dimensional turbulence more strongly than pairwise-merger intuition would suggest.

## 3. Bose–Einstein condensates: bifurcations, crystals, and decay

In rotating trapped Bose–Einstein condensates, multi-vortex states arise by symmetry-breaking bifurcation from radially symmetric vortices. For the two-dimensional Gross–Pitaevskii equation
\[
i\,u_t = -(\partial_x^2+\partial_y^2)u + (x^2+y^2)u + |u|^2u + i\,\Omega\,(x\partial_y-y\partial_x)u,
\]
a charge-\(m_0\) vortex \(U(r,\theta)=e^{i m_0\theta}\psi_{m_0}(r)\) has a Hessian with \(m_0(m_0+1)/2\) pairs of negative eigenvalues near \(\Omega=0\), and as \(\Omega\) increases there are \(1+m_0(m_0-1)/2\) global bifurcations [1701.01494]. These branches include an asymmetric single vortex for \(m_0=1\), an asymmetric vortex pair for \(m_0=2\), and polygonal vortex arrays for \(m_0\ge 2\). The polygonal states are read off from zeros of Hermite–Gauss combinations and can include a central antivortex together with a ring of charge-one vortices [1701.01494].

In the many-vortex regime, a separate reduction from the rotating GP equation produces an effective ODE for co-rotating vortex crystals in isotropic or anisotropic traps [1708.04336]. In the isotropic case, the continuum limit yields an explicit areal density and a fold in the crystal branch, with maximum admissible vortex number
\[
N_{\max} = \frac{1}{\nu} \left\{ (\Omega \nu + 2)\left[\frac{1}{2}\ln(\Omega \nu + 2) - \ln 2 - \frac{1}{2}\right] + 2 \right\},
\]
where \(\nu = 1/\log(1/\epsilon)\) [1708.04336]. In strong anisotropy, vortices align on the long axis and the system reduces to an effective one-dimensional density problem. For two vortices, a pair on the long axis is linearly stable while a pair on the short axis is unstable, consistent with both the ODE reduction and full GP simulations [1708.04336].

A different BEC multi-vortex mechanism appears when a multiply charged vortex decays. In a cigar-shaped \({}^{87}\mathrm{Rb}\) condensate, a quadruply charged vortex imprinted by topological phase engineering decays into four singly charged lines that twist around one another as helical Kelvin waves [1505.00616]. The governing model is the three-dimensional GP equation
\[
i\hbar \frac{\partial \psi}{\partial t}
=
\left[
-\frac{\hbar^2}{2m}\nabla^2 + V(r) + g|\psi|^2
\right]\psi.
\]
The decay is not a simple radial splitting: axial inhomogeneity causes different \(z\)-slices to unwind at different times, so the emergent singly charged vortices acquire a helical structure [1505.00616]. The resulting multi-vortex state is then used numerically to generate an almost isotropic turbulent tangle, whose velocity statistics display power-law tails rather than Gaussian ones [1505.00616].

The stability of multiply charged vortices can also be analyzed variationally in quasi-two-dimensional condensates [1505.06227]. There the linearized collective-mode problem contains a vortex-coupled quadrupole mode \(Q_v\) whose frequency can become imaginary, providing a decay channel into singly charged vortices. Static Gaussian pinning with width comparable to the vortex-core radius can stabilize the state, and a sinusoidally modulated Gaussian potential produces a Mathieu-type system whose parametric-resonance tongues delimit additional stability windows [1505.06227]. This suggests that some multi-vortex configurations in BECs are not only interaction-selected but also dynamically engineered through mode control.

## 4. Multicomponent condensates and spin–orbit coupled vortex matter

In multicomponent homogeneous condensates, the decisive object is often the pairwise vortex interaction rather than the single-vortex spectrum [1504.03203]. For \(N_c\) components with coherent couplings,
\[
i\hbar \frac{\partial\psi_i}{\partial t}
=
\left[
-\frac{\hbar^2}{2m_i}\nabla^2 + V_i + \sum_j g_{ij}|\psi_j|^2
\right]\psi_i
-
\sum_{j\neq i} w_{ij}\psi_j,
\]
the energy contains density couplings \(g_{ij}\) and Rabi couplings \(w_{ij}=\hbar\Omega_{ij}\) [1504.03203]. In the miscible regime, intra-component vortices repel, but inter-component interactions can be attractive if \(\gamma_{12}<0\), and the competition between density repulsion and coherent phase locking can produce non-monotonic pair potentials with a minimum at finite separation [1504.03203]. This generates vortex dimers and, in three-component systems, equilateral trimers with bond length obeying the empirical scaling
\[
d_{\min}(\omega) \approx A(\gamma)\,\omega^p, \qquad p \approx -0.32
\]
for the parameter range studied [1504.03203]. The same pairwise mechanisms explain non-triangular many-vortex lattices and dimerized arrays under rotation, even though the rotating lattice problem itself is not solved directly in that work [1504.03203].

Spin–orbit coupling plus an optical lattice yields another route to exotic multi-vortex order [1609.05489]. In a two-component BEC with Rashba SOC in a square lattice, the single-particle minima move with SOC strength \(\lambda\) through four-fold, eight-fold, and Brillouin-zone-edge or corner degeneracies. When the minima touch the boundary of the first Brillouin zone, the interacting ground state can support simultaneously half-quantum vortex lattices, vortex–antivortex pair lattices, and fundamental vortex lattices [1609.05489]. The phase with this coexistence occurs when the corner minima \((\pm \pi,\pm \pi)\) are occupied in a parity-symmetric superposition and \(g_{\uparrow\downarrow}>g\), so that the two spin components populate complementary sublattices [1609.05489]. Here “multi-vortex” does not mean merely many vortices, but several topologically distinct vortex-lattice species embedded in the same ground state.

A related but distinct notion of layered vortex structure appears in generalized Maxwell–Higgs theories [1908.07871]. There the visible sector obeys BPS equations such as
\[
g' = \frac{a g}{r}, \qquad -\frac{a'}{r} = \frac{1-g^2}{f(\chi)},
\]
or, in a two-\(U(1)\) model,
\[
h' = \frac{c h}{r}, \qquad - \frac{c'}{q r} = q(w^2-h^2),
\]
with a dielectric function \(f\) controlled by an auxiliary neutral or hidden field [1908.07871]. The result is not a lattice of separate cores but a single radial vortex endowed with multiple magnetic rings. The number of rings is controlled by the dielectric modulation rather than by the winding \(n\), while the total flux remains quantized [1908.07871]. This suggests a useful distinction: in field theory, “multi-vortex” can refer either to several separate topological cores or to a single topological object with multilayered internal vorticity structure.

## 5. Superconductors, plasmas, and curved-surface vortex ensembles

In multiply connected superconducting structures, multi-vortex states can be encoded not as visible Abrikosov cores but as sets of loop winding numbers constrained by fluxoid quantization [1905.06303]. For a partitioned ring, the gauge-invariant phase differences \(\gamma_k\) satisfy
\[
\sum_{j:\,i\text{-th loop}} \gamma_j = 2\pi(n_i - \phi_i),
\]
and the free energy is a quadratic form in \(\tilde\phi_i = \phi_i - n_i\) [1905.06303]. The resulting stability regions form hexagonal honeycombs in flux space, each cell corresponding to a distinct multi-vortex configuration \(\hat n=(n_1,n_2,\dots)\). A micron-size magnetic particle on a cantilever locally tunes \((\phi_1,\phi_2,\dots)\), allowing stabilization of specific multi-vortex states and even braiding of individual vortices within a larger ensemble [1905.06303]. The braiding demonstrated is geometric rather than topologically protected, but it realizes deterministic exchange operations in a superconducting network [1905.06303].

A different superconducting use of “multi-vortex” concerns vortex matter with multi-scale inter-vortex interactions [1605.00524], [1404.4305]. A phenomenological pair potential of the form
\[
\frac{V(r)}{V_0}
=
e^{-r/\lambda}
-
c_2 e^{-r/\xi}
+
c_3\frac{\lambda\{\tanh[a(r-b)] + 1\}}{r+\delta}
\]
encodes short-range repulsion, intermediate-range attraction, and long-range Pearl-like repulsion [1605.00524]. Phase diagrams in the \((\xi/\lambda,\rho)\) plane then contain not only triangular lattices, but dimer, trimer, tetramer, stripe, honeycomb, kagome, and cluster phases [1605.00524]. Closely related layered-superconductor modeling shows that honeycomb, square, hexagonal, and kagomé vortex lattices can be stabilized by engineering multiple screening and core-overlap scales [1404.4305]. A plausible implication is that “multi-vortex” in superconductors increasingly refers to designer vortex matter whose unit objects may themselves be molecular clusters rather than isolated flux lines.

In dusty plasmas, the term describes several co-rotating vortices self-organized in an extended dust cloud [1701.03233]. The mechanism is a charge-gradient-driven vorticity source,
\[
S_\omega \propto \frac{1}{m_d}\nabla Q_d \times \mathbf{E}_{\text{total}} + \frac{1}{m_d}\nabla \times \mathbf{F}_i,
\]
generated by an inhomogeneous plasma background [1701.03233]. The steady rotation frequency follows the Vaulina scaling
\[
\omega \approx \left|\frac{F_i}{m_d \nu_{fr}}\right| \left|\frac{\beta}{e Z_0}\right|,
\]
and the characteristic vortex diameter is
\[
D_0 = \alpha \sqrt{\frac{\eta_k}{\omega^* + \nu_{fr}}}.
\]
As RF power decreases, the dust column shortens and the characteristic vortex size increases, driving a transition from three vortices to two and then one [1701.03233]. Here the multi-vortex state is a dissipative pattern of macroscopic circulation cells rather than a set of phase singularities.

On curved surfaces, point-vortex clusters on a torus introduce yet another meaning [2506.13158]. The toroidal geometry generates a geometric self-force through the Robin function, and the vortex dynamics is written in terms of the Schottky–Klein prime function and its logarithmic derivative
\[
K(\zeta)
=
\frac{1}{1-\zeta}
-
\frac{1}{2\pi A}\Psi_\rho\!\left(\frac{\log \zeta}{2\pi A}\right)
+
\frac{1}{2\pi A}\Psi_\rho\!\left(-\frac{\log \zeta}{2\pi A}\right),
\]
with \(\Psi_q\) the \(q\)-digamma function [2506.13158]. Chiral same-sign clusters remain geometrically confined and undergo collective toroidal drift, whereas achiral zero-net-circulation clusters disperse over the surface [2506.13158]. This suggests that on curved manifolds, confinement or dispersion of multi-vortex matter may be controlled as much by geometry and topology as by pairwise circulation signs.

## 6. Broader meanings, applications, and recurring mechanisms

Across photonics, condensed matter, fluid mechanics, plasma physics, and quantum foundations, several recurring mechanisms organize multi-vortex behavior. One is balance or competition among multiple physical scales. In dissipative optical media, stable multi-vortex solitons require balance between diffraction, Kerr self-focusing, gain, and two-photon absorption [1109.0139]. In superconductors and multicomponent condensates, unconventional multi-vortex matter emerges from competing inter-vortex forces acting on different length scales [1605.00524], [1504.03203], [1404.4305]. In layered Maxwell–Higgs models, auxiliary fields reshape a single vortex into a multilayer object through a tunable dielectric profile [1908.07871].

Another recurring mechanism is symmetry breaking from a simpler parent state. Rotating BEC multi-vortex polygons bifurcate from radially symmetric vortices as negative Hessian directions are removed by increasing rotation [1701.01494]. Multiply charged BEC vortices decay into multi-vortex sets when quadrupolar modes or axial inhomogeneity destabilize the original core [1505.00616], [1505.06227]. In SOC optical lattices, several vortex-lattice species coexist because corner degeneracies in the single-particle spectrum enforce particular coherent superpositions [1609.05489]. In superconducting networks, distinct multi-vortex states appear as neighboring cells of a discrete fluxoid honeycomb in control-parameter space [1905.06303].

A third theme is the tension between pairwise and collective descriptions. Some systems are well described by explicit pairwise interactions: multicomponent BEC vortex molecules [1504.03203], superconducting multi-scale vortex matter [1605.00524], and co-rotating vortex polygons in viscous fluids [1506.05227]. Other systems require intrinsically collective formulations: Hermite-moment closures for deformable viscous vortices [1010.2475], continuum crystal densities in rotating BECs [1708.04336], or toroidal Hamiltonians with global geometric regularization [2506.13158]. This suggests that “multi-vortex” is less a single object than a hierarchy of descriptions ranging from few-body topological molecules to many-body vortex matter.

Applications likewise reflect this breadth. Multi-ring gain landscapes permit deterministic generation of robust topological light structures [1109.0139]. Concentric integrated resonators generate and sort multiple OAM channels on chip [1303.4487], while reconfigurable metasurfaces convert several coaxial vortex-wave modes into spatially separated near-field focal spots for low-interference multi-user communication [2502.17479]. Superconducting multi-vortex control offers a route toward programmable braiding primitives [1905.06303]. Layered-superconductor vortex lattices have been proposed as magnetic templates for ultracold-atom quantum emulators [1404.4305], and dusty-plasma multi-vortices provide a macroscopic testbed for driven dissipative vorticity organization [1701.03233].

A common misconception is that “multi-vortex” always means many identical, separate vortex cores. The literature does not support that restriction. It may denote multi-core dissipative solitons [1109.0139], multi-valued vortex-filament wavefunctions [1911.00499], molecularly bound clusters [1504.03203], multilayer internal ring structures [1908.07871], fluxoid configurations in multiply connected superconductors [1905.06303], or several concurrent OAM channels in photonics [1303.4487], [2502.17479]. A plausible implication is that the most robust definition is not geometric multiplicity alone, but the presence of multiple coupled topological or vortical degrees of freedom within one physical state or dynamical system.

In that broad sense, multi-vortex research traces how topology, interaction range, geometry, dissipation, and external structuring conspire to create composite vortex matter. The specific equations differ—from dissipative nonlinear Schrödinger models to Gross–Pitaevskii, Ginzburg–Landau, point-vortex Hamiltonians, and constrained multicomponent GP systems—but the central problem remains the same: how multiple vortices coexist, interact, stabilize, bifurcate, or reorganize into higher-order structures.

Source: https://www.emergentmind.com/topics/multi-vortex