Borromean Hydrodynamics

Updated 7 July 2026

Borromean hydrodynamics is the theory of multicomponent supercounterfluids (M ≥ 3) where only relative-phase transport is active while total mass flow remains zero.
It hinges on compact-gauge invariance that forces the hydrodynamics to depend solely on phase differences, yielding M – 1 counterflow modes and M distinct elementary vortices.
This framework predicts a unique universality class, revealing modified BKT criticality and novel vortex interactions that distinguish it from ordinary two-component superfluids.

Borromean hydrodynamics is primarily the long-wavelength, finite-temperature theory of multicomponent supercounterfluids with $M\ge 3$ in which dissipationless transport survives only in relative-flow channels while the total mass current is arrested. Its defining structure is a compact-gauge invariance under equal shifts of all component phases by a compact field, so the hydrodynamics depends only on phase differences, supports $M-1$ counterflow Nambu–Goldstone modes, and admits $M$ elementary vortex species defined modulo the uniform vector $(1,1,\ldots,1)$ , a mismatch that distinguishes the Borromean regime from ordinary and two-component superfluids (Golic et al., 2024, Babaev et al., 2023). A distinct later usage applies the same adjective to a covariant intersection-theoretic formulation of fluid dynamics in which thermodynamics, spacetime geometry, and differential topology enter as an irreducible triple structure (Nekrasov et al., 31 Dec 2025).

1. Definition and distinguishing content

A Borromean supercounterfluid is a multicomponent superfluid state with $M\ge 3$ components in which supertransport is possible only in counterflow. Net mass flow is arrested, while dissipationless transport occurs in relative channels. In this state, composite vortices with equal winding in all components proliferate and cost no macroscopic energy; only counterflow, or relative-phase, excitations remain hydrodynamically active (Golic et al., 2024). This is the physical origin of the term “Borromean” in the condensed-matter literature: no single component, and no uniform co-flow of all components, furnishes an autonomous hydrodynamic transport channel.

The contrast with lower-component cases is sharp. In a single-component superfluid, supertransport occurs in the total phase, vortices carry a single $U(1)$ charge, and the two-dimensional transition is the standard XY/BKT transition with the usual Nelson–Kosterlitz jump. In a two-component counterflow system with net flow arrested, the theory reduces to an effective single-component superfluid of a relative phase, so the vortex plasma remains scalar-charged and the universal properties are those of a standard XY/BKT theory. For $N\ge 3$ , by contrast, the vortex plasma is multicomponent or vector-charged, and the supercounterfluid-to-normal transition belongs to a distinct universality class (Golic et al., 2024).

A second hallmark is the mismatch between the number of elementary defects and the number of hydrodynamic sound branches. For $N\ge 3$ , the long-wavelength theory has $N-1$ independent phonon modes but $N$ distinct elementary vortex excitations. For $M-1$ 0, the compact-gauge identification collapses the two unit vortices into a single elementary type, so the counting of defects and modes coincides. The $M-1$ 1 versus $M-1$ 2 mismatch is therefore a defining structural feature of the Borromean regime rather than a generic property of counterflow hydrodynamics (Babaev et al., 2023).

2. Compact-gauge invariance and hydrodynamic variables

The organizing principle of Borromean hydrodynamics is compact-gauge invariance,

$M-1$ 3

where $M-1$ 4 is a compact phase field defined modulo $M-1$ 5 and is allowed to carry vortices. Because the phase fields $M-1$ 6 are themselves compact, the gauge transformation includes topological content rather than merely smooth reparameterizations. Gauge-invariant long-wave expressions must therefore depend only on differences $M-1$ 7, or, equivalently, on a stiffness matrix $M-1$ 8 satisfying

$M-1$ 9

so that the uniform vector $M$ 0 is annihilated and the action is restricted to the subspace orthogonal to the common phase (Golic et al., 2024).

In matrix form the effective action is

$M$ 1

In the component-symmetric case, the action depends only on pairwise phase-gradient differences, and a convenient gauge fixing introduces relative phases

$M$ 2

Compact-gauge invariance eliminates one redundant degree of freedom, leaving $M$ 3 physical relative phases, forbids any term sensitive to the uniform gradient $M$ 4, and enforces vanishing total mass current identically in the long-wave theory (Golic et al., 2024).

The constitutive relation follows from variation with respect to $M$ 5,

$M$ 6

hence

$M$ 7

Counterflow currents remain finite and are driven by phase-gradient differences. At long wavelengths, the hydrodynamic equations take the superfluid form

$M$ 8

with physical Josephson relations written for the relative phases,

$M$ 9

The gauge redundancy removes the uniform $(1,1,\ldots,1)$ 0 mode and leaves $(1,1,\ldots,1)$ 1 gapless counterflow modes. In the symmetric case, one such mode has stiffness $(1,1,\ldots,1)$ 2, while the remaining $(1,1,\ldots,1)$ 3 modes have stiffness $(1,1,\ldots,1)$ 4; with suitable compressibility terms, all acquire linear dispersion $(1,1,\ldots,1)$ 5 (Golic et al., 2024).

3. Vortices, modular charges, and algebraic order

Away from vortex cores, extremal configurations obey Laplace equations, and localized topological sectors can be written as

$(1,1,\ldots,1)$ 6

The integer vector $(1,1,\ldots,1)$ 7 is defined modulo the uniform vector $(1,1,\ldots,1)$ 8, so the gauge-invariant content lies in relative charges

$(1,1,\ldots,1)$ 9

This modular arithmetic expresses the fact that equal winding of all phases is a gauge defect with no long-distance energy cost (Golic et al., 2024, Babaev et al., 2023).

In two dimensions, or per unit length in three dimensions, the vortex energy is

$M\ge 3$ 0

with

$M\ge 3$ 1

normalized so that elementary vortices have $M\ge 3$ 2. Pair interactions form a multicomponent Coulomb gas. Intracomponent interactions are logarithmically repulsive,

$M\ge 3$ 3

whereas intercomponent interactions are logarithmically attractive and weaker by a factor $M\ge 3$ 4,

$M\ge 3$ 5

This vector-charged plasma is qualitatively different from the scalar Coulomb gas of the standard XY model (Golic et al., 2024).

The basic off-diagonal correlator is the composite density matrix

$M\ge 3$ 6

which is symmetry-independent of $M\ge 3$ 7 and decays algebraically in two dimensions,

$M\ge 3$ 8

At the BKT point,

$M\ge 3$ 9

which is larger than the single-component value $U(1)$ 0 for all $U(1)$ 1, directly reflecting the Borromean universality class (Golic et al., 2024).

A later dynamical field theory reformulates the same defect structure through inter-species bosonic fields $U(1)$ 2. In that description, condensation breaks $U(1)$ 3 to $U(1)$ 4, and the ordered phase supports $U(1)$ 5 distinct flavors of energetically stable elementary vortex solutions even though

$U(1)$ 6

Thus only $U(1)$ 7 vortex generators are topologically independent, but $U(1)$ 8 elementary vortex flavors remain locally stable and are related by modular arithmetic (Kuklov et al., 29 Jul 2025).

4. Universality class and phase transitions

The two-dimensional finite-temperature transition retains the Nelson–Kosterlitz universal jump in form,

$U(1)$ 9

but the stiffness here governs the counterflow sector rather than total superflow. The renormalization-group flow of the Borromean vortex gas can be written in Kosterlitz–Thouless form,

$N\ge 3$ 0

provided the effective pair fugacity is rescaled by the Borromean factor

$N\ge 3$ 1

This factor, together with the counterflow exponent $N\ge 3$ 2, is what distinguishes Borromean criticality from both single-component and two-component counterflow systems, even though the flow of $N\ge 3$ 3 is formally identical after the fugacity rescaling (Golic et al., 2024).

At criticality, the finite-size stiffness obeys

$N\ge 3$ 4

and the correlator acquires the corresponding logarithmic corrections. The two-dimensional case is especially instructive because the Borromean nature of the system is strongly manifested while still admitting an asymptotically exact analytic description (Golic et al., 2024).

For dimensions greater than two, a later Landau-Ginzburg field theory based on the inter-species order parameters $N\ge 3$ 5 predicts that the transition into the Borromean super-counterfluid state is generically first-order. The mechanism is the presence of cubic invariants,

$N\ge 3$ 6

which phase-lock triplets and tend to drive a first-order transition in $N\ge 3$ 7 (Kuklov et al., 29 Jul 2025). This is consistent with the earlier observation that three-dimensional criticality may be nontrivial and that compact-gauge invariance likely plays a crucial role (Golic et al., 2024).

5. Statistical models, lattice formulations, and dynamical field theory

A compact-gauge-invariant lattice realization is provided by the Borromean XY model,

$N\ge 3$ 8

Unlike the standard XY Hamiltonian, this depends only on differences of phase gradients across components and is invariant under local shifts $N\ge 3$ 9. Fourier expansion yields an integer-current or loop representation,

$N\ge 3$ 0

on every bond. Local current conservation and bondwise arrest of the algebraic current sum are the discrete counterparts of compact-gauge invariance and vanishing net flow (Golic et al., 2024).

The minimal $N\ge 3$ 1 Borromean loop model allows only empty bonds, with weight $N\ge 3$ 2, or bonds carrying two counterflows, with weight $N\ge 3$ 3. Worm-algorithm simulations in this representation measure the composite density matrix and the stiffness via winding-number fluctuations. The numerical results validate the analytic theory in several ways: the aspect-ratio dependence of $N\ge 3$ 4 at $N\ge 3$ 5 agrees with theory without fitting parameters; the BKT critical exponent is recovered, with $N\ge 3$ 6 for $N\ge 3$ 7; and joint fits of $N\ge 3$ 8 and $N\ge 3$ 9 near criticality give

$N-1$ 0

while the rescaled pair density $N-1$ 1 remains small, validating the KT flow (Golic et al., 2024).

A complementary continuum formulation introduces Hermitian inter-species fields $N-1$ 2, with $N-1$ 3 and $N-1$ 4, transforming as

$N-1$ 5

Condensation of these fields breaks $N-1$ 6 to $N-1$ 7, so only relative phases are ordered. In the phase-only limit, the gradient energy becomes

$N-1$ 8

with $N-1$ 9 possessing the null vector $N$ 0. The hydrodynamic equations are

$N$ 1

and the linearized relative-phase sector yields $N$ 2 sound modes through the generalized eigenproblem $N$ 3 (Kuklov et al., 29 Jul 2025).

This dynamical field theory also predicts a counterflow AC Josephson effect. For a species pair $N$ 4,

$N$ 5

so a chemical-potential bias drives antisymmetric AC currents with zero net current (Kuklov et al., 29 Jul 2025).

The hydrodynamic framework extends beyond the component-symmetric case by replacing the single stiffness parameter with a symmetric matrix $N$ 6 satisfying $N$ 7 and positivity on the subspace orthogonal to the uniform vector. The action remains

$N$ 8

and the currents remain $N$ 9 with $M-1$ 00. Diagonalization on the $M-1$ 01-dimensional relative-phase subspace yields a spectrum of stiffness eigenvalues that modifies mode velocities, vortex energetics, and critical exponents. Asymmetry splits degeneracies but preserves the Borromean character: there is still no uniform superflow mode, and the vortex plasma remains vector-charged (Golic et al., 2024).

Off-diagonal intercomponent couplings can convert the Borromean supercounterfluid into a Borromean insulator. Adding $M-1$ 02-periodic even functions $M-1$ 03 to the hydrodynamic Hamiltonian pins all relative phases if the resulting mass matrix is positive-definite on the relative subspace, thereby gapping all $M-1$ 04 counterflow modes. When the minima occur at nontrivial phase differences, the pinned state can break time-reversal symmetry. For $M-1$ 05, one explicit chiral pattern is

$M-1$ 06

which supports domain walls carrying persistent counterflow textures (Babaev et al., 2023).

The principal experimental settings discussed in the literature are multicomponent Bose gases in optical lattices or the continuum, ring traps or annuli, two-dimensional films, and materials with frustrated Josephson couplings and broken time-reversal symmetry. The robust signatures are absence of net superflow, dissipationless counterflow transport, $M-1$ 07 gapless counterflow phonons, quantized counterflow circulation with vortex charges defined modulo the uniform vector, and the logarithmic force pattern of repulsive intracomponent versus attractive intercomponent vortex interactions. In elongated two-dimensional geometries, supercurrent states modify the phase-twist response and produce measurable aspect-ratio dependence of $M-1$ 08 (Golic et al., 2024). Species-selective weak links provide a direct route to the counterflow AC Josephson effect (Kuklov et al., 29 Jul 2025).

A related many-body manifestation appears in three-component ultracold Bose gases through the “Borromean droplet,” a self-bound state that exists only for the full ternary mixture while binary subsystems remain unbound. In that setting, the ternary induced attraction

$M-1$ 09

is stronger than the binary induced attraction

$M-1$ 10

and Gross–Pitaevskii hydrodynamics augmented by the Lee–Huang–Yang term governs continuity, Euler dynamics, compressibility, and multi-component sound within the droplet regime (Ma et al., 2021). This suggests a broader family of Borromean many-body phenomena in which transport, binding, or stability emerge only at the genuinely three-component level.

7. Alternative geometric usage of the term

A mathematically distinct usage appears in the work “Fluid dynamics as intersection problem,” where “Borromean hydrodynamics” denotes a triple interplay of equation of state, spacetime geometry, and differential topology rather than multicomponent counterflow order (Nekrasov et al., 31 Dec 2025). The phase space is the infinite-dimensional symplectic manifold

$M-1$ 11

with canonical symplectic form

$M-1$ 12

The coisotropic constraint manifold is

$M-1$ 13

while the equation of state and background metric define a Lagrangian submanifold through the generating functional

$M-1$ 14

On-shell fluid configurations are the intersections

$M-1$ 15

In this language, topology and symmetry reside in the coisotropic moment-map constraints, geometry enters through $M-1$ 16 and $M-1$ 17, and thermodynamics enters through $M-1$ 18 and its Legendre transform (Nekrasov et al., 31 Dec 2025).

This framework also accommodates anomalies and inflow. A deformation of the symplectic structure yields the anomalous flux

$M-1$ 19

and corresponds to a five-dimensional Chern–Simons action

$M-1$ 20

In the authors’ terminology, the structure is “Borromean” because nontrivial transport effects, such as the anomaly-induced flux $M-1$ 21 and the rocket term $M-1$ 22, require the simultaneous presence of topology, geometry, and thermodynamics; removing any one sector trivializes the effect (Nekrasov et al., 31 Dec 2025). This usage is conceptually separate from Borromean supercounterfluid hydrodynamics, but both employ the term to denote phenomena that exist only through an irreducible three-way structure.