Gravity as Boundary-Layer Ebbing

Every orbital system in this framework is wrapped in counter-rotating boundary layers — co- and counter-rotating flows that nearly cancel, forming the barriers that make stable matter possible. These boundaries do three things, depending on how they are disturbed:

1. Push Back

Quantum Potential
Reaction force against internal flow gives particles their wave nature.

2. Leak

Gravity
A tiny ebbing current of substrate particles transits the boundary.

3. Eject

Photon
A counter-rotating vortex dipole (modon) breaks away from the layer.

All three use the same boundary layer, the same dc1 substrate, the same counter-rotating mechanics. The difference is the mode of interaction: reaction force, leak current, or ejection. This chapter is about the second mode — the leak current that we experience as gravity, and why it is so extraordinarily weak.

The Mechanism

Gravity as a Leak Current A source mass on the left drives a dense dc1 current rightward against a counter-rotating boundary layer. The layer is a near-perfect mirror: almost the entire incident current is reflected back, and only a single thin thread — about one dc1 particle in ten quadrillion — transits to reach a distant boundary on the right, delivering the net momentum we feel as gravity. The transit fraction is f_cross = 4 pi G / v_rot,outer, approximately 1.1e-15, which is why gravity is so weak. Gravity as a Leak Current Inside a boundary layer two opposing currents nearly cancel — only about one dc1 particle in a quadrillion transits, and that thread is gravity. M source mass trapped rotational energy incident dc1 current ≈ all of it reflected back counter-rotating boundary layer ≈ neutral zone the leak: net dc1 current carries momentum fcross ≈ 1.1 × 10−15 of the flux gets through m distant boundary F = GMm / r² Of every ~1015 dc1 particles that strike a boundary, about one transits. The layers are near-perfect mirrors — that is the whole reason gravity is weak.

Gravity is a physical flow, not a force at a distance. It arises from the net dc1 current that leaks through the counter-rotating boundary layers of orbital system complexes.

Within a boundary layer, the co- and counter-rotating systems nearly cancel, creating an approximately neutral zone. But a tiny fraction f_\text{cross} of dc1 particles transit the boundary per unit time, carrying momentum from the gravitational source. This “ebbing current” applies force to each boundary as a whole — and since mass is the total rotational energy enclosed by those boundaries (see Mass as Leaking Rotational Kinetic Energy), the force is proportional to the enclosed mass.

Mathematical Form

The River of Space Substrate flows radially inward toward a central mass, speeding up as v_ebb(r) equals the square root of 2GM over r — slow blue arrows far out, fast red arrows near the centre. Where the inflow reaches the speed of light it crosses the horizon at r_s = 2GM/c squared. Substituting this self-consistent flow into the Unruh-Visser acoustic metric yields the exact Painleve-Gullstrand form of the Schwarzschild solution, from which gravitational redshift, light deflection, Shapiro delay, and perihelion precession all follow exactly. The River of Space Space itself flows inward at vebb = √(2GM/r). Feed that flow into the acoustic metric and you get the Schwarzschild solution — exactly, not approximately. M horizon: vebb = c rs = 2GM/c² inflow speed grows toward M far / slow faster v → c From flow to spacetime the self-consistent steady inflow… vebb(r) = √(2GM / r) vebb dvebb/dr = −GM/r²  from substrate Euler + continuity …dropped into the acoustic metric… ds² = −(c² − vebb2) dt² − 2vebb dr dt + dr² + r² dΩ² = Painlevé–Gullstrand Schwarzschild exact — not linearized, not approximate where vebb = c the river outruns sound: an acoustic horizon, rs = 2GM/c² every static GR test follows exactly redshift · light-bending · Shapiro delay · perihelion precession · GPS

The gravitational ebbing current density is:

j_\text{grav} = f_\text{cross} \cdot n_1 \cdot m_1 \cdot v_\text{drift}

where v_\text{drift} is the net drift velocity of dc1 particles between massive systems, driven by the asymmetry created by a mass M at distance r. For this to reproduce Newtonian gravity (F = GMm/r^2), the drift velocity must scale as M/r^2 — which follows if the dc1 current is sourced by the total orbital system energy of the source and falls off as 1/r^2 due to geometric dilution in 3D.

The constraint system gives a simplified form for the gravitational constant:

\boxed{G = \frac{f_\text{cross} \cdot v_\text{rot,outer}}{4\pi}}

The velocity here is the outer-scale lattice rotation v_\text{rot,outer} = \omega_0 \xi \approx 0.0025\,c \approx 7.6 \times 10^5 m/s — the macroscopic flow speed relevant to the gravitational leak current, not the inner-scale orbital velocity. With this velocity:

f_\text{cross} = \frac{4\pi G}{v_\text{rot,outer}} \approx 1.1 \times 10^{-15}

Note

As written, f_\text{cross} = 4\pi G / v_\text{rot,outer} has units [\text{m}^3/(\text{kg}\cdot\text{s}^2)] / [\text{m/s}] = [\text{m}^2/(\text{kg}\cdot\text{s})], not dimensionless as labeled (P8). The numerical value 1.1 \times 10^{-15} is correct in SI. See open problems WIP-15 for the path to repair.

Intermediate steps (full → simplified): The full gravitational ebbing force between masses M and m separated by r is:

F = f_\text{leak} \cdot n_1 \cdot m_1 \cdot v_\text{drift} \cdot A_\text{boundary}

Setting F = GMm/r^2, with v_\text{drift} \propto M/r^2 from geometric dilution, and identifying n_1 m_1 = \rho_\text{DM}, the density and area factors (\rho_\text{DM}, A_\text{boundary}, geometric prefactors from the 2D sheet structure) absorb into f_\text{cross} when the expression is reduced to a single-parameter form using v_\text{rot,outer}. The simplified G = f_\text{cross} \cdot v_\text{rot,outer} / (4\pi) is a numerical recipe valid in MKS — it gives the correct value of G — but the hidden dimensional factors (involving \rho_\text{DM}, the inter-sheet spacing d, and the chirality-coherent 2D→3D projection) have been absorbed into f_\text{cross}, making it appear dimensionless when it is not.

This is the same family of issue as the bridge equation and the old C1 modon condition: the substrate’s vortex lattice is organized into chirality-coherent 2D sheets, and the correct 3D→2D projection requires the thermodynamic calculation of the inter-sheet spacing — the same calculation needed to derive the Higgs VEV. The numerical results and physical interpretation are unaffected; the issue is presentational completeness. See WIP-15.

Gravity’s extraordinary weakness — G \sim 10^{-11} in SI units — is a direct consequence of f_\text{cross} being \sim 10^{-15}: only about one in a quadrillion dc1 particles transits a boundary per interaction time. The counter-rotating layers are nearly perfect barriers.

Recovering General Relativity

The substrate does not merely approximate Newtonian gravity — it reproduces general relativity exactly. The self-consistent steady-state dc1 inflow velocity is:

v_\text{ebb}(r) = \sqrt{\frac{2GM}{r}}

Substituting this into the Unruh-Visser-Volovik acoustic metric gives the exact Painlevé-Gullstrand form of the Schwarzschild solution — not approximate, not linearized. The substrate’s Euler and continuity equations produce this flow self-consistently (v_\text{ebb} \cdot dv_\text{ebb}/dr = -GM/r^2 exactly), closing the loop through a fixed-point argument.

All classical static GR tests — gravitational redshift, light deflection, Shapiro delay, perihelion precession, GPS corrections — are exact consequences of this acoustic metric. The substrate density adjusts hydrostatically as \rho(r) = \rho_0 \exp(-\Phi(r)/c^2), which is automatic for the barotropic equation of state P = \rho c^2.

Rotating bodies: frame-dragging from the azimuthal flow

Frame-Dragging from the Azimuthal Flow A spinning mass with angular momentum J entrains the surrounding substrate into azimuthal flow that weakens outward as 1 over r squared and varies as sine theta. Frame-dragging is just the local angular velocity of this flow. A free draining vortex would conserve circulation, giving v_phi proportional to 1 over r and a dragging rate of 1 over r squared, which is not the Lense-Thirring law. The substrate is instead entrained: the force-free azimuthal flow obeys an Euler equation with roots plus one and minus two, so the decaying solution gives v_phi proportional to 1 over r squared and a dragging rate of 2GJ over c squared r cubed, exactly Kerr to linear order. Gravity Probe B confirms both the geodetic and frame-dragging precessions. Frame-Dragging is the Spinning Substrate A spinning mass entrains the substrate into azimuthal flow. Frame-dragging is nothing but the local angular velocity of that flow — an entrained dipole, not a free-draining vortex. spin axis J vφ  ∝  sinθ / r² strongest at the equator frame-dragging rate = local flow rotation: ωdrag = vφ / (r sinθ) = 2GJ / c²r³ Why 1/r³, not the bathtub's 1/r² free vortex (bathtub): Γ = 2πr vφ const → vφ ∝ 1/r → ω ∝ 1/r² — not the Lense–Thirring law entrained boundary: f″ + (2/r)f′ − (2/r²)f = 0, roots +1, −2 → vφ ∝ 1/r² → ωdrag ∝ 1/r³ — exactly Kerr / Lense–Thirring to linear order Gravity Probe B Precession Predicted GP-B measured Geodetic radial ebb (mass sector) 6604 6601.8 ± 18.3 Frame-dragging azimuthal vφ (spin sector) 41 37.2 ± 7.2 mas/yr · geodetic matches GR to 0.03%; same G as the mass sector, applied to J

Real masses spin, and their gravitational fields carry angular momentum. In the substrate this must show up as an azimuthal entrainment v_\phi added to the radial ebb — the spinning boundary layers dragging the surrounding dc1 into helical flow. The acoustic metric makes the connection exact: the cross term -2\,\vec v\cdot d\vec x\,dt that already carries the gravitomagnetic sector has an azimuthal piece -2\,v_\phi\,(r\sin\theta)\,d\phi\,dt, so g_{t\phi}=-v_\phi\,r\sin\theta and g_{\phi\phi}=(r\sin\theta)^2 (the conformal factor cancels). The frame-dragging (zero-angular-momentum-observer) rate is then

\omega_\text{drag} = -\frac{g_{t\phi}}{g_{\phi\phi}} = \frac{v_\phi}{r\sin\theta}.

Frame-dragging is the local angular velocity of the substrate flow. The radial ebb v_\text{ebb} lives in g_{tt} and gives Schwarzschild; the azimuthal v_\phi is the only new field needed for the rotating (Kerr) sector.

Why the falloff is 1/r^3, not the bathtub’s 1/r^2. The tempting picture — a draining vortex — is wrong. A free vortex conserves circulation \Gamma = 2\pi r\,v_\phi=\text{const}, giving v_\phi\propto 1/r and \omega_\text{drag}\propto 1/r^2, which is not the Lense–Thirring law. The substrate around a spinning body is not freely draining; it is entrained by the rotating boundary. A purely azimuthal, force-free flow v_\phi=f(r)\sin\theta obeys the \phi-component of the vector Laplacian,

f'' + \frac{2}{r}f' - \frac{2}{r^2}f = 0,

an Euler equation with indicial roots n=+1 (the rigid co-rotating interior) and n=-2 (the decaying exterior). The unique decaying entrainment is therefore

v_\phi(r,\theta) = \frac{2GJ}{c^2}\,\frac{\sin\theta}{r^2} \quad\Longrightarrow\quad \omega_\text{drag}(r) = \frac{2GJ}{c^2 r^3},

which is exactly Kerr / Lense–Thirring to linear order. This is the spin-dipole (\ell=1) analog of the mass-monopole’s Newtonian 1/r: a mass sources a monopole (the radial ebb \sqrt{2GM/r}), an angular momentum sources a dipole (the azimuthal rotlet). The radial law is fixed by geometry alone — zero parameters, the same status as the Schwarzschild profile; only the amplitude carries a coupling, and it is the same G as SC1, applied to the boundary’s angular-momentum current J rather than its mass M.

Gravity Probe B. Feeding Earth’s parameters into this picture, the geodetic precession (from the radial SC1 sector and orbital motion) and the frame-dragging precession (from the new v_\phi) are, for the polar GP-B orbit (a=7027.4 km):

Effect Substrate sector Predicted GP-B measured
Geodetic (de Sitter) radial ebb (SC1) 6604 mas/yr 6601.8\pm18.3
Frame-dragging (Lense–Thirring) azimuthal v_\phi 41 mas/yr 37.2\pm7.2

The geodetic value matches GR’s 6606.1 mas/yr to \sim0.03\%; the frame-dragging value (a simple circular orbit-average; the full GP-B model’s 39.2 mas/yr adds eccentricity and Schwarzschild-coupling corrections) sits well within the measured 37.2\pm7.2. The calculation is in scripts/kerr_frame_dragging.py.

Where v_\phi goes supersonic (|v_\phi|=c) the substrate has an acoustic ergosphere — the rotating analog of Unruh’s “dumb hole,” inside which no static substrate parcel can exist. What remains for a full, non-linearized Kerr match (the standard the Schwarzschild sector already meets) is the combined radial-plus-azimuthal acoustic metric and the explicit boundary-layer entrainment fixing the amplitude from first principles; see open problems WIP-24.

The Cosmological Constant Problem — and Its Resolution

The counter-rotating boundary layers between all orbital systems contain energy. This energy is gravitationally invisible in equilibrium because co- and counter-rotating contributions nearly cancel, it doesn’t couple to electromagnetic probes (dark particles only), and it shows up only through its gravitational effects.

In quantum field theory, the vacuum energy is calculated to be \sim 10^{120} times larger than observed — the worst prediction in physics. The substrate framework resolves this through Volovik’s thermodynamic identity.

In a superfluid at zero temperature and complete thermodynamic equilibrium, the Gibbs-Duhem relation gives:

\varepsilon + P = 0

This is the equation of state for dark energy (w = -1). But in equilibrium, \varepsilon and P are both exactly zero — not just their sum. The superfluid self-tunes: any attempt to add vacuum energy changes the density, which changes the chemical potential, which drives flows that relax the energy back to zero. No fine-tuning is needed.

The residual: an order-unity disequilibrium

One Disequilibrium, Two Rulers The observed cosmological constant is one residual disequilibrium of the substrate, rho_Lambda = rho_ref times (delta T over T_c) squared. Measured against the substrate's own density rho_DM, the disequilibrium delta T over T_c is the square root of 5.8 over 2.4, about 1.6 — order unity, no fine-tuning. Measured against the gravitational Planck density rho_Pl, the same residual reads as 3.4e-62, about 10 to the minus 61.5, the famous worst prediction in physics. The two readings differ by exactly the square root of rho_DM over rho_Pl, equal to (m1 over M_Planck) squared, equal to (Planck length over xi) squared, equal to f_cross omega_0 hbar over 4 pi c cubed xi, about 2.9e-62. Small Lambda is the same number as weak gravity. One Disequilibrium, Two Rulers The same residual ρΛ reads as order-unity or as 10−61.5 — depending only on which density you measure it against. The residual disequilibrium ρΛ = ρref · (δT/Tc)² → δT/Tc = √(ρΛref) ρΛ¼ = 2.24 meV ≈ m1c² = 2.07 meV — the dark-energy scale is the dc1 scale ruler 1 ruler 2 Ruler 1 · the substrate's own density ρref = ρDM ≈ 2.4 × 10−27 kg/m³ δT/Tc = √(5.8 / 2.4) ≈ 1.6 O(1) An order-unity disequilibrium — the vacuum sits one step from equilibrium. No fine-tuning at all. Ruler 2 · the gravitational Planck density ρref = ρPl = c5/ℏG² ≈ 5 × 1096 kg/m³ δT/Tc = √(ρΛPl) ≈ 3.4 × 10−62 ≈ 10−61.5 The famous “worst prediction in physics” — but that number lives only against this ruler. The gap between the two rulers is exactly the gravitational hierarchy: √(ρDMPl) = (m1/MPl)² = (ℓPl/ξ)² = fcrossω0ℏ / 4πc³ξ ≈ 2.9 × 10−62 Small Λ is the same number as weak gravity.

The observed \Lambda comes from the universe not being in perfect equilibrium. Cosmic expansion prevents the substrate from fully relaxing, and the residual vacuum energy is quadratic in the departure from equilibrium:

\rho_\Lambda = \rho_\text{ref} \cdot \left(\frac{\delta T}{T_c}\right)^2

The disequilibrium fraction \delta T/T_c this implies depends entirely on which density \rho_\text{ref} the residual is measured against — and that choice, usually left implicit, is the whole story.

Against the substrate’s own density, the disequilibrium is order unity. The substrate’s natural energy density is its ground-state value \rho_\text{ref} \approx n_1 m_1 \approx \rho_\text{DM} (close-packing). Measured against it,

\frac{\delta T}{T_c}\bigg|_\text{substrate} = \sqrt{\frac{\rho_\Lambda}{\rho_\text{DM}}} \approx \sqrt{\frac{5.8\times10^{-27}}{2.4\times10^{-27}}} \approx 1.6 = \mathcal{O}(1).

There is no fine-tuning of the substrate’s state at all — it sits an order-unity fraction away from full relaxation. This is what the dark-energy data measure directly: the local density today is f(0) = 1.25, and the moraine’s harmonic edge z_\text{harm} = -0.25 lies in our future (see Dark Energy and the Crust). We are still inside the downstream wake of the previous bubble; the substrate has not yet re-equilibrated. The nonzero \Lambda today is that un-drained disequilibrium. (This is constraint C7.)

The famous 10^{-61.5} appears only against the gravitational Planck density. The “worst prediction in physics” measures \rho_\Lambda against \rho_\text{Pl} = c^5/\hbar G^2 \approx 5\times10^{96} kg/m³. That reference gives

\frac{\delta T}{T_c}\bigg|_\text{Planck} = \sqrt{\frac{\rho_\Lambda}{\rho_\text{Pl}}} = 3.4\times10^{-62} \approx 10^{-61.5}.

The two readings differ by exactly the ratio of the two Planck scales, \sqrt{\rho_\text{DM}/\rho_\text{Pl}} = (m_1/M_\text{Pl})^2, so the entire “10^{-61.5} fine-tuning” is the order-unity disequilibrium times the substrate-to-gravitational hierarchy factor:

\frac{\delta T}{T_c}\bigg|_\text{Planck} = \mathcal{O}(1)\times\left(\frac{m_1}{M_\text{Pl}}\right)^2 = 2.9\times10^{-62}.

Small \Lambda is the same number as weak gravity. (m_1/M_\text{Pl})^2 is not an independent small number. With M_\text{Pl}^2 = \hbar c/G and the Volovik speed m_1 = \hbar/c\xi (C1) alone, it is simply the squared ratio of the Planck length \ell_\text{Pl} = \sqrt{\hbar G/c^3} to the substrate cell \xi; the induced-gravity relation G = f_\text{cross}\,\omega_0\xi/4\pi (C3) then rewrites that same ratio through the boundary-transit probability:

\frac{\delta T}{T_c}\bigg|_\text{Planck} = \left(\frac{m_1}{M_\text{Pl}}\right)^2 = \left(\frac{\ell_\text{Pl}}{\xi}\right)^2 = \frac{\hbar G}{c^3\xi^2} = \frac{f_\text{cross}\,\omega_0\,\hbar}{4\pi c^3\xi} \approx 2.9\times10^{-62}.

The two middle forms make plain that the number is geometric — the cell sits some 10^{31} Planck lengths across (\ell_\text{Pl}/\xi \approx 1.6\times10^{-31}, just m_1/M_\text{Pl} read as lengths, since \xi and \ell_\text{Pl} are the reduced Compton wavelengths of m_1 and M_\text{Pl}) — while the last shows its source is the same tiny boundary-transit probability f_\text{cross}\approx10^{-15} that makes gravity weak. The cosmological constant and Newton’s G are the same problem: deriving f_\text{cross}/\omega_0 (see Open Problems WIP-15 item 2, WIP-16) would predict both at once — the same collapse the framework found between the bridge 4\pi and the Higgs 8\pi.

This same hierarchy is the one hand-tuned input of the closest published log-EOS dark fluid. Chavanis’s logotropic model carries a single dimensionless constant B = 1/\ln(\rho_P/\rho_\Lambda) = 1/283, which he notes is quantum in origin (B\to0 as \hbar\to0); its argument is identically the substrate’s weak-gravity number, \ln(\rho_P/\rho_\Lambda) = 2\,|\ln(m_1/M_\text{Pl})^2| = 283, so B_\text{Chavanis} = 1/(2\,|\ln(m_1/M_\text{Pl})^2|). The number Chavanis inputs to make his logarithm cosmological is the one the substrate routes through the marginal point — see the logotropic comparison for the full three-way map.

The dark-energy scale is the dc1 scale — and the two lengths differ by a known factor. The close-packing relation \rho_\text{DM}c^2 = (m_1c^2)^4/(\hbar c)^3 fixes the substrate’s infrared quantum from the matter density alone, m_1c^2 = \rho_\text{DM}^{1/4} = 1.77\ \text{meV}, whose Compton length is the lattice cell \xi_\text{DM} = \hbar c/\rho_\text{DM}^{1/4} = 111.8\ \mu\text{m} — Route 1 of the bridge equation. The dark-energy density defines its own length by the same fourth root,

\xi_\text{DE} = \frac{\hbar c}{\rho_\Lambda^{1/4}}, \qquad \rho_\Lambda^{1/4} = 2.24\ \text{meV}, \qquad \xi_\text{DE} = 88.1\ \mu\text{m}.

This length is not our invention: \xi_\text{DE} = (\hbar c/\rho_\Lambda)^{1/4} \approx 85\ \mu\text{m} is the canonical dark-energy length (Beane 1997; Kapner & Adelberger 2007 title their sub-millimetre-gravity paper exactly “Tests of the Gravitational Inverse-Square Law Below the Dark-Energy Length Scale”), and it is why Eöt-Wash torsion-balance experiments probe \sim100\ \mu\text{m} at all. The two lengths are not equal, and their ratio is not free:

\frac{\xi_\text{DM}}{\xi_\text{DE}} = \frac{\rho_\Lambda^{1/4}}{\rho_\text{DM}^{1/4}} = \left(\frac{\Omega_\Lambda}{\Omega_\text{DM}}\right)^{1/4} = 1.269.

So the long-noted “\rho_\Lambda^{1/4}\sim meV \sim m_1c^2” coincidence, sharpened, is an identity carrying one dark-sector number: the lattice cell exceeds the dark-energy length by exactly (\Omega_\Lambda/\Omega_\text{DM})^{1/4}. The substrate has a single energy-density scale, and dark matter and dark energy are two fourth roots of it — the coincidence problem (\rho_\Lambda\sim\rho_\text{DM} today) dissolves into one density, and the residual 27\% between “the medium” and “its vacuum energy” is the dark-sector ratio. This is the framework’s honest replacement for the broken electroweak cube root (C-08/C-09): not a second, independent determination of a length, but a single dimensionless identity whose measured side, \Omega_\Lambda/\Omega_\text{DM}, comes from cosmology.

One number, three hats — and a caution against counting it thrice. That same ratio \Omega_\Lambda/\Omega_\text{DM} = 2.59 is the only dark-sector datum in what follows, and it reappears at three powers. The moraine crest we sit on (Dark Energy and the Crust) is its fourth root, f(0) = (\Omega_\Lambda/\Omega_\text{DM})^{1/4} = 1.27; the substrate-referenced disequilibrium above is its square root, \delta T/T_c = \sqrt{\rho_\Lambda/\rho_\text{DM}} = 1.61; and the flatness-normalized Volovik base density is \Omega_{\Lambda,\text{base}} = \Omega_\Lambda^{3/4}\Omega_\text{DM}^{1/4} = 0.540 (the DESI self-consistent fit quotes 0.548, using its fitted f(0)=1.25 rather than the predicted 1.27; WIP-18). These obey f(0)^2 = \delta T/T_c and \Omega_{\Lambda,\text{base}} = \Omega_\Lambda/f(0) identically — one quantity wearing three hats, not three independent confirmations, and the framework must not present them as such. Two honesty notes follow. First, the DESI moraine fit constrains only the shape f(z)/f(0) of the crust — the flatness-normalized expansion law E^2(z) = \Omega_m(1+z)^3 + \Omega_\text{tot}\,f(z)/f(0) depends on the ratio, never on the absolute crest f(0) — so the crest’s agreement with 1.27 is a consistency check, not an independent second read of the ratio. Second, this base density (0.548) is not the close-packing equilibrium (\rho_\Lambda = \rho_\text{DM}, \Omega = 0.26): the base is a weighted geometric mean tilted toward today’s \Omega_\Lambda, and the two must not be conflated.

The Planck-referenced number is a frozen photon mass. The split above is bookkeeping; underneath it is a mechanism, and the substrate’s logarithmic equation of state supplies it. The photon is the substrate’s sound mode (see Emergent Speed of Light); the logarithm makes that mode exactly massless only when the vacuum sits at its marginal, perfectly Lorentz-invariant density — the \mu\to0 critical point that, on the log’s Bogoliubov spectrum, is the literal edge of stability the substrate self-tunes toward (massless light \Leftrightarrow marginal stability; see Substrate Particles § The Marginal Point, grounded in § The Logarithmic Equation of State). The approach slows critically: the rate that relaxes the residual super-criticality (the vacuum sitting just above marginality) vanishes at the critical point, so cosmic expansion outruns it and freezes a small residual in place. Read as a photon mass, that frozen residual is the Hubble energy itself,

m_\gamma c^2 \;\sim\; \hbar H_0 \;\approx\; 1.4\times10^{-33}\ \text{eV}, \qquad \lambda_\text{Compton} \;\sim\; \frac{c}{H_0} = R_\text{Hubble}

— a Compton wavelength of order the horizon, the smallest photon mass the observable universe can resolve, and some 10^{15} below the observational bound m_\gamma\lesssim10^{-18}\,\text{eV}. Expressing this rate-scale residual as an energy density is exactly what the Friedmann equation H_0^2 = (8\pi G/3)\,\rho_\text{tot} does, and carrying it through returns the Planck-referenced disequilibrium (m_1/M_\text{Pl})^2, the rate-to-density conversion factor being the dark-energy fraction 8\pi/3\Omega_\Lambda — Friedmann itself. So the famous 10^{-61.5} is not merely a choice of reference frame: it is the frozen super-criticality of a relaxing vacuum, and the celebrated hierarchy is cosmology’s own dictionary between a rate and a density.

The mechanism selects the history, not the number. The relaxation ODE above is a tracker: written out, \dot{(\delta T)} = -\delta T/\tau + \alpha H with \rho_\Lambda = \rho_\ast(\delta T)^2 and critical slowing \tau\to\infty as \delta T\to0, coupled to Friedmann, has been integrated (Open Problems WIP-16). Two results are robust and one negative is load-bearing. It reproduces the correct dark-energy history — \rho_\Lambda\to0 at high redshift, rising through \Omega_\Lambda\sim0.7 today and freezing toward de Sitter — which is exactly the behaviour a naive \rho_\Lambda\propto H^2 tracker gets wrong (that one holds \Omega_\Lambda constant and fails at high z); critical slowing is what freezes the residual late instead. And it is a genuine attractor in initial conditions: the starting \delta T is forgotten, so “now” is not a fine-tuned instant but a generic point on a universal freeze-out. But the frozen value \rho_\Lambda/\rho_m today is set one-to-one by the drive strength (the inherited relaxation depth) — no drive law tested, self-limiting or not, selects it. So the mechanism explains why \rho_\Lambda\sim\rho_\text{DM} is natural and late without predicting the ratio 2.59. This is a sharper competitive position than the closest published twin, Chavanis’s logotropic dark fluid, which fixes the analogous number (\Omega_\text{de}/\Omega_\text{dm}=e) only through an admitted “cosmic coincidence”: the substrate trades that coincidence for a mechanism (attractor + correct history) at the cost of the precise value. Chavanis himself suggests his e-relation “may correspond to a fixed point in a more sophisticated theory” — the substrate’s tracker attractor is a candidate for exactly that (full comparison).

What remains open. Consistent with the tracker, the value of \Lambda — equivalently the de Sitter horizon entropy S_\text{dS} = 1/(\delta T/T_c)^2 \approx 2\times10^{122}, or the bubble’s size at nucleation — is an inherited initial condition: the residual relaxation depth of the previous cycle, the same status as the crust amplitude B. Volovik self-tuning plus the horizon-entropy fluctuation law \delta T/T_c = 1/\sqrt{S_\text{dS}} make \rho_\Lambda = \rho_\text{Pl}/S_\text{dS} a marginal, scale-free self-consistency — it fixes the form of the relation but not the value, exactly what one expects at the substrate’s critical (\mu\to0) point. So the framework predicts the dark-energy scale (m_1) and reduces the apparent tuning to the gravitational hierarchy; the value is set by the nucleation dynamics of \mathcal{B}^{-1} (see A Universe That Boils breadcrumb 1, and Open Problems WIP-17).

Gravity’s Place in the Framework

Gravity, the quantum potential, and photon propagation all arise from the same boundary physics operating at different scales. The quantum potential (Two Fluids → Quantum Potential) is the reaction force of counter-rotating layers on co-rotating flow. Gravity is the macroscopic leak current through those same layers. And photons are modons — counter-rotating vortex dipoles ejected when boundaries reorganize — which is the subject of the next chapter.