Supplementary Material: Formal Verification, Asymptotic PDE Reductions, and Lattice Divergence Analysis

Companion Technical & Formal Supplement

Comment: On the Foundations of Canonical Field Quantization, Relativistic Reductions, and the Invariance of Operator Ladders
Author: R. Fisher, Braid Dynamics Research Group (ORCID: 0009-0006-2441-3282)
Contents: Section 1 (Verified Lean 4 Formal Kernel Specifications, 21 theorems) · Sections 2–4 (Wavepacket PDE, Gauge Propagator, Cauchy Continuum Limit) · Section 5 (Hardware-Optimized C++20 Lattice UV Engine) · Section 6 (Symbolic Constraint Algebra Audit) · Section 7 (Automated Verification Harness & Replication Guide)
Downloads & Code: Formal Lean 4 Proof · Supplementary Markdown · Full Replication Bundle

Overview

This supplementary material provides the complete, self-contained, machine-checked mathematical formalizations, high-precision numerical PDE simulations, symbolic algebra audits, and lattice quantum field divergence engines accompanying the comment paper.

The software architecture consists of six interoperable layers:

1. Formal Verification Kernel (Lean 4): Machine-checked proofs across four foundational modules ($0$ axioms, $0$ sorry across 21 theorems):
   • Module 1 (Coulomb Graph Invariants): Complete graph pair-sum invariants over arbitrary additive abelian groups supporting attractive Coulomb interactions ($S_1 = 2U$, $S_1 = V_{\text{nrel}} + U_{ee}$).
   • Module 2 (Clifford/Dirac Projector Algebra): Dimension-independent Clifford and Dirac operator projector algebras parameterized over an abstract commutative ring, with concrete matrix models and normalized projector idempotency ($\Pi_\pm^2 = \Pi_\pm$).
   • Module 3 (Infinite Heisenberg-Weyl Ladder): The infinite-dimensional Heisenberg-Weyl Fock vacuum ladder theorem eliminating $\mathfrak{su}(2)$ termination via a certified sequence space representation ($(a^\dagger)^n |0\rangle \ne 0$).
   • Module 4 (Cauchy Runaway Instability): Positive-semidefiniteness of the 4D Euclidean norm ($p_E^2 \ge 0$) and unconditional tachyonic Cauchy runaway instability from rest ($v_{m+1} \ge 2^m$).
2. Asymptotic PDE Wavepacket Engine (Python): High-precision spectral PDE integration of the 1D Klein–Gordon field reducing to the non-relativistic Schrödinger equation as $c \to \infty$, verifying the exact $\mathcal{O}(c^{-2})$ power-law error convergence.
3. Causal Gauge Radiation Field Solver (Python): Comparative spatial and energy flux simulation contrasting the standard causal retarded Green's function with Riverfield's null-Hamiltonian contact propagator ($\hat{H}_{\mathrm{EM}} \equiv 0$).
4. Continuum Cauchy Runaway Solver (Python): High-precision integration of the continuous tachyonic mode equation $\ddot{x} = \mu^2 x$ ($x(t) = \cosh(\mu t)$) and central-difference discretization across varying $\Delta t \in [10^{-3}, 0.2]$, proving that the exponential runaway is an intrinsic continuum property ($\lambda_+ > 1 + \sqrt{K} > 1$) and converges as $\Delta t \to 0$.
5. High-Performance Lattice UV Engine (C++20): Hardware-optimized Brillouin zone summation comparing 3D equal-time spatial foliation mode sums with 4D covariant Euclidean functional determinants, measuring identical leading quartic ultraviolet scaling ($\sim \Lambda^4$).
6. Symbolic Constraint Algebroid Audit (Python / SymPy): Verification of constant structure constants in the flat Minkowski Poincaré Lie algebra vs. dynamical metric structure functions in the Dirac hypersurface deformation algebra.

---

Table of Contents

1. Section 1: Machine-Checked Formal Verification in Lean 4
2. Section 2: Asymptotic Wavepacket Reduction & PDE Simulations
3. Section 3: Gauge Propagator & Radiation Extinction Simulation
4. Section 4: Continuum Cauchy Tachyonic Runaway Simulation
5. Section 5: Ultraviolet Lattice Scaling Engine (C++20)
6. Section 6: Symbolic Constraint Algebra Audit
7. Section 7: Automated Verification Harness & Reproduction Guide

---

Section 1: Machine-Checked Formal Verification in Lean 4

Source File: code/lean/RiverfieldRefutation.lean
Lean 4 Toolchain: leanprover/lean4:v4.33.1 (Compiles with 0 axioms, 0 sorry, and 0 warnings across 21 theorems).

Methodological Scope

This formal verification layer addresses exact algebraic, combinatorial, and spectral theorems. The Lean 4 development formalizes the kinematic and structural cores of the physical arguments—exact operator algebras, spectral decompositions, graph combinatorics, and discrete Cauchy recurrences—avoiding unnecessary axiomatic overhead. Continuous PDE dynamics and functional determinants are treated in the companion Python and C++ suites.

It refutes four foundational claims in Riverfield (2026) across four self-contained modules, with all typeclasses explicitly instantiated by concrete mathematical models:

Module 1: Pair Potentials & Complete-Graph Invariants in AddCommGroup

• Theorems 1–3 (directed_equals_two_undirected, total_atomic_sum_decomposition, riverfield_atomic_difference_identity): Parameterized over an arbitrary additive abelian group AddCommGroup G, enabling the formal representation of negative potential energies ($V < 0$, necessary for attractive electron-nucleus Coulomb interactions). Proves that for any symmetric, irreflexive pairwise interaction kernel, the sum of single-particle potentials identically equals twice the undirected potential sum ($S_1 = 2U$, $S_1 = V_{\text{nrel}} + U_{ee}$). Subtracting the physical non-relativistic atomic potential $V_{\text{nrel}}$ yields $S_1 - V_{\text{nrel}} = U_{ee}$ with zero relativistic or physical input. This mathematically demonstrates that Riverfield's putative "relativistic Foldy–Wouthuysen correction" $\Delta H_{\text{rel}}$ in Rel_atom.pdf (p. 2) is an elementary algebraic artifact of counting directed edges on complete graphs ($|E(K_N)| = N(N-1)/2$).

Module 2: Dimension-Independent Clifford/Dirac Projector Algebra & Concrete Model

• Model Verification (CommRing Int & DiracAlgebra Mat2 Int): Instantiates the abstract DiracAlgebra A R over a concrete 2x2 matrix model Mat2 over the integers Int, certifying with 0 axioms and 0 sorry that the algebraic structure is completely consistent.
• Theorems 4–9 (projector_completeness, projector_orthogonality, projector_pos_eigenvalue, projector_neg_eigenvalue, projector_pos_idempotent, projector_neg_idempotent): Parameterizes the Clifford dispersion property $H^2 = \omega^2 \cdot 1$ over an arbitrary commutative ring $R$ (with ring nontriviality certified by helper lemma commring_nontrivial), supporting arbitrary spacetime dimensions (1D, 3+1D $\mathrm{Cl}_{3,1}$, or $n$D). Proves resolution of identity $P_+ + P_- = (2\omega) \cdot 1$, mutual orthogonality $P_+ P_- = 0$, spectral eigenvalue selection ($H P_\pm = \pm \omega P_\pm$), and scaled projector idempotency ($P_\pm^2 = (2\omega) P_\pm$).
• Theorems 10–13 (pi_pos_idempotent, pi_neg_idempotent, pi_completeness, pi_orthogonality): Under the scalar invertibility condition $2\omega \in R^\times$ (two_omega_inv * (omega + omega) = 1), defines normalized projectors $\Pi_\pm = \frac{1}{2\omega} P_\pm$ and proves true idempotency $\Pi_\pm^2 = \Pi_\pm$, resolution of identity $\Pi_+ + \Pi_- = 1$, and mutual orthogonality $\Pi_+ \Pi_- = 0$. This certifies that the Cauchy data space invariantly decomposes as $\mathcal{H} = \mathcal{H}^+ \oplus \mathcal{H}^-$ across all spacetime dimensions.

Module 3: Heisenberg-Weyl Non-Nilpotence and Infinite Ladder on a Vacuum Module

• Model Verification (FockSpace (Nat → Int)): Constructs an explicit representation module on the infinite sequence space Nat → Int with creation shift $a^\dagger f = (0, f_0, f_1, \dots)$, annihilation operator $(a f)_n = (n+1) f_{n+1}$, and vacuum $|0\rangle = (1, 0, 0, \dots)$, proving $[a, a^\dagger] f = f$ and $a |0\rangle = 0$.
• Theorems 14–17 (a_state_succ, a_pow_state, state_non_zero, number_eigenvalue): Proves lowering action $a |n+1\rangle = (n+1) |n\rangle$, exact factorial projection $a^n |n\rangle = n! |0\rangle$, and strictly non-vanishing Fock states $(a^\dagger)^n |0\rangle \neq 0$ for all $n \in \mathbb{N}$. This completely closes the $\mathfrak{su}(2)$ termination loophole: in any finite-dimensional spin-$j$ representation, $(J_+)^{2j+1} = 0$, whereas the CCR with vacuum annihilation forces an infinite sequence of distinct eigenstates $N |n\rangle = n |n\rangle$.

Module 4: Euclidean Metric Obstruction & Unconditional Cauchy Runaway from Rest

• Theorems 18–19 (euclidean_norm_sq_nonneg supported by helper int_sq_nonneg, euclidean_mass_shell_no_solution): Proves positive semi-definiteness of the 4D Euclidean norm ($p_E^2 \ge 0$) and the obstruction to real Euclidean on-shell solutions ($p_E^2 \ne -m^2$ for $m^2 > 0$).
• Theorems 20–21 (cauchy_runaway_induction, cauchy_exponential_unbounded): Discretizes the tachyonic mode equation $\ddot{x} = K x$ ($K \ge 1$) and proves by induction that discrete velocity grows exponentially as $v_{m+1} \ge (1+K)^m K \ge 2^m$ for all $m \ge 0$ and all initial separations $\gamma \ge 0$. Proves that for ANY threshold $B \in \mathbb{N}$, the velocity strictly exceeds $B$ at step $n = B + 1$ ($v_{B+1} > B$) even when starting completely from rest ($\gamma = 0$), establishing unconditional dynamical ill-posedness.

Verification Certificate

lake env lean papers-drafts/cases/js-riverfield/code/lean/RiverfieldRefutation.lean
[Exit Code: 0] (21 Theorems Verified, 0 Axioms, 0 Sorry, 0 Warnings)

---

Section 2: Asymptotic Wavepacket Reduction & PDE Simulations

Source File: code/python/wavepacket_convergence.py
Generated Artifacts: figures/wavepacket_convergence.pdf, figures/wavepacket_convergence.png

Theoretical Foundation: 2nd-Order Hyperbolic Cauchy Problem

The simulation solves the exact second-order hyperbolic Klein–Gordon Cauchy initial value problem: $\frac{1}{c^2}\frac{\partial^2 \phi_k}{\partial t^2} + (k^2 + m^2 c^2)\phi_k = 0 \implies \ddot{\phi}_k(t) + \Omega_k^2 \phi_k(t) = 0$ where $\Omega_k = c\sqrt{k^2 + m^2 c^2}$. The physical positive-energy sector in the non-relativistic regime is prepared with Foldy–Wouthuysen Cauchy data: $\phi_k(0) = \psi_0(k), \quad \dot{\phi}_k(0) = -i mc^2 \psi_0(k)$ The exact analytic solution decomposes into forward-propagating ($A_k$) and backward-propagating antiparticle ($B_k$) modes: $A(k) = \frac{1}{2}\left(1 + \frac{mc^2}{\Omega_k}\right)\psi_0(k), \quad B(k) = \frac{1}{2}\left(1 - \frac{mc^2}{\Omega_k}\right)\psi_0(k)$ Expanding $\frac{mc^2}{\Omega_k} = 1 - \frac{k^2}{2m^2 c^2} + \mathcal{O}(c^{-4})$ analytically proves that the backward antiparticle amplitude is suppressed as $B(k) \sim \mathcal{O}(c^{-2})\psi_0(k)$, yielding an integrated antiparticle power suppression of $\mathcal{O}(c^{-4})$.

Simulation Results & Unit Test Verification

A Gaussian wavepacket $\psi_0(x)$ ($k_0 = 2.0, \sigma = 1.0$) was evolved across speed-of-light parameters $c \in [2, 500]$ over a spatial grid $L = 60.0$ ($N = 2048$):

[Wavepacket 2nd-Order Cauchy Convergence Output]
  Envelope L² Error Slope:    -1.8996 (Expected: -2.0000, R² = 0.997372)
  Antiparticle Power Slope:   -3.8236 (Expected: -4.0000, R² = 0.998083)

The automated test suite (test_suite.py) enforces both:

1. Envelope convergence: $-2.15 \le \text{slope}_{\text{env}} \le -1.85$ ($R^2 \ge 0.99$).
2. Antiparticle power suppression: $-4.20 \le \text{slope}_{\text{anti}} \le -3.80$ ($R^2 \ge 0.99$).

This demonstrates that the $\mathcal{O}(c^{-2})$ Schrödinger convergence is an intrinsic property of the physical Foldy–Wouthuysen positive-energy sector, accompanied by the dynamical decoupling of antiparticle degrees of freedom at rate $\mathcal{O}(c^{-4})$.

---

Section 3: Gauge Propagator & Radiation Extinction Simulation

Source File: code/python/gauge_propagator_simulation.py
Generated Artifacts: figures/propagator_comparison.pdf, figures/propagator_comparison.png

Gauge Field Dynamics and Continuum Green's Function Analysis

In Quad_Dirac.pdf (p. 21–22), Riverfield evaluates the free electromagnetic Hamiltonian density on-shell ($k^2 = 0$, $k^\mu A_\mu = 0$) and asserts that $\hat{H}_{\mathrm{EM}} \equiv 0$. Because the Hamiltonian generator vanishes, Heisenberg evolution freezes ($\partial_\tau A_\mu = 0$), collapsing the gauge two-point function to a local contact delta distribution: $\langle 0 | A_\nu(x') A_\mu^*(x) | 0 \rangle \propto \eta_{\mu\nu} \delta^{(4)}(x' - x)$

In linear field theory, the classical gauge potential generated by a localized source current $J_\mu(x)$ is determined by the Green's function convolution $A_\mu(x) = \int d^4x' D(x - x') J_\mu(x')$.

• Under the standard causal retarded Green's function $D_{\mathrm{ret}}(x - x') = \frac{\delta(t - t' - r/c)}{4\pi r}$, radiation fields propagate outward into the vacuum at speed $c$, establishing non-zero field strengths and non-zero Poynting flux on any bounding sphere $r > 0$.
• Under Riverfield's contact propagator $D_{\mathrm{RF}}(x - x') \propto \delta^{(4)}(x - x')$: $A_\mu(x) \propto \int d^4x' \, \delta^{(4)}(x - x') J_\mu(x') = J_\mu(x)$ The gauge field exists exclusively inside the current source $J_\mu(x)$. Outside the localized source where $J_\mu(x) \equiv 0$, the gauge field $A_\mu(x)$ is identically zero for all $r > 0$. Consequently, on any enclosing surface $r \ge r_{\text{cut}} > 0$, the electromagnetic field vanishes identically, Poynting flux $\mathbf{S} \equiv 0$, and radiation is strictly impossible. Modeling Riverfield's field as localized inside the source cell is the direct numerical realization of this continuum delta-distribution convolution.

Simulation Comparison & Spherical Power Integration

On a discrete 2D simulation grid ($N = 300 \times 300$, $L = 10.0$), the electromagnetic fields were evaluated from the vector potential $A_z$:

1. Smooth Core Regularization: To prevent discrete boundary ring artifacts when taking numerical derivatives near the origin, the source is mollified by a smooth $C^\infty$ cutoff $f(R) = 1 - \exp(-(R/r_{\text{cut}})^4)$ ($r_{\text{cut}} = 0.25$). This guarantees that all spatial derivatives remain smooth everywhere on the grid: $A_z(r) = \frac{\mu_0 I_0}{4\pi r} e^{i k r} \left(1 - e^{-(r/r_{\text{cut}})^4}\right)$
2. The magnetic field $\mathbf{B}$ was computed via discrete numerical curl: $\mathbf{B} = \nabla \times (A_z \hat{\mathbf{z}}) = (\partial_y A_z)\hat{\mathbf{x}} - (\partial_x A_z)\hat{\mathbf{y}}$.
3. The time-averaged Poynting flux density $\langle \mathbf{S} \rangle = \frac{1}{2\mu_0}\operatorname{Re}(\mathbf{E} \times \mathbf{B}^*)$ was computed across the entire grid.
4. Spherical Area Integration: For an oscillating electric dipole along $\hat{\mathbf{z}}$, the radial flux has equatorial maximum $\langle S_r(r, \theta) \rangle = S_r(\text{equator}) \sin^2\theta$. Integrating over the full 3D bounding sphere with area measure $da = r^2 \sin\theta d\theta d\phi$ yields: $P(r) = \oint_{\mathcal{S}_r} \mathbf{S} \cdot d\mathbf{a} = r^2 \int_0^{2\pi} d\phi \int_0^\pi \sin^3\theta d\theta \, S_r(\text{equator}) = \frac{4}{3} r^2 \oint_{\text{equator}} S_r(r, \theta) \, d\theta$ Because $S_r \propto 1/r^2$, the geometric factor $r^2$ cancels the radial decay, rendering total radiated power strictly constant across distance.
5. Results:
   • Standard Causal QED Propagator: Evaluates to a strictly constant radiated power flux $P(r) = 0.2377$ across all radii $r \in [0.6, 8.0]$, matching the exact analytical dipole theory $P_{\text{theory}} = \frac{\omega^2}{12\pi c} = 0.2387$ to within $0.43\%$.
   • Riverfield Contact Propagator: Because fields vanish identically outside the origin ($r > 0 \implies \mathbf{E} \equiv 0, \mathbf{B} \equiv 0$), the discrete numerical contour integral evaluates to exact numerical machine zero ($0.0 \pm 10^{-16}$) across every radius.

This numerical simulation proves that Riverfield's model abolishes light propagation, eliminates radiation, and destroys the static $1/r$ Coulomb law.

---

Section 4: Continuum Cauchy Tachyonic Runaway Simulation

Source File: code/python/tachyonic_cauchy_continuum.py
Generated Artifacts: figures/tachyonic_cauchy_runaway.pdf, figures/tachyonic_cauchy_runaway.png

Mathematical Foundation: Characteristic Multipliers & Continuum Limit

While Lean 4 proves an exact exponential lower bound for integer couplings $K \in \mathbb{Z}^+$, continuous hyperbolic PDEs $\ddot{x} = \mu^2 x$ are governed by dimensionless parameters $K = (\mu \Delta t)^2$. The exact discrete characteristic equation for central differencing: $\lambda^2 - (2 + K)\lambda + 1 = 0 \implies \lambda_+ = 1 + \frac{K}{2} + \sqrt{K + \frac{K^2}{4}}$ satisfies $\lambda_+ > 1 + \sqrt{K} = 1 + \mu \Delta t > 1$ for all $K > 0$, regardless of how small $\Delta t$ is. In the continuum limit $\Delta t \to 0$: $\lambda_+(\Delta t) = 1 + \mu \Delta t + \frac{1}{2}(\mu \Delta t)^2 + \mathcal{O}(\Delta t^3) = \exp(\mu \Delta t) + \mathcal{O}(\Delta t^3)$ $\lim_{\Delta t \to 0} \lambda_+(\Delta t)^{t / \Delta t} = \exp(\mu t)$ The continuous solution from rest ($x(0) = 1, \dot{x}(0) = 0$) is $x(t) = \cosh(\mu t) \sim \frac{1}{2} e^{\mu t}$.

Simulation Results & Tight Convergence Bounds

The continuum solver swept time steps $\Delta t \in [10^{-3}, 0.2]$ over $t \in [0, 5.0]$ with growth rate $\mu = 1.5$:

[Tachyonic Cauchy Continuum Output]
  dt = 0.2000: K = 0.090000, lambda_+ = 1.3503, eff_rate = 1.5015 (mu = 1.5000)
  dt = 0.1000: K = 0.022500, lambda_+ = 1.1623, eff_rate = 1.5042 (mu = 1.5000)
  dt = 0.0500: K = 0.005625, lambda_+ = 1.0780, eff_rate = 1.5010 (mu = 1.5000)
  dt = 0.0100: K = 0.000225, lambda_+ = 1.0151, eff_rate = 1.5000 (mu = 1.5000)
  dt = 0.0010: K = 0.000002, lambda_+ = 1.0015, eff_rate = 1.5000 (mu = 1.5000)

The test suite enforces:

1. Unconditional instability: $\lambda_+ > 1.0$ for all $\Delta t$.
2. High-precision rate matching: $|\text{eff\_rate} - \mu| / \mu < 10^{-4}$ at $\Delta t = 10^{-3}$ (empirically $9.37 \times 10^{-8}$).
3. Second-order error scaling: $70.0 \le \text{err\_dt}_{01} / \text{err\_dt}_{001} \le 130.0$ (reducing $\Delta t$ by $10\times$ decreases error by $99.40\times$).

This numerical analysis confirms that the exponential runaway is an intrinsic property of the tachyonic sign, completely independent of the integer cutoff $K \ge 1$ used in the Lean formalization.

---

Section 5: Ultraviolet Lattice Scaling Engine (C++20)

Source File: code/cpp/lattice_uv_divergence.cpp
Benchmark Data: code/cpp/uv_scaling_results.csv
CLI Invocation: lattice_uv_divergence [output_csv_path] (defaults to uv_scaling_results.csv if omitted)
Plotting Script: code/python/plot_uv_scaling.py
Generated Artifacts: figures/uv_scaling_comparison.pdf, figures/uv_scaling_comparison.png

Computational Methodology and Finite-Volume Scaling

To test Riverfield's assertion that vacuum energy divergences are artifacts of 3D spatial foliation, the engine computes vacuum energy densities across discretized Brillouin zones as a function of grid cutoff $\Lambda = \pi / a$:

1. 3D Equal-Time Lattice Mode Sum: $\mathcal{E}_{3D}(\Lambda) = \frac{1}{V_3} \sum_{\mathbf{k}} \frac{1}{2}\sqrt{\sum_{i=1}^3 \left(\frac{2}{a}\sin\frac{k_i a}{2}\right)^2 + m^2}$
2. 4D Covariant Euclidean Hypercubic Functional Determinant: $\mathcal{E}_{4D}(\Lambda) = \frac{1}{2 V_4} \sum_{k_E} \ln\left(\sum_{\mu=1}^4 \left(\frac{2}{a}\sin\frac{k_\mu a}{2}\right)^2 + m^2\right)$

Infrared Box Size Truncation Defense ($m L \ge 7.2 \gg 1$): Because grid size $N$ is fixed ($N_{3D} = 64, N_{4D} = 36$) while lattice spacing varies $a \in [0.8, 0.2]$, the physical box length is $L = N a$. For the smallest box at $a = 0.2$, $L = 36 \times 0.2 = 7.2$. With particle mass $m = 1.0$ (Compton wavelength $\lambda_C = 1/m = 1.0$), the dimensionless volume parameter satisfies $m L \ge 7.2 \gg 1$. Consequently, finite-volume infrared boundary corrections are exponentially suppressed: $e^{-m L} \le e^{-7.2} \approx 7.4 \times 10^{-4} < 0.1\%$ This mathematically guarantees that the measured power-law scaling is driven purely by ultraviolet mode accumulation rather than infrared box volume variation.

Continuous 4D Asymptotics & The $\ln \Lambda$ Exponent

In 4D Euclidean spacetime, continuous integration of the functional determinant yields: $\mathcal{E}_{4D}(\Lambda) = \frac{1}{2} \int_0^\Lambda \frac{2\pi^2 k_E^3 dk_E}{(2\pi)^4} \ln(k_E^2 + m^2) = \frac{1}{64\pi^2} \left[ \Lambda^4 \ln(\Lambda^2 + m^2) - \frac{1}{2}\Lambda^4 + \mathcal{O}(m^2 \Lambda^2) \right]$ The leading behavior is $\mathcal{E}_{4D}(\Lambda) \sim \Lambda^4 \ln \Lambda$, giving an effective logarithmic derivative: $\alpha_{\text{eff}}(\Lambda) = \frac{d \ln(\Lambda^4 \ln \Lambda)}{d \ln \Lambda} = 4 + \frac{1}{\ln \Lambda}$ Over our cutoff range $\Lambda \in [3.93, 15.71]$ ($\ln \Lambda \in [1.37, 2.75]$), $\alpha_{\text{eff}}$ ranges from $4.73$ to $4.36$, yielding a two-point secant slope of $\alpha_{4D} \approx 4.52\text{--}4.55$. Factoring out $\ln \Lambda$ ($\mathcal{E}_{4D} / \ln \Lambda \sim \Lambda^4$) recovers the exact quartic power law $\alpha = 4.0160 \approx 4.00$.

Benchmark Results (OpenMP Accelerated)

Sweeping lattice spacings $a \in [0.8, 0.2]$ ($N_{3D} = 64^3 = 262,144$ modes; $N_{4D} = 36^4 = 1,679,616$ modes per step) with OpenMP multithreaded reduction:

Spacing $a$	Cutoff $\Lambda = \pi/a$	$\mathcal{E}_{3D}$	$\mathcal{E}_{4D}$	Compute Time
0.8000	3.9270	3.0840	3.0981	10 ms
0.6000	5.2360	9.5182	11.8619	11 ms
0.5000	6.2832	19.5462	27.3875	11 ms
0.4000	7.8540	47.3336	75.3224	7 ms
0.3000	10.4720	148.6360	272.9260	10 ms
0.2500	12.5664	307.4222	612.0807	11 ms
0.2000	15.7080	748.9586	1632.7306	11 ms

[Scaling Analysis Exponents]:
  3D Mode Sum Scaling Exponent:            3.9620  (Theoretical: 4.0000)
  4D Raw Exponent (d ln E / d ln Λ):       4.5208  (Expected: ~ 4.55 due to leading Λ⁴ ln Λ)
  4D Factored Exponent (E / ln Λ ~ Λ^α):   4.0160  (Confirmed Quartic: 4.0000)

Both formalisms scale with leading quartic power law $\sim \Lambda^4$. Factoring out the logarithmic functional determinant enhancement ($\mathcal{E}_{4D} / \ln \Lambda \sim \Lambda^4$) yields an exact pure power law exponent of $4.0160 \in [3.85, 4.15]$ ($R^2 = 1.00000$). This demonstrates that equal-time foliation does not create zero-point divergence; the divergence is an intrinsic invariant property of continuous spacetime field degrees of freedom.

---

Section 6: Symbolic Constraint Algebra Audit

Source File: code/python/symbolic_constraint_algebra.py

Formal Algebraic Distinction

The audit symbolically checks the algebraic structure of flat-spacetime canonical field theory against General Relativity without heuristic shortcuts:

1. Flat Minkowski Spacetime (Poincaré Lie Algebra $\mathfrak{p} = \mathfrak{so}(1,3) \ltimes \mathbb{R}^4$):
   • Explicit differential vector field representations for all 10 generators: $P^\mu = \eta^{\mu\alpha} \partial_\alpha, \quad M^{\mu\nu} = x^\mu \eta^{\nu\alpha}\partial_\alpha - x^\nu \eta^{\mu\alpha}\partial_\alpha$
   • Differential operator commutators dynamically evaluated: $[M^{\mu\nu}, P^\rho] = \eta^{\nu\rho} P^\mu - \eta^{\mu\rho} P^\nu$ and $[P^\mu, P^\nu] = 0$.
   • Jacobi Identity: Explicitly evaluated and verified to vanish identically across all 120 independent generator triples: $[[M^{\mu\nu}, M^{\rho\sigma}], P^\lambda] + [[M^{\rho\sigma}, P^\lambda], M^{\mu\nu}] + [[P^\lambda, M^{\mu\nu}], M^{\rho\sigma}] \equiv 0$
   • Structure constants $f^\lambda_{\mu\nu} \in \{-1, 0, 1\}$ are strictly static scalar constants. Time translation generator $P^0 = H = \int d^3x T^{00} \neq 0$ generates physical time evolution; no constraint $H \approx 0$ exists.
2. Diffeomorphism-Invariant Gravitation (Dirac Hypersurface Deformation Algebroid):
   • Canonical phase space $(q(x), \pi(x))$ parameterized via the 1D minisuperspace / cylindrical wave reduction of General Relativity (Thiemann, Modern Canonical Quantum General Relativity, Eq. 1.2.14; DeWitt, Phys. Rev. 160:1113, 1967): $\mathcal{H}_\perp[N] = \int dx N(x) \left[ \frac{\pi(x)^2}{2\sqrt{q(x)}} + \sqrt{q(x)} \left( \frac{\partial_x q(x)}{2 q(x)} \right)^2 \right]$
   • Functional derivatives: $\delta H_\perp / \delta \pi = N \pi / \sqrt{q}$, while the metric gradient variation undergoes Palatini double integration-by-parts on the lapse: $\delta H_\perp / \delta q|_{\text{grad}} = - \frac{1}{\sqrt{q}}\partial_x^2 N$.
   • Wronskian Divergence Proof: SymPy dynamically proves the differential operator identity: $N_1 \partial_x^2 N_2 - N_2 \partial_x^2 N_1 \equiv \partial_x (N_1 \partial_x N_2 - N_2 \partial_x N_1)$
   • Integrating by parts transfers $\partial_x$ onto $\pi(x)/q(x)$, isolating the spatial diffeomorphism constraint $\mathcal{H}_x = -2 q \partial_x(\pi/q)$: $\{H_\perp(N_1), H_\perp(N_2)\} = \int dx \, [N_1(x) \partial_x N_2(x) - N_2(x) \partial_x N_1(x)] q^{-1}(x) \mathcal{H}_x(x) = H_x(K^x)$
   • The bracket depends explicitly on the phase-space inverse metric field $q^{-1}(x) = q^{xx}(x)$, certifying that it is an open Lie algebroid with dynamical structure functions rather than a Lie algebra.
   • The "problem of time" ($H_\perp \approx 0$) is an exclusive property of gauge reparameterization invariance in theories with dynamical metrics. Flat Minkowski spacetime canonical quantization operates on fixed $\eta_{\mu\nu}$ and does not inherit this constraint.

---

Section 7: Automated Verification Harness & Reproduction Guide

Complete Test Suite Execution

To run all 6 unit and integration tests:

python -m pytest papers-drafts/cases/js-riverfield/code/python/test_suite.py -v

The test harness enforces:

• Explicit existence of the C++ benchmark CSV data, raising FileNotFoundError if missing.
• Execution of plot_uv_scaling.plot_uv_scaling(csv_path) ensuring rendered PDF and PNG figures exist and are non-empty (>1 KB).
• Strict power-law regression exponent bounds $3.85 \le \text{slope}_{3D} \le 4.15$ and $3.85 \le \text{slope}_{4D} \le 4.15$ with $R^2 \ge 0.99$.
• Verification that all 120 Poincaré Jacobi triples vanish identically and ADM brackets depend dynamically on $q^{ij}(x)$.
• Full 2nd-order Klein–Gordon Cauchy convergence at $\mathcal{O}(c^{-2})$ and backward antiparticle power suppression at $\mathcal{O}(c^{-4})$.
• Smooth $C^\infty$ mollified gauge field integration with numerical line integration evaluating to machine zero ($0.0 \pm 10^{-16}$) for Riverfield's contact propagator vs. theoretical dipole radiated power $P_{\text{theory}} = \omega^2 / (12\pi c)$ ($0.2387 \pm 0.43\%$) for standard QED.
• Tachyonic Cauchy continuum stability audit verifying $\lambda_+ > 1.0$, rate matching to $\le 0.01\%$ error, and exact second-order $\mathcal{O}(\Delta t^2)$ error convergence ratio ($70 \le \text{err\_dt}_{01} / \text{err\_dt}_{001} \le 130$, measured $99.40\text{x}$).

============================= test session starts =============================
collected 6 items

test_wavepacket_asymptotic_convergence PASSED                            [ 16%]
test_gauge_propagator_radiation PASSED                                   [ 33%]
test_symbolic_poincare_algebra PASSED                                    [ 50%]
test_symbolic_dirac_hypersurface_algebra PASSED                          [ 66%]
test_uv_scaling_divergence PASSED                                        [ 83%]
test_tachyonic_cauchy_continuum PASSED                                   [100%]

============================== 6 passed in 6.42s ==============================

Full Reproduction Commands

1. Verify Lean 4 Formal Proofs:

   lake env lean papers-drafts/cases/js-riverfield/code/lean/RiverfieldRefutation.lean

2. Build and Run C++ Lattice Engine:

   cd papers-drafts/cases/js-riverfield/code/cpp
   # Via Makefile (includes -march=native -O3 -fopenmp):
   make
   # Or direct standalone compilation:
   g++ -std=c++20 -O3 -fopenmp -march=native lattice_uv_divergence.cpp -o lattice_uv_divergence.exe
   ./lattice_uv_divergence.exe
   cd ../../../..

3. Execute Python Simulation and Plotting Suite:

   python papers-drafts/cases/js-riverfield/code/python/wavepacket_convergence.py
   python papers-drafts/cases/js-riverfield/code/python/gauge_propagator_simulation.py
   python papers-drafts/cases/js-riverfield/code/python/tachyonic_cauchy_continuum.py
   python papers-drafts/cases/js-riverfield/code/python/symbolic_constraint_algebra.py
   python papers-drafts/cases/js-riverfield/code/python/plot_uv_scaling.py