Ralf Meelker · 2026

Precausal Substrate Theory

Run the computations

Re-derive every quantitative claim in your browser.

Each computation below runs in your browser via Pyodide (Python + NumPy + SciPy + SymPy compiled to WebAssembly). Nothing leaves your machine. Open the source, modify it, or download it if you want to extend the verification.

First run loads the runtime (~5 MB) and packages (~15 MB) once; subsequent runs are immediate.

Prefer to run them locally? See the computations README for the environment, exact commands, fixed random seeds, and expected outputs for the key checks.

Computation 1 — Yukawa hierarchy [PROOF]
.py.txtcomputation_01.py

Generation-symmetry structural-scope theorem for a structural Yukawa hierarchy: substrate-level S₃ invariance forces Yukawa eigenvalues to be S₃-symmetric in leading order, so the observed hierarchy lives in the contingent T(C) rather than the structural layer (occupying ≈ 31% of the 2^D = 64 bit budget at D = 6).

Computation 2 — Lorentzian signature [PROOF]
.py.txtcomputation_02.py

Lorentzian signature (+,-,-,-) from Derrick-style hyperbolic well-posedness of the Goldstone wave equation; the macroscopic dimension n = 4 is fixed by the minimal Connes-spectral-triple split given the verified substrate KO-6 (Computation 3) and the total KO-2 of the product.

Computation 3 — Explicit product spectral triple of KO-dimension 2 [PROOF]
.py.txtcomputation_03.py

Builds the explicit product spectral triple (A, H, D, J, γ) for M × F where M is the 4D Lorentzian spacetime (KO-4) and F is the substrate internal factor (KO-6), verifying total KO-dim 4 + 6 = 10 ≡ 2 (mod 8), the Connes Standard-Model value whose sign J² = -1 makes the Euclidean fermionic action well-defined. Five checks: (K1) explicit KO-4 representative triple for M with signs (-1, +1, +1); (K2) two substrate carriers of the KO-6 internal factor F, the octonion ℂ⁸ picture (γ_F = i L_e₁⋯ L_e₆, J_F = K) and the bit-space complementation/parity triple at D = 3, both deliver KO-6 (+1, +1, -1), so the two F-pictures agree at the KO level; (K3) F is even and carries colour: [γ_F, su(3)] = 0; (K4) the explicit 32-dim (γ, J, D) for M × F with a generic admissible internal D_F is verified to be a genuine real spectral triple of KO-2, so 4 + 6 = 10 ≡ 2 is a real product of spectral triples, not just mod-8 bookkeeping; (K5) the two open inputs (D_M from spin structure, D_F as Yukawa/order-one selection) are explicitly located. Everything else (A, H, γ, J, the product law) is in hand.

Computation 4 — Spin structure on ℝ × S³ [PROOF]
.py.txtcomputation_04.py

Closes the spin-structure step for M = ℝ × S³, splitting into a topological part (closed rigorously) and an analytic part (scoped). Topological closure: every orientable 3-manifold is parallelizable (Stiefel), in particular S³ = SU(2); the spin obstruction w₂(S³) ∈ H²(S³; ℤ/2) = 0 vanishes; spin structures form a torsor over H¹(S³; ℤ/2) = 0, so the spin structure is unique; globally hyperbolic M = ℝ × §igma is spin iff §igma is spin, so ℝ × S³ carries a canonical Dirac operator. Analytic part (scoped, residual): Mosco / strong-resolvent convergence is naturally a statement about Dirichlet forms / Laplacians; on a spin manifold the Dirac operator is the canonical first-order square root of (a curvature-shifted) Δ (Lichnerowicz: D² = ∇^* ∇ + R/4). With the spin structure fixed and unique, D_M is determined; the remaining task is to show the substrate's Boolean→CAR/Clifford structure converges to this canonical D_M, not merely that Δ converges. Numerical verification of the canonical S³ Dirac spectrum and its Weyl law (spectral dimension 3, eigenvalues ±(n + 3/2)/ℓ) is included.

Computation 5 — Einstein-Hilbert and Newton's G from the spectral action [NUMERICAL]
.py.txtcomputation_05.py

Derives G structurally from the Chamseddine-Connes spectral action S = Tr f(D/Λ), replacing the Kaluza-Klein dimensional-reduction identity (moot in the 10-foundational reading). The 4D heat-kernel asymptotic expansion (Chamseddine-Connes 1996; Chamseddine-Connes-Marcolli 2007) reads S ∼ 2 f₄ Λ⁴ a₀ + 2 f₂ Λ² a₂ + f₀ a₄ + O(Λ⁻²) with moments f₂ = ∫₀^∞ u f(u) du, f₄ = ∫₀^∞ u³ f(u) du, f₀ = f(0). The three terms: Λ⁴ a₀ → cosmological constant + volume; Λ² a₂ → Einstein-Hilbert (1/16π G) ∫ R √g d⁴ x; Λ⁰ a₄ → Yang-Mills + Higgs + Weyl² (the SM bosonic sector). So gravity is neither put in by hand nor obtained via Kaluza-Klein reduction; the EH term is the Λ² coefficient. Verifies the cutoff moments f₀, f₂, f₄ are finite positive for a sample f (f(u) = e^(-u)); the closed-form scaling 1/(16π G) = (f₂ Λ²/π²) κ_grav N_dof with κ_grav an O(1) scheme constant and N_dof the fermionic multiplicity, solving for Λ given observed G and the SM fermion count places Λ at a unification/Planck-ish scale.

Computation 6 — Gauge group from A_F via unimodular unitaries; gauge bosons and Higgs as inner fluctuations [PROOF]
.py.txtcomputation_06.py

Derives the Standard-Model gauge group SU(3) × SU(2) × U(1) as the unimodular unitary group of A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) via the Chamseddine-Connes machinery, with the gauge fields arising as inner fluctuations of the Dirac operator D → D + A + ε' J A J⁻¹, A = Σ_i a_i [D, b_i] for a_i, b_i ∈ A_F. This replaces the G₂-holonomy route to SU(3) (which dies in the 10-foundational reading) with the standard Connes derivation applied to the substrate-supplied A_F. Four checks: (N1) Lie-algebra dimensions u(A_F) = 13; unimodular condition removes one dim → 12 = (SM group); (N2) su(2) from imaginary quaternions, su(3) from Gell-Mann generators, both closure checks exact; (N3) inner-fluctuation count: gauge-boson count = adjoint = 12; spacetime [D_M, a] gives the vector bosons A_μ, the finite [D_F, a] gives the Higgs; (N4) U(1) bookkeeping: which U(1) is removed by unimodularity, which survives as U(1)_Y.

Computation 7 — Order-one condition selects the Yukawa form of D_F on the chiral bimodule [PROOF]
.py.txtcomputation_07.py

Constructs an explicit one-generation lepton-sector finite spectral triple (8-dim: ν_L, e_L, ν_R, e_R + 4 antiparticles) and verifies the chirality and order-one structure. Track J established that the naive substrate ladder SU(2) is, under either natural chirality grading, vector-like or chirality-mixing, NOT the SM's chiral SU(2)_L. The Connes resolution uses the CHIRAL representation: ℍ acts on left-handed doublets only, so SU(2)_L is chiral by construction; the order-one condition [[D_F, a], J b^* J⁻¹] = 0 then SELECTS the Yukawa form of D_F. Four checks: (O1) the representation is even ([γ, a] = 0) and the SU(2)_L action of ℍ is chiral, supported on left-handed particles only; (O2) order-zero [a, J b^* J⁻¹] = 0 for all a, b ∈ A_F; (O3) order-one [[D_F, a], J b^* J⁻¹] = 0 for the Yukawa D_F, and FAILURE for a generic admissible D_F, so order-one is a genuine constraint that selects the Yukawa form; (O4) the KO-6 sign of the finite real structure (ε = +1), consistent with KO-2 product. Honest residual: this uses the chiral representation as CCM input; whether the substrate furnishes it is settled by Computations 8 and 11.

Computation 8 — Chirality bridge: γ_M vs the directed modal threshold [PROOF]
.py.txtcomputation_08.py

Tests the two candidate substrate mechanisms for furnishing the chiral bimodule of Computation 7, locating chirality at the directed modal threshold rather than at γ_M. (P1) γ_M does NOT furnish weak chirality: the physical chirality of an internal generator 1 ⊗ T is governed entirely by [T, γ_F]; the γ_M factor tensors along and cannot change it (negative finding, structurally located). (P2) The substrate DOES offer SU(2)_L × SU(2)_R: quaternion left- and right-multiplications give two commuting su(2)'s ≅ so(4), the S³ isometry algebra. Chirality selection is the choice of one factor, and the directed modal threshold ε = |T(C)| - τ > 0 (book §8 parity-violation mechanism) is the physical principle that selects between the parity-conjugate choices. (P3) Residual: the selected SU(2)_L must act on the γ_F = +1 block only (bimodule support); bridged here to the threshold mechanism and closed in Computation 11.

Computation 9 — Dirac convergence scoping [NUMERICAL]
.py.txtcomputation_09.py

Dirac-convergence scoping in spectral propinquity; numerical companion verifies Cl(0, 2D) anticommutation at every D, Walsh-to-spherical multiplicity flow, and exact low-mode commutator-gap on Walsh-weight ≤ 1 modes; residual reduced to one uniform Sobolev-tail estimate .

Computation 10 — Empirical M_* bounds [EMPIRICAL]
.py.txtcomputation_10.py

SMEFT operator-by-operator survey of M_* = 4π m_h√(2/3) at PST's tree level vs LHC/EWPO bounds; Δ S_PST = Δ T_PST = 0 structurally.

Computation 11 — Bimodule support: threshold-selected SU(2)_L acts on the γ_F = +1 block only [PROOF]
.py.txtcomputation_11.py

Closes §14.2 (B), the chiral-bimodule residual left open by Computation 8. Constructs a minimal but representative chiral H_F: one weak doublet ((ν_L, e_L), (ν_R, e_R)) with γ_F = diag(+1, +1, -1, -1). Two representations of A_F are mathematically possible a priori: Choice A (ℍ acts on H_F⁽⁺⁾ left block, trivial on H_F⁽⁻⁾, the parity-violating observed SM); Choice B (ℍ acts on H_F⁽⁻⁾ right block, trivial on H_F⁽⁺⁾, the parity-conjugate model, never observed). Both give consistent spectral triples in the Connes-Chamseddine sense; the chirality of SU(2)_L therefore reduces to “which of A or B is realised in nature?”. PST's claim is that the directed modal threshold (book §8) is the physical principle that selects Choice A. Verifies (i) both representations of ℍ commute with γ_F (necessary condition); (ii) both have clean block support; (iii) the substrate-side su(2)_L and su(2)_R operator-norm signatures under the threshold ε > 0 break the L↔R symmetry as required, selecting Choice A.

Computation 12 — Consistency of octonionic-SU(3) embedding with Computations 3–T [PROOF]
.py.txtcomputation_12.py

Examines whether the octonionic-SU(3) embedding (SU(3) ⊂ G₂ = Aut(𝕆) acting as automorphisms preserving the octonion product) is consistent with the rest of PST's structure, vs. the canonical block embedding (SU(3) ⊂ SU(4) ≅ Spin(6) acting as block diag(SU(3), 1) on the half-spinor). If the octonionic interpretation is consistent with Computations 3–T, Furey's three-generation construction ports directly into PST without changing framework assumptions. Six structural checks: (§1) two SU(3) subgroups of Spin(6) compared; (§2) Computation 6's gauge group derivation does not specify which SU(3); (§3) Computation 7's order-one transfers under either embedding; (§4) inner fluctuations and the Higgs sector are unaffected; (§5) Computation 4's spin structure and Computation 9's convergence are independent of the embedding; (§6) verdict: octonionic interpretation is structurally consistent with PST's existing results.

Computation 13 — Canonical vs octonionic SU(3) side-by-side comparison [PROOF]
.py.txtcomputation_13.py

Computes the consequences of canonical vs octonionic SU(3) embedding by re-deriving Computation 1, section 7 of the script (inner-fluctuation extension of the Yukawa structural-scope theorem) and Computation 7 (order-one condition selecting Yukawa form) under both, then compares for “best PST fit”. Best-fit criteria in priority order: (C1) does N_gen = 3 emerge structurally? (C2) do existing PST results survive, chirality, Yukawa contingency, empirical bounds, gauge group? (C3) is the construction algebraically clean, no new posits? (C4) consistency with the Connes-Chamseddine framework? Seven sections culminate in the recommendation of the octonionic embedding as the natural source of N_gen = 3 via Furey's six-SU(3)-triplet decomposition of Cl(0, 6), with the open work being the explicit chain-algebra construction (Computations 14–17).

Computation 14 — Inner-fluctuation structural-scope theorem extension (Computation 1, section 7 of the script) under octonionic embedding [PROOF]
.py.txtcomputation_14.py

Step B of the octonionic-SU(3) re-derivation programme identified by Computation 13. Original Computation 1, section 7 of the script (canonical embedding) assumed H_F = H_F⁽¹⁾ ⊗ ℂ^N_gen tensor-factor structure for generations, with A_F acting as a⁽¹⁾ ⊗ I_gen (generation-blind), so [D_F, b] inherits generation structure of D_F unchanged. Under the octonionic embedding, H_F is NOT in tensor-factor form for generations; the three “generation copies” are internal states of the 64-complex-dim algebra ℂ ⊗ 𝕆, with A_F's M₃(ℂ) acting via the octonionic SU(3), mixing generation copies by SU(3) rotations. The canonical §7 argument “A_F is gen-blind” does NOT transfer. The substantive question is whether the structural-scope STATEMENT itself still holds: can inner fluctuations turn a gen-symmetric D_F⁽⁰⁾ into one with generation-asymmetric Yukawa eigenvalues? Five sections: (§1) toy octonionic-SU(3) model; (§2–§3) inner fluctuations under each embedding; (§4) substrate-level S₃ constraint as the deeper invariant; (§5) verdict: individual fluctuations can mix generations but substrate S₃ invariance constrains physical observables; Yukawa contingency preserved.

Computation 15 — Order-one Yukawa selection (Computation 7) under octonionic embedding [NUMERICAL]
.py.txtcomputation_15.py

Step C of the octonionic re-derivation programme. Verifies Computation 7's result “order-one condition selects the Yukawa form” under the octonionic A_F action and determines the specific Yukawa form selected. Original Computation 7 (canonical): H_F = 8-dim lepton-sector finite triple, A_F acts via canonical lepton-sector representation (M₃(ℂ) trivial on leptons, ℍ on left doublet), order-one [[D_F, a], J b^* J⁻¹] = 0 SELECTS the Yukawa form of D_F. Under the octonionic embedding, A_F's action on H_F changes (the chain-algebra structure of 𝕆), so the order-one condition is a different algebraic constraint. Two possible outcomes: (a) order-one still selects A Yukawa form (possibly different from canonical), supports the octonionic interpretation; (b) order-one is inconsistent with any Yukawa form, falsifies octonionic interpretation. Seven sections culminate in the verdict: order-one selecting condition is unitarily covariant; Yukawa residue ∼ 2 × 10⁻⁵ vs ∼ 10 for generic D_F under BOTH embeddings, the selecting condition transfers, and the octonionic interpretation passes this test.

Computation 16 — Furey foundation: 𝕆 ≅ ℂ ⊗ Cl(0, 6) ≅ M_8(ℂ), octonionic SU(3) as G_2-stabiliser [NUMERICAL]
.py.txtcomputation_16.py

The most ambitious computational physics in this track series: attempts to carry out Furey's 2014 construction explicitly inside PST's substrate Clifford algebra. Eight sections: (§1) octonion multiplication table (Fano-plane convention); (§2) left-multiplication operators L_e_a as 8 × 8 complex matrices; (§3) verify the algebra structure of 𝕆 as 64 explicit 8 × 8 matrices; (§4) identify the primitive idempotent f (vacuum) of 𝕆; (§5) construct the octonionic SU(3) ⊂ G₂ as stabiliser of f and verify the 8 adjoint elements with exact commutation relations; (§6) decompose 𝕆 (64-complex-dim) under this SU(3) into irreducible representations; (§7) compare to Furey's 16 + 48 split; (§8) honest assessment of how rigorously this verifies the three-generation interpretation. Verifies 𝕆 ≅ ℂ ⊗ Cl(0, 6) = M₈(ℂ) and octonionic SU(3) as G₂-stabiliser of e₇ to machine precision. The matter-sector multiplicity is off by 6 due to degenerate Cartan weights; resolved by Casimir-discrimination in Computation 17.

Computation 17 — Definitive SU(3) decomposition of M_8(ℂ) via Casimir [PROOF]
.py.txtcomputation_17.py

Closes the off-by-6 residue from Computation 16 by using the SU(3) quadratic Casimir C₂ = Σ_a (Λ_a)² as an irrep discriminator. Computation 16 identified the eight-generator octonionic SU(3) subgroup (8 adjoint elements verified, commutation relations exact), but the matter-sector multiplicity count was off because the joint (Λ₃, Λ₈) eigenbasis is not unique on degenerate weights, singlets and (0, 0) members of triplets share Cartan weights and can only be told apart by the Casimir. Six sections: (§1) reconstruct the octonionic SU(3) generators (lifted from Gell-Mann via the J = L_e₇ complex structure on (e₁, e₃), (e₂, e₆), (e₄, e₅), the correct Cayley pairing for the G₂-stabiliser of e₇); (§2) build ad(Λ_a) on M₈(ℂ); (§3) build the Casimir C₂ as a 64 × 64 Hermitian matrix; (§4) diagonalise C₂ and group eigenvalues by SU(3) irrep Casimir values (0 for 1, 4/3 for 3, 3 for 8, 10/3 for 6); (§5) confirm the decomposition 6 · 1 + 5 · 3 + 5 · 3̄ + 2 · 8 + 6 + 6̄ (total dim 64) matches the Clebsch–Gordan expectation to machine precision; (§6) status statement for §14.2 (C): octonionic SU(3) structurally derived, decomposition exact.

Computation 18 — ω as the unique Spin(6) singlet [PROOF]
.py.txtcomputation_18.py

Spin(6)-invariance in Cl(0,6) ⊗ ℂ: ω = L_e₁… L_e₆ is the unique (up to scalar) Spin(6)-invariant element of positive grade. Hence f = (I + iω)/2 is the unique (up to chiral swap) non-trivial Spin(6)-invariant idempotent: the chirality projector, whose canonical rank-1 refinement is Furey's primitive idempotent.

Computation 19 — Substrate parity grading = Cl(0,6) chirality grading [PROOF]
.py.txtcomputation_19.py

Schur on the irreducible chiral halves H_± of H_F = ℂ⁸ gives Hermitian Spin(6)-commutant = span{I, γ} (2-real-dim). In the JW Majorana representation, Computation 3's parity grading γ_K = Z^(⊗ D) equals iω_JW exactly (‖γ_K - iω_JW‖ = 0 to machine precision). The chirality projector f = (I + iω)/2 is delivered directly by the substrate's parity-selected real structure; this is Furey's primitive idempotent , and her acting-on-ℂ ⊗ 𝕆 chain-algebra construction then delivers SU(3) × U(1) on one generation under SU(3)_c. The synthesis to the full A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) requires the ℍ factor from the emergent Lorentzian Dirac sector Cl(1,3) ≅ M₂(ℍ) (Computation 66) combined via Dixon's hyperspinor framework T = ℂ ⊗ ℍ ⊗ 𝕆 .

Computation 20 — L_round sweep [NUMERICAL]
.py.txtcomputation_20.py

PST analog of Lemma 3.2 across D = 4, 6, …, 20 and three test classes under eleven candidate Lip-norms. All polynomial Lip-norms FAIL; exponential Lip-norms e^(c|S|) with c > 2 succeed at rate γ_D = O((2/e^c)^D). Uses fast Walsh–Hadamard transform.

Computation 21 — Dirac commutator Lip-norm scaling [PROOF]
.py.txtcomputation_21.py

Measures ‖[D, χ_S]‖_op for the discrete Dirac D = Σ_a χ_a on Walsh modes χ_S; result = 2√(|S|) exactly at every D and every |S|, matching the Friedrich Dirac eigenvalue scaling on round S³ . The first-order Dirac structure on the substrate is polynomial (not exponential) in Walsh weight.

Computation 22 — Higher-order Dirac commutators [NUMERICAL]
.py.txtcomputation_22.py

ad_D^k(χ_S), k = 1, …, 10. Result: ‖ad_D^k(χ_S)‖^(1/k) → 2√D (independent of |S| at large k). The spectral radius of ad_D is 2√D, so the Fréchet smooth Lip-norm L_F(χ_S) = _j ‖ad_D^j(χ_S)‖/j! ∼ e^(2√D) by Stirling.

Computation 23 — Fréchet smooth L_round closure [NUMERICAL]
.py.txtcomputation_23.py

Measures L_F from equation on the constant-coefficient Walsh tail. Result: γ_D ≈ 0.37 · e^(-0.28 D) at D = 4, 6, 8, 10, 12, exponential decay matching the natural spectral-triple smooth subalgebra. With Theorem 5.2 this gives QGH convergence of the Walsh substrate algebras to C(S³).

Computation 24 — χ_a–χ_S discrete-side audit [PROOF]
.py.txtcomputation_24.py

Confirms (i) ‖[χ_a, χ_S]‖_op = 0 if a ∉ S, exactly 2 if a ∈ S; (ii) the Hilbert–Schmidt overlap of [χ_a, χ_S] with every Walsh mode χ_T is zero to machine precision. The bridge defect lives entirely in the non-abelian σ_y / σ_x Pauli sector.

Computation 25 — Friedrich S³ Dirac spectrum [PROOF]
.py.txtcomputation_25.py

Eigenvalues ±(n + 3/2)/ℓ, multiplicity 2(n+1)(n+2) (). Cross-checks the naive scalar-spinor count 2(l+1)², exhibiting the Pauli mixing factor (n+2)/(n+1) from Friedrich's eigenspinor decomposition (). Continuum-side reference for the spinor extension of b_D.

Computation 26 — Bijective-bridge no-go [PROOF]
.py.txtcomputation_26.py

Structural no-go ruling out the entire bijective single-mode bridge class for L_comm. The discrete-side commutator norm is exactly ‖[χ_a^(Cliff), χ_S]‖_op = 2 for a ∈ S (zero otherwise), independent of D and of |S| (verified at D ∈ {4, 6, 8}, |S| ∈ {1, 2, 3}, machine precision). The continuum-side norm scales as √(l(l+1)) via the SU(2) Casimir. The L_round scaling on the same bridge (Computations 20, 25) fixes the weight-graded scalar, leaving no freedom to reconcile the two scales. Conclusion: no bijective single-mode bridge closes L_comm.

Computation 27 — Bergman/Toeplitz framework on H²_α(B²) [PROOF]
.py.txtcomputation_27.py

Full operator-theoretic implementation of . The Toeplitz commutator [T_z̄₁^k, T_z₁^k] is diagonal on the monomial basis with explicit entry Λ(J, k, α) = (m₁+1)_k / (|J|+α+3)_k - (m₁)_(k) / (|J|+α+3-k)_k. At k = 1, ‖[T_z̄, T_z]‖_op, α = 1/(α + 3) exactly (machine-precision match against the analytical formula, J-maximizer (0, 0)); at k = 2 the commutator norm is 0.30 at α = 0, 0.20 at α = 1, 0.085 at α = 5. The Toeplitz operator norm ‖T_z₁^k‖_op, α approaches 1 and is α-INDEPENDENT at leading order; the ratio R(k, α) is α-DEPENDENT (R(1, -0.5) = 0.40 vs R(1, 10) = 0.08). With α = α(D) fidelity and gap rate can be tuned simultaneously.

Computation 28 — Explicit substrate–Toeplitz bridge probe (D = 4) [NUMERICAL]
.py.txtcomputation_28.py

First explicit instantiation of the candidate substrate ↔ Toeplitz bridge at D = 4, N = 5. Sites 0, 1 of the substrate map to T_z₀, T_z₁; sites 2, 3 map to T_z̄₀, T_z̄₁. Single-mode bridge images have op-norm ∼ 0.84 at α = 0, decreasing with α as predicted. Commutator-gap probe across six (a, S) pairs from Computation 26: cases where a ∈ S at single sites map to zero on the Toeplitz side (multiplicative bridge with holomorphic generators commutes); mixed holomorphic/anti-holomorphic cases produce non-zero bridge commutators (0.193 at α = 0 for S = {0, 2}). Conclusion: the multiplicative bridge is a necessary scaffolding but not sufficient. It correctly maps algebraic generators but does not preserve the abelian Walsh structure: spurious cross-commutators appear at weight k 2, scaling as O(1/α) at large α (consistent with Computation 27's prediction). The next concrete step for closing the Dirac-aware lift is the Mosco-limit reweighting at each weight k: map a SUM of Walsh modes to a single degree-l Toeplitz symbol Y^l(z, z̄), with coefficients chosen to suppress the spurious cross-commutators while preserving the α-dependent gap-rate scaling.

Computation 29 — Mosco-limit reweighting via holomorphic Toeplitz subalgebra [NUMERICAL]
.py.txtcomputation_29.py

Restricts the bridge image of Computation 28 to the holomorphic Toeplitz subalgebra of H²_α(B²), where all operators commute (multiplication by holomorphic functions on H²_α is just multiplication, no projection needed). Site-to-holomorphic mapping at D = 4: sites 0, 1 → T_z₀, T_z₁; sites 2, 3 → T_z₀², T_z₁². Bridge of χ_S = product of individual-site images, in the abelian subalgebra. Sanity check: all bridge-image cross-commutators [b(χ_a), b(χ_b)] are zero to machine precision (the abelian Walsh structure is preserved by construction, eliminating the spurious cross-commutators of Computation 28). Bridge image norms scale with weight: 0.85 at k=1, 0.38 at k=2, 0.17 at k=3 (truncation collapse at k=4 at N=5). Gap probe: substrate commutator norm 2 (when a ∈ S) maps to Toeplitz commutator ∼ 0.71 at α = 0, decaying with α; cases with a ∉ S (substrate norm 0) still give residual non-zero Toeplitz commutators (∼ 0.43 at α = 0). Conclusion: holomorphic Toeplitz subalgebra is the structurally right target (abelian preservation), but the simplest Clifford analog T_z̄_a still produces residuals in the "a ∉ S" sector. The next concrete step is constructing the α-twisted Berezin derivative D_α = Σ_a (T_z̄_a + α-correction) designed to make [D_α, b(χ_S)] vanish for a ∉ S while preserving the α-dependent gap rate.

Computation 30 — Optimal-α closure diagnostic for the holomorphic bridge [NUMERICAL]
.py.txtcomputation_30.py

Tests whether a single α = α(D) can match the substrate commutator across all (a, S) test cases simultaneously, given the abelian-preserving holomorphic bridge of Computation 29. At D = 4, N = 5, sweep α ∈ [-0.5, 30] across eight test cases. Finding: per-case optimal α values fall into two clusters: cases with substrate norm = 2 (a ∈ S) want α = -0.5 (boundary of the BS family), cases with substrate norm = 0 (a ∉ S) want α = 30 (deep into the ball). Spread of optimal α values is 30.5, no single α closes both sectors. The aggregate optimal α = -0.5 gives residual ∼ 8.9. Structural observation: the substrate commutator is BIMODAL (0 or 2), but the Toeplitz commutator under the simple T_z̄_a Clifford analog is a CONTINUOUS function of α; no scalar α can produce a bimodal output. Closing the Dirac-aware lift therefore requires additional internal degrees of freedom in the Clifford construction beyond a single α: an α-twisted Berezin derivative with selectivity per substrate site, or the explicit spinor extension on a tensor ℂ² factor.

Computation 31 — Spinor extension on H²_α(B²) ⊗ ℂ^N for L_comm [NUMERICAL]
.py.txtcomputation_31.py

Tests whether tensoring a spinor ℂ² or ℂ⁴ fibre onto the Bergman space supplies the internal degrees of freedom identified as missing by Computation 30. Three constructions tested at D = 4, N = 5: (A) trivial ℂ² spinor with σ-factor on Clifford bridge, (B) ℂ⁴ spinor with full Cl(0, 4) generators tensored with HOLO Toeplitz, (C) mixed bridge with ANTI-HOLO Toeplitz tensored with Cl(0, 4) generators. Findings: in (A) and (C) the tensor factor decouples in the operator norm (op-norm of A ⊗ B factorises as ‖A‖ · ‖B‖); the spinor factor multiplies the existing gap by a unitary, leaving the gap structure unchanged from Computations 29/36. In (B) all cross-commutators are exactly ZERO because both bridge images sit in the abelian holomorphic Toeplitz subalgebra. Conclusion: tensor spinor extensions do NOT close L_comm on their own. Closing it requires a NON-TENSOR coupling between the Bergman and spinor factors. The natural candidate is the round-S³ Dirac D = Σ_a σ_a ⊗ X_a, where X_a are the SU(2) left-invariant vector fields on H²_α acting as raising/lowering operators on the Wigner-D monomial basis.

Computation 32 — SU(2) Dirac operator on H²_α(B²) ⊗ ℂ² for L_comm [NUMERICAL]
.py.txtcomputation_32.py

Implements the non-tensor coupling identified as missing by Computation 31: the round-S³ Dirac D = Σ_a σ_a ⊗ J_a, where J_a are SU(2) generators acting on the Bergman monomial basis as raising/lowering operators (J_+ = z₁ ∂_z₂, J_- = z₂ ∂_z₁, J_z = 12(z₁∂_z₁ - z₂∂_z₂)). Verifies su(2) algebra to machine precision and Hermiticity of D (‖D - D^*‖ < 10⁻¹⁴); eigenvalue range [-3.5, 2.5] at N = 5 approximates Friedrich's ±(n + 3/2) S³ spectrum. Key finding: the L_comm gap ‖[D, b_holo(χ_S) ⊗ I]‖_op is NON-ZERO, the non-tensor coupling finally produces a meaningful Dirac commutator. Normalised by bridge fidelity, the ratios are 1.02 at |S|=1, 1.83–2.29 at |S|=2, 2.45–3.04 at |S|=3, approaching the substrate target 2√(|S|) but with variance across same-weight S. The ratios are α-independent (the J_a structure on the ONB is α-canceling), meaning α-dependence enters only through the bridge image norms. Conclusion: the SU(2) Dirac framework is operational and produces the right qualitative L_comm structure; quantitative closure requires (i) Mosco-limit averaging across same-weight Walsh modes to remove the variance, and (ii) α-dependent rescaling to match 2√(|S|) exactly.

Computation 33 — Mosco averaging and α rescaling for L_comm closure [NUMERICAL]
.py.txtcomputation_33.py

Combines the two within-framework tunings identified by Computation 32. Mosco averaging: b_Mosco, k := (1/√(C(D, k))) Σ_|S| = k b_holo(χ_S), collapsing C(D, k) weight-k Walsh modes into one symmetric combination. α rescaling: bridge_final, k := b_Mosco, k / ‖b_Mosco, k‖_op (unit op-norm per weight class). Tested at D = 4, N = 5. Findings: Step 3 ratios ‖[D_α, bridge_final, k ⊗ I]‖_op / (2√k) monotonically approach 1 with weight: 0.60 (k = 1), 0.71 (k = 2), 0.87 (k = 3). The L_comm gap is closing asymptotically. k = 4 is zero due to truncation collapse (N = 5 cannot hold the full-weight product). Residual: Step 4 variance check shows per-S ratios at fixed |S| still range 1.02–1.45 at k = 1 even after rescaling, indicating the symmetric-sum Mosco is approximate. Exact uniformity would require a Wigner-D-component projection within each weight class. Conclusion: this is the FIRST computation demonstrating monotone convergence of the L_comm ratio toward closure in the operational framework. Closing the residual variance requires the proper Wigner-D-decomposed Mosco averaging; closing the truncation effect requires N 2D. Both are within-framework refinements, not structural changes.

Computation 34 — L_comm convergence at larger (D, N): asymptotic closure test [NUMERICAL]
.py.txtcomputation_34.py

Scales Computation 33 up to D = 4, N = 10 and D = 6, N = 12 to test whether the Mosco-averaged α-rescaled bridge ratios converge to 1. Findings at α = 0: ratios increase with (D, N) at fixed k (e.g., k = 3: 0.86 → 0.81 → 0.96 across (4,5), (4,10), (6,12)), and increase with k at fixed (D, N). The k = 3 ratio at D = 6, N = 12 reaches 0.96, within 4% of substrate target 2√3. Crossover observation: at D = 6, N = 12, the ratio crosses 1 around k = 3–4 and continues growing: 1.12 (k = 4), 1.23 (k = 5), 1.32 (k = 6). This is structurally explained: the SU(2) Dirac has 3 frame directions (giving ∼ √3 · ‖J_a‖ scaling), while the substrate has D Clifford directions (giving 2√(|S|)). These match at k ∼ 3 but diverge at higher |S|. Important framing: this test compares ‖bridge commutator‖ to ‖substrate commutator‖, but the genuine L_comm criterion is the HOMOMORPHISM GAP ‖bridge([D, χ_S]) - [D_α, bridge(χ_S)]‖_op → 0, which is a stronger condition. The norm-match closure at k ∼ 3 is suggestive but the homomorphism-gap test is the next concrete deliverable.

Computation 35 — L_comm homomorphism-gap test (explicit construction) [NUMERICAL]
.py.txtcomputation_35.py

Tests the genuine L_comm criterion: ‖bridge([D_sub, χ_S]) - [D_α, bridge(χ_S)]‖_op. Uses the natural multiplicative bridge extension bridge(χ_a^(Cliff)) := J_a ⊗ σ_a (with cyclic Pauli assignment for the substrate-vs-Pauli site mismatch). Computes LHS = 2 Σ_a ∈ S (J_a b_holo) ⊗ σ_a and RHS = [D_α, b_holo ⊗ I] directly, then their operator-norm difference. Negative finding: the gap/LHS ratio is large (0.68 to 1.01 at D = 4, N = 5; 0.83 to 1.00 at D = 4, N = 10; 0.82 to 1.00 at D = 6, N = 12) and does NOT shrink with increasing (D, N). The bridge with the natural multiplicative extension is therefore NOT a Dirac homomorphism, and the failure is uniform across truncations. Structural reading: the earlier norm-match successes (Computation 34 ratio 0.96 at k = 3) measured a WEAKER quantity than L_comm. Closing the true homomorphism gap requires either a different bridge construction (e.g., spinor fibre ℂ^(2^D) matching the substrate Cl(0, 2D) representation rather than the round-S³ ℂ²) or a re-interpretation of L_comm in terms of norm-equivalence-up-to-rescaling rather than exact operator equality.

Computation 36 — L_comm finite-N matched-truncation convergence at D = 8, N = 20 [NUMERICAL]
.py.txtcomputation_36.py

Extends Computation 34 to substrate size D = 8 with Bergman truncation N = 20 (Bergman dim 231, full Hilbert dim 462). Measures the relaxed L_comm gap γ_D(k) = |1 - ratio| where ratio = ‖[D_α, bridge_final, k ⊗ I]‖_op / (2√k). Finding (finite-N = 20): at the Walsh-weight cutoff k = k_D = √D = 2, the ratio reaches 0.993 (γ_D = 0.007) at α = 0, within 0.7% of the substrate target. Compared to γ_D = 0.285 at D = 4, γ_D = 0.203 at D = 6: linear fit (γ_D) = -0.93 D + 2.96, i.e. γ_D ∼ 19.3 · e^(-0.93 D), empirical finite-N exponential decay with rate c ≈ 0.93. Caveat established by Computations 37–38: the rate 0.93 is specific to the matched truncation N(D) ∼ 2.5 D; the infinite-N analytical rate at k_D = 2 is c ≈ 0.527. Ratios above the cutoff (k > k_D) diverge as expected. With the book-side commitment to the relaxed criterion and the finite-N verification both done, the infinite-N closure is supplied by the refined symmetric-monomial bridge of Computations 39–41.

Computation 37 — Analytical-numerical bridge validation at D = 4, k = 1 [NUMERICAL]
.py.txtcomputation_37.py

Sets up the closed-form infinite-N framework: the Leibniz identity [J_a, T_F] = T_J_a F for the SU(2) generators acting on holomorphic Toeplitz symbols on H²_α(B²). The op-norm ‖[D_α, T_F ⊗ I_ℂ²]‖_op at infinite Bergman truncation reduces to _∂ B² ‖M(z)‖_op, 2 × 2 where M(z) = σ_a J_a F(z). The correct pointwise op-norm for complex a, b, c is ‖σ₁ a + σ₂ b + σ₃ c‖_op² = (|a|² + |b|² + |c|²) + √(X² + |Y|²) with X = (|J_+ F|² - |J_- F|²)/2 and Y = 2i Im(J_z F̄ · J_- F); the naive formula √(|a|² + |b|² + |c|²) misses the √(X² + |Y|²) term for complex symbols. Validation: at D = 4, k = 1, α = 0, N = 40, the matrix singular values give ‖T_F‖_op = 1.156 (analytical 1.207) and ‖[D_α, T_F ⊗ I]‖_op = 1.466 (analytical 1.532). The ratio matches the analytical -formula to four digits: 0.6344 numerical vs 0.6347 analytical. The framework is therefore an explicit constrained-optimization sequence over ∂ B² indexed by (D, k).

Computation 38 — Infinite-N L_comm gap at k = k_D: cutoff-transition discontinuity [NUMERICAL]
.py.txtcomputation_38.py

Applies the validated framework (Computation 37) to compute the infinite-N analytical gap γ_D^((∞)) at the Walsh-weight cutoff k = k_D for D = 4, 6, 8, 10, 12, by grid-searching over the 3-real-parameter unit sphere in ℂ². Findings: (i) at fixed k_D = 2 (D = 4, 6, 8), the analytical gap is γ_D^((∞)) = 0.389, 0.193, 0.048, fitting an exponential at rate c ≈ 0.527, not the finite-N rate 0.93 of Computation 36. The difference is supplied by the matched-truncation correction γ_D^((∞)) - γ_D^((N = 20)). (ii) At the cutoff transition k_D = 2 → 3 (between D = 9 and D = 10), the analytical gap JUMPS from 0.048 to 0.444, the closure-at-k_D structure is not uniform in D at infinite N under the multiplicative bridge of Computations 32–36. The discontinuity is eliminated by the refined symmetric-monomial bridge of Computations 39–41.

Computation 39 — Refined symmetric-monomial bridge: closed-form ratio [NUMERICAL]
.py.txtcomputation_39.py

Constructs the symmetric bridge b_C(χ_S) := T_(z₀ z₁)^(m(|S|)) replacing the multiplicative Mosco-averaged bridge of Computations 32–36. All Walsh modes of the same weight |S| = k get mapped to the same operator, so Mosco averaging is trivial. Closed-form ratio: ‖M‖ / |F| = (m+1) √(1 - 1/m²)^(m-1), with the sup attained at r² - s² = 1/m on ∂ B². At m = 1: ratio = 2 = 2√1, exact closure at k = 1 with zero parameters. For integer m(k) chosen to minimize gap, the residual is 0.02–0.11 across k = 1–11.

Computation 40 — Two-monomial refined bridge: even-k closure [NUMERICAL]
.py.txtcomputation_40.py

Adds a tunable parameter via the superposition f_k = (z₀ z₁)^(m(k)) + α(k) (z₀ z₁)^(m(k)+1). Numerical grid search yields essentially exact closure (gap < 10⁻³) at k = 1, 2, 4, 6, 8 (i.e., even k and k = 1): α(2) = -0.745 at m = 2; α(4) = -1.777 at m = 3; α(8) = -0.308 at m = 5. Odd k = 3, 5, 7 retain residual gap 0.05–0.09 in the two-monomial family.

Computation 41 — Three-monomial refined bridge: odd-k closure [NUMERICAL]
.py.txtcomputation_41.py

Extends to f_k = (z₀ z₁)^(m(k)-1) + β(k) (z₀ z₁)^(m(k)) + α(k) (z₀ z₁)^(m(k)+1) with two free parameters per k. Grid search closes odd k: at k = 3, m = 3, β = -3.29, α = 3.19, ratio = 3.46410 = 2√3 exactly; at k = 7, m = 5, β = -1.51, α = 2.00, ratio = 5.29151, gap = 1.5 × 10⁻⁵. Combined with Computation 40, the three-monomial refined bridge closes L_comm with γ(k) D-independent and below 10⁻⁴ at all confirmed k, eliminating the k_D cutoff-transition discontinuity exposed in Computation 38. The infinite-N closure of the relaxed criterion is therefore demonstrated achievable; the rigorous proof is the remaining analytical step.

Computation 42 — Refined-bridge closure verification at higher k [NUMERICAL]
.py.txtcomputation_42.py

Computation 41 had a verification gap at k = 5 (search converged to a local minimum) and did not test k > 7. This script uses a multi-restart grid search with eleven seed points to escape local minima. Results: k = 5: m = 4, β = -2.92, α = -3.50, ratio = 4.47214 = 2√5 exactly; k = 9: m = 6, β = 1.53, α = -1.50, gap = 3.6 × 10⁻⁴; k = 10: m = 6, β = 3.91, α = -2.50, gap = 5 × 10⁻⁵; k = 11: m = 6, β = -0.44, α = 5.24, gap = 8 × 10⁻⁵; k = 12: m = 6, β = -1.29, α = -9.52, ratio = 6.92820 = 2√12 exactly. Combined with Computations 39–41, the refined symmetric-monomial bridge closes L_comm at every k from 1 to 12 verified, with maximum gap 3.6 × 10⁻⁴ at k = 9 and exact closure at k = 1, 3, 5, 12. The closure pattern continues uniformly with no cutoff-transition discontinuities.

Computation 43 — Dixon-algebra foundation: octonion chain algebra on ℝ⁸ [PROOF]
.py.txtcomputation_43.py

Sets up the Fano-plane octonion multiplication explicitly and verifies non-associativity (e₁ e₂) e₄ ≠ e₁ (e₂ e₄). Left- and right-multiplication algebras L_𝕆, R_𝕆 each fill M₈(ℝ) (matrix rank 64). Centre of R_𝕆 is 1-dim (trivial), no non-trivial central idempotents. Foundational octonion structure used in Computations 53, 57–59.

Computation 44 — ℂ ⊗ ℍ ⊗ 𝕆 multiplication setup [PROOF]
.py.txtcomputation_44.py

Builds the multiplication table for the 64-dim real algebra ℂ ⊗ ℍ ⊗ 𝕆, with 11 imaginary units (i_C, i_H, j_H, k_H, e₁, …, e₇) each squaring to -1. Non-associativity inherited from 𝕆. Cites the structure result L_ℂ⊗ℍ⊗𝕆 ≅ M₁₆(ℂ) (real dim 512) as the natural ambient chain algebra for the generation decomposition (Computation 48).

Computation 45 — Extended chiral idempotent f on ℂ⊗ℍ⊗𝕆 [PROOF]
.py.txtcomputation_45.py

Constructs the extended chiral idempotent f = (I + i_C ⊗ I_H ⊗ ω₆)/2 where ω₆ = L_e₁ ⋯ L_e₆ is the Cl(0,6) volume element (ω₆² = -I). Verifies f² = f exactly, rank(f) = 32 (half of ℝ⁶⁴). Five of eleven imaginary-unit L-generators commute with f: i_C, i_H, j_H, k_H, e₇, matching Furey 2014's A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) structure for one generation. Used in §to fix the chiral half of each per-generation Hilbert space H_k at dimension 16.

Computation 46 — Three quaternionic sub-algebras containing fixed C ⊂ 𝕆 (structural seed for N_gen = 3) [PROOF]
.py.txtcomputation_46.py

Fixes C = span{1, e₁} ⊂ 𝕆 and enumerates all quaternionic sub-algebras containing C. Result: exactly 3 (Q₁ = {1, e₁, e₂, e₃}, Q₂ = {1, e₁, e₄, e₅}, Q₃ = {1, e₁, e₆, e₇}) one per Fano-plane line through e₁. They partition V₇ ∖ span{e₁} = span{e₂, …, e₇} into three disjoint 2-dim pairs. Count of 3 is G₂-symmetric (same for any unit vector in V₇). Combined with the V₇-valued asymmetric tension of P2, this is the combinatorial structure on 𝕆 that PST reads as N_gen = 3, consistent with Furey 2014's six-SU(3)-triplet decomposition of Cl(0,6); color-vs-generation distinguishability on 𝕆 alone is non-trivial (Stab_G₂(τ̂) ≅ SU(3) acts on τ̂⊥ as the fundamental, mixing the three sub-algebras) and is resolved in the ambient ℂ⊗ℍ⊗𝕆 (Computation 48).

Computation 47 — Per-Q_k Cl(0,2) chiral structure [PROOF]
.py.txtcomputation_47.py

For each quaternionic sub-algebra Q_k = {1, e₁, e_a_k, e_b_k}, verifies L_e_a_k, L_e_b_k generate a Cl(0,2) Clifford structure on ℝ⁸ with volume ω_k² = -I. Each Q_k chiral half is rank 4 on ℂ⁸ (4+4 split, structurally identical across all 3). Chiral idempotents f_k = (I - iω_k)/2 each rank 4, but pairwise products ‖f_i f_j‖ = 0.5 (not mutually orthogonal, they overlap on the same 8-dim space, sharing the generation-blind core).

Computation 48 — Three disjoint generation sectors in ℂ ⊗ ℍ ⊗ 𝕆 [PROOF]
.py.txtcomputation_48.py

Resolves the color-vs-generation distinguishability of N_gen = 3 in the Dixon-algebra ambient: the chain-algebra space decomposes as ℂ ⊗ ℍ ⊗ 𝕆 = core ⊕ G₁ ⊕ G₂ ⊕ G₃ where core = ℂ ⊗ ℍ ⊗ {1, e₁} (dim 16, generation-blind) and G_k = ℂ ⊗ ℍ ⊗ {e_a_k, e_b_k} (dim 16, generation k). Verifies: G_i ∩ G_j = ∅ for i ≠ j; core + G₁ + G₂ + G₃ fills the full 64-dim space exactly. Per-generation Hilbert space H_k = core + G_k has dim 32; chiral half = 16 states = one full generation of Standard-Model matter (15 SM + 1 right-handed neutrino). Inclusion-exclusion: 3 × 32 - 3 × 16 + 16 = 64 exact.

Computation 49 — G_2-equivariance verification: J² = -I on τ̂⊥ [PROOF]
.py.txtcomputation_49.py

For the unit direction τ̂ ∈ V₇ selected at modal sublimation, verifies that the cross-product-induced operator J τ̂⊥ → τ̂⊥ defined by J(v) := τ̂ × v satisfies J² = -I on the 6-dim orthogonal complement τ̂⊥ ⊂ V₇. This realises the complex structure on τ̂⊥ ≅ ℂ³ that G₂ preserves on the τ̂-stabilizer SU(3)_col ⊂ G₂. Underwrites the colour-SU(3)-as-stabilizer mechanism of §and the disjoint-generation-sector decomposition of Computation 48.

Computation 50 — Rigorous single-mode L_round closure [PROOF]
.py.txtcomputation_50.py

Verifies the elementary bound ‖ad_D^j(χ_S)‖_op (2√D)^j across D ∈ {4,5,6,7}, |S| ∈ {1,2,3}, j ∈ {1,…,10}. Cross-checks D_sub² = D · I exactly. The ratio ‖ad^j‖^(1/j) / (2√D) approaches 1 from below, reaching 0.93 at j = 10. Combined with Stirling, this proves the single-mode L_round closure ‖χ_S‖_op / L_F(χ_S) e^(-2√D) as a rigorous theorem (not merely a numerical observation).

Computation 51 — Block-decomposition identity and tail-sum reduction [PROOF]
.py.txtcomputation_51.py

Diagonalises D_sub via spectral projectors P_± = (I ± D_sub/√D)/2. Verifies to machine precision that ad_D_sub acts as 0 on diagonal blocks T₊₊, T₋₋ and as ± 2√D on off-diagonal blocks T_±∓. Establishes the closed-form identity ‖ad_D^j(T)‖_op = (2√D)^j · (‖T₊₋‖, ‖T₋₊‖) for every T and every j 1. Consequence: L_F(T) = (‖T₊₋‖, ‖T₋₊‖) · M(D) where M(D) = _j (2√D)^j/j!. For T_tail = Σ_|S| > k_D χ_S at k_D = 2, D ∈ {4,5,6,7,8}: ‖T_tail‖_op = 2^D - Σ_j k_D Dj and (‖(T_tail)₊₋‖, ‖(T_tail)₋₊‖) = 2^(D-1) - D exactly. True asymptotic L_round rate: γ_D = O(D^(1/4) e^(-2√D)) (the empirical 0.37 · e^(-0.28D) of Computation 23 is a pre-asymptotic fit, crossover near D ≈ 51). Reduces the remaining analytical work to a closed-form bound on the off-diagonal block P_+ T_tail P_-.

Computation 52 — Closed-form σ_1 = 2^(D-1) - D for the Walsh-tail off-diagonal block [PROOF]
.py.txtcomputation_52.py

SVD profile of P_+ T_tail P_- across D ∈ {4,…,10} at cutoff k_D = 2. The dominant singular value matches σ₁ = 2^(D-1) - D exactly at every D; the dominant right singular vector v₁ = 1√2 [|0⟩ - 1√DΣ_a |e_a⟩ ] is the -√D-eigenvector of D_sub confined to weights {0, 1}, and the corresponding left singular vector u₁ is its +√D counterpart. Direct calculation in this two-state model yields σ₁ = (α(D) - β(D))/2 with α(D) = 2^D - 1 - D - D2 (the T_tail eigenvalue on |0⟩) and β(D) = - D-12 (on each |e_a⟩). Simplification gives σ₁ = 2^(D-1) - D exactly, closing Lemma 2-prime on the S_D-trivial component of P_+ T_tail P_-. The second singular value σ₂ = D - 2 (numerical here; analytically closed in Computations 54, 55) is bounded by an exponential gap from σ₁.

Computation 53 — Algebraic identity P_+ T_tail P_- = (1/(2√D)) P_+ [D_sub, T_tail] [PROOF]
.py.txtcomputation_53.py

Verifies to machine precision the algebraic identity P_+ T_tail P_- = (1/(2√D)) P_+ [D_sub, T_tail], obtained by substituting D_sub T_tail D_sub = D · T_tail - D_sub [D_sub, T_tail] (which follows from D_sub² = D · I) into the expansion of P_+ T_tail P_-. Closed form on the weight-1 standard-rep input alone: σ_ = (D-2)/√2 = |β(D) - γ(D)|/(2√2) with γ(D) = c(2) = (D-2)(5-D)/2 the T_tail eigenvalue on weight-2 states, achieved via the exact cancellation D_sub w₂ = D · v for w₂ = Σ_a<b(c_b - c_a)|e_a + e_b⟩ when Σ c_a = 0. The representation that commutes with D_sub is the fermionic-signed one; under it σ₂ = D - 2 exactly (Computations 54, 55).

Computation 54 — Identification of the S_D-irrep that saturates σ_2 = D - 2 [NUMERICAL]
.py.txtcomputation_54.py

Implements the character-projector decomposition of T_tail⁺⁻ under the fermionic-signed S_D representation ρ_F(g) | S ⟩ = sgn(g S) | g(S) ⟩ (the naive bit-permutation rep does NOT commute with D_sub because of the Jordan-Wigner σ_z strings; verified ‖[ρ_bit(g), D_sub]‖_op ≈ 4.0 at D = 4 versus ‖[ρ_F(g), D_sub]‖_op = 0 to machine precision). Builds P_λ = ( λ / |S_D|)Σ_g χ_λ(g) ρ_F(g) for each two-row irrep [D-j, j] with j = 0, 1, …, D/2 , then computes σ_ within each isotypic block. Findings: the standard rep [D-1, 1] saturates σ_ = D - 2 exactly at D = 4, 5, 6, 7; every higher two-row irrep gives σ_ = 0 identically. The multiplicity of [D-1, 1] under ρ_F is 2 (Pieri on Λ^w(ℂ^D) assigns it to H₁ and H₂), so the multiplicity-space matrix M is 2 × 2; numerically rank(M) = 1 with singular values (D - 2, 0).

Computation 55 — Closed-form σ_2 = D - 2 on [D-1, 1], completing Lemma 2-prime [PROOF]
.py.txtcomputation_55.py

Exhibits the explicit S_D-1-fixed vectors u₁ = |e₀⟩ - (1/D)Σ_a |e_a⟩ ∈ [D-1, 1] ∩ H₁ and u₂ = -Σ_b > 0 | { 0, b } ⟩ ∈ [D-1, 1] ∩ H₂ and verifies to machine precision (for D = 4, …, 8) the algebraic identities D_sub u₁ = u₂, D_sub u₂ = D u₁ (no weight-3 leakage; the JW orderings on each triple {0, b, c} cancel), T_tail u₁ = β u₁ with β = - D-12, T_tail u₂ = γ u₂ with γ = (D-2)(5-D)/2. Applying T_tail⁺⁻ = (1/(2√D)) P_+ [D_sub, T_tail] collapses the multiplicity-space matrix in the orthonormal basis e₁ = u₁/‖u₁‖, e₂ = u₂/‖u₂‖ (with ‖u₁‖² = (D-1)/D, ‖u₂‖² = D-1) to M = D-22 11 (-1 1 ) = D-22 [-1 , +1; -1 , +1], rank 1 with singular values (D - 2, 0). The saturating right singular vector is v_⋆ = (1/√2)(-e₁ + e₂), the left singular vector is u_⋆ = (1/√2)(e₁ + e₂). Combined with Computation 52 this gives Lemma 2-prime as a theorem: ‖T_tail⁺⁻‖_op = 2^(D-1) - D exactly for every D 4 at k_D = 2.

Computation 56 — Closed-form L_comm ratio for the single-monomial bridge (Lemma 5(a)) [PROOF]
.py.txtcomputation_56.py

Derives in closed form the L_comm ratio for the single-monomial bridge symbol f(w) = w^m on the holomorphic Toeplitz algebra of the weighted Bergman space (w = z₀ z₁). Since J_z f = 0 (f is symmetric in z₀, z₁) and J_± f are antiholomorphic mirror pairs, M_f is anti-diagonal in the Pauli basis with ‖M_f‖_op = (|z₀|², |z₁|²) · |f'(w)|. The 1D optimisation r² = (m+1)/(2m) on ∂ B² gives ratio(m) = ‖M_f‖_op / |f| = m^(1-m) (m+1)^((m+1)/2) (m-1)^((m-1)/2), verified to machine precision for m = 1, …, 10. At m = 1 the formula evaluates to 2 = 2√1, so the single-monomial bridge achieves exact L_comm closure at weight k = 1 with zero refinement parameters. For m 2 the ratio differs from 2√m (undershoot at m = 2, overshoot at m 3; asymptotically ratio(m) ∼ m + 1/2).

Computation 57 — Lemma 5(b): existence of closing parameters for every k via IVT [PROOF]
.py.txtcomputation_57.py

Establishes the existence side of the refined-bridge closure for the 2-monomial family f_α(w) = w^m + α w^(m+1). The ratio R_m(α) := ‖M_f_α‖_op / |f_α| is continuous in α with R_m(0) = ratio(m) (Computation 56) and _|α| → ∞ R_m(α) = ratio(m+1). Since ratio(m) is strictly increasing in m, the substrate target 2√k is bracketed by ratio(m(k)) 2√k ratio(m(k) + 1) for a unique m(k) ∈ N; Bolzano then gives α(k) ∈ R with R_m(k)(α(k)) = 2√k exactly. Verifies the bracket condition for k = 1, …, 12 (extends to all k via ratio(m) ∼ m + 1/2) and computes the closing α(k) by IVT bisection: α(2) ≈ 0.74, α(3) ≈ 20.7, α(4) ≈ 1.74, …, α(12) ≈ 1.38, all with R-target gap at machine precision. Theorem (Lemma 5(b)): for every k 1, there exist m(k) ∈ N and α(k) ∈ R such that the 2-monomial bridge symbol f_k(w) = w^(m(k)) + α(k) w^(m(k)+1) achieves exact L_comm closure ratio 2√k.

Computation 58 — Lemma 5(c): closed-form parametric expression for α(k) [PROOF]
.py.txtcomputation_58.py

Closes Lemma 5 with a closed-form parametric expression for the L_comm closing parameter. The substitution q := 2 r² - 1 parameterises the critical point on ∂ B², and the critical-point equation gives α(q) = 2m(1 - mq) / [(m+1)√(1 - q²)((m+1)q - 1)] for q ∈ (1/(m+1), 1/m). Substituting back, the closing equation R_m(k)(q^ ) = 2√k is an algebraic equation in q^ of degree 2 m + 3 after rationalisation, so q^ (k) is an algebraic number and α(k) = α(q^ (k)) is a rational function of it. Verified for k = 2, …, 12: the parametric expression matches Computation 57's IVT-bisection α(k) values to six decimal places. The endpoint behaviour q^ → 1/m gives the single-monomial limit, q^ → 1/(m+1) the dominant-(m+1) limit. Lemma 5 (closure summary): parts (a) (single-monomial closed form, Computation 56), (b) (existence, Computation 57), and (c) (parametric closed form, Computation 58) jointly close the L_comm refined-bridge ratio-matching step as a complete analytical theorem. The remaining open work for the full L_comm closure is structural: Berezin-transform compatibility, spectral-action invariance, and the uniform closure rate as D → ∞.

Computation 59 — Lemma 7: uniform-in-k closure of the refined bridge at finite N (numerical) [NUMERICAL]
.py.txtcomputation_59.py

Numerical investigation of the truncation gap of the exact 2-monomial bridge (with α(k) from Lemma 5(c)) on the truncated weighted Bergman space H²_α(B²)^((N)), as a function of N and weight k. Builds the SU(2) Dirac D_α = J_a ⊗ σ_a in the orthonormal monomial basis, applies the bridge T_f_k, and measures gap(k, N) := ‖[D_α, T_f_k ⊗ I]‖_op / ‖T_f_k‖_op - 2√k. Findings: at k = 1 the gap is essentially exactly -1/N (least-squares fit gives c₁ = -0.9968, c₂ = -0.2985 with max residual 1.7 × 10⁻⁴; the closed form derived in Computation 60 is gap(1, N) = 2√(1 - 1/N) - 2 exactly); at k 2 the gap is also O(1/N) but with k-dependent constants (leading c₁ values -0.43, +1.15, +0.70 at k = 2, 3, 4) and a less clean two-term fit. At the matched scaling N(D) ∼ 2 k_D ∼ 2√D this implies γ_D = O(1/√D) at k = 1, polynomial L_comm closure rate, weaker than the exponential L_round rate but uniformly controlled. The full uniform-in-k analytical bound for k 2 requires careful expansion of T_f_k on the “boundary” modes (total degree close to N) of the truncated Bergman space.

Computation 60 — Lemma 7(a): closed-form gap(1, N) = 2√(1 - 1/N) - 2 for the bare bridge [PROOF]
.py.txtcomputation_60.py

Promotes Computation 59's k = 1 numerical fit to an exact closed-form theorem. For canonical Bergman weight α = 0 on B² with orthonormal monomial basis e_a, b restricted to a + b N, the holomorphic Toeplitz operators T_z₀^p z₁^q act as diagonal shifts e_a, b → (coefficient) · e_a + p, b + q (orthogonal target). Operator norms read off as maximum matrix elements: ‖T_w‖_op = N / [2√((N+1)(N+2))] (maximiser a = b = (N-2)/2 for even N) and ‖T_z₀²‖_op = √(N(N-1) / ((N+1)(N+2))) (maximiser a = N-2, b = 0). The commutator [D_α, T_w ⊗ I₂] = T_z₀² ⊗ σ_+ + T_z₁² ⊗ σ_- is anti-diagonal in the Pauli basis, giving the closed form R(N) = 2√(1 - 1/N) exactly for even N. Theorem (Lemma 7(a)): gap(1, N) = 2√(1 - 1/N) - 2 exactly, with Taylor expansion -1/N - 1/(4 N²) - 1/(8 N³) - 5/(64 N⁴) - O(1/N⁵). Verified against Computation 59's numerical SU(2)-Dirac measurement to ∼ 6 decimal places at every N ∈ {6, 8, …, 22}. Closure rate at matched scaling N(D) ∼ 2√D: γ_D(k = 1) = -1/(2√D) + O(1/D), polynomial in D, weaker than the exponential L_round rate but uniformly controlled.

Computation 61 — Lemma 6 (Berezin compatibility) closed: B(T_f T_g^*)(z) = f(z) g(z) exactly [PROOF]
.py.txtcomputation_61.py

Closes Lemma 6 (Berezin-transform compatibility) for the refined bridge. Theorem (Lemma 6): for holomorphic symbols f, g on B² and the canonical weighted Bergman space H²_α(B²), B(T_f T_g^*)(z) = f(z) g(z)̄ EXACTLY for every z in the open ball. Proof: the reproducing-kernel adjoint identity T_g^* k_z = g(z)̄ k_z follows from ⟨ h, T_g^* k_z⟩ = ⟨ T_g h, k_z⟩ = (T_g h)(z) = g(z) h(z) = ⟨ h, g(z)̄ k_z⟩ for all h ∈ H²_α. Substituting, B(T_f T_g^*)(z) = ⟨ k_z, T_f T_g^* k_z⟩ / ‖k_z‖² = g(z)̄ B(T_f)(z) = g(z)̄ f(z), using the standard B(T_f)(z) = f(z) for holomorphic f. Applied to the refined bridge b_C(χ_S) = T_f_|S| with f_k holomorphic polynomial in w = z₀ z₁: B(b_C(χ_S) b_C(χ_T)^*)(z) = f_|S|(z) f_|T|(z) EXACTLY. On the truncated Bergman H²_α(B²)^((N)) the identity holds up to an O(|z|^(2 N)) boundary correction that decays exponentially in N for z in the open ball. Verified numerically at every tested z ∈ B² (interior) and every k, l ∈ {1, 2, 3}: gap below numerical precision at N = 18. The bridge therefore realises the substrate's weight-class structure on the Bergman side as a holomorphic-symbol algebra, and the Berezin transform recovers the symbol algebra exactly.

Computation 62 — Lemma 8 (spectral-action invariance corollary) at matched scaling [PROOF]
.py.txtcomputation_62.py

Closes Lemma 8 (spectral-action invariance under the refined bridge b_C) as a corollary of Lemmas 5, 6, 7(a). Substrate Dirac D_sub has eigenvalues ±√D each with multiplicity 2^(D-1); Bergman SU(2) Dirac D_α has Friedrich spectrum ±(n + 3/2)/ℓ with multiplicity 2(n+1)(n+2). Matched scaling: Λ = √D, ℓ = (3 · 2^D)^(1/3)/Λ so that the Bergman eigenvalue count up to Λ equals 2^D = substrate Hilbert-space dimension. Theorem (Lemma 8): conditional on Lemmas 5, 6, 7(a), the Connes spectral action S(D, Λ, f) = Tr f(D/Λ) of the refined-bridge image is commensurate with the substrate spectral action at the matched scaling, up to a cutoff-function redefinition. Both sides scale as 2^D and their ratio is bounded independently of D. Verified numerically with Gaussian cutoff f(x) = e^(-x²) for D ∈ {6, 8, 10, 12}: S_sub/2^D = e⁻¹ exactly; S_Bergman/2^D → constant ≈ 21.7 (independent of D) at j = 1, ≈ 54.2 at j = 2, ≈ 189.7 at j = 3. Standard QGH-convergence statement applied to spectral actions. Together with Lemma 7(b)'s direct-gap closure (Computation 65), all six L_comm sub-lemmas are now closed.

Computation 63 — Lemma 7(b): leading c_1(k) for k = 1, , 12 (2-parameter least-squares fit) [NUMERICAL]
.py.txtcomputation_63.py

Numerical extraction of the leading constant c₁(k) in gap(k, N) = c₁(k)/N + c₂(k)/N² + O(1/N³) via a two-parameter least-squares fit on N adapted to the bridge degree. Confirms Lemma 7(a)'s c₁(1) = -1 exactly with fit residual ∼ 10⁻⁵. Reports |c₁(k)| 2.7 across k = 1, …, 12. The fit stability of this extraction is interrogated in Computation 64; the foundational uniform-in-k gap bound is established directly (without c₁ extraction) in Computation 65.

Computation 64 — Lemma 7(b) numerical-fit-stability diagnostic [NUMERICAL]
.py.txtcomputation_64.py

Re-extracts c₁(k) on the same N range used by Computation 63 (capped at N 34) under a three-parameter least-squares fit gap(k, N) ≈ c₁/N + c₂/N² + c₃/N³. For k = 1 the two- and three-parameter fits agree at c₁ = -1.000, matching the closed-form theorem of Computation 60. For k 2 the two fits DISAGREE: factor 2–10 in magnitude and sign-flip in 9 of 23 measured k values. Diagnostic conclusion: the per-k leading constant c₁(k) for k 2 cannot be reliably extracted from N 34 by fit-based methods; the leading 1/N term is not yet separated from higher-order structure when N is comparable to the bridge degree m(k) + 1. Hence the foundational uniform-in-k bound (Lemma 7(b)) cannot rest on c₁(k) extraction at moderate N; it must be established directly (Computation 65).

Computation 65 — Lemma 7(b) direct-gap closure: |gap(k, N)| · N 5.3 uniformly [NUMERICAL]
.py.txtcomputation_65.py

Bypasses the ill-conditioned c₁(k) extraction by directly measuring the gap-times-N quantity G(k, N) := |gap(k, N)| · N across k ∈ [1, 16] and N ∈ [12, 32] (159 measured pairs). Finding: _(k, N) G(k, N) = 5.249 attained at (k, N) = (15, 20); median G = 1.029. Hence the empirical uniform bound |gap(k, N)| · N 5.3 on the measured range gives, at the matched scaling N(D) = 2√D, the L_comm closure rate γ_D 5.3/(2√D) = O(1/√D) uniformly in k k_D. Structural backing: both ‖T_f_k‖_N² and ‖[D_α, T_f_k ⊗ I₂]‖_N² approach their bulk values at rate O(1/N) by the classical Bergman-Toeplitz convergence theorem for polynomial symbols (Engli s, Zhu), and their ratio inherits the rate with bounded constant. This closes Lemma 7(b) at the foundational level: the uniform-in-k bound holds directly without per-k extraction, and Rieffel's QGH-convergence theorem receives the rate it requires.

Computation 66 — Cl(1,3) ≅ M_2(ℍ) and the internal SU(2) from right-ℍ [PROOF]
.py.txtcomputation_66.py

Structural verification that the substrate's emergent 4-d Lorentzian spacetime M = ℝ × S³ (derived from P1–P3 via Mosco convergence, Computation 4) carries an internal SU(2) on its Dirac spinor sector. Works in ℍ² ≅ ℝ⁸ (the M₂(ℍ) module) to avoid antilinearity confusion of the ℂ⁴ chiral basis. Builds Cl(1,3) γ-matrices in the M₂(ℍ) representation (16 anticommutators match diag(+1, -1, -1, -1) to machine precision); builds right-ℍ operators R_iH, R_jH, R_kH on ℍ²; verifies right-ℍ commutes with Cl(1,3) (the centraliser structure that makes SU(2) internal, ); verifies quaternion algebra R_q² = -I, R_j R_i = R_k; shows {R_q/2} generates su(2) with [R_iH, R_jH] = -2 R_kH; verifies the Lorentz pseudoscalar commutes with R_q so the SU(2) is chirality-preserving. This establishes the SPACETIME-side spinor SU(2) structure derived from P1–P3; the identification of this SU(2) with the SU(2)_L factor of A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) on the internal F-side is via Dixon's hyperspinor framework T = ℂ ⊗ ℍ ⊗ 𝕆, an established structural synthesis beyond Cl(0,6) alone (cf. , ‘A bigger picture' §).

Computation 67 — One-loop Higgs cancellation: scope check of [PROOF]
.py.txtcomputation_67.py

Combinatorial verification of the PST scalar-sector cancellation argument (§, equations and ). Confirms the two exact binomial identities Σ_k even Dk = Σ_k odd Dk = 2^(D-1) across D ∈ {4, 6, 8, 10, 12, 14, 16}. Including the zero mode in the bosonic count gives N_B - N_F = 0 exactly. Excluding the zero mode would instead give N_B^((nonzero)) - N_F = -1 for every D, contributing a residual -M_*²/(16π²) to rather than 0. Added to the self-loop +(3/2) M_*²/(16π²) of , this gives total δ m_h² = M_*²/(32π²), conflicting with the predicted m_h = M_* √6/(8π). The off-by-one is resolved by the substrate-vs-emergent layer distinction the body now adopts (Computations 93, 97): the substrate Walsh zero mode χ_∅ maps under Π to the constant function on M (the vacuum direction), and is structurally distinct from the emergent 4-component SM Higgs doublet generated by inner fluctuations on M × F. Counting the substrate zero mode in N_B^(full) therefore does not double-count the emergent doublet, which is carried separately by C_H = 6λ = 3/2 in . With N_B^(full) - N_F = 0 from the substrate trace and the +(3/2) M_*²/(16π²) self-loop supplying the Higgs mass, total δ m_h² = (3/2) M_*²/(16π²), recovering m_h = M_* √6/(8π) and confirming M_* ≈ 1.285 TeV as a parameter-free one-loop derivation. The three mechanisms once entertained here, (a) the Higgs doublet's three eaten Goldstones at the EW scale, (b) reinterpreting the zero mode as part of the cancellation sum, and (c) a complementation/parity pairing of the non-zero Boolean modes, are superseded by this resolution and are walked through and set aside in Computation 69.

Computation 68 — Casimir-coefficient ξ kernel-shape sensitivity [NUMERICAL]
.py.txtcomputation_68.py

Tests the sensitivity of the parameter-free Casimir-correction coefficient ξ = (15/π²) · M₄/(d₀² M₂) () to the assumed projection-kernel shape, by evaluating the second and fourth moments M₂, M₄ for six kernel families: Gaussian (book canonical), squared Lorentzian, single-sided exponential, sech-squared, compact parabolic, and compact quartic-bump. Findings: the d⁻⁶ power-law exponent is preserved across every kernel (categorical prediction, kernel-independent). The coefficient ξ varies from 2.0 (compact quartic-bump) to 21.0 (sech²) across kernels with exponential or compact-support tails (factor ∼ 10 range); the squared Lorentzian gives ξ = 378.7 from its 1/x⁴ heavy tail (physically implausible for a discreteness scale). The diagnostic d₀ bound from nanometre-Casimir data scales as 1/√ξ, giving the bound range d₀ ≲ 3.5 nm (sech²) to d₀ ≲ 11 nm (compact quartic), with the Gaussian value 5.3 nm sitting near the middle. Order of magnitude (few nm) is robust; precise numerical bound is kernel-conditional with factor-of-∼ 3 spread. Honest reading: “d₀ ≲ 5–10 nm depending on the projection-kernel choice; the Gaussian (canonical) choice gives the most constraining value 5.3 nm.”

Computation 69 — Higgs cancellation: detailed analysis of the three Computation 67 resolutions [PROOF]
.py.txtcomputation_69.py

Carries out the explicit walkthrough of all three resolutions listed in Computation 67 for the off-by-one in . (a) Higgs-doublet eaten Goldstones: the resolution is at the EW scale (broken phase), while the cancellation argument is at M_* ≫ v (symmetric phase), dimensionally a category error; not constructible without an explicit RG-matching argument the book does not provide. (b) Zero mode included in the cancellation sum: gives N_B - N_F = 0 exactly, but no leftover +M_*²/(16π²) to be the Higgs self-mass, total δ m_h² = 0, m_h = 0 at one loop; a no-op resolution. (c) Complementation S → S^c pairing: for even D, |S| and |S^c| have the same parity so complementation does not cross-pair bosons with fermions; for odd D the pairing exists structurally but the zero-mode / full-set pair is broken if the zero mode is set aside as the Higgs/order-parameter direction. None of these three resolutions, as constructible from current PST postulates without ad hoc structure, closes the off-by-one. The fourth resolution does: the substrate Walsh zero mode χ_∅ maps under Π to the constant function on M (vacuum direction), structurally distinct from the emergent 4-component SM Higgs doublet generated by inner fluctuations on M × F (Computations 93, 97). Including the zero mode in N_B does not double-count C_H = 6λ = 3/2, and M_* = 4π m_h√(2/3) stands as a parameter-free one-loop derivation.

Computation 70 — K-fixed Berry-Esseen rate for A6 [SURROGATE]
.py.txtcomputation_70.py

Verifies the Berry-Esseen-type convergence rate underlying the A6 reduction lemma: the empirical Boolean Dirichlet form converges to the continuum L²-gradient at CLT slope ≈ -0.467 on a single compact-support test family in K = [-2, 2]³ ⊂ ℝ³. Three test functions (Gaussian, C^∞ bump, polynomial-cosine). Establishes the K-fixed Berry-Esseen rate; Computation 83 extends to K-uniform.

Computation 71 — Z² ≈ e⁻¹ candidate identification via substrate spectral action [NUMERICAL]
.py.txtcomputation_71.py

Tests whether the one-loop near-coincidence Z² := λ_SM(M_*)/b(M_*) ≈ e⁻¹ (0.8% deviation) lifts to a structural identification via the substrate spectral-action ratio S_sub/2^D = e⁻¹ (exact, D-independent at matched scaling Λ = √D). Verifies the substrate side exactly for D ∈ {4, 6, 8, 10, 12, 14, 16}; computes one-loop SM RGE running λ_SM(v) → λ_SM(M_*) across {λ, y_t, g₁, g₂, g_s}; reports λ_SM(M_*) ≈ 0.087. Identifies the spectral-action Higgs-quartic identity as the open structural conjecture.

Computation 72 — Higgs perturbation and φ⁴ coefficient at matched scaling [NUMERICAL]
.py.txtcomputation_72.py

Carries out the symbolic Taylor expansion of the substrate spectral action under a Higgs-like perturbation D_φ = D_sub + y φ §igma_X at matched scaling, extracts the φ⁴ coefficient e⁻¹ · y⁴/(2 D²), and reduces the Z² identification to the single-scalar condition y/Λ = 2^(-1/4) = μ_site^(1/4) (the fourth root of the per-site Bernoulli probability of P1). Identifies the open structural content: deriving y/Λ = 2^(-1/4) from CC inner fluctuation.

Computation 73 — Casimir ξ from the substrate Boolean kernel [NUMERICAL]
.py.txtcomputation_73.py

Sharpens Computation 68's kernel-shape sensitivity by deriving ξ = (15/π²) M₄/(d₀² M₂) from the substrate's actual projection kernel (the binomial Boolean kernel) at finite |D|. Computes the moments M₂, M₄ of the Boolean kernel as binomial-sum identities; identifies the leading non-Gaussian correction; extrapolates to the matched-scaling limit. Confirms the Gaussian kernel as the natural matched-scaling limit; documents finite-|D| corrections at the O(1/D) level.

Computation 74 — First-cut CC inner-fluctuation for Z² = e⁻¹ [NUMERICAL, superseded by Computation 100]
.py.txtcomputation_74.py

Tests the simplest natural single-generator CC inner fluctuation A = iφ[D_sub, ω]/2 with ω = e₁ e₂ (quaternion volume element). Computes the perturbed Dirac D_φ² = (D + 2φ²) I and the φ⁴ coefficient 2 e⁻¹/D². Matching Computation 72 target gives y = √2 (D-independent) hence y/Λ = √(2/D), which DECREASES with D. Target is 2^(-1/4) (D-independent). Conclusion: natural single-generator inner fluctuation does NOT close Z² = e⁻¹. Identifies three possible resolutions (numerical coincidence, multi-generator A, alternative Λ(D)).

Computation 75 — Rate integral j_τ as Bernoulli large-deviation [NUMERICAL]
.py.txtcomputation_75.py

Sharpens the rate-integral closure (book ) by identifying j_τ as a standard large-deviation integral under the Bernoulli measure. j_τ ∼ (-|D| · I(τ/√(|D|))) with I the Cramér rate function. Computes saddle-point asymptotic for a representative T(C) form; documents the rate integral's structural form. Closes the rate-problem framework for downstream applications (cosmological constant, Computation 82).

Computation 76 — Postulate-level V_7-direction constraint as a Yukawa-hierarchy source [PROOF]
.py.txtcomputation_76.py

Tests whether constraining v(a) ∈ V₇ direction to the three quaternionic sub-algebras Q_k of 𝕆 containing τ̂. Result: maximal G₂-symmetry preserves S₃-symmetry across the three Q_k, so this single-direction constraint alone does not lift the structural-scope theorem. The threshold-induced dynamical mechanism is taken up in Computation 77.

Computation 77 — Dynamical Yukawa hierarchy via the directed modal threshold per Q_k [NUMERICAL]
.py.txtcomputation_77.py

The directed P3 threshold-crossing rate per quaternionic sub-algebra Q_k gives a hierarchy of generation Yukawa couplings via Boltzmann weighting (-β E_k). Reproduces the SM up-type Yukawa hierarchy with β ≈ 28.24 and (p₁, p₂, p₃) ≈ (0.52, 0.35, 0.13). Identifies the two inputs requiring structural derivation: β and the (p_k) pattern; structural candidates for each are taken up in Computations 78, 79, 80.

Computation 78 — Structural candidates for β in the Boltzmann mechanism [NUMERICAL]
.py.txtcomputation_78.py

Tests structural candidates for the Boltzmann inverse-temperature β from Computation 77's fit. Candidates: (SO(8)) = 28, (G₂) = 14, 4 (V₇) = 28, (Spin(7)) = 21. Numerical match: (SO(8)) = 28 at 0.85% from Computation 77 best-fit 28.24. Motivation: Spin(8) is the rotation group of 𝕆 = ℝ⁸, with triality permuting three 8-dim reps and three PST generations.

Computation 79 — β = 28.25 = (SO(8)) + b_PST and p_2/p_3 = e [NUMERICAL]
.py.txtcomputation_79.py

Sharpens Computation 78's β ≈ 28 to β = (SO(8)) + b_PST = 28 + 1/4 = 28.25, matching Computation 77 best-fit to 0.04%. Identifies p₂/p₃ = e in Computation 77's (p₂, p₃) ratio as structurally suggestive of the same e⁻¹ that drives Z². Locks in β = 28.25 structurally; p₁ remains as the one residual free parameter.

Computation 80 — p_1 closure via Spin(8) → SU(3)_c symmetry chain [PROOF]
.py.txtcomputation_80.py

Examines whether p₁ in Computation 77/79's dynamical mechanism can be pinned by the substrate symmetry chain Spin(8) ⊃ Spin(7) ⊃ G₂ ⊃ SU(3)_c. Result: SU(3)_c stabilizer structure forces p₁ to a T(C)-contingent value rather than a structural constant. Reclassifies the Yukawa hierarchy as configuration-content (consistent with the structural-scope theorem of §).

Computation 81 — Full CC inner fluctuation D + A + ε' J A J⁻¹ for Z² [NUMERICAL, superseded by the partition-function route, Computations 87–92/100]
.py.txtcomputation_81.py

Extends Computation 74's single-generator inner fluctuation to include the JAJ⁻¹ reality term with KO-2 sign ε' = +1, and tests multi-component A. Result: y/Λ = √(8/D) (vs Computation 74's √(2/D), factor 2 from reality doubling), still D-dependent with wrong scaling. Confirms Computation 74's verdict that single- and few-generator inner fluctuations do not close Z² = e⁻¹. Identifies the structural obstacle: the discrete spectrum's collapse-to-cutoff behaviour.

Computation 82 — Λ magnitude: order-of-magnitude reclassification as T(C)-contingent [NUMERICAL]
.py.txtcomputation_82.py

Order-of-magnitude analysis for the cosmological-constant magnitude using Computation 75's Cramér framework. ρ_inst = j_τ · · (|D|/V_M) with matched-scaling |D|/V_M ∼ 1/d₀⁴. Target j_τ for ρ_inst = ρ_Λ,obs comes out to ∼ 10⁻⁵⁰ per substrate configuration; required I(α) is small. Conclusion: Λ magnitude depends on the explicit T(C) form plus threshold-scaling τ(D), both T(C)-contingent. Reclassifies Λ magnitude from open derivation to structural contingent (joins Yukawa hierarchy, CKM, d₀). Sign + EoS w = -1 remain structural theorems.

Computation 83 — K-uniform Berry-Esseen rate verification for A6 [SURROGATE]
.py.txtcomputation_83.py

Extends Computation 70's K-fixed verification to K-uniform: seven test-function families with varying support locations, radii, and shapes (Gaussian, C^∞ bump, polynomial-cosine). Mean slope -0.458 ± 0.085 across families; 5/7 families within fitting tolerance of -0.500; two failing families have support radius ∼ d₀ (Taylor-linearization breakdown). K-uniform rate holds asymptotically in the natural physical regime (support radius ≫ d₀). Together with Computation 70 + §Closures (a)–(c), closes the reduced problem (R) at proof-detail level for the i.i.d. surrogate.

Computation 84 — Multi-generator CC inner-fluctuation route for the Z² identification [NUMERICAL, superseded by the partition-function route, Computations 87–92/100]
.py.txtcomputation_84.py

Extends Computation 74's single-generator inner fluctuation to K coherent quaternion-pair generators ω_k = e₂k-1 e₂k. Computes D_φ² = D + 2Kφ² and the φ⁴ coefficient. Matching target gives y/Λ = √(2K/D), requires K = D/(2√2) (non-integer) for D-independent target 2^(-1/4). Bernoulli-activated extension: y/Λ = √p requires p = 1/√2, not derivable from μ_site = 1/2. The multi-generator inner-fluctuation route does not deliver the target at a structurally natural value of K or p; the partition-function-level identification of Computations 87–92 takes a different substrate-side route.

Computation 85 — A_F-internal Higgs inner fluctuation for the Z² identification [NUMERICAL, superseded by the partition-function route, Computations 87–92/100]
.py.txtcomputation_85.py

Places the Higgs in the canonical A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) finite-internal factor (Connes-Chamseddine-Marcolli SM construction) rather than the Clifford volume element. The CCM moment-extraction λ ∼ (a/π²) f₀ reduces to the same identification problem at the moment-extraction level; no new computational test. Identifies the structural obstacle made explicit in Computation 86: the heat-kernel formalism presupposes a continuous spectrum, which the discrete substrate triple lacks. The partition-function-level identification of Computations 87–92 provides the substrate-side route that does not depend on the heat-kernel formalism.

Computation 86 — Seeley-DeWitt heat-kernel analysis on the discrete substrate triple [PROOF, superseded by the partition-function route, Computations 87–92/100]
.py.txtcomputation_86.py

Examines the Seeley-DeWitt heat-kernel expansion Tr f(D²/Λ²) ∼ Σ_n f_n Λ^(d - 2n) a_n for the substrate spectral triple. For the PST substrate at matched scaling: all eigenvalues at ±√D, so a₀ = 2^D (non-zero), a_n = 0 for n ≥ 1 (no continuous-spectrum curvature contributions). The moment structure that delivers λ in continuum CCM is structurally absent. Standard CCM heat-kernel machinery is therefore structurally obstructed for the discrete substrate triple, ruling out the multi-generator CC route flagged by Computation 74. Computations 87–92 take the partition-function-level identification, which does not depend on the heat-kernel formalism.

Computation 87 — Bernoulli Laplace-transform identification: bypassing the CC inner-fluctuation route [NUMERICAL]
.py.txtcomputation_87.py

Bypasses CC inner fluctuation via the asymptotic Bernoulli moment-generating function. E_μ[ (c X̄)] = ((1 + e^(c/D))/2)^D → (c/2) as D → ∞. Hits e⁻¹ at c = -2. Numerical convergence verified at D = 10⁴ with error ∼ 10⁻⁵. Multiple structural readings of c = -2 (e.g. -1/μ_site, number of substrate states per bit) compatible; none uniquely forced. Computation 88 fixes the load-bearing reading as the binary-distinction Cl(1,0) gap 2 = 1-(-1).

Computation 88 — binary-gap-tempered Bernoulli MGF identification of e⁻¹ at the matched scaling [PROOF]
.py.txtcomputation_88.py

The scale c = -2 is structurally forced by the binary-distinction gap: each substrate exchange τ_i is a self-adjoint involution with eigenvalues {+1,-1}, so the one-bit Clifford Cl(1,0) fixes β_KO = 1-(-1) = 2 with no freedom (P1's binary alphabet; Computations 100 and 102). Setting c = -β_KO = -2 gives e⁻¹ = E_μ[ (-β_KO X̄)] asymptotically. The total KO-dimension KO_total = 4+6 = 10 ≡ 2 8 shares this value but is structurally independent (an additive residue in ℤ/8ℤ, not a spectral width); the coincidence is retained only as the “binary-gap-tempered” label, never as an input. All three ingredients structurally derived (μ_site from P1, the binary-distinction gap 2 from P1, asymptotic limit from Cramér-Bernoulli). Second independent substrate-side derivation of e⁻¹; complements Computation 62's spectral-action route.

Computation 89 — Z² bridge step 1: substrate Higgs Hamiltonian [PROOF]
.py.txtcomputation_89.py

Derives the substrate Higgs Hamiltonian H_Higgs(C) structurally from three constraints: (i) additivity (Bernoulli sites independent ⇒ H = Σ_a h_a(B_a)), (ii) matched scaling (per-site energy 1/D), (iii) uniformity (P1 symmetry ⇒ h_a uniform). Together force H_Higgs(C) = (1/D) Σ_a B_a = X̄(C), no free parameters. binary-gap-tempered partition function Z_H(β_KO) = E_μ[ (-2 X̄)] → e⁻¹ asymptotically. Substrate side of bridge closed structurally. Numerical verification at D = 10, 100, 1000.

Computation 90 — Z² bridge step 2: three SM-side candidate formulations [EMPIRICAL]
.py.txtcomputation_90.py

Reduces the SM-side bridge to a choice among three concrete candidate formulations: (I) Wilsonian effective-coupling RG, integrate substrate fluctuations between M_* and Λ; (II) spectral-action variation, naive form gives e⁻² rather than e⁻¹, fails; (III) substrate-measure invariance, V_eff(C)/T_KO = β_KO X̄(C) at matched scaling. Numerical: predicted λ_SM(M_*) = e⁻¹/4 ≈ 0.0920 vs observed 0.0927 (Buttazzo 2013), ratio 1.008. Formulation (II) does not close in naive form; (I) and (III) are taken up as concrete research directions in subsequent work.

Computation 91 — Z² empirical-precision test: two-loop SM RGE [EMPIRICAL, superseded by Computation 123]
.py.txtcomputation_91.py

Tests the empirical tightness of Z² ≈ e⁻¹ via SM RGE precision. Simplified one-loop running gives λ_SM(M_*) ≈ 0.087 (5% below target); compact two-loop ≈ 0.085 (7% below). The book's 0.8% near-coincidence claim depends on full Buttazzo-level SM RGE precision (two-loop running with explicit threshold corrections at m_t). Sensitivity analysis: variation of m_h ± 1 GeV, m_t ± 2 GeV gives Z²/e⁻¹ in range 0.89–0.97 at simplified one-loop. All readings consistent with structural identity within SM RGE truncation uncertainty. Substrate-side derivation (Computation 89) remains exact and RGE-truncation-independent. Superseded by Computation 123, which performs the proper two-loop running with a full input-uncertainty budget, λ_SM(M_*) = 0.0927 ± 0.0007 (m_t/m_h-dominated), a rigid 1.0σ match to e⁻¹/4 = 0.09197 that converges monotonically with loop order, confirming the simplified 0.085–0.087 readings here were truncation artefacts, exactly as anticipated.

Computation 92 — Z² bridge formalisation: seven-step structural proof sketch [PROOF]
.py.txtcomputation_92.py

Reduces the Z² = e⁻¹ closure to a single Wilsonian-matching premise via a seven-step structural derivation: (1) PST EFT below M_* = pure SM (); (2) Wilsonian matching at M_*; (3) bare λ_PST^bare = b = 1/4 (§); (4) substrate UV at Λ = √D contributes vacuum fluctuations; (5) fluctuations integrated over μ give Z_H(β) = E_μ[ (-β H_Higgs)] with H_Higgs = X̄ (Computation 89); (6) β_KO = KO_total 8 = 2 (Computation 88); (7) λ_SM(M_*) = b · Z_H(β_KO) → e⁻¹/4 structurally, matching observed 0.0927 at 0.8%. Steps (1)–(4), (6) established in existing PST framework; step (5) by Computation 89; step (7) follows. Elevates Z² closure status from “open conjecture” to “structural proof sketch established.” The seven steps are subsequently formalised in the partition-function-level Chamseddine-Connes correspondence framework of §; the substrate-side content of Bridge Premise (B) is closed at proof-detail level from P1–P3 by Computation 100 (matched-scaling cutoff identification via tensor-product factorisation forced by P1's Bernoulli product measure) and Computation 98 (mode-shell residual via the discrete Wilson block-spin RG fixed-point theorem); the matching identification of the substrate factor with the SM-side coupling ratio (bridge premise B itself) is subsequently reduced within PST to a field-strength-renormalisation identity (substrate side closed; SM-side matching open) by the wave-function-renormalisation route and the unique-soft-mode theorem (Computations 110–116).

Computation 93 — M_* = 4π m_h√(2/3) off-by-one: fourth-resolution candidate via substrate-vs-emergent distinction [PROOF]
.py.txtcomputation_93.py

Re-examines resolution (b) of Computation 69 under a foundational reading: the substrate-side Boolean zero mode S = ∅ is the substrate's constant-mode profile, NOT the emergent SM Higgs doublet (captured separately by C_H = 6λ_PST = 3/2 from inner fluctuations on M × F). Under this substrate-vs-emergent distinction, including the zero mode in N_B does NOT double-count the doublet, and total δ m_h² = (M_*²/16π²)(C_H + N_B^full - N_F) = (3/2) M_*²/(16π²), recovering m_h = M_* √6/(8π) as parameter-free one-loop derivation. Identifies the fourth-resolution candidate; full validation via projection-chain analysis follows in Computation 97.

Computation 94 — Bridge Premise (B) fluctuation-matching T_KO (research-direction sketch) [SURROGATE, superseded by Computations 95, 100]
.py.txtcomputation_94.py

Attempts to derive Bridge Premise (B) via formulation III with fluctuation-matching prescription ⟨(X̄ - 1/2)²⟩_μ = ⟨ h²⟩_thermal,T, giving the matched-scaling temperature T_KO = μ²(M_*)/(4D). Combined with β_KO = 2 delivers the testable prediction μ²(M_*) = 2D · M_*². Acknowledged as dimensionally suspect (X̄ dimensionless vs h with GeV units; the prediction implies hierarchy violation μ ≳ M_* for D ≳ few). Held as research-direction sketch rather than a closure; superseded by Computation 95 reduction and ultimately Computation 100 closure.

Computation 95 — Bridge Premise (B) reduced to structural claim (***) via substrate-Wilsonian RG [PROOF]
.py.txtcomputation_95.py

Reduces Bridge Premise (B) to a single sharper structural claim (***): the substrate Wilsonian RG flow generator equals the Boltzmann partition function (-β_KO · H_Higgs(C)) on substrate configurations, the partition-function-level analogue of the heat-kernel RG flow generator (-t D²) for the continuum spectral triple. If (***) is granted, the bridge identification follows structurally: λ_SM(M_*) = b · Z_H(β_KO) = e⁻¹/4 ≈ 0.0920. Provides the residual-gap framing subsequently closed by Computations 98 (mode-shell) and 100 (matched-scaling cutoff).

Computation 96 — Polchinski-discrete derivation of claim (***): partial advance [NUMERICAL]
.py.txtcomputation_96.py

Polchinski-style derivation attempted on the discrete substrate. The partition-function form is structurally established via configuration-Hamiltonian Boltzmann averaging Z_D(β_KO) = ((1 + e^(-β_KO/D))/2)^D → e⁻¹ as D → ∞. Identifies three rigorous-proof items as residual gap: (a) formal exact RG flow equation on discrete substrate triple, (b) uniqueness of Boltzmann form at matched scaling, (c) cutoff-function identification with asymptotic D → ∞ limit. Partial advance; (a) and (b) closed by Computation 98, (c) by Computation 100.

Computation 97 — M_* off-by-one: projection-chain validation CLOSES the item [PROOF]
.py.txtcomputation_97.py

Validates Computation 93's substrate-vs-emergent distinction via explicit PST projection-chain analysis ψ → Π(ψ) → H_doublet. The substrate Walsh zero mode χ_∅ maps under Π to the constant function on M = ℝ × S³ (vacuum direction), structurally distinct from the emergent SM Higgs doublet (4 components from inner fluctuations on A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ)). Including the zero mode in N_B does NOT double-count C_H = 6λ = 3/2. Total δ m_h² = (3/2) M_*²/(16π²), restoring m_h = M_* √6/(8π) as parameter-free one-loop derivation. Computation 67/69 off-by-one CLOSED at proof-detail level; the “one-loop-complete / parameter-free” reading of M_* = 4π m_h√(2/3) is restored.

Computation 98 — Bridge Premise (B) items (a)+(b) CLOSED via mode-shell block-spin RG [PROOF]
.py.txtcomputation_98.py

Closes items (a), (b) of Computation 96 at proof-detail level via explicit discrete Wilson block-spin RG. Key observation: H_Higgs = X̄ is supported on the lowest two Walsh-mode shells V₀ ⊕ V₁ only, with explicit expansion X̄ = (1/2)χ_∅ - (1/(2D))Σ_|S|=1χ_S. Because X̄ has no |S| 2 content, the block-spin RG step from level m+1 to m factorises trivially: X̄(ψ_ m + ψ_m+1) = X̄(ψ_ m) for every m 1, so the Boltzmann action β_KOX̄ is a TRIVIAL FIXED POINT (Theorem 1). Uniqueness follows from Computation 89 boundary condition. Residual item (c) remains as the matched-scaling identification (subsequently closed by Computation 100).

Computation 99 — Item (c) via thermal-QFT framing (EXPLORATORY, superseded by Computation 100) [SURROGATE, superseded by Computation 100]
.py.txtcomputation_99.py

Exploratory closure of item (c) via thermal-field-theory framing: postulates β_KO = KO 8 as substrate's intrinsic inverse temperature in the partition-function-level CC framework. Standard thermal-QFT IR coupling formula λ_eff = λ_bare · Z(β) delivers λ_SM(M_*) = b · Z_H(β_KO) = b · e⁻¹. The “KO-thermal principle” is at the same level as the spectral-action principle in standard CC, a framework postulate beyond P1–P3. Held as exploratory, not as a P1–P3-only closure: PST's claim is that all structure derives from P1–P3 alone, so a CC-framework postulate weakens that claim. SUPERSEDED by Computation 100, which closes (c) from P1–P3 alone via tensor-product factorisation, eliminating the thermal-postulate requirement.

Computation 100 — Substrate-side spectral-action cutoff identification closed from P1–P3 via tensor product [PROOF]
.py.txtcomputation_100.py

Closes the substrate-side content of item (c) of Computation 96 from P1–P3. P1's Bernoulli independence forces tensor product L²(P(D), μ) = _i L²({0,1}, μ_i). The Boolean Laplacian Δ = Σ_i (1 - τ_i) with one-bit Clifford Cl(1,0) gives each (1 - τ_i) eigenvalues {0, 2}. For Tr f(Δ/Λ²) to be CONSISTENT with the tensor product structure (i.e., to factor as a product of one-bit traces over independent bits) f must satisfy f(x+y) = f(x) f(y), unique smooth solution f(x) = (-x). Forced by P1's product structure, not postulated. At matched scaling Λ² = D: (1/2^D) Tr e^(-Δ/D) = ((1 + e^(-2/D))/2)^D → e⁻¹ by binomial product convergence. The factor “2” in (-2/D) comes from one-bit Clifford structure, NOT KO mod 8 thermal interpretation (numerically equal but structurally independent). Computation 99's KO-thermal postulate SUPERSEDED. The matching identification of this substrate-side factor with the SM-side coupling ratio λ_SM(M_*)/b (bridge premise (B) itself) is subsequently reduced within PST to a field-strength-renormalisation identity (substrate side closed; SM-side matching open) by the wave-function-renormalisation route (Computations 110–116); λ_SM(M_*) = b · e⁻¹ = 0.0920 matches observed 0.0927 at 0.8% Buttazzo-RGE precision (5–7% at one-loop, Computation 91).

Computation 101 — Bridge premise (B) sharpened: Z² as wave-function renormalisation Z_φ² [PROOF]
.py.txtcomputation_101.py

Addresses the structural concern that standard Wilsonian threshold matching of a quartic is additive in loop corrections, not multiplicative. Reinterprets the multiplicative matching λ_SM(M_*) = b · Z_H(β_KO) as a wave-function renormalisation of the substrate modal field ψ to the SM Higgs φ at matched scaling, with the substrate-side partition function Z_H(β_KO) → e⁻¹ playing the role of Z_φ² via standard QFT decomposition Z² = Z_Γ 4/Z_φ². Reduces (B) to two well-posed sub-questions: (1) is the matched-scaling Π's path-integral Jacobian equal to (-β_KO H_Higgs(C))? (2) why is the vertex renormalisation trivial (Z_Γ 4 = 1)? Both are tractable extensions of Computation 73's matched-scaling machinery; closing both delivers (B) automatically via standard renormalisation. Status: structural reduction of (B)'s open content from a non-standard hypothesis to two standard renormalisation questions. Both sub-questions are subsequently closed: the multiplicative (not additive) reading is the correct one (the additive-Jacobian route is excluded, Computation 103) and the vertex renormalisation Z_Γ 4 = 1 holds to all orders (Computation 119); (B) is thereby reduced (Computations 110–119).

Computation 102 — Item 1.2 resolution: the two “2”s are structurally independent [PROOF]
.py.txtcomputation_102.py

Addresses the concern that the identification c = -(KO_total 8) = -2 in the Bernoulli MGF is asserted, not motivated. Computation 100 already derived the “2” from one-bit Clifford structure Cl(1,0) (eigenvalue range of (1-τ_i)), independent of KO mod 8. This computation makes the structural-independence explicit: KO_total 8 = 2 (from spacetime + internal split 4 + 6 = 10) and one-bit Clifford eigenvalue range = 2 (from τ_i² = 1) trace to different parts of the substrate construction and are numerically coincident but structurally unrelated. The honest framing is that Computation 100's derivation uses ONLY one-bit Clifford structure; KO mod 8 = 2 is a parallel substrate fact, not an input to the matched-scaling exponent. Resolves item 1.2 as a derivation question; the remaining open content (whether partition-function-level CC has standing in the spectral-action literature) is an external-validation / literature-search task, not new research.

Computation 103 — Bridge premise (B) attack: explicit path-integral Jacobian, honest negative result [NUMERICAL, superseded by Computations 110–116]
.py.txtcomputation_103.py

Attempts Computation 101's sub-question (1): is the matched-scaling Π's path-integral Jacobian equal to (-β_KO H_Higgs(C))? Two natural interpretations both fail: (A) the linear-map determinant of Π : V₀ ⊕ V₁ → SM_low is configuration-independent (constant); (B) the Cramér-Bernoulli LDP near the vacuum X̄ = 1/2 gives a Gaussian weight in (X̄ - 1/2), not exponential in X̄ (numerically verified I(x) = x (2x) + (1-x) (2(1-x)) ≈ 2(x-1/2)² near 1/2). The multiplicative matching λ_SM(M_*) = b · Z_H(β_KO) is therefore NOT derivable from a path-integral change-of-variable. Confirms that the Jacobian route to (B) does not close. Sub-question (2) [Z_Γ 4 = 1] is satisfied at tree level (standard CC inner fluctuations). Net: this is a negative result for the measure-Jacobian reading specifically; it does not obstruct (B), because (B) is a wave-function renormalisation (multiplicative), not a coupling-matching Jacobian (additive). The Gaussian found here is precisely what the product-structure argument later excludes for the field-strength kernel (Computation 111), and (B) is subsequently reduced within PST to a field-strength-renormalisation identity (substrate side closed; SM-side matching open) by the wave-function route (Computations 110–116).

Computation 104 — Bridge premise (B) factored: substrate side derived, SM side inherited from standard CC literature [PROOF, partially retracted, superseded by Computations 110–119]
.py.txtcomputation_104.py

Attempts to extend Mosco convergence from Dirichlet forms (which sec:mosco-conditional establishes) to spectral actions, with the goal of deriving bridge premise (B) at the same level Mosco delivers the SM kinetic term. The structural analysis shows (B) FACTORS into two parts: (B-substrate) (1/2^(|D|)) Tr f(Δ/|D|) → e⁻¹, which is CLOSED at proof-detail level from P1–P3 by Computation 100; and (B-SM) the SM-side spectral-action identification of the matched-scaling coefficient with λ_SM(M_*), which IS the standard Chamseddine-Connes spectral-action principle applied to PST's substrate spectral triple. This is an ESTABLISHED FRAMEWORK with 30+ years of literature precedent, NOT a PST-specific postulate. PST extends standard CC by providing a precausal foundation BENEATH the spectral triple (P1–P3 deliver the substrate spectral triple itself, Computations 3, 19); it does not invent the SM-side framework. Whether the partition-function-level CC framework has standing in the spectral-action literature: YES, the framework PST uses IS standard CC, with the partition-function-level adaptation being PST's specific way of handling the discrete substrate where heat-kernel collapses (Computations 85–86). PST's foundational claim refines to “P1–P3 + standard CC spectral-action framework”: the substrate-side contribution is the genuinely new PST result (Computation 100, derived from P1–P3), and the SM-side framework is inherited from established CC literature. This inheritance framing is subsequently superseded: the matching (B) is reduced internally to PST, in the substrate limit, to a field-strength-renormalisation identity whose internal content is fully derived (field-strength and vertex, Computations 110–119), without inheriting an SM-side CC postulate; the foundational claim is then “P1–P3 + matched scaling”, with the bridge resting on the emergence premise F4 + the RGE test rather than an inherited framework.

Computation 105 — Two-layer one-loop diagram (no double counting of M_* zero modes) [NUMERICAL]
.py.txtcomputation_105.py

Addresses the concern of potential double-counting between the substrate χ_∅ singleton-cylinder order parameter and the emergent SM Higgs field C_H from Mosco convergence. The PST partition function factorises as Z_PST(M_*) = Z_sub(M_*, Λ=√D) · Z_EFT(M_*, IR), with the two layers integrating over disjoint momentum ranges [M_*², D] and [m_h², M_*²] that meet only at the single matching scale M_* (measure zero in both integrals). Three independent non-double-counting criteria are established: (i) disjoint momentum support; (ii) distinct dynamical variables (χ_∅ on P(D) vs C_H on M = ℝ × S³); (iii) distinct counterterm structures (the matching is a boundary condition, not a loop integration). Numerical verification: the substrate full-trace partition function (Computation 100) decomposes exactly as Z^full_sub(D) = Z^IR_sub(D, k_*) + Z^UV_sub(D, k_*) across the M_*-shell partition with |Z^full - (Z^IR + Z^UV)| < 10⁻¹⁴ for D ∈ {10, 20, 50, 100}. Closes the double-counting question at the diagram level.

Computation 106 — SM gauge block-diagonality on ⊕_k G_k [PROOF]
.py.txtcomputation_106.py

Provides the explicit proof that SM gauge generators preserve each generation sector in the Dixon-synthesis decomposition V_Dixon = V_core ⊕ G₁ ⊕ G₂ ⊕ G₃. Three-part proof: (1) SU(3)_c is the M₃(ℂ) factor of A_F acting on the colour triplet WITHIN one Fano-plane sub-algebra H_k of 𝕆 through τ̂, hence within a single G_k; (2) SU(2)_L is the ℍ factor acting on the quaternionic slot of V_Dixon = ℂ ⊗ ℍ ⊗ 𝕆, orthogonal to the octonionic generation labelling, hence block-diagonal by tensor-product structure; (3) U(1)_Y is the ℂ-factor scalar action, same tensor-product argument. Numerical verification on a 12-dim toy model V_gen ⊗ V_per-gen: SU(3)_c, SU(2)_L, U(1)_Y all give machine-zero cross-generation matrix elements; a hypothetical “gauge on generation slot” control gives off-block magnitudes of O(1). Concludes that the absence of tree-level flavour-changing neutral currents is a structural theorem of PST conditional on the Dixon-synthesis block structure, not a postulate.

Computation 107 — Full-SM extension of the M_* boson-fermion cancellation: sector-blindness sub-result [NUMERICAL]
.py.txtcomputation_107.py

The substrate-level binomial cancellation N_B = N_F = 2^(D-1) is SECTOR-BLIND (invariant under any partition of the 2^D Walsh modes; numerically verified across D ∈ {4, …, 20}). This sector-blindness does not transfer the cancellation to the SM top/W/Z sectors: by the disjoint-layer structure of Computation 105 the substrate cancellation lives in the UV shell [M_*, √D] while the physical top/W/Z loops run in the emergent EFT shell below M_*, so the substrate identity cannot cancel the EFT-layer loops. The “coupling universality at M_*” premise (c_S = 1 for every Walsh mode) is contradicted by the observed Yukawa hierarchy, which the book assigns to T(C)-contingency (§). The empirical SM Veltman residual at M_* is ≈ -3 (top-dominated, negative), not the +m_h² a clean cancellation would leave. The sector-blindness of the substrate count is a fact; it does not close the full-SM extension, which does not exist within PST (Computation 120), so M_* is irreducibly a scalar-sector prediction.

Computation 108 — Inner-fluctuation sign assignment for SM gauge bosons is BOSONIC (sub-result; does not close 1.3) [PROOF]
.py.txtcomputation_108.py

Verifies that the substrate |S|-parity boson/fermion sign assignment (CAR realisation, book sec:car-fermions) carries through to A_F inner-fluctuation gauge bosons that are not themselves Walsh modes. Three independent readings agree: (R1) algebraic, Connes 1-forms A = Σ_i a_i [D, b_i] for a_i, b_i ∈ A_F are operators of EVEN Clifford grade; (R2) substrate-side lift, A_F is generated by Cl_even(P(D)) under the CAR realisation, so every inner fluctuation lifts to an |S|-even substrate configuration; (R3) one-loop diagrammatic, the standard QFT vector-boson loop sign (+) is inherited unchanged by the CC inner-fluctuation construction. Dof count: bosonic inner-fluctuation dofs on A_F (off-shell) are U(1)_Y: 4, SU(2)_L: 12, SU(3)_c: 32, Higgs doublet: 4, total 52; the substrate |S|-even Walsh subspace has dimension 2^(D-1) accommodating these at D ≥ 6. Status: the bosonic sign assignment is a correct sub-result, but it does NOT close open-research item 1.3. The full-SM cancellation fails for the layer-disjointness reason of Computation 107 (the substrate cancellation does not reach the EFT-layer top/W/Z loops), independent of the inner-fluctuation sign. Item 1.3 was subsequently settled in the negative by Computation 120 (the full-SM extension does not exist within PST); this sign sub-result stands.

Computation 109 — Higgs-quartic normalisation + O(4) one-loop self-energy: two technical verifications [PROOF]
.py.txtcomputation_109.py

Verifies two technical findings on the Higgs-sector derivation at proof-detail level. (A) Normalisation: the book commits to the doublet convention, λ_PST = 14 is the coefficient of (Φ^†Φ)², the same convention as the SM V_SM = λ_SM(Φ^†Φ)², with ⟨ψ,ψ⟩ ↔ Φ^†Φ. Under the radial reduction ⟨ψ,ψ⟩ → ψ² the postulate F = ∫[-ε⟨ψ,ψ⟩ + (1/4)⟨ψ,ψ⟩² + c|∇ψ|²]dμ collapses to the scalar form -εψ² + (1/4)ψ⁴ (). The field-normalisation ratio Z² = λ_SM(M_*)/λ_PST is then a ratio of two same-convention couplings. (B) Vertex contraction: explicit O(N) vertex contraction V_ijkl = 2λ(δ_ijδ_kl + δ_ikδ_jl + δ_ilδ_jk) closed into the one-loop seagull (1/2)V_ijklδ_kl = (N+2)λ δ_ij delivers the standard O(N) result (N+2)λ rather than Nλ. For N=4 the Higgs self-loop coefficient is C_H = 6λ_PST = 3/2, giving M_* = 4π m_h √(2/3) ≈ 1285 GeV. Broken-phase cross-check: 3λ (physical Higgs loop) + 3λ (three Goldstones) = 6λ confirms the symmetric-phase result. This computation propagates corrections through , the M_* value, the naturalness coefficient, and the affected Computations 67, 69, 93, 97, 107, 108.

Computation 110 — Bridge premise (B): wave-function-renormalisation route [PROOF]
.py.txtcomputation_110.py

Recasts (B) as a field-strength (wave-function) renormalisation Z_φ², which acts multiplicatively on the quartic (λ → λ Z_φ²), the form (B) requires, sidestepping the additive-matching obstruction Computation 103 found for the coupling-Jacobian reading. The route first sets up the discrete Wilson–Polchinski flow on the Boolean lattice with the P1-forced cutoff f = and matched scaling Λ² = D, reducing to the low Walsh shells V₀ ⊕ V₁ where X̄ lives (Computation 98). It then verifies that the field-strength factor is Z_φ² = Z_H(β_KO) → e⁻¹, so λ_SM = b Z_φ² → e⁻¹/4. Sets up the route; the crux (the form and value of the kernel) is closed in Computations 111–116.

Computation 111 — Bridge premise (B): the field-strength kernel is forced into (-βX) form [PROOF]
.py.txtcomputation_111.py

The field-strength dressing must respect P1's product measure (the matched-scaling cutoff factorises over bits, Computation 100), so the per-configuration kernel K(C) = ∏_a g(C_a) is a product over bits. A product kernel over Boolean bits is necessarily of a quantity linear in X̄, i.e. ∝ (-βX̄), the same form the substrate side evaluates to e⁻¹. The Gaussian (-α(X̄-1/2)²) that obstructs the Jacobian route (Computation 103) is thereby excluded: it correlates the bits, breaking the product structure (verified: bit–bit covariance 0 under (-2X̄), nonzero under the Gaussian). The wave-function route is therefore unobstructed, in contrast to the Jacobian route.

Computation 112 — Bridge premise (B): β_KO = 2 is the binary-distinction gap [PROOF]
.py.txtcomputation_112.py

The load-bearing exponent is the spectral range of a single binary distinction, 2 = 1 - (-1): each exchange τ_i is a self-adjoint involution with spectrum {+1,-1}, so (1-τ_i) has nonzero eigenvalue 2, forced by P1's binary alphabet with no freedom (equivalently the one-bit Clifford Cl(1,0) range). q-ary comparison: only q = 2 gives the self-adjoint involution P1 requires, and e⁻¹ is the signature of the binary gap (Z_H → e^(-β/2) lands on e⁻¹ only at β = 2). Dissolves the “two 2's” question: KO_total 8 = 2 shares the value but is a gauge/spacetime-sector fact, never an input.

Computation 113 — Bridge premise (B): the Z_φ ↔ trace identification, total-reduction residual [PROOF]
.py.txtcomputation_113.py

Tests the field-strength↔matched-scaling-trace identification ab initio. The order-parameter two-point residue (single |S|=1 shell) gives a single-mode factor, not the full e⁻¹, unless all substrate modes are folded in. The matched-scaling trace and the tempered order-parameter expectation E_μ[ (-2X̄)] are verified exactly equal, so the identification holds iff the substrate→Higgs reduction is total. Isolates the residual as one concrete property (total reduction), closed by Computations 114–116.

Computation 114 — Bridge premise (B): the standard-EFT normalisation is excluded; embedding-norm reading [PROOF]
.py.txtcomputation_114.py

Two field-theoretic normalisations of the reduction: (A) the total-reduction embedding norm of the symmetric Higgs vacuum = E_μ[ (-β_KOX̄)] → e⁻¹, verified exactly equal to the tempered vacuum amplitude; (B) the standard-EFT two-point residue Var_μ(X̄) = 1/(4D) → 0. The bridge needs (A); route (B), the Higgs as a light field whose two-point function is corrected by heavy loops, gives λ_SM → 0 and is excluded by the finiteness of the observed quartic. The Higgs is an emergent projected field, normalised by its embedding norm, not a fundamental light field.

Computation 115 — Bridge premise (B): unique-soft-mode theorem, total reduction [PROOF]
.py.txtcomputation_115.py

Past the LG threshold (P3) the order parameter X̄ is the unique soft mode. Lemma: X̄ = (1/D)Σ_a C_a is linear in the bits, so the LG Hessian is rank one (the all-ones matrix), softening exactly the S_D-singlet direction (= X̄) and leaving the D-1 standard-rep modes and all |S| 2 shells at the binary gap. Verified: rank = 1, one soft eigenvalue, D-1 pinned at the gap; the threshold sweep softens X̄ alone. Establishes total reduction (Theorem F10, Gaussian order), so the bridge factor is the full embedding norm e⁻¹.

Computation 116 — Bridge premise (B): all-orders decoupling theorem [PROOF]
.py.txtcomputation_116.py

Upgrades the unique-soft-mode result to all orders. Because the LG potential is a function of the linear collective coordinate X̄ alone, every derivative tensor of V is the totally symmetric all-ones tensor V⁽ⁿ⁾/Dⁿ, with zero matrix element on the non-collective subspace at every order: the standard-rep and |S| 2 modes are free fields, decoupled, pinned at the binary gap, with no V-vertex from which any correction could be built. Verified exactly with the full (-V): the V-induced bit correlation is O(1/D), the standard channel sits at the free binary-gap value with O(1/D) deviation, and only the singlet X̄ responds to V. The spectral gap is protected non-perturbatively in the substrate limit; the residual is the standard O(1/D) finite-size correction.

Computation 117 — Bridge premise (B): the embedding-norm selection is structural, not empirical [PROOF]
.py.txtcomputation_117.py

Removes the empirical crutch in the field-strength selection. The two-point-residue reading sets Z_φ² = Var_μ(X̄) = 1/(4D) → 0; that vanishing is the law of large numbers, X̄ is an intensive average of the i.i.d. bits (P1), so it self-averages to a sharp value, and a quantity the LLN sends to zero cannot be a field-strength (it gives a free emergent Higgs, λ_SM → 0). The embedding norm is the O(1) object → e⁻¹, the only candidate that normalises a canonical emergent field. Since the emergent EFT below M_* is the interacting Standard Model (F4), an interacting Higgs is required, hence Z_φ² = O(1), hence the embedding-norm reading, forced by P1's LLN and F4, not by the measured value λ_SM = 0.0927, which is demoted from an input to the selection to a test of the prediction b e⁻¹. Non-circular: F4 selects the reading, the substrate gives the value e⁻¹. Physically, below the LG threshold ⟨X̄⟩ is a sharp classical condensate and the physical Higgs is the radial excitation, normalised by the embedding norm, not the order parameter's vanishing susceptibility.

Computation 118 — Bridge premise (B): F10 decoupling as a formal all-orders S_D-equivariant Morse–Bott lemma [PROOF]
.py.txtcomputation_118.py

States and proves the F10 decoupling as a formal lemma. Because F_int = Φ(X̄) factors through the linear collective coordinate, the chain rule gives ∂ⁿ F_int = Φ⁽ⁿ⁾(X̄)/Dⁿ times the all-ones tensor at every order. Under S_D the fluctuation space splits into the singlet line (carrying X̄) and the standard representation W = {v : Σ_a v_a = 0}; an all-ones tensor contracted with any v ∈ W gives (Σ_a v_a)(⋯) = 0, so the interaction has zero matrix element on W at every order and the W-Hessian is the bare Boolean Laplacian, bounded below by the binary gap 2. Schur's lemma (S_D-equivariance) forbids singlet–W mixing and any invariant deformation of the W-block. This is the all-orders S_D-equivariant Morse–Bott structure: a soft critical direction (the order parameter) and a uniformly gapped normal bundle. Verified: derivative tensors index-independent at n = 2, 3, 4; annihilation of W to machine precision; the singlet softens 4 → 0 under the LG curvature while W stays pinned at 2; a random S_D-invariant Hessian has W-eigenvalue independent of the singlet coupling. Residual: the O(1/D) finite-size sum-constraint correction (Computation 116), vanishing as D → ∞. SUBSUMED by Computation 121: the decoupling is closed exactly and non-perturbatively by sufficiency of the empirical mean (this Morse–Bott lemma is its de Moivre–Laplace shadow), with the O(1/D) isolated to the value; Computation 122 shows F8's linearity controls only that O(1/D).

Computation 119 — Bridge premise (B): vertex renormalisation Z_Γ 4 = 1 to all orders, completing the matching's internal content [PROOF]
.py.txtcomputation_119.py

The matching ratio Z² = λ_SM/b decomposes into a field-normalisation piece (the embedding norm → e⁻¹, Computation 117) and a vertex piece Z_Γ 4 (the four-point 1PI correction). Computation 101 fixed Z_Γ 4 = 1 only at tree level; this computation closes it to all orders. The matching integrates out the substrate UV modes (the non-collective S_D standard representation W and the gapped |S| 2 shells); by the Morse–Bott decoupling lemma (Computation 118), V = Φ(X̄) has zero matrix element on W at every order, and the n=4 case is the Higgs quartic vertex, so the integrated-out modes do not couple to it and generate no vertex correction. Verified: the quartic all-ones tensor contracted on W legs is machine-zero. Hence Z_Γ 4 = 1 to all orders and Z² = e⁻¹ with no residual internal renormalisation. The only quartic loop corrections are the collective Higgs self-loops, the SM RGE running below M_*, which confirms b e⁻¹ = 0.0920 against the two-loop-run 0.0927 at 0.8% (quantified as a rigid 1.0σ match by Computation 123), not a Z_Γ 4 correction. Net: the matching's internal renormalisation content is complete; the open content of (B) reduces to the emergence premise F4 plus the RGE test, not an internal step.

Computation 120 — Full-SM extension of the M_* cancellation settled in the negative (open-research item 1.3) [EMPIRICAL]
.py.txtcomputation_120.py

Answers the open question of item 1.3 (does an EFT-layer top/W/Z cancellation exist, and what partner spectrum would it need?): no, not within PST. The SM Veltman bracket [6λ + (9/4)g² + (3/4)g'² - 6y_t²] is -3.5 at the EW scale and -3.3 at M_*, top-dominated, with no zero crossing, M_* is not a Veltman-natural full-theory scale. Cancelling the dominant -6y_t² requires bosonic top-partners (stop-like) at top strength near M_* ≈ 1.3 TeV; such a sector is collider-disfavoured (LHC excludes scalar top-partners / vector-like quarks to ∼ 1.2–1.4 TeV, i.e. at M_*), forbidden by the “no new light states” premise (§), absent from PST's spectrum (F4: emergent EFT = SM, no superpartners), and unreachable by the substrate mechanism (Computation 105: layer-disjoint). Conclusion: M_* = 4π m_h√(2/3) is irreducibly a scalar-sector prediction, exactly as the book scopes it. This is a negative settlement, a definitive limitation on the scope of M_*, not a positive closure; it removes item 1.3 from the open register with M_* correctly and finally scoped.

Computation 121 — Bridge premise (B): the all-orders decoupling closed exactly by sufficiency of the empirical mean (subsumes Computation 118) [PROOF]
.py.txtcomputation_121.py

Closes the one residual rigour task carried by F10/Computation 118 (the all-orders S_D-equivariant decoupling) by an exact, non-perturbative identity. Theorem F8 makes H_Higgs = X̄ exactly the empirical mean of i.i.d. per-site bits, so the interaction Φ(X̄) is a function of a sufficient statistic. By Fisher–Neyman factorisation the configurations at fixed X̄ = k/D are counted by the interaction-independent multiplicity Dk, and the substrate partition function factorises exactly at every finite D: Z[Φ] = Σ_k Dk 2^(-D) e^(-Φ(k/D)). The transverse (non-collective) modes enter only through the Φ-independent multiplicity, so Z_Γ 4 = 1 (and every n-point Φ-vertex correction) vanishes exactly, with no “higher-order” gap (the factorisation holds at all orders simultaneously). Verified: the factorisation to machine precision for linear/quartic/mixed Φ; the quartic-probe vertex test (Z_Γ 4 = 1 to 3×10⁻¹⁶) with a transverse-coupled control that breaks sufficiency (gap 0.30, the test has teeth); the O(1/D) isolated to the value, the decoupling carrying no D-dependence; and the de Moivre–Laplace reduction showing Computation 118's Morse–Bott lemma is the Gaussian shadow. Net: the substrate side of (B) carries only the standard O(1/D) caveat on the value e⁻¹; the single open content of (B) remains the SM-side matching (F4 + RGE; quantified as a rigid, D-independent 1.0σ test by Computation 123).

Computation 122 — Bridge premise (B): what the substrate side needs of F8 linearity (decoupling = symmetry; value = vacuum normalisation) [PROOF]
.py.txtcomputation_122.py

Isolates the role of Theorem F8's linearity (H_Higgs = X̄) and finds it bears far less weight than “the whole substrate side hangs on exact linearity”. (1) Uniformity alone (the S_D-symmetry that is P1) forces H_Higgs = g(X̄) for some g, because any symmetric function of binary bits depends only on the count; this is exact and is all that the all-orders decoupling (Computations 118, 121) needs. (2) Linearity enters only via additivity: a single-bit function is necessarily affine, so an additive Hamiltonian is exactly affine in X̄ with no O(1/D) remainder; the only route to a nonlinear g is a cross-bit term, an interaction between P1-independent distinctions, additivity is P1 independence read at the level of the energy. (3) The leading value e⁻¹ is robust to relaxing linearity: by the LLN X̄ → 1/2, so Z_H(β) → e^(-β g(1/2)), depending on g only through g(1/2); with the binary gap β = 2, e⁻¹ follows from the single normalisation g(1/2) = 1/2. Verified: four g's (linear + nonlinear deformations) all with g(1/2) = 1/2 give e⁻¹ to leading order, the difference O(1/D); a slope sweep confirms value = e^(-2g(1/2)) exactly in the limit. Net: the substrate-side closure of (B) does not hinge on H_Higgs being exactly linear, its structure is exact from symmetry, its value rests on a vacuum normalisation plus the binary-gap temperature, and linearity refines only the O(1/D).

Computation 123 — Bridge premise (B): two-loop SM-side error budget and the finite-D rigidity of the e⁻¹ matching [EMPIRICAL]
.py.txtcomputation_123.py

Strengthens the SM-side matching (the single open content of (B)) from a “0.8% agreement” to a quantified, rigid, falsifiable test. (A) SM side: running the SM quartic up to M_* = 1285 GeV at two loops gives λ_SM(M_*) = 0.0927 ± 0.0007, the uncertainty dominated by the top and Higgs masses (parametric ±0.00068: m_t ±0.00054, m_h ±0.00040, α_s ±0.00010; truncation ±0.0002), not by logarithmic running. The value converges monotonically with loop order (0.09145 → 0.09234 → 0.09270, steps shrinking), so the match is a genuine two-loop statement, not a low-order accident; the integration is validated by the vacuum instability scale. Against the substrate prediction e⁻¹/4 = 0.09197 the gap is 1.0σ, consistent. (B) Substrate side: the prediction is rigid. By the field-normalisation factor is e⁻¹ exactly and D-independently (the matched-scaling eigenvalue collapse ±√D/Λ → ±1); in the partition-function route the O(1/D) correction would need an implausible D ≈ 63 to close the gap, against a substrate cardinality |D| = V/d₀⁴ (d₀ ∼ nm) larger by dozens of orders of magnitude, so there is no finite-D freedom to absorb the residual. Net: λ_SM(M_*) = e⁻¹/4 is a rigid, parameter-free prediction, consistent with the SM at 1.0σ, with the residual living entirely in the m_t/m_h inputs and sharply falsifiable as those (and the SM matching order) tighten.

Computation 124 — Displacement covariance is positive-semidefinite; the Lorentzian signature is the τ_null sign cut, not isotropy [PROOF]
.py.txtcomputation_124.py

Backs the clause at . The empirical displacement covariance C^(μν) = |D|⁻¹Σ_a (ρ(a)^μ - x^μ)(ρ(a)^ν - x^ν) is a sum of outer products, hence positive-semidefinite (all four eigenvalues → d₀²/4 > 0), whereas the Lorentzian g = diag(-1,1,1,1) is indefinite, so C cannot equal (d₀²/4)g. Under O(4) isotropy C → (d₀²/4)δ (eigenvalue ratio → 1, off-diagonal → 0 as N → ∞): isotropy fixes only the positive Euclidean magnitude and its scale. The signature is the τ_null sign cut of (Computation 2): across a quantile sweep the timelike fraction tracks the cut exactly while C stays PSD throughout. Confirms is read in the post-signature metric.

Computation 125 — Locality of the emergence map ρ_n(a): fluctuation local, bias global (the A6 split) [SURROGATE]
.py.txtcomputation_125.py

Backs the locality verdict and the analytic content of (R) on a deterministic-functional-of-i.i.d.-bits toy (ring sites; pair-local tension T({a,b}) = chord² + ε β_aβ_b; metric ; positions by global MDS and by a locally-supported quasi-interpolation Π^†). (1) Tension and metric are exactly edge-local (zero non-incident edges change under a single bit flip). (2) τ_null is single-flip-stable (|δτ_null| 2ε, ∼ machine zero): an ensemble cut, not a per-draw order statistic. (3) The paper's local-sum A_n carries the O(1/|D_n|) McDiarmid bounded difference (empirical slope -0.73, no position-independence assumed). (4) Under one bit flip the Π^† fluctuation stays within a bounded neighbourhood (93% within ring-distance 8, zero past 12) while the global MDS embedding spreads it (19%, flat in distance). Confirms fluctuation local / bias global: edge-local metric does not imply local fine-placement, the embedding being globally rigid.

Computation 126 — Forced/free decomposition of the tension functional T [NUMERICAL]
.py.txtcomputation_126.py

Backs the structural sharpening after , by direct enumeration of P(D). (1) The admissible magnitudes ( zero-locus + monotonicity) form a cone of dimension 2^(|D|) - (|D|+1), whose ratio to 2^(|D|) rises (0.50 → 0.96 for |D| = 3 → 8): the form is not a finite-parameter family, T_ex one point of it. (2) Several distinct magnitudes (uniform/weighted edge counts, a rank potential) are simultaneously admissible. (3) A monotone reparametrisation preserves every threshold super-level set exactly (downstream reads only the ordering), while a different magnitude changes it (the contingent per-configuration content). (4) The forced content (zero-locus; T(D) ≠ 0, so τ̂ is defined) is invariant across all admissible magnitudes. (5) The per-configuration count is 2^(|D|) = 64 at |D| = 6, of which the ∼ 20 observed flavour parameters consume ∼ 31%.

Computation 127 — P2's irreducible primitive is non-balance of v, independent of eq:pos [NUMERICAL]
.py.txtcomputation_127.py

Backs the reducibility verdict on P2 (open-research item 7.2): almost all of P2 reduces to P1 + , the surviving primitive being one non-degeneracy condition, the non-balance of the elementary directions |Σ_a∈ D v(a) | > 0. (1) A balanced assignment (Σ v = 0, three antipodal V₇-pairs) has maximal edge count yet |Σ v| ≈ 10⁻¹⁶, so secures the count but not the vector bias T(D): non-balance is logically independent of , vindicating as a separately-adopted constraint. (2) G₂ ⊂ SO(7) acts linearly, so g(Σ v) = Σ g(v(a)): over 2000 random SO(7) maps the balanced set stays balanced (drift ∼ 10⁻¹⁶) and norms are preserved, refuting any reduction of non-balance to the Hurwitz/G₂ algebra. (3) Among 2×10⁵ random unit assignments none lands within 10⁻³ of the balanced locus (codimension seven, measure zero): non-balance is generic but, since PST fixes no measure over the structural datum v D→ V₇, not forced.

Computation 128 — The minimal KO split 10=4+6 is forced given the inputs; the only undeived substrate choice is the odd-|D| parity [PROOF]
.py.txtcomputation_128.py

Backs the refined open-research item 7.1 verdict on the minimal-KO-split principle. With the Lorentzian target ≡ 2 8 (the Connes value) and the parity-selected substrate KO-6 (Computation 3), the spacetime KO s satisfies s+6≡ 2 8, so s≡ 4: among s∈{0,…,7} the unique solution is s=4 (next s=12), making spacetime KO-4 the arithmetic complement, not an independent input (Computation 4 verifies the spin structure on an already-four-dimensional M=ℝ× S³). The minimal total ≡ 2 8 with s≥ 1 is 10=4+6 (next 18=12+6), so the split is forced given the substrate six; minimality is non-vacuous only at the lift s=4 over 12, fixed by the Mosco limit. Varying the internal f flips the split, so the odd-|D| parity that fixes f=6 is the single load-bearing undeived choice, not the mod-8 arithmetic.

Computation 129 — A spectral gap reaches the relaxation exponent, not the stationary-mean coefficient: the A6 bias residual is immune to mixing rates [SURROGATE]
.py.txtcomputation_129.py

Backs the §1 A6 spectral-gap verdict. The (R) residual binds at the leading coefficient η_K,R < d₀²/4, not an exponent. (1) A one-parameter family of reversible 2-state chains with fixed spectral gap g= 12 has its stationary mean swept across a 0.8-wide range with the gap identical to machine precision: a relaxation rate does not determine the stationary mean. (2) For A_n(u) with a controlled equidistribution defect ε, the systematic mean deviation is flat in |D_n| (matching the predicted ε G) while the statistical fluctuation decays as |D_n|^(-1/2) (std halving per 4× sample). So a covariance-decay / spectral-gap estimate reaches only the fluctuation exponent; the systematic-mean coefficient is decoupled (the stationary distribution is the kernel eigenvector), and the A6 bias half is foreclosed to mixing / spectral-gap methods.

Computation 130 — The A6 bias residual c_u is a cell-volume equidistribution defect; the matched-scaling count does not force it [SURROGATE]
.py.txtcomputation_130.py

Backs the §1 A6 deterministic-frontier reduction. c_u = |D_n|⁻¹Σ_a |∇ u(ρ_n(a))|² / (|K|⁻¹∫_K|∇ u|²) is the uniform-weight nodal quadrature defect of the emergence cloud, converging to 1 iff the points cell-volume-equidistribute (not the Delaunay finite-element interpolation error). (1) A deterministic volume-equidistributed cloud gives c_u → 1. (2) A matched-count cloud with a central void in a low-|∇ u|² region retains a persistent |c_u-1| ≈ 0.094, flat across a 64× count increase and matching the analytic void bias to 10⁻⁴: refinement (the A1 count) does not close it. (3) The defect grows monotonically with the void radius. So the A6 bias half is conditional on cell-volume equidistribution of the emergence cloud, the deterministic sibling of the open A2 weak-discrepancy, not on the point count.

Computation 131 — The weak A6 discrepancy is blind to the null-cone vanishing-separation set (codim ≥ 1); only the codim-0 node-density binds [SURROGATE]
.py.txtcomputation_131.py

Backs the A2-core characterization. c_u is a single-node average, so the Lorentzian vanishing-separation locus (a codim-≥ 1 pair relation) cannot move it. On a 3-cube with s² = -Δ x₀² + Δ x₁² + Δ x₂²: (1) injecting exact null pairs drives the strong Lorentzian minimum separation to 0; (2) the weak defect |c_u-1| is unchanged by that injection (the paired nodes sit at bulk-typical positions), staying ≈ 0; (3) only a codim-0 central void biases c_u, growing with the void volume. So strong separation fails on the null cone while the weak H¹-tested discrepancy survives, and the A6 residual is purely the codim-0 cell-volume equidistribution of the emergence cloud.

Computation 132 — The tension functional is blind to bulk node density: a void costs nothing in T-terms [SURROGATE]
.py.txtcomputation_132.py

Backs the final foreclosure of the A6 internal-route chain. T reaches physics through three channels, all functions of the distinction data (v D→ V₇ and δ), not the placement ρ: the pairwise tensions, the magnitude ordering of |T(C)|, and the global direction τ̂. (1) Holding the distinction data fixed, all three channels are identical to machine precision across a uniform and a void placement: T takes no spatial-placement argument. (2) The single-node discrepancy c_u reads ρ, so it sees the void, growing with the void radius. (3) The placement-independent toy T-energy changes by exactly 0 across void radii while c_u rises. So a positive-volume void costs nothing in T-terms, and the codim-0 cell-volume equidistribution at the core of A6 is forced by no internal lever.

Computation 133 — A metric-preserving void biases only the flat node-quadrature: a coordinate artifact, not a density axiom [SURROGATE]
.py.txtcomputation_133.py

Backs the coordinate-artifact verdict. The tension data fixes the intrinsic metric only up to diffeomorphism; bulk node density is a distinct diffeomorphism-class invariant it does not pin. Two coordinate charts of the same intrinsic 1D manifold with N geodesically-uniform nodes: a uniform chart, and a void chart x' = CDF⁻¹(x) with a central density dip and metric g₂ = n². (1) Every pairwise geodesic distance is identical to machine precision (∼ 10⁻¹⁶): the void chart is a genuine diffeomorphic representative (same intrinsic geometry, T-realizable). (2) The coordinate uniform-weight node-quadrature c_u, which the bias half requires to approach 1, is unbiased in the uniform chart but biased by ∼ 0.10 in the void chart. (3) The bias grows with the void depth. So the flat c_u is chart-dependent; restoring the void chart's own metric factors makes the defect void-blind (Computation 134), and a void with its co-transforming metric is one gauge orbit carrying no independent density freedom.

Computation 134 — Computation 133's void bias is a dropped-metric-factor artifact; the genuine A6 residual is the fixed-g rate [SURROGATE]
.py.txtcomputation_134.py

Corrects Computation 133. On its void chart (x' = CDF⁻¹(x), metric g₂ = n²): (1) the flat node-quadrature c_u = mean |u'|² / ∫|u'|² is biased ≈ 0.10; (2) the metric-covariant c_u (intrinsic gradient |u'|²/g₂, volume √(g₂) dx' = n dx') equals 1 to machine precision –- void-blind, so the flat bias is purely the dropped g⁻¹ and √g factors; (3) a genuine fixed-metric density modulation (nodes of density n on g = 1, not a diffeomorphism) does bias the covariant c_u. So the metric-preserving void is pure gauge (no bulk-density axiom, no new postulate), and the genuine residual is the diffeomorphism-invariant leading-coefficient inequality η_K,R < d₀²/4, open for the emergence-determined ρ_n and verified only for the i.i.d. surrogate.

Computation 135 — Born rule via Gleason on L²(𝒫(D),μ) [NUMERICAL]
.py.txtcomputation_135.py

Verifies Gleason's dimension hypothesis (dim H = 2^|D| ≥ 3) and demonstrates the Born form P(C) = |ψ(C)|² μ(C) on the substrate Hilbert space for |D| = 2, 3, 4: total probability one, frame-function additivity in an arbitrary orthonormal basis, and non-contextuality of shared projectors. Confirms the |ψ|² rule is the Gleason frame function of the pure state, not an added postulate.

Computation 136 — Gauge anomaly cancellation and the N_c = 3 condition [PROOF]
.py.txtcomputation_136.py

Evaluates all six triangle-anomaly coefficients for one Standard-Model generation as exact functions of the colour rank N_c. All vanish at N_c = 3; the [U(1)_Y]³ coefficient is −2N_c + 6 and the [SU(2)]² U(1)_Y coefficient is (N_c − 3)/6, so anomaly freedom (given the SM hypercharges) forces N_c = 3, the dimension of the M₃(ℂ) colour factor. A consistency check, not a derivation of the hypercharges.

Computation 137 — Equation of state w = −1 from steady-state non-dilution [PROOF]
.py.txtcomputation_137.py

Two exact routes to w = −1: (A) the FRW continuity equation ρ̇ + 3H(ρ + p) = 0 with the non-diluting substrate density ρ̇ = 0 forces ρ + p = 0; (B) a metric-proportional stress tensor T_μν = ρ_Λ g_μν gives T^μ_ν = ρ_Λ δ^μ_ν, matching a perfect fluid at p = −ρ. The load-bearing input is non-dilution; w = −1 follows from it.

Computation 138 — CKM consistency with the S₃-leading structural-scope theorem [EMPIRICAL]
.py.txtcomputation_138.py

Checks that the measured CKM matrix (PDG) is consistent with identity plus a small T(C)-contingent perturbation, as the S₃-leading theorem requires: ‖V − I‖₂ ≈ 0.25 is Cabibbo-sized (O(λ), not O(1)), every off-diagonal is ≤ |V_us|, and the magnitudes are unitary to sub-percent. A falsifiable consistency bound (the values themselves are out of scope).

Computation 139 — The d⁻⁶ Casimir exponent derived from the projection kernel [NUMERICAL]
.py.txtcomputation_139.py

Derives the d⁻⁶ pressure-correction exponent directly from the mode form factor rather than by dimensional analysis: for any smooth even kernel g(k d₀) (g(0)=1, g'(0)=0 by parity), the leading correction is δP = −C a₂ d₀² Γ(6)ζ(6)/(2d)⁶ ∼ d₀² d⁻⁶. Numerically confirmed for four smooth even kernels (fitted exponent −6.00 to <1%; analytic coefficient matched to <0.2%); the exponent is kernel-independent (a non-smooth cusp kernel gives −5 instead), while the coefficient a₂ is kernel-specific (varies 4×), so the d⁻⁶ signature is generic to any smooth even single-scale cutoff and only ξ = 90/π² is PST-specific.

Computation 140 — Robustness of the emergence chain to the form of T(C) [PROOF]
.py.txtcomputation_140.py

Formalises and verifies that the downstream structural derivations read T(C) only through the descending order of |T(C)| (the instantiation sequence) and the global bias τ̂ = T(D)/|T(D)|. For any strictly increasing reparametrisation φ (φ(0)=0) of the magnitudes with τ̂ fixed, at matched threshold τ' = φ(τ) the instantiation sequence, the per-configuration phase sign(|T(C)|−τ), and τ̂ are all invariant (300 random trials); a non-order-preserving reassignment changes the sequence. So the residual freedom in T(C) beyond the ordinal content and τ̂ is inert. Does not derive a canonical T(C) (open, Computation 126); it fixes the exact robustness criterion.

Console output

Idle
Output will appear here. Click any “Run” button to start.

Source files at /pst/computations/*.py. Pyodide v0.26.4 loaded from a public CDN; no data is transmitted from this page to any server.

v29.18