Two-Scale Vortex Construction in Multiscale Systems

• Two-Scale Vortex Construction is a strategy that organizes vortex phenomena by coupling localized vortex cores with larger-scale geometric or dynamical structures.
• It utilizes competing interaction ranges and nested geometries in systems such as layered superconductors, vortex-beam scattering, quantum droplets, and turbulence to achieve non-Abrikosov lattice configurations.
• This approach demonstrates practical multiscale control through simulation, experimental design, and effective theories, unifying diverse settings like vortex matter, quantum fluids, and continuum field models.

Searching arXiv for the cited papers and closely related work to ground the article. Searching "two-scale vortex construction superconducting layered vortices honeycomb square kagome" Searching "(Meng et al., 2014) Honeycomb square kagomé vortex lattices layered superconductors" “Two-Scale Vortex Construction” (Editor’s term) denotes a class of constructions in which vortex phenomena are organized simultaneously at two distinct scales, or are produced by coupling a local vortex structure to a larger geometric, dynamical, or collective scale. In the literature, the two scales may be two interaction lengths in a layered superconductor, two correlated transverse-momentum scales in a beam collision, an inner and outer annular scale in a nested quantum droplet, a forcing scale and a condensate scale in two-dimensional turbulence, or a filament core thickness and an order-one vortex geometry in Euler or Ginzburg–Landau theory [(Meng et al., 2014); (Ivanov, 2011); (Lin et al., 2020); (Parfenyev, 2022); (Davila et al., 2018); (Román, 2017)]. Taken together, these works treat two-scale vortex construction not as a single formalism, but as a recurring multiscale pattern in vortex matter, vortex beams, quantum fluids, and continuum field theory.

1. Conceptual range

A consistent theme across the literature is that the “two scales” are not always the same physical objects. In some settings they are two preferred interaction distances; in others they are a microscopic core and a macroscopic support curve; in others they are two correlated cone radii or a small forcing scale feeding a system-sized coherent vortex [(Meng et al., 2014); (Ivanov, 2011); (Parfenyev, 2022); (Averkiou et al., 15 Nov 2025)].

Setting	Two scales	Construction principle
Layered superconductors	Multiple $\lambda_i,\xi_j$ and lattice spacing	Competing repulsive and attractive intervortex ranges (Meng et al., 2014)
Vortex-beam scattering	$\varkappa_1,\varkappa_2$	Triangle-constrained triple-twisted amplitude (Ivanov, 2011)
Quantum droplets	Inner and outer annular radii	Nested vortex QD inside large-hole vortex QD (Lin et al., 2020)
2D turbulence and BECs	Forcing or pair scale and cluster scale	Evaporative heating, bounded-domain clustering, or forcing localization (Groszek et al., 2015, Gauthier et al., 2018, Parfenyev, 2022)
Vortex-lattice EFT	Coarse-grained lattice and long-wavelength mode	Coset construction and inverse Higgs constraints (Głódkowski et al., 17 Jul 2025)
Euler and GL desingularization	$O(\varepsilon)$ core and $O(1)$ filament, sheet, or ring geometry	Gluing, minimal connections, and thin-ring asymptotics (Davila et al., 2018, Román, 2017, Cao et al., 2021, Meyer et al., 2024)

This suggests that the term is best understood structurally rather than materially. The unifying feature is a controlled relation between a local vortex-bearing object and a larger-scale organizing framework.

2. Interaction-engineered vortex matter in superconductors

A direct construction of two-scale vortex matter was proposed for layered superconducting systems with multi-scale inter-vortex interaction (Meng et al., 2014). The physical starting point is the known possibility of nonmonotonic intervortex forces in multicomponent or layered superconductors: short-range repulsion, intermediate-range attraction, and, in layered geometries, an additional long-range repulsive tail. In a two-band system the paper quotes

$$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$

and for layered systems

$$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$

The layered form is the key engineering statement: each layer can contribute a repulsive magnetic term with its own $\lambda_i$, while each superconducting component or layer can contribute an attractive core-overlap term with its own $\xi_j$ (Meng et al., 2014).

The practical route is a layered heterostructure composed of different superconducting materials and/or different interlayer spacings. The paper identifies several control parameters: choice of superconducting materials, interlayer spacing, layer thickness, temperature, and distinct roles for upper and lower layers. The effective phenomenological potential used in simulations,

$$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$

encodes short-range repulsion, intermediate-range attraction, and long-range power-law repulsion, with the long-range tail associated with the demagnetization field and Pearl-type physics (Meng et al., 2014).

Within this framework, the construction target is not the conventional triangular Abrikosov lattice but non-Abrikosov symmetries. The paper reports, for one explicit potential, a honeycomb lattice at density $\rho=1.50$, a hexagonal lattice at $\varkappa_1,\varkappa_2$0, and a square lattice at $\varkappa_1,\varkappa_2$1; for a second potential it reports a kagomé lattice at $\varkappa_1,\varkappa_2$2 (Meng et al., 2014). The simulations use Langevin dynamics with approximately $\varkappa_1,\varkappa_2$3 to $\varkappa_1,\varkappa_2$4 vortices, random initial configurations, and simulated annealing to $\varkappa_1,\varkappa_2$5. The reported perfection values are close to unity, with coordination numbers within $\varkappa_1,\varkappa_2$6 in the honeycomb, hexagonal, square, and defect-free kagomé examples. In this construction, two-scale or multiscale behavior appears as interaction engineering: the preferred spacing is not single-scale, and lattice symmetry is selected by the detailed balance of competing ranges and by the vortex density.

3. Kinematical, nested, and compound constructions

A different two-scale mechanism appears in vortex-beam scattering. In a collision between one incoming Bessel vortex beam and one counterpropagating plane wave, the outgoing state is a pair of twisted beams entangled both in orbital angular momentum and in cone opening angles (Ivanov, 2011). A twisted one-particle state has fixed transverse-momentum magnitude $\varkappa_1,\varkappa_2$7 and definite $\varkappa_1,\varkappa_2$8-projection of OAM $\varkappa_1,\varkappa_2$9,

$O(\varepsilon)$0

The relevant scale is the conical momentum $O(\varepsilon)$1, or equivalently the opening angle

$O(\varepsilon)$2

For the triple-twisted amplitude, the two outgoing cone radii $O(\varepsilon)$3 are not independent: they must satisfy

$O(\varepsilon)$4

The final state is therefore not of the form $O(\varepsilon)$5; its OAM content and opening-angle content are jointly constrained by collision geometry (Ivanov, 2011). In this setting, two-scale vortex construction is kinematical rather than geometric: two outgoing transverse scales are created together, and postselection on one opening-angle sector sharpens the OAM correlations of the partner beam.

In quasi-2D Bose–Bose mixtures under thick transverse confinement, the same phrase describes a concentric composite object rather than a correlated scattering state (Lin et al., 2020). The effective equation is

$O(\varepsilon)$6

with cubic self-attraction and LHY-induced quartic repulsion. Stationary vortex droplets are sought as

$O(\varepsilon)$7

The paper reports stable vortex QDs for $O(\varepsilon)$8, and then constructs a nested state by embedding one vortex QD with a smaller topological number into another vortex QD with a larger number of topological charge (Lin et al., 2020). The explicit examples are an inner droplet $O(\varepsilon)$9 nested inside an outer droplet $O(1)$0, and an inner droplet $O(1)$1 nested inside the same outer droplet. Stability evidence is dynamical rather than spectral: both configurations persist beyond $O(1)$2, with “small turbulence” but without destruction of the composite state (Lin et al., 2020). Here the two scales are the inner droplet scale and the outer annular scale, and the construction relies on the large-hole geometry of high-$O(1)$3 droplets together with the stronger quartic LHY repulsion of the thick-confinement model.

4. Turbulent, negative-temperature, and forcing-mediated organization

In two-dimensional quantum turbulence and two-dimensional superfluids, two-scale vortex construction often appears as the coexistence of small-scale vortex structures and large-scale coherent organization. In highly oblate Bose–Einstein condensates in power-law traps, one paper identifies a small scale governed by vortex–antivortex pairs, vortexonium, annihilation, and sound emission, and a large scale given by same-sign Onsager clusters (Groszek et al., 2015). The main result is that Onsager vortex formation is strong in steep-walled near-uniform traps and suppressed in harmonic traps; the paper explicitly argues that this suppression is energetic, not a failure of evaporative heating. It also states that vortex–antivortex annihilation at zero temperature requires at least three vortices, and often involves three- and four-body processes (Groszek et al., 2015). In that setting, the two scales are linked by many-vortex dynamics rather than by a prescribed interaction potential.

An experimental realization of giant vortex clusters in a planar $O(1)$4 BEC made this bounded-domain picture concrete (Gauthier et al., 2018). The trap is a hard-walled ellipse, the central healing length is $O(1)$5, and a double-paddle stir directly injects a clustered high-energy state. The order parameter is the major-axis dipole moment

$O(1)$6

and the paper reports that the clustered state persists for up to $O(1)$7, approximately $O(1)$8 turnover times, remaining in the negative-temperature clustered regime. In this case the small scale is the individual quantized vortex core, while the large scale is the system-spanning dipolar arrangement selected by the elliptical boundary (Gauthier et al., 2018).

A related forcing-based construction appears in forced 2D turbulence on a $O(1)$9 periodic domain (Parfenyev, 2022). The universal condensate profile

$$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$0

is shown to be only the asymptotic small-forcing-scale limit. When the forcing correlation length $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$1 is finite, homogeneous forcing makes the condensate profile steeper than $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$2, while localized forcing can make it flatter than $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$3. The paper interprets this through a radius-dependent effective pumping intensity $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$4 and gives the diagnostic relation

$$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$5

Here the small scale is the forcing scale $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$6, and the large scale is the domain-sized coherent vortex condensate (Parfenyev, 2022).

Instability-driven turbulence with nonlinear damping furnishes yet another multiscale construction (Kan et al., 2023). The paper identifies four stationary branches: large-scale vortices, hybrid states with embedded shielded vortices, a dense shielded-vortex gas, and a hexagonal vortex crystal. The fundamental object is the shielded vortex, described as a tripolar vortex with an elliptical core and two opposite-sign satellites $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$7 apart, such that the circulation generated by any given vortex becomes small beyond a finite radius comparable to the largest forcing scale. The crystal emerges from the gas at $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$8, with nearest-neighbor spacing near $$E_\mathrm{int}=C_B^2K_0\left( \frac{r}{\lambda}\right) - C_1^2{2\pi} K_0 \left(\frac{r}{\xi_1}\right) - C_2^2K_0 \left(\frac{r}{\xi_2}\right),$$9, where $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$0 (Kan et al., 2023). This is a clear compound two-scale construction: each shielded vortex already has a core-plus-shield structure, and the collection of such vortices can then self-organize into a larger lattice.

5. Coarse-grained and long-wavelength effective descriptions

Some of the most explicit multiscale formulations are effective theories that begin from a coarse-grained medium and then reduce further to its gapless sector. For a rotating trapped condensate with a dense vortex lattice, one paper constructs an effective field theory by embedding the system in Newton–Cartan spacetime, identifying the magnetic Bargmann algebra, and introducing a compact phase field $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$1 together with two scalar fields $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$2 labeling vortices in the coarse-grained solid (Głódkowski et al., 17 Jul 2025). Around the homogeneous equilibrium

$$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$3

the coset construction yields covariant derivatives and inverse Higgs constraints

$$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$4

These eliminate the boost Goldstone, the rotational Goldstone, and the Kohn sector, leaving a single gapless mode. The leading low-energy action becomes

$$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$5

and the long-wavelength excitation is the Tkachenko mode with $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$6, while the gapped Kohn mode satisfies $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$7 (Głódkowski et al., 17 Jul 2025). The two scales here are the coarse-grained vortex crystal and the still longer-wavelength Tkachenko EFT.

An effective-string analogue appears for vortex rings (Gubser et al., 2014). The action

$$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$8

encodes local string tension, coupling to a background $$E_\mathrm{int} = \sum_i C_B{}_i^2K_0\left( \frac{r}{\lambda_i}\right) - \sum_jC_j^2{2\pi}K_0\left(\frac{r}{\xi_j}\right).$$9, and bilocal interaction. The short-distance divergence is absorbed into the dynamically generated scale

$\lambda_i$0

For a circular ring the translational velocity is

$\lambda_i$1

and the instability bands are centered near $\lambda_i$2. In head-on collisions, the late-time minimum separation is

$\lambda_i$3

while the most unstable wavelength satisfies $\lambda_i$4 as $\lambda_i$5 (Gubser et al., 2014). This is a particularly explicit two-scale separation: the core scale and its renormalized surrogate $\lambda_i$6 enter the effective tension, but the dominant instability lives at a much larger wavelength.

6. Desingularization and gluing in Euler and Ginzburg–Landau models

A large mathematical literature constructs two-scale vortex objects by desingularizing point vortices, filament measures, or thin rings into smooth or sheet-like structures. In two-dimensional Euler flow on a bounded domain, smooth $\lambda_i$7-vortex solutions can be built by gluing an outer Green-function field to inner Liouville bubbles (Davila et al., 2018). The core profile is

$\lambda_i$8

and the vorticity concentrates as

$\lambda_i$9

The two scales are explicit: the domain/intervortex scale is $\xi_j$0, while the vortex-core scale is $\xi_j$1 (Davila et al., 2018).

For the 3D Ginzburg–Landau functional, a quantitative vortex approximation constructs a polyhedral one-current $\xi_j$2 representing vortex filaments and proves the $\xi_j$3-level lower bound

$\xi_j$4

together with

$\xi_j$5

This is again a core-versus-filament construction: the energy pays core cost per unit filament length, while the approximating object is a macroscopic one-current (Román, 2017).

For vortex sheets, one paper starts from a macroscopic point-vortex relative equilibrium and replaces each point vortex by a microscopic closed vortex sheet of size $\xi_j$6, producing co-rotating $\xi_j$7-gons with angular velocity

$\xi_j$8

and translating dipoles with

$\xi_j$9

The support curves are analytic and strictly convex, and the construction is carried out directly in contour-dynamics variables through the Birkhoff–Rott operator (Cao et al., 2021). In this case the two scales are the $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$0 point-vortex skeleton and the $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$1 sheet geometry.

Several three-dimensional constructions make the scale separation geometrically explicit. Thin ring-shaped vortex sheets in the two-phase Euler equations are built with large ring radius $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$2 and small cross-sectional radius $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$3, with slenderness

$$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$4

and translation speed

$$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$5

in rescaled variables (Meyer et al., 2024). Smooth helical vortex filaments in $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$6 are constructed by gluing a compactly supported core profile solving $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$7 to a macroscopic helix, with inner variable $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$8 and limiting measure $$E_\mathrm{int}= c_1 e^{-r/\lambda} - c_2 e^{-r/\xi} + c_3 \frac{\lambda \{\tanh[\alpha(r-\beta)]+1\}}{r+\delta},$$9 as $\rho=1.50$0 (Averkiou et al., 15 Nov 2025). Time-periodic leapfrogging vortex rings are then realized for 3D Euler by coupling a slow Hamiltonian two-filament dynamics to fast internal transport of finite-thickness cores, using degenerate KAM-type analysis, pseudodifferential operator techniques, and Nash–Moser iteration; the solution is periodic in an almost uniformly translating frame and exists for parameters in a Cantor-like set of asymptotically full measure (García et al., 23 Mar 2026). Across these works, the recurring architecture is the same: an order-one filament, sheet, or ring geometry is resolved by an $\rho=1.50$1 core or cross-section.

7. Limits, misconceptions, and recurring principles

A common misconception is that two-scale vortex construction always means two independent scales. The beam-scattering example shows the opposite: $\rho=1.50$2 are constrained by a triangle inequality, and the OAM content is non-factorized (Ivanov, 2011). The superconducting example likewise requires simultaneous tuning of interaction profile and density; the same fabrication principle does not automatically produce every lattice at every density (Meng et al., 2014).

Another misconception is that all such constructions are exact analytic families. Some are numerical proof-of-principle studies. The superconducting lattice paper does not solve the full layered Ginzburg–Landau system directly and notes future work on vortex bending and non-pairwise forces (Meng et al., 2014). The nested quantum-droplet states are supported by long-time real-time evolution rather than by a separate BdG spectrum (Lin et al., 2020). The dense vortex-gas and vortex-crystal states in instability-driven turbulence are stationary branches identified from long-time DNS and continuation in $\rho=1.50$3, not Newton solves of exact steady states (Kan et al., 2023).

Conversely, not every reduced model includes explicit core physics. The vortex-lattice EFT from coset construction is not a microscopic derivation from individual vortex-core dynamics, but a symmetry-based coarse-grained theory of the long-wavelength sector (Głódkowski et al., 17 Jul 2025). The effective string treatment of vortex rings integrates out short-distance structure into $\rho=1.50$4, and the paper is explicit that reconnection requires microscopic physics outside the EFT (Gubser et al., 2014). In bounded BEC turbulence, suppression of Onsager clustering in harmonic traps is explained by the energy of the clustered branch, not by the absence of annihilation-driven heating (Groszek et al., 2015).

Taken together, these works suggest three recurrent principles. First, two-scale constructions are usually controlled by a small parameter, a bounded domain, or a symmetry reduction that separates local and global behavior. Second, the local object is rarely discarded: it reappears as a degree, a core profile, a shield, a cone radius, or a renormalized length such as $\rho=1.50$5. Third, large-scale order is typically selected either by interaction competition, by geometry, or by a reduced variational or Hamiltonian structure. In that sense, two-scale vortex construction is best understood as a family of multiscale strategies for connecting vortex cores, vortex interactions, and collective organization across very different physical theories.