Appendix B: Reproducibility map

Every quantitative claim is reproduced by the companion package: thirty-four standalone Python scripts and two real datasets, the SPARC rotation curve of NGC 2403 and the Pantheon+ supernovae. The scripts use symplectic integrators, fixed and un-tuned constants except where a single fit is noted, and no random seeds. The vacuum-dispersion tension is made reproducible, not hidden.

Every quantitative claim is reproduced by the companion package: thirty-four standalone Python scripts and two real datasets (the SPARC rotation curve of NGC 2403 and the Pantheon+ supernovae). The scripts use symplectic integrators, fixed (un-tuned) constants except where a single fit is explicitly noted, and no random seeds for the deterministic results.

Every quantitative claim is reproduced by the companion package (VP_reproducibility): thirty-four standalone Python scripts and two real datasets (SPARC NGC2403_rotmod.dat; Pantheon+ Pantheon+_extract.tsv). The scripts use symplectic integrators, fixed (un-tuned) constants except where a single fit is explicitly noted, and no random seeds for the deterministic results. The vacuum-dispersion tension (Chapter 2) is made reproducible (the transverse-mode estimate, the §10.9.1 angle reclassification, the collective-disturbance mechanism, the Goldstone/Maxwell synchronization account, the collective-stiffness ceiling with the gamma residual at scale a, and the relativistic energy–velocity law), not hidden.

ChapterScriptHeadline expected output
1ch1_inflow_rates.pyν_H=293.227; per-body Q table
2ch2_light.pyc∝√(K) (R²=1); E_(QG) vs Fermi (8–15 orders); fcc coord. 12, isotropic
2ch2_lightangle.pyangle table sinχ=λ/(mD) at canon D: visible 550nm→89.90^(∘); pre-registered pair 632.99/532→89.9378^(∘)/89.8248^(∘); gamma 1pm→11.9^(∘), 1fm→0.012^(∘) (quasi-longitudinal)
2ch2_grb.pycollective disturbance vs high-k photon: indep. packet spreads ×11, collision pulse coherent ×1.0 (dispersionless, toy)
2ch2_goldstone.pysynced rotors lock (r→0.96); Goldstone mode ω=2√(J/I)|sin(k/2)|: dispersionless at small k, curves at high k
2ch2_stiffness.pyceiling v_g≤c (collective stiffness); gamma residual at scale a: E_(rm QG,2)≈882 GeV, 8 orders below Fermi
2ch2_relativistic.pytwo laws: relativistic Δ v/c=tfrac12(E₀/E)² (shrinks, v_g≤c) vs phonon (E/E_(rm QG))²; massless ⇒ dispersionless
2ch2_gammashot.py2D launch: gamma packet disperses (peak ≈0.43c) vs visible (≈ c); gamma spreads LESS transversely (no 90^∘ fan-out)
2ch2_collision.pycollision test: speed is amplitude-independent (stiffness; gentle=violent) but wavelength-dependent (ccos(ka/2)); isolates the open item to the wavelength axis
2ch2_shake.pyshake test: a long-wave vibration relays at c and stays coherent (×1.1); a gamma-wavelength vibration disperses (v_g) and smears (×16), so a literal gamma packet cannot survive cosmic distance
2ch2_gammacontent.pygamma wavelengths live inside lattice waves: a churning disturbance is broadband (reaches the gamma band, amplitude-independent fraction); burst energy V/a³ quanta ( 1m³ holds one burst)
2ch2_burstprop.pydecisive test: the gamma (high-k) content separates from the low-k (GW-like) front at v_g=ccos(ka/2), weak=strong—a burst does not carry its gamma coherently
2ch2_jammed2d.pydensity=1 (jammed/no-void) freezes the longitudinal channel but leaves the transverse (light) dispersion cos(ka/2) unchanged: jamming rescales c, not the deficit
2ch2_upconvert.pyup-conversion test: a gentle (low-strain) wave cannot shake the lattice to the gamma (high-k) scale; the threshold is per-bond strain 1, the same in 1D and 3D, so a wave with strain h is 20 orders short—the local-production escape fails too
3ch3_gravity.py1/r² bounded vs 1/r⁵ runaway; EP 1.9×10⁻¹⁴; mass–G degeneracy
3ch3_gate.pysaturation g(eₐ)→1 (g_*); choke radius r_(ch)=√(K_F/Bc)=1; |S|≤ceₐ
4ch4_solar_system.pyperiods ≤0.73%; T²/a³=1.00001
5ch5_spin_tidal.pyspin-up ω→ v_(sw)/a; tidal lock ω/n:4→1; Δ g
6ch6_galaxy_rar.pya₀=1.08×10⁻¹⁰; NGC 2403 Υ=0.567, χ²/dof=1.99
7ch7_lattice_optics.pyPantheon+ χ²/dof=0.50 vs ΛCDM 0.44; z_(min)=1.72 vs 1.61
8ch8_deficit.pydM/dr=4π (flat); core ρ=0 (absolute zero); rays bend
8ch8_bullet.pycluster lensing/gas offset, 0.15–0.86L_(off) (onset τ_Δ τ_(coll))
9ch9_lattice_cmb.pylangleKE⟩/langlePE⟩=0.998; slope→ c; n_(max) 4×10⁵⁴; L 6×10⁴⁵ erg/s
5ch_galactic_spin.pygalactic tide on planets 10⁻¹⁸–10⁻¹³ of Sun's; 5e16× to matter at Earth ⇒ negligible (spins accommodated)
10ch10_post_newtonian.pylight bending 1.751” (→Eddington 1.75); Mercury perihelion Newton 0.0 vs PPN(γ=β=1) 42.98”/cy (degenerate; Yoshida-4 cross-check +42.4”)
11ch11_grav_waves.pylattice pulse speed 1.000c; chirp (M_c=28M_odot) 35→250 Hz in 0.19s; Hulse–Taylor dot P_b=-2.40×10⁻¹² (degenerate)
12ch12_puzzles.pyOlbers integral diverges without attenuation, → nL/κₒₚₜ (finite) with κₒₚₜ=H₀/c ⇒ dark sky (distinguishing)
13ch13_bbn.pyfreeze-out n/p1/6→1/7, Yₚ=0.250 (target met); stellar D/H 2.5→1.0→0.4×10⁻⁵ destroyed, wrong way (conflicting)
14ch14_cmb_aniso.pyfirst peak ℓ=220⇒θ=π/220=0.82^∘; single-scale toy gives one bump (harmonic ℓ 540/810 need a coherent mechanism — see sim_acoustic_peaks) (\textsf{HYP}, open)
15ch15_lss_bao.pybroadband ξ(r) reproduced; 150Mpc BAO bump from a fixed inflow length (sim_acoustic_peaks); deriving that length is the open input (\textsf{HYP}, open)
14/15sim_acoustic_peaks.py(1) harmonic series needs phase coherence (1:2:3 from a synchronized clock; random phase → flat); (2) one fixed inflow-shell length Rₛ gives the BAO bump at Rₛ and P(k) zeros at nπ/Rₛ; (3) random driving → smooth P(k) (corr 0.99), peaks need a resonant cavity. Open program: derive L≈150Mpc (\textsf{HYP})
Ech12_sunspot.pyinflow-sink: B_z100→ 3000 G; T 3800 K; Wilson 300 km (\textsf{HYP}/\textsf{SPEC}, degenerate)
16ch10_ledger.pyhonest-ledger synthesis figure: per-chapter status (solid/degenerate/distinguishing/tension) + three boxed predictions (renderer, no new sim)
16ch_accuracy_dashboard.pyper-scale accuracy + input tags (M/F/D): a₀ 90% (distinguishing); Kepler/Moon/Mercury degenerate; Venus accommodated; galactic tide negligible