The Vacuum’s Superfluid Lattice

Particles as speed-of-light vortices sealed by counter-rotating shear layers: one boundary model, from the electron to the cosmos

Author
Affiliation

Jeff Vroom

Independent

Published

November 8, 2026

Abstract

I found a model that fits broadly, in which the vacuum is a superfluid condensate of a light (~2 meV) particle and every particle is a vortex in it whose rotation is set by the signal speed c. Two adjustments to established superfluid-vacuum physics are needed: the Bohm quantum potential is the reaction force of a counter-rotating boundary layer (following Simeonov’s two-fluid derivation), and the photon is a modon, a self-propelled counter-rotating vortex pair. A motivated choice for the modon’s lowest Bessel mode in one cell, together with the measured dark-matter density, fixes a lattice cell of ξ ≈ 97 μm with no continuously fitted parameters. Measured in units of an effective quantum set by the electron mass and the Weinberg angle, the same cell gives the Higgs vacuum expectation value as 246.1 ± 0.3 GeV, within the 0.13% uncertainty carried by the dark-matter density, and the MOND scale follows as c√(Gρ_DM) to within 3% of McGaugh’s estimate. The model recovers the Painlevé–Gullstrand form of the Schwarzschild metric and sees mass as rotational energy leaking through quantized boundaries. Inputs, open factors, and Python scripts reproducing the main numbers are included.

Introduction

In 1861 Maxwell built electromagnetism on a sea of molecular vortices, and in 1867 Thomson proposed his vortex atom. Neither had a medium that could supply inertia or keep atoms stable. A self-binding superfluid fills those gaps, and the idea that empty space is such a superfluid has since become an established thread of research. This paper continues that thread and adds two adjustments:

  1. The Bohm quantum potential is the reaction force of a superfluid’s counter-rotating boundary layer, building on Simeonov’s two-fluid derivation (Simeonov 2025). His two fluids map directly onto the two components of a superfluid: the co-rotating bulk is the condensate, and his counter-flowing second fluid is the counter-rotating component of the Hall–Vinen–Bekarevich–Khalatnikov (HVBK) mutual-friction force.
  2. The photon is a modon (Figure 1) in that same fluid, using the vortex-dipole solutions of Larichev and Reznik (Larichev and Reznik 1976): the self-propelled pair of counter-rotating vortices that oceanographers and atmospheric physicists already know from Gulf Stream rings and the Madden–Julian oscillation. A modon carries energy and momentum but transports no net mass, moves at the speed its medium allows, and needs a matching energy to absorb or disrupt — every property a photon must have. A modon that travels at the speed of light requires a substrate whose vortices turn at that speed.

Both adjustments come from the shear-layer dynamics of frictionless vortices moving at the speed of light, and together they turn established superfluid physics into equations that span quantum mechanics and general relativity. The boundary layer behaves like a superfluid rotating in a cylinder, described by helium’s mutual-friction equations, which model lossless boundary interactions through a reactive and a dissipative channel. Combine that with the smallest modon the fluid can carry (set by Planck’s constant), and a single length falls out: the lattice cell size \xi. The measured dark-matter density fixes it with no continuously fitted parameters,

\xi_\mathrm{CP}=\left(\frac{\hbar}{\rho_\mathrm{DM}\,c}\right)^{1/4}\approx112\;\mu\mathrm{m},

and the geometry gives the cell occupancy, the fraction of one wavefunction each cell holds once neighboring cores overlap: f=4\pi/(K\sqrt2)=0.5666. The 4\pi lets the lattice gravitate across the full solid angle, the 1/K comes from the lowest Bessel-function mode of a modon held in one cell, and the 1/\sqrt2 comes from the balance that sets the size of a superfluid vortex core. Using this fraction instead of 1 gives \xi=\xi_\mathrm{CP}f^{1/4}\approx97\,\mum (Equation 11).

Particle physics reaches the same size from the other direction. The cell is the Compton wavelength of a ~2 meV vacuum particle, nicknamed dc1. The effective quantum m_\mathrm{eff}=1.70 MeV is one unit of vortex energy shared by the electron and the proton, fixed by the electron mass and the Weinberg angle. The cell is \nu=m_\mathrm{eff}/m_1\approx8.35\times10^8 of its reduced Compton wavelengths wide. This condensation number gives the Higgs vacuum expectation value. With the boundary geometry and the measured dark-matter density, the headline results are:

Quantity Framework Observed Off by Inputs
Higgs vacuum expectation value 246.1\pm0.3 GeV 246.22 GeV -0.04\%\pm0.13\% (inputs only) \rho_\mathrm{DM}, m_e (pole), \sin^2\theta_W (\overline{\mathrm{MS}}, M_Z) via \nu
MOND scale c\sqrt{G\rho_\mathrm{DM}} 1.16\times10^{-10} m/s² 1.20\times10^{-10} m/s² \sim3\% from McGaugh (20% systematic) \rho_\mathrm{DM}, G
Koide relation, read geometrically A=\sqrt2\;\Leftrightarrow\;Q=2/3 0.666664(5) within the m_\tau error lepton masses

The \approx97\,\mum lattice fills the space between matter, transmits light without loss, has no hard spectral floor, and couples seamlessly to chemistry.

The framework builds on five established results:

  • Volovik, The Universe in a Helium Droplet (Volovik 2003) — superfluid ^3He reproduces fundamental properties of the vacuum at both atomic and galactic scales, including an emergent relativistic (Dirac) spectrum.
  • Zloshchastiev and collaborators (Zloshchastiev 2020; Avdeenkov and Zloshchastiev 2011) — a self-binding condensate whose logarithmic equation of state (the superfluid-vacuum model) has a Lorentz-invariant ground state.
  • Bush, Oza and collaborators (Bush and Oza 2020) — a droplet bouncing on a vibrating bath, a particle in self-generated resonance with a responsive substrate, experimentally reproduces quantized orbits, tunneling, and single-particle diffraction.
  • Simeonov (Simeonov 2025) — the Madelung/Bohm quantum potential follows from the material derivative of the osmotic velocity in a two-fluid model.
  • Khoury and Berezhiani (Khoury 2015; Berezhiani and Khoury 2015) — if dark matter is a superfluid, its phonons mediate a MOND-like force on galactic scales, with a sharp transition to ordinary cold-dark-matter behavior above a critical velocity.
Two counter-rotating vortex cores stacked vertically inside one dashed closed envelope; the upper core is +omega counter-clockwise, the lower is -omega clockwise. A bold horizontal arrow to the right marks self-propulsion at speed c, and exterior streamlines loop from front to back, illustrating zero net mass transport.
Figure 1: The photon as a modon. A modon is a bound pair of counter-rotating vortex cores — here +\omega (counter-clockwise) above -\omega (clockwise) — sharing a single closed envelope of streamlines of scale \xi. Each core sits in the velocity field of the other, so the pair advects itself forward (the induced jet between the cores) and travels through the substrate at the signal speed c without the substrate doing net work. Because the two spins are equal and opposite, the substrate mass each core would carry cancels exactly (\int\psi\,\mathrm{d}A=0): the modon transports no net mass, propagates at c carrying energy E=h\nu and momentum p=E/c — the properties a photon must have.

Understanding this model

This section gives an intuitive picture before the math of §1.

The lattice size, \xi \approx 97\;\mum, feels at first glance far too large for the structure that fills the vacuum. Picture \xi not as a rigid grid but as the perturbation envelope of a stiff superfluid spinning at \approx 0.776\,c (the substrate’s own inner rotation speed, c\sqrt{2\alpha_{mf}}, §1.2). It mimics a vacuum but is nothing like one. It is an energetic, closely packed, strongly coupled condensate formed from balanced vortices breathing in anti-phase, tied together by topologically protected vortex lines. Its particles move too fast and sit too close to be resolved as individuals; instead they form lines connected to a vortex core, which disconnect and reconnect. As disconnected lines they bounce and rebound at the Compton wavelength. With an anti-phase partner they form counter-spinning nested eddies, while retaining a quantized momentum that returns them to the connected state along a path that is very hard to break. That stiffness is what lets them bend quickly enough to transmit modons — radiation without energy loss. Counter-spinning boundary layers form around the electron’s orbit and keep it from collapsing. Those same stiff layers limit the leak we feel as gravity, and between them the stream falls.

The stiff substrate also provides the energy that supports photons as self-advecting modons. Neutrinos are untwisted vortices wrapped in a modon dressing and move the same way. Everything else pushes a bow wave and trails a balancing stern wave. At steady speed the two cancel, and the particle coasts without drag, so inertia behaves just as it would in a true vacuum. Only near the signal speed, c, does the bow wave stiffen into ram pressure.

The lattice cell size is large for a simple reason: through the Compton wavelength, lighter particles have longer wavelengths, and dc1 appears very light because the cell’s leak is mostly contained. Its cell vortex’s lines start crowded at a core, disconnect, and oscillate outward at \sim0.776\,c — frictionlessly, so every collision is perfectly elastic — before reconnecting a wavelength later. A lighter, less dense core sends its ripples farther per cycle and overlaps more deeply with its anti-phase breathing partner, so the interior pairs need a larger envelope before a containing counter-spinning boundary can close around them.

That envelope is measured by the condensation number \nu=m_\mathrm{eff}/m_1\approx8.35\times10^8. This one number shows up in two ways. As the visible mass ratio, \nu is the number of dc1 quanta that condense in concert to make one effective quantum — the 1.70 MeV vortex unit inside every particle. As a length count, the same number is the Compton lift \xi/\bar\lambda_\mathrm{eff} — how many effective-quantum reach-lengths (116 fm) fit across one lattice cell. That sheer size is where the vacuum’s energy hides.

All of these disconnecting and reconnecting vortex lines make the substrate an energetic seamless glue, underlying everything and imposing a large, blurry window on observation. It hides further still by nesting inside stable, balanced layers, concealing a subtle scale-invariant organization formed by the substrate energy.

Next comes the substrate’s reading of the Standard Model and how it connects to cosmology.

The substrate overview

The framework posits a new particle, nicknamed dc1, of mass m_1\approx2\;\mathrm{meV}/c^2. It condenses into a strongly coupled superfluid on the BEC side of the BCS–BEC crossover and forms the quantized vortices that carry the substrate’s internal weather. Because its self-interaction is logarithmic, it also self-binds into a lattice of cell size \xi.

These identifications follow:

  • The vacuum is the dc1 condensate, organized into a lattice of cell vortices, six stacked anti-phase pairs per cell. Its rest-frame density is the cosmological dark-matter density, n_1 m_1 = \rho_\mathrm{DM}.
  • Mass is the held rotational kinetic energy, seen through the fraction of a vortex’s spin that escapes its outermost counter-rotating boundary, the rest staying trapped and unweighable (§1.2).
  • A fermion (electron, quark) is a persistent particle vortex of dc1, a topologically protected excitation, with an odd number of counter-rotating boundary layers separating its co-rotating interior from the surrounding substrate. The odd number of layers leaves it unbalanced, leaking mass-energy as vortex lines. The heavier generations (muon and tau) carry more of these wrapping layers.
  • A boson (photon, W, Z) has even boundary parity, a balanced counter-rotating pair. The photon is the fully balanced case, a modon with zero net momentum and no net leak of vortex lines, hence massless. It forms when an electron drops to a lower orbital and its two boundary layers roll up together, and the plane of the pair gives light its polarization.
  • Antimatter is opposite-spin vortices. When matter and antimatter meet they roll up into a modon, stop leaking, and appear massless, their angular momenta canceling.
  • A proton is a Borromean ring, three locked vortices. Colors come from topological interlocking phases, and flavors from braids of the knot.
  • Charge is the twist of vortex strands, following the helon model (Bilson-Thompson 2005). Each quark has three strands; the up quark has two of them twisted and the down quark only one. The proton’s vortex topology combines them into a charge of +1 at its boundary, while the twists cancel for the neutron.
  • An electron and neutrino are the twisted and untwisted vortices spun out to conserve topology when a proton or neutron changes.
  • Spin statistics come from counter-rotating boundary-layer angular momentum: the two-layer gyroscope of a fermion gives a double-cover topology, hence the 720° return.
  • The CPT theorem is the conservation of topologically protected vortex lines: P flips a particle’s motion relative to its spin, and T reverses the flow direction itself rather than running time backwards.
  • Symmetries are preserved through one continuous system of equations built on the topology of vortex cores and their boundary layers.

An illustrated dictionary of the framework on a white background, a header naming the one species dc1 of mass about 2 millielectronvolts, then six framed tiles. The vacuum: a triangular lattice of small blue co-rotating cells, one per site, with the cell spacing marked xi and the caption that its rest density equals the dark matter, n1 m1 = rho_DM. Mass: a single blue co-rotating spinner labelled plus omega, net spin uncancelled, bleeding three grey arrows of rotational kinetic energy outward to a tag reading leak rate equals m. A fermion: a blue co-rotating core wrapped in one red dashed counter-rotating layer, a dashed arrow entering as plus psi and flipping to minus psi inside, marked one layer, odd, and (minus one) to the first power equals minus one, antisymmetric, Pauli. A boson: a blue plus-omega core stacked over a red minus-omega core inside one dashed envelope, an induced jet and a bold v equals c arrow, even parity (minus one) squared equals plus one, net mass zero. A proton: three interlocked Borromean rings colored blue, green and red, labelled u, u and d, locked around a small gold central core, noted as three colors making one singlet, remove one and the rest fall free. Spin: an outer blue and inner red counter-rotating ring track with a big 720 degrees to return at center, one turn 360 degrees flips psi to minus psi and two turns 720 degrees return psi to plus psi, spin one-half. A footer shows the gear analogy: two blue plus-omega gears flanking a small red minus-omega idler, with text explaining that two like-spinning regions force a counter-spinning layer between them, that every red boundary above is one of these, and that counting them gives parity while locking three gives the proton, beside a legend keying blue to co-rotating plus omega and red to counter-rotating minus omega.

The same boundary model helps explain boundary problems at all scales. One boundary response function, the mutual-friction model of superfluid helium rotating in a cylinder, gives the Higgs VEV, a geometric reading of the Koide lepton relation, the charged-lepton mass ratios from a single phase angle, and a relation for the fine-structure constant that the measured couplings satisfy at the electroweak scale (§5).

Gravity is the rare leak through a particle’s outer rim, which then free-falls to the next boundary. Where the leaks from two masses meet, they form a self-binding channel that exerts gravity’s force. At low accelerations the lattice boundary responds quadratically, which gives the MOND scale a_0=c\sqrt{G\rho_\mathrm{DM}}; above the outer-rim speed, as in colliding clusters, that response breaks down. a_0 tracks the substrate’s mean density, which was higher in the early universe (the lattice scale itself is fixed), so the effect was stronger then, explaining galaxies that formed too early. Lensing and time dilation are effects of the flow and of how matter crosses lattice cells: moving fast through the substrate means more cell crossings per second, which slows the Compton clock, not time itself.

On cosmological scales, the Big Bang is Our Big Bubble, a bubble nucleated in a substrate that boils. The CMB comes from opposite-spin matter and antimatter vortices rolling up into modons until all the antimatter is gone. The evenly spaced harmonics of the CMB’s acoustic peaks are the ripples of density variation in the newly nucleated bubble, frozen after the plasma’s sound-speed transition. At recombination the photons rebound off those ripples and record a snapshot of the frozen pattern.

The paper builds the picture in four steps:

  1. The quantum potential as the reaction force of the counter-rotating boundary layer.
  2. The photon’s minimum energy as a modon, which fixes the bridge equation.
  3. Gravity, as a slow leak of substrate through the same boundaries.
  4. The lattice’s vertical geometry, and the patterns of underlying organization it reveals, including how chemistry shapes and uses its boundary layers.

Closing sections present the framework’s predictions, Python scripts to reproduce the numbers, and its open problems.

1 Counter-rotating layers and the quantum potential

In a frictionless substrate of vortex lines rotating around vortex cores at high speed, shear layers form at quantized boundaries. They organize the leak and seal it in: a wrapping layer forms once a quantum of leaking energy is large enough to support one. The lowest-energy way to absorb the shear of speed-of-light vortex lines is a thin counter-rotating boundary layer of eddies (Figure 2).

These two parts of the boundary — a co-rotating bulk and a counter-rotating correction — are the two-fluid structure of HVBK superfluid hydrodynamics and the structure Simeonov’s model needs.

Three parts on a white background. Left, a cross-section: a blue co-rotating vortex core marked plus omega counter-clockwise, sitting inside a faint dashed circle of resting substrate; at the core's rim a translucent red annulus holds six small clockwise minus-omega rollers, labelled the counter-rotating skin or shield (fluid 2) around the co-rotating core (fluid 1), captioned a shielded vortex with net circulation zero. Center, a magnifier lens zooms into one point on the rim: the seam is drawn as a horizontal dashed line with a grey still region above and a blue bulk-flow arrow pointing left below, and a single red minus-omega roller straddling the seam, captioned two speeds, one seam makes a roller spun the opposite way to the bulk. Right, a graph of azimuthal speed v-theta against distance from center: a blue straight line rises through the core to a peak at the rim (vorticity omega equals plus two Omega), then a red curve falls back toward zero faster than a dashed slate 1/r reference curve labelled the unshielded tail with omega equals zero; an annotation notes that where v-theta drops faster than 1/r the vorticity flips sign to omega less than zero, the skin, and the net circulation returns to zero.
Figure 2: The shear zone: where counter-rotation comes from. Left: around a single co-rotating vortex core (+\omega, blue), a thin counter-rotating skin (-\omega, red) forms to balance the energy. Center: zoom into the seam — two regions in relative motion put a single roller between them, spun the opposite way to the bulk. Right: the velocity profile shows the azimuthal speed v_\theta rises through the core to a peak at the rim, then drops back toward zero faster than the 1/r tail an unshielded vortex would keep. Wherever v_\theta falls faster than 1/r the vorticity flips sign. This is the mechanism behind the quantum potential, and core to core it is the photon.

Let fluid 1 (co-rotating, the “particle”) have density \rho_1 and velocity \mathbf v_1, and let fluid 2 (counter-rotating boundary) respond to gradients in fluid 1 with an osmotic velocity

\mathbf v_2 = -D\,\nabla \ln\rho_1, \qquad D = \frac{\hbar}{2m}, \tag{1}

where D is the diffusion constant of the counter-rotating layer. Simeonov (Simeonov 2025) shows that the reaction force of fluid 2 back on fluid 1 is, with R=\sqrt{\rho_1},

\mathbf F_\mathrm{reaction} = -\nabla Q, \qquad Q = -\frac{\hbar^2}{2m}\,\frac{\nabla^2 R}{R}, \tag{2}

which is the Bohm quantum potential.1 It is the back-pressure of the counter-rotating boundary on the curvature of the co-rotating density, shown for the hydrogen orbital in Figure 3. Three behaviors follow: near a density maximum \nabla^2 R<0 so Q>0 — a repulsive “quantum pressure” that keeps the electron from spiralling in; near a node Q dips sharply negative and its gradient pushes flow away from the node, maintaining it; in a uniform region \nabla^2R=0 and Q=0, no boundaries, no quantum effects. This is the substrate content of Bush’s hydrodynamic analogy: the electron is a real structure resonating with a real responsive substrate, and the “pilot wave” is the dc1 flow that the counter-rotating layer organizes around it.

Left: a hydrogen 1s electron as a co-rotating blue vortex core (+omega) with the proton at center and density peaking there, ringed by a thin layer of small counter-rotating red eddies (-omega) that absorb the edge shear, labelled with the osmotic law v2 = -D grad ln rho1, D = hbar/2m. An arrow labelled 'reaction force F = -grad Q' connects to a graph on the right of energy versus radius: the quantum potential Q (red) rises steeply and is positive near the center, the Coulomb potential V = -e^2/r (slate) dives down, and their sum V+Q=E is a flat constant line, showing the outward quantum pressure exactly cancelling the inward Coulomb pull.
Figure 3: The quantum potential as a boundary-layer back-pressure. Left: a hydrogen 1s electron is a co-rotating vortex of dc1 (+\omega, blue) whose density \rho_1 peaks at the proton and falls off outward — a slice of the n=1 orbital. Where the co-rotating flow shears against the resting substrate, the cheapest response is a thin counter-rotating boundary layer of eddies (-\omega, red — fluid 2), which tracks gradients in \rho_1 with the osmotic velocity \mathbf v_2=-D\,\nabla\ln\rho_1, D=\hbar/2m. Right: the reaction force of that layer back on the bulk is the Bohm quantum potential Q=-\frac{\hbar^2}{2m}\nabla^2R/R (with R=\sqrt{\rho_1}). Near the density peak \nabla^2R<0, so Q>0: an outward “quantum pressure.” For the ground state it exactly cancels the inward Coulomb potential V=-e^2/r, leaving the sum V+Q=E constant — the substrate reason the electron neither spirals in nor radiates.

1.1 The effective quantum

Equation 1 used Planck’s constant in the diffusion constant D=\hbar/(2m). Here it comes back as the action of the counter-rotating boundary layer, and that action is tied to the effective quantum’s mass-energy.

The HVBK mutual-friction force is built with a dimensionless coupling \alpha_{mf} — the fraction of each boundary interaction that dissipates energy rather than deflecting the flow. For a superfluid with quantized circulation \kappa_q = h/m_\mathrm{eff}, the effective diffusivity of vortex-mediated transport is D_{sf} = \kappa_q/(4\pi\,\alpha_{mf}).2 Setting this equal to the quantum diffusion constant D=\hbar/(2m) of Equation 1 — where m is the mass of the particle the orbital builds (the electron, for hydrogen) — and substituting \kappa_q = 2\pi\hbar/m_\mathrm{eff} gives the central mass relation

\boxed{\;m_\mathrm{eff}\,\alpha_{mf} = m\;} \tag{3}

and, equivalently, \hbar = 2mD: Planck’s constant reframed as the action of the counter-rotating boundary layer. The quantum comes from the discrete minimum vortex circulation, m_\mathrm{eff}\kappa_q, where m_\mathrm{eff} is the mass of one effective quantum of dc1 circulation and \alpha_{mf} is the mutual-friction coupling for the relevant scale. With the electroweak value \alpha_{mf} = \tan^2\theta_W = \sin^2\theta_W/(1-\sin^2\theta_W) = 0.30083 and m=m_e, the electron mass, the effective quantum is m_\mathrm{eff} = m_e/\alpha_{mf} = 1.70\;\mathrm{MeV}/c^2, heavier than the electron by 1/\alpha_{mf}\approx3.3. Applying the same relation in the nuclear sector, with this one shared quantum, gives \alpha_{mf}^{(N)} = m_p/m_\mathrm{eff}\approx552. That number is not a larger leak fraction (no single boundary can leak more than half, §1.2) but a count: the proton presents N_p\approx1836 boundary seams, each leaking the electron’s same \alpha_{mf}=0.3008, so \alpha_{mf}^{(N)}=N_p\,\alpha_{mf}. It stores \sim1836 effective quanta (3.12 GeV) and shows 30\% of them, 938 MeV.4

1.2 E=mc^2 — vortex energy fully expressed

All particles are vortices spinning in place with ordinary rotational kinetic energy. Only part of that is visible, because the inner energy is trapped behind the counter-rotating boundary layer that forms in the region right around the particle. Only the leaked part shows up in observations. The leaking fraction is \alpha_{mf}: about 30\% of the effective quantum’s spin energy dissipates outward. The other \approx70\% stays reactive and invisible to weighing — but not invisible to a hard enough probe, where it shows up as the scattering phase and the magnetic moment (§5.4).

Together, though, the flywheel energy of the spin and the leak share one speed, the speed of the medium, c, and the rim is locked to a fixed fraction of it, v_\mathrm{rot}^2 = 2\alpha_{mf}c^2 (v_\mathrm{rot}=0.776\,c, the substrate’s own inner rotation speed). This is consistent with

\underbrace{\tfrac12\,m_\mathrm{eff}\,v_\mathrm{rot}^2}_{\text{flywheel energy stored}} \;=\; \tfrac12\,m_\mathrm{eff}\,(2\alpha_{mf}c^2) \;=\; \underbrace{(\alpha_{mf}\,m_\mathrm{eff})}_{\text{the part that leaks}}c^2 \;=\; m_e c^2 . \tag{4}

The \tfrac12 and the 2 cancel; \alpha_{mf} converts the effective quantum into the observed electron mass; what is left is E=mc^2. The c^2 is the flywheel’s rim speed, pinned to c because c is the fastest the substrate carries anything. And mass is not converted into energy: mass is that energy, seen through the finite aperture of a counter-rotating boundary. A nuclear reactor does not convert mass to energy; it tears open the stiff barrier that fixes the aperture, exposing the vast rotational energy that is already there.

The factor \alpha_{mf}=\tan^2\theta_W comes from the Weinberg angle, measured in particle accelerators, and the same per-boundary leak applies in the nucleus and in the electron’s inner core. It is close to as large as it can be: because \alpha_{mf}=\tfrac12\sin2\delta_0 (footnote to §1.1), no single boundary can leak more than half. The electron, with one boundary, lets out under a third of one quantum’s spin, which is why it is light. The proton gets heavy not by leaking harder but by presenting more boundaries — \sim1836 seams, each leaking the same 30\% and hiding the same 70\%. Merged boundaries are also why a bound nucleus weighs less than its parts: two merged boundaries present a smaller leaking surface than two separate ones.

The shared effective quantum m_\mathrm{eff} is fixed by Equation 3 from m_e and \alpha_{mf}. It will matter in Step 2, because the same circulation quantum \kappa_q = 2\pi\hbar/m_\mathrm{eff} sets the speed of light and the modon’s action.

2 The photon as a modon, and the bridge equation

2.1 The modon and the emergent speed of light

A modon is a steadily propagating dipole of two counter-rotating vortex cores that advect each other forward: each core sits in the velocity field of the other, and their mutual induction drives the pair through the substrate without the substrate doing net work (Larichev and Reznik 1976; Saffman 1992). Its energy is entirely internal, and its net mass transport is zero. These are the properties a photon must have: massless, self-propelled, and energy E=h\nu. A modon in the ocean has friction and is still difficult to disrupt. The photon’s energy is held in topologically protected vortex lines, so disrupting it takes a matching counter-rotating boundary that can absorb or transform it.

Those same counter-rotating layers show how the atom emits one (Figure 4). A bound orbital is built from the standing wave’s counter-spinning shells. When the electron drops from n\to n', the released energy \Delta E=E_n-E_{n'} has to go somewhere: it thins the seam between two shells, that sheared layer rolls up into a counter-rotating pair, and the pair pinches off as a free modon carrying \Delta E=h\nu. The photon is born massless and moving at c, set by the orbital’s own rotation speed — each direction rotating at 0.776\,c.

Four panels left to right. I: a hydrogen orbital drawn as two concentric counter-spinning shells, inner +omega blue counter-clockwise and outer -omega red clockwise, proton at center. II: the atom drops from level n to a lower n-prime, releasing energy Delta-E, and the seam between the two shells thins at one spot where the shear concentrates. III: that sheared layer rolls up into a small counter-rotating blue-and-red vortex pair at the rim. IV: the pair detaches as a free modon with a blue +omega core above a red -omega core, an induced jet between them, propelling itself rightward at v=c and carrying Delta-E = h-nu, while the atom relaxes into a smaller n-prime orbital.
Figure 4: How an atom emits a photon: a modon pinches off. (I) A bound orbital is built of the standing wave’s counter-spinning shells — an inner co-rotating shell (+\omega, blue) and an outer counter-rotating shell (-\omega, red), the same two-fluid pairing that supplies the quantum potential. (II) When the electron drops n\to n', the released energy \Delta E=E_n-E_{n'} has nowhere to hide: the seam between the two shells thins where their shear concentrates. (III) That sheared layer rolls up, Kelvin–Helmholtz fashion, into an embryonic counter-rotating vortex pair. (IV) The pair pinches off as a free modon and self-advects at the substrate speed c, carrying \Delta E=h\nu away while the atom settles into the smaller n' orbital. The freed pair is the orbital’s own two counter-spinning shells, now detached — +\omega above, -\omega below — which is why a photon is born already massless (\int\psi\,\mathrm{d}A=0) and already moving at c.

In the strong-coupling BEC regime the dc1 quasiparticle spectrum is Dirac-like,

E^2 = \mu^2 + c^2 p^2, \qquad c = \frac{\hbar}{m_1\,\xi} \tag{5}

with a single isotropic speed set by the dc1 mass and the coherence length \xi (Volovik 2003).5 Every low-energy excitation — phonon, modon, and the tensor metric perturbation of Step 3 — inherits this one speed, so the framework predicts c_\mathrm{GW}=c exactly; gravity should travel at the speed of light, as LIGO measured with GW170817: |c_\mathrm{GW}/c-1| < 6\times10^{-15}. The speed of light is the maximum group velocity of organized disturbances in a superfluid of 2\;\mathrm{meV} particles with a hundred-micron coherence length.

The symbol \xi is the coherence length of the dc1 condensate, the perturbation envelope of the photon/modon. The relation comes from Volovik’s superfluid helium, with the speed of light as the medium speed. It is also the Compton wavelength of the dc1 particle, and because this condensate is close-packed, it is also the lattice cell size.

The modon has a precise solution from Larichev–Reznik (Larichev and Reznik 1976) that matches the interior and exterior boundary energies, and long-lived eddy pairs in the ocean and atmosphere show it in nature. Inside a critical radius the flow circulates and oscillates; outside it dies away exponentially into the surrounding substrate; and the two must join smoothly at the edge of the cell, r=\xi. Picture a self-contained swirl with a sharp skin — trapped circulation within, a fading halo without. That is the modon in flight: the photon.

A cell standing in the lattice is held differently. Its anti-phase partner cancels the dipole beyond the cell (§4.1), so the paper treats the cell’s edge as a node rather than a smooth join onto a halo. The interior pattern is a Bessel function, \psi\sim J_1(pr), for the same reason a circular drumhead vibrates in Bessel-function modes: it is the natural shape of a wave confined to a disk. A node at the edge, with the swirl completing exactly one lobe across the cell, pins the interior wavenumber to the first zero of J_1, the pure number j_{11}=3.8317 (the same constant that sets the fundamental dipole tone of a clamped circular drum). The medium adds its gap: the dc1 spectrum is gapped (Equation 5), so the mode obeys (-\nabla^2+\kappa^2)\,\psi=(K/\xi^2)\,\psi with screening wavenumber \kappa=m_1c/\hbar=1/\xi, the same rate at which any static disturbance outside a cell dies away. In units of the cell size the eigenvalue is the dimensionless structure constant

K = (p\xi)^2 + (\kappa\xi)^2 = j_{11}^2 + 1 = 15.68 . \tag{6}

K is a pure geometric number, the square of a Bessel zero plus one. The term j_{11}^2 is the lowest dipole mode of a clamped cell. The term +1=(\kappa\xi)^2 is the medium’s gap, fixed like the cell itself by the single scale \xi=\hbar/(m_1c), so the “+1” is a consequence of the cell size. For a clamped disk this closed form is exact: the gap adds exactly 1 to every clamped eigenvalue and leaves the mode unchanged. A free modon instead joins smoothly onto a halo that decays at the same rate, and its Larichev–Reznik matching gives K=16.39. The difference is a boundary condition, a node or a smooth tail. For the lattice, the node for K is selected from the structure and is supported by the predictions (§5.1, §5.2).

The same p^2+\kappa^2 is the denominator of the mode’s linear Rossby-wave dispersion, U=-\beta/(p^2+\kappa^2) with \beta the background vorticity gradient; setting that speed equal to c recovers the Volovik relation c=\hbar/(m_1\xi). This is an identity, not a new constraint: it restates why a counter-rotating dipole in the substrate travels at exactly the speed of light.

Two existence conditions follow. First, the matching has no solution below scale \xi: the smallest modon is exactly one lattice cell wide. Second, that smallest modon has wavelength \lambda=\xi; since a modon carries energy E=hc/\lambda, this fixes a minimum energy

E_\mathrm{min} = \frac{hc}{\xi} = 2\pi\,m_1 c^2 \approx 12.8\;\mathrm{meV}, \tag{7}

the minimum energy of a single localized modon — one whose whole dipole fits in one cell — at \lambda\sim100\,\mum, a frequency of \sim3 THz. This is a floor on one elementary modon, not a cutoff on the spectrum: lower-frequency light propagates as a coherent superposition of many cells, the collective limit of the same modon field (§5.1).6

How, then, does ordinary light live above this floor? Visible light (\sim2 eV) sits far above threshold, \lambda\ll\xi: it keeps the full \xi-scale envelope but concentrates its energy in a compact sub-cell dipole core — a boat in a harbor, a small hull trailing a basin-wide wake. At the floor the hull grows to fill the harbor, and core and envelope coincide at \xi. Below it (\lambda>\xi) the dipole can no longer localize inside one cell, and the excitation spreads into a delocalized, many-cell collective mode — the way radio, microwave, and the CMB propagate, each photon a single quantum of collective winding spread over many cells, just as a long-wavelength phonon is one quantum of collective motion far below any single atom’s vibration. So 12.8 meV (\sim3 THz) is a crossover in modon structure, not an edge in the spectrum: nothing goes dark there, only the shape of the packet changes.

2.2 The bridge equation

The next step is to find the lattice size from the dark-matter density, assuming one wavefunction per cell, and then the geometry that pulls it into place: the bridge equation.

Together they fix \xi, \nu, and the 2 meV dc1 particle behind the predictions of §5.

Route 1 — from cosmology. Three relations describe the dc1 superfluid: the Volovik speed Equation 5, the composition n_1 m_1 = \rho_\mathrm{DM}, and close-packing n_1\xi^3 \approx 1 (one condensate wavefunction per cell, the Bose-condensation degeneracy condition n\lambda_{dB}^3\sim1).7 Eliminating n_1 and m_1 leaves a single equation,

\xi_\mathrm{CP} = \left(\frac{\hbar}{\rho_\mathrm{DM}\,c}\right)^{1/4} \approx 111.8\;\mu\mathrm{m}, \tag{8}

with inputs \hbar, c, and the dark-matter density \rho_\mathrm{DM}=2.254\times10^{-27}\;\mathrm{kg/m}^3 from Planck 2018 (Planck Collaboration et al. 2020).8 (Equivalently n_1\xi^3=1 and Equation 5 reduce, with no numbers, to \rho_\mathrm{DM}\xi^3=m_1 — close-packing is the same statement as emergent light.)

Route 2 — the electroweak side, through the condensation number. From particle physics, three conditions — the modon constant, a gravitational self-consistency condition called SC2, and the Gross–Pitaevskii core balance — fix the cell occupancy, and with it the condensation number that the electroweak measurements land on independently.

SC2 is the requirement that the lattice not merely carry disturbances but gravitate — that a lump of substrate energy bend the emergent (acoustic) metric exactly as Newton says a mass should, \nabla^2\Phi = 4\pi G\rho (Step 3). The 4\pi comes from the surface area of a unit sphere, the Gauss solid angle over which flux from a point spreads in three dimensions. For the substrate’s vortices to source gravity with that coefficient, one circulation quantum must satisfy

\kappa_q\,\Omega_v = 4\pi c^2, \qquad \Omega_v = \frac{2m_\mathrm{eff}c^2}{\hbar}, \tag{9}

where \kappa_q=2\pi\hbar/m_\mathrm{eff} is the quantum’s circulation — its swirl, in m²/s — and \Omega_v is its Compton zitterbewegung clock, its internal rate in 1/s. The striking part is that their product is 4\pi c^2 identically: m_\mathrm{eff} cancels. So SC2 adds no adjustable dynamics — it is the statement that the same quantum which carries circulation also sources gravity, and the two roles can only lock together through the geometric 4\pi. Both sides have units m²/s², so SC2 is dimensionally clean and unit-invariant, and supports the cell occupancy equation’s leading 4\pi term (Equation 10).9 From the lattice’s parameters it reads \Omega_v\,\xi/c = 2\nu — a pure number set equal to a pure number.

The third ingredient is the Gross–Pitaevskii core balance. At a vortex core the condensate’s kinetic energy \hbar^2/(2m_1\xi^2) balances its interaction energy m_1c^2; equating the two fixes the healing length \xi_\mathrm{GP}=\hbar/(\sqrt2\,m_1c)=\xi/\sqrt2 (Fetter 2009). The \sqrt2 comes from the 2 in the kinetic operator \hbar^2/2m — an elementary fact of non-relativistic quantum mechanics.

The three assemble into a fraction. SC2 contributes its Gauss’s-law 4\pi; the cell’s modon mode contributes 1/K with K=j_{11}^2+1; the core balance contributes 1/\sqrt2. Their product is the cell occupancy, the number of dc1 wavefunctions the substrate actually puts in one cell:

\boxed{\;f = \frac{4\pi}{K\sqrt2} = 0.5666\;} \tag{10}

4\pi is a solid angle, K is a Bessel zero squared plus one, \sqrt2 is the 2 in \hbar^2/2m. This is the bridge equation. A possible fourth factor, a three-dimensional correction \eta, is set to 1, consistent with the randomly oriented domain texture of §5.3.

Route 1’s close-packing “\approx1” is the f\!=\!1 idealization: one wavefunction exactly filling each cell. But binding energy pulls the cores into overlap, and the modon envelope must close inside the same width. The geometry shows the real occupancy as f, and since \rho_\mathrm{DM}\,c\,\xi^4/\hbar=f makes \xi\propto f^{1/4}, the true lattice cell width is

\xi = \xi_\mathrm{CP}\,f^{1/4} = 111.8\;\mu\mathrm{m}\times0.8676 = 97.0\;\mu\mathrm{m}. \tag{11}

In a nutshell, the cell occupancy 0.5666 holds three conditions at once: heal a vortex core (1/\sqrt2, Gross–Pitaevskii), close the modon envelope (1/K, Bessel mode), and present the correct flux to gravity (4\pi, Gauss).

The other thing the electroweak side supplies is the effective quantum, m_\mathrm{eff}=m_e/\alpha_{mf}=1.699 MeV/c^2 (Equation 3), whose reduced Compton length \bar\lambda_\mathrm{eff}=\hbar/(m_\mathrm{eff}c)=116.2 fm is the ruler the cell is measured in. The one number left over is the count of rulers — the condensation number,

\boxed{\;\nu \equiv \frac{m_\mathrm{eff}}{m_1} = \frac{\xi}{\bar\lambda_\mathrm{eff}} = \frac{\omega_\mathrm{eff}}{\omega_1} \approx 8.35\times10^8\;} \tag{12}

— a ratio of two masses, two lengths, or two frequencies, and so a pure number.

2.2.1 The agreement

Three measurements from unrelated experiments constrain \nu: a cosmological density, a collider VEV, and a heliospheric flow. The first two legs also share the effective quantum, built from \sin^2\theta_W and m_e, so what each leg adds is its own measurement:

Leg Relation Measured input \nu
Cosmology \times geometry \nu=\xi_\mathrm{CP}f^{1/4}/\bar\lambda_\mathrm{eff}, f=4\pi/(K\sqrt2) \rho_\mathrm{DM} (Planck 2018), m_\mathrm{eff} (8.349\pm0.021)\times10^8
Electroweak \nu=v^2/(8\pi\,m_\mathrm{eff}^2c^4) v=246.22 GeV (collider), m_\mathrm{eff} (8.356\pm0.006)\times10^8
Heliospheric \nu=4\pi\,(c/v_L)^3 v_L=751.5 km/s (Ulysses) 7.98\times10^8

Cosmological density, collider electroweak physics, and heliospheric flow converge on one \nu\approx8.2\times10^8 to \sim5\%. The two tightest legs agree to 0.08\%, well inside the 0.25\% that the 1\% uncertainty in \rho_\mathrm{DM} puts on the cosmology leg (\nu\propto\rho_\mathrm{DM}^{-1/4}). The electroweak leg’s error comes from \sin^2\theta_W in m_\mathrm{eff}. The two sides use the effective quantum with opposite powers — \nu\propto m_\mathrm{eff} on the length side (through 1/\bar\lambda_\mathrm{eff}), \nu\propto m_\mathrm{eff}^{-2} on the VEV side — so m_\mathrm{eff} cannot cancel out of the comparison. Their agreement is a constraint tying \rho_\mathrm{DM}, v, and \sin^2\theta_W together.

3 Gravity as a leak, and the acoustic metric

3.1 The waterfall and the Schwarzschild metric

The substrate leaks through stiff counter-rotating boundary layers — an ebbing flow that only rarely breaks through a boundary’s rapids, but once through falls freely as a gravitational waterfall until it reaches the next boundary, every 97 μm. When streams from two opposing masses meet, they create a self-binding channel that exerts the gravitational force. The rare boundary crossing is felt as weight, while the free-falling waterfall between boundaries is the velocity field that shapes the metric, giving the illusion of curved spacetime.

Unruh, Visser, and Volovik showed that sound in a flowing fluid obeys an effective Lorentzian (acoustic) metric (Unruh 1981; Barceló et al. 2005; Volovik 2003).10 For a steady radial inflow v_\mathrm{ebb}(r) the acoustic line element is

ds^2 = -\!\left(c^2 - v_\mathrm{ebb}^2\right)dt^2 - 2\,v_\mathrm{ebb}\,dr\,dt + dr^2 + r^2 d\Omega^2 . \tag{13}

A dc1 particle falling from rest at infinity reaches radius r with v_\mathrm{ebb}(r) = \sqrt{2GM/r} from the substrate’s Euler equation (the convective acceleration v\,dv/dr = -GM/r^2 is the free-fall law, with zero pressure gradient, as the equivalence principle requires). Substituting into Equation 13 gives

ds^2 = -\!\left(1 - \frac{2GM}{rc^2}\right)c^2 dt^2 - 2\sqrt{\tfrac{2GM}{r}}\,dr\,dt + dr^2 + r^2 d\Omega^2 , \tag{14}

which is the Painlevé–Gullstrand form of the Schwarzschild metric — the “river model” of a black hole (Hamilton and Lisle 2008). Here the river is the substrate waterfall.11

Two panels. Left, the mechanism: three stiff red counter-rotating -omega boundary bands stacked toward a mass below, each with a narrow central gap; a blue substrate-ebb arrow descends, squeezing through each gap (the leak, felt as weight) and lengthening as it free-falls in the open space between bands (the waterfall, which shapes the metric), with the reading 'through a boundary, the rare leak felt as weight; between boundaries, free fall which shapes the metric.' Right, the consequence: a central mass M ringed by a bold horizon circle where v_ebb = c at r = 2GM/c squared; blue arrows converge inward from all sides as the river v_ebb(r) = root(2GM/r); along a ray to the right, black outgoing-light arrows grow with distance (net speed c minus v_ebb), shrink to near-zero just outside the horizon (marked frozen), and reverse to point inward just inside it (v_ebb greater than c); a boxed line element ds squared = minus (c squared minus v_ebb squared) dt squared minus 2 v_ebb dr dt plus dr squared plus r squared d-Omega squared is labelled the Painleve-Gullstrand form of Schwarzschild, the river taken literally.
Figure 5: Gravity as a leak, and the Painlevé–Gullstrand river. Left — the mechanism: the substrate ebb falls toward a mass through a stack of stiff counter-rotating -\omega boundary layers; it can only rarely transit a boundary — that leak is the force we feel as weight — while between boundaries it free-falls, the waterfall whose velocity field shapes the metric. Right — the consequence: coarse-grained, the inflow is a smooth radial river v_\mathrm{ebb}(r)=\sqrt{2GM/r} — free fall from rest at infinity, fixed by the substrate’s own Euler equation with no pressure gradient (the equivalence principle). Put into the acoustic metric it is exactly the Painlevé–Gullstrand form of the Schwarzschild line element, the river taken literally. The horizon is simply where the river reaches c: outgoing light, of net speed c-v_\mathrm{ebb}, slows to zero there and is swept inward beyond it.

These classical tests follow: the full nonlinear redshift \nu(r)/\nu_0 = \sqrt{1-2GM/rc^2}, light deflection \Delta\theta = 4GM/(bc^2), the Shapiro delay, and perihelion precession 6\pi GM/(ac^2(1-e^2)).12 The effect on atomic clocks, both gravitational and kinematic, is one mechanism — a clock’s signal runs at c relative to the substrate. Every way we measure time is the ticking of such a clock, and it runs on lattice-cell crossings per second. Movement through the substrate adds to that count and slows the clock — and every cell in your body — while time itself rolls steadily forward. With motion, the clock’s signals cross cells along a diagonal, so its ticks run slower by \sqrt{1-|\mathbf v_\mathrm{rel}|^2/c^2}, with \mathbf v_\mathrm{rel} its velocity relative to the local river.

At the dynamical level the linearized substrate equations reproduce the linearized Einstein equations, with the GR factor of 4 (the difference between 4\pi G and 16\pi G) supplied automatically because pressure gravitates: the stiff equation of state P=\rho c^2 gives an effective source \rho_\mathrm{eff}=\rho+3P/c^2 = 4\rho, and the substrate’s pressure is the rotational kinetic energy of its vortices.13 The full nonlinear completion proceeds order by order with coefficients fixed by the fluid equations; a formal all-orders proof is an open problem (§7).

Although the substrate is stiff locally, with both pressure and energy in its equation of state, that energy also balances locally. It affects how mass and radiation move through the substrate, but beyond that the anti-phase action cancels. Photons move masslessly, and massive objects form self-canceling disturbances. From observations at a distance it looks like a pressureless, temperatureless vacuum in space. These two properties, stiff up close and balanced at a distance, give accurate redshift predictions from a locally energetic substrate.14

3.2 The weakness of gravity

The river fixes the shape of gravity, and the gravitational constant sets the strength. The substrate relies on the measured G, expressed through a tiny crossing fraction:

G=\frac{\hbar c}{m_1^2}\,\mathcal F,\qquad \mathcal F\equiv\frac{Gm_1^2}{\hbar c}=\left(\frac{\ell_\mathrm{Pl}}{\xi}\right)^2\approx2.8\times10^{-62}. \tag{15}

\hbar c/m_1^2 is the coupling the substrate would have if its quanta passed the boundaries freely; \mathcal F is how far the counter-rotating boundaries hold real gravity below that. Writing the leak as a rate \Gamma and matching the units of G and the density leads to 4\pi G\rho_\mathrm{DM}=\Gamma^2. Measuring the rate against the cell’s Compton clock, \Gamma=f_\mathrm{cross}\,\omega_1 with \omega_1=c/\xi, and using close-packing \rho_\mathrm{DM}=f\hbar/c\xi^4 (the fixed cell density, which equals the cosmic mean today) gives the crossing probability per breath of the cell,

f_\mathrm{cross}=\sqrt{4\pi f}\;\frac{\ell_\mathrm{Pl}}{\xi}\approx4.4\times10^{-31}. \tag{16}

One dc1 quantum gets through about once in 10^{31} breaths.

The same number appears twice more. The MOND scale of §3.4 is the cell’s acceleration diluted by this probability, a_0=c\sqrt{G\rho_\mathrm{DM}}=f_\mathrm{cross}\,c^2/(\sqrt{4\pi}\,\xi). And with the Friedmann equation it gives the e-fold count N_*=\ln(c/H_0\xi)\approx69, compared with the \approx50–60 usually inferred for inflation. The cell is the geometric mean of the Planck length and the Hubble radius, \xi^2\approx4.2\,\ell_\mathrm{Pl}\,c/H_0, the long-noted dark-energy-scale coincidence. For now \mathcal F is taken from observation; deriving it is an open problem (§7).

3.3 The small cosmological constant

On cosmological scales the same model gives the Friedmann equations. From Volovik’s self-tuning (Volovik 2003) the equilibrium superfluid vacuum obeys the Gibbs–Duhem identity \varepsilon + P = 0 with \varepsilon as the vacuum energy density and P as the pressure. In full equilibrium both relax to zero — the vacuum is gravitationally inert, contributing nothing to the acceleration equation, with no fine-tuning. The observed cosmological constant is then a small residual disequilibrium, the vacuum prevented from fully relaxing by the ongoing expansion. The dark-energy scale is the substrate scale, \rho_\Lambda^{1/4}\approx m_1 c^2\approx2\;\mathrm{meV}, which dissolves the “why is \rho_\Lambda\sim\rho_\mathrm{DM} today” coincidence into a single number. The confounding mismatch with 10^{120} becomes an order-unity disequilibrium times the same weak-gravity hierarchy (m_1/M_\mathrm{Pl})^2=\mathcal F (Equation 15) that makes G small — reframing the cosmological-constant problem and Newton’s G as one problem.15 (The value of \Lambda remains an inherited initial condition.)

Two panels. Left, the mechanism: a wide red bar at the top marks the naive quantum-field-theory vacuum energy, rho_vac approximately ten to the 120 times rho_Lambda, the catastrophe. A bold blue arrow drains straight down, labelled epsilon plus P equals zero — the Gibbs–Duhem self-tuning that relaxes both epsilon and P to zero, so the ten-to-the-120 does not gravitate. The arrow drains into a tiny blue residual nub sitting on a baseline marked epsilon equals P equals zero, gravitationally inert; an arrow from the right notes that the ongoing expansion leaves this small residual rho_Lambda un-relaxed, and a box gives its scale, rho_Lambda to the one-quarter approximately m1 c squared approximately 2 millielectronvolts, the substrate scale itself. Right, why it looks tiny: a central pill holds the single residual disequilibrium delta-T over T-c and branches to two boxes. The blue box reads it against the substrate density rho_DM as the square root of rho_Lambda over rho_DM, about 1.6 — order unity, nothing tuned. The slate box reads the same residual against the Planck density rho_Pl as the square root of rho_Lambda over rho_Pl, equal to 3.4 times ten to the minus 62 — the notorious tiny Lambda. A gold bar beneath both states that the two readings differ by exactly the ruler ratio, the square root of rho_DM over rho_Pl, equal to m1 over M-Planck squared — the same hierarchy that makes G weak, so small Lambda and weak gravity are one number, not two.
Figure 6: A small cosmological constant. Left: in full equilibrium a superfluid vacuum obeys the Gibbs–Duhem identity \varepsilon+P=0 and self-tunes both to exactly zero (Volovik (Volovik 2003)), so the naive 10^{120} zero-point energy does not gravitate; the observed \Lambda is only the small residual the ongoing expansion leaves un-relaxed, and its scale is the substrate’s own, \rho_\Lambda^{1/4}\approx m_1c^2\approx2\,meV. Right — why it looks tiny: that one residual disequilibrium \delta T/T_c is order-unity (\approx1.6) seen against the substrate’s own density \rho_\mathrm{DM} — nothing tuned — and is the notorious 3.4\times10^{-62} only against gravity’s Planck density \rho_\mathrm{Pl}. The two rulers differ by exactly (m_1/M_\mathrm{Pl})^2, the same hierarchy (M_\mathrm{Pl}^2=\hbar c/G) that makes G weak — so a small \Lambda and weak gravity are one number, not two.

Through the lens of the substrate framework, the universe makes more sense as Our Big Bubble — a nucleation event in a substrate that boils. In this picture dark energy is the residual of the substrate’s incomplete relaxation after the transition (Volovik self-tuning); its time dependence is not developed in this paper.

3.4 MOND from the superfluid

Because the substrate is dark matter, the anomalous observations around galaxies are explained by a collective response of the substrate, not by the gravity of an unseen halo. Following Khoury and Berezhiani (Khoury 2015; Berezhiani and Khoury 2015), the substrate behaves as a coherent superfluid below the Landau critical velocity and as a collisionless gas above it. That critical velocity, v_L\approx2.5\times10^{-3}c\approx750 km/s, falls between galactic and cluster velocity scales, which is why rotation curves go MOND-like while clusters stay CDM-like.16 That MOND behavior is a property of the substrate: the counter-rotating boundary is parity-even — it holds +\omega_0 and -\omega_0 equally, the anti-phase breath of §4.1 — so its response has no linear term and is quadratic at leading order, which coarse-grained is exactly the Bekenstein–Milgrom deep-MOND field equation, the Newtonian limit appearing only once the Hubble flow’s DC bias breaks that parity. This reaches Khoury and Berezhiani’s superfluid dark matter from boundary symmetry rather than from a chosen P(X)\propto X^{3/2}.

This cosmological number uses only the dark-matter density that sets the lattice (Equation 8) and the measured G: the MOND acceleration scale is

a_0 = c\sqrt{G\rho_\mathrm{DM}} = 1.16\times10^{-10}\;\mathrm{m/s^2}, \tag{17}

matching the measured g_\dagger = (1.20\pm0.02_\mathrm{stat}\pm0.24_\mathrm{sys})\times10^{-10}\;\mathrm{m/s^2} (McGaugh et al. 2016) to \sim3\% — with no continuously fitted parameter.17 The same v_L\approx750 km/s is the terminal speed of the Sun’s fast polar wind — the polar jet of the Sun’s own loop (§4.4) saturating at that Landau velocity, a cross-scale match at the 1.5\% level (Ulysses mean 751.5 km/s) — one velocity seen at two scales. It follows from v_L=c\,(4\pi/\nu)^{1/3}=740 km/s (1.5\% low) using the same \nu that sizes the lattice cell and, through v=\sqrt{8\pi\,m_\mathrm{eff}^2c^4\nu}, the Higgs VEV, another coarse match for \nu.

This is a close reframing of Milgrom’s coincidence. He noted that the MOND scale is close to the speed of light times the expansion rate, which by the Friedmann equation scales as the square root of G times the total density. Here only the dark-matter density — the lattice density — enters, together with the signal speed c, and the result lands just 3\% below McGaugh’s measured scale.

That makes the MOND scale a property of the lattice: a boundary phase that holds in smooth substrate weather and is disrupted when flows are pushed past the outer-rim speed.

Two panels. Left, a vertical speed axis in kilometres per second crossed by a single horizontal gold threshold at the Landau critical velocity v_L equals omega-zero times xi, about 750 kilometres per second, which is also the Sun's fast polar-wind speed (Ulysses 751.5), marked by a gold dot riding the line. Below the threshold the band is shaded blue and labelled coherent superfluid leading to deep MOND; the Milky Way at less than about 200 kilometres per second and dwarf galaxies at 30 to 80 kilometres per second sit there as blue bars. Above the threshold the band is shaded red and labelled collisionless gas leading to CDM; galaxy clusters at 800 to 1500 kilometres per second sit there as a red bar. The caption reads: one substance, gated by velocity not position. Right, the radial-acceleration relation, a log-log plot of observed acceleration g_obs against the Newtonian acceleration g_bar from visible matter. A dashed grey one-to-one line is the Newtonian regime g_obs equals g_bar at high acceleration; the blue curve follows it at the top right but bends above it at low acceleration toward the deep-MOND form g_obs equals the square root of g_bar times a-zero. The transition scale a-zero is marked by gold dashed guide lines on both axes meeting at a gold knee dot. A boxed result in the empty upper-left states a-zero equals c times the square root of G times rho_DM equals 1.16 times ten to the minus 10 metres per second squared, against the measured g-dagger equals 1.20 plus or minus 0.02 times ten to the minus 10, agreeing to about 3 percent.
Figure 7: MOND from the superfluid. Left — one speed sorts the phases: the substrate stays a coherent superfluid below the Landau critical velocity v_L=\omega_0\xi\approx750 km/s and becomes a collisionless gas above it, so dwarf galaxies (30–80 km/s) and the Milky Way (\lesssim200 km/s) sit below it and go deep-MOND, while galaxy clusters (800–1500 km/s) sit above it and stay CDM — one substance, gated by velocity, not position; the same v_L is, strikingly, the Sun’s fast polar-wind speed (Ulysses mean 751.5 km/s). Right — the scale falls out: with today’s \rho_\mathrm{DM} fixed by the bridge equation, the radial-acceleration relation bends from the Newtonian g_\mathrm{obs}=g_\mathrm{bar} to the deep-MOND g_\mathrm{obs}=\sqrt{g_\mathrm{bar}\,a_0} at the scale a_0=c\sqrt{G\rho_\mathrm{DM}}=1.16\times10^{-10}\,\mathrm{m/s^2} — within \sim3\% of the measured g_\dagger=(1.20\pm0.02_\mathrm{stat}\pm0.24_\mathrm{sys})\times10^{-10} (McGaugh et al. 2016).

4 The shape of the lattice

The axis and plane of the modon give light its polarization. They also require each lattice cell to be organized as sheets of anti-phase oscillating pairs.

The cells themselves organize into a triangular (Abrikosov) lattice, held together by a residual quadrupole that also hides the leaking vortex lines between cells.

Together these give the energy of the vacuum a shape and a texture.

4.1 Sheets and the anti-phase breath

The lattice’s cell vortices organize into coherent sheets of aligned co-rotating vortices, where adjacent ones have the opposite Compton phase. This allows the vortex lines to almost fully close inside of each cell.

Some vortex lines run perpendicular to the sheets. The vertical period comes from the Glaberson–Johnson–Ostermeier instability wavelength of the vortex lines (Glaberson et al. 1974):

d_\mathrm{GJO} = \xi\sqrt{\frac{\ln(\xi/\xi_\mathrm{GP})}{4\pi}} \approx 0.166\,\xi \approx 16\;\mu\mathrm{m}, \tag{18}

with \xi_\mathrm{GP}=\xi/\sqrt2.18 Because the stack alternates handedness, like-handed sheets repeat at \approx16\,\mum but an anti-phase boundary falls every half-period, at \approx8\,\mum — the spacing for the boundary-locked structure. So those two values, 8 and 16\,\mum, form the substrate’s octave.

Each counter-rotating intermediate layer holds the anti-phase partner of the sheet before it. Together they bind like a Cooper pair in a superconductor — expansion and contraction in opposite phase, tied by a counter-spinning vortex in the middle. This anti-phase breathing trades energy across the shared seam without loss, and cancels outside the cell. The vacuum hums everywhere at \omega_1 = m_1c^2/\hbar\approx3\times10^{12}\;\mathrm{rad/s}, but the hum is silent at every scale we can directly observe.

Three panels. Left, in-plane: a triangular array of blue +omega vortex cores, each with a faint central density peak and a small open core ring, with the nearest-neighbour spacing marked xi approximately 100 micrometres on the edge of a highlighted triangular cell; a small Gaussian density-profile inset and a legend note that each cell is a self-bound gausson with no pinning species (held by dc1's logarithmic equation of state, width xi_GP = xi over root 2), all one handedness, a chirality-coherent sheet. Middle, vertical side view: blue +omega sheets and red -omega intermediate layers alternating along a vertical axis, threaded by vertical vortex-line filaments, with the like-handed period marked 16 micrometres (d_GJO) and the half-period 8 micrometres, the substrate's octave. Right, the seam: a large blue +omega core with outward arrows labelled 'expands' above a small red -omega core with inward arrows labelled 'contracts', separated by a dashed seam, with text explaining the anti-phase breath at omega_1 = m_1 c^2 / hbar, that the hum cancels above one cell leaving the infrared floor, and that the pair is a Cooper pair at rest and a photon in flight.
Figure 8: The shape of the substrate lattice. Left (in-plane): the lattice’s cell vortices form a triangular (Abrikosov) array of a single handedness — a chirality-coherent sheet — each a co-rotating +\omega vortex, one per cell, with nearest-neighbour spacing equal to the bridge-equation cell size \xi\approx100\,\mum; each cell is a self-bound gausson (width \xi_\mathrm{GP}=\xi/\sqrt2) — held by dc1’s logarithmic equation of state. Middle (vertical): these sheets stack along a perpendicular axis; like-handed +\omega sheets (blue) repeat at the Glaberson–Johnson–Ostermeier period d_\mathrm{GJO}\approx16\,\mum, with a counter-rotating -\omega intermediate layer (red) at the half-period, so a boundary falls every \approx8\,\mum — the substrate’s own octave — while the vortex lines thread the stack perpendicular to the sheets. Right (the seam): each +\omega vortex and its -\omega partner breathe in anti-phase — one expanding as the other contracts, trading energy across their shared seam without loss, at the dc1 hum frequency \omega_1=m_1c^2/\hbar\approx3\times10^{12}\,rad/s; because the pair is anti-phase the hum cancels above one cell. Held standing, the pair is a Cooper pair; set propagating, it is the photon (Figure 1).

4.2 The quadrupole honeycomb

The cancellation is not perfect, and the residual nests in the spaces in between, further hiding the energy of these vortex lines. Each cell is an anti-phase pair wrapped in one \xi-wide envelope; the envelopes tile on the triangular array. The sliver the breath can’t cancel leaks through the lattice boundary, and this surviving quadrupole holds the lines in the gaps between them. Those gaps are a single connected honeycomb of hollows, the lattice dual to the triangular array, one hollow between each trio of cells. The residual threads that honeycomb and self-screens to \sim1\% within one cell (the e^{-4.44} falloff of a source whose dipole has been arranged away), so a probe standing even one \xi away sees nothing. In Fourier language this is S(\mathbf q\to0)\to0: killing the dipole removes the long-range, small-\mathbf q weight and leaves only the quadrupole at the cell scale, the ring at q\approx2\pi/\xi. It is why a medium at the dark-matter density presents as a perfectly fluid soft boundary — the substrate hides in the seams of its own fabric.

Three panels. Left, the tiling: a triangular array of blue circular cell envelopes, each holding an anti-phase pair of a blue plus-omega core over a red minus-omega core, with the cores varying in size while the envelopes stay the same, spacing marked xi approximately 100 micrometres; the gaps between each trio of envelopes are filled by a gold honeycomb net of hollows with a small four-lobed quadrupole glyph at each node. Middle, a zoom into one hollow: three cell envelopes meet around a curved-triangle gold gap holding a four-lobed quadrupole, labelled that anti-phase cancels the dipole and the quadrupole survives, a perfectly fluid soft boundary. Right, a log plot of leak strength versus distance in lattice constants: a gold anti-phase quadrupole curve falling to about one percent within one cell, a grey dashed aligned-dipole curve reaching much farther, and a caption that removing the dipole gives S of q going to zero as q goes to zero with the surviving weight at a single ring at q about 2 pi over xi, beside a small dark ring glyph.
Figure 9: The soft boundary — the honeycomb between the cells. Left: the cells tile on a triangular array, each an anti-phase pair (+\omega over -\omega) wrapped in one \xi\approx100\,\mum envelope; the cores vary in size from cell to cell while the envelopes stay identical — varying-energy modons sharing one wake (the boat-in-the-harbor scale separation of §2.1). The gaps between each trio of envelopes form a honeycomb of hollows, dual to the array. Middle: a zoom into one hollow — the anti-phase breath cancels the dipole (the long-range, small-\mathbf q leak); what can’t cancel is the quadrupole, bottled in the gap between the trio, a perfectly fluid soft boundary that is nearly, but never quite, closed. Right: why it hides — the anti-phase quadrupole leak self-screens to \sim1\% within one lattice constant (e^{-4.44}), while an aligned dipole would reach macroscopically; removing the dipole is S(\mathbf q\to0)\to0, the surviving weight pushed to the single diffuse ring at q\approx2\pi/\xi — the stealth vacuum, read from its seams.

4.3 The ladder: a √2 lock and a φ anti-lock

At boundaries on every scale, two patches of the lattice meet across a seam — two of the chirality-coherent sheets of §4.1, or a structure’s ordered interior meeting the resting substrate outside it. Each side is held by its own chemistry. In ordinary fluid mechanics the seam between two flows is a featureless shear layer that smears together by diffusion, with no length or shape to organize it. But the substrate’s seam is a counter-rotating boundary layer between stacked lattice sheets, with a definite spacing — a layered attraction/repulsion that chemistry uses to shape boundaries. How those sheets line up defines a ratio — something the next layer can either match or avoid. The ratio is mediated and adjusted by chemistry.

When the boundary layers of two sheets meet, they choose lock or anti-lock as the surrounding chemistry dictates. The next layers above and below can follow, tiling like a crystal and locking or anti-locking in cascades, the family of patterns that chemistry already explains well. The lattice is otherwise scale-free, emitting light at one speed with no other break. All it has is this mediated ratio that repeats: a boundary at one scale seeds a like boundary at the ratio times that scale, and that one seeds the next, up and down. The pattern propagates the way a crystal grows — start from a single rung and it tiles outward in both directions, start from a handful and they tile together. What is copied from rung to rung is not a length but a ratio. At the limit it could even create discrete scale invariance, the mechanism behind the Efimov effect’s universal 22.7 tower (Braaten and Hammer 2006) and Sornette’s log-periodic precursors to rupture.19 The result looks fractal, but it is one repeating energy pattern at a fixed ratio, influencing chemistry to build the next ratio out of it.

So when a structure locks on a rung, it can sit in register on a tooth: an integer ratio or the bare rung. This is how parts bind, nest and hand energy back and forth. Examples include the cochlea’s octave layout, the entorhinal grid modules’ measured \sim1.4 spacing, and vesicle-coat and microtubule nesting.

On the other hand, when it anti-locks, it sits as far from every tooth as a ratio can, refusing to resonate. This is how parts stay independent and avoid overlapping: when chemistry chooses the gap, they repel through the dissonance.

The teeth (\sqrt2, the octave) allow binding and the gap (\varphi) keeps the boundary separate. These come from the substrate’s pairing factor: \xi^2 = 2\,\xi_\mathrm{GP}^2.

Chemistry chooses the moves and how they roll up across the boundary. Do the locks tile evenly, or are they interleaved with anti-locks at intervals? The pure anti-lock pattern spirals up the scale ladder.

The substrate’s tower is formed from half-octaves of \sqrt2; the octave takes two rungs — the 8\to16\,\mum vertical span between sheets. The gap is the one number that refuses every tooth at once, the most-irrational \varphi=1.618 — the ladder’s shadow.20 The same gap takes a different arithmetic in each space it lives in: \varphi on a circle (phyllotaxis’s golden angle), disordered hyperuniform “blue noise” on a plane (the retinal cone mosaic), and mutually-prime periods in time (the 13- and 17-year periodical cicada) — one principle, maximal incommensurability, three geometries.

A two-panel schematic of a self-similar tower. Left, a side view: a vertical dashed axis with the same motif — a teal disc and an indigo disc joined by a short dashed seam with arrowheads, each disc wrapped by a counter-rotating arc — drawn four times at sizes a factor of root two apart, labelled from the bottom up as the vortex core xi-GP, xi equals root two xi-GP, times two one octave, and times two root two, with a dotted continuation to organelle and cell. Brackets at the far left mark one octave as times two and one half-octave as times root two, with an example span of 8 to 16 micrometres. Right, a top-down view: concentric rings spaced by root two seen down the tower's axis, with the same teal-and-indigo paired motif sitting on every ring at a size matching that ring. Caption beneath: zoom the whole picture by root two and it lands on itself.
Figure 10: The self-similar tower: one motif, stamped at every \times\sqrt2 rung. The tower is a single motif — a teal vortex breathing against its indigo counter-rotating partner, energy traded across the dashed seam — repeated at a geometric ladder of sizes a factor of \sqrt2 apart. Left, side view: the staircase climbs from the vortex core \xi_\mathrm{GP} by \times\sqrt2 to \xi, by another to the octave \times2, and on up toward organelle and cell — two half-octave rungs to each octave, the 8\to16\,\mum span among them. Right, top-down view: the same tower seen down its axis is a nest of concentric rings spaced \times\sqrt2, the identical paired breath sitting on every ring. Rescale the whole picture by \sqrt2 and it lands on itself: that invariance is the tower. The rungs are discrete rather than a smooth power law because the critical substrate’s renormalisation-group flow closes into a limit cycle, one lap of which is \times\sqrt2.

The boundary can thus organize, not merely diffuse. When two lattices meet, their energy and chemistry set up a lock/anti-lock potential at every scale, with a variety of subtle effects.

It helps flows persist far past what viscosity should allow — a coherent ocean current, a vortex cell, a cell’s membrane network sit on rungs and are held there, trading energy along the lossless seam instead of dissipating it. The ladder is a hidden potential surface lying under structure at every scale.

The same boundary dynamics also set angles. When the tiling wraps a plane or a center, the lock/anti-lock choice becomes an angle. The cleanest lock is three-fold, 120^\circ — the lowest-frustration closed cycle, the angle of the hexagonal sheet and of the cell’s three-way membrane junctions. The cleanest anti-lock is the golden angle, 137.5^\circ = 360^\circ/\varphi^2, whose rational near-misses (1/3\to120^\circ, 3/8\to135^\circ, …) are the Fibonacci spiral arms a growing front slips through on its way to never locking at all. One ladder, two ends: the angle that nests and the angle that refuses.

A single wide panel with three sub-figures. Left, a golden Fibonacci spiral of dots, teal at the center grading to amber at the rim, labelled the golden spiral and anti-lock, phi, 137.5 degrees. Centre, a three-fold junction: three green arms meeting at 120 degrees inside a faint dashed hexagon, each arm carrying a small counter-rotating vortex glyph, with the 120 degree angle marked by an arc, labelled the 120 degree three-fold and lock, hexagons, ER junctions. Right, a fan of grey divergence rays opening clockwise from a reference ray marked previous, with the golden 137.5-degree ray drawn heavy in amber, labelled and every angle between, from 120 to 180 degrees the golden one never repeats. Beneath all three, a single line: mix the two rung by rung, arrow, even nesting, intermediate spacings, golden spirals, hexagonal sheets, branching networks, blue-noise mosaics, prime cycles.
Figure 11: The repertoire of the lock/anti-lock choice. From the ladder of ratios, chemistry uses the rung-by-rung choice to set the shape. Left: pure anti-lock — successive elements placed at the golden angle 137.5^\circ=360^\circ/\varphi^2, which never repeats and so packs without ever aligning, the Fibonacci spiral of phyllotaxis and the blue-noise mosaic on a plane. Center: clean lock — the three-fold 120^\circ junction, the lowest-frustration closed cycle, the angle of the hexagonal sheet and of the cell’s three-way membrane junctions. Right: the two are the ends of one continuum, not a binary — the fan of divergence angles from 120^\circ to 180^\circ, every one of them available, with the golden 137.5^\circ marked as the single member that never repeats. Mixed rung by rung, the same two moves compose even nesting, intermediate spacings, hexagonal sheets, branching networks and mutually-prime cycles in time.

Locked boundaries are stiff and only crossed by a punch-through. It’s more expensive to cross the locked boundary, so flow will be direct, not at an angle — like subducting slabs that punch through the mantle’s 660-km discontinuity or stall against it, or flux lines threading a type-II superconductor one quantum at a time.

One clear example from cellular dynamics is how the substrate energy activates the endoplasmic reticulum: a membrane network that grows and re-knits itself like a crystal in real time, its three-way junctions at 120^\circ, its tubule lengths clustering on the ladder’s rungs rather than spreading smoothly, its junctions migrating continuously to find the substrate-energy minimum. It is the ladder made dynamical — lock-pole geometry running a live search for its rungs across the inside of a cell.

4.4 The same loop at every scale

With mass as leaking rotational energy, a vortex, a vortex knot, an orbiting pair, a whole complex of orbiting structures — anything carrying more rotational energy than its counter-rotating layers can cancel — has a net spin it cannot hide, and that uncancelled remainder is felt as mass. A fermion, with odd parity, is unbalanced: a persistent vortex with one rotational layer too many. Wrapping boundary layers must contain that energy, and the way they do forms a common feedback topology.

Unbalanced systems at all scales show the same pattern: a co-rotating equatorial disk, two opposing polar jets along the spin axis, and an enclosing counter-rotating boundary layer. The disk shape has the lowest energy — least moment of inertia per unit boundary. At the same time, the lowest-energy exit for the excess momentum is straight up or down the axis — the jet. The cheapest return is a counter-rotating sheath wrapping the disk’s rim, conveying flow back inward while matching the substrate’s own boundary. Disk, jets, counterflow, and a thin remainder radiated outward as waves. It is boundary-minimization made dynamical: the reason a closed torus is the cheapest envelope for a standing loop becomes, for a system that must shed, this disk-jet-counterflow loop.

Twelve small line-drawings on a white background, the same canonical loop redrawn in its own form at each scale and arranged in two rows of six, from the electron at 10 to the minus 12 metres up to the cosmic bubble at 10 to the 26 metres: electron, hydrogen atom, codon, DNA, mitochondrion, neuron, cortical column, human being, Earth, Sun, galaxy, and bubble, color-banded blue for the quantum scales, teal for the biological scales, and amber for the cosmic scales. A key at lower left reads the loop three ways — a blue polar jet labelled the signal it radiates outward, a red counter-rotating layer labelled the hidden balance that wraps it, and a slate spinning mass labelled the feedback loop at the core — beside a blue-and-red yin-yang emblem and the closing line 'the shape of the unbalanced: what a leaking spinner does with the energy it cannot keep.'
Figure 12: The same loop at every scale. A net-spinning mass dropped into the stiff, low-dissipation substrate organizes the lattice around it identically at every scale: a co-rotating equatorial disk (slate), an enclosing counter-rotating boundary layer (red), and two polar jets along the spin axis (blue), with a thin remainder radiated outward as waves. The geometry is forced, not fitted — least moment of inertia per unit boundary for the disk, the cheapest exit straight up the axis for the jet, the cheapest return a counter-rotating sheath. Because the substrate is near-lossless the loop persists, so the same three-part machine recurs across ~25 orders of magnitude: the electron, a hydrogen atom, a codon, DNA, a mitochondrion, a neuron, a cortical column, a human being, the Earth (geodynamo with auroral funnels), the Sun (convective disk and polar wind), a galaxy (galactic disk with AGN jets), and the cosmic bubble. This is the shape of the unbalanced — what a leaking spinner does with the rotational energy it cannot keep — and the complement of the balanced modon (Figure 13).

The excess energy in the substrate balances two ways at once. The jets offload angular momentum up the axis; the disk weaves it outward into the lattice’s repeating fabric, where it is stored, twisted, and fed back. Any disturbance to a lattice sheet pushes energy into the axis or pulls it away from it — modulating the jets. The jets’ pull controls the rotational speed of the vortex. A leaking spinning mass creates a feedback loop: in through the disk, out through the jets, traded along the lossless counter-rotating seam, a fraction radiated — two coupled loops in one lattice containment field.

The loop has no viscosity and no measurable decay because it cycles the energy along topologically protected vortex lines, locked into frictionless motion. This three-part machine recurs across all scales: the substrate vortex cell at \xi\approx97\,\mum, an accretion disk with its relativistic jets, the geodynamo with its auroral funnels for jets, a galactic disk, the sheets of the cosmic web — and, turned inward, the nested loops that run a living cell.

Because the substrate prefers locally oriented planar order, its vortices lying in chirality-coherent sheets (§4.1), the disk plane is the substrate’s plane: orbital lobes, planetary rings, the ecliptic, accretion and galactic disks, and cosmic-web filaments are one geometry expressed in whatever local material is doing the rotating.

Frame-dragging observations are explained as the loop’s azimuthal component: a rotating mass entrains the azimuthal part of the gravitational inflow (§3.1) through the same mutual friction that couples the components of rotating superfluid helium, reproducing the Lense–Thirring rate for a slow rotator — GR predicts 39.2 mas/yr for Gravity Probe B, which measured 37.2\pm7.2 mas/yr.21

Now let’s move to the massless topology, where the spinning energy is not leaking but balanced and self-sustaining.

4.5 The modon topology

When there is a balanced pair in the substrate — or a larger group organized as nested balanced loops, like a cell with its mitochondria — it behaves in many ways like the modon. The energy balances out, effectively hiding the structure from observation. In motion, the substrate energy appears massless like the photon, its momenta offsetting. At rest, the pair becomes counter-spinning layers of coherent energy circulation, obeying the same Bessel boundary-matching as the modon but held in place.

These are the coherence layers, the standing fabric the leaking loops weave into. The benzene ring forms a toroidal vortex to share electrons — a modon in the substrate. Aromatic stacks of vortices are composite modons of coherent energy.

A white-background figure with two rows of six small line-drawings, each a counter-rotating dipole with a blue positive-vorticity core and a red negative-vorticity core. The top row, labelled MOBILE, shows the self-advecting modons: the photon, a Gulf Stream ring, the MJO twin cyclones, a Neptune storm pair, the counter-rotating galaxy NGC 4550, and a binary black hole. The bottom row, labelled STATIONARY and NESTED, shows the held-in-place modons: a Cooper pair with anti-phase breath arrows, a benzene ring with two counter-spinning faces, the DNA double helix with one strand blue and one red, a cell nested with a red mitochondrion, the brain's blue and red hemispheres joined by the corpus callosum, and Jason and Tuzo — Earth's two antipodal mantle superplumes drawn as a cutaway of the planet with two opposite-spinning blue and red feet standing on the core–mantle boundary. A key at lower left draws one dipole — equal and opposite cores, one Bessel match, net spin zero, net mass zero, radiates nothing — beside a blue-and-red yin-yang emblem and the closing line 'the shape of the balanced: energy with nothing to shed.'
Figure 13: The modon topology, at every scale. Replace the one net spinner with two counter-rotating cores — +\omega (blue) over -\omega (red), sharing one envelope across a clean counter-rotating seam — and the bookkeeping inverts: equal and opposite spins, net angular momentum and net mass transport both zero (\int\psi\,\mathrm dA=0), nothing uncancelled, nothing to shed. The pair is held by the same Bessel structure at every scale (J_1 inside the envelope, with a K_1 halo outside when it travels), in two forms. Top — mobile, when the pair advects itself: the photon at the substrate scale, and the identical solution in seawater, planetary air, and spacetime scaled up by ~10^{13} — a Gulf Stream ring, the Madden–Julian oscillation’s twin cyclones, Neptune’s dark-spot pair, NGC 4550’s counter-rotating disks, a binary black hole. Bottom — stationary / nested, when the dipole is held in place: a Cooper pair breathing anti-phase, a benzene ring’s two counter-spinning faces, DNA’s two counter-wound strands, a eukaryotic cell nested with its mitochondrion, the brain’s two hemispheres across the corpus callosum, and — one rung up in scale, to the whole planet — Jason and Tuzo, Earth’s two antipodal mantle superplumes (the LLSVPs), opposite-spinning feet standing on the core–mantle boundary, the slowest stationary modon the planet supports, held antipodal across Wilson cycles for \gtrsim 300 Myr. This is the shape of the balanced: the standing fabric the unbalanced loops (Figure 12) weave into.

Coherence scales by nesting topological modon layers. Each passes its boundary to the next by boundary matching, and the number and smoothness of those matches set the structure’s coherence — its stability and longevity.

Each modon holds a coherent, stable pattern — the modon’s coin. A photon is the simplest modon — a single frequency stored as persistent energy. But the topology also allows a structured or composite modon — an aromatic stack, or a nested cell — to hold and transmit a long vector of information, using the layered lock/anti-lock pattern along a stack of lattice sheets. It enables structures to perform a deep coherence match from these pure energy patterns. This offers a new lens on cellular dynamics: organized by substrate energy, not by diffusion alone.

The lattice cell size is also about the size of a living cell. The conjecture is that a cell is a many-layered nested topological modon, built from the plasma membrane, cortex, nuclear double-wrap and mitochondrial double-membrane, and that a neuron’s myelin adds more wrapping layers for more coherence.

Three cards in a left-to-right spectrum. Left card, ONE RUNG, the photon: one counter-rotating blue-over-red dipole with a c arrow, and beneath it a single blue bar labelled psi equals a-one, captioned one frequency, one number, E equals h nu. Centre card, A FEW RUNGS, an aromatic stack or chord: a vertical stack of three counter-rotating rings, and beneath it three teal bars, psi equals a-one a-two a-three. Right card, MANY RUNGS, a brain modon: a composite of four counter-rotating cores in one envelope with a faint motion arrow, and beneath it eight teal bars of varying length, psi equals a-one a-two up to a-N, captioned the full comb aboard, a long vector sent whole. A spectrum arrow beneath reads one object, more rungs to the right, basis log-spaced by root-two not harmonic. A closing gold band shows a blue-red modon yin-yang labelled THE MODON COIN beside the equality energy equals pattern, and the line: the modon's energy is its long vector, one packet carries power and meaning, passed boundary to boundary by the coherence-match, at c, without loss.
Figure 14: The modon carries the long vector. A modon’s complexity is the dimension of the long vector it carries. Left: the photon is the one-rung special case — a single counter-rotating dipole carrying one number, its energy E=h\nu but no pattern. Center: a modon with structure — an aromatic stack, a chord — binds multiple rungs into one travelling object, a short long vector. Right: a composite modon — like a brain — has a full comb, a long vector \psi=(a_1,\dots,a_N) log-spaced by \sqrt2 (the ladder of §4.3). Bottom: the modon’s energy is a long vector, so a single substrate packet transmits energy and information together, passed boundary to boundary by the coherence match \langle a|b\rangle at c without loss.

These ideas about cells are conjecture, not results of this paper. They are developed at https://lightfluid.org.

5 What the framework predicts

The framework’s backbone derives the lattice cell size \xi and ties it to the condensation number \nu=m_\mathrm{eff}/m_1, which gives the core predictions with no additional parameters.

The Tier 1a predictions use the Weinberg angle, the electron’s rest mass, Planck’s constant \hbar, the speed of light c, the gravitational constant G, the dark matter density \rho_\mathrm{DM}, one geometric backbone (f=4\pi/(K\sqrt2)), and one dynamical thread (\alpha_{mf} = \tan^2\theta_W = \sin^2\theta_W/(1-\sin^2\theta_W) = 0.3008, the mutual-friction coupling of §1.1).

The Weinberg angle used, 0.23122\pm0.00006 (Takahashi et al. 2026), is the \overline{\mathrm{MS}} value at M_Z, where it is defined from renormalized couplings, matching \alpha_{mf}=g'^2/g^2 as a ratio of gauge couplings (§1.1). The electron’s mass at M_Z absorbs the top-quark and Higgs loop corrections, so the rest mass is preferred when combining it with \alpha_{mf} to compute the effective quantum’s mass-energy.

5.1 Tier 1a — backbone predictions

Quantity Expression Predicted Observed Discrepancy
Higgs VEV v \sqrt{8\pi\,m_\mathrm{eff}^2 c^4\nu}22 246.1\pm0.3 GeV 246.22 GeV -0.04\%\pm0.13\% (inputs only)
Koide relation, read geometrically Q=\tfrac13+A^2/6, A=\sqrt223 2/3 0.666664\pm0.000005 0.4\sigma
MOND scale a_0 c\sqrt{G\rho_\mathrm{DM}} 1.16\times10^{-10} m/s² (1.20\pm0.02\pm0.24)\times10^{-10} m/s² \sim3\%
Landau critical velocity v_L (same \nu) c\,(4\pi/\nu)^{1/3} \approx740 km/s Ulysses fast wind 751.5 km/s 1.5\%
Min localized-modon energy hc/\xi = 2\pi m_1c^224 12.8 meV — crossover at \lambda=\xi
Cosmic coincidence (from a_0) \sqrt{3\Omega_\mathrm{DM}/8\pi} 0.178 0.183 \sim3\%

The VEV and a_0 are the two strongest independent tests. The Landau velocity reuses the VEV’s \nu to find a weaker correspondence to the fast solar wind. The cosmic coincidence is a_0 rewritten through the Friedmann equation.

In the Koide row, the three generations are seen as the three branches of one three-fold junction (the \mathbb{Z}_3 behind color and the \pm\tfrac23 quark charges). The phase then cancels, and Koide’s ratio depends only on the amplitude of the mass pattern, Q=\tfrac13+A^2/6. Because of this, A=\sqrt2 implies Q=\tfrac23. This is a geometric reading of Koide’s 1981 relation where \sqrt2 comes from the pairing factor \xi^2=2\xi_\mathrm{GP}^2, the same factor that appears in the cell occupancy f. It also accounts for why there are three generations, since of N equal phases only N=3 puts the balanced value on \tfrac23.

5.2 Tier 1b — predictions with an open factor

Quantity Expression Predicted Observed Discrepancy Open Factor
Muon/electron ratio m_\mu/m_e three-phase clock, \delta=2/9 rad25 206.77 206.768 +0.001\% 2/9
B-DNA bp/turn 2\pi r f/h26 10.47 10.5\pm0.1 0.3\% \tan\alpha_\text{pitch}=f
Microtubule wall 3/(2f)27 2.647 2.594 2.1\% \tan\alpha_\text{3start}=f/\pi

Both m_\mu/m_e and m_\tau/m_e use the same \delta=2/9 rad matched from the data. The 2/9 is read geometrically as pairing-two over three-fold-squared, but that is not a formal derivation.

Both the B-DNA and microtubule rows tie a helix tangent to f, motivated by the chiral/boundary projection but not derived. DNA splits its motion into \sin\alpha across substrate sheets and \cos\alpha within them, so \tan\alpha_\text{pitch}=f. The microtubule uses the same tangent definition, axial rise per cycle over perimeter per cycle, but its closed wall locks to f/\pi. That extra \pi is part of the open factor.

For consistency, the same scattering phase used for \alpha_{mf} should give the fine structure constant as a boundary relation among couplings, \alpha=\sin^2\!\delta_0\sin^2\theta_W/\pi28, as seen through the energy of a probe. Weaker probes are screened and so measure a smaller value. For \theta_W measured at the Z pole, with a probe \mu=M_Z=91 GeV, the relation gives \alpha=1/135.1, too weak by 4.5% against the observed value of 1/128. The coupling matches with \mu^*=168–182 GeV, which brackets both the top-quark mass and a candidate based on the geometry: v/\sqrt2=174 GeV, the Higgs VEV divided by the recurring \sqrt2. Walking that value back to zero momentum using the observed screening gives 1/\alpha(0)=136.7–137.5, consistent with the measured 137.036.

5.3 The hidden lattice

One reason a hundred-micron superfluid filling space can be so stealthy is the way it transmits radiation. Photons move through the free substrate along the least-energy path — Lorentz invariance — a property that emerges because every vortex, every excitation of the substrate shares the same cone c=\hbar/m_1\xi and the same timing.

The photon stays coherent, meaning the lattice must have a structure factor that prevents spatial scattering: S(\mathbf q) is the overlap of its density with the probe wave. The function might predict diffuse or even sharp scattering, but the photon keeps its shape, revealing a property of the lattice. The texture must be a domain glass of randomly oriented triangular Abrikosov crystallites, the same assumption made with the bridge equation’s \eta=1. It’s a disordered hyperuniform solid: S(\mathbf q\to0)\to0. There’s no comb, no dispersion. Instead the model predicts a single diffuse ring at the cell scale |\mathbf q|\simeq2\pi/\xi, a powder ring with an edge that coincides with the modon-localization crossover E_\mathrm{min}=2\pi m_1c^2. This is essentially the lattice size reached a third way.

Feeding that S(\mathbf q) into a first-order (Born) scattering integral fixes where the vacuum could first scatter far-IR light. A photon can transfer momentum only up to 2k=4\pi/\lambda, so every channel for \lambda>2\xi\approx200\,\mum (\nu<1.5 THz) lies inside the stealth window and the vacuum is strictly transparent there, as stealthy hyperuniform media are (Leseur et al. 2016). As the wavelength shortens below 2\xi, a faint scattering edge turns on and rises to the ring at \lambda=\xi (\sim3 THz). The observed far-IR/submm sky is transparent: Herschel and ALMA resolve dusty star-forming galaxies out to z\sim4–6 (Pilbratt et al. 2010; Casey et al. 2014). That bounds the per-cell cross-section below \sim10^{-31}\xi^2, while a random gas of cells at the same density would be opaque by \sim30 orders of magnitude. The sky’s far-IR transparency therefore requires the hyperuniform suppression of the domain glass. The edge itself, pinned to the cell scale and nothing below it, is the prediction still to be tested.

The same stealth answers an obvious question about gravity. A \sim100\,\mum superfluid lattice filling space sits squarely inside the sub-mm window that Eöt-Wash torsion balances have swept for deviations from Newtonian gravity (roughly 40\,\mum to 11 mm) — yet they see none (Lee et al. 2020). The cell scale does not present itself to laboratory masses as an extra gravitational structure: ordinary matter couples to the lattice only through the same boundary leak that makes gravity weak — a crossing probability of \sim4\times10^{-31} per breath of the cell (Equation 16) — so the cell adds no separate oscillatory sub-mm force for a torsion balance to find. The same universal weakness satisfies the equivalence-principle tests — one gravitational channel, common to all matter, produces no composition-dependent force — consistent with MICROSCOPE’s \eta\lesssim10^{-15} (Touboul et al. 2022). The one macroscopic scale above the cell — the coherence ceiling \ell_L\approx3.9 cm — is a smoothness null for the same reason, not a signal.

5.4 Tier 2 — physics re-derived

This section re-derives accurately measured physics with a substrate interpretation. It shows the model’s consistency and in some cases adds a physical reading to a known relation, like the W/Z mass ratio, frame-dragging, and (g-2)/2=\alpha/2\pi.

  • The Bohm quantum potential (§1)
  • The hydrogen Rydberg spectrum from pilot-wave standing-wave matching29
  • Entanglement and Bell’s E(\theta)=-\cos\theta, CHSH 2\sqrt2 — from a twist wave on the vortex channel joining the split pair30
  • The Born rule for a spin projection, \cos^2(\theta/2) — the single-quantum Malus law — from the reactive gear-reduction of the counter-rotating boundary layer as the analyzer axis tilts31
  • The Aharonov–Bohm phase e\Phi/\hbar — real and topological, because the vector potential is the substrate’s co-rotating flow the charge swims in, not a bookkeeping device, so a charge circling shielded flux picks up its enclosed circulation even where the field vanishes32
  • The Sagnac phase \Delta\varphi=\tfrac{2m}{\hbar}\vec\Omega\cdot\vec A — a rotation imprints the enclosed substrate circulation and the matter-wave fringe count is the number of circulation quanta h/m threading the loop, the frame-dragging term read for a laboratory apparatus rather than a gravitating mass33
  • The Casimir force -\pi^2\hbar c/240\,d^4 comes from a boundary-shaped modon pressure, a rearrangement of the vacuum’s field, not its energy, with no zero-point sea, per Schwinger/Jaffe, exact where measured (10 nm–μm)
  • Vacuum magnetic birefringence \Delta n\propto B^2 shows that field-aligned dc1 flow is the optic axis and the modon’s splitting is the phase delay with the field’s own flow as the lattice rows. It reproduces the Euler–Heisenberg form and scale (B/B_c)^2; magnetar RX J1856 (Mignani 2016), PVLAS bound.
  • The electron g=2 and the anomalous moment from two layers with different moments of inertia34
  • (g-2)/2=\alpha/2\pi; quark charges +\tfrac23,-\tfrac13 from the junction solid angle, matching the three helon strands of the preon-braid model (Bilson-Thompson 2005), where an up quark has two twisted strands to a down quark’s one. \sim99\% of the proton mass comes from boundary-layer energy (matching lattice QCD)35
  • The quantum Hall effects — fractional charge e/3 as a Laughlin vortex’s winding fraction (the same charge-as-winding ontology as the quark \tfrac13, here measured directly by shot noise), integer quantization \sigma_{xy}=\nu e^2/h from topological protection of quantized circulation (exact to <10^{-9}, the resistance standard), and the odd/even-denominator rule as the boundary-parity fermion/boson split — even denominators the paired breath at \nu=\tfrac5236
  • The W/Z mass ratio M_W/M_Z=\cos\theta_W — the standard gauge relation, but as the equatorial spin velocity \beta_\mathrm{eq} of the counter-rotating boundary shell: \cos\theta_W=0.877 against the measured 80.4/91.2=0.882 (0.5\%)37
  • The CKM-small / PMNS-large mixing asymmetry — quark mixing is small and lepton mixing is large because a mixing matrix is the relative rotation of two \mathbb{Z}_3 generation clocks: up- and down-type quarks are color-loaded to the same side of the generation lock (co-aligned \to small CKM), while charged leptons and neutrinos straddle it (anti-aligned \to large PMNS); the atmospheric angle sits at the pairing-\sqrt2 tilt, \theta_{23}\approx45° (\cos45°=1/\sqrt2), and a large CP-violating Dirac phase is generic in both sectors as the imaginary part of the complex clock \omega=e^{2\pi i/3} — the sign of the asymmetry and near-maximal \theta_{23} reproduced, the individual angles and the CP value still owed38
  • Frame-dragging (Lense–Thirring) at the Gravity Probe B rate — the entrained azimuthal component of the same gravitational leak-inflow, v_\phi=2GJ\sin\theta/(c^2r^2) (§4.4): geodetic 6604 vs GP-B 6601.8\pm18.3, frame-drag 40.7 (simple orbit-averaged GJ/2c^2a^3; the full GR mission prediction is 39.2) vs 37.2\pm7.2 mas/yr
  • Rest energy E=mc^2 (§1.2) — the flywheel energy of the internal vortex, \tfrac12m_\mathrm{eff}v_\mathrm{rot}^2 with the rim locked at v_\mathrm{rot}^2=2\alpha_{mf}c^2, with an aperture \alpha_{mf} that decides how much leaks out and is weighable; the \tfrac12 and the 2 cancel and m_ec^2 is what remains, with c derived (Equation 5) and \alpha_{mf}=\tan^2\theta_W fixed, so nothing is left free to choose
  • Special relativity — time dilation \gamma, length contraction, and E=\gamma mc^2 — from a moving light-clock: the electron’s internal Compton breath runs on c-bounded substrate signals, so a moving clock slows by the transverse-light-clock factor \gamma=1/\sqrt{1-v^2/c^2} and the de Broglie wave is that same clock seen sideways, unifying the Michelson–Morley null with the pilot wave (exact, zero-parameter)39
  • The baryonic Tully–Fisher law M_b\propto v^440
  • Kleiber’s metabolic law B\propto M^{3/4} from the same boundary-leakage theorem read on a branching vascular network41
  • The inflationary spectral index n_s\approx0.967 from a sound-speed phase transition42
  • The nuclear binding curve showing iron as the peak.43
  • Mass defect \leftrightarrow EMC effect. Mass is a boundary leak, not a bulk count, so binding merges two counter-rotating boundaries into one internal seam — dropping the leaked ledger (the mass defect) and reshaping the reactive ledger (in-medium quark structure, quenched moments/g_A). The sign and universality of the defect are structural, following from a total of -E_\text{bind}/c^2 (which both frameworks match), with EMC strength \propto binding per nucleon.
  • Superconductivity’s Cooper pair as the substrate’s anti-phase breath made macroscopic — two same-chirality electrons \pi out of phase in their Compton breathing, one contracted while the other expands, bound by the shared counter-rotating vortex between them (the BCS singlet ↑↓ read as opposite phase, not spin), the charged twin of superfluid ^3He and the lab realization of the very +\omega/-\omega sheet pair the substrate is built from (§4.1)44
  • The Hawking temperature T_H=\hbar c^3/8\pi GM, exact — the sonic horizon of the same Painlevé–Gullstrand inflow (§3.1), with 8\pi=2\times4\pi_\mathrm{SC2} the gravitational normalization shared with the Higgs VEV45
  • The black-hole entropy area law S_\mathrm{BH}=A/4\ell_\mathrm{Pl}^2 — entropy as boundary-crossing leak counted on the one-way horizon seam, hence area not volume46
  • Neutron stars show the substrate’s vortex unpinning and mutual friction in a neutron superfluid: the HVBK/\alpha_{mf} boundary coupling matches the standard glitch model, with recovery set by mutual friction. A canonical 1.4\,M_\odot, 12-km star has a surface that sits at r/r_s\approx3, so the substrate inflow there is v_\text{ebb}=c\sqrt{r_s/R}\approx0.6\,c — the framework’s counter-rotating boundary skin at more than half the signal speed and still closing the seam that keeps matter matter. The dipole is the geodynamo, the canonical loop scaled a billionfold. The pulsar/magnetar split shows a consistent birth and evolution pattern where the rotational budget is partitioned between a coherent-jet channel (the beamed spin-down) and a stored-field channel (released in bursts and giant flares). The glitches are mutual friction (\alpha_{mf}) in that superfluid. The Crab and Vela wind nebulae are the canonical disk–jet–counterflow loop photographed, and the pulsar\tomagnetar sequence is a B/B_c ladder of natural birefringence laboratories climbing from \sim0.3\,B_c to \sim45\,B_c.
  • The baryon asymmetry \eta_B\approx\varepsilon_\mathrm{chirality}^9=5.8\times10^{-10} — all three Sakharov conditions are met, with the chiral-vacuum bias frozen across nine junction interfaces: three strands on each of three quarks. The ninth power of that bias gives an approximate match, 5\% low, which is about 8\sigma from the measured value47
  • The sonoluminescence flash width stays wavelength-independent into the far-UV / soft X-ray as the substrate predicts. It comes from the bubble’s boundary collapse, from a maximum R_\text{max}\approx 40–50\;\mum, about half the lattice cell size, to R_\text{min}\approx 0.5\;\mum at supersonic speeds. The same process that emits a spectral line when an atom’s boundary collapses here causes a stream of photon/modons spread evenly across the spectrum, without an orbital boundary to form the sharp line. Noble gases allow the substrate’s pure energy to be released as modons. Other atoms absorb energy and disrupt the flash, which offers an explanation for this poorly understood phenomenon.
  • The vacuum itself is a time-crystal formed from the anti-phase \omega_1 breathing. It is hidden in time, the way the lattice is hidden in space, beating a half-cycle out of step, summing to zero above one cell. Lab time crystals surface the vacuum’s breath through chemistry.

6 Reproducing the numbers

The code is at lightfluid.org and in the additional resources with this paper.

Most numerical predictions are reproduced by three Python scripts:

python3 scripts/substrate_atomic.py     # particle-physics backbone + lattice
python3 scripts/substrate_gravity.py    # gravity sector
python3 scripts/substrate_galactic.py   # dark energy and cosmology

7 Open problems

  • The 4\pi in SC2 and in the bridge equation, and the results that follow, rely on the substrate’s emergent Lorentz invariance being exact — supported by observation (Michelson–Morley, GW170817) but not derived. Once it is shown, the induced (Einstein–Hilbert) action supplies the 4\pi automatically.

  • The nonlinear completion of the Einstein equations requires a formal proof.

  • Derive the gravity number \mathcal F=Gm_1^2/\hbar c=(\ell_\mathrm{Pl}/\xi)^2\approx2.8\times10^{-62} (Equation 15), equivalently the per-breath crossing probability f_\mathrm{cross}\approx4.4\times10^{-31}. G is used here as measured. A derivation likely comes from the dissipative fraction \alpha_{mf} counted over a number of boundaries, together with a model of the channel created between two masses (§3.2). It would also derive a_0 in cell units and the e-fold count.

  • Close continuity for the gravitational ebb (§3.1). The two-fluid counterflow closes it exactly, provided gravitational waves ride the superfluid’s own phase mode; if instead they ride the total-density mode, they would feel no Shapiro delay, and GW170817’s arrival within 1.7 s of its gamma rays would rule the candidate out.

  • For dark matter, the quantitative step from the Hubble flow’s DC parity-bias to the coefficient of a_0=c\sqrt{G\rho_\mathrm{DM}} (§3.4) — equivalently, a phonon-mediated force profile reproducing the full radial-acceleration relation at galactic scales — is still open.

  • Finding the nuclear surface scale profile a_S. Only the iron-peak has been predicted so far.

  • Derive the proton’s seam count N_p\approx1836 (equivalently \alpha_{mf}^{(N)}=N_p\,\alpha_{mf}\approx552, §1.1) independently — from the topology of the nuclear vortex junction, the way \tan^2\theta_W fixes the per-seam leak \alpha_{mf}=0.3008. At present N_p is back-solved from the measured proton mass through the shared quantum m_\mathrm{eff}, so the proton-to-electron mass ratio is a re-parametrization rather than a prediction; an independent count would turn m_p/m_e\approx1836 into a genuine output.

  • Derive the condensation number \nu=m_\mathrm{eff}/m_1=8.35\times10^8 (Equation 12). It is now over-determined — cosmology \times the geometric cell occupancy, the measured Higgs VEV, and the fast solar wind all land on it, the first two to 0.08\% (inside the 0.25\% set by \rho_\mathrm{DM}, at the reference choices of §5.1) — but no bottom-up calculation produces it from first principles. It is the framework’s single sharpest open number: the lattice size, the Higgs VEV, and v_L all ride on it.

  • Derive the VEV relation v^2=8\pi\,m_\mathrm{eff}^2c^4\nu itself, from a chirality energy functional. The 8\pi=2\times4\pi and the \sqrt\nu are currently readings of a numeric result.

  • Derive the matching scale of the fine-structure relation (§5.2). The measured couplings satisfy \alpha=\sin^2\delta_0\sin^2\theta_W/\pi at \mu^*=168–182 GeV; a derivation placing the boundary relation at v/\sqrt2 would turn the 1/\alpha(0)=136.7–137.5 consistency check into a prediction good to about 0.3\%. Two related items: above M_Z the substrate’s own running law has \alpha_{mf} falling with energy while the Standard Model’s mixing angle rises, and the paper uses the M_Z value \alpha_{mf}=0.3008 throughout (m_\mathrm{eff}, the Higgs VEV), whereas the relation is satisfied with \alpha_{mf}=0.3065 at \mu^*. Which scale the leak fraction belongs to is not yet settled, and it matters: at \alpha_{mf}=0.3065 the Higgs VEV comes out at 239 GeV (-2.8\%) instead of 246.1 GeV. The electron mass carries the same open question: the baseline uses its pole mass, and the \overline{\mathrm{MS}} mass at M_Z would put the VEV at 229 GeV (-7.1\%).

  • Settle the modon constant’s boundary condition. K=j_{11}^2+1 is exact for a clamped cell, a node at r=\xi plus the medium’s gap; a free Larichev–Reznik modon with the same gap gives K=16.39, and a hexagonal clamped cell of the same area gives 15.99. The paper uses the clamped value for the standing lattice. Settling it means solving the gapped dipole mode on the actual lattice cell, which first needs two things the paper does not yet fix: whether \xi is the modon radius or the nearest-neighbour spacing (the separatrix sits at r=\xi, but neighbours are drawn at spacing \xi), and which neighbours are anti-phase, which decides whether each face of the cell is a node. Between the clamped and free values the predicted VEV shifts by 0.6\%, the DNA pitch and the microtubule wall by about 4.5\%, and \eta_B by a third.

  • Show that the substrate ladder is a genuine discrete scale invariance (a renormalization-group limit cycle).

  • Formal proof of the lattice spacing from the logarithmic equation of state’s maximum packing result.

Conclusion

I am not an expert in any of the domains this paper covers, but I have been reading science papers for a long time. I have spent the time on due diligence to make sure these results are consistent, to the best of my ability, and I think they deserve attention: this simple model reproduces the best existing theories and adds new relations among measured numbers, all built from one response function.

Volovik found key properties of the universe in a supercooled droplet of He-3, the one atom that stays fluid cold enough to tune into the anti-phase breathing of the substrate, and so essentially a mirror of the substrate itself. Bush, Zloshchastiev, Khoury, and many others supplied more pieces of the puzzle.

When the Madelung/Bohm and Simeonov fluid is treated as a superfluid, the quantum potential becomes the reaction force of the boundary layer. The photon makes sense as a modon — wrapping up these counter-spinning layers and using the substrate’s energy to launch itself. Once Planck’s constant, light’s minimum packet size, sets the modon’s minimum size, the lattice size is fixed and these predictions follow. The behavior of superfluids unites these results under one physics.

The stealthy superfluid turns quantum uncertainty into a measurement artifact. Entanglement is explained by the split pair’s topologically protected vortex lines, which still join the two particles across long distances. When you measure one end, a twist wave runs down that cleared, laminar channel faster than the speed of light, reproducing the Bell correlations without allowing signaling (§5.4). The unique properties of the superfluid — its stiff local behavior, its cold appearance from a distance, and its phase transitions — explain gravitational lensing and time dilation as a signal speed pinned to a flowing medium.

The wider investigation at https://lightfluid.org explores how vortices help explain covalent bonds, the energy behind cellular dynamics, and more.

AI Disclosure

Claude, from Anthropic, contributed major advances to this work and at the same time has been an obstacle I have had to manage carefully. I mostly use single-prompt sessions with a focused research goal and the relevant papers, to keep the science honest. I know enough not to trust Claude’s judgment on its own, which is why I spent the last year on due diligence, making sure I understand the results, checking the equations with independent Python scripts, and reviewing each one carefully many times using the best sources.

Acknowledgements

This paper was supported by the work of many scientists.

In particular, Volovik, Bush, Oza, Simeonov, Zloshchastiev, Khoury, Larichev and Reznik provided the foundational math and inspiration.

I’m also grateful to my amazing employer Posit PBC, supporters of open source software for data science. This is an unaffiliated side project, self-funded and released as open source material.

References

Avdeenkov, A. V., and K. G. Zloshchastiev. 2011. “Quantum Bose liquids with logarithmic nonlinearity: self-sustainability and emergence of spatial extent.” Journal of Physics B: Atomic, Molecular and Optical Physics 44: 195303. https://arxiv.org/abs/1108.0847.
Bancal, J.-D., S. Pironio, A. Acín, Y.-C. Liang, V. Scarani, and N. Gisin. 2012. “Quantum nonlocality based on finite-speed causal influences leads to superluminal signalling.” Nature Physics 8: 867–70. https://arxiv.org/abs/1110.3795.
Barceló, C., S. Liberati, and M. Visser. 2005. “Analogue gravity.” Living Reviews in Relativity 8: 12.
Bekarevich, I. L., and I. M. Khalatnikov. 1961. “Phenomenological derivation of the equations of vortex motion in He II.” Soviet Physics JETP 13: 643–46.
Bell, J. S. 1987. “How to teach special relativity.” In Speakable and Unspeakable in Quantum Mechanics. Cambridge University Press.
Berezhiani, L., and J. Khoury. 2015. “Theory of dark matter superfluidity.” Physical Review D 92: 103510.
Bialynicki-Birula, I., and J. Mycielski. 1976. “Nonlinear wave mechanics.” Annals of Physics 100: 62–93.
Bilson-Thompson, S. O. 2005. A topological model of composite preons. https://arxiv.org/abs/hep-ph/0503213.
Braaten, E., and H.-W. Hammer. 2006. “Universality in few-body systems with large scattering length.” Physics Reports 428: 259–390.
Bush, J. W. M., and A. U. Oza. 2020. “Hydrodynamic quantum analogs.” Annual Review of Fluid Mechanics 52: 235–80.
Casey, C. M., D. Narayanan, and A. Cooray. 2014. “Dusty star-forming galaxies at high redshift.” Physics Reports 541 (2): 45–161. https://doi.org/10.1016/j.physrep.2014.02.009.
Donnelly, R. J. 1991. Quantized Vortices in Helium II. Cambridge University Press.
Fetter, A. L. 2009. “Rotating trapped Bose–Einstein condensates.” Reviews of Modern Physics 81: 647.
Glaberson, W. I., W. W. Johnson, and R. M. Ostermeier. 1974. “Instability of a vortex array in He II.” Physical Review Letters 33: 1197.
Goles, E., O. Schulz, and M. Markus. 2001. “Prime number selection of cycles in a predator–prey model.” Complexity 6 (4): 33–38.
Hall, H. E., and W. F. Vinen. 1956. “The rotation of liquid helium II. II. The theory of mutual friction in uniformly rotating helium II.” Proceedings of the Royal Society A 238: 215–34.
Hamilton, A. J. S., and J. P. Lisle. 2008. “The river model of black holes.” American Journal of Physics 76: 519.
Khoury, J. 2015. “A dark matter superfluid.” https://arxiv.org/abs/1507.03013.
Larichev, V. D., and G. M. Reznik. 1976. “Two-dimensional solitary Rossby waves.” Doklady Akademii Nauk SSSR 231: 1077–79.
Lee, J. G., E. G. Adelberger, T. S. Cook, S. M. Fleischer, and B. R. Heckel. 2020. “New Test of the Gravitational 1/r^2 Law at Separations down to 52 \mum.” Physical Review Letters 124: 101101. https://doi.org/10.1103/PhysRevLett.124.101101.
Leseur, O., R. Pierrat, and R. Carminati. 2016. “High-density hyperuniform materials can be transparent.” Optica 3 (7): 763–67. https://doi.org/10.1364/OPTICA.3.000763.
Levitov, L. S. 1991. “Phyllotaxis of flux lattices in layered superconductors.” Physical Review Letters 66: 224.
McGaugh, S. S., F. Lelli, and J. M. Schombert. 2016. “Radial acceleration relation in rotationally supported galaxies.” Physical Review Letters 117: 201101.
Mehdi, Z., J. J. Hope, S. S. Szigeti, and A. S. Bradley. 2023. “Mutual friction and diffusion of two-dimensional quantum vortices.” Physical Review Research 5: 013184. https://arxiv.org/abs/2205.04065.
Pilbratt, G. L., J. R. Riedinger, T. Passvogel, et al. 2010. “Herschel Space Observatory — An ESA facility for far-infrared and submillimetre astronomy.” Astronomy & Astrophysics 518: L1. https://doi.org/10.1051/0004-6361/201014759.
Planck Collaboration, N. Aghanim, et al. 2020. “Planck 2018 results. VI. Cosmological parameters.” Astronomy & Astrophysics 641: A6. https://arxiv.org/abs/1807.06209.
Pletzer, Belinda, Hubert Kerschbaum, and Wolfgang Klimesch. 2010. “When frequencies never synchronize: the golden mean and the resting EEG.” Brain Research 1335: 91–102. https://doi.org/10.1016/j.brainres.2010.03.074.
Rosen, N. 1968. “Dynamics of a particle as a solution of a nonlinear field theory.” Journal of Mathematical Physics 9: 996.
Saffman, P. G. 1992. Vortex Dynamics. Cambridge University Press.
Sakharov, A. D. 1967. “Vacuum quantum fluctuations in curved space and the theory of gravitation.” Doklady Akademii Nauk SSSR 177: 70–71.
Simeonov, L. 2025. “Quantum potential from the material derivative of the osmotic velocity: a two-fluid Madelung framework.” https://arxiv.org/abs/2509.02868.
Sornette, D. 1998. “Discrete scale invariance and complex dimensions.” Physics Reports 297: 239–70.
Stensola, H., T. Stensola, T. Solstad, K. Frøland, M.-B. Moser, and E. I. Moser. 2012. “The entorhinal grid map is discretized.” Nature 492: 72–78.
Takahashi, F. et al. 2026. “Review of Particle Physics.” International Journal of Modern Physics A 41: 2630011. https://doi.org/10.1142/S0217751X26300115.
Touboul, P., G. Métris, M. Rodrigues, et al. 2022. “MICROSCOPE Mission: Final Results of the Test of the Equivalence Principle.” Physical Review Letters 129: 121102. https://doi.org/10.1103/PhysRevLett.129.121102.
Unruh, W. G. 1981. “Experimental black-hole evaporation?” Physical Review Letters 46: 1351.
Volovik, G. E. 2003. The Universe in a Helium Droplet. Oxford University Press.
Zloshchastiev, K. G. 2020. “Superfluid vacuum theory and deformed dispersion relations.” International Journal of Modern Physics A 35: 2040032.

Footnotes

  1. Simeonov’s derivation (Simeonov 2025) takes a single fluid of density \rho obeying the diffusion law \mathbf J=-D\nabla\rho; its velocity is the osmotic velocity \mathbf u=\mathbf J/\rho=-D\nabla\ln\rho. Form its material (convective) derivative — the osmotic acceleration \mathrm d\mathbf u/\mathrm dt=\partial_t\mathbf u+(\mathbf u\cdot\nabla)\mathbf u — and substitute the continuity equation \partial_t\rho=-\nabla\!\cdot\mathbf J=D\nabla^2\rho for the time part. The time and convective pieces collapse into a single gradient, \mathrm d\mathbf u/\mathrm dt=\nabla\!\left(-2D^2\,\nabla^2\!\sqrt{\rho}/\sqrt{\rho}\right), whose bracket is the Bohm quantum potential per unit mass, Q/m, the instant one sets D=\hbar/2m (then 2D^2=\hbar^2/2m^2 and Q=-\frac{\hbar^2}{2m}\,\nabla^2\!\sqrt{\rho}/\sqrt{\rho}). So Q is the inertia of the osmotic flow, and the convective term (\mathbf u\cdot\nabla)\mathbf u is what produces it. To turn this identity into a force on a physical density, Simeonov uses two fluids: fluid 2 mutually diffuses down fluid 1’s gradient, \mathbf J_2=-D\nabla\rho_1, and is re-equalised to it on a fast timescale, so its osmotic acceleration is the expression above with \rho\to\rho_1; by Newton’s third law its reaction on fluid 1 is exactly \mathbf F_\mathrm{reaction}=-\nabla Q (Equation 2). Fed into fluid 1’s Euler equation beside an ordinary potential U, this returns the full Madelung system — hence, through \Psi=\sqrt{\rho_1}\,e^{iS/\hbar}, the Schrödinger equation; a relativistic diffusion law J^\mu=D\,\partial^\mu\rho runs the identical steps to Klein–Gordon. For Step 1 the framework reads fluid 1 as the co-rotating dc1 bulk and fluid 2 as a counter-rotating HVBK boundary layer, so that his diffusion constant becomes the circulation-set diffusivity D_{sf}=\kappa_q/(4\pi\alpha_{mf}) of §1.1 and Simeonov’s fast density-equalising “jumps” become the boundary layer continuously re-tracking the bulk gradient. The quantum potential and \hbar=2mD are then a single statement read at two levels: Simeonov’s kinematics, and the substrate’s hydrodynamics.↩︎

  2. Here \kappa_q/4\pi is the standard local-induction self-velocity coefficient of a quantized vortex filament (Saffman 1992; Donnelly 1991), \kappa_q = h/m_\mathrm{eff} the Onsager–Feynman circulation quantum, and \alpha_{mf} the HVBK mutual-friction coupling (Hall and Vinen 1956; Bekarevich and Khalatnikov 1961). The expression is a dimensional identification of the friction-limited transport scale — its content is the scaling D_{sf}\propto\kappa_q/\alpha_{mf} with the 4\pi local-induction prefactor — not a quoted closed-form coefficient. (The distinct thermal vortex-diffusion coefficient has the Einstein form D\propto\alpha_{mf}\,k_BT (Mehdi et al. 2023); what enters here is the athermal, circulation-set transport scale.)↩︎

  3. The value \alpha_{mf}=0.3008 comes from the Weinberg angle using a friction ratio equation. The HVBK mutual-friction force has two pieces: a reactive part that deflects back into the core and a dissipative part that crosses. Both are loss-free interactions. Identifying the reactive channel with the SU(2)_L weak coupling g and the dissipative channel with the U(1)_Y hypercharge coupling g' makes the friction ratio the coupling ratio, \alpha_{mf}=\frac{g'^2}{g^2}=\tan^2\theta_W,\qquad\text{equivalently}\qquad \sin^2\theta_W=\frac{\alpha_{mf}}{1+\alpha_{mf}}, so the measured \sin^2\theta_W=0.2312 fixes \alpha_{mf}=0.2312/0.7688=0.3008. This sits in the range real two-fluid superfluids show in their active regime (He-II passes \alpha\approx0.3 near T\approx2.0 K, T/T_\lambda\approx0.93): \approx30\% of each boundary interaction dissipates, \approx70\% deflects. The same parameter also represents a vortex-scattering phase, \alpha_{mf}=\tfrac12\sin2\delta_0 with \delta_0=18.48^\circ (the weak-scattering branch, selected by \alpha\ll1) — the geometric input behind the fine-structure relation of §5.2. Currently \alpha_{mf} comes from a measured input. It might be possible to compute it from the equilibrium spin rate of the counter-rotating boundary.↩︎

  4. m_\mathrm{eff} is taken to be a property of the substrate, not of the particle it builds — one quantum of dc1 circulation, the same in the electronic and nuclear sectors. Quoting m_e = 0.511\;\mathrm{MeV}/c^2, Equation 3 fixes it once and for all: m_\mathrm{eff} = m_e/\alpha_{mf} = 0.511/0.3008 = 1.70\;\mathrm{MeV}/c^2. Applying Equation 3 in the nuclear sector with this same quantum then defines the proton coupling \alpha_{mf}^{(N)} = m_p/m_\mathrm{eff} = 1836\times0.3008 \approx 552 (552\times1.70\;\mathrm{MeV} \approx 938\;\mathrm{MeV} = m_p c^2). Seen as a single coupling, 552 would break the ceiling \alpha_{mf}\le\tfrac12; seen as a count, it is N_p\,\alpha_{mf} with N_p\approx1836 boundary seams, each with the electron’s per-seam leak. This is a re-parametrization of m_p, not a prediction of it: N_p is back-solved from the measured proton mass, so the identity N_p/N_e = m_p/m_e is algebraic. It would become a genuine prediction only if the seam count could be fixed independently from the topology of the nuclear vortex junction, the way \tan^2\theta_W fixes the per-seam leak. That derivation is left open (§7).↩︎

  5. The observation about the spectrum comes from Volovik, The Universe in a Helium Droplet (2003), §7.4.3–7.4.4. The Bogoliubov quasiparticle energy of a pair-correlated fermionic vacuum is E^2 = M^2(\mathbf p) + c^2 p^2 with c=\Delta_0/p_F (his eqn 7.49); in the strong-coupling deep-BEC limit mc^2\gg|\mu| this deforms continuously into the relativistic Dirac form E^2=\mu^2+c^2p^2 (eqn 7.51), with the chemical potential \mu playing the role of the rest mass. The framework’s one added identification is that the dc1 coherence length is the dc1 Compton wavelength, \xi=\hbar/(m_1 c) — equivalently \Delta_0\simeq\mu=m_1c^2 and p_F\simeq\hbar/\xi, so c=\Delta_0/p_F=\hbar/(m_1\xi). This same \xi=\hbar/(m_1c) is the close-packing input of Route 1 (Equation 8). Strong coupling is also what makes the cone isotropic: the spectrum depends on |\mathbf p| alone, so c_\parallel=c_\perp — unlike weak-coupling ^3He-A, whose point-node spectrum is strongly anisotropic (c_\perp/c_\parallel\sim10^{-5}; Volovik §7.4.6). A single isotropic cone shared by every sector is the origin of c_\mathrm{GW}=c.↩︎

  6. Two steps. (i) Existence condition 1 fixes the threshold wavelength at one cell, \lambda_\mathrm{min}=\xi (the perturbation envelope has no smooth J_1–K_1 match below \xi). A modon carries photon energy E=h\nu=hc/\lambda, so the one-cell modon has E_\mathrm{min}=hc/\xi. (ii) Substituting \hbar/\xi=m_1c (the Volovik relation, Equation 5) gives the clean identity E_\mathrm{min}=hc/\xi=2\pi(\hbar/\xi)c=2\pi\,m_1c^2: exactly 2\pi times the dc1 rest energy m_1c^2 — the 2\pi is one full circulation loop of the cell (the dc1 circulation quantum being \kappa_1=2\pi\hbar/m_1), m_1c^2 the rest energy of the single dc1 wavefunction per cell. Numerically m_1=m_\mathrm{eff}/\nu=2.03 meV gives E_\mathrm{min}=12.8 meV at \lambda=\xi=97\,\mum, i.e. a crossover frequency c/\xi\approx3.1 THz.↩︎

  7. The elimination proceeds with three steps. Volovik’s speed (Equation 5) gives m_1=\hbar/(c\xi); the composition n_1m_1=\rho_\mathrm{DM} then gives n_1=\rho_\mathrm{DM}\,c\,\xi/\hbar; and close-packing n_1\xi^3=1 gives \rho_\mathrm{DM}\,c\,\xi^4/\hbar=1, i.e. \xi=(\hbar/\rho_\mathrm{DM}c)^{1/4}. The only physical content beyond \hbar,c,\rho_\mathrm{DM} is the quantum-degeneracy condition n\lambda_{dB}^3\sim1 that defines Bose condensation — one coherent wavefunction per de Broglie volume — read here at the lattice cell, \lambda_{dB}\to\xi. Numerically \hbar/(\rho_\mathrm{DM}c) = 1.055\times10^{-34}/(2.254\times10^{-27}\cdot2.998\times10^8) = 1.561\times10^{-16}\,\mathrm{m}^4, whose fourth root is 111.8\,\mum. The “\approx1” is exact only for one wavefunction per cell; the true occupancy is the cell occupancy f, so Route 1 is the f\!=\!1 idealization and the real cell is smaller by f^{1/4}.↩︎

  8. Planck reports the dimensionless cold-dark-matter density \Omega_c h^2 = 0.1200 \pm 0.0012 (Planck Collaboration et al. 2020); the mass density follows from the critical density, \rho_\mathrm{DM} = \Omega_c h^2 \cdot \rho_\mathrm{crit}/h^2 = 0.1200 \times 1.878\times10^{-26}\;\mathrm{kg/m}^3 = 2.254\times10^{-27}\;\mathrm{kg/m}^3, with \rho_\mathrm{crit}/h^2 = 3H_0^2/(8\pi G h^2) = 1.878\times10^{-26}\;\mathrm{kg/m}^3. The \pm1\% uncertainty on \Omega_c h^2 carries directly to \rho_\mathrm{DM}.↩︎

  9. The circulation \kappa_q=2\pi\hbar/m_\mathrm{eff} is the ordinary quantum of circulation for a quantum of mass m_\mathrm{eff} (§1.1); the Compton frequency \Omega_v=2m_\mathrm{eff}c^2/\hbar is twice its rest energy over \hbar, the zitterbewegung rate. Explicitly \kappa_q\Omega_v=(2\pi\hbar/m_\mathrm{eff})(2m_\mathrm{eff}c^2/\hbar)=4\pi c^2, with m_\mathrm{eff} cancelling — which is why SC2 is a normalization identity, not a dynamical equation: it fixes which 4\pi the lattice must present in order to gravitate. The field-theory reason that coefficient is precisely the induced-Newton 4\pi is that one-loop induced gravity (Seeley–DeWitt) and Jacobson’s horizon thermodynamics both land on the same number.↩︎

  10. The acoustic metric is Unruh’s (Unruh 1981): linearizing the Euler + continuity equations of a barotropic fluid about a background flow \mathbf v shows that sound — here the modon — obeys a Lorentzian metric whose null cones are tilted and dragged by the flow; Barceló–Liberati–Visser made it general and rigorous (Barceló et al. 2005) and Volovik identified it with superfluid quasiparticles (Volovik 2003). The overall conformal factor \rho/c_s that multiplies the bracket in the general expression has been dropped because it rescales ds^2 but leaves the null cones alone: every light path and the entire causal structure is unchanged, and it divides out exactly when the metric is read in the rain-frame form Equation 14. The substrate identifications are c_s\to c (the modon speed Equation 5) and \mathbf v\to v_\mathrm{ebb}(r)\,\hat r (the radial dc1 inflow). This is the kinematic half — the metric given the flow; the dynamical half, that the substrate’s own equations produce exactly this flow, is carried by Equation 14 below.↩︎

  11. Two independent statements agree. Kinematics: substituting the free-fall profile v_\mathrm{ebb}=\sqrt{2GM/r} into Equation 13 returns the Painlevé–Gullstrand “rain” form of the Schwarzschild line element — the 1921–22 coordinates of Painlevé and Gullstrand, related to the textbook Schwarzschild coordinates by a redefinition of t alone (no curvature invariant is changed), read physically as Hamilton & Lisle’s river model (Hamilton and Lisle 2008). Dynamics: that profile is not an assumed force law but the steady state of the substrate’s own Euler equation. For v_\mathrm{ebb}=\sqrt{2GM/r} the convective acceleration is v_\mathrm{ebb}\,dv_\mathrm{ebb}/dr=-GM/r^2 with dP/dr=0 — free fall with no pressure gradient, exactly as the equivalence principle requires of a comoving frame. The two halves close a fixed point: the flow that sets the metric (kinematics, Equation 13) is the flow the substrate dynamics produce, so the identification is self-consistent rather than circular. Continuity is the one piece not yet closed. At uniform density the steady ebb needs a distributed sink \propto r^{-3/2}. The leading candidate reads the ebb as the superfluid half of a zero-net-flux two-fluid counterflow, in which only the total mass current is conserved and the two components convert locally at exactly that rate. It holds only if gravitational waves see the superfluid velocity, which the GW170817 arrival time tests (§7).↩︎

  12. Each is a geodesic of Equation 14. Redshift is the g_{00} of a static clock, \nu/\nu_0=\sqrt{1-2GM/rc^2}. The deflection \Delta\theta=4GM/(bc^2) carries the GR factor of two over the Newtonian-corpuscular 2GM/(bc^2) because the modon is both refracted by the gravitationally-graded signal speed c(r)=c\,(1-2GM/rc^2) and dragged by the inflow v_\mathrm{ebb} — two equal contributions, where the corpuscular estimate counts only the first. The Shapiro delay \Delta t=(2GM/c^3)\ln(4r_1r_2/b^2) is the same reduced signal speed integrated along the path; the perihelion precession 6\pi GM/(ac^2(1-e^2)) is the non-closure of bound geodesics. The combined moving-clock shift, \nu/\nu_0=\sqrt{1-2GM/rc^2-v^2/c^2}\simeq1-GM/rc^2-v^2/2c^2, is the Painlevé–Gullstrand proper time of a clock moving across the river: the radial inflow and the tangential motion add as one vector, |\mathbf v_\mathrm{rel}|^2=v_\mathrm{ebb}^2+v^2, and one \gamma follows — the GPS correction, +45.7-7.2=+38.5\;\mus/day. The boundary layer’s cost of holding coherence against the flow is a separate line of the energy ledger, \gamma mc^2=mc^2/\gamma+\gamma mv^2 (the beat plus the bow wave); it does not draw down the tick.↩︎

  13. The linearized substrate equations give \Box\,\bar h_{00}=-(16\pi G/c^2)\,\rho_\mathrm{matter} — the 00 component of the linearized Einstein equations — where the naïve fluid-Poisson route would give only 4\pi G. The factor of four is that pressure gravitates: the standard general-relativistic source is the Tolman combination \rho_\mathrm{eff}=\rho+3P/c^2, and the substrate’s stiff equation of state P=\rho c^2 (sound speed c_s=c) makes \rho_\mathrm{eff}=\rho+3\rho=4\rho, so 4\rho\times4\pi G=16\pi G. The substrate’s pressure is the rotational kinetic energy of its vortices. The weight is equation-of-state specific and applies to the substrate’s own self-gravitation; matter embedded in it sources gravity through its own \rho_\mathrm{eff} — radiation 2\rho (P=\rho c^2/3), dust \rho (P\simeq0) — the standard kinematics/dynamics split of analog gravity (Barceló et al. 2005). The 4\pi here is the same Gauss’s-law solid angle that normalizes SC2 (Equation 9) and heads the cell occupancy (Equation 10).↩︎

  14. The substrate’s stiff fluid nature, w=1, does not wreck BBN and the CMB because Gibbs–Duhem pins its background value \bar P\to0 at the marginal point (§3.3) while only the slope \delta P/\delta\rho=c^2 stays stiff — the logarithmic equation of state decouples the two. So the substrate gravitates to Friedmann as w\approx0, c_{s,\mathrm{bg}}\approx0 cold dark matter with clustering driven by a tachyonic instability above the marginal density rather than a horizon-scale Jeans pressure, and with the internal 0.776\,c vortex circulation counting as rest-mass energy, not bulk kinetic pressure.↩︎

  15. In a superfluid at T=0 in full equilibrium the Gibbs–Duhem relation forces \varepsilon+P=0, and the vacuum self-tunes both \varepsilon and P to exactly zero — any added vacuum energy shifts the density, hence the chemical potential, driving a flow that relaxes it away; this is Volovik’s reason the bare 10^{120} does not gravitate (Volovik 2003). The observed \Lambda is the small residual the expansion leaves un-relaxed, and its apparent size depends entirely on the reference density. Against the substrate’s own density \rho_\mathrm{DM} the disequilibrium is order unity, \delta T/T_c=\sqrt{\rho_\Lambda/\rho_\mathrm{DM}}\approx1.6 — there is no tuning of the substrate’s state. The notorious 10^{-61.5} appears only against the gravitational Planck density \rho_\mathrm{Pl}=c^5/\hbar G^2: \delta T/T_c=\sqrt{\rho_\Lambda/\rho_\mathrm{Pl}}=3.4\times10^{-62}, and the two readings differ by exactly \sqrt{\rho_\mathrm{DM}/\rho_\mathrm{Pl}}=(m_1/M_\mathrm{Pl})^2 — the same hierarchy (M_\mathrm{Pl}^2=\hbar c/G) that makes G weak, so small \Lambda and weak gravity have the same origin. Independently, the close-packing relation \rho_\mathrm{DM}c^2=(m_1c^2)^4/(\hbar c)^3 fixes the dark-energy scale with no new input, \rho_\Lambda^{1/4}=2.24\;\mathrm{meV}\approx m_1c^2, so \rho_\Lambda\sim\rho_\mathrm{DM} is the same substrate scale.↩︎

  16. The Landau critical velocity is the substrate’s outer-rim lattice rotation v_L=\omega_0\xi, with \omega_0 the angular rate at the cell vortex’s outer rim; from the condensation number, v_L=c\,(4\pi/\nu)^{1/3}\approx740\;\mathrm{km/s}\approx0.0025\,c — the threshold above which flow through a boundary excites vortices and dissipates. It gates the phase by velocity, not by position: dwarf galaxies (30–80 km/s) and the Milky Way (\lesssim200 km/s) sit below it and stay coherently superfluid (deep MOND), while galaxy clusters (800–1500 km/s) sit above and revert to collisionless CDM — the same substance in two regimes, which is why MOND phenomenology appears in galaxies but not clusters, with no spatial boundary invoked. It is the substrate-ladder sign rule read gravitationally: a system’s speed relative to v_L fixes its phase before it is measured.↩︎

  17. The boundary’s parity-even current–phase relation has no linear term, J\propto|\nabla\Phi|\,\nabla\Phi, and coarse-grained over the N=r/\xi boundaries between source and test point gives the Bekenstein–Milgrom field equation \nabla\!\cdot(|\nabla\Phi|\,\nabla\Phi)\propto4\pi G\rho_b; the Newtonian linear term is induced only once the Hubble flow’s DC bias breaks the parity. So the scale of a_0 comes out where every factor is a substrate parameter — c=\hbar/(m_1\xi) (Equation 5), G from the boundary ebbing current (measured; its unit-free content \mathcal F of Equation 15 is not yet derived), \rho_\mathrm{DM}=n_1m_1. In cell units, today, a_0=f_\mathrm{cross}\,c^2/(\sqrt{4\pi}\,\xi) (Equation 16), so deriving a_0 and deriving G are the same problem. (The identity is present-epoch only: f_\mathrm{cross} and G are fixed by the cell, while a_0 tracks the cosmic-mean \rho_\mathrm{DM}, which dilutes.) Equivalently a_0=c/t_\mathrm{ff} with t_\mathrm{ff}=1/\sqrt{G\rho_\mathrm{DM}} the substrate’s gravitational free-fall time. Through the Friedmann relation this restates Milgrom’s long-standing coincidence 2\pi a_0\approx cH_0 with a definite prefactor, a_0/cH_0=\sqrt{3\Omega_\mathrm{DM}/8\pi}=0.178 (Gauss’s 8\pi, geometry’s 3, and \Omega_\mathrm{DM}), against g_\dagger/cH_0=0.183 for Planck’s H_0 — the same 3\% as a_0 itself, since it is the same relation. At z=0 this is a relabeling — the cH_0/2\pi and cH_0/6 forms also sit inside the \pm20\% systematic band of McGaugh et al. (2016) — and the discriminating content is the redshift law: a_0\propto (1+z)^{3/2} tracks \rho_\mathrm{DM}, reaching \approx2.8\,a_0(0) at z=1 where an a_0\propto cH reading gives \approx1.8\,a_0(0).↩︎

  18. The vertical period is derived from the wavelength of the Glaberson–Johnson–Ostermeier instability (Glaberson et al. 1974) of the vortex lines that thread the stack. Axial flow along a rotating vortex array destabilises its Kelvin (vortex-wave) modes once it exceeds the slowest wave speed; the GJO dispersion \omega=2\Omega+\nu_s k^2 has its phase-velocity minimum at k_c=\sqrt{2\Omega/\nu_s}, so the unstable mode selects a single wavelength d=1/k_c=\sqrt{\nu_s/(2\Omega_\mathrm{sheet})}. Two standard classical-fluid inputs close it with no new coefficient: Feynman’s relation at one vortex per cell, \Omega_\mathrm{sheet}=\kappa_q/(2\xi^2) (Donnelly 1991), and Saffman’s vortex-line tension \nu_s=(\kappa_q/4\pi)\ln(\xi/\xi_\mathrm{GP}) (Saffman 1992). The circulation quantum \kappa_q cancels, leaving the pure geometric ratio of Equation 18. The short-distance cut-off is the close-packing length \xi_\mathrm{GP}=\xi/\sqrt2 — cores touching, the same pairing-two as the ladder below and the bridge equation’s 1/\sqrt2 — so \ln(\xi/\xi_\mathrm{GP}) =\tfrac12\ln2 and d_\mathrm{GJO}=\xi\sqrt{\ln2/(8\pi)}=0.166\,\xi. On the cell size, \xi=97.0\,\mum (Equation 11), that is d_\mathrm{GJO}\approx16\,\mum; the naive f\!=\!1 length \xi_\mathrm{CP}=112\,\mum would give \approx19\,\mum, so the sheet spacing inherits the same f^{1/4} that separates the two.↩︎

  19. Discrete scale invariance is what a scale-free system shows once it is handed a short-distance cutoff: the continuous rescaling symmetry breaks to a discrete subgroup, and a geometric tower of states spaced by a fixed ratio appears. The textbook instance is the Efimov effect — three particles at the unitary limit feel the unique scale-free 1/r^2 potential and would carry a continuum of bound states; a three-body cutoff discretises it into a tower with the universal ratio \lambda_0=e^{\pi/s_0}\approx22.7, set by a renormalisation-group limit cycle (Braaten and Hammer 2006). Sornette finds the same log-periodic signature in classical rupture and earthquakes (Sornette 1998). The substrate is built to sit at such a point — the gapless, isotropic \mu\to0 end of its Bogoliubov spectrum, the same condition under which it emits light at one speed — cut off below by the core \xi_\mathrm{GP}. The full derivation behind the tower effect is still open. The substrate’s ratio is found from the pairing geometry rather than computed from a limit cycle as Efimov’s 22.7 is. If we assume the tower exists: the fundamental rung \xi_\mathrm{GP}\to\xi sits on the marginal dispersion’s E\sim1/r^2 branch, so rescaling length by \sqrt2 rescales energy by 2 — the pairing/breathing octave — and the log-period is \ln\sqrt2=\tfrac12\ln2, the same \tfrac12\ln2 the packing factor carries. What is owed is that the marginal point is a genuine limit cycle at all (the anti-phase pairing driving the effective coupling past the “fall-to-the-center” threshold) — the same marginal-point Bogoliubov problem listed with the discrete-scale-invariance item in the open problems (§7).↩︎

  20. The “most irrational” status of \varphi=(1+\sqrt5)/2 is precise: its continued-fraction expansion is all 1s, so its rational convergents approach it more slowly than for any other number — it is the value hardest to approximate by a ratio, hence the orientation a structure adopts to avoid every rung at once. Levitov derived exactly this for a lattice of magnetic flux tubes with a competing length scale: the energetically stable configuration locks onto the golden mean (Levitov 1991) — the same anti-resonance selection the framework reads into phyllotaxis’s golden angle and the resting cortex’s band spacing (Pletzer et al. 2010). The gap appears as \varphi only on a circle (angles mod 2\pi); the identical avoid-overlap principle takes a different arithmetic in a different space — disordered-hyperuniform “blue noise” for points on a plane (the retinal cone mosaic), and mutually-prime periods for cycles in time, which is why 13- and 17-year cicadas minimise predator co-emergence (Goles et al. 2001). One principle (maximal incommensurability), three geometries set by what the variable lives in.↩︎

  21. A mass with angular momentum J makes a gyroscope at radius r precess at the Lense–Thirring rate \omega_\mathrm{fd}=2GJ/(c^2r^3). In the substrate this is the azimuthal part of the gravitational inflow (Step 3) entrained by the rotating boundary layer, through the same mutual friction that couples the normal and superfluid components of rotating helium (Hall and Vinen 1956). The mechanism is microscopic, not a numerical correction, so for slow rotators it reproduces GR exactly — Gravity Probe B measured Earth’s frame-dragging at 37.2\pm7.2 mas/yr against GR’s 39.2 mas/yr. The framework departs from GR only where the substrate’s own critical (Landau) velocity begins influencing gravity: the azimuthal entrainment cannot exceed v_L\approx750 km/s, so the inferred frame-dragging frequency should plateau rather than diverge as J/M approaches the maximal-spin regime — testable in the inner accretion flow of high-spin AGN via iron-line spectroscopy.↩︎

  22. Physically the Higgs VEV (v) is the baseline strength of the Higgs field that permeates all space. It’s the amplitude of the aligned ground state, v=\sqrt{8\pi}\,m_\mathrm{eff}c^2\sqrt{\nu}, built from the effective quantum m_\mathrm{eff}=m_e/\alpha_{mf}=1.699 MeV and the condensation number \nu=m_\mathrm{eff}/m_1=8.349\times10^8 — the count of dc1 quanta (m_1=2.03 meV) in one effective quantum (the same \nu as Equation 12). The \sqrt\nu is the relativistic Bose-field amplitude law: a complex field’s amplitude-squared counts quanta linearly (|\phi|^2=n/2\omega, so v^2=n_\chi/\omega_\chi is the in-quadrature sum of \nu chirality fluctuations), the same reason every condensate order parameter scales as \sqrt{\text{density}}. The prefactor comes from 8\pi=2\times4\pi_\mathrm{SC2} — the Gauss 4\pi behind \nabla^2\Phi=4\pi G\rho (the Einstein–Hilbert normalization, with emergent Lorentz invariance) times the radiation-equation-of-state weight \rho_\mathrm{eff}=\rho+3P/c^2=2\rho of the massless chirality Goldstones. The chain m_\mathrm{eff}\!\leftarrow\!\sin^2\theta_W and \nu=m_\mathrm{eff}/m_1 never references v, and every step is dimensionally clean and unit-invariant. The \nu it uses is the one built from cosmology and the geometric cell occupancy, \nu=\xi_\mathrm{CP}f^{1/4}/\bar\lambda_\mathrm{eff}=8.349\times10^8 (Equation 12), not from the electroweak side. So the inputs are: measured \rho_\mathrm{DM}, measured \sin^2\theta_W (\overline{\mathrm{MS}} at M_Z), measured m_e (pole mass), one Bessel zero (in the clamped-cell constant K=j_{11}^2+1); out comes 246.12\pm0.33 GeV against a measured 246.22. The error bar is the propagated 1\sigma on \Omega_ch^2 and \sin^2\theta_W. Since v\propto m_\mathrm{eff}^{3/2}\rho_\mathrm{DM}^{-1/8}K^{-1/8}, the result is far more sensitive to the electroweak inputs (\sin^2\theta_W and m_e) than to the cosmological one or to the modon geometry; the scheme, scale and modon-constant dependence is discussed in §5.1.↩︎

  23. Koide’s 1981 relation folds the three charged-lepton masses into Q=(m_e+m_\mu+m_\tau)/(\sqrt{m_e}+\sqrt{m_\mu}+\sqrt{m_\tau})^2, measured 0.666664\pm0.000005 (PDG 2026 (Takahashi et al. 2026), m_\tau=1776.93\pm0.09 MeV, which dominates the error). A generation is read as which of the three branches of one three-fold (Y-junction) vortex node a knot threads — the same \mathbb{Z}_3 that gives color and the \pm\tfrac23 quark charges (§5.4) — so the three \sqrt{m_k} sit at a common mean times (1+A\cos\theta_k) with \theta_k=\delta+\tfrac{2\pi}{3}k. Three unit vectors 120^\circ apart give \sum\cos\theta_k=0 and \sum\cos^2\theta_k=\tfrac32, so the phase cancels and Q=\tfrac13+A^2/6 for any \delta. The lattice’s pairing-\sqrt2 — the same factor in \xi^2=2\xi_\mathrm{GP}^2 — quantizes the deviation amplitude at A=\sqrt2, giving Q=\tfrac13+(\sqrt2)^2/6=\tfrac23 exactly. Since Q and A are one-to-one, this is a reading of Koide’s relation, not a separate prediction. The measured Q corresponds to A=1.414209\pm0.000011, against \sqrt2=1.414214. Equivalently the triad tilts 45^\circ (\cos45^\circ=1/\sqrt2) from the democratic axis (1,1,1) — the exact midpoint between mass degeneracy (Q=\tfrac13) and full hierarchy (Q=1) — and of N equal phases only N=3 lands the balanced Q_N=2/N on \tfrac23, which is also why three generations.↩︎

  24. Derived at Equation 7: the modon matching condition (Equation 6) has no interior/exterior solution below one cell, so the smallest localized modon is exactly \xi wide and carries a floor energy E_\mathrm{min}=hc/\xi; the Volovik relation \hbar/\xi=m_1c (Equation 5) collapses this to the parameter-free 2\pi m_1c^2, with the 2\pi one full cell circulation (\kappa_1=2\pi\hbar/m_1) and m_1c^2 the rest energy of the single dc1 quantum per cell (m_1=2.03 meV \Rightarrow12.8 meV, \nu_\mathrm{floor}=c/\xi\approx3.1 THz, \lambda=\xi\approx97\,\mum). The floor is grounded by m_1 as the frequency of a crossover, not a spectral edge. The speed of light stays at c while the quantization flips from an ordinary modon above to a delocalized collective winding spread over \lambda\gg\xi below, a photon stretched across many cells.↩︎

  25. Fixing A=\sqrt2 (hence Q=2/3) and the electron scale m_e=\alpha_{mf}m_\mathrm{eff} leaves exactly one undetermined number: the phase \delta where the rigid triad sits on its circle. Each lepton’s \sqrt{m} is then 1+\sqrt2\cos\theta_k, so all charged-lepton ratios ride on that single angle. The data pick \delta=0.222222\,rad =2/9 rad to six figures (pairing-two over three-fold-squared; the empirical 2/9 was first noticed by Brannen in 2006), and that one rational lands two independent ratios at once: m_\mu/m_e=206.77 (observed 206.768, +0.001\%) and m_\tau/m_e=3477.5 (observed 3477.37, +0.003\%). In the tripled phase the identity is clean, 3\delta=\tfrac23=Q, so \delta=Q/3=\tfrac19+A^2/18=\tfrac29. What the geometry supplies is the amplitude A=\sqrt2 (equivalently Q=2/3, Koide’s empirical relation given a reason), “why three,” and the massless edge \delta_\mathrm{edge}=\pi/12 forced by A=\sqrt2; what remains a bet is that the residual phase is exactly 2/9 rad — read off the data, not yet derived from a boundary calculation. A sharpened m_\tau pulling the best-fit phase off 2/9 would break the bet while leaving Q=2/3 intact.↩︎

  26. The cell occupancy locks the pitch angle of the counter-rotating duplex to itself, \tan\alpha_\text{pitch}=f — the strand splits its motion into \sin\alpha across substrate sheets and \cos\alpha within them, and equilibrium sits where that ratio matches the substrate’s own chiral/boundary ratio f. With the chemistry-fixed backbone radius r\approx10.0 Å and base rise h=3.4 Å, the turns-per-pitch ratio is N=2\pi r f/h=10.47 bp/turn with no fitted parameter, versus the canonical 10.5\pm0.1 measured in solution. Equivalently the pitch angle itself: f=0.5666 against the observed \tan\alpha_\text{pitch}=35.7/(2\pi r)=0.5681, a 0.26\% match. The numerical agreement is parameter-free; what remains conjectural is the functional form \tan\alpha_\text{pitch}=f (motivated by the chiral/boundary projection but not yet derived from the substrate Lagrangian) and deriving r itself from substrate constants.↩︎

  27. The closed-cylinder counterpart of the DNA pitch, using the same tangent definition — axial rise per cycle over perimeter per cycle. The 13-protofilament wall’s 3-start lateral helix advances three monomer rises h_\text{mon} axially per full azimuthal turn 2\pi R, so \tan\alpha_\text{3start}=3h_\text{mon}/(2\pi R). The substrate locks this to f/\pi — one Gauss factor of \pi absorbed by the cylinder’s closed topology, the perimeter-to-diameter ratio between the open helix and the closed wall — giving R/h_\text{mon}=3/(2f)=2.647. Observed: R/h_\text{mon}=10.6/4.087=2.594 from cryo-EM, a 2.1\% match well inside the scatter on the wall radius R across structures. The same open derivation as the DNA case applies, though the substrate-Lagrangian backing (the cylindrical reduction 4\pi\to4 of the bridge equation’s Gauss factor) is sharper here.↩︎

  28. \delta_0=18.48^\circ is the s-wave scattering phase of a dc1 quasiparticle off the quantized vortex boundary — the same boundary that supplies \hbar (§1.1) — fixed by the Weinberg angle through \alpha_{mf}=\tfrac12\sin2\delta_0, with the weak-scattering branch selected by the requirement \alpha\ll1. The same phase sets the SU(2) coupling (g^2=4\sin^2\delta_0), so \alpha=g^2\sin^2\theta_W/4\pi follows with no new parameter. At M_Z the chain gives g = 0.634 against the measured 0.652, and 1/\alpha=135.1 against the measured 127.95 at that scale. Run upward with the Standard Model’s two-loop equations, which involve only measured particles, the measured couplings satisfy the relation exactly at \mu^*, where \delta_0=18.9^\circ. Matched at v/\sqrt2 and run to zero momentum with the measured vacuum polarization, it gives 1/\alpha(0)=136.7–137.5, and the result moves by about 0.45\% for every 10 GeV of matching scale. Below \mu^* no substrate correction is needed, and none is allowed: a substrate vacuum polarization at the boundary doublet’s energies (\sim0.04–1.7 MeV) would shift the hydrogen and muonic-hydrogen Lamb shifts and the electron g-2 by 10^4–10^8 times their measured precision. The anomalous moment takes \alpha at its own scale: \alpha(0)/2\pi is the Schwinger term (+0.15\% from the measured moment, the size of the two-loop correction).↩︎

  29. The electron’s internal structure oscillates at the Compton frequency \omega_c=m_ec^2/\hbar, pumping pressure ripples of wavelength \lambda_c=h/(m_ec) into the substrate; the electron’s motion Doppler-compresses them into a pilot-wave envelope of de Broglie wavelength \lambda_B=h/(m_ev). Levels are quantized by requiring the envelope to close on itself constructively around the orbit, 2\pi r=n\lambda_B — the same standing-wave matching that fixes the Larichev–Reznik modon. With Coulomb force balance this gives r_n=n^2a_0 and E_n=-13.6/n^2 eV; the ground state locks since at v=\alpha c the de Broglie wavelength 137\,\lambda_c equals the Bohr circumference 2\pi a_0.↩︎

  30. The split pair stays joined by a half-quantum vortex channel in the substrate’s SU(2)\to U(1) order parameter — a topologically protected line, like those seen in superfluid ^3He-A, whose half-integer (\pi) winding encodes the singlet constraint J_A+J_B=0 and cannot be undone by any local fluctuation. As the two particles separate they sweep the bulk’s counter-rotating layers aside, so the channel interior is a laminar corridor, not the obstacle course that holds modons to c. A torsional twist wave bound to that corridor’s core is sub-emergent structure and runs at v_\mathrm{ch}\gg c — bounded below at 5\times10^6\,c by aligned Bell tests, and plausibly at the mass-hierarchy scale m_e/m_1\approx2.5\times10^8 (an ordinary Kelvin bending wave on the line is capped below c). Measuring A along \hat{\mathbf a} snaps its axis and launches that twist wave carrying the rotation to B; the geometric identity R_A(-\hat{\mathbf s}_0)=-\hat{\mathbf a} erases the hidden axis \hat{\mathbf s}_0 from B’s outcome, leaving E(\theta\mid\hat{\mathbf s}_0)=\sin^2(\theta/2)-\cos^2(\theta/2)=-\cos\theta for every \hat{\mathbf s}_0, not just on average — the local-hidden-variable \tfrac13 dilution is gone — and CHSH then reaches the Tsirelson bound |S|=4/\sqrt2=2\sqrt2. Because A’s outcome is random, B’s marginal averages to 50/50, so no two-party experiment can signal. That does not extend to more parties: any model with finite-speed hidden influences that reproduces quantum mechanics must signal faster than light in aligned multipartite geometries (Bancal et al. 2012). The framework accepts this. The signal runs forward in the CMB frame, so it closes no causal loops, and it appears only for events aligned in that frame to within L/v_\mathrm{ch}. The finite v_\mathrm{ch} is the falsifiable edge — Bell correlations should relax from -\cos\theta toward the classical -\tfrac13\cos\theta for pairs whose two measurements are simultaneous in the substrate frame to within L/v_\mathrm{ch} — a reach L_\mathrm{max}=v_\mathrm{ch}\,\Delta t for an experiment aligned to \Delta t — so an aligned test on a 50–100 km baseline probes the mass-hierarchy scale directly. Full derivation at lightfluid.org.↩︎

  31. A spin measurement reads the projection of the electron’s counter-rotating boundary-layer angular momentum onto the analyzer axis. As the axis tilts by \theta the reactive coupling between boundary layer and detector falls like a gear mesh taken off-axis, and the transmitted weight follows \cos^2(\theta/2) — the single-quantum Malus law that is the Born rule for spin. The half-angle is the same two-layer, double-cover topology that gives the 720° return; carried on the vortex channel joining a split pair, the identical twist-wave geometry returns Bell’s E(\theta)=-\cos\theta and CHSH 2\sqrt2, so measurement statistics and non-locality are one construction seen head-on and at a distance.↩︎

  32. In the substrate the electromagnetic vector potential \mathbf A is the co-rotating flow-velocity field itself — a physically real circulation, not a gauge convenience. A charge carried around a region threaded by magnetic flux \Phi accumulates a phase e\Phi/\hbar from the enclosed circulation even where the local field \mathbf B vanishes (outside an ideal solenoid), exactly as a walker circling a vortex picks up its quantized circulation regardless of the empty core. The shift is therefore expected to be real, topological, and parameter-free, matching the measured Aharonov–Bohm effect; what the substrate reading removes is the “action at a distance” puzzle, since the potential is simply the local flow the charge swims in.↩︎

  33. A rotating apparatus imprints the enclosed substrate circulation 2\vec\Omega\cdot\vec A on any wave that traverses the loop, and the matter-wave fringe count is the number of circulation quanta h/m threading the enclosed area, so the accumulated phase is \Delta\varphi=\tfrac{2m}{\hbar}\vec\Omega\cdot\vec A. This is the frame-dragging g_{t\phi} term (§4.4) read for a laboratory apparatus rather than for a gravitating mass — the same entrained azimuthal flow — so light interferometers, atom gyroscopes, and ring lasers are all reading one substrate circulation. A consistency check on the shared gravitational/quantum flow, not a new number.↩︎

  34. The electron is a dual-spin gyroscope: a co-rotating core and a counter-rotating boundary shell of nearly equal moment of inertia, whose net angular momentum L_\mathrm{core}-L_\mathrm{boundary}=\hbar/2 is the residual. An external field couples to both with opposite sign, so the net coupling is proportional to 2L_\mathrm{spin} — hence g=2 at leading order. A small inertial asymmetry \eta=(I_1-I_2)/(I_1+I_2) shifts this to g_e=2/(1-\eta^2)\approx2(1+\eta^2), supplying the anomalous moment.↩︎

  35. The fractional charges follow from the solid-angle geometry of the co-rotating flow at the three-fold Y-junction: the up- and down-quark orbital orientations couple to the confinement boundary with monopole strengths +\tfrac23 and -\tfrac13, a count reproduced independently by ribbon-twist counting in the preon-braid picture. The Schwinger term is the g=2 asymmetry made quantitative: the electron’s own electromagnetic self-energy perturbs the core/boundary mass split at order \alpha, with a 1/2\pi shell-averaging factor (l=1 Fourier component over the spherical shell), so \eta^2=\alpha/2\pi and (g-2)/2=\alpha/2\pi.↩︎

  36. A Laughlin quasiparticle in the fractional quantum Hall effect is a vortex in the two-dimensional electron fluid, and its charge is its winding fraction — so the measured e/3 (shot noise, 1997) is charge as a vortex winding number, the same ontology the framework uses for the quark’s +\tfrac23,-\tfrac13 from the three-fold junction solid angle (see the quark-charge note in the Tier 2 list): the fraction \tfrac13 appears twice for one reason, once minted at a vortex junction and once measured in a 2D electron gas. The integer plateaus \sigma_{xy}=\nu e^2/h are the topological protection of quantized circulation — why the von Klitzing resistance standard holds to <10^{-9}. The same resistance fixes \alpha as an impedance ratio, \alpha=Z_0/2R_K with Z_0=\sqrt{\mu_0/\varepsilon_0} and R_K=h/e^2: it compares how the lattice carries a wave with how it carries one quantized twist, and its running is the screening of that twist by short-lived twist–antitwist pairs, ordinary vacuum polarization. The odd/even split in the fractional sequence is the framework’s boundary-parity fermion/boson rule (§4.5) read on flux attachment: odd denominators are the Jain composite-fermion sequence, while the even-denominator \nu=\tfrac52 state is the paired breath (Moore–Read / Read–Green p-wave BCS) — the same anti-phase pairing that binds a Cooper pair (see the Cooper-pair note), here in the quantum Hall fluid. A shared ontology (charge = winding, parity = pairing), not an identical microscopic mechanism.↩︎

  37. The equation M_W/M_Z=\cos\theta_W follows from the gauge structure shared by the Standard Model and the substrate. The substrate adds a meaning behind \cos\theta_W. Writing the Weinberg angle as \sin^2\theta_W=1-\beta_\mathrm{eq}^2=1/\gamma_\mathrm{eq}^2 identifies \cos\theta_W=\beta_\mathrm{eq} as the equatorial spin velocity of the fermion’s counter-rotating boundary shell (§1.1) — an oblate spheroid of polar-to-equatorial axis ratio R_\mathrm{polar}/R_\mathrm{eq}=0.548 and eccentricity 0.837, the deformation expected of a shell spinning near the Compton frequency, with g (stronger) coupling to the tighter polar curvature and g' (weaker) to the gentler equatorial curvature. This \beta_\mathrm{eq}=0.877\,c is distinct from — and larger than — the inner rotation speed \approx0.776\,c the substrate itself turns at (the c\sqrt{2\alpha_{mf}} inner-scale orbital velocity): the boundary shell spins faster than the effective quantum it confines, as a confining structure should. Numerically \cos\theta_W=\sqrt{1-0.2312}=0.877 against the measured M_W/M_Z=80.4/91.2=0.882, a 0.5\% match. The small residual is a renormalization-scheme difference: \sin^2\theta_W=0.2312 is the \overline{\mathrm{MS}} value at M_Z, while the on-shell definition \sin^2\theta_W\equiv1-M_W^2/M_Z^2\approx0.223 makes the relation exact by construction.↩︎

  38. A mixing matrix is the relative rotation between two families’ \mathbb{Z}_3 generation clocks — the same three-fold junction whose branches are the three generations (§5.1 and its Koide note). In the quark sector the up-type and down-type towers are both color-loaded onto the same side of the generation lock, so their clocks are nearly co-aligned and the residual rotation — the CKM matrix — is small. In the lepton sector the charged leptons and neutrinos sit on opposite sides and straddle the lock, so their clocks are anti-aligned and the residual rotation — the PMNS matrix — is large. The sign of that asymmetry (small CKM, large PMNS) is the parameter-free content. Two further features follow: the atmospheric angle sits at the lattice’s pairing-\sqrt2 tilt, \theta_{23}\approx45° with \cos45°=1/\sqrt2 — the same \sqrt2 that fixes the Koide amplitude (§5.1) — reproduced within the measured 43–49°; and because each clock is the complex cube-root-of-unity \omega=e^{2\pi i/3}, a large CP-violating Dirac phase (the imaginary part of the clock) is generic in both sectors, so the framework expects a near-maximal \delta_\mathrm{PMNS}. What remains owed is the individual mixing angles and the numerical CP value, which need the full neutrino-oscillation calculation (§7).↩︎

  39. The electron carries a real internal oscillation at the Compton frequency \omega_c=m_ec^2/\hbar (see the hydrogen note) — a substrate breath, not a bookkeeping clock — and every signal that keeps it ticking is bounded by the one substrate signal speed c. A clock built on c-bounded signals and carried at speed v is the textbook transverse light-clock: it slows by \gamma=1/\sqrt{1-v^2/c^2}, lengths along the motion contract by 1/\gamma, and the energy of the moving breath is E=\gamma mc^2. The de Broglie wave is that same internal oscillation viewed from the lab frame — the Compton breath beat sideways — so relativistic kinematics and the pilot wave are one structure. This makes special relativity a consequence of an internal clock in a medium with a single signal speed, and unifies the Michelson–Morley null (no preferred frame is detectable because every clock and rod is built from the same c-bounded substrate) with the pilot-wave electron. Exact and zero-parameter; the lone substrate assumption — that the internal clock is physically real — is itself a live prediction (electron-channeling / Gouanère-type searches for the Compton oscillation). The real rest frame survives reciprocity: observers who synchronize their clocks with light signals set the forward clock behind by vD/c^2, and with those dials each crew reads the other’s clocks slow by the same \gamma (Bell 1987). The moving particle’s energy splits as \gamma mc^2=mc^2/\gamma+\gamma mv^2: the slowed beat plus the bow wave the pusher paid for, not a budget the tick draws down.↩︎

  40. The counter-rotating boundary has even parity (equal +\omega_0 and -\omega_0 content), which forces a quadratic current–phase response, J\propto|\delta\phi|\,\delta\phi. In the continuum gravitational limit this is the deep-MOND field equation \nabla\!\cdot[|\nabla\Phi|\,\nabla\Phi]\propto4\pi G\rho_b, whose spherically symmetric solution gives flat rotation curves with v^4=a_0GM_b — the baryonic Tully–Fisher law. The acceleration scale a_0=c\sqrt{G\rho_\mathrm{DM}} is fixed by substrate parameters, not fitted to the relation.↩︎

  41. Minimizing a coherent pulse’s reflection at each vascular junction derives the area-preserving branching rule that West, Brown & Enquist had to assume, and carrying it through the tree returns Kleiber’s 3/4 — one instance of the same boundary-leakage theorem that also fixes the lattice’s geometric ladder.↩︎

  42. Inflation here is driven by latent-heat release during a first-order superfluid transition: compression nucleates more superfluid, suppressing the sound speed (c_s^2=\varepsilon_sc^2, \varepsilon_s\ll1) and acting as negative pressure. The tilt then follows the generic slow-roll/k-inflation form n_s-1\approx-3/2N_*-s with s\approx1/2N_*, so n_s\approx1-2/N_*, but the framework makes the e-fold count natural — N_*\approx\ln(c/H_0\xi)\approx69 from the coherence length and Hubble scale, no potential tuning, eased toward 60 by the subluminal bubble wall. That gives n_s\approx0.967 at N_*=60, within 1\sigma of 0.965\pm0.004 (0.971 for the uncorrected 69). The same N_* is the weak-gravity number of §3.2 read cosmologically.↩︎

  43. Binding is a “residual seam”: when nucleon confinement boundaries overlap, their tails fuse into a shared counter-rotating seam, a saturating contact interaction (~8 MeV/nucleon) that is the strong force rather than an exchanged particle. The curve’s shape is the competition of volume (seam saturation), surface (missing edge contacts), and Coulomb (co-rotating boundary repulsion) terms; setting \mathrm d(B/A)/\mathrm dA=0 gives a peak at A_\mathrm{peak}\approx2a_S/a_C, i.e. the iron/nickel region — the largest nucleus where surface-seam binding still outpaces accumulated Coulomb repulsion.↩︎

  44. All electrons share one circulation handedness, so what pairs two electrons is opposite Compton phase: when phonon-mediated distortion brings two electrons into proximity \pi out of phase in their breathing, one is contracted to its inner scale while the other is expanded, and the inside/outside flow gradient between them self-organizes into a shared counter-rotating vortex — structurally the same seam that binds nucleons (§4.1), at a millionth the energy and a billion times the size. The BCS coherence length is not the electron size but the spatial envelope of that anti-phase give-and-take; the two phase-flips make the pair even-parity, hence a boson that condenses. The pair flows around a defect because its net disturbance averages to zero (topological protection), and the macroscopic SQUID-measurable phase is the condensate’s collective breathing clock — the same limit the ladder’s lock pole reaches when a whole ensemble locks breath-for-breath into one clock (§4.3). Two textbook results follow in the same language: the isotope effect T_c\propto M^{-1/2} is the phonon channel-distortion rate \omega_D\propto M^{-1/2}, and the Meissner effect with both London equations is the reactive (B\to0) limit of the same HVBK mutual friction that fixes \alpha_{mf} (§1.1) — the superconductor as the dissipationless corner of the framework’s two-fluid dynamics. The neutral twin, with the charge stripped off, is superfluid helium tuning in to the vacuum’s breath.↩︎

  45. A mass draws a steady inward substrate current with the free-fall profile v_\mathrm{ebb}(r)=\sqrt{2GM/r} of §3.1; fed into the acoustic metric (Equation 13) it is exactly the Painlevé–Gullstrand Schwarzschild line element, so the horizon is the physical surface where the inflow reaches the signal speed, v_\mathrm{ebb}=c\Rightarrow r_s=2GM/c^2. Because the metric is exact, Unruh’s 1981 result follows exactly rather than by analogy: the surface gravity is \kappa=\tfrac12|\mathrm dv_\mathrm{ebb}^2/\mathrm dr|_{r_s}=c^4/4GM, and k_BT_H=\hbar\kappa/2\pi c=\hbar c^3/8\pi GM (for a solar mass, \sim6\times10^{-8} K). The radiated quanta are modons shed by the failing horizon boundary skin — the same pinch-off as atomic emission (Figure 4), driven by a flow gradient instead of an orbital drop. The prefactor 8\pi=2\times4\pi_\mathrm{SC2} is the framework’s gravitational normalization once more: the Gauss 4\pi of \nabla^2\Phi=4\pi G\rho times the radiation-EOS weight 2, the same 8\pi as the Higgs VEV (§5) and the Friedmann H^2=\tfrac{8\pi G}{3}\rho. A re-derivation of a known result, exact because the metric is.↩︎

  46. In this framework entropy is lost pairing coherence counted in boundary crossings: the dissipative fraction \alpha_{mf} of each coherent “breath” handed across a counter-rotating seam disperses incoherently among the lattice cells and cannot be re-gathered. A horizon is the ultimate one-way boundary — everything crossing it takes the dissipative channel to completion (nothing can be re-collected from outside, 1-\alpha_{mf}\to0), so the irreversibly lost information is registered on the seam it had to cross, the horizon area, not the interior volume. That is the substrate reason the law is area, not volume. The coefficient S_\mathrm{BH}=A/4\ell_\mathrm{Pl}^2 (\ell_\mathrm{Pl}^2=\hbar G/c^3) is the same 1/4 the de Sitter horizon carries — the cosmological and black-hole horizons being one kind of surface seen from its two sides, so the two are locked, and it ties to the \Lambda accounting of §3.3. What the framework adds is the reason there is an area law; the value 1/4 is shared, not yet counted independently from substrate cells (the same open calculation as \Lambda’s value).↩︎

  47. The three Sakharov conditions are native. (i) Baryon-number violation: knots are minted wholesale only at the boil — the first-order transition where normal substrate orders into the superfluid — and are topologically locked (unbreakable three-quark knots, hence also proton stability) ever after. (ii) C/CP violation: an antiparticle is a knot with every rotation reversed, and the vacuum condensed into a single chirality (the same fact that makes the weak force left-handed), so swapping co- for counter-rotating changes the knot’s relation to the handed background — CP is broken by a definite amount, the chirality packing factor \varepsilon_\mathrm{chirality}=\sqrt{\pi\ln2}/K=0.0941, fixed by the same modon factor K=j_{11}^2+1 that sets the cell occupancy (Equation 10), the \sqrt{\pi\ln2} the entropy of the one fillable Majorana zero mode per paired vortex. (iii) Departure from equilibrium: the boil proceeds by bubble nucleation, and a wall sweeping outward is equilibrium-breaking by definition. A baryon is a bound three-quark knot, each quark a three-armed Y-junction, so freezing a net handedness is a choice repeated across 3\times3=9 chirality-bearing interfaces, one factor of \varepsilon_\mathrm{chirality} each: \eta_B\approx\varepsilon_\mathrm{chirality}^9=(0.0941)^9=5.8\times10^{-10} versus the observed (6.1\pm0.04)\times10^{-10} (5\%, zero tunable parameters, about 8\sigma). The mechanism is solid, but that the per-interface bias is exactly one factor and the count exactly 9 is a motivated conjecture, not a stress-tensor calculation. The exponent is sharply diagnostic (n=8 would overshoot tenfold), which makes it falsifiable rather than fitted.↩︎

Citation

BibTeX citation:
@misc{vroom2026,
  author = {Vroom, Jeff},
  title = {The {Vacuum’s} {Superfluid} {Lattice}},
  date = {2026-08-12},
  url = {https://doi.org/10.5281/zenodo.23873867},
  doi = {10.5281/zenodo.23873867},
  langid = {en},
  abstract = {I found a model that fits broadly, in which the vacuum is
    a superfluid condensate of a light (\textasciitilde2 meV) particle
    and every particle is a vortex in it whose rotation is set by the
    signal speed c. Two adjustments to established superfluid-vacuum
    physics are needed: the Bohm quantum potential is the reaction force
    of a counter-rotating boundary layer (following Simeonov’s two-fluid
    derivation), and the photon is a modon, a self-propelled
    counter-rotating vortex pair. A motivated choice for the modon’s
    lowest Bessel mode in one cell, together with the measured
    dark-matter density, fixes a lattice cell of ξ ≈ 97 μm with no
    continuously fitted parameters. Measured in units of an effective
    quantum set by the electron mass and the Weinberg angle, the same
    cell gives the Higgs vacuum expectation value as 246.1 ± 0.3 GeV,
    within the 0.13\% uncertainty carried by the dark-matter density,
    and the MOND scale follows as c√(Gρ\_DM) to within 3\% of McGaugh’s
    estimate. The model recovers the Painlevé–Gullstrand form of the
    Schwarzschild metric and sees mass as rotational energy leaking
    through quantized boundaries. Inputs, open factors, and Python
    scripts reproducing the main numbers are included.}
}
For attribution, please cite this work as:
Vroom, Jeff. 2026. “The Vacuum’s Superfluid Lattice.” Preprint, August 12. https://doi.org/10.5281/zenodo.23873867.