The Resonant Theory: A Hydrodynamic Framework for the Quantum Plenum (v10.5)

Author: Justin Time Miller
Version: Version 10.5 (October 2026 — Continuous Non-Linear Acoustic Pressure PDE, Two-Fluid Ram-Impedance Derivation & Soliton Boundary Closure)
Date: October 2026

Abstract

Contemporary theoretical physics relies on non-baryonic dark-sector parameterizations—dark matter and dark energy—to reconcile macroscopic cosmological observations with the Standard Model. This paper formulates an axiomatic, publication-grade physical framework that replaces the empty inviscid vacuum with a continuous, compressible, two-fluid hydrodynamic substrate termed the Quantum Plenum, modeled on the Tisza–Landau two-fluid formalism of superfluid $^{3}\text{He}$ and $^{4}\text{He}$ [1, 2, 5]. From the electromagnetic fine-structure constant $\alpha$ and natural logarithmic attenuation $e$, we define two strictly dimensionless hydrodynamic coupling invariants: the Universal Acoustic Susceptibility $A_{\text{RT}} = \alpha^{-1}/e = 50.412727$ and its reciprocal damping ratio, the Acoustic Venting Coefficient $\beta_{\text{TRT}} = e / A_{\text{RT}} = 0.05392055$, which defines the equilibrium normal-fluid mass fraction $\chi_n \equiv \rho_{n,0}/\rho_0 = \beta_{\text{TRT}}$. Conditioning on six explicitly declared constitutive postulates (P1–P6) and discrete structural indices (Table 1), and anchoring to the local present-epoch acoustic wave speed $c_0 = 2.99792458 \times 10^8\text{ m s}^{-1}$ and carrier frequency $f_0 = 7.64496 \times 10^7\text{ Hz}$ ($L_0 = c_0 / f_0 = 3.92144\text{ m}$), we establish a dimensionally homogeneous system of characteristic Plenum scales (Table 2) for kinematic viscosity ($\nu_0 = c_0 L_0 / A_{\text{RT}} = 2.3320 \times 10^7\text{ m}^2\text{s}^{-1}$), angular relaxation time ($\tau_0 = (2\pi f_0)^{-1} = 2.0818 \times 10^{-9}\text{ s}$), and minimum cavitation energy ($E_{\text{cav}} = A_{\text{RT}} h f_0 = 2.5537 \times 10^{-24}\text{ J}$).

To ensure complete methodological and claim-calibration transparency, all quantitative results are categorized up front into four distinct epistemological tiers (Table 4): the framework offers one prospective, out-of-sample falsification prediction—the 2028–2030 Europa Clipper outer-cruise Doppler residual ($\bar{a}_{\text{plenum}} = (8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$, derived directly from two-fluid momentum flux across an axisymmetric 3D dipole wake and falsifiable if $|\mathbf{a}_{\text{res}}| < 2.0 \times 10^{-10}\text{ m s}^{-2}$ under an exhaustive four-regime decision protocol)—while the subatomic mass ratio, apparent Hubble shift, continuous non-linear galactic pressure PDE, solar tachocline torque balance ($45.9855\text{ laps}$, $M_{\text{Oort}} = 13.98 \pm 1.18\,M_\oplus$), and lunar recession are retrodictive consistency benchmarks conditioned on explicitly declared constitutive postulates and discrete structural integers (Table 1).

Conditioned on six explicitly declared constitutive postulates (P1–P6) and discrete structural indices (Table 1), the framework yields:

1. Lepton–Baryon Ground-State Soliton Rest-Energy Functional & Screened Bjerknes Nuclear Binding: Under the declared ground-state topological soliton energy functional $\mathcal{H}_{\text{sol}}(p, q, D) \equiv E_e \left[ (p \cdot q)\pi^D + \Delta M_{\text{visc}}(p, \beta_{\text{TRT}}) \right]$ on the 5-torus $\mathbb{T}^5$ (Postulate P3) with pairwise poloidal strand shear increment $\Delta M_{\text{visc}}(p, \beta_{\text{TRT}}) \equiv \binom{p}{2}\frac{2}{\pi}\beta_{\text{TRT}} = \frac{p(p-1)}{\pi}\beta_{\text{TRT}}$, the unknotted electron reference state is the $(p_e, q_e, D) = (1, 1, 0)$ ground state with zero internal standing-wave modes ($\pi^0 \equiv 1$, $\binom{1}{2}=0 \implies \Delta M_{\text{visc}} = 0$, giving $\mathcal{H}_{\text{sol}}(1,1,0) = E_e = m_e c_0^2$), while the proton is the knotted $(p, q, D) = (2, 3, 5)$ Trefoil knot $3_1$ whose torus covering map $\Pi_{2,3}: \mathbb{T}^2 \to \mathbb{T}^2$ has Jacobian degree $d_{\text{cover}} = p \cdot q = 6$, sweeping across $D = 5$ dimensionless radian-normalized half-wave standing modes ($\mathcal{V}_5 = (\pi\text{ rad}/1\text{ rad})^5 = \pi^5 = 306.019685$) and coupling across its single unordered poloidal strand pair ($\binom{2}{2}=1 \implies \Delta M_{\text{visc}} = \frac{2}{\pi}\beta_{\text{TRT}} = 0.034327$). The resulting constitutive proton-to-electron rest-energy ratio $\mu_{\text{TRT}} \equiv \mathcal{H}_{\text{sol}}(2,3,5)/\mathcal{H}_{\text{sol}}(1,1,0) = 6\pi^5 + \frac{2}{\pi}\beta_{\text{TRT}} = 1836.152436$ evaluates the CODATA 2022 value $1836.15267343(11)$ to $0.13\text{ ppm}$ ($|\Delta\mu| = 2.37 \times 10^{-4}$). Phase-locked breathing modes between pulsating solitons obey a Yukawa-screened secondary Bjerknes potential whose exponential factor $\exp(-r/\lambda_{\text{screen}})$ arises from evanescent mode decay over the finite carrier relaxation length $\lambda_{\text{screen}} = c_0/\omega_0 = L_0 / (2\pi) = 0.62412\text{ m}$.
2. Mach-1 Sonic-Horizon Characteristic Degeneracy & Cavitation Mass Floor: Under the steady compressible potential-flow equation $(\delta_{ij} - u_i u_j / c_0^2)\partial_i u_j = 0$, the longitudinal principal symbol $1 - |\mathbf{u}|^2/c_0^2 = 1 - M^2$ vanishes at the Mach-1 sonic horizon $|\mathbf{u}| = c_0$ ($\rho_{\text{sonic}} = \rho_0 e^{-1/2} = 0.60653\,\rho_0$, $\Delta P_{\text{sonic}} = (1 - e^{-1/2})P_0 = 0.39347\,P_0$), exhibiting parabolic characteristic degeneracy at the Mach-1 sonic horizon ($|\mathbf{u}| = c_0$, with hyperbolic flow for $|\mathbf{u}| > c_0$). Conditioned on the declared core nucleation energy scale $E_{\text{cav}} \equiv A_{\text{RT}} h f_0 = 2.55373 \times 10^{-24}\text{ J}$ (Postulate P4), this defines the characteristic topological core mass floor $\mu_0 = E_{\text{cav}} / c_0^2 = 2.84140 \times 10^{-41}\text{ kg}$ without vacuum density singularity ($\rho \to 0$), while sub-threshold perturbations disperse as massless acoustic radiation ($m = 0$).
3. Zero-Mean Solenoidal Regularized 3D Vorticity PDE: Introducing the zero-mean solenoidal Leray–Helmholtz projected cubic venting term $-\beta_{\text{TRT}}\tau_0 \mathbb{P}_0(|\boldsymbol{\omega}|^2\boldsymbol{\omega})$ into the 3D vorticity equation preserves $\nabla \cdot \boldsymbol{\omega} = 0$ and the zero spatial mean $\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$ identically on both $\mathbb{R}^3$ and $\mathbb{T}^3$, establishing an a priori finite-horizon Grönwall enstrophy bound on $\mathbb{R}^3$ and time-uniform absorbing-ball bound on $\mathbb{T}^3$ for smooth solutions under Navier–Stokes with cubic Hall–Vinen damping (Theorem 1).
4. Horizon Quarter-Wave Attenuation, CMB Acoustic Peak Separation & Solar Tachocline Balance: Formulating photon propagation along a forward line-of-sight ray $s \in [0, d_L]$ from emitter to observer in the 1D quarter-wave causal horizon resonator ($k_1 R_H = \pi/2$) under the declared constitutive cavity action-transfer coupling $\gamma_{\text{cav}} = k_1 \beta_{\text{TRT}} = \frac{\pi}{2 R_H}\beta_{\text{TRT}}$ (Postulate P5) evaluates the strictly positive $+8.47\%$ apparent local Hubble shift ($H_{\text{late, TRT}} = H_{\text{early}}(1 + \frac{\pi}{2}\beta_{\text{TRT}}) = 73.11\text{ km s}^{-1}\text{Mpc}^{-1}$ from $H_{\text{early}} = 67.40\text{ km s}^{-1}\text{Mpc}^{-1}$) and geometric-mean BAO scale $\lambda_{\text{BAO}} = 152.53\text{ Mpc}$, while explicitly separating the transverse comoving sound horizon $r_s(z_*)$ frozen into the pre-recombination plasma ($\ell_1 = 220.04$) from post-recombination radial line-of-sight carrier attenuation. In the heliosphere, exact solar differential rotation over the $11.0\text{-yr}$ Schwabe cycle accumulates $N_{\text{laps}} = 4017.75\text{ d} \times (24.4754^{-1} - 34.000^{-1})\text{ d}^{-1} = 45.9855\text{ laps}$ ($\Delta\Phi_{\text{Sun}} = 16,554.8^\circ$), yielding an outer Oort boundary mass balance of $M_{\text{Oort}} = 13.98 \pm 1.18\,M_\oplus$.
5. Static Hydrostatic Realization of Newtonian Gravity, Positive Soliton Hamiltonian & Universal Free-Fall: In the static background Plenum ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$, identically satisfying mass continuity $\partial_t\rho + \nabla\cdot(\rho\mathbf{u}) = 0$ everywhere in the exterior vacuum), the continuous carrier radiation reaction body force $\mathbf{f}_{\text{rad}} = +\nabla P$ maintains hydrostatic equilibrium in the unknotted ambient continuum ($\rho_0 \partial_t \mathbf{u} = -\nabla P + \mathbf{f}_{\text{rad}} = \mathbf{0}$). By contrast, a topological matter soliton is a localized, phase-locked vortex singularity whose closed phase boundary scatters incident carrier phonons asymmetrically and whose integrated two-fluid Hamiltonian energy $m_i c_0^2 = \int_{\mathcal{V}_{\text{sol}}} (\frac{1}{2}\rho_s |\mathbf{u}_s|^2 + \frac{1}{2}\rho_n |\mathbf{u}_n|^2 + \Delta\mathcal{U}_{\text{pot}})\,d^3x > 0$ is strictly positive because its quantized toroidal circulation kinetic energy vastly exceeds the shallow core cavitation deficit. Consequently, solitons experience the inward wave-gradient Bernoulli force $\mathbf{F}_{\text{grad}} = -(m_i/\rho_0)\nabla P$, yielding universal inward free-fall acceleration $\mathbf{a} = -(1/\rho_0)\nabla P = -GM/r^2\hat{\mathbf{r}}$ strictly independent of body composition or density (Weak Equivalence Principle).
6. Continuous Non-Linear Acoustic Pressure PDE & Emergent Galactic Rotation Plateau: Both solar-system and galactic gravity are governed by a single continuous non-linear acoustic pressure PDE, $\nabla \cdot [\mu(|\nabla P|/(\rho_0 a_{\text{gal}}))\nabla P] = 4\pi G \rho_0 \rho_M$ with $\mu(x) = x/\sqrt{1+x^2}$ (Postulate P6), where the cosmological acoustic floor $a_0 \equiv c_0 H_{\text{early}}\beta_{\text{TRT}} = 3.531 \times 10^{-11}\text{ m s}^{-2}$ scaled by the circular disc boundary projection factor $\mathcal{C}_{\text{disc}} = \pi$ sets the galactic shear threshold $a_{\text{gal}} \equiv \pi a_0 = 1.109 \times 10^{-10}\text{ m s}^{-2}$ (within $7.6\%$ of the empirical MOND scale). For $x \gg 1$, $\mu(x) \to 1$ recovers the linear Poisson law $\nabla^2 P = 4\pi G \rho_0 \rho_M$ and Newtonian acceleration $a = -GM/r^2$; for $x \ll 1$, Gauss divergence integration yields $|\nabla P|/\rho_0 = \sqrt{GM a_{\text{gal}}}/r$ and flat rotation curves $v_{\text{rot}}^4 = GM a_{\text{gal}}$ directly from the PDE without ad-hoc algebraic interpolation.
7. Two-Fluid Momentum-Flux Derivation of Spacecraft Ram-Impedance & Exhaustive Decision Protocol: From the background metric cavity expansion acceleration $a_H = c_0 H_{\text{early}}$, the single-cycle acoustic impedance $Z_{\text{cyc}} = 2\pi + \beta_{\text{TRT}}$, and the 3D dipole wake meridian projection $A_{\text{eff}} = A_{\text{RT}}/\pi^2$, two-fluid momentum flux derives the inner-heliospheric hyperbolic cruise deceleration $\mathbf{a}_{\text{plenum}}^{(0)} = -a_H (Z_{\text{cyc}}/A_{\text{eff}})\hat{\mathbf{r}} = -c_0 H_{\text{early}}\frac{2\pi+\beta_{\text{TRT}}}{A_{\text{RT}}/\pi^2}\hat{\mathbf{r}} = -8.124 \times 10^{-10}\text{ m s}^{-2}\hat{\mathbf{r}}$ and outer-cruise midpoint $\bar{a}_{\text{plenum}} = (8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$. Table 4 compiles the master epistemological ledger and the 100% mathematically exhaustive four-regime decision protocol for the 2028–2030 Europa Clipper Doppler test.

1. Introduction & Literature Positioning

1.1 Physical Picture: Space as an Acoustic Two-Fluid Medium

To grasp the core physical intuition of The Resonant Theory (TRT) before entering its formal differential equations, consider the difference between an empty geometric stage and a submerged acoustic chamber. In standard cosmology and quantum theory, empty space is treated as a frictionless void; when galaxies rotate faster than visible matter allows, or when distant supernovae appear dimmer than early-universe expansion predicts, invisible "dark matter" particles and repulsive "dark energy" fields must be added to balance the equations. TRT investigates a mechanical alternative: what we call "empty space" is modeled as a continuous, compressible two-fluid medium—the Quantum Plenum—analogous to laboratory superfluid helium ($^{3}\text{He}$ and $^{4}\text{He}$). At equilibrium, $94.61\%$ of this medium acts as an inviscid superfluid condensate ($\rho_{s,0} / \rho_0 = 1 - \beta_{\text{TRT}}$) that transmits electromagnetic and gravitational waves without shear resistance, while the remaining $5.39\%$ ($\chi_n = \rho_{n,0}/\rho_0 = \beta_{\text{TRT}} = 0.05392055$) acts as a viscous "normal fluid" boundary component. Just as a vortex ring in liquid helium is a self-sustaining topological structure of the fluid itself, elementary particles in TRT are modeled as quantized toroidal vortex solitons in the Plenum, while macroscopic orbital and cosmological residuals emerge from the $5.39\%$ normal-fluid attenuation and standing-wave boundary geometry of the medium.

1.2 The Cosmological Crisis and Theoretical Lineage

The persistent tension between quantum field theory and general relativity has required auxiliary dark-sector parameterizations across cosmological scales. Late-time cosmic acceleration and the $5\sigma$ Hubble tension between early-universe Cosmic Microwave Background (CMB) inferences ($H_{\text{early}} = 67.40 \pm 0.50\text{ km s}^{-1}\text{Mpc}^{-1}$ [8]) and local Type Ia supernova distance-ladder measurements ($H_{\text{late}} = 73.04 \pm 1.04\text{ km s}^{-1}\text{Mpc}^{-1}$ [9]) are conventionally attributed to dark energy ($\Lambda$), while flat galactic rotation curves [7] and outer heliospheric angular-momentum residuals are attributed to cold dark matter halos.

Rather than introducing new particle species, The Resonant Theory (TRT) investigates whether these anomalies can be organized within a single constitutive hypothesis: that the physical vacuum acts as a compressible, two-fluid electro-hydrodynamic substrate. This inquiry synthesizes three established physics lineages while introducing an explicit normal-fluid dissipation closure:

1. Superfluid Vacuum & Analogue Gravity Theories: Following Tisza (1938) [1] and Landau (1941) [2], Volovik (2003) [5] showed that quasiparticles in superfluid $^{3}\text{He-A}$ reproduce emergent Lorentz invariance, gauge fields, and effective acoustic metrics [4]. However, standard superfluid-vacuum models treat the vacuum ground state as strictly inviscid ($\nu = 0$) below the Landau critical velocity, providing no quantitative mechanism for secular orbital drag or cosmological redshift attenuation.
2. Topological Soliton Models of Baryons: Following Kelvin's vortex-atom hypothesis and Skyrme's (1961) nonlinear meson theory [3], Faddeev and Niemi (1997) [6] demonstrated that knotted toroidal vortex solitons—specifically the $(2,3)$ Trefoil knot $3_1$—emerge as stable energy minima in three-dimensional field theories, though standard Skyrmion models require tunable coupling parameters to match empirical baryon-to-lepton mass ratios.
3. Galactic Phenomenology Scope & Separation of Program Papers: While this foundational paper focuses on subatomic rest energies, cosmological horizon attenuation, and heliospheric drag, Section 7 formulates the continuous non-linear acoustic pressure PDE $\nabla \cdot [\mu(|\nabla P|/(\rho_0 a_{\text{gal}}))\nabla P] = 4\pi G \rho_0 \rho_M$ connecting solar-system Newtonian gravity to the background acoustic acceleration floor $a_0 \equiv c_0 H_{\text{early}}\beta_{\text{TRT}} = 3.531 \times 10^{-11}\text{ m s}^{-2}$ and effective galactic disc boundary threshold $a_{\text{gal}} \equiv \pi a_0 = 1.109 \times 10^{-10}\text{ m s}^{-2} \approx 1.2 \times 10^{-10}\text{ m s}^{-2}$ underlying flat rotation curves; full parameter-by-parameter fitting across the 175-galaxy SPARC database [7] is formally designated for the companion paper (Paper 2) in this research series.
4. Methodological Transparency & Separation of Postulates from Predictions: TRT models the vacuum as a viscous-normal plus inviscid-superfluid two-component substrate governed by two dimensionless electro-hydrodynamic invariants, $A_{\text{RT}} = \alpha^{-1}/e = 50.412727$ and $\beta_{\text{TRT}} = e / A_{\text{RT}} = 0.05392055$. To maintain strict distinction between constitutive axioms, discrete topological/harmonic integers, retrodictive benchmarks, and prospective out-of-sample predictions, Section 2 declares all six foundational postulates (P1–P6) in Table 1, and Section 8 classifies every quantitative claim in Table 4.

---

2. Foundational Postulates, Two-Fluid Constitutive Equations & Closure Relations for $\beta_{\text{TRT}}$

Physical Picture: Why does a single dimensionless number, $\beta_{\text{TRT}} = 0.05392055$, appear across the boundary-layer skin of a subatomic proton, the damping of a spinning fluid vortex, and the attenuation of light across cosmological distances? Within the Tisza–Landau two-fluid framework, all three effects represent standard constitutive closures of the equilibrium normal-fluid mass fraction ($\chi_n = 5.392\%$). At the medium's characteristic relaxation resonance ($\omega\tau_0 = 1$), the $5.392\%$ normal-fluid fraction sets the per-radian wave attenuation rate; in a spinning vortex tube, centrifugal phase separation expels that same $5.392\%$ normal-fluid fraction into an annular boundary sheath of leading-order thickness $\delta^*/R \approx \frac{1}{2}\beta_{\text{TRT}} = 2.696\%$ (exact nonlinear root $1 - \sqrt{1-\beta_{\text{TRT}}} = 2.733\%$); and in high-vorticity cores, Hall–Vinen mutual friction proportional to $\chi_n$ vents excess spin into acoustic waves. Proposition 1 formalizes these three constitutive closure relations.

2.1 State Variables and Tisza–Landau Two-Fluid Decomposition

Postulate P1 (Two-Fluid Quantum Plenum): Physical space is occupied by a continuous, compressible two-fluid substrate described at every spacetime point $(\mathbf{x}, t)$ by the Tisza–Landau mass and momentum decomposition [1, 2]: \begin{equation} \rho(\mathbf{x}, t) = \rho_s(\mathbf{x}, t) + \rho_n(\mathbf{x}, t), \qquad \mathbf{j}(\mathbf{x}, t) = \rho_s \mathbf{u}_s + \rho_n \mathbf{u}_n \tag{1} \end{equation} where:

• $\rho_s$ $[\text{kg m}^{-3}]$ and $\mathbf{u}_s$ $[\text{m s}^{-1}]$ denote the density and velocity field of the superfluid condensate component, which is irrotational outside quantized vortex cores ($\nabla \times \mathbf{u}_s = 0$) and has zero intrinsic shear viscosity ($\eta_s = 0$);
• $\rho_n$ $[\text{kg m}^{-3}]$ and $\mathbf{u}_n$ $[\text{m s}^{-1}]$ denote the density and velocity field of the normal viscous component, which carries thermal/acoustic excitations and dynamic shear viscosity $\eta_n = \rho_n \nu_0$.

At the present cosmological epoch ($t = t_0$, redshift $z = 0$), the undisturbed background Plenum has uniform total density $\rho_0 = \rho_{s,0} + \rho_{n,0}$, equilibrium normal-fluid mass fraction $\chi_n \equiv \rho_{n,0} / \rho_0 = \beta_{\text{TRT}} = 0.05392055$, and barotropic equation of state $P(\rho) = c_0^2 \rho$. Linearizing the mass continuity equation $\partial_t \rho + \nabla \cdot (\rho \mathbf{u}) = 0$ and momentum conservation equation $\rho_0 \partial_t \mathbf{u} = -\nabla P + \mathbf{f}_{\text{rad}}$ (where $\mathbf{f}_{\text{rad}}$ denotes the acoustic carrier radiation reaction body-force density) around the unperturbed vacuum state ($\nabla P = \mathbf{0}$, $\mathbf{f}_{\text{rad}} = \mathbf{0}$, $\rho = \rho_0$) yields the standard wave equation $\partial_t^2 \delta\rho - c_0^2 \nabla^2 \delta\rho = 0$, where the present-epoch local longitudinal acoustic wave speed is the invariant vacuum speed of light: \begin{equation} c_0 = \sqrt{\left.\frac{\partial P}{\partial \rho}\right|_{\rho_0}} = 2.99792458 \times 10^8\text{ m s}^{-1} \quad [\text{L T}^{-1}] \tag{2} \end{equation} In the presence of a static baryonic mass distribution, the outward radiation reaction body force $\mathbf{f}_{\text{rad}} = +\nabla P$ balances the inward pressure gradient across the ambient unknotted continuum, maintaining the background substrate in exact static hydrostatic equilibrium ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), as formulated in Section 7.1.

Operational Distinction Between Present-Epoch Local $c_0$ and Cosmological Lookback Attenuation: Throughout this framework, $c_0$ is a fixed, invariant local constant at the present epoch ($z = 0$) used in all atomic, laboratory, and solar-system equations (preserving the constancy of $\alpha$ in Section 3.1). Over cosmological lookback paths ($z > 0$), normal-fluid Stokes–Kirchhoff damping along the line of sight attenuates propagating photon frequency $\omega(s)$ without altering local present-day $c_0$, as formulated in Section 5.1.

2.2 Proposition 1: Constitutive Closure Relations Unifying the Operational Roles of $\beta_{\text{TRT}}$

In TRT, the dimensionless constant $\beta_{\text{TRT}} = 0.05392055$ is defined fundamentally as a single physical state ratio: the equilibrium normal-fluid mass fraction of the Plenum, \begin{equation} \chi_n \equiv \frac{\rho_{n,0}}{\rho_0} \equiv \beta_{\text{TRT}} = 0.05392055 \quad [\text{dimensionless: } 1] \tag{3} \end{equation} Proposition 1 (Two-Fluid Constitutive Closures): Under the standard two-fluid closure hypotheses of (i) unit relaxation quality factor $Q_n(\omega_0\tau_0 = 1) = 1$ at the angular relaxation resonance $\tau_0 \equiv \omega_0^{-1}$, (ii) centrifugal phase expulsion of the normal fluid into a peripheral annular vortex sheath, and (iii) Hall–Vinen cubic mutual friction at high strain rates $\tau_0|\boldsymbol{\omega}| \gtrsim 1$, the equilibrium normal-fluid fraction $\chi_n = \beta_{\text{TRT}}$ in Eq. (3) determines three secondary transport coefficients:

1. Per-Radian Acoustic Energy & Single-Quantum Frequency Attenuation Rate ($\Gamma_{\text{rad}} = \beta_{\text{TRT}}$): Consider a longitudinal wave mode $\mathbf{u}(\phi) = \mathbf{u}_0 e^{i\phi}$ ($\phi = k r - \omega t$) coupled to the Plenum's angular relaxation resonance $\omega_0\tau_0 = 1$. In the two-fluid mixture (1), the superfluid fraction $\rho_{s,0}/\rho_0 = 1 - \chi_n$ oscillates without viscous loss, whereas the normal-fluid fraction $\chi_n = \rho_{n,0}/\rho_0$ undergoes Debye/Stokes–Kirchhoff relaxation damping with unit loss tangent $Q_n^{-1}(\omega_0\tau_0=1) = 1$. Thus the fraction of classical wave energy $E = \frac{1}{2}\rho_0 |\mathbf{u}|^2$ dissipated per radian of phase advance ($d\phi = 1\text{ rad}$) is: $\frac{1}{E}\frac{dE}{d\phi} = -\frac{\rho_{n,0}}{\rho_0} Q_n^{-1}\big|_{\omega_0\tau_0=1} = -\chi_n = -\beta_{\text{TRT}}$

   (Scope & Distinction: Proposition 1(1) governs microscopic longitudinal acoustic waves and single-quantum damping at the relaxation resonance $\omega_0\tau_0 = 1$. For optical electromagnetic wave packets traversing macroscopic cosmological distances ($s \gg L_0$, $\omega \gg \omega_0$), wave-action transfer does not occur through resonant single-quantum Stokes absorption of the photon's own optical phase $\omega t$; rather, under Postulate P5, it couples to the background metric horizon cavity mode via the constitutive quarter-wave cavity action-transfer rate $\gamma_{\text{cav}} = k_1 \beta_{\text{TRT}} = \frac{\pi}{2 R_H}\beta_{\text{TRT}}$, as formulated in Section 5.1.) For a single coherent electromagnetic/acoustic wave packet (photon/phonon) of energy $E_\gamma = \hbar\omega$ propagating along a forward number-conserving ray ($N_\gamma = E / (\hbar\omega) = 1$, $\hbar = \text{const}$), adiabatic action transfer to the normal-fluid background satisfies $\frac{dE_\gamma}{E_\gamma} = \frac{d(\hbar\omega)}{\hbar\omega} = \frac{d\omega}{\omega}$, yielding the exact per-radian frequency attenuation law: $\frac{1}{\omega}\frac{d\omega}{d\phi} = \frac{1}{E_\gamma}\frac{dE_\gamma}{d\phi} = -\beta_{\text{TRT}} \quad \left[\text{rad}^{-1}\right]$

2. Annular Vortex Boundary-Layer Thickness Ratio ($\delta^*/R$ — Exact vs. Leading-Order Linearization): In a cylindrical or toroidal vortex tube of outer cross-sectional radius $R$, centrifugal pressure expels the normal-fluid fraction $\chi_n = \beta_{\text{TRT}}$ into a peripheral annular shear sheath of radial thickness $\delta^*$ surrounding the inviscid superfluid core of radius $R - \delta^*$. Equating the annular cross-sectional area fraction $A_{\text{sheath}}/A_{\text{total}}$ to $\chi_n = \beta_{\text{TRT}}$ yields the quadratic equation: $\frac{A_{\text{sheath}}}{A_{\text{total}}} = \frac{\pi R^2 - \pi (R - \delta^*)^2}{\pi R^2} = 2\left(\frac{\delta^*}{R}\right) - \left(\frac{\delta^*}{R}\right)^2 = \chi_n = \beta_{\text{TRT}}$ Solving this quadratic exactly and comparing with its leading-order thin-sheath linearization ($\delta^*/R \ll 1$) gives: $\left(\frac{\delta^*}{R}\right)_{\text{exact}} = 1 - \sqrt{1 - \beta_{\text{TRT}}} = 0.0273339, \qquad \left(\frac{\delta^*}{R}\right)_{\text{lin}} = \frac{1}{2}\beta_{\text{TRT}} = 0.0269603$ The quadratic truncation residual is $\left(\frac{\delta^*}{R}\right)_{\text{exact}} - \frac{1}{2}\beta_{\text{TRT}} = 3.74 \times 10^{-4}$ ($\mathcal{O}(\beta_{\text{TRT}}^2)$), confirming that the linear thin-sheath closure $\delta^*/R \approx \frac{1}{2}\beta_{\text{TRT}}$ is accurate to $1.37\%$ relative ($0.037\%$ absolute).

3. Cubic Vorticity Mutual-Friction Coefficient ($\beta_{\text{TRT}}\tau_0$): Under the Hall–Vinen mutual-friction closure for a two-fluid vortex core of vorticity $\boldsymbol{\omega}$ $[\text{s}^{-1}]$, coupling between the normal-fluid fraction $\chi_n = \beta_{\text{TRT}}$ and quantized vortex lines transfers rotational energy into longitudinal Kelvin/acoustic modes at a nonlinear rate proportional to $\chi_n \tau_0^{-1} (\tau_0|\boldsymbol{\omega}|)^2 \boldsymbol{\omega} = \beta_{\text{TRT}}\tau_0 |\boldsymbol{\omega}|^2\boldsymbol{\omega}$ $[\text{s}^{-2}]$. Applying the zero-mean solenoidal Leray–Helmholtz projector $\mathbb{P}_0$ yields the curl-compatible damping term $-\beta_{\text{TRT}}\tau_0 \mathbb{P}_0(|\boldsymbol{\omega}|^2\boldsymbol{\omega})$ analyzed in Section 6.1. $\blacksquare$

2.3 Summary of Declared Constitutive Postulates & Discrete Structural Indices

To ensure complete methodological transparency and avoid conflating continuous parameter fitting with discrete topological/harmonic indexing, Table 1 lists all six foundational postulates and discrete structural indices used in this framework.

Table 1: Declared Constitutive Postulates and Discrete Structural Indices of the TRT Framework

---

3. Dimensionless Coupling Constants & Characteristic Plenum Scales

Physical Picture: In engineering and fluid mechanics, dimensionless numbers (like the Reynolds or Mach number) describe how a fluid behaves regardless of scale, whereas dimensional scales (meters, seconds, kilograms) set the physical size of the system. In TRT, the fine-structure constant $\alpha^{-1} \approx 137.036$ measures how many electromagnetic wave cycles the Plenum sustains per natural logarithmic attenuation factor $e \approx 2.718$, giving the dimensionless Acoustic Susceptibility $A_{\text{RT}} = \alpha^{-1}/e = 50.412727$. Once we anchor a single macroscopic carrier frequency $f_0 = 76.4496\text{ MHz}$ (corresponding to a fundamental standing wavelength of $L_0 = 3.92144\text{ meters}$), every dimensional transport property of the vacuum—its cycle period $T_0$, its angular relaxation time $\tau_0 = \omega_0^{-1}$, its kinematic viscosity $\nu_0$, and the minimum energy $E_{\text{cav}}$ required to nucleate a cavitated vortex core—is fixed with exact SI units.

We maintain a strict separation between dimensionless coupling ratios (Section 3.1) and dimensional physical scales carrying SI units of length $[\text{L}]$, time $[\text{T}]$, mass $[\text{M}]$, and kinematic viscosity $[\text{L}^2\text{T}^{-1}]$ (Section 3.2).

3.1 The Two Strictly Dimensionless Coupling Invariants ($A_{\text{RT}}$ and $\beta_{\text{TRT}}$)

Using the CODATA 2022 inverse fine-structure constant $\alpha^{-1} = 4\pi \varepsilon_0 \hbar c_0 / e_c^2 = 137.035999177$ and Euler's number $e = 2.71828182846$, Postulate P2 defines the Universal Acoustic Susceptibility ($A_{\text{RT}}$) as the dimensionless ratio of electromagnetic coupling capacity to single-neper logarithmic attenuation: \begin{equation} A_{\text{RT}} \equiv \frac{\alpha^{-1}}{e} = \frac{137.035999177}{2.71828182846} = 50.412727 \quad [\text{dimensionless: } 1] \tag{4} \end{equation} Its reciprocal damped ratio defines the Acoustic Venting Coefficient ($\beta_{\text{TRT}}$): \begin{equation} \beta_{\text{TRT}} \equiv \frac{e}{A_{\text{RT}}} = \frac{e^2}{\alpha^{-1}} = 0.05392055 \quad [\text{dimensionless: } 1] \tag{5} \end{equation} Both $A_{\text{RT}}$ and $\beta_{\text{TRT}}$ are strictly dimensionless pure numbers everywhere in this work (never carrying degree, radian, or dimensional units).

3.2 Characteristic Dimensional Scales of the Plenum

Heuristic binary-octave scaling down from the cyclic Planck frequency $f_P \equiv \frac{\omega_P}{2\pi} = \frac{1}{2\pi}\sqrt{\frac{c_0^5}{\hbar G}} = 2.95209 \times 10^{42}\text{ Hz}$ across $N_0 = 114.85$ octaves plus the hexagonal close-packed ($\theta_{\text{HCP}} = \pi/6 = 30^\circ$) normal-fluid mass-loading shift $\Delta N_{\text{visc}} = \beta_{\text{TRT}} \cos(\pi/6) = 0.046697\text{ octaves}$ gives: \begin{equation} f_{\text{oct}} \equiv \frac{f_P}{2^{114.85 + \beta_{\text{TRT}}\cos(\pi/6)}} = \frac{2.95209 \times 10^{42}\text{ Hz}}{2^{114.8967}} = 7.63454 \times 10^7\text{ Hz}, \qquad f_0 = 7.64496 \times 10^7\text{ Hz} \quad [\text{T}^{-1}] \tag{6} \end{equation} (where $N_0 = 114.85 = 117.5 - \log_2(2\pi)$ formulates the octave chain directly in cyclic frequency units without any implicit $2\pi$ conversion factor). Because $\omega_P \propto G^{-1/2}$, the $2.2 \times 10^{-5}$ experimental uncertainty in $G$ propagates to $\delta\omega_P / \omega_P = \frac{1}{2}(\delta G / G) = 1.1 \times 10^{-5}$, whereas the relative offset between the heuristic octave estimate $f_{\text{oct}} = 7.63454 \times 10^7\text{ Hz}$ and the dimensional carrier anchor $f_0 = 7.64496 \times 10^7\text{ Hz}$ is $\frac{f_0 - f_{\text{oct}}}{f_{\text{oct}}} = \frac{1.042 \times 10^5\text{ Hz}}{7.63454 \times 10^7\text{ Hz}} = 1.365 \times 10^{-3} \; (0.1365\% \approx 0.14\%)$. This narrow $0.14\%$ residual reflects unmodeled higher-order geometric packing and boundary dispersion between the idealized Planck-scale lattice and macroscopic carrier modes; accordingly, $f_0 = 7.64496 \times 10^7\text{ Hz}$ is treated explicitly as a declared empirical dimensional anchor (Postulate P4), with the binary-octave scaling $f_{\text{oct}}$ providing a heuristic cosmological scale motivation aligned to within $0.14\%$. From $f_0 = 7.64496 \times 10^7\text{ Hz}$, $c_0 = 2.99792458 \times 10^8\text{ m s}^{-1}$, and Planck's constant $h = 6.62607015 \times 10^{-34}\text{ J s}$, the characteristic dimensional scales of the Quantum Plenum follow directly from two-fluid kinetic transport theory:

1. Fundamental Carrier Wavelength ($L_0$): The wavelength of a longitudinal carrier wave propagating at speed $c_0$ with frequency $f_0$ is: \begin{equation} L_0 \equiv \frac{c_0}{f_0} = \frac{2.99792458 \times 10^8\text{ m s}^{-1}}{7.64496 \times 10^7\text{ s}^{-1}} = 3.92144\text{ m} \quad [\text{L}] \tag{7} \end{equation}
2. Fundamental Carrier Period ($T_0$) & Characteristic Angular Relaxation Time ($\tau_0$): The single-cycle period of the carrier wave is: \begin{equation} T_0 \equiv \frac{1}{f_0} = \frac{L_0}{c_0} = 1.30805 \times 10^{-8}\text{ s} \quad [\text{T}] \tag{8} \end{equation} The characteristic angular relaxation timescale associated with carrier angular frequency $\omega_0 \equiv 2\pi f_0 = 4.80347 \times 10^8\text{ rad s}^{-1}$ is: \begin{equation} \tau_0 \equiv \frac{1}{\omega_0} = \frac{1}{2\pi f_0} = \frac{T_0}{2\pi} = 2.08183 \times 10^{-9}\text{ s} \quad [\text{T}] \tag{9} \end{equation} guaranteeing that at the carrier angular frequency, the Debye relaxation resonance condition $\omega_0 \tau_0 \equiv 1$ is satisfied identically.
3. Plenum Kinematic Viscosity ($\nu_0$): By Chapman–Enskog and Landau–Khalatnikov kinetic transport theory, a compressible fluid with sound speed $c_0$ and fundamental wave-packet mean free path $L_0$ possesses bare acoustic diffusivity $D_0 = c_0 L_0$. Dividing $D_0$ by the dimensionless acoustic susceptibility quality factor $A_{\text{RT}}$ defines the equilibrium kinematic viscosity of the Plenum: \begin{equation} \nu_0 \equiv \frac{D_0}{A_{\text{RT}}} = \frac{c_0 , L_0}{A_{\text{RT}}} = \frac{(2.99792458 \times 10^8\text{ m s}^{-1})(3.92144\text{ m})}{50.412727} = 2.33199 \times 10^7\text{ m}^2\text{ s}^{-1} \quad [\text{L}^2\text{T}^{-1}] \tag{10} \end{equation}
4. Hydrodynamic Cavitation Energy Threshold ($E_{\text{cav}}$) and Mass Floor ($\mu_0$): A single carrier phonon carries quantum energy $E_1 = h f_0$. Nucleating a persistent cavitated vortex core against the Plenum's acoustic susceptibility barrier $A_{\text{RT}}$ requires the coherent superposition of $A_{\text{RT}}$ carrier quanta within a single wavelength volume, defining the minimum cavitation energy $E_{\text{cav}}$ and its Einstein–de Broglie mass floor $\mu_0 = E_{\text{cav}}/c_0^2$: \begin{equation} E_{\text{cav}} \equiv A_{\text{RT}} h f_0 = 2.55373 \times 10^{-24}\text{ J}, \qquad \mu_0 \equiv \frac{E_{\text{cav}}}{c_0^2} = \frac{A_{\text{RT}} h f_0}{c_0^2} = 2.84140 \times 10^{-41}\text{ kg} \tag{11} \end{equation}

Table 2: Characteristic Scale System of the Quantum Plenum (Exact Arithmetic Verification)

---

4. Topological Soliton Particles & The Universal Rest-Energy Functional

Physical Picture: Why is a proton $1836.15267$ times heavier than an electron, and why do atomic nuclei destabilize when too many nucleons crowd into a shell? In TRT, elementary fermions are not point particles; they are self-locked vortex solitons in the Plenum. An electron is the minimal unknotted vortex filament ($0_1$)—the ground state of the medium carrying zero internal spatial-phase standing modes. A proton, by contrast, is the simplest self-locking 3D torus knot: the $(2,3)$ Trefoil knot ($3_1$), winding twice through its core ($p=2$) and three times around its perimeter ($q=3$). Those $p \times q = 2 \times 3 = 6$ sheets sweep through $D = 5$ independent spatial-phase half-wave standing modes ($\pi$ radians each), generating a geometric volume factor of $\pi^5 = 306.0197$. Multiply that by the 6 interwoven sheets: $6 \times \pi^5 = 1836.118$. Coupling to the $5.392\%$ normal-fluid boundary sheath adds the tiny remaining skin friction ($\frac{2}{\pi}\beta_{\text{TRT}} = 0.034327$), yielding $1836.152436$—matching experiment to $0.13\text{ ppm}$. When nucleons pack into clusters (like Carbon-12), their boundary sheaths reduce empty space to $\approx 20\%$; twelve spheres surround one with $>60^\circ$ ($>1\text{ rad}$) between centers, while forcing a $13\text{th}$ sphere violates the 1-radian boundary, triggering acoustic shear cavitation and nuclear fission.

4.1 The Lepton–Baryon Ground-State Soliton Functional $\mathcal{H}_{\text{sol}}(p, q, D)$

Under Postulate P3 (Table 1), the rest energy of the stable charged lepton and baryon ground-state toroidal solitons in the Quantum Plenum is governed by the closed-form constitutive functional parameterized by the coprime torus winding integers $(p, q)$ and the number $D \in \{0, 5\}$ of active internal half-wave standing-wave phase modes: \begin{equation} \mathcal{H}{\text{sol}}(p, q, D) \equiv E_e \left[ (p \cdot q) \cdot \pi^D + \Delta M{\text{visc}}(p, \beta_{\text{TRT}}) \right] \tag{12} \end{equation} where $E_e \equiv m_e c_0^2$ is the canonical rest energy of the electron reference state, and the pairwise poloidal strand boundary-layer shear increment is defined by: \begin{equation} \Delta M_{\text{visc}}(p, \beta_{\text{TRT}}) \equiv \binom{p}{2}\frac{\beta_{\text{TRT}}}{\pi}\int_0^\pi \sin\theta,d\theta = \frac{p(p-1)}{2}\left(\frac{2}{\pi}\beta_{\text{TRT}}\right) = \frac{p(p-1)}{\pi}\beta_{\text{TRT}} \tag{12a} \end{equation} (Scope of Postulate P3: Eq. (12) is formulated explicitly for the two topologically stable charged fermion ground states: the unknotted lepton ground state $(p_e, q_e, D) = (1, 1, 0)$ and the minimal knotted baryonic ground state $(p_p, q_p, D) = (2, 3, 5)$ [3, 6]. Excited unstable resonances such as $\mu^\pm, \tau^\pm, \Delta$ involve non-equilibrium radial/multipole acoustic excitations outside the ground-state scope of this foundational paper.)

The three mathematical factors of Eq. (12) are rigorously defined as follows:

1. Torus Self-Covering Map & Jacobian Degree ($d_{\text{cover}} = p \cdot q$): A $(p, q)$ toroidal vortex phase state induces the canonical product self-covering map $\Pi_{p,q}: \mathbb{T}^2 \to \mathbb{T}^2$ on the surface 2-torus $\mathbb{T}^2 = S^1_\theta \times S^1_\varphi$ defined by: $\Pi_{p,q}(\theta, \varphi) = (p\theta \bmod 2\pi, \; q\varphi \bmod 2\pi), \qquad \det(D\Pi_{p,q}) = \det\!\begin{pmatrix} p & 0 \\ 0 & q \end{pmatrix} = p \cdot q$ Pulling back the normalized area 2-form $\omega_{\mathbb{T}^2} = \frac{1}{(2\pi)^2}d\theta \wedge d\varphi$ yields the exact topological covering degree $d_{\text{cover}} \equiv \deg(\Pi_{p,q}) = \int_{\mathbb{T}^2} \Pi_{p,q}^*\omega_{\mathbb{T}^2} = p \cdot q$.
2. Dimensionless Radian-Normalized Phase-Space Volume Factor ($\mathcal{V}_D = \pi^D$): In an extended acoustic resonator, standing waves oscillate from node to antinode across half-wave phase intervals $\Delta\vartheta_k = \pi\text{ rad}$ ($k = 1, \dots, D$). Normalizing each angular coordinate by the natural radian unit $\vartheta_0 \equiv 1\text{ rad}$ defines the strictly dimensionless $D$-dimensional standing-wave phase volume: \begin{equation} \mathcal{V}D \equiv \int{[0, \pi]^D} \prod_{k=1}^D \frac{d\vartheta_k}{\vartheta_0} = \prod_{k=1}^D \left(\frac{\pi\text{ rad}}{1\text{ rad}}\right) = \pi^D \quad [\text{dimensionless: } 1] \tag{13} \end{equation} For $D = 0$ (zero internal standing-wave modes), the empty product is identically $\mathcal{V}_0 = \pi^0 \equiv 1$.
3. Pairwise Poloidal Strand Boundary-Layer Drag ($\Delta M_{\text{visc}}(p, \beta_{\text{TRT}})$): Along a meridian cross-section $\theta \in [0, \pi]$, a torus knot with poloidal winding number $p$ threads $p$ distinct vortex strands through the core. Mutual normal-fluid boundary-layer shear occurs between every unordered pair of distinct poloidal strands $\{a, b\} \subset \{1, \dots, p\}$, of which there are $\binom{p}{2} = \frac{p(p-1)}{2}$ pairs. Multiplying the number of interacting strand pairs $\binom{p}{2}$ by the single-pair meridian dipole shear integral $\frac{\beta_{\text{TRT}}}{\pi}\int_0^\pi \sin\theta\,d\theta = \frac{2}{\pi}\beta_{\text{TRT}}$ yields Eq. (12a) identically for all $p \ge 1$.

4.2 Uniform Evaluation of the Lepton Ground State and Baryon Mass Ratio

Evaluating Eq. (12) across the lepton and baryon ground states yields the proton-to-electron mass ratio $\mu_{\text{TRT}} \equiv m_p / m_e$:

1. The Electron Reference State as the $(1, 1, 0)$ Lepton Ground State: The electron is the fundamental unknotted vortex loop ($0_1$) with trivial covering degree $\deg(\Pi_{1,1}) = 1 \cdot 1 = 1$, quantized superfluid circulation $\Gamma_e = h / m_e$, zero internal standing-wave modes ($D = 0 \implies \mathcal{V}_0 = \pi^0 \equiv 1$), and only a single poloidal strand ($p=1 \implies \binom{1}{2} = 0$ strand pairs $\implies \Delta M_{\text{visc}}(1, \beta_{\text{TRT}}) = 0$ identically). Evaluating Eq. (12) for $(1, 1, 0)$ gives: \begin{equation} \mathcal{H}_{\text{sol}}(1, 1, 0) = E_e \left[ (1 \cdot 1) \cdot \pi^0 + 0 \right] = E_e = m_e c_0^2 \tag{14} \end{equation}
2. The Proton as the Knotted $(2, 3, 5)$ Baryonic Soliton: The proton is modeled as the minimal non-trivial torus knot—the $(p, q) = (2, 3)$ Trefoil knot ($3_1$) [3, 6], whose torus covering degree is $d_{\text{cover}} = \deg(\Pi_{2,3}) = 2 \times 3 = 6$. Unlike the unknotted electron filament, the proton is an extended 5-torus cavity supporting standing half-wave resonances across all $D = 5$ phase dimensions (toroidal azimuth, poloidal twist, radial core acoustic standing-wave phase, superfluid condensate phase, and normal-fluid shear phase), with dimensionless phase volume $\mathcal{V}_5 = \pi^5 = 306.019685$. Its geometric rest-energy ratio is: \begin{equation} M_{\text{geom}}(2, 3, 5) \equiv \deg(\Pi_{2,3}) \cdot \mathcal{V}_5 = 6\pi^5 = 6(306.019685) = 1836.118109 \tag{15} \end{equation}
3. Pairwise Poloidal Meridian Boundary-Layer Increment ($\Delta M_{\text{visc}}(2, \beta_{\text{TRT}}) = \frac{2}{\pi}\beta_{\text{TRT}}$): For the Trefoil knot ($p=2$), the two poloidal strands form exactly $\binom{2}{2} = \frac{2(1)}{2} = 1$ interacting strand pair across the meridian dipole profile $\theta \in [0, \pi]$. Evaluating Eq. (12a) gives: \begin{equation} \Delta M_{\text{visc}}(2, \beta_{\text{TRT}}) = \binom{2}{2}\frac{\beta_{\text{TRT}}}{\pi}\int_0^\pi \sin\theta,d\theta = 1 \cdot \frac{2}{\pi}\beta_{\text{TRT}} = \frac{2}{\pi}(0.05392055) = 0.0343269 \tag{16} \end{equation}

Summing the geometric and boundary-layer ratios yields the constitutive TRT mass-ratio formula: \begin{equation} \mu_{\text{TRT}} \equiv \frac{\mathcal{H}{\text{sol}}(2, 3, 5)}{\mathcal{H}{\text{sol}}(1, 1, 0)} = 6\pi^5 + \frac{2}{\pi}\beta_{\text{TRT}} = 1836.118109 + 0.034327 = 1836.152436 \tag{17} \end{equation} Precision Calibration vs. CODATA 2022: Comparing $\mu_{\text{TRT}} = 1836.152436$ to the CODATA 2022 experimental value $\mu_{\text{exp}} = 1836.15267343(11)$ gives an absolute residual of $|\mu_{\text{TRT}} - \mu_{\text{exp}}| = 2.37 \times 10^{-4}$, or a relative approximation accuracy of $0.13\text{ ppm}$ ($99.999987\%$). Because this $0.13\text{ ppm}$ residual exceeds the $6 \times 10^{-11}$ CODATA experimental uncertainty, Eq. (17) is classified in Table 4 as a leading-order constitutive soliton benchmark (prior to $\mathcal{O}(\beta_{\text{TRT}}^3)$ toroidal aspect-ratio curvature corrections) rather than an exact identity.

4.3 Secondary Bjerknes Nuclear Binding, Carrier Relaxation Screening ($\lambda_{\text{screen}}$), Void Fraction ($\epsilon_{\text{eff}}$) & The 1-Radian ($57.2958^\circ$) Limit

In a compressible two-fluid medium of equilibrium density $\rho_0$, two spherical solitons pulsating in phase at angular frequency $\omega_0 = 2\pi f_0 = 2\pi c_0 / L_0$ with volume oscillation amplitudes $\dot{V}_1, \dot{V}_2$ interact through their phase-locked acoustic breathing modes. In an idealized inviscid, unattenuated fluid, the classical secondary Bjerknes potential $\Phi_{\text{Bj},0}(r) \propto -1/r$ produces a long-range $1/r^2$ attraction. In the two-fluid Quantum Plenum, however, near-field breathing-mode pressure perturbations couple to the normal-fluid relaxation resonance ($\omega_0 \tau_0 = 1$, Eq. 9); over a single radian relaxation cycle $\tau_0 = \omega_0^{-1}$, evanescent phase-locked modes satisfy the screened Helmholtz dispersion equation $(\nabla^2 - \lambda_{\text{screen}}^{-2})\Phi_{\text{Bj}}(r) = 0$, where the characteristic carrier relaxation screening length is: \begin{equation} \lambda_{\text{screen}} \equiv c_0 \tau_0 = \frac{c_0}{\omega_0} = \frac{L_0}{2\pi} = \frac{3.92144\text{ m}}{2\pi} = 0.62412\text{ m} \quad [\text{L}] \tag{18a} \end{equation} (At the macroscopic carrier scale $L_0$, $\lambda_{\text{screen}} = L_0/(2\pi) = 0.62412\text{ m}$ sets the single-radian carrier screening length, while at the femtometer baryonic core scale $R_{\text{nuc}} = L_0 / N_{\text{sub}}$, the corresponding sub-harmonic vortex breathing mode decays over the femtometer nuclear screening length $\lambda_{\text{nuc}} = \lambda_{\text{screen}}/N_{\text{sub}} \sim 1.4\text{ fm}$.) Multiplying the $1/r^2$ secondary Bjerknes interaction by this Yukawa relaxation envelope $\exp(-r/\lambda_{\text{screen}})$ yields the screened secondary Bjerknes acoustic radiation force: \begin{equation} \mathbf{F}_{\text{Bj}}(r) = -\frac{\rho_0 \langle \dot{V}1 \dot{V}2 \rangle}{4\pi r^2} \exp!\left(-\frac{r}{\lambda{\text{screen}}}\right)\hat{\mathbf{r}}, \qquad \lambda{\text{screen}} \equiv \frac{c_0}{\omega_0} = \frac{L_0}{2\pi} = 0.62412\text{ m} \tag{18} \end{equation} In spherical nuclear coordination geometry, 12 equal nucleon spheres of radius $R$ arranged on the coordination shell ($R_{\text{coord}} = 2R$) around a central nuclear barycentric potential void ($R_{\text{core}} = 0$) form the Carbon-12 ($^{12}\text{C}$) close-packed coordination cluster:

1. Effective Interstitial Void Fraction ($\epsilon_{\text{eff}}$ — Both Linearized and Exact Sheath Roots): By Kepler–Hales close-packing, bare rigid spheres of radius $R$ have maximum FCC/HCP packing density $\phi_{\text{FCC}} = \frac{\pi}{3\sqrt{2}} = 0.740480$ (bare void fraction $1 - \phi_{\text{FCC}} = 0.259520$). Including the normal-fluid boundary-layer sheath derived in Proposition 1 expands the effective hydrodynamic exclusion radius to $R_{\text{eff}} = R(1 + \delta^*/R)$, yielding the Effective Interstitial Void Fraction for both the leading-order linear closure ($\delta^*/R = \frac{1}{2}\beta_{\text{TRT}} = 0.026960$) and the exact quadratic root ($\delta^*/R = 1 - \sqrt{1-\beta_{\text{TRT}}} = 0.027334$): \begin{equation} \epsilon_{\text{eff}}^{\text{(lin)}} = 1 - \phi_{\text{FCC}}\left(1 + \frac{1}{2}\beta_{\text{TRT}}\right)^3 = 0.19800 \approx 0.200, \qquad \epsilon_{\text{eff}}^{\text{(exact)}} = 1 - \phi_{\text{FCC}}\left(2 - \sqrt{1 - \beta_{\text{TRT}}}\right)^3 = 0.19712 \approx 0.200 \tag{19} \end{equation} Both roots agree to within $8.8 \times 10^{-4}$ ($0.4\%$ relative), establishing an effective interstitial void fraction of $\approx 20\%$.
2. Spatial Hard-Core Chord Distance ($d_{\text{chord}}$) vs. Geodesic Surface Wave Path ($s(\theta)$): On a spherical coordination shell of radius $R_{\text{coord}} = 2R$, two distinct separation conditions govern physical stability:
   • Spatial Hard-Core Euclidean Chord Distance ($d_{\text{chord}}$): The straight-line Euclidean distance through 3D space between two sphere centers separated by angle $\theta$ is $d_{\text{chord}}(\theta) = 2 R_{\text{coord}}\sin(\theta/2) = 4R\sin(\theta/2)$. Hard-sphere non-penetration requires $d_{\text{chord}}(\theta) \ge 2R \iff \sin(\theta/2) \ge 1/2 \iff \theta \ge \pi/3 = 60.0000^\circ$, which is satisfied exactly by close-packed FCC kissing configurations.
   • Geodesic Surface Acoustic Boundary-Wave Path ($s(\theta)$ & The 1-Radian Packing Closure): Along the coordination shell surface, the acoustic boundary-layer shear wave propagates along the geodesic arc $s(\theta) = R_{\text{coord}}\,\theta = 2R\,\theta$. Under the declared geodesic surface-wave exclusion closure, independent secondary Bjerknes pressure wells require the geodesic surface arc $s(\theta)$ to be at least one full nucleon exclusion diameter $d_{\text{ex}} = 2R$: $s(\theta) \ge d_{\text{ex}} \iff 2R\,\theta \ge 2R \iff \theta \ge 1\text{ radian} = \frac{180^\circ}{\pi} = 57.2958^\circ$ Because the scale factor $2R$ cancels identically, $\theta \ge 1\text{ rad}$ defines a scale-free geodesic angular lower bound. For $N = 12$ outer spheres, both the Tammes icosahedral maximum $\theta_{\text{Tammes}}(12) = 63.4349^\circ$ and FCC kissing angle $60.0000^\circ$ exceed both $60^\circ$ and $1\text{ rad}$ ($57.2958^\circ$). Forcing a $13\text{th}$ sphere compresses the maximum angular separation to $\theta_{\text{Tammes}}(13) = 57.1367^\circ < 1\text{ rad} < 60^\circ$, violating both the Euclidean hard-core chord bound and the geodesic surface-wave bound and triggering acoustic shear fission.

4.4 Mach-1 Sonic-Horizon Characteristic Degeneracy & Topological Core Mass-Floor Scale

Under the linear barotropic equation of state $P(\rho) = c_0^2 \rho$ (Postulate P1), the static isothermal bulk modulus $\kappa \equiv \rho(\partial P/\partial\rho) = \rho c_0^2 > 0$ is strictly positive at all non-zero densities, so the fluid exhibits no thermodynamic bulk modulus instability ($\partial P/\partial\rho = 0$). Rather, the physical boundary separating smooth, subsonic potential flow from a cavitated vortex core is a kinematic PDE characteristic degeneracy (the Mach-1 sonic horizon) [4, 5].

In a barotropic compressible fluid with ambient boundary conditions $u(\infty) = 0$ and $\rho(\infty) = \rho_0$, steady irrotational flow ($\mathbf{u} = \nabla\Phi_v$) satisfies the compressible Bernoulli invariant: \begin{equation} \frac{1}{2}|\mathbf{u}|^2 + \int_{\rho_0}^\rho \frac{dP}{\rho'} = \frac{1}{2}|\mathbf{u}|^2 + c_0^2 \ln!\left(\frac{\rho}{\rho_0}\right) = 0 \implies \rho(\mathbf{u}) = \rho_0 \exp!\left(-\frac{|\mathbf{u}|^2}{2 c_0^2}\right) \tag{20a} \end{equation} Taking the spatial gradient $\nabla\rho = -\frac{\rho}{c_0^2}(\mathbf{u}\cdot\nabla)\mathbf{u}$ and substituting it into the steady mass continuity equation $\nabla\cdot(\rho\mathbf{u}) = \rho\nabla\cdot\mathbf{u} + \mathbf{u}\cdot\nabla\rho = 0$ yields the exact steady compressible potential-flow PDE: \begin{equation} \left(\delta_{ij} - \frac{u_i u_j}{c_0^2}\right)\partial_i u_j = 0 \tag{20b} \end{equation} Along the streamline direction $\hat{\mathbf{s}} = \mathbf{u}/|\mathbf{u}|$, the longitudinal principal symbol of Eq. (20b) is $\sigma_L(\mathbf{u}) = 1 - |\mathbf{u}|^2/c_0^2 = 1 - M^2$ (where $M \equiv |\mathbf{u}|/c_0$ is the local Mach number). For subsonic flow ($M < 1$), Eq. (20b) is strictly elliptic; at the Mach-1 sonic horizon $|\mathbf{u}| = c_0$ ($M = 1$), $\sigma_L(\mathbf{u})$ vanishes identically ($\sigma_L = 0$), rendering the PDE parabolically degenerate and forming an acoustic sonic horizon [4, 5] at the finite sonic-horizon core density and pressure depression: \begin{equation} \rho_{\text{sonic}} = \rho_0 \exp!\left(-\frac{c_0^2}{2 c_0^2}\right) = \rho_0 e^{-1/2} = 0.60653 , \rho_0 > 0, \qquad \Delta P_{\text{sonic}} \equiv P_0 - P_{\text{sonic}} = \left(1 - e^{-1/2}\right) P_0 = 0.39347 , P_0 \tag{20c} \end{equation} Thus core phase separation occurs at the finite Mach-1 sonic horizon $\rho_{\text{sonic}} = 0.60653\,\rho_0$ without ever requiring a logarithmic vacuum density singularity ($\rho \to 0$). Conditioned on Postulate P4 (Table 1), which declares the characteristic cavitation energy scale $E_{\text{cav}} \equiv A_{\text{RT}} h f_0 = 2.55373 \times 10^{-24}\text{ J}$ as the constitutive activation energy required to sustain a cavitated core at the Mach-1 sonic horizon, the framework classifies localized states into two distinct dynamical regimes:

1. Sub-Threshold Wave Perturbations ($E_k < E_{\text{cav}}$): Perturbations carry energy below the declared constitutive threshold $E_{\text{cav}}$; they remain governed by the linear subsonic potential equation ($M < 1$) and disperse as linear, massless longitudinal acoustic radiation propagating at speed $c_0$ ($m = 0$).
2. Super-Threshold Vortex Nucleation ($E_k \ge E_{\text{cav}}$): For wave packets exceeding the constitutive cavitation energy threshold ($E_k \ge E_{\text{cav}}$), dynamic enstrophy forms a stable, persistent cavitated core boundary at the Mach-1 sonic horizon ($|\mathbf{u}| = c_0$), establishing the invariant characteristic core nucleation mass floor $\mu_0$: \begin{equation} \mu_0 \equiv \frac{E_{\text{cav}}}{c_0^2} = \frac{A_{\text{RT}} h f_0}{c_0^2} = 2.84140 \times 10^{-41}\text{ kg} \tag{20d} \end{equation} For an arbitrary super-threshold excitation of energy $E_k \ge E_{\text{cav}}$, the associated rest mass is parameterized by $m(E_k) \equiv \frac{E_k}{c_0^2}\,\Theta(E_k/E_{\text{cav}} - 1) \ge \mu_0$, with exact dimensional homogeneity ($[E_k / E_{\text{cav}}] = \text{J}/\text{J} = 1$, $[\mu_0] = \text{kg}$). (Epistemological Classification: Under Table 4, $\mu_0$ is classified as a Category A/B declared constitutive core nucleation scale anchored by Postulate P4 at the Mach-1 sonic horizon, rather than an ab initio theorem derived from packet volume integrals under an unmodified linear barotrope.)

---

5. Cosmological Wave-Attenuation Model & Resolution of the Hubble Tension

Physical Picture: Why do measurements of the early universe (via the Cosmic Microwave Background) infer a baseline expansion rate of $H_{\text{early}} = 67.40\text{ km s}^{-1}\text{Mpc}^{-1}$, while local supernova telescopes (SH0ES) measure $H_{\text{late}} = 73.04\text{ km s}^{-1}\text{Mpc}^{-1}$? In an empty vacuum, photon redshift can only arise from metric expansion, so an $+8.4\%$ jump in apparent redshift requires invoking a late-time "dark energy" acceleration. In the two-fluid Plenum, a 1D line-of-sight acoustic waveguide extending from the observer ($r = 0$, pressure antinode $\delta P'(0) = 0$) to the causal horizon ($r = R_H$, pressure node $\delta P(R_H) = 0$) supports the exact fundamental cosine standing mode $\delta P(r) = \delta P_0 \cos(k_1 r)$ whose first zero occurs at $k_1 R_H = \pi/2\text{ radians}$ ($90^\circ$ quarter-wave cavity). Across that $\pi/2\text{ rad}$ horizon phase, the $5.392\%$ normal-fluid component contributes an integrated acoustic attenuation susceptibility $\eta_{\text{TRT}} = \frac{\pi}{2}\beta_{\text{TRT}} = +0.0847$ ($+8.47\%$). When local distance-ladder surveys measure the total redshift slope $dz_{\text{obs}}/dd_L$, they observe the sum of baseline expansion and horizon-coupled acoustic attenuation: $67.40 \times (1 + 0.084698) = 73.11\text{ km s}^{-1}\text{Mpc}^{-1}$—matching the SH0ES measurement to $0.07\sigma$ without dark energy.

5.1 Exact Quarter-Wave Horizon Eigenmode ($k_1 R_H = \pi/2$), Line-of-Sight Ray Propagation, CMB First Acoustic Peak Separation ($\ell_1 = 220.04$) and the $+8.47\%$ Hubble Shift

We now derive the mathematical link between normal-fluid acoustic damping $\beta_{\text{TRT}}$ and the fractional discrepancy between early-universe CMB inferences ($H_{\text{early}} = 67.40 \pm 0.50\text{ km s}^{-1}\text{Mpc}^{-1}$ [8]) and local distance-ladder measurements ($H_{\text{late}} = 73.04 \pm 1.04\text{ km s}^{-1}\text{Mpc}^{-1}$ [9]).

1. Exact 1D Line-of-Sight Helmholtz Boundary-Value Problem ($k_1 R_H = \pi/2$): Along a 1D line-of-sight acoustic waveguide ray $r \in [0, R_H]$ connecting the observer at $r = 0$ to the causal Hubble horizon at $r = R_H = c_0 / H_{\text{early}} = 4,447.96\text{ Mpc}$, longitudinal pressure perturbations $\delta P(r)$ obey the 1D Helmholtz equation: $\frac{d^2\,\delta P(r)}{dr^2} + k^2\,\delta P(r) = 0, \qquad r \in [0, R_H]$ subject to mixed Neumann–Dirichlet quarter-wave resonator boundary conditions: $\text{(i) Observer pressure antinode: } \left.\frac{d\,\delta P}{dr}\right|_{r=0} = 0, \qquad \text{(ii) Causal horizon pressure node: } \delta P(R_H) = 0$ The unique general solution of $\delta P''(r) + k^2\delta P(r) = 0$ satisfying the Neumann condition $\delta P'(0) = -k\,\delta P_0 \sin(0) = 0$ is the exact cosine eigenfunction: $\delta P(r) = \delta P_0 \cos(k r)$ Imposing the Dirichlet horizon condition $\delta P(R_H) = \delta P_0 \cos(k R_H) = 0$ yields the exact eigenvalue spectrum $k_n R_H = \left(n - \frac{1}{2}\right)\pi$ for $n = 1, 2, 3, \dots$, whose fundamental mode ($n = 1$) has its exact first root at $k_1 R_H = \pi/2$: $k_1 R_H = \frac{\pi}{2} \implies \phi(r) \equiv k_1 r = \frac{\pi}{2}\frac{r}{R_H} \in \left[0, \frac{\pi}{2}\right], \qquad \Delta\phi = \phi(R_H) - \phi(0) = \frac{\pi}{2}\text{ rad}$ (Equivalently, in 3D spherical coordinates where $\frac{1}{r^2}\frac{d}{dr}(r^2 \frac{d\psi}{dr}) + k^2\psi = 0$, the canonical radial displacement potential $u(r) \equiv r\,\psi(r)$ satisfies $u''(r) + k^2 u(r) = 0$ with origin regularity $u(0) = 0$, giving $u(r) = u_0 \sin(kr)$; imposing the open-horizon Neumann condition $u'(R_H) = k u_0 \cos(k R_H) = 0$ likewise yields the exact fundamental eigenvalue $k_1 R_H = \pi/2$.)

2. Line-of-Sight Ray Propagation Coordinate & Action-Transfer Coupling: To eliminate any coordinate-direction ambiguity, let $s \in [0, d_L]$ denote the forward physical ray propagation distance measured from the emitting source ($s = 0$, where $\omega(0) = \omega_{\text{emit}}$) to the observer ($s = d_L$, where $\omega(d_L) = \omega_{\text{obs}}$).

   Along the propagation path $ds$, the photon traverses the macroscopic cosmological standing cavity mode $\delta P(r) = \delta P_0 \cos(k_1 r)$. Because the cavity mode carries fundamental spatial wavenumber $k_1 = \frac{\pi}{2 R_H}$, traversing distance $ds$ advances the traversed cavity eigenphase by $d\theta_{\text{cav}} = k_1 ds$ radians. Under Postulate P5 (Cymatic Shell & Waveguide Closures), the normal-fluid fraction $\chi_n = \beta_{\text{TRT}}$ sets the action-transfer coupling between the propagating electromagnetic wave packet and the background metric cavity: wave action is transferred at the rate $\beta_{\text{TRT}}$ per radian of cavity eigenphase traversed ($d\theta_{\text{cav}} = k_1 ds$), yielding the constitutive quarter-wave horizon action-transfer rate $\gamma_{\text{cav}} \equiv \beta_{\text{TRT}} k_1 = \beta_{\text{TRT}}\frac{\pi}{2 R_H}$: \begin{equation} \frac{dE_\gamma}{E_\gamma} = \frac{d\omega}{\omega} = -\beta_{\text{TRT}} , d\theta_{\text{cav}} = -\beta_{\text{TRT}} (k_1 ds) = -\left(\beta_{\text{TRT}}\frac{\pi}{2 R_H}\right) ds \tag{21a} \end{equation} (Epistemological Classification: Under Table 4, Eq. (21a) is classified as a Category B retrodictive benchmark governed by the Postulate P5 macro-waveguide action-transfer closure. It is distinct from the microscopic acoustic resonance damping of Proposition 1(1), ensuring strictly achromatic optical transmission without introducing chromatic Rayleigh scattering or line-broadening.) Dividing by $ds$ yields the strictly constant, achromatic spatial attenuation rate declared under Postulate P5: \begin{equation} \gamma_{\text{cav}} \equiv -\frac{1}{\omega}\frac{d\omega}{ds} = k_1 \beta_{\text{TRT}} = \frac{\pi}{2 R_H}\beta_{\text{TRT}} \tag{21b} \end{equation} Across the entire causal horizon depth $s = R_H$, the integrated horizon attenuation bias is: \begin{equation} \eta_{\text{TRT}} \equiv \int_0^{R_H} \gamma_{\text{cav}},ds = k_1 R_H \beta_{\text{TRT}} = \frac{\pi}{2}\beta_{\text{TRT}} = \frac{\pi}{2}(0.05392055) = +0.084698 \quad (+8.47%) \tag{21} \end{equation} Including both the cosmological metric expansion rate $H_{\text{early}}/c_0$ and the cavity action-transfer attenuation $\gamma_{\text{cav}}$, the photon angular frequency $\omega(s)$ evolves along the forward line of sight according to: \begin{equation} \frac{1}{\omega}\frac{d\omega}{ds} = -\frac{H_{\text{early}}}{c_0} - \gamma_{\text{cav}} = -\frac{H_{\text{early}}}{c_0}\left(1 + \frac{\pi}{2}\beta_{\text{TRT}}\right), \qquad \omega(0) = \omega_{\text{emit}} \tag{21c} \end{equation} Integrating Eq. (21c) from the emitting source ($s = 0$) to the observer at distance $s = d_L$ gives: \begin{equation} \omega(d_L) = \omega_{\text{emit}} \exp!\left[ -\frac{H_{\text{early}} d_L}{c_0}\left(1 + \frac{\pi}{2}\beta_{\text{TRT}}\right) \right] \tag{21d} \end{equation} Defining the observed cosmological redshift $1 + z_{\text{obs}}(d_L) \equiv \frac{\omega_{\text{emit}}}{\omega(d_L)}$ yields: \begin{equation} 1 + z_{\text{obs}}(d_L) = \exp!\left[ \frac{H_{\text{early}} d_L}{c_0}\left(1 + \frac{\pi}{2}\beta_{\text{TRT}}\right) \right] \approx 1 + \frac{H_{\text{early}} d_L}{c_0}\left(1 + \frac{\pi}{2}\beta_{\text{TRT}}\right) + \mathcal{O}(d_L^2) \tag{21e} \end{equation}

3. Resolution of Redshift Law vs. CMB First Acoustic Peak ($\ell_1 = 220.04$) & Exact Evaluation of $H_{\text{late, TRT}}$: A critical physical distinction separates the transverse comoving sound horizon $r_s(z_*)$ imprinted on the Cosmic Microwave Background at photon decoupling ($z_* \approx 1090$) from the radial line-of-sight carrier attenuation accumulated by free-streaming optical photons at late times ($z \ll z_*$):

   • Pre-Recombination Transverse Acoustic Horizon ($r_s(z_*)$ and $\ell_1 = 220.04$): Prior to recombination ($t \le t_*$, $z \ge z_* \approx 1090$), photons are tightly coupled by Thomson scattering to the local relativistic baryon–photon plasma, where acoustic waves propagate locally at the plasma sound speed $c_s = c_0 / \sqrt{3(1 + R_b)}$. At decoupling ($z_* \approx 1090$), the comoving sound horizon $r_s(z_*) = \int_0^{t_*} c_s(t)\,dt/a(t)$ is frozen across the last-scattering surface as a transverse physical ruler. Because a uniform achromatic line-of-sight frequency shift (Eq. 21d) rescales photon energy $\hbar\omega$ without bending transverse ray angles $d\theta$, the observed angular acoustic scale $\theta_* \equiv r_s(z_*) / D_A(z_*)$—and thus the CMB first acoustic peak multipole $\ell_1 = \pi / \theta_* = 220.04$—is governed purely by the quarter-wave causal cavity geometry ($k_1 R_H = \pi/2$) and the true unattenuated background expansion rate $H_{\text{early}} = 67.40 \pm 0.50\text{ km s}^{-1}\text{Mpc}^{-1}$ [8], without double-counting $H_0$ or mixing in the line-of-sight attenuation bias $\eta_{\text{TRT}}$.
   • Post-Recombination Radial Line-of-Sight Estimator ($H_{\text{late, TRT}} = 73.11\text{ km s}^{-1}\text{Mpc}^{-1}$): Conversely, local Type Ia supernova distance-ladder surveys (SH0ES [9]) do not measure transverse early-plasma angles; they measure the radial line-of-sight frequency shift $z_{\text{obs}}(d_L)$ of free-streaming optical photons traversing the post-recombination macroscopic cavity mode $\delta P(r) = \delta P_0 \cos(k_1 r)$ along $s \in [0, d_L]$. Differentiating Eq. (21e) with respect to $d_L$ at $d_L \to 0$ yields the strictly positive local Hubble estimator shift: \begin{equation} H_{\text{late, TRT}} \equiv c_0 \left.\frac{dz_{\text{obs}}}{dd_L}\right|{d_L \to 0} = H{\text{early}} \left(1 + \frac{\pi}{2}\beta_{\text{TRT}}\right) = 67.40 \times 1.084698 = 73.11\text{ km s}^{-1}\text{Mpc}^{-1} \tag{22} \end{equation} This predicted value ($73.11\text{ km s}^{-1}\text{Mpc}^{-1}$) agrees with the SH0ES measurement ($73.04 \pm 1.04\text{ km s}^{-1}\text{Mpc}^{-1}$ [9]) to within $0.07\sigma$ ($0.10\%$ relative error) without dark energy ($\Lambda$) or any modification to present-epoch local $c_0$.

5.2 Cymatic Shell Quantization and the Geometric-Mean BAO Scale ($\lambda_{\text{BAO}} = 152.53\text{ Mpc}$)

Along the propagation path $s \in [0, R_H]$, the exponential action-attenuation factor is $\exp(\gamma_{\text{cav}} s) = \exp(k_1 \beta_{\text{TRT}} s)$. At the causal horizon boundary $s = R_H$ ($\theta_{\text{cav}} = \pi/2$), the total attenuation factor is $\exp(\frac{\pi}{2}\beta_{\text{TRT}}) \approx 1 + \frac{\pi}{2}\beta_{\text{TRT}} = 1 + \eta_{\text{TRT}}$. At the cavity spatial midpoint $s = R_H / 2$ ($\theta_{\text{mid}} = k_1 (R_H/2) = \pi/4$), the cumulative attenuation is $\exp(\frac{1}{2}\eta_{\text{TRT}}) = \sqrt{\exp(\eta_{\text{TRT}})}$. To leading order in $\eta_{\text{TRT}} = \frac{\pi}{2}\beta_{\text{TRT}} = 0.0847$, this midpoint expansion corresponds to the leading-order geometric-mean approximation between the unattenuated early rate $H_{\text{early}}$ and the attenuated late rate $H_{\text{late, TRT}}$ (differing from the exact exponential midpoint $\exp(\eta_{\text{TRT}}/2) = 1.04326$ by the leading-order Taylor difference $\frac{1}{4}\eta_{\text{TRT}}^2 = 0.179\% \approx 0.18\%$): \begin{equation} \bar{H}{\text{cav}} \approx \sqrt{H{\text{early}},H_{\text{late, TRT}}} = H_{\text{early}}\sqrt{1 + \frac{\pi}{2}\beta_{\text{TRT}}} = 67.40 \times \sqrt{1.084698} = 70.196\text{ km s}^{-1}\text{Mpc}^{-1} \tag{23a} \end{equation} which defines the effective cavity length $L_{\text{cav}} \equiv c_0 / \bar{H}_{\text{cav}} = 4,270.77\text{ Mpc}$. Under Postulate P5, standing pressure waves inside $L_{\text{cav}}$ are quantized by the canonical spherical harmonic shell closure sequence: \begin{equation} N_{\text{shell}} \in {2, 8, 20, 28, 50, 82, 126} \tag{23b} \end{equation} corresponding to the Mayer–Jensen 3D spherical well eigenmode sequence with spin-orbit boundary coupling. In this sequence, the 5th shell closure ($N_5 = 50$) corresponds to the integer floor of $A_{\text{RT}} = 50.4127$ (governing the bimodal Earth-Moon waveguide $2N_5 = 100$), and the 4th spherical shell closure ($N_4 = 28$) defines the primary baryonic acoustic peak wavelength: \begin{equation} \lambda_{\text{BAO}} \equiv \frac{L_{\text{cav}}}{N_4} = \frac{c_0}{N_4 H_{\text{early}}\sqrt{1 + \frac{\pi}{2}\beta_{\text{TRT}}}} = \frac{4,270.77\text{ Mpc}}{28} = 152.53\text{ Mpc} \tag{23} \end{equation} (Evaluated at the unattenuated early-universe endpoint $R_H = c_0 / H_{\text{early}} = 4,447.96\text{ Mpc}$ and attenuated late-universe endpoint $L_{\text{late}} = c_0 / H_{\text{late, TRT}} = 4,100.64\text{ Mpc}$, the $N_4 = 28$ harmonic spans $[146.45, 158.86]\text{ Mpc}$ with geometric mean $152.53\text{ Mpc}$ ($\sqrt{146.45 \times 158.86} = 152.53\text{ Mpc}$), matching the empirical BAO comoving sound horizon $r_d \approx 147\text{–}153\text{ Mpc}$.)

5.3 Dimensionally Consistent Solar Tachocline Torque Balance & Calibrated 8-Planet Summation

Over an 11.0-year Schwabe magnetic cycle ($T_{\text{cyc}} = 4,017.75\text{ d}$), solar differential rotation between the equatorial tachocline ($P_{\text{eq}} = 24.4754\text{ d}$) and polar regions ($P_{\text{pole}} = 34.000\text{ d}$) accumulates an exact continuous torsional lap count of: $N_{\text{laps}} = 4,017.75\text{ days} \times \left(24.4754^{-1} - 34.000^{-1}\right)\text{ day}^{-1} = 4,017.75 \times (0.04085735 - 0.02941176) = 45.9855\text{ laps}$ (where $N_{\text{mode}} = 46$ is the nearest integer harmonic mode number, while the exact continuous phase wind-up angle $\Delta\Phi_{\text{Sun}}$ is evaluated directly from $N_{\text{laps}} = 45.9855\text{ laps}$): \begin{equation} \Delta\Phi_{\text{Sun}} = 360^\circ \times T_{\text{cyc}}\left(\frac{1}{P_{\text{eq}}} - \frac{1}{P_{\text{pole}}}\right) = 360^\circ \times 45.9855\text{ laps} = 16,554.8^\circ ; (288.9355\text{ rad}) \tag{24} \end{equation} To keep $A_{\text{RT}} = 50.412727$ strictly dimensionless everywhere while ensuring complete arithmetic reproducibility from the text alone, let $\Phi_0 \equiv 1^\circ = \frac{\pi}{180}\text{ rad}$ denote the dimensional unit phase angle so that $A_{\text{RT}}\Phi_0 = 50.412727^\circ = 0.879868\text{ rad}$ defines the characteristic torsional phase scale per Earth mass $M_\oplus = 5.972 \times 10^{24}\text{ kg}$. Using the calibrated planetary orbital-impedance phase absorptions $\Phi_{p,i} = A_{\text{RT}}\Phi_0 (M_i / M_\oplus) w_i$, where $w_i = (a_\oplus / a_i)^{1/2} \eta_{\text{syn},i}$ is the product of each planet's Keplerian lever arm $(a_\oplus/a_i)^{1/2}$ and its empirical synodic lap coupling coefficient $\eta_{\text{syn},i} \equiv |T_{\text{syn},i}^{-1} / T_{\text{tach}}^{-1}|$ with the solar tachocline (explicitly listed below as calibrated solar-system inputs):

• Jupiter ($M_J = 317.83\,M_\oplus$, $a_J = 5.2044\text{ AU}$, $(a_\oplus/a_J)^{1/2} = 0.43835$, $\eta_{\text{syn},J} = 1.7003 \implies w_J = 0.74532$): $\Phi_{p,\text{Jup}} = 50.412727^\circ \times 317.83 \times 0.74532 = 11,942.0^\circ$ ($75.34\%$ of planetary total)
• Saturn ($M_S = 95.159\,M_\oplus$, $a_S = 9.5826\text{ AU}$, $(a_\oplus/a_S)^{1/2} = 0.32304$, $\eta_{\text{syn},S} = 1.7191 \implies w_S = 0.55534$): $\Phi_{p,\text{Sat}} = 50.412727^\circ \times 95.159 \times 0.55534 = 2,664.0^\circ$ ($16.81\%$)
• Neptune ($M_N = 17.147\,M_\oplus$, $a_N = 30.070\text{ AU}$, $(a_\oplus/a_N)^{1/2} = 0.18236$, $\eta_{\text{syn},N} = 3.9200 \implies w_N = 0.71485$): $\Phi_{p,\text{Nep}} = 50.412727^\circ \times 17.147 \times 0.71485 = 618.0^\circ$ ($3.90\%$)
• Uranus ($M_U = 14.536\,M_\oplus$, $a_U = 19.218\text{ AU}$, $(a_\oplus/a_U)^{1/2} = 0.22811$, $\eta_{\text{syn},U} = 3.1347 \implies w_U = 0.71506$): $\Phi_{p,\text{Ura}} = 50.412727^\circ \times 14.536 \times 0.71506 = 524.0^\circ$ ($3.31\%$)
• Inner Terrestrial Group ($\sum M_{\text{inner}} = 1.989\,M_\oplus$, $(a_\oplus/a_{\text{eff}})^{1/2} = 1.0000$, $\eta_{\text{syn,in}} = 1.0173 \implies w_{\text{in}} = 1.0173$): $\Phi_{p,\text{inner}} = 50.412727^\circ \times 1.989 \times 1.0173 = 102.0^\circ$ ($0.64\%$)

Summing all 8 planetary contributions gives the exact total $\Phi_p = 11,942.0^\circ + 2,664.0^\circ + 618.0^\circ + 524.0^\circ + 102.0^\circ = 15,850.0^\circ$ ($276.635\text{ rad}$), leaving an uncompensated torsional phase residual $\Delta\Phi = \Delta\Phi_{\text{Sun}} - \Phi_p = 16,554.8^\circ - 15,850.0^\circ = 704.8^\circ$ ($12.3011\text{ rad}$). Balancing $\Delta\Phi$ at the outer Oort standing-wave boundary across the $\pm \frac{\pi}{2}\beta_{\text{TRT}}$ impedance band yields the consistency estimate: \begin{equation} M_{\text{Oort}} = M_\oplus \left(\frac{\Delta\Phi}{A_{\text{RT}},\Phi_0}\right)\left(1 \pm \frac{\pi}{2}\beta_{\text{TRT}}\right) = M_\oplus \left(\frac{704.8^\circ}{50.412727 \times 1^\circ}\right)(1 \pm 0.0847) = 13.98 \pm 1.18,M_\oplus \tag{25} \end{equation} In Eq. (25), the angular units cancel explicitly ($\text{deg}/\text{deg} = \text{rad}/\text{rad} = 1$), keeping $A_{\text{RT}}$ strictly dimensionless and yielding an outer Oort anchor mass consistency range of $12.80\text{ to }15.16\,M_\oplus$ (classified in Table 4 as a calibrated solar-system consistency model rather than an ab initio prediction).

---

6. Zero-Mean Regularized Vorticity PDE, A Priori $L^2$ Enstrophy Bound & Exhaustive Orbital Falsifiability

Physical Picture: In classic fluid mechanics (the unmodified Navier–Stokes equations), a 3D fluid lacks an intrinsic upper bound on spin: in principle, vortex stretching can amplify vorticity ($\boldsymbol{\omega}$) toward infinity at a single point ("finite-time blowup"). In the two-fluid Quantum Plenum, however, a vortex core spinning near the medium's relaxation frequency $\omega_0 = \tau_0^{-1}$ experiences Hall–Vinen mutual friction against the $5.392\%$ normal-fluid component, venting excess rotational energy into longitudinal sound waves. Mathematically, this acts like a zero-mean solenoidal cubic pressure-relief valve ($-\beta_{\text{TRT}}\tau_0 \mathbb{P}_0(|\boldsymbol{\omega}|^2\boldsymbol{\omega})$): because the cubic valve dissipation ($|\boldsymbol{\omega}|^4$) grows faster than quadratic vortex stretching ($|\boldsymbol{\omega}|^2\|\nabla\mathbf{u}\|_{L^2}$), global $L^2$ enstrophy cannot blow up in finite time for smooth solutions. In the solar system, this same two-fluid medium produces two mechanical signatures: (1) a sunward spacecraft ram-impedance deceleration spanning $[8.124, 8.812] \times 10^{-10}\text{ m s}^{-2}$ (midpoint $\bar{a}_{\text{plenum}} = (8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$), which NASA's solar-powered Europa Clipper will test prospectively in 2028–2030, and (2) an Earth–Moon waveguide recession ceiling of $3.9214\text{ cm yr}^{-1}$ reduced by $2.696\%$ boundary-layer slip ($-0.1057\text{ cm yr}^{-1}$) to yield a net lunar recession rate of $3.816\text{ cm yr}^{-1}$.

6.1 Zero-Mean Solenoidal Leray–Helmholtz Regularized 3D Vorticity Equation

In classical 3D incompressible fluid mechanics ($\nabla \cdot \mathbf{u} = 0$), the vorticity $\boldsymbol{\omega} = \nabla \times \mathbf{u}$ $[\text{s}^{-1}]$ satisfies $\nabla \cdot \boldsymbol{\omega} = 0$ and on the periodic 3-torus $\mathbb{T}^3 = [0, L]^3$ of finite volume $V = L^3$, must satisfy the zero spatial mean condition $\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$ identically (since the integral of any curl over a closed boundaryless manifold vanishes).

When local vorticity approaches $\tau_0^{-1} = \omega_0$, two-fluid mutual friction sheds rotational energy into acoustic modes at rate $\beta_{\text{TRT}}\tau_0 |\boldsymbol{\omega}|^2\boldsymbol{\omega}$. To guarantee that the damping term simultaneously preserves: (i) zero divergence ($\nabla \cdot \boldsymbol{\omega} = 0$) and (ii) the zero-mean condition on $\mathbb{T}^3$ ($\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$), we introduce the Zero-Mean Solenoidal Projector $\mathbb{P}_0$: \begin{equation} \mathbb{P}0 \mathbf{F} \equiv \mathbb{P}\mathbf{F} - \frac{1}{V}\int{\mathbb{T}^3} \mathbb{P}\mathbf{F},d^3x \tag{26a} \end{equation} (where on $\mathbb{R}^3$ with rapid decay at infinity, the mean is already zero and $\mathbb{P}_0 \equiv \mathbb{P}$). The resulting Curl-Compatible TRT-Regularized 3D Vorticity Equation is: \begin{equation} \frac{\partial \boldsymbol{\omega}}{\partial t} + (\mathbf{u} \cdot \nabla) \boldsymbol{\omega} = (\boldsymbol{\omega} \cdot \nabla) \mathbf{u} + \nu_0 \nabla^2 \boldsymbol{\omega} - \beta_{\text{TRT}},\tau_0 , \mathbb{P}_0!\left(|\boldsymbol{\omega}|^2 \boldsymbol{\omega}\right), \qquad \mathbf{u} = \nabla \times (-\Delta)^{-1}\boldsymbol{\omega} \tag{26} \end{equation} Verification of Curl Invariance, Mean Preservation, and SI Dimensions:

1. Preservation of $\nabla \cdot \boldsymbol{\omega} = 0$ and $\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$: Because the constant vector $\frac{1}{V}\int_{\mathbb{T}^3}\mathbb{P}\mathbf{F}\,d^3x$ is spatially uniform, $\nabla \cdot (\mathbb{P}_0 \mathbf{F}) = \nabla \cdot (\mathbb{P}\mathbf{F}) - 0 = 0$ identically. Furthermore, integrating $\mathbb{P}_0 \mathbf{F}$ over $\mathbb{T}^3$ gives $\int_{\mathbb{T}^3}\mathbb{P}\mathbf{F}\,d^3x - \frac{V}{V}\int_{\mathbb{T}^3}\mathbb{P}\mathbf{F}\,d^3x = \mathbf{0}$ identically. Thus, if $\boldsymbol{\omega}_0$ has zero divergence and zero spatial mean, $\boldsymbol{\omega}(\cdot, t)$ remains strictly divergence-free and zero-mean for all $t \ge 0$, guaranteeing that $\mathbf{u} = \nabla \times (-\Delta)^{-1}\boldsymbol{\omega}$ is uniquely defined on $\mathbb{T}^3$ for all time.
2. Dimensional Homogeneity: Since $[\beta_{\text{TRT}}] = 1$, $[\tau_0] = \text{s}$, $[\mathbb{P}_0] = 1$, and $[|\boldsymbol{\omega}|^2\boldsymbol{\omega}] = \text{s}^{-3}$, the damping term has exact SI units of $[\text{s} \cdot \text{s}^{-3}] = [\text{s}^{-2}]$, matching $\partial_t \boldsymbol{\omega}$.

6.2 Theorem 1: A Priori Global $L^2$ Enstrophy Bound for Smooth Solutions

(Scope Qualification: The Clay Mathematics Institute Millennium Problem concerns the unmodified continuum Navier–Stokes system with $\beta_{\text{TRT}} = 0$. Theorem 1 establishes a rigorous a priori $L^2$ enstrophy bound and $L^4_t L^4_x$ vorticity spacetime integrability for smooth solutions of the physically regularized TRT vorticity system (26) with $\beta_{\text{TRT}}\tau_0 > 0$ on $\mathbb{R}^3$ and the periodic 3-torus $\mathbb{T}^3$.)

Theorem 1 (A Priori Enstrophy Boundedness on $\mathbb{R}^3$ and $\mathbb{T}^3$): Let $\Omega$ be either the whole space $\mathbb{R}^3$ (with rapid decay at infinity) or the periodic 3-torus $\mathbb{T}^3 = [0, L]^3$ (of finite volume $V = L^3 < \infty$), and let $\boldsymbol{\omega}_0 \in L^2_\sigma(\Omega)$ with $\nabla \cdot \boldsymbol{\omega}_0 = 0$, and on $\mathbb{T}^3$, with zero spatial mean $\int_{\mathbb{T}^3}\boldsymbol{\omega}_0(\mathbf{x})\,d^3x = \mathbf{0}$. For any smooth solution $\boldsymbol{\omega}(\cdot, t)$ of Eq. (26) existing on $[0, T]$ with $\nu_0 > 0$ and $\beta_{\text{TRT}}\tau_0 > 0$:

(a) On both $\Omega = \mathbb{R}^3$ and $\Omega = \mathbb{T}^3$, the enstrophy $\|\boldsymbol{\omega}(t)\|_{L^2(\Omega)}^2$ satisfies the finite-horizon Grönwall a priori bound: \begin{equation} |\boldsymbol{\omega}(t)|{L^2(\Omega)}^2 + \int_0^t \left(2\nu_0 |\nabla\boldsymbol{\omega}(s)|{L^2(\Omega)}^2 + \beta_{\text{TRT}}\tau_0 |\boldsymbol{\omega}(s)|{L^4(\Omega)}^4\right) ds \le |\boldsymbol{\omega}0|{L^2(\Omega)}^2 \exp!\left(\frac{t}{\beta{\text{TRT}}\tau_0}\right) < \infty \quad \forall t \in [0, T] \tag{27} \end{equation} (b) Moreover, on the periodic 3-torus $\Omega = \mathbb{T}^3 = [0, L]^3$, if the smooth solution continues for all $t \ge 0$, the enstrophy is uniformly bounded for all $t \ge 0$ by the time-independent absorbing ball: \begin{equation} |\boldsymbol{\omega}(t)|{L^2(\mathbb{T}^3)}^2 \le \max!\left(|\boldsymbol{\omega}0|{L^2(\mathbb{T}^3)}^2, ; \frac{V}{\beta{\text{TRT}}^2 \tau_0^2}\right) < \infty \quad \forall t \ge 0 \tag{28} \end{equation}

Proof: Step 1 (Zero-Mean Projector Orthogonality & Exact $L^4$ Dissipation): Taking the $L^2(\Omega)$ inner product of Eq. (26) with $\boldsymbol{\omega}$, the advection term vanishes by incompressibility ($\int_\Omega \boldsymbol{\omega}\cdot(\mathbf{u}\cdot\nabla)\boldsymbol{\omega}\,d^3x = 0$), and the viscous term integrates by parts to $-\nu_0\|\nabla\boldsymbol{\omega}\|_{L^2(\Omega)}^2$. For the damping term, because $\boldsymbol{\omega}$ has zero spatial mean on $\mathbb{T}^3$ ($\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$), the spatially constant mean-subtraction in $\mathbb{P}_0$ drops out identically: $\int_{\mathbb{T}^3} \boldsymbol{\omega} \cdot \left[ \frac{1}{V}\int_{\mathbb{T}^3}\mathbb{P}(|\boldsymbol{\omega}|^2\boldsymbol{\omega})\,d^3x' \right] d^3x = \left(\frac{1}{V}\int_{\mathbb{T}^3}\mathbb{P}(|\boldsymbol{\omega}|^2\boldsymbol{\omega})\,d^3x'\right) \cdot \int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = 0$ Because $\mathbb{P}$ is self-adjoint and $\mathbb{P}\boldsymbol{\omega} = \boldsymbol{\omega}$, we obtain: $\int_\Omega \boldsymbol{\omega} \cdot \mathbb{P}_0(|\boldsymbol{\omega}|^2\boldsymbol{\omega})\,d^3x = \int_\Omega \boldsymbol{\omega} \cdot \mathbb{P}(|\boldsymbol{\omega}|^2\boldsymbol{\omega})\,d^3x = \int_\Omega (\mathbb{P}\boldsymbol{\omega})\cdot(|\boldsymbol{\omega}|^2\boldsymbol{\omega})\,d^3x = \int_\Omega |\boldsymbol{\omega}|^4\,d^3x = \|\boldsymbol{\omega}\|_{L^4(\Omega)}^4$ Thus the exact enstrophy balance equation is: \begin{equation} \frac{1}{2}\frac{d}{dt}|\boldsymbol{\omega}|{L^2(\Omega)}^2 + \nu_0 |\nabla \boldsymbol{\omega}|{L^2(\Omega)}^2 + \beta_{\text{TRT}}\tau_0 |\boldsymbol{\omega}|{L^4(\Omega)}^4 = \int\Omega \boldsymbol{\omega} \cdot (\boldsymbol{\omega} \cdot \nabla)\mathbf{u},d^3x \tag{29} \end{equation}

Step 2 (Degree-3 $L^4 \times L^4 \times L^2$ Hölder Estimate & Plancherel Identity): By Hölder's inequality with exponents $(p, q, r) = (4, 4, 2)$ satisfying $\frac{1}{4} + \frac{1}{4} + \frac{1}{2} = 1$: \begin{equation} \left|\int_\Omega \boldsymbol{\omega} \cdot (\boldsymbol{\omega} \cdot \nabla)\mathbf{u},d^3x\right| \le |\boldsymbol{\omega}|{L^4(\Omega)}^2 , |\nabla \mathbf{u}|{L^2(\Omega)} \tag{30} \end{equation} On $\mathbb{T}^3$, because $\int_{\mathbb{T}^3}\boldsymbol{\omega}\,d^3x = \mathbf{0}$, the Fourier zero mode $\hat{\boldsymbol{\omega}}(\mathbf{0})$ vanishes identically, ensuring $\hat{\mathbf{u}}(\mathbf{k}) = i \mathbf{k} \times \hat{\boldsymbol{\omega}}(\mathbf{k})/|\mathbf{k}|^2$ is completely non-singular. By Plancherel's theorem, $\|\nabla \mathbf{u}\|_{L^2(\Omega)} = \|\boldsymbol{\omega}\|_{L^2(\Omega)}$ with exact constant $1$, yielding: \begin{equation} \left|\int_\Omega \boldsymbol{\omega} \cdot (\boldsymbol{\omega} \cdot \nabla)\mathbf{u},d^3x\right| \le |\boldsymbol{\omega}|{L^4(\Omega)}^2 , |\boldsymbol{\omega}|{L^2(\Omega)} \tag{31} \end{equation}

Step 3 (Elementary Quadratic Young Inequality & Absorption): Applying the quadratic Young inequality $a b \le \frac{\epsilon}{2}a^2 + \frac{1}{2\epsilon}b^2$ with $a = \|\boldsymbol{\omega}\|_{L^4(\Omega)}^2$, $b = \|\boldsymbol{\omega}\|_{L^2(\Omega)}$, and $\epsilon = \beta_{\text{TRT}}\tau_0 > 0$: \begin{equation} |\boldsymbol{\omega}|{L^4(\Omega)}^2 , |\boldsymbol{\omega}|{L^2(\Omega)} \le \frac{\beta_{\text{TRT}}\tau_0}{2}|\boldsymbol{\omega}|{L^4(\Omega)}^4 + \frac{1}{2\beta{\text{TRT}}\tau_0}|\boldsymbol{\omega}|{L^2(\Omega)}^2 \tag{32} \end{equation} Substituting Eq. (32) into Eq. (29), absorbing $\frac{\beta_{\text{TRT}}\tau_0}{2}\|\boldsymbol{\omega}\|_{L^4}^4$, and multiplying by $2$: \begin{equation} \frac{d}{dt}|\boldsymbol{\omega}|{L^2(\Omega)}^2 + 2\nu_0 |\nabla\boldsymbol{\omega}|{L^2(\Omega)}^2 + \beta{\text{TRT}}\tau_0 |\boldsymbol{\omega}|{L^4(\Omega)}^4 \le \frac{1}{\beta{\text{TRT}}\tau_0}|\boldsymbol{\omega}|_{L^2(\Omega)}^2 \tag{33} \end{equation} Integrating Eq. (33) via Grönwall's inequality establishes Part (a) (Eq. 27) on $[0, T]$.

Step 4 (Logistic Riccati Absorbing Ball on $\mathbb{T}^3$): On $\mathbb{T}^3$ of volume $V$, Cauchy–Schwarz gives $\|\boldsymbol{\omega}\|_{L^2}^2 \le V^{1/2}\|\boldsymbol{\omega}\|_{L^4}^2$, hence $-\|\boldsymbol{\omega}\|_{L^4}^4 \le -\frac{1}{V}\|\boldsymbol{\omega}\|_{L^2}^4$. Dropping the viscous term $2\nu_0\|\nabla\boldsymbol{\omega}\|_{L^2}^2 \ge 0$, the enstrophy $Y(t) \equiv \|\boldsymbol{\omega}(t)\|_{L^2(\mathbb{T}^3)}^2$ obeys: \begin{equation} \frac{dY}{dt} \le \frac{1}{\beta_{\text{TRT}}\tau_0} Y - \frac{\beta_{\text{TRT}}\tau_0}{V} Y^2 = \frac{\beta_{\text{TRT}}\tau_0}{V} Y \left(\frac{V}{\beta_{\text{TRT}}^2\tau_0^2} - Y\right) \tag{34} \end{equation} By the standard ODE comparison principle, $Y(t) \le \max\!\left(Y(0), \frac{V}{\beta_{\text{TRT}}^2\tau_0^2}\right)$ for all $t \ge 0$, establishing Part (b) (Eq. 28). $\blacksquare$

6.3 Two-Fluid Momentum-Flux Derivation of the Spacecraft Ram-Impedance Law, Europa Clipper Error Budget & Exhaustive Decision Protocol

1. Four-Step Two-Fluid Momentum-Flux Derivation of the Spacecraft Ram-Impedance Force Law (Eq. 35): We derive the hyperbolic spacecraft deceleration law directly from two-fluid acoustic momentum flux. Consider a macroscopic spacecraft of inertial rest mass $m_{\text{sc}}$ moving radially outward at hyperbolic cruise speed $\mathbf{v} = v_r \hat{\mathbf{r}}$ ($v_r \ll c_0$) through the expanding Hubble–Plenum background:

   • Step 1 (Cosmological Metric Cavity Expansion Acceleration $a_H$): In the expanding cosmological cavity of radius $R_H = c_0 / H_{\text{early}}$, the background acoustic wave frame expands radially with characteristic cosmological acceleration: $a_H \equiv c_0 H_{\text{early}} = (2.99792458 \times 10^8\text{ m s}^{-1})(2.1843 \times 10^{-18}\text{ s}^{-1}) = 6.5483 \times 10^{-10}\text{ m s}^{-2}$
   • Step 2 (Characteristic Acoustic Cycle Impedance $Z_{\text{cyc}}$): As the spacecraft's constituent baryonic solitons traverse the unperturbed two-fluid substrate over one full carrier wave cycle ($\Delta\phi_{\text{cyc}} = 2\pi\text{ rad}$), the coupled wave-packet envelope encounters the inviscid superfluid cycle phase $2\pi$ plus the dissipative boundary normal-fluid mass fraction $\chi_n = \beta_{\text{TRT}}$ (Proposition 1), defining the dimensionless cycle impedance: $Z_{\text{cyc}} \equiv 2\pi + \beta_{\text{TRT}} = 6.283185 + 0.053921 = 6.33711 \quad [\text{dimensionless: } 1]$
   • Step 3 (3D Dipole Wake Meridian Projection $A_{\text{eff}}$): A translating macroscopic body in a compressible acoustic medium generates an axisymmetric 3D acoustic dipole wake whose standing half-wave nodes span two orthogonal meridian cross-sections ($\Delta\theta_1 \times \Delta\theta_2 = \pi \times \pi = \pi^2\text{ rad}^2$). Projecting the bulk acoustic susceptibility $A_{\text{RT}}$ across this orthogonal meridian phase area defines the effective dipole wake susceptibility: $A_{\text{eff}} \equiv \frac{A_{\text{RT}}}{\pi^2} = \frac{50.412727}{9.869604} = 5.10788 \quad [\text{dimensionless: } 1]$
   • Step 4 (Net Inertial Drag Per Unit Mass from Carrier Back-Reaction): By two-fluid momentum conservation across the dipole wake control volume, the net sunward acoustic carrier back-reaction force $\mathbf{F}_{\text{plenum}} = m_{\text{sc}}\mathbf{a}_{\text{plenum}}^{(0)}$ opposing hyperbolic radial motion scales with $m_{\text{sc}} a_H$ amplified by the cycle-to-wake impedance ratio $Z_{\text{cyc}}/A_{\text{eff}}$ (Postulate P5). Dividing by $m_{\text{sc}}$ yields the exact inner-heliospheric uncompressed baseline deceleration $a_{\text{min}} = |\mathbf{a}_{\text{plenum}}^{(0)}|$: \begin{equation} \mathbf{a}{\text{plenum}}^{(0)} = -a_H \left(\frac{Z{\text{cyc}}}{A_{\text{eff}}}\right)\hat{\mathbf{r}} = -c_0 H_{\text{early}} \left(\frac{2\pi + \beta_{\text{TRT}}}{A_{\text{RT}} / \pi^2}\right) \hat{\mathbf{r}} = -\left(6.5483 \times 10^{-10}\text{ m s}^{-2}\right)\left(\frac{6.33711}{5.10788}\right)\hat{\mathbf{r}} = -8.124 \times 10^{-10}\text{ m s}^{-2},\hat{\mathbf{r}} \tag{35} \end{equation} In the outer heliospheric waveguide regime, inclusion of the quarter-wave horizon compression factor $(1 + \frac{\pi}{2}\beta_{\text{TRT}}) = 1.0847$ defines the outer-heliospheric compressed ceiling $a_{\text{max}} \equiv a_{\text{min}}(1 + \frac{\pi}{2}\beta_{\text{TRT}}) = 8.812 \times 10^{-10}\text{ m s}^{-2}$. Across the radial transition regime $[a_{\text{min}}, a_{\text{max}}] = [8.124, 8.812] \times 10^{-10}\text{ m s}^{-2}$, the exact midpoint $\bar{a}_{\text{plenum}} = \frac{1}{2}(a_{\text{min}} + a_{\text{max}})$ and model half-width $\Delta a_{\text{band}} = \frac{1}{2}(a_{\text{max}} - a_{\text{min}})$ are: $\bar{a}_{\text{plenum}} = a_{\text{min}}\left(1 + \frac{\pi}{4}\beta_{\text{TRT}}\right) = 8.468 \times 10^{-10}\text{ m s}^{-2} \approx 8.47 \times 10^{-10}\text{ m s}^{-2}, \qquad \Delta a_{\text{band}} = \pm 0.344 \times 10^{-10}\text{ m s}^{-2}$ yielding the exact theoretical model band $(8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$ (and when combined linearly with the $\pm 0.65 \times 10^{-10}\text{ m s}^{-2}$ $3\sigma$ spacecraft systematic floor from Table 3, a total observational acceptance interval of $[7.47, 9.47] \times 10^{-10}\text{ m s}^{-2}$, consistent with the historical Pioneer 10/11 residual of $(8.74 \pm 1.33) \times 10^{-10}\text{ m s}^{-2}$ [10]).

2. Europa Clipper 2028–2030 Navigation Systematics & Doppler Error Budget: Turyshev et al. (2012) [10] showed that on Pioneer 10/11, anisotropic thermal radiation recoil from four nuclear SNAP-19 RTGs could account for $\sim 8 \times 10^{-10}\text{ m s}^{-2}$ of sunward acceleration. NASA's Europa Clipper (outer cruise $4.0\text{ AU} \le r \le 5.2\text{ AU}$, 2028–2030) eliminates RTG thermal recoil completely by carrying zero RTGs and using symmetric $100\text{ m}^2$ solar wings. Table 3 documents the quantitative non-gravitational acceleration error budget during quiet outer-cruise X/Ka-band Doppler tracking windows (grounded in standard deep-space radiometric navigation flight standards; e.g., Moyer 2000 [11]), demonstrating a combined $3\sigma$ systematic uncertainty floor of $\sigma_{3\sigma} = 0.65 \times 10^{-10}\text{ m s}^{-2}$—more than an order of magnitude below the predicted TRT signal ($\bar{a}_{\text{plenum}} = 8.47 \times 10^{-10}\text{ m s}^{-2}$).

Table 3: Europa Clipper Outer-Cruise ($4.0\text{–}5.2\text{ AU}$) Non-Gravitational Acceleration & Doppler Error Budget

Pre-Registered Exhaustive Operational Decision Protocol for Europa Clipper Outer Cruise: To ensure unambiguous, mathematically rigorous falsifiability, all possible prospective Doppler acceleration residual outcomes $a_{\text{res}} \in (-\infty, +\infty)$ from Europa Clipper's 2028–2030 outer-cruise Doppler tracking ($4.0\text{–}5.2\text{ AU}$) are partitioned into four mutually exclusive and collectively exhaustive operational regimes:

• Coordinate & Sign Convention: The observable $a_{\text{res}}$ is defined as the signed radial line-of-sight acceleration anomaly, where positive values ($a_{\text{res}} > 0$) denote an anomalous sunward (inward) deceleration force opposing outward spacecraft motion, and negative values ($a_{\text{res}} < 0$) denote an anomalous anti-sunward (outward) acceleration force.

1. Regime 1 — Definitive Null Falsification ($|a_{\text{res}}| < 2.00 \times 10^{-10}\text{ m s}^{-2}$): Refutes the predicted Plenum ram-impedance deceleration (Eq. 35) at $>9\sigma$ statistical confidence relative to the $1\sigma = 0.22 \times 10^{-10}\text{ m s}^{-2}$ navigation uncertainty floor (i.e. $-2.00 \times 10^{-10} < a_{\text{res}} < 2.00 \times 10^{-10}\text{ m s}^{-2}$), falsifying the Postulate P5 spacecraft impedance closure.
2. Regime 2 — Confirmatory TRT Detection Band ($7.47 \times 10^{-10} \le a_{\text{res}} \le 9.47 \times 10^{-10}\text{ m s}^{-2}$): Statistically consistent with the predicted TRT model band $(8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$ within the combined $3\sigma$ navigation systematic floor ($\pm 0.65 \times 10^{-10}\text{ m s}^{-2}$).
3. Regime 3 — Intermediate Transition / Inconclusive Zone ($2.00 \times 10^{-10} \le a_{\text{res}} < 7.47 \times 10^{-10}\text{ m s}^{-2}$): Requires joint covariance re-estimation of solar radiation pressure, optical surface degradation, and bus thermal louvers before definitive classification.
4. Regime 4 — Off-Band Falsification ($a_{\text{res}} \le -2.00 \times 10^{-10}\text{ m s}^{-2}$ OR $a_{\text{res}} > 9.47 \times 10^{-10}\text{ m s}^{-2}$): Any negative (anti-sunward) residual exceeding noise, or any positive sunward residual exceeding $9.47 \times 10^{-10}\text{ m s}^{-2}$, decisively refutes the specific TRT impedance coupling closure (Eq. 35).

Every real value $a_{\text{res}} \in (-\infty, +\infty)$ is assigned to exactly one of the four operational regimes without omission or ambiguity.

1. Reconciled Secular Lunar Recession Law & Dynamical Boundary Governor Model: In the TRT Earth–Moon waveguide model (Postulate P5), the bimodal $2 N_5 = 100$ harmonic of the carrier wavelength $L_0 = 3.92144\text{ m}$ over the terrestrial orbital period $\tau_{\text{yr}} = 3.15576 \times 10^7\text{ s}$ defines the unattenuated inviscid waveguide expansion ceiling $\dot{r}_0$: \begin{equation} \dot{r}0 \equiv \frac{L_0}{2 N_5 , \tau{\text{yr}}} = \frac{3.92144\text{ m}}{100\text{ yr}} = 3.92144\text{ cm yr}^{-1} \quad [\text{SI: } 1.2426 \times 10^{-9}\text{ m s}^{-1}] \tag{36} \end{equation} In celestial mechanics, terrestrial ocean-tidal friction is the dissipative fluid mechanism through which angular momentum is physically transferred. However, integrating standard ocean-tidal dissipation models backward in time produces the well-known Gerstenkorn timing anomaly (collapsing the Moon onto Earth within just $\sim 1.5\text{–}2.0\text{ Gyr}$, in sharp conflict with $3.8\text{–}4.4\text{ Gyr}$ zircon geological records). In TRT, the macroscopic $2N_5 = 100$ acoustic standing waveguide acts as a macroscopic dynamical boundary governor / resonant attractor: ocean-tidal dissipation self-organizes to match the allowed phase-locked cavity expansion rate governed by the Plenum boundary-layer slip (Proposition 1).

   Evaluating the net secular recession rate under both the leading-order linear thin-sheath approximation ($\delta^*/R = \frac{1}{2}\beta_{\text{TRT}} = 0.026960$) and the exact quadratic root ($\delta^*/R = 1 - \sqrt{1 - \beta_{\text{TRT}}} = 0.027334$) yields: \begin{align} \dot{r}{\text{net}}^{\text{(lin)}} &= \dot{r}0\left(1 - \frac{1}{2}\beta{\text{TRT}}\right) = 3.92144\text{ cm yr}^{-1} \times (1 - 0.026960) = 3.8157\text{ cm yr}^{-1} \approx 3.816\text{ cm yr}^{-1} \tag{37a} \ \dot{r}{\text{net}}^{\text{(exact)}} &= \dot{r}0\left(1 - (1 - \sqrt{1 - \beta{\text{TRT}}})\right) = 3.92144\text{ cm yr}^{-1} \times (1 - 0.027334) = 3.8143\text{ cm yr}^{-1} \approx 3.814\text{ cm yr}^{-1} \tag{37b} \end{align} The exact quadratic and linear thin-sheath closures span the truncation interval $[3.8143, 3.8157]\text{ cm yr}^{-1}$ (midpoint $3.8150 \pm 0.0007\text{ cm yr}^{-1}$, $0.037\%$ truncation spread), which agrees with the present-day Lunar Laser Ranging (LLR) measurement ($\dot{r}_{\text{obs}} = 3.82 \pm 0.07\text{ cm yr}^{-1}$) to within $0.07\sigma$. As an illustrative constant-rate kinematic baseline (noting that detailed paleotidal evolution involves variable ocean-basin resonance geometry over geological epochs), linear backward extrapolation of $\dot{r}_{\text{net}} = 3.815\text{ cm yr}^{-1}$ over $4.50 \times 10^9\text{ yr}$ gives $\Delta r = 171,675\text{ km}$ ($r_0 = 384,400 - 171,675 = 212,725\text{ km} \gg r_{\text{Roche}} \approx 9,500\text{ km}$), illustrating how waveguide-governed recession avoids the unphysical $1.5\text{ Gyr}$ Gerstenkorn collapse.

---

7. Hydrostatic Realization of Newton's Gravitational Law, Positive Soliton Hamiltonian & The Continuous Non-Linear Acoustic Pressure PDE

Physical Picture: In a compressible fluid, what we experience as "gravity" is modeled as a static acoustic pressure-depletion gradient rather than curvature of an empty void. Because a baryonic mass distribution $\rho_M(\mathbf{x})$ continuously cycles and scatters coherent acoustic wave energy at its vortex core boundaries, it establishes a steady, static radiation-pressure depression below the ambient far-field pressure $P_0$. Crucially, the unknotted ambient Quantum Plenum substrate itself remains in static hydrostatic equilibrium ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), held balanced by carrier radiation reaction $\mathbf{f}_{\text{rad}} = +\nabla P$ so there is no material drainage of the vacuum medium into the body. Why, then, does a matter soliton accelerate inward if the ambient substrate is static, and why doesn't its cavitated core float outward like a bubble? First, a topological soliton is not a passive parcel of unknotted ambient fluid: it is a localized, phase-locked vortex singularity whose closed sonic-horizon boundary scatters incident carrier phonons asymmetrically, generating the wave-gradient Bernoulli force $\mathbf{F}_{\text{grad}} = -(m_i/\rho_0)\nabla P$. Second, its inertial and gravitational mass $m_i$ is determined not by its static core density deficit, but by its total integrated two-fluid Hamiltonian energy $m_i c_0^2 = \int_{\mathcal{V}_{\text{sol}}} (\frac{1}{2}\rho_s|\mathbf{u}_s|^2 + \frac{1}{2}\rho_n|\mathbf{u}_n|^2 + \Delta\mathcal{U}_{\text{pot}})d^3x > 0$: because the relativistic toroidal circulation kinetic energy vastly exceeds the shallow core cavitation deficit, $m_i > 0$ is strictly positive, forcing the soliton inward toward lower acoustic pressure ($\mathbf{a} = -(1/\rho_0)\nabla P$). Finally, a single continuous non-linear acoustic pressure PDE $\nabla \cdot [\mu(|\nabla P|/(\rho_0 a_{\text{gal}}))\nabla P] = 4\pi G \rho_0 \rho_M$ governs the pressure field across all scales—recovering $\nabla^2 P = 4\pi G \rho_0 \rho_M$ ($a = -GM/r^2$) in the solar system and transitioning smoothly to $|\nabla P|/\rho_0 = \sqrt{GM a_{\text{gal}}}/r$ ($v_{\text{rot}}^4 = GM a_{\text{gal}}$) in the deep galactic disc regime.

7.1 Weak-Field Acoustic Gravity, Continuum Radiation Reaction & Positive Soliton Hamiltonian Mass

We first examine the gravitational sector in the high-gradient solar-system and laboratory regime ($|\nabla P|/(\rho_0 a_{\text{gal}}) \gg 1$ and $|\Phi(\mathbf{x})|/c_0^2 = GM/(r c_0^2) \ll 1$), where fractional density variations satisfy $|\delta\rho|/\rho_0 = |\Delta P|/P_0 = GM/(r c_0^2) \ll 1$ (e.g., $\sim 10^{-6}$ at the solar surface and $\sim 10^{-9}$ at Earth's surface) and $\rho(\mathbf{x}) = \rho_0[1 + \mathcal{O}(|\Phi|/c_0^2)] \approx \rho_0$.

Under Postulate P6 (Table 1), in the high-gradient limit $|\nabla P|/(\rho_0 a_{\text{gal}}) \gg 1$ where the constitutive flux function $\mu(x) = x/\sqrt{1+x^2} \to 1$, a stationary baryonic mass distribution $\rho_M(\mathbf{x})$ [$\text{kg m}^{-3}$] of total mass $M = \int_{\mathbb{R}^3} \rho_M(\mathbf{x})\,d^3x$ in the static background Quantum Plenum ($\mathbf{u}_0 = \mathbf{0}$, $\mathbf{j}_0 = \rho_s \mathbf{u}_s + \rho_n \mathbf{u}_n = \mathbf{0}$) establishes a static acoustic pressure depletion governed by the linearized Poisson equation: \begin{equation} \nabla^2 P(\mathbf{x}) = 4\pi G \rho_0 , \rho_M(\mathbf{x}), \qquad \lim_{|\mathbf{x}|\to\infty} P(\mathbf{x}) = P_0 \tag{38} \end{equation}

1. Exact Source-Free Mass Continuity Harmony

In the exterior vacuum region outside the compact baryonic source ($r > R_{\text{core}}$, where $\rho_M = 0$), Eq. (38) reduces to Laplace's equation for the acoustic pressure field: \begin{equation} \nabla^2 P(\mathbf{x}) = 0 \quad \text{for } |\mathbf{x}| > R_{\text{core}} \tag{38a} \end{equation} Because the physical substrate medium is static ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), the steady mass continuity equation (Postulate P1, Eq. 1) satisfies: \begin{equation} \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}_{\text{substrate}}) = 0 + \nabla \cdot (\rho_0 \cdot \mathbf{0}) = 0 \tag{38b} \end{equation} identically and trivially at every exterior point $r > R_{\text{core}}$. Unlike ad-hoc hydrodynamic sink geometries that assert continuous fluid inflow into a non-singular body (which would violate source-free continuity), TRT's gravitational field is a pure static acoustic radiation-pressure depression with zero material substrate velocity ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), requiring no unphysical mass sinks or fluid drainage in the exterior vacuum.

Applying Gauss's divergence theorem to Eq. (38) over a sphere of radius $r > R_{\text{core}}$, $\oint_{\partial B_r} \nabla P \cdot d\mathbf{A} = 4\pi r^2 \frac{dP}{dr} = 4\pi G \rho_0 M$, yields the exact outward-increasing radial pressure gradient and static pressure profile: \begin{equation} \frac{dP(r)}{dr} = +\frac{\rho_0 G M}{r^2} > 0, \qquad P(r) = P_0 - \int_r^\infty \frac{\rho_0 G M}{{r'}^2},dr' = P_0 - \frac{\rho_0 G M}{r} \implies \Delta P(r) \equiv P(r) - P_0 = -\frac{\rho_0 G M}{r} < 0 \tag{39} \end{equation}

2. Hydrostatic Momentum Balance of the Unknotted Ambient Continuum

In Section 2.1, the linearized momentum conservation equation of the ambient two-fluid Plenum continuum couples the pressure gradient $-\nabla P$ to the body-force density $\mathbf{f}_{\text{rad}}$ exerted by acoustic carrier radiation: \begin{equation} \rho_0 \frac{\partial \mathbf{u}{\text{substrate}}}{\partial t} = -\nabla P + \mathbf{f}{\text{rad}} \tag{39a} \end{equation} For the ambient unknotted substrate to remain in exact static hydrostatic equilibrium ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$, $\partial_t \mathbf{u}_{\text{substrate}} = \mathbf{0}$), Eq. (39a) requires exact mechanical force cancellation everywhere in the exterior continuum: \begin{equation} -\nabla P(\mathbf{x}) + \mathbf{f}{\text{rad}}(\mathbf{x}) = \mathbf{0} \implies \mathbf{f}{\text{rad}}(\mathbf{x}) = +\nabla P(\mathbf{x}) = +\frac{\rho_0 G M}{r^2}\hat{\mathbf{r}} \tag{39b} \end{equation} Physical Mechanism: Baryonic matter consists of phase-locked topological vortex solitons (Postulate P3) that continuously scatter and absorb the isotropic high-frequency acoustic carrier modes ($f_0$) of the Quantum Plenum. The momentum transferred from the scattered carrier wave flux exerts an outward-directed acoustic radiation reaction body-force density $\mathbf{f}_{\text{rad}}(\mathbf{x}) = +\frac{\rho_0 G M}{r^2}\hat{\mathbf{r}}$ upon the surrounding unknotted continuum. This outward radiation reaction force precisely counterbalances the inward pressure-gradient body force $-\nabla P(\mathbf{x}) = -\frac{\rho_0 G M}{r^2}\hat{\mathbf{r}}$ at every radius $r > R_{\text{core}}$, holding the ambient continuum in exact static hydrostatic equilibrium ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$).

3. Topological Soliton Wave-Gradient Radiation Force & Universal Free-Fall (Weak Equivalence Principle)

We now derive the gravitational acceleration of a localized physical test body. By Postulate P3, elementary fermions and composite nuclei are localized topological vortex solitons carrying total rest energy $E_i = m_i c_0^2$. In a barotropic acoustic medium where $P = c_0^2 \rho$, the acoustic enthalpy perturbation per unit mass is $w(\mathbf{x}) = \int_{P_0}^{P(\mathbf{x})} \frac{dP'}{\rho'} \approx \frac{\Delta P(\mathbf{x})}{\rho_0}$. A localized soliton of integrated inertial mass $m_i = E_i / c_0^2$ immersed in the background pressure variation $\Delta P(\mathbf{x})$ possesses the acoustic interaction potential energy: \begin{equation} U_i(\mathbf{x}) = \int_{\mathcal{V}{\text{sol}}} \frac{\mathcal{H}{\text{sol}}(\mathbf{x}')}{c_0^2} w(\mathbf{x}'),d^3x' = m_i \frac{\Delta P(\mathbf{x})}{\rho_0} \tag{39c} \end{equation} Under wave-action conservation and acoustic radiation theory, a localized wave packet or topological soliton of rest energy $E_i = m_i c_0^2$ in a background pressure gradient experiences the wave-gradient acoustic radiation force (the generalized Bernoulli–Bjerknes wave-gradient force): \begin{equation} \mathbf{F}_{\text{grad}} = -\nabla U_i(\mathbf{x}) = -\frac{E_i}{\rho_0 c_0^2}\nabla P(\mathbf{x}) = -\frac{m_i c_0^2}{\rho_0 c_0^2}\nabla P(\mathbf{x}) = -\frac{m_i}{\rho_0}\nabla P(\mathbf{x}) \tag{39d} \end{equation} Applying Newton's second law $\mathbf{F}_{\text{grad}} = m_i \mathbf{a}$ to determine the physical center-of-mass acceleration $\mathbf{a}$ of the soliton: \begin{equation} m_i \mathbf{a} = -\frac{m_i}{\rho_0}\nabla P(\mathbf{x}) \implies \mathbf{a} = -\frac{1}{\rho_0}\nabla P(\mathbf{x}) \tag{40} \end{equation} Because the body's inertial mass $m_i$ appears linearly in both the acoustic wave-gradient force $\mathbf{F}_{\text{grad}}$ and the inertial resistance $m_i \mathbf{a}$, the inertial mass $m_i$ cancels identically. Substituting the radial pressure gradient from Eq. (39) into Eq. (40) yields: \begin{equation} \mathbf{a}(\mathbf{x}) = -\frac{1}{\rho_0}\left(+\frac{\rho_0 G M}{r^2}\hat{\mathbf{r}}\right) = -\frac{G M}{r^2}\hat{\mathbf{r}} \quad \left[\text{SI: } \frac{1}{\text{kg m}^{-3}}\frac{\text{N m}^{-2}}{\text{m}} = \text{m s}^{-2}\right] \tag{40a} \end{equation} Consequently, the Eötvös parameter $\eta(A, B) \equiv 2\frac{|\mathbf{a}_A - \mathbf{a}_B|}{|\mathbf{a}_A + \mathbf{a}_B|} \equiv 0$ vanishes identically, satisfying the Weak Equivalence Principle (WEP) and the sub-femto-g empirical bounds of the MICROSCOPE space experiment ($\eta(\text{Ti}, \text{Pt}) < 10^{-15}$) regardless of a body's mass, baryon-to-lepton ratio, or nuclear binding fraction.

4. Resolution of Soliton Boundary Scattering vs. Continuum Hydrostatics and Positive Hamiltonian Mass

Two foundational questions arise in any hydrodynamic theory of gravity: (a) If the carrier radiation reaction body force $\mathbf{f}_{\text{rad}} = +\nabla P$ holds the ambient substrate static ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), why does $\mathbf{f}_{\text{rad}}$ not act equally on a soliton formed of the substrate and cancel its gravitational acceleration to zero? (b) Since the soliton's Mach-1 sonic core has a static density below ambient ($\rho_{\text{sonic}} = 0.60653\,\rho_0 < \rho_0$, Eq. 20c), why does the soliton not experience outward Archimedes buoyancy (anti-gravity) in the outward-increasing pressure gradient $dP/dr > 0$? Both questions are resolved rigorously by the two-fluid boundary and Hamiltonian structure of a topological vortex soliton:

1. Continuous Carrier Hydrostatics vs. Localized Soliton Boundary Scattering: The radiation reaction body-force density $\mathbf{f}_{\text{rad}} = +\nabla P$ in Eq. (39b) represents the weak, coherent acoustic radiation pressure acting on the unknotted, irrotational ambient continuum ($\nabla \times \mathbf{u}_s = \mathbf{0}$, $|\mathbf{u}_{\text{substrate}}| = 0$), which transmits background carrier modes smoothly without topological phase singularities. A topological matter soliton, by contrast, is not a passive volume element of ambient fluid: it is a localized, quantized vortex singularity ($\oint_{\mathcal{C}} \mathbf{u}_s \cdot d\mathbf{l} = n \kappa_0 \ne 0$) enclosed by the parabolically degenerate Mach-1 sonic horizon ($|\mathbf{u}| = c_0$, Eq. 20b). Incident carrier phonons cannot pass undisturbed through the Mach-1 sonic horizon; instead, they undergo asymmetric refractive scattering around the closed circulating vortex loop (the acoustic analogue of the Aharonov–Bohm/Magnus wave-scattering cross section). Because the background acoustic energy density $u_{\text{ac}}(\mathbf{x}) = P(\mathbf{x})$ is higher on the exterior side ($r + \delta r$) than on the interior barycentric side ($r - \delta r$), asymmetric carrier phonon scattering across the closed sonic-horizon boundary imparts a net inward wave-momentum transfer equal to the Bernoulli–Bjerknes wave-gradient force $\mathbf{F}_{\text{grad}} = -(m_i/\rho_0)\nabla P$ (Eq. 39d). Thus, the continuum hydrostatic balance $-\nabla P + \mathbf{f}_{\text{rad}} = \mathbf{0}$ holds the unknotted background medium stationary while the localized vortex singularity translates freely down the pressure gradient.
2. Positive Integrated Two-Fluid Hamiltonian Energy ($m_i > 0$) vs. Core Cavitation Deficit: A static gas bubble in water floats upward toward lower pressure only if its total mass-energy density is less than the displaced ambient fluid ($\Delta M_{\text{eff}} < 0$). For a relativistic topological vortex soliton in the Quantum Plenum, however, the effective inertial and active/passive gravitational mass $m_i$ is governed by the total integrated two-fluid Hamiltonian energy functional over the soliton volume $\mathcal{V}_{\text{sol}}$: \begin{equation} m_i c_0^2 \equiv \int_{\mathcal{V}{\text{sol}}} \mathcal{H}{\text{sol}},d^3x = \int_{\mathcal{V}_{\text{sol}}} \left( \frac{1}{2}\rho_s |\mathbf{u}_s|^2 + \frac{1}{2}\rho_n |\mathbf{u}n|^2 + \Delta \mathcal{U}{\text{pot}}(\rho) \right) d^3x > 0 \tag{40b} \end{equation} where $\frac{1}{2}\rho_s |\mathbf{u}_s|^2 + \frac{1}{2}\rho_n |\mathbf{u}_n|^2$ is the strictly positive circulating kinetic energy density of the quantized superfluid and normal-fluid velocity fields (reaching $|\mathbf{u}| = c_0$ at the sonic horizon and sweeping $p \cdot q \, \pi^D$ standing-wave phase modes across $\mathbb{T}^5$, Eq. 12), and $\Delta \mathcal{U}_{\text{pot}}(\rho)$ is the local compressional potential energy relative to the ambient state. Because the quantized toroidal circulation kinetic energy integrated over the extended $1/r_{\perp}$ vortex velocity field vastly exceeds the shallow static core pressure deficit ($\Delta P_{\text{sonic}} = 0.39347\,P_0$ confined strictly to the thin sonic-horizon filament), the total integrated Hamiltonian mass perturbation $m_i$ is strictly positive ($m_i > 0$) for all super-threshold solitons ($E_k \ge E_{\text{cav}} \implies m_i \ge \mu_0 > 0$). By Eq. (39c), a localized wave-energy packet with strictly positive Hamiltonian mass $m_i > 0$ immersed in an outward-increasing pressure gradient ($dP/dr = +\rho_0 GM/r^2 > 0$) minimizes its total acoustic enthalpy $U_i(r) = m_i \Delta P(r)/\rho_0 = -G M m_i / r$ by accelerating radially inward toward the low-pressure center ($\mathbf{a} = -\frac{1}{\rho_0}\nabla P = -\frac{GM}{r^2}\hat{\mathbf{r}}$), guaranteeing universal gravitational attraction rather than outward buoyancy.

If an individual test soliton accelerates along the inward acceleration field (40a) from rest at infinity ($r \to \infty$, $v = 0$), its particle trajectory speed $v(r)$ satisfies the convective energy integral: \begin{equation} v\frac{dv}{dr} = -\frac{GM}{r^2} \implies \frac{1}{2}v^2(r) = \frac{GM}{r} \implies v(r) = \sqrt{\frac{2GM}{r}} \tag{40c} \end{equation} which is the kinematic velocity of the moving test soliton along its trajectory, not a bulk flow velocity of the static Plenum continuum ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$).

7.2 The Continuous Non-Linear Acoustic Pressure PDE & Flat Galactic Rotation Curves

In astrophysical disc galaxies (such as the 175 systems cataloged in the SPARC database [7]), observed orbital velocities at large radii $r \gg R_{\text{disc}}$ flatten to an asymptotically constant velocity plateau $v_{\text{rot}}(r) \to v_c = \text{const}$, defying the Keplerian falloff $v_{\text{Kep}}(r) = \sqrt{GM/r} \propto r^{-1/2}$ expected from visible baryonic mass alone. Rather than attributing this deficit to an unobserved cold dark matter halo or imposing an ad-hoc algebraic interpolation between regimes, TRT derives both solar-system and galactic dynamics from a single continuous non-linear acoustic pressure PDE (Postulate P6).

We first identify the characteristic cosmological and galactic acceleration scales of the Plenum. From the background expansion rate $H_{\text{early}}$ and normal-fluid mass fraction $\chi_n = \beta_{\text{TRT}}$, the Universal Plenum Acoustic Acceleration Floor ($a_0$) is: \begin{equation} a_0 \equiv c_0 H_{\text{early}} \beta_{\text{TRT}} = (6.5483 \times 10^{-10}\text{ m s}^{-2}) \times 0.05392055 = 3.531 \times 10^{-11}\text{ m s}^{-2} \tag{41} \end{equation} Along a thin planar disc galaxy embedded in the 3D spherical cavity, the acoustic shear boundary layer scales with the circular disc transverse boundary projection factor—the ratio of the disc circumference to its diameter: \begin{equation} \mathcal{C}{\text{disc}} \equiv \frac{2\pi R{\text{disc}}}{2 R_{\text{disc}}} = \pi \tag{41a} \end{equation} yielding the effective galactic disc boundary acceleration threshold declared under Postulate P6: \begin{equation} a_{\text{gal}} \equiv \pi a_0 = \pi (c_0 H_{\text{early}} \beta_{\text{TRT}}) = 1.109 \times 10^{-10}\text{ m s}^{-2} \approx 1.2 \times 10^{-10}\text{ m s}^{-2} \tag{42} \end{equation} which matches the empirical Milgrom acceleration constant $a_0^{\text{(MOND)}} \approx 1.2 \times 10^{-10}\text{ m s}^{-2}$ determined from SPARC rotation curves [7] to within $7.6\%$.

Formulation & Exact Asymptotic Integration of the Continuous Non-Linear Acoustic Pressure PDE

In a compressible two-fluid substrate where radial acoustic radiation pressure couples to the normal-fluid shear boundary layer at ultra-low pressure gradients ($|\nabla P|/\rho_0 \le a_{\text{gal}}$), the acoustic susceptibility flux is governed by the dimensionless gradient ratio $x \equiv \frac{|\nabla P|}{\rho_0 a_{\text{gal}}}$. Under Postulate P6, the static acoustic pressure field $P(\mathbf{x})$ across all gravitational regimes satisfies the continuous non-linear divergence-form acoustic pressure PDE: \begin{equation} \nabla \cdot \left[ \mu!\left(\frac{|\nabla P(\mathbf{x})|}{\rho_0 a_{\text{gal}}}\right) \nabla P(\mathbf{x}) \right] = 4\pi G \rho_0 , \rho_M(\mathbf{x}), \qquad \mu(x) \equiv \frac{x}{\sqrt{1 + x^2}}, \qquad \lim_{|\mathbf{x}|\to\infty} P(\mathbf{x}) = P_0 \tag{43} \end{equation} Because $\mu(x) = x(1+x^2)^{-1/2}$ is smooth, strictly positive, and monotonically increasing for all $x > 0$ (with $\frac{d}{dx}[x\mu(x)] = \frac{2x + x^3}{(1+x^2)^{3/2}} > 0$), Eq. (43) is a well-posed quasilinear elliptic PDE throughout $\mathbb{R}^3 \setminus \{\mathbf{0}\}$ that continuously connects the solar-system and deep-galactic regimes without piecewise stitching:

1. Newtonian / Solar-System High-Gradient Limit ($x \equiv \frac{|\nabla P|}{\rho_0 a_{\text{gal}}} \gg 1$): In the solar system and laboratory regime where $|\nabla P|/\rho_0 \gg a_{\text{gal}}$, the constitutive susceptibility function expands as $\mu(x) = (1 + x^{-2})^{-1/2} = 1 - \frac{1}{2}x^{-2} + \mathcal{O}(x^{-4}) \to 1$. Substituting $\mu(x) \to 1$ into Eq. (43) recovers the exact linear Poisson pressure-depletion law (Eq. 38): $\nabla^2 P(\mathbf{x}) = 4\pi G \rho_0 \, \rho_M(\mathbf{x}) \implies \frac{1}{\rho_0}\frac{dP}{dr} = \frac{GM}{r^2} \implies \mathbf{a}(r) = -\frac{1}{\rho_0}\nabla P = -\frac{GM}{r^2}\hat{\mathbf{r}}$
2. Deep-Galactic Normal-Fluid Boundary Shear Limit ($x \equiv \frac{|\nabla P|}{\rho_0 a_{\text{gal}}} \ll 1$): At large galactocentric radii $r \gg R_{\text{disc}}$ outside a galaxy of total baryonic mass $M_{\text{baryon}} = \int \rho_M(\mathbf{x})\,d^3x$ where $|\nabla P|/\rho_0 \ll a_{\text{gal}}$, $\mu(x) = x(1+x^2)^{-1/2} \to x = \frac{|\nabla P|}{\rho_0 a_{\text{gal}}}$, reducing Eq. (43) to the degenerate $p$-Laplacian ($p = 3$) acoustic flux equation: $\frac{1}{\rho_0 a_{\text{gal}}}\nabla \cdot \left( |\nabla P| \nabla P \right) = 4\pi G \rho_0 \, \rho_M(\mathbf{x})$ Integrating this divergence equation over a sphere $B_r$ of radius $r$ enclosing the baryonic mass $M_{\text{baryon}}$ via Gauss's divergence theorem gives: $\oint_{\partial B_r} \mu\!\left(\frac{|\nabla P|}{\rho_0 a_{\text{gal}}}\right) \nabla P \cdot d\mathbf{A} = 4\pi r^2 \left(\frac{|\nabla P(r)|}{\rho_0 a_{\text{gal}}}\right) |\nabla P(r)| = 4\pi G \rho_0 M_{\text{baryon}}$ Dividing both sides by $4\pi r^2 \rho_0$ and taking the positive square root yields the exact asymptotic deep-field pressure gradient and inward soliton acceleration (with $g_{\text{Newton}}(r) \equiv G M_{\text{baryon}}/r^2$): \begin{equation} \frac{1}{\rho_0}|\nabla P(r)| = \sqrt{\frac{G M_{\text{baryon}} a_{\text{gal}}}{r^2}} = \sqrt{g_{\text{Newton}}(r),a_{\text{gal}}} \implies \mathbf{a}(r) = -\frac{1}{\rho_0}\nabla P(r) = -\frac{\sqrt{G M_{\text{baryon}} a_{\text{gal}}}}{r}\hat{\mathbf{r}} \tag{43a} \end{equation} Integrating $dP/dr = \rho_0 \sqrt{G M_{\text{baryon}} a_{\text{gal}}} / r$ shows that the acoustic pressure transitions smoothly to a logarithmic normal-fluid shear profile $P(r) - P(R_0) = \rho_0 \sqrt{G M_{\text{baryon}} a_{\text{gal}}} \ln(r/R_0)$. Equating the inward gravitational acceleration $|\mathbf{a}(r)| = \sqrt{G M_{\text{baryon}} a_{\text{gal}}}/r$ from Eq. (43a) to the centripetal acceleration $v_{\text{rot}}^2 / r$ yields the exact asymptotically flat galactic rotation plateau and baryonic Tully–Fisher law: \begin{equation} \frac{v_{\text{rot}}^2}{r} = \frac{\sqrt{G M_{\text{baryon}} a_{\text{gal}}}}{r} \implies v_{\text{rot}}^4 = G M_{\text{baryon}} a_{\text{gal}} \implies v_c = \left(G M_{\text{baryon}} a_{\text{gal}}\right)^{1/4} = \text{const} \tag{44} \end{equation} Thus, both the inverse-square Newtonian solar-system law and the flat galactic rotation curve plateau emerge as exact asymptotic integrals of the single continuous PDE (43). (Full galaxy-by-galaxy SPARC rotation curve fitting and radial acceleration relation (RAR) residuals are formally developed in Paper 2 of this series.)

---

8. Conclusion & Master Classification Ledger of Claims and Predictions

By separating strictly dimensionless electro-hydrodynamic couplings ($A_{\text{RT}} = 50.412727$, $\beta_{\text{TRT}} = 0.05392055$) from dimensional SI Plenum scales (Table 2) and explicitly declaring all six constitutive postulates and discrete structural indices (Table 1), The Resonant Theory yields a self-contained mathematical and phenomenological architecture. To ensure complete clarity regarding epistemological status, Table 4 explicitly classifies every quantitative result in this paper into (A) Declared Constitutive Postulates & Empirical Anchors, (B) Retrodictive Benchmarks Against Known Constants, (C) Calibrated Solar-System Consistency Models, and (D) Prospective Out-of-Sample Experimental Falsification Tests.

Table 4: Master Classification Ledger of Constitutive Inputs, Retrodictive Benchmarks, and Prospective Predictions

---

9. References

1. L. Tisza, "Transport Phenomena in Helium II," Nature, vol. 141, p. 913, 1938. doi: 10.1038/141913a0
2. L. D. Landau, "Theory of the Superfluidity of Helium II," Physical Review, vol. 60, no. 4, pp. 356–358, 1941. doi: 10.1103/PhysRev.60.356
3. T. H. R. Skyrme, "A Non-Linear Field Theory," Proceedings of the Royal Society of London A, vol. 260, no. 1300, pp. 127–138, 1961. doi: 10.1098/rspa.1961.0018
4. W. G. Unruh, "Experimental Black-Hole Evaporation?" Physical Review Letters, vol. 46, no. 21, pp. 1351–1353, 1981. doi: 10.1103/PhysRevLett.46.1351
5. G. E. Volovik, The Universe in a Helium Droplet. Oxford University Press, 2003. doi: 10.1093/acprof:oso/9780199564842.001.0001
6. L. Faddeev and A. J. Niemi, "Stable Knot-Like Structures in Classical Field Theory," Nature, vol. 387, pp. 58–61, 1997. doi: 10.1038/387058a0
7. F. Lelli, S. S. McGaugh, and J. M. Schombert, "SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves," The Astronomical Journal, vol. 152, no. 6, p. 157, 2016. doi: 10.3847/0004-6256/152/6/157
8. Planck Collaboration (N. Aghanim et al.), "Planck 2018 Results. VI. Cosmological Parameters," Astronomy & Astrophysics, vol. 641, p. A6, 2020. doi: 10.1051/0004-6361/201833910
9. A. G. Riess et al., "A Comprehensive Measurement of the Local Value of the Hubble Constant with $1\text{ km s}^{-1}\text{Mpc}^{-1}$ Uncertainty from the Hubble Space Telescope and the SH0ES Team," The Astrophysical Journal Letters, vol. 934, no. 1, p. L7, 2022. doi: 10.3847/2041-8213/ac5c5b
10. S. G. Turyshev et al., "Support for the Thermal Origin of the Pioneer Anomaly," Physical Review Letters, vol. 108, no. 24, p. 241101, 2012. doi: 10.1103/PhysRevLett.108.241101
11. T. D. Moyer, "Formulation for Observed and Computed Values of Deep Space Network Data Types for Navigation," JPL Publication 00-7, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, 2000.