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

Author: Justin Time Miller
Version: 10.4
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}$, 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, galactic acceleration floor, 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: 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}$) under a single closed-form functional.
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 & Apparent Hubble Shift: 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}$ with zero coordinate ambiguity.
5. Static Hydrostatic Realization of Newtonian Gravity & Universal Free-Fall: In the static background Plenum ($\mathbf{u} = \mathbf{0}$, where the mass continuity equation $\partial_t\rho + \nabla\cdot(\rho\mathbf{u}) = 0$ is identically satisfied everywhere in the exterior vacuum), a Poisson pressure-depletion law $\nabla^2 P(\mathbf{x}) = 4\pi G\rho_0\rho_M(\mathbf{x})$ establishes an exterior harmonic depression $\nabla^2 P = 0 \implies P(r) = P_0 - \frac{\rho_0 GM}{r}$ held in exact static hydrostatic equilibrium ($\rho_0 \partial_t \mathbf{u} = -\nabla P + \mathbf{f}_{\text{rad}} = \mathbf{0}$) by acoustic carrier radiation reaction $\mathbf{f}_{\text{rad}} = +\nabla P = +\frac{\rho_0 GM}{r^2}\hat{\mathbf{r}}$ without fluid drainage. Topological solitons of inertial mass $m_i$ experience the wave-gradient radiation force $\mathbf{F}_{\text{grad}} = -\frac{m_i}{\rho_0}\nabla P$, yielding the universal free-fall acceleration $\mathbf{a} = \frac{\mathbf{F}_{\text{grad}}}{m_i} = -\frac{1}{\rho_0}\nabla P = -\frac{GM}{r^2}\hat{\mathbf{r}}$ strictly independent of body composition or density, identically satisfying the Weak Equivalence Principle and MICROSCOPE bounds.
6. Emergent Galactic Acceleration Threshold: The cosmological acoustic acceleration 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}} \equiv 2\pi R / (2R) = \pi$ evaluates the galactic 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}$ (within 7.6% of the empirical MOND scale), under which flat rotation curves and the baryonic Tully–Fisher scaling follow algebraically via a deep-field interpolation (with parameter-by-parameter SPARC galaxy dataset fitting formally scoped for Paper 2).
7. Exhaustive Decision Protocol & Epistemological Ledger: Table 4 compiles the master ledger of retrodictive benchmarks (including lunar recession spanning $[3.8143, 3.8157]\text{ cm yr}^{-1}$ between the exact quadratic and linear thin-sheath closures) and formalizes the 100% mathematically exhaustive four-regime decision protocol for the prospective 2028–2030 Europa Clipper Doppler test ($\bar{a}_{\text{plenum}} = (8.47 \pm 0.34) \times 10^{-10}\text{ m s}^{-2}$), classifying every real signed residual $a_{\text{res}} \in (-\infty, +\infty)$ without omission.

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 derives 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, maintaining the substrate in exact static hydrostatic equilibrium ($\mathbf{u} = \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, Void Fraction ($\epsilon_{\text{eff}}$) & The 1-Radian ($57.2958^\circ$) Limit

In a compressible fluid of density $\rho_0$, two spherical solitons pulsating in phase at angular frequency $\omega_0$ with volume amplitudes $\dot{V}_1, \dot{V}_2$ experience an attractive 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 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. Exact Evaluation of the Local Distance-Ladder Estimator ($H_{\text{late, TRT}} > H_{\text{early}}$): Early-universe CMB analysis [8] infers the background expansion rate $H_{\text{early}} = 67.40 \pm 0.50\text{ km s}^{-1}\text{Mpc}^{-1}$ by measuring the angular acoustic scale $\theta_* \equiv r_s(z_*) / D_A(z_*)$ at decoupling ($z_* \approx 1090$)—which is a geometric ratio of acoustic scales at recombination—whereas local distance-ladder surveys (SH0ES [9]) measure the local apparent flux-distance redshift slope $cz/d_L$, defining the late-universe Hubble estimator from the low-redshift linear slope $H_{\text{late}} \equiv c_0 \left.\frac{dz_{\text{obs}}}{dd_L}\right|_{d_L \to 0}$. Differentiating Eq. (21e) at $d_L = 0$ yields the strictly positive 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 a total torsional wind-up of $N_{\text{laps}} = 4,017.75 \times (24.4754^{-1} - 34.000^{-1}) = 46.0000\text{ laps}$, corresponding to the exact wind-up angle $\Delta\Phi_{\text{Sun}}$: \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 46.0000\text{ laps} = 16,560.0^\circ ; (289.0265\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,560.0^\circ - 15,850.0^\circ = 710.0^\circ$ ($12.3918\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{710.0^\circ}{50.412727 \times 1^\circ}\right)(1 \pm 0.0847) = 14.08 \pm 1.19,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.89\text{ to }15.27\,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 Spacecraft Impedance Force Model, Europa Clipper Error Budget & Exhaustive Decision Protocol

1. Constitutive Spacecraft Ram-Impedance Force Law & Exact Midpoint/Half-Width Band: A spacecraft of mass $m_{\text{sc}}$ moving radially outward at hyperbolic velocity $\mathbf{v} = v_r \hat{\mathbf{r}}$ ($v_r \ll c_0$) through the expanding Hubble-Plenum background experiences a cosmological acoustic impedance force $\mathbf{F}_{\text{plenum}} = m_{\text{sc}}\mathbf{a}_{\text{plenum}}$. Under the constitutive waveguide closure declared in Postulate P5, during one full carrier wave cycle ($\Delta\phi_{\text{cyc}} = 2\pi\text{ rad}$) the dimensionless cycle impedance is $Z_{\text{cyc}} = 2\pi + \beta_{\text{TRT}} = 6.33711$, while projecting the bulk susceptibility $A_{\text{RT}}$ across the two orthogonal half-wave toroidal cross-sections ($\pi \times \pi = \pi^2$) of a 3D dipole wake defines the reduced cross-sectional susceptibility $A_{\text{eff}} = A_{\text{RT}}/\pi^2 = 5.10788$. Multiplying the baseline Hubble expansion acceleration $-c_0 H_{\text{early}}\hat{\mathbf{r}}$ ($c_0 H_{\text{early}} = 6.5483 \times 10^{-10}\text{ m s}^{-2}$) by the ratio $Z_{\text{cyc}}/A_{\text{eff}}$ defines the inner-heliospheric uncompressed baseline deceleration $a_{\text{min}}$: \begin{equation} \mathbf{a}{\text{plenum}}^{(0)} = -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 & Emergent Galactic Acceleration Floor

Physical Picture: In a compressible fluid, what we experience as "gravity" in the weak-field Newtonian regime ($|\Phi|/c_0^2 = GM/(r c_0^2) \ll 1$) 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 absorbs coherent acoustic wave energy at its vortex core boundaries, it establishes a steady, static radiation-pressure depression $P(r) = P_0 - \frac{\rho_0 G M}{r}$ below the ambient far-field pressure $P_0$. Crucially, the background Quantum Plenum substrate itself remains in static hydrostatic equilibrium ($\mathbf{u}_{\text{substrate}} = \mathbf{0}$), so there is no continuous material drainage of the vacuum medium into the body. Because pressure increases outward from the low-pressure mass toward $P_0$ ($\frac{dP}{dr} = +\frac{\rho_0 G M}{r^2} > 0$), any test parcel or physical body immersed in that outward-increasing pressure field experiences an inward acoustic radiation buoyancy acceleration $\mathbf{g} = -\frac{1}{\rho_0}\nabla P = -\frac{GM}{r^2}\hat{\mathbf{r}}$ directed toward the barycentric low-pressure center.

7.1 Static Linearized Poisson Pressure-Depletion Law, Hydrostatic Momentum Balance & Universal Free-Fall

We formulate the gravitational sector explicitly in the weak-field linearized acoustic regime where $|\Phi(\mathbf{x})|/c_0^2 = GM/(r c_0^2) \ll 1$, so that 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), the undisturbed background Quantum Plenum substrate is static ($\mathbf{u}_0 = \mathbf{0}$, $\mathbf{j}_0 = \rho_s \mathbf{u}_s + \rho_n \mathbf{u}_n = \mathbf{0}$), and 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$ acts as a localized acoustic wave scatterer and absorption center that 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} = \mathbf{0}$), the steady mass continuity equation (Postulate P1, Eq. 1) satisfies: \begin{equation} \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 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} = \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 Plenum Substrate

In Section 2.1, the linearized momentum conservation equation of the two-fluid Plenum substrate 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}}{\partial t} = -\nabla P + \mathbf{f}{\text{rad}} \tag{39a} \end{equation} For the substrate to remain in exact static hydrostatic equilibrium ($\mathbf{u} = \mathbf{0}$, $\partial_t \mathbf{u} = \mathbf{0}$), Eq. (39a) requires exact mechanical force cancellation everywhere in the exterior vacuum: \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 substrate. 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}}$. The substrate is thus held in exact static hydrostatic equilibrium ($\mathbf{u} = \mathbf{0}$) without any fluid acceleration, inflow, or material drainage—in precise physical analogy to how outward radiation pressure balances inward gravitational pressure in a stellar atmosphere in hydrostatic balance ($-\nabla P + \mathbf{f}_{\text{rad}} = \mathbf{0}$).

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

We now derive the gravitational acceleration of a physical test body. In TRT, matter particles are not foreign, rigid objects displacing substrate fluid; by Postulate P3, elementary fermions and composite nuclei are topological vortex solitons formed of the Quantum Plenum substrate itself.

A localized soliton or composite body of total inertial rest mass $m_i$ represents an integrated substrate energy perturbation $E_i = m_i c_0^2 = \int_{\Omega_i} \mathcal{H}_{\text{sol}}\,d^3x$. In the presence of a background acoustic pressure variation $\Delta P(\mathbf{x})$, the local energy density of the acoustic mode couples to the background pressure field. In a barotropic acoustic medium where $P = c_0^2 \rho$, the acoustic enthalpy per unit mass is $w = \int \frac{dP}{\rho} \approx \frac{\Delta P}{\rho_0}$. A localized soliton of integrated mass $m_i$ possesses the interaction potential energy: \begin{equation} U_i(\mathbf{x}) = \int_{\Omega_i} \delta\rho(\mathbf{x}') w(\mathbf{x}'),d^3x' = m_i \frac{\Delta P(\mathbf{x})}{\rho_0} \tag{39c} \end{equation} Equivalently, under wave-action conservation and acoustic radiation theory, a localized wave packet or stationary topological soliton of wave energy $E_i = m_i c_0^2$ in a background acoustic field experiences the wave-gradient acoustic radiation force (the generalized Bjerknes radiation 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 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} Eq. (40a) establishes that:

1. The resulting acceleration is strictly attractive (directed radially inward toward the center of mass).
2. The acceleration is strictly universal: it is independent of the body's total mass $m_i$, constituent baryon-to-lepton ratio, nuclear binding energy fraction, or material density $\rho_{\text{body}}$.
3. 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}$) as an exact theorem of topological soliton dynamics in the Plenum.

4. Kinematic Free-Fall Speed as a Particle Trajectory Consequence

If an individual test parcel of matter accelerates along the inward acceleration field (40a) from rest at infinity ($r \to \infty$, $v = 0$), its particle velocity $v(r)$ along the trajectory satisfies the convective energy relation: \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{40b} \end{equation} This $\sqrt{2GM/r}$ is the kinematic velocity of an accelerating test particle along its trajectory, not a bulk material velocity of the static Plenum substrate, preserving $\mathbf{u}_{\text{substrate}} = \mathbf{0}$ and exact source-free continuity everywhere in the vacuum.

7.2 The Emergent Acoustic Acceleration Floor & 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(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. In standard cosmology, this deficit is attributed to an unobserved spherical cold dark matter halo.

In The Resonant Theory, a galaxy of total baryonic mass $M_{\text{baryon}}$ embedded in the cosmological two-fluid Plenum is subject to the background acoustic horizon condensation. When the local Newtonian gravitational acceleration drops below the characteristic cosmological Plenum acceleration scale $a_0$, the normal-fluid shear boundary layer transitions to a logarithmic acoustic pressure profile. We define the Universal Plenum Acoustic Acceleration Floor ($a_0$) directly from the background expansion rate $H_{\text{early}}$ and normal-fluid fraction $\beta_{\text{TRT}}$: \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} Eq. (42) evaluates to within 7.6% of 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]. Under Postulate P6, when the local Newtonian gravitational acceleration drops far below the galactic boundary threshold ($g_{\text{Newton}} \ll a_{\text{gal}}$), the acoustic pressure gradient transitions from the linear Poisson regime to the nonlinear acoustic shear boundary-layer regime via the standard deep-field interpolation $g(r) \to \sqrt{g_{\text{Newton}} a_{\text{gal}}} = \frac{\sqrt{G M_{\text{baryon}} a_{\text{gal}}}}{r}$. Equating this inward acceleration to the centripetal acceleration $v^2/r$ recovers the baryonic Tully–Fisher relation algebraically: \begin{equation} \frac{v_c^2}{r} = \frac{\sqrt{G M_{\text{baryon}} a_{\text{gal}}}}{r} \implies v_c = \left(G M_{\text{baryon}} a_{\text{gal}}\right)^{1/4} = \text{const} \tag{43} \end{equation} demonstrating that the asymptotically flat rotation curve plateau and baryonic Tully–Fisher scaling follow algebraically from the Plenum's boundary threshold. (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.