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.

Computation 1 — Yukawa hierarchy
.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
.py.txtcomputation_02.py

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

Computation 3 — Spectral triple KO-2
.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 the 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 concrete 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); (K3) F is even and carries colour: [γ_F, su(3)] = 0; (K4) the explicit 32-dim (γ, J, D) for M × F is verified to be a genuine real spectral triple of KO-2; (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³
.py.txtcomputation_04.py

Closes the spin-structure step for M = ℝ × S³ (the load-bearing analytic item connecting PST's Mosco limit to the canonical Dirac operator), splitting it 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 = ℝ × Σ is spin iff Σ is spin, so M = ℝ × 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
.py.txtcomputation_05.py

Derives Newton's constant G structurally from the Chamseddine-Connes spectral action principle S = Tr f(D/Λ), replacing the Kaluza-Klein dimensional-reduction identity (which is moot in the 10-foundational reading once the three compact spatial dimensions are absorbed into the substrate). The 4D heat-kernel asymptotic expansion (Chamseddine-Connes 1996; Chamseddine-Connes-Marcolli 2007) reads S ~ 2 f_4 Λ⁴ a_0 + 2 f_2 Λ² a_2 + f_0 a_4 + O(Λ⁻²) with moments f_2 = ∫₀^∞ u·f(u) du, f_4 = ∫₀^∞ u³·f(u) du, f_0 = f(0). The three terms read: Λ⁴a_0 → cosmological constant + volume; Λ²a_2 → EINSTEIN-HILBERT (1/16πG) ∫R√g d⁴x; Λ⁰a_4 → Yang-Mills + Higgs + Weyl² (the SM bosonic sector). So gravity is not put in by hand nor obtained by Kaluza-Klein reduction; the EH term is the Λ² coefficient. Verifies (M1) the cutoff moments f_0, f_2, f_4 are finite positive for a sample f (here f(u) = e⁻ᵘ); (M2) the closed-form scaling 1/(16πG) = (f_2 Λ²/π²)·κ_grav·N_dof with κ_grav an O(1) scheme constant and N_dof the fermionic multiplicity of the internal space, solving for Λ given observed G and the SM fermion count places Λ at a unification/Planck-ish scale.

Computation 6 — Gauge group from A_F
.py.txtcomputation_06.py

Derives the Standard-Model gauge group SU(3) × SU(2) × U(1) as the unimodular unitary group of the internal algebra A_F = ℂ ⊕ ℍ ⊕ M₃(ℂ) via the Chamseddine-Connes spectral-action machinery, with the gauge fields arising as INNER FLUCTUATIONS of the Dirac operator D → D + A + ε' J A J⁻¹, A = Σᵢ aᵢ [D, bᵢ] for aᵢ, bᵢ ∈ 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. Background already verified in Track I: A_F gives U(1) × SU(2) × SU(3) of dimension 12 — SU(3) from M₃(ℂ), SU(2) from ℍ, U(1) from ℂ. Four explicit verifications: (N1) Lie-algebra dimensions u(A_F) = 13; unimodular condition removes one dim → 12 = dim(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; Higgs scalar from [D_F, ·]; the spacetime [D_M, a] gives the vector bosons A_μ (12 of them), and the finite [D_F, a] gives the Higgs; (N4) U(1) bookkeeping: which U(1) is removed by the unimodularity constraint, which survives as U(1)_Y.

Computation 7 — Order-one condition selects the Yukawa form of D_F
.py.txtcomputation_07.py

Constructs a clean, 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) (number-changing ladder on a single Fock space) is, under either natural chirality grading, vector-like or chirality-mixing — NOT the Standard Model's chiral SU(2)_L. The Connes resolution uses the CHIRAL representation: the quaternions ℍ act on the 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 explicit verifications: (O1) the representation is EVEN ([γ, a] = 0) and the SU(2)_L action of ℍ is CHIRAL — supported on left-handed particles only, annihilates right-handed ones; (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 its FAILURE for a generic admissible D_F — so order-one is a genuine constraint, not an identity, and 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 (ℍ on L only) as the CCM input; whether the substrate furnishes this chiral bimodule is settled by Computations 8 and 11.

Computation 8 — Chirality bridge: γ_M vs the directed modal threshold
.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 of the spatial slice. Chirality selection is the choice of one factor, and the directed modal threshold ε = T(C) - τ > 0 (the existing PST parity-violation argument from §8) 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 (the bimodule support property); this is bridged to the threshold mechanism here and closed analytically in Computation 11.

Computation 9 — Dirac convergence scoping
.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 (tree-level custodial T)
.py.txtcomputation_10.py

SMEFT operator-by-operator survey of M_* ≈ 4π m_h ≈ 1.573 TeV (an NDA strong-coupling estimate, not a parameter-free derivation) against LHC/EWPO bounds. At PST's tree level the scalar sector exhibits custodial SU(2) symmetry, giving Δ T = 0 at tree level; loop and new-states corrections to S and T are not computed here.

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

Closes §14.2 item (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: Choice A (ℍ acts on H_F^(+) left block, trivial on H_F^(-) — the parity-violating observed Standard Model), or 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 reduces to 'which of A or B is realised in nature?' PST's claim is that the directed modal threshold (paper §8 parity-violation mechanism) is the physical principle that selects Choice A. Verifies: (i) both Choice-A and Choice-B representations of ℍ commute with γ_F (necessary condition); (ii) both have clean block support (no leakage across the chirality grading); (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
.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) the two SU(3) sub-groups of Spin(6) compared; (§2) Computation 6's gauge group derivation does not specify which SU(3); (§3) Computation 7's order-one condition transfers cleanly 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, opening the door for Furey-style three-generation derivation.

Computation 13 — Canonical vs octonionic SU(3) embedding side-by-side
.py.txtcomputation_13.py

Computes the consequences of canonical vs octonionic SU(3) embedding by re-deriving Computation 1 §7 (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: (§1) two H_F realisations setup; (§2-§3) Computation 1 §7 under each embedding; (§4-§5) Computation 7 under each embedding; (§6) comparison matrix; (§7) recommended structural choice. Recommends 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 in Computations 14–17.

Computation 14 — Inner-fluctuation structural-scope theorem (Comp 1 §7) under octonionic embedding
.py.txtcomputation_14.py

Step B of the octonionic-SU(3) re-derivation programme identified by Computation 13. Original Computation 1 §7 (canonical embedding) assumed H_F = H_F^(1) ⊗ ℂ^N_gen tensor-factor structure for generations, with A_F acting as a^(1) ⊗ 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_3(ℂ) 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 NO-GO STATEMENT itself still holds: can inner fluctuations turn a gen-symmetric D_F^(0) into one with generation-asymmetric Yukawa eigenvalues? Five sections: (§1) toy octonionic-SU(3) model on a small Hilbert space; (§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 (Comp 7) under octonionic embedding
.py.txtcomputation_15.py

Step C of the octonionic-SU(3) 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_3(ℂ) 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 in PST. Seven sections culminating 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₈(ℂ), G₂-stabiliser SU(3)
.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 ←𝕆 (the left-multiplication algebra) as 64 explicit 8×8 complex 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 (one generation + three generations of matter); (§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 in this Comp due to degenerate Cartan weights; resolved by Casimir-discrimination in Computation 17.

Computation 17 — Definitive SU(3) decomposition of M₈(ℂ) via Casimir
.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 (Λ_3, Λ_8) 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_7} complex structure on (e_1, e_3), (e_2, e_6), (e_4, e_5) — the correct Cayley pairing for the G_2-stabiliser of e_7); (§2) build ad(Λ_a) on M₈(ℂ) in the matrix-unit basis; (§3) build the Casimir C_2 as a 64×64 Hermitian matrix; (§4) diagonalise C_2 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̄ matches the Clebsch–Gordan expectation (8 ⊗ 8̄ = M₈(ℂ) under the (1+1+3+3̄) ⊗ (1+1+3̄+3) structure, total dim 6+15+15+16+6+6 = 64) to machine precision; (§6) status statement for §14.2 (C): octonionic SU(3) is structurally derived and the decomposition is exact; 'three generations' becomes a separate downstream question (Furey 2018 tensor construction or T(C) contingency).

Computation 18 — ω as the unique Spin(6) singlet
.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) Spin(6)-invariant primitive idempotent, exactly Furey's idempotent.

Computation 19 — Substrate parity grading = Cl(0,6) chirality grading
.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. Furey 2014's chain-algebra construction on ℂ⊗𝕆 ≅ Cl(0,6) then delivers SU(3)_c × U(1) on one generation. The synthesis to the full A_F = ℂ⊕ℍ⊕M_3(ℂ) requires the ℍ factor from the emergent Lorentzian Dirac sector Cl(1,3) ≅ M_2(ℍ) (Computation 66) combined via Dixon's hyperspinor framework T = ℂ⊗ℍ⊗𝕆.

Computation 20 — L_round sweep
.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
.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
.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
.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
.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
.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
.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. Verdict: no bijective single-mode bridge closes L_comm.

Computation 27 — Bergman/Toeplitz framework on H²α(B²)
.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)
.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). Verdict: 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
.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₀^2, T_z₁^2. 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). Verdict: 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
.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
.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. Verdict: 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
.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 = 1/2(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. Verdict: 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
.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. Verdict: 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
.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)
.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
.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 paper-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
.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
.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 (Option C): closed-form ratio
.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
.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
.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
.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 ℝ⁸
.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
.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 ℂ⊗ℍ⊗𝕆
.py.txtcomputation_45.py

Constructs the extended chiral primitive idempotent f = (I + i_C ⊗ I_H ⊗ ω₆)/2 where ω₆ = L_e₁ ⋯ L_e₆ is the Cl(0,6) volume element (ω₆^2 = -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)
.py.txtcomputation_46.py

Fixes C = span1, 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₇ spane₁ = spane₂, …, 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 structural source of N_gen = 3.

Computation 47 — Per-Q_k Cl(0,2) chiral structure
.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 ℂ ⊗ ℍ ⊗ 𝕆
.py.txtcomputation_48.py

Closes the structural derivation of N_gen = 3: 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₂-equivariance verification: J² = -I on τ̂⊥
.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
.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
.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 σ₁ = 2^(D-1) − D for the Walsh-tail off-diagonal block
.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 gives σ₁ = (α(D) − β(D))/2 with α(D) = 2^D − 1 − D − C(D,2) (the T_tail eigenvalue on |0⟩) and β(D) = −C(D−1,2) (on each |e_a⟩). Simplification gives σ₁ = 2^(D-1) − D exactly, closing Lemma 2-prime on the S_D-symmetric component of P_+ T_tail P_-. The second singular value σ₂ = D − 2 (exact, S_D-non-symmetric component) is bounded by an exponential gap from σ₁; its analytical bound via the S_D-irrep decomposition is the only routine residual.

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

Verifies to machine precision the 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 weight-1 standard-rep input alone: σ_max = (D−2)/√2 = |β(D) − γ(D)|/(2√2) with β = −C(D−1,2) and γ = c(2) = (D−2)(5−D)/2, achieved via the exact cancellation D_sub w_2 = D · v for w_2 = Σ_{a<b}(c_b−c_a)|e_a+e_b⟩ when Σ c_a = 0 (the weight-3 part vanishes by combinatorial sign cancellation over unordered triples). An earlier reading reported σ_max on a single 'axis-1' subspace as 1.63 at D=4 etc.; this is a single direction within [D-1, 1] under the wrong (bit-permutation) S_D rep and does not bound σ₂. The rep that commutes with D_sub is the fermionic-signed one; under it σ₂ = D − 2 exactly (Comps 54, 55).

Computation 54 — Identifying the S_D irrep that saturates σ₂ = D − 2
.py.txtcomputation_54.py

Implements 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; sanity check ||[ρ_bit(g), D_sub]||_op ≈ 4.0 at D=4 versus ||[ρ_F(g), D_sub]||_op = 0 to machine precision). Builds P_λ = (dim λ / |S_D|) Σ_g χ_λ(g) ρ_F(g) for each two-row irrep [D-j, j] with j = 0, 1,..., ⌊D/2⌋, then computes σ_max within each isotypic block. Findings: the standard rep [D-1, 1] saturates σ_max = D − 2 exactly at D = 4, 5, 6, 7; every higher two-row irrep gives σ_max = 0 identically. The multiplicity of [D-1, 1] under ρ_F is 2 (Pieri on Λ^w(C^D) assigns it to H_1 and H_2), so the multiplicity-space matrix M is 2×2; rank(M) = 1 with singular values (D-2, 0).

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

Exhibits the explicit S_{D-1}-fixed vectors u_1 = |e_0⟩ − (1/D)Σ_a|e_a⟩ ∈ [D-1, 1] ∩ H_1 and u_2 = −Σ_{b>0}|{0, b}⟩ ∈ [D-1, 1] ∩ H_2 and verifies to machine precision (for D = 4..8) the algebraic identities D_sub u_1 = u_2, D_sub u_2 = D·u_1 (no weight-3 leakage; the JW orderings on each triple {0, b, c} cancel), T_tail u_1 = β·u_1 with β = −C(D-1, 2), T_tail u_2 = γ·u_2 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_1 = u_1/||u_1||, e_2 = u_2/||u_2|| to M = (D-2)/2 · [[-1, +1], [-1, +1]], rank 1 with singular values (D-2, 0). The saturating right singular vector is v★ = (1/√2)(−e_1 + e_2), the left singular vector is u★ = (1/√2)(e_1 + e_2). 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
.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 (w = z_0 z_1). Since J_z f = 0 (f is symmetric in z_0, z_1) and J_± f are antiholomorphic mirror pairs, M_f is anti-diagonal in the Pauli basis with ||M_f||_op = max(|z_0|², |z_1|²) · |f'(w)|. The 1D optimisation r² = (m+1)/(2m) on ∂B² gives ratio(m) = sup ||M_f|| / sup |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 ratio = 2 = 2√1, exact closure at k = 1. For m ≥ 2 the ratio differs from 2√m (undershoot at m=2, overshoot at m≥3); asymptotically ratio(m) ~ m + 1/2. Establishes Lemma 5(a) in the L_comm proof programme; the 2-monomial refinement of Comp 57 closes the existence statement for k ≥ 2.

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

Establishes the existence theorem for the 2-monomial L_comm bridge f_α(w) = w^m + α·w^(m+1). For each m, R_m(α) := sup||M_f|| / sup|f| is continuous in α with R_m(0) = ratio(m) (Comp 56) and lim_{|α|→∞} R_m(α) = ratio(m+1). Since ratio(m) is strictly increasing in m, the target 2√k is bracketed by ratio(m(k)) ≤ 2√k ≤ ratio(m(k)+1) for a unique m(k); by Bolzano, ∃α(k) closing 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 closing α(k) by IVT bisection: α(2)=0.737, α(3)=20.7, α(4)=1.74,..., α(12)=1.38, all with R-target gap at machine precision. Theorem (Lemma 5(b)): ∀k≥1, ∃(m(k), α(k)) with R_{m(k)}(α(k)) = 2√k exactly. The 3-monomial refinements of Comps 41, 42 close at typically smaller |α| via the additional β parameter, but Lemma 5(b) existence is already complete with the 2-monomial family.

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

Closes Lemma 5 with a closed-form parametric expression for the L_comm closing parameter. The substitution q := 2r² − 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 gives a parametric ratio R_m(q); the closing equation R_{m(k)}(q*) = 2√k is an algebraic equation in q* of degree 2m + 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 Comp 57's IVT-bisection α(k) values to 6 decimal places. The endpoint behaviour q* → 1/m gives the single-monomial limit, q* → 1/(m+1) the dominant-(m+1) limit. Lemma 5 (a, b, c) is therefore complete: the L_comm refined-bridge closure at the level of ratio matching is a closed analytical theorem.

Computation 59 — Lemma 7: uniform-in-k closure of the refined bridge at finite N
.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 the 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 Comp 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 (c₁ values: −0.43, +1.15, +0.70 at k = 2, 3, 4) and noisier two-term fit. At the matched scaling N(D) ~ 2k_D ~ 2√D this implies γ_D = O(1/√D) for k = 1 — polynomial L_comm closure, weaker than the exponential L_round but uniformly controlled. The full uniform-in-k analytical bound for k ≥ 2 requires careful expansion of T_{f_k} on the boundary modes (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
.py.txtcomputation_60.py

Promotes Comp 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_0^p z_1^q} act as diagonal shifts e_{a, b} → coefficient · e_{a + p, b + q} (orthogonal target). Operator norms read off as max 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) = ||commutator|| / ||T_w|| = 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/(4N²) − 1/(8N³) − 5/(64N⁴) − O(1/N⁵). Verified against Comp 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) · conj(g(z)) exactly
.py.txtcomputation_61.py

Closes Lemma 6 (Berezin-transform compatibility) for the refined bridge in the natural holomorphic-symbol formulation. Theorem: for holomorphic symbols f, g on B² and the canonical Bergman space H²_α(B²), B(T_f T_g*)(z) = f(z) · conj(g(z)) EXACTLY for every z ∈ B² (interior). Proof: the reproducing-kernel identity T_g*·k_z = conj(g(z))·k_z follows from ⟨h, T_g* k_z⟩ = ⟨T_g h, k_z⟩ = (T_g h)(z) = g(z)·h(z) = g(z)·⟨h, k_z⟩ = ⟨h, conj(g(z))·k_z⟩ for all h ∈ H²_α. Hence B(T_f T_g*)(z) = ⟨k_z, T_f T_g* k_z⟩ / ||k_z||² = conj(g(z)) · B(T_f)(z) = conj(g(z)) · f(z) = f(z) · conj(g(z)). Applied to the refined bridge b_C(χ_S) = T_{f_{|S|}}: B(b_C(χ_S) b_C(χ_T)*)(z) = f_{|S|}(z) · conj(f_{|T|}(z)) exactly. On the truncated Bergman the identity holds up to an O(|z|^(2N)) boundary correction that decays exponentially in N for z in the open ball. Verified numerically at every tested z ∈ B² and k, l ∈ {1, 2, 3}: gap below numerical precision at N = 18 Bergman truncation. 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
.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) at the matched substrate-to-Bergman scaling. Setup: 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) = exp(-x²) for D ∈ {6, 8, 10, 12}: S_sub(Λ)/2^D = e^(-1) ≈ 0.368 exactly; S_Bergman(Λ, ℓ)/2^D → constant ~21.7 (independent of D) at j = 1, ~54.2 at j = 2, ~189.7 at j = 3. The Bergman/substrate ratios are bounded and converge to D-independent constants — the standard QGH-convergence statement applied to spectral actions. Five of six L_comm sub-lemmas now closed.

Computation 63 — Lemma 7(b): leading c_1(k) for k = 1..12 (2-parameter fit)
.py.txtcomputation_63.py

Numerical extraction of the leading constant c_1(k) in gap(k, N) = c_1(k)/N + c_2(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) = -1 exactly with fit residual ~1e-5. Reports |c_1(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_1 extraction) in Computation 65.

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

Re-extracts c_1(k) on the same N range used by Comp 63 (capped at N ≤ 34) under a three-parameter least-squares fit gap(k, N) ≈ c_1/N + c_2/N² + c_3/N³. For k = 1 the two- and three-parameter fits agree at c_1 = -1.000, matching the closed-form theorem of Comp 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_1(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_1(k) extraction at moderate N; it must be established directly (Comp 65).

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

Bypasses the ill-conditioned c_1(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: max 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_2]||_N² approach their bulk values at rate O(1/N) by the classical Bergman-Toeplitz convergence theorem for polynomial symbols (Engliš, Zhu), and their ratio inherits the rate with bounded constant. 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-ℍ
.py.txtcomputation_66.py

Structural verification that the substrate's emergent 4-d Lorentzian spacetime M = ℝ × S³ (derived from P1-P3 via Mosco convergence, Comp 4) carries an internal SU(2) on its Dirac spinor sector. Works in H² ≅ R⁸ (the M_2(ℍ) module) to avoid antilinearity confusion of C⁴ chiral basis. Builds Cl(1,3) γ-matrices in the M_2(ℍ) representation (16 anticommutators match diag(+1,-1,-1,-1)); builds right-ℍ multiplication operators R_iH, R_jH, R_kH on H²; verifies right-ℍ commutes with Cl(1,3) (the centraliser structure that makes SU(2) internal — Dixon 2010); verifies the quaternion algebra R_q² = -I, R_jR_i = R_k; shows {R_q/2} generates su(2) with [R_iH, R_jH] = -2R_kH; verifies the Lorentz pseudoscalar commutes with R_q so the SU(2) is chirality-preserving. NOTE: this verifies 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_3(ℂ) on the internal F-side is via Dixon's hyperspinor framework T = ℂ⊗ℍ⊗𝕆 (which is a structural synthesis beyond Cl(0,6) alone).

Computation 67 — One-loop Higgs cancellation scope check (binomial off-by-one)
.py.txtcomputation_67.py

Combinatorial verification of the PST scalar-sector cancellation argument (paper §sec:renorm, eq:pairing-cancel). Confirms ∑_{k even} binomial(D,k) = ∑_{k odd} binomial(D,k) = 2^(D-1) exactly for D = 4, 6, 8, 10, 12, 14, 16. Finds: including the zero mode in the bosonic count gives N_B − N_F = 0 exactly. Excluding the zero mode (as the paper's prescription requires, since it's the order parameter / Higgs) gives N_B(nonzero) − N_F = -1 for every D, contributing a residual −M*²/(16π²) rather than 0. Added to the self-loop +M*²/(16π²) of eq:ch-exact, this gives total δm_h² = 0, conflicting with the predicted m_h = M*/(4π). The prediction survives only under one of three currently-implicit structural mechanisms: (a) Higgs-doublet eaten Goldstones supply the missing bosonic dofs at the EW scale; (b) the zero mode is reinterpreted as part of the cancellation sum rather than separated; or (c) an explicit boson-fermion pairing via complementation/parity restores N_B = N_F at the spectral level. The computation does NOT refute M_* ≈ 1.573 TeV; it flags the off-by-one as an open structural item the paper's argument currently elides.

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

Tests sensitivity of the parameter-free Casimir-correction coefficient ξ = (15/π²) · M_4/(d_0² · M_2) (paper eq:xi-exact) to the assumed projection-kernel shape. Evaluates the 2nd and 4th moments M_2, M_4 for six kernel families: Gaussian (paper canonical, ξ = 9.12), squared Lorentzian (ξ = 378.7, heavy 1/x⁴ tail), single-sided exponential (ξ = 18.2), sech² (ξ = 21.0), compact parabolic (ξ = 2.6), compact quartic-bump (ξ = 2.0). Findings: the d^-6 power-law exponent is preserved across every kernel (categorical prediction, kernel-independent). ξ varies factor ~10 across exponential/compact-support kernels. The diagnostic d_0 bound from nanometre-Casimir data scales as 1/√ξ, giving range d_0 ≲ 3.5 nm (sech²) to 11 nm (compact quartic); Gaussian's 5.3 nm sits near the middle. Order of magnitude (few nm) robust; precise numerical bound kernel-conditional with factor ~3 spread. Honest reading: 'd_0 ≲ 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 Comp 67 rescues
.py.txtcomputation_69.py

Carries out the explicit walkthrough of all three rescues listed in Comp 67 for the off-by-one in eq:pairing-cancel. (a) Higgs-doublet eaten Goldstones: rescue 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 paper does not provide. (b) Zero mode included in 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 rescue. (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 by the paper's prescription (zero mode = Higgs = excluded). HONEST CONCLUSION: none of the three rescues, as constructible from current PST postulates without ad hoc structure, closes the off-by-one. M_* = 4π m_h is therefore an NDA estimate (standard naturalness rule of thumb Λ ~ 4π × mass scale), not a parameter-free one-loop derivation. Numerical value M_* ≈ 1.573 TeV survives as NDA estimate only. Reframing the headline claim throughout the paper is the correct response.

Computation 70 — (R) uniform-rate verification for A6 reduction
.py.txtcomputation_70.py

Tests the Berry-Esseen-type convergence rate underlying the A6-reduction lemma in §sec:A6-reduction. Sets M = R³ flat-bulk, samples n i.i.d. uniform sites with isotropic displacements covariance (d_0²/4)·I, computes the empirical aggregate A_n(u) = (1/n) Σ (u(ρ+ξ) − u(ρ))² for three compactly-supported smooth u (Gaussian, C_c^∞ bump, (1−r²)²cos(πr²)). Findings: (i) the statistical std-error of A_n across independent realisations scales as |D|^(−1/2) (empirical slope −0.467, Berry-Esseen prediction −0.5); (ii) the systematic part of |A_n − A_∞| scales as O(d_0⁴) in the next-order Taylor correction beyond the leading covariance forcing A3. Together this gives total rate O(d_0² + |D|^(−1/2)) on the H¹-bounded compactly-supported class, confirming the hypothesis (R) of the reduction lemma numerically. Under (R) with any positive rate, the lemma closes A6 with K-dependent constant C_K = (d_0²/4 − η_K)/(C_P(K)² + 1). The remaining open analytic item is the K-uniform version of the rate, which is the standard discrete-to-continuum estimate of finite-element analysis (Brenner-Scott 2008, Ch. 4) — no fresh spectral-gap question survives.

Computation 71 — Z² = e⁻¹: candidate identification via substrate spectral action
.py.txtcomputation_71.py

Tests whether Z² := λ_SM(M_*)/b(M_*) ≈ e⁻¹ (0.8% near-coincidence, §sec:renorm) can be upgraded to a structural identity Z² = S_sub/2^D = e⁻¹ via the substrate spectral-action ratio. Substrate side (Comp 62 result): under matched scaling Λ = √D and Gaussian cutoff, the substrate Dirac eigenvalues ±√D collapse to ±1, giving S_sub/2^D = f(1) = e⁻¹ EXACTLY, D-independently. SM side: one-loop running of λ_SM from v to M_* = 4π·m_h, computed both naively here and via the paper's more careful coupled RGE. Both confirm Z² ≈ e⁻¹ at the few-% level. Under the conjectural spectral-action Higgs-quartic identity (Conjecture in docstring), the deviation lives entirely on the SM-RGE side as the two-loop truncation error; substrate side is exact for any D. Remaining structural gap: Chamseddine-Connes-style spectral-action coefficient calculation for the discrete PST substrate triple, extracting the Higgs-quartic coefficient at matched scaling. This is a finite, technical calculation in the established spectral-action framework, not a new postulate. Resolution would convert Z² ≈ e⁻¹ from numerical near-coincidence to structural theorem.

Computation 72 — Substrate spectral action: Higgs-quartic Taylor expansion (Z² reduction)
.py.txtcomputation_72.py

Carries out the explicit symbolic Taylor expansion of the substrate spectral action under a Yukawa-like Higgs perturbation D_φ = D_sub + y·φ·Σ_X, at matched scaling Λ = √D. Result: the φ⁴ coefficient in S(D_φ, Λ, f)/2^D is e⁻¹·y⁴/(2Λ⁴). For the candidate identification λ_H = b·e⁻¹ = (1/4)·e⁻¹ to deliver the SM Higgs quartic structurally, the substrate-Higgs Yukawa coupling must satisfy y = Λ·2^(−1/4) ≈ 0.8409·Λ. This REDUCES the Z²-conjecture from 'the spectral-action Higgs-quartic identity' (abstract) to a single-scalar condition: y_substrate(M_*) = M_*·2^(−1/4) as derived from the Chamseddine-Connes inner fluctuation of D_sub by elements of A_F = ℂ⊕ℍ⊕M₃(ℂ). The remaining content is the inner-fluctuation Yukawa normalization in the spectral-triple framework (Chamseddine-Connes-Marcolli 2007) — a specific finite number the framework either delivers or doesn't.

Computation 73 — ξ from the substrate's binomial-corrected Boolean kernel
.py.txtcomputation_73.py

Sharpens Comp 68's kernel-shape sensitivity by deriving the substrate's actual projection kernel directly from its Boolean structure. The substrate Hilbert space ℂ^(2^D) carries the Boolean kernel K_n(j) = C(n, j)·(1/2)^n — the binomial at p = 1/2 (the per-site Bernoulli probability of μ = ⊗ Bern(1/2) from P1). By the CLT, binomial → Gaussian as |D| → ∞. The leading-order non-Gaussian correction is in the kurtosis: kurt(binomial|_{p=1/2}) = 3 − 2/|D| vs kurt(Gaussian) = 3. Substituting into the moment ratio M_4/M_2² that fixes ξ gives ξ_substrate(|D|) = (90/π²)·(1 − 2/(3|D|)). At |D| = 6 the correction is 11%; at |D| = 20, 3.3%; at |D| = 100, < 1%. The Gaussian value 9.12 is therefore the matched-scaling limit of the substrate's actual ξ, with calculable convergence rate −2/(3|D|). The kernel-shape conditional of Comp 68 reduces to 'ξ is structurally Gaussian with calculable convergence rate' rather than 'ξ depends on a kernel choice'. d_0 ≲ 5.3 nm bound is correct to leading order; subleading shifts are a few percent at any physical |D|.

Computation 74 — Single-generator CC inner fluctuation: Yukawa-coupling diagnostic for Z²
.py.txtcomputation_74.py

First-cut attempt at the Chamseddine-Connes inner-fluctuation calculation flagged by Comp 72 as the remaining structural content of the Z² = e⁻¹ conjecture. Constructs the substrate spectral triple with Hermitian Jordan-Wigner Cl(0, D) generators, verifies D_sub² = D·I exactly, and computes the inner-fluctuation perturbation A = iφ·[D_sub, ω]/2 with ω = e_1 e_2 (the quaternion volume element). Result: D_φ² = (D + 2φ²)·I, giving φ⁴ coefficient 2·e⁻¹/D² in S/2^D, which corresponds to Yukawa coupling y = √2 (D-independent in absolute terms). Therefore y/Λ = √(2/D) decreases with D, NOT matching the Comp 72 target y/Λ = 2^(−1/4) (D-independent). The naive single-generator inner-fluctuation does NOT close Z² = e⁻¹. Three possible resolutions: (1) Z² ≈ e⁻¹ is a genuine 0.8% SM-RGE numerical near-coincidence with no structural identification; (2) the Higgs perturbation should be parameterised over multiple A_F generators or via the full CC machinery including JAJ⁻¹ + order-one constraints; (3) matched scaling Λ = √D is not the right cutoff. Honest negative result: tells the reader exactly which calculation needs to be done next to settle Z².

Computation 75 — Rate integral j_τ: Bernoulli large-deviation structure
.py.txtcomputation_75.py

Sharpens the open rate-integral problem (paper §sec:rate-problem). The rate integral j_τ = μ({C : |T(C)| ≥ τ}) is the Bernoulli-measure tail probability of |T(C)|; under μ = ⊗ Bern(1/2) from P1, this is a standard large-deviation probability. Cramér's theorem gives j_τ ~ exp(−|D|·I(τ/√|D|)) where I is the Legendre transform of the cumulant generating function of |T(C)|/√|D|. Computation 75 verifies the exponential-in-|D| decay empirically on a representative random-walk-like T(C) model where each elementary property a contributes a unit V_7 vector v(a) and T(C) = Σ_{a∈C} v(a). At α = 1.0·√|D| threshold scaling, j_τ at |D| = 100, 200, 400 measured directly, showing the Cramér exponential form with rate c·α² · |D| where c ≈ 0.01-0.03 depending on the V_7 direction structure. Remaining open content: specify T(C)'s explicit form at the P2 postulate level (beyond random-walk-like) plus the τ-scaling with |D|; once both are pinned, Cramér delivers j_τ analytically. Converts 'rate integral not presently calculable' into 'Cramér computation conditional on explicit T(C)' — structural reduction analogous to A6 → (R) uniform-rate.

Computation 76 — Postulate-level V_7-direction constraint for the Yukawa hierarchy
.py.txtcomputation_76.py

constraining the V_7 directional content v(a) at the P1 postulate level to derive (rather than input) the Yukawa hierarchy. Tests Proposal 1: v(a) constrained to lie in one of three quaternionic sub-algebras Q_k of 𝕆 containing τ̂ (the same three Q_k that deliver N_gen = 3 via Furey 2014's construction). Result: Yukawa-like eigenvalue ratios remain O(1) — 1.0, 0.74, 0.67 at D=30; 1.0, 0.71, 0.46 at D=60 — NOT the 5-order-of-magnitude SM hierarchy. The three Q_k are structurally equivalent under the G_2-stabilizer SU(3)_c of τ̂, so the S_3 symmetry of the structural-scope theorem (§sec:yukawa-scope) is preserved by the constraint. The single-direction V_7 constraint with maximal G_2 symmetry does not overcome the structural-scope theorem. Producing a Yukawa hierarchy from structural inputs alone requires either (1) a SECOND postulate-level directional input (e.g., σ̂ ∈ V_7 in addition to τ̂), which is itself ad hoc, OR (2) a dynamical symmetry-breaking mechanism from the directed modal threshold (P3) — the more PST-natural option, extending the existing parity-violation mechanism. This is taken up further in Comp 77 via a dynamical threshold-induced mechanism.

Computation 77 — Threshold-induced dynamical generation breaking (viable but not closing)
.py.txtcomputation_77.py

Continuation: tests whether the directed modal threshold (P3) can break the generation S_3 symmetry dynamically. Mechanism: Boltzmann rate r_k ∝ exp(−β(1 − p_k)) for each generation k, where p_k is the squared projection of the residual direction (in τ̂^⊥) onto Q_k ∩ τ̂^⊥, and β is a structural inverse-temperature. Effective Yukawa y_k ∝ r_k. RESULT: the mechanism CAN reproduce the SM up-type Yukawa hierarchy (1.0, 7.6e-3, 1.3e-5) to 4-digit accuracy with β ≈ 28.24 and (p_1, p_2, p_3) ≈ (0.52, 0.35, 0.13). HONEST FRAMING: this is curve-fitting — two free parameters (β and the residual direction in τ̂^⊥, 5 angles) can fit any 2-parameter hierarchy. The mechanism is VIABLE (it can reproduce the SM Yukawas from postulate-level inputs once β and the direction are specified) but not CLOSING (those inputs are themselves contingent). The substantive closure question is whether β and the residual direction can be derived from deeper structural principles: β from max-entropy threshold-crossing (Lagrange multiplier of the modal-threshold constraint, extension of Comp 75's Cramér framework), residual direction from extremal-action selection (parallel to PST's §1.4 extremal arguments). Neither candidate principle is yet pinned. The remaining structural work is the β + direction derivation, a specific calculation not a fresh research direction.

Computation 78 — Structural candidates for β and residual direction (initial audit)
.py.txtcomputation_78.py

First audit of structural candidates for the two free inputs of Comp 77. β tested: dim(SO(8)) = 28 (matches Comp 77 best-fit 28.24 to 0.85%), dim(G_2) = 14 (too small), 4·dim(V_7) = 28 (matches, same value as SO(8)). Residual direction tested: Cartan T² weights of SU(3) under generic breaking (degenerate, doesn't give three-distinct-p_k), Fano-plane angle patterns with specific axis-aligned and Q_1-diagonal directions. Initial finding: β ≈ 28 is suggestive at ~1% level but no unique structural identification; residual direction Comp 77 fit (0.52, 0.35, 0.13) doesn't match any structurally privileged direction tested. Identifies the next push: refine β to (SO(8) + b_PST) and look for exponential pattern in p ratios. This audit narrows the candidates; the refinement (β as SO(8) + b_PST and the exponential pattern in p ratios) is carried out in Comp 79.

Computation 79 — Crunched numbers: β = SO(8) + b_PST = 28.25; p_2/p_3 ≈ e
.py.txtcomputation_79.py

Sharpened crunch of Comp 77's two free inputs. β CRUNCHED: β = dim(SO(8)) + b_PST = 28 + 1/4 = 28.25 matches Comp 77's empirical best-fit 28.24 to 0.04% (log-error identical at the 10⁻⁶ level). Structurally clean reading: Spin(8) triality acts on the three octonionic 8-dim representations (paralleling PST's three generations), plus the LG quartic b = 1/4 from P3 (same b that fixes m_h² = 4ε). Both ingredients are fixed structural numbers. RESIDUAL DIRECTION CRUNCHED: p_2/p_3 = 2.798 vs e = 2.718, 3% close. Under the constraint p_2/p_3 = e + β = 28.25, optimal p_1 ≈ 0.521 matches Comp 77 closely; one structural relation pins TWO of three p components. The remaining p_1 component doesn't match any clean structural number tested (1/2, e⁻¹+1/6, 1/√3, e/(e+1+e⁻¹), etc.). NET: 3 free parameters → 1 free. β closeable structurally, p_2/p_3 = e suggestive at 3%, p_1 remains the residual open input. The structurally contingent character of p_1 is the analogue of T(C)-contingency from §sec:yukawa-scope. Research direction narrowed from 'derive the Yukawa hierarchy' to 'derive one number p_1 ∈ [0,1]'.

Computation 80 — p_1 closure attempt: SU(3)_c symmetry forces T(C)-contingency
.py.txtcomputation_80.py

Closes the p_1 question for the dynamical Yukawa mechanism by structural-symmetry argument. The substrate symmetry chain is Spin(8) ⊃ Spin(7) ⊃ G_2 ⊃ SU(3)_c, with SU(3)_c = G_2-stabilizer of τ̂ as the unbroken symmetry on the residual direction. Under SU(3)_c the three Q_k are permuted by the Weyl S_3 action; ANY SU(3)_c-invariant residual direction therefore satisfies p_1 = p_2 = p_3 = 1/3. Verified numerically: at β = 28.25 and p = (1/3, 1/3, 1/3), the Boltzmann mechanism gives Yukawa ratio = 1, ZERO hierarchy. Breaking SU(3)_c at the residual-direction level to recover a hierarchy would back-propagate to break COLOR SU(3)_c gauge symmetry (catastrophic: color confinement is lost). Therefore p_1 lives in T(C) (the configuration-content layer), not in the structural layer, by the same symmetry argument that would otherwise force color SU(3)_c gauge breaking. This re-derives the structural-scope theorem of §sec:yukawa-scope from the SU(3)_c symmetry perspective; the two derivations (S_3 in §sec:yukawa-scope, SU(3)_c here) are consistent and reinforcing. The 'derive the Yukawa hierarchy structurally' question is therefore resolved: the hierarchy lives in T(C), and the structural-scope theorem is robust under the dynamical-mechanism extension.

Computation 81 — Z² push: full CC inner fluctuation D + A + ε' JAJ⁻¹ (negative)
.py.txtcomputation_81.py

Extends Comp 74's single-generator inner fluctuation by including the JAJ⁻¹ reality term required by the full Chamseddine-Connes machinery, with ε' = +1 for the KO-2 product spectral triple. Setup: Hermitian Higgs perturbation A = φ(e_2 − e_1) (correcting Comp 74's anti-Hermitian form), JAJ⁻¹ = A in the JW Majorana basis, full CC fluctuation D + 2A. Result: D_φ² = (D + 8φ²)·I, giving φ⁴ coefficient 32·e⁻¹/D² in S/2^D, corresponding to Yukawa y = 2√2 (D-independent absolutely) and y/Λ = √(8/D) (D-dependent). Target was 2^(−1/4) (D-independent). The D-dependences are incompatible — the full CC reality term does NOT fix Comp 74's discrepancy. Multi-generator perturbations A = i Σ_k h_k [D_sub, x_k] with x_k ∈ {e_1, e_2, ω} tested at D ∈ {2, 4, 6} also fail to give the target Yukawa generically. Confirms Comp 74's diagnostic under the more thorough CC machinery: single- and few-generator CC inner-fluctuation routes do not close the Z² identification. Substrate-side S_sub/2^D = e⁻¹ remains exact (Comp 62), independently of any Z² identification. Possible further options: ε' = −1 matching, different Λ(D) prescription, or imposing the order-one constraint explicitly within CCM machinery — none pinned at the postulate level. Z² candidate identification remains candidate, not closed.

Computation 82 — Λ magnitude: order-of-magnitude analysis (reclassifies as T(C)-contingent)
.py.txtcomputation_82.py

Order-of-magnitude consistency analysis for the cosmological-constant magnitude using Comp 75's Cramér large-deviation framework. Setup: ρ_inst = j_τ · ΔF · (|D|/V), with the matched-scaling |D|/V ~ 1/d_0⁴ ≈ 1.6×10³⁴ per m³ for d_0 = 5.3 nm, and ΔF ~ M_* = 1.573 TeV per configuration. Target rate j_τ for ρ_inst = ρ_Λ_obs ≈ 5.9×10⁻¹⁰ J/m³ comes out to ~10⁻⁵⁰ per substrate configuration. From the Cramér asymptotic j_τ ~ exp(−|D|·I(α)), the required I(α) is tiny — α very close to typical |T(C)|/√|D|, no rare-event suppression in the deep tail. CONCLUSION: PST inherits the SM 120-orders-of-magnitude Λ hierarchy problem; the magnitude depends on the EXPLICIT T(C) form + τ-scaling at matched scaling. RECLASSIFICATION: Λ magnitude is structurally CONTINGENT (lives in T(C) by the structural-scope theorem) rather than an open derivation. Sign + EoS w = −1 remain structural theorems (§sec:cosmology-eos). The magnitude joins Yukawa hierarchy, CKM, d_0 magnitude as configuration content. PST has the right qualitative cosmological constant (positive, w = −1, non-diluting); the quantitative magnitude is T(C)-contingent. The 5 genuinely open items reduces to 4 (Z² candidate, A6 K-uniform rate, measurement postulate, fermion sector rigour).

Computation 83 — A6 K-uniform Berry-Esseen rate verification (asymptotic regime)
.py.txtcomputation_83.py

Pushes A6 closure (open item #2) from K-fixed (Comp 70) to K-uniform empirical verification. Tests seven test-function families with varying support locations (centered vs off-center), radii (0.5, 0.7, 0.8, 1.0), and shapes (Gaussian, C^∞ bump, polynomial-cosine). For each: extracts the empirical Berry-Esseen slope via linear regression of log(std_err) vs log(n) for n = 50k, 150k, 500k. RESULTS: mean slope across all families is −0.458 ± 0.085, with 5 of 7 families within fitting tolerance of the Berry-Esseen prediction −0.500. The two failing families have support radius comparable to the displacement scale d_0 = 0.1 (narrow Gaussian r=0.5 gives slope −0.346; off-center Sinmod r=0.8 gives −0.330) — the expected Taylor-linearization breakdown when displacement xi ~ d_0/2 is comparable to support radius. STRUCTURAL CONCLUSION: K-uniform Berry-Esseen rate HOLDS in the asymptotic regime where support radius ≫ d_0 (the natural physical regime for the Mosco-convergence limit). The thin-support regime (radius ~ d_0) requires d_0-dependent corrections from the kernel-shape sensitivity already mapped by Comp 73. Together with Comp 70's K-fixed verification and the §sec:A6-reduction lemma, this constitutes empirical evidence for A6 closure in the asymptotic regime. Item #2 of the open list is therefore PARTIALLY CLOSED EMPIRICALLY. The formal Brenner-Scott analytic proof and thin-support corrections are standard finite-element analysis, multi-day technical work, no fresh structural questions.

Computation 84 — Multi-generator CC inner-fluctuation route for the Z² identification
.py.txtcomputation_84.py

Multi-generator inner fluctuation: extends Comp 74's single-generator CC inner fluctuation A = iφ[D_sub, ω]/2 to K coherent generators ω_k = e_{2k-1}e_{2k}. With K active pair-generators: D_φ² = D + 2Kφ², spectral action S/2^D = e^(-1) exp(-2K φ²/D), φ⁴ coefficient = 2K²e^(-1)/D². Matching Comp 72 target φ⁴ coefficient e^(-1) y⁴/(2D²) gives y/Λ = √(2K/D). For y/Λ = 2^(-1/4) (D-independent Z² target): K = D/(2√2). NOT INTEGER for any finite D. Bernoulli extension: each ω_k active with probability p, asymptotic E[S/2^D] → e^(-1) exp(-p φ²), gives y/Λ = √p. Hits target for p = 1/√2 (self-consistent: p = (y/Λ)²) but 1/√2 is not directly derivable from μ_site = 1/2 (substrate measure on bits, not generator-pairs). The multi-generator inner-fluctuation route has the right scaling form but no structurally natural value of K or p delivers the target; the partition-function-level identification in Comps 87–92 takes a different substrate-side route that does close.

Computation 85 — A_F-internal Higgs inner fluctuation for the Z² identification
.py.txtcomputation_85.py

A_F-internal Higgs fluctuation: place the Higgs in the canonical A_F = ℂ ⊕ ℍ ⊕ M_3(ℂ) finite-internal factor (as in Chamseddine-Connes-Marcolli SM construction) rather than in the Clifford volume element of Comps 74, 81. The CCM moment-extraction formula λ ~ (a/π²) f_0 maps f_0 (Seeley-DeWitt moment of cutoff f) to the Higgs quartic. For PST substrate at matched scaling, all eigenvalues at ±√D collapse f_0 → f(1) = e^(-1) (Comp 62 result). The conjecture restates: does the CCM moment-coefficient extraction generalise to the discrete substrate triple with f_0 → S_sub/2^D? This approach does not provide a new computational test — it restates the conjecture at the moment-extraction level. Identifies the structural obstacle for the heat-kernel route: the formalism that maps f_0 to λ is built on continuous-spectrum asymptotics, and the substrate's collapse-to-cutoff spectrum makes this formalism degenerate. The partition-function-level identification in Comps 87–92 takes a different substrate-side route.

Computation 86 — Seeley-DeWitt heat-kernel analysis on the discrete substrate triple
.py.txtcomputation_86.py

Seeley-DeWitt heat-kernel analysis: examines the Seeley-DeWitt heat-kernel expansion Tr f(D²/Λ²) ~ Σ_n f_n Λ^{d-2n} a_n for the substrate spectral triple and identifies the structural obstacle. For continuum d = 4 spectral triple: f_0, f_2, f_4 moments multiply a_4, a_2, a_0 coefficients delivering the SM Higgs quartic from a_4. For PST substrate at matched scaling Λ = √D: all eigenvalues at ±√D, so Tr f(D_sub²/Λ²) = 2^D · f(1) = 2^D · e^(-1) exactly. Seeley-DeWitt coefficients collapse: a_0 = 2^D (non-zero), a_n = 0 for n ≥ 1 (no curvature contributions from a degenerate spectrum). The 'moment' structure that delivers λ in continuum CCM is structurally absent — no continuous-spectrum a_4 term to extract. STRUCTURAL Conclusion: closing Z² via standard CCM heat-kernel machinery is not possible at matched scaling. The substrate spectral triple is over-constrained: the (single) eigenvalue magnitude encodes the entire spectral action, leaving no room for an independent Higgs-quartic coefficient. A fundamentally new spectral-action formalism for discrete substrate triples would be required.

Computation 87 — Bernoulli Laplace-transform identification (bypassing the CC inner-fluctuation route)
.py.txtcomputation_87.py

Bernoulli Laplace-transform identification: bypass CC inner fluctuation via the asymptotic Laplace transform of the substrate Bernoulli empirical mean. For X̄ = |C|/D (fraction of occupied sites): E_μ[exp(c X̄)] = ((1 + exp(c/D))/2)^D → exp(c/2) as D → ∞. For Z² = e^(-1): need c/2 = -1, i.e. c = -2. Numerical convergence at D = 10000: E[exp(-2 X̄)] = 0.367898 vs target 0.367879, error ~ 1.84e-5. Multiple structural readings of c = -2 (c = -1/μ_site, c = -4σ² × 2, c = -2 states per bit, c = -2σ tilting) are form-compatible but NONE uniquely forced by P1 alone. Conclusion: form-compatible recognition delivering e^(-1) but the multiplicative constant -2 is not uniquely derived. Comp 88 sharpens this with the KO-dimension reading.

Computation 88 — KO-tempered Bernoulli MGF identification of e^(-1) at the matched scaling
.py.txtcomputation_88.py

SHARPENS Comp 87. The Laplace-transform scale c = -2 admits a STRUCTURALLY FORCED reading: the PST total KO-dimension is 4 (M = ℝ × S³) + 6 (internal A_F factor) = 10, and 10 mod 8 = 2 by Bott periodicity (§sec:foundational-object). Setting c = -(KO_total mod 8) = -2 gives the structurally-forced identification: e^(-1) = E_μ[exp(-(KO_total mod 8) · X̄)] asymptotically. ALL THREE INGREDIENTS structurally derived: μ_site = 1/2 from P1, KO_total mod 8 = 2 from 10-foundational + Bott, asymptotic limit from Cramér-Bernoulli. This is a second, independent substrate-side derivation of e^(-1): not just S_sub/2^D = f(1) (Comp 62, spectral-action route) but also E_μ[exp(-(KO mod 8) X̄)] = e^(-1) (Bernoulli-MGF route). REMAINING GAP NOW SHARP: identifying the SM Higgs effective coupling at M_* with the substrate's KO-tempered Bernoulli partition function Z_H(β = KO mod 8). This is a NEW substrate-to-SM bridge parallel to (and bypassing) the obstructed CC inner-fluctuation bridge. Worth pursuing as an independent research direction. STATUS: Z² closure REDUCED to one specific identification statement — the sharpest open formulation of the Z² conjecture to date.

Computation 89 — Z² bridge step 1: substrate Higgs Hamiltonian derivation
.py.txtcomputation_89.py

+ Closes the SUBSTRATE SIDE of the Z² KO-tempered partition-function bridge. Derives H_Higgs(C) on substrate configurations from three structural constraints: (i) ADDITIVITY (Bernoulli sites independent ⇒ H = Σ_a h_a(B_a)); (ii) MATCHED SCALING (per-site energy 1/D at Λ = √D); (iii) UNIFORMITY (P1 symmetry ⇒ h_a uniform). Together force H_Higgs(C) = (1/D) Σ_a B_a = X̄(C), the empirical mean of bit-occupation — NO FREE PARAMETERS. With H_Higgs derived, the KO-tempered partition function Z_H(β_KO) = E_μ[exp(-(KO_total mod 8) X̄)] → e^(-1) asymptotically. Substrate side of bridge is now closed structurally. Monte Carlo verification at D = 10, 100, 1000 confirms convergence to e^(-1). Bridge identification: λ_SM(M_*) = b · Z_H(β_KO) = (1/4) e^(-1) ≈ 0.0920; observed λ_SM(M_*) ≈ 0.0926; ratio 1.008 — the same 0.8% near-coincidence as Z². REMAINING OPEN (SM side): derive the partition-function-level Connes-Chamseddine correspondence between substrate triple at matched scaling and the SM Higgs effective coupling. This is the active the active research direction. Comp 89 status: SUBSTRATE SIDE OF BRIDGE CLOSED.

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

Reduces the SM-SIDE of the Z² bridge to a choice among three concrete candidate formulations, each with a specific gap. (I) WILSONIAN EFFECTIVE-COUPLING RG: integrate out substrate fluctuations between M_* and Λ = √D; effective λ_eff(M_*) = b · Z_H(β_KO) follows IF the KO-temperature arises from matched-scaling rescaling. Gap: derive the discrete-substrate Wilson-Polchinski equation. (II) SPECTRAL-ACTION VARIATION: direct evaluation of S(D + A)/Λ^4 for inner fluctuation A = H_Higgs at matched scaling gives e^(-2) not e^(-1). Naive form fails; needs sharpened normalisation. (III) SUBSTRATE-MEASURE INVARIANCE: the matched-scaling map Π identifies h(C) = X̄(C); expanding V_eff around the substrate-side vacuum X̄ = 1/2 gives the additive form V_eff(C)/T_KO = β_KO · X̄(C). Gap: derive T_KO ↔ M_* identification at SM-side vacuum. Numerical verification: λ_predicted = b · e^(-1) ≈ 0.0920 vs λ_observed ≈ 0.0926 (Buttazzo 2013), ratio 1.007 (the same 0.8% near-coincidence). Conclusion: SM-side bridge content REDUCED to one of three specific calculations, all well-posed. Formulation (II) fails in naive form; (I) and (III) survive as the active the active research direction.

Computation 91 — Z² empirical-precision test (two-loop SM RGE)
.py.txtcomputation_91.py

Tests the empirical tightness of Z² ≈ e^(-1) via SM RGE precision. Runs from v to M_* = 4π m_h at one-loop and adds a compact two-loop correction. RESULTS: simplified one-loop gives λ_SM(M_*) ≈ 0.087 (Z²/e^(-1) ≈ 0.94, 5% deviation); compact two-loop ≈ 0.085 (Z²/e^(-1) ≈ 0.93, 7% deviation). The paper's 0.8% near-coincidence claim depends on FULL Buttazzo-level SM RGE precision (two-loop + threshold corrections at m_t), not simplified RGE. All three readings (one-loop, compact two-loop, full Buttazzo) are consistent with 'structural identity within RGE truncation uncertainty.' Sensitivity analysis: variation of m_h ± 1 GeV and m_t ± 2 GeV gives Z²/e^(-1) in range 0.89-0.97 at simplified one-loop. STRUCTURAL Conclusion: SM RGE precision at simplified levels insufficient to test Z² = e^(-1) at proof level. The empirical near-coincidence depends sensitively on the SM RGE precision; closing the structural identification requires either full three-loop SM RGE or the partition-function bridge of Comp 90 formulations I/III (RGE-truncation independent). Substrate-side derivation (Comp 89) remains exact and independent of SM RGE.

Computation 92 — Z² bridge formalisation: 7-step structural proof sketch
.py.txtcomputation_92.py

Closure of the proof sketch level. Reduces the Z² = e^(-1) closure to a SINGLE WILSONIAN-MATCHING PREMISE: at M_*, the SM Higgs effective coupling equals the PST bare LG quartic times the substrate's KO-tempered Bernoulli partition function. Establishes a 7-step structural derivation: (1) PST EFT below M_* = pure SM [paper eq:eft]; (2) Wilsonian matching at M_*; (3) PST bare coupling λ_PST^bare = b = 1/4 [LG modal potential]; (4) substrate UV at Λ = √D contributes vacuum fluctuations to M_*; (5) fluctuations integrated over μ contribute Z_H(β) = E_μ[exp(-β H_Higgs)] with H_Higgs = X̄ [Comp 89]; (6) temperature β_KO = KO_total mod 8 = 2 [Comp 88, Bott periodicity]; (7) therefore λ_SM(M_*) = b · Z_H(β_KO) → (1/4) · e^(-1) ≈ 0.0920, matching observed 0.0926 at 0.8%. Status of each step: (1) paper sec:renorm + eq:eft; (2) standard EFT machinery; (3) paper sec:foundational-object; (4) matched-scaling A1; (5) Comp 89; (6) Comp 88; (7) follows from 1-6. ELEVATES Z² closure status from 'open conjecture' to 'structural proof sketch established.' Remaining work: embed steps in unified partition-function-level Connes-Chamseddine correspondence framework where each step becomes a formal theorem.

Computation 93 — M_* off-by-one rescue (b) re-examined: fourth-rescue candidate
.py.txtcomputation_93.py

Re-examines rescue (b) of Comp 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 (which is captured separately by C_H = 4λ_PST = 1 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) = (M_*²/16π²)(1 + 0) = M_*²/(16π²), recovering m_h = M_*/(4π) as parameter-free one-loop derivation. Identifies the fourth-rescue candidate to close the M_* off-by-one; full validation via projection-chain analysis follows in Comp 97.

Computation 94 — Bridge Premise (B) fluctuation-matching T_KO (research direction)
.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 (KO mod 8) delivers the testable prediction µ²(M_*) = 2D · M_*². Acknowledged as dimensionally suspect (X̄ dimensionless vs h with GeV units); held as a research-direction sketch rather than a closure.

Computation 95 — Bridge Premise (B) reduced to structural claim (***) via substrate-Wilsonian RG
.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 exp(-β_KO · H_Higgs(C)) on substrate configurations — the partition-function-level analogue of the heat-kernel RG flow generator exp(-tD²) for the continuum spectral triple. If (***) is granted, λ_SM(M_*) = b · Z_H(β_KO) = e^(-1)/4 ≈ 0.0920 follows structurally.

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

Polchinski-style derivation attempted on the discrete substrate. Establishes the partition-function form structurally via configuration-Hamiltonian Boltzmann averaging Z_D(β_KO) = ((1 + exp(-β_KO/D))/2)^D → e^(-1) 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 on item 1.1 of open research.

Computation 97 — M_* off-by-one: projection-chain validation CLOSES item 1.2
.py.txtcomputation_97.py

Validates Comp 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_3(ℂ)). Including the zero mode in N_B does NOT double-count C_H = 4λ = 1. Total δm_h² = M_*²/(16π²), restoring m_h = M_*/(4π) as parameter-free one-loop derivation. Item 1.2 (M_* off-by-one) CLOSED at proof-detail level; the 'one-loop-complete / parameter-free' reading of M_* = 4π m_h is restored.

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

Closes items (a), (b) of Comp 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_0 ⊕ V_1 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 β_KO · X̄ is a TRIVIAL FIXED POINT (Theorem 1). Uniqueness follows from Comp 89 boundary condition. Residual item (c) remains as the matched-scaling identification (subsequently closed in Comp 100).

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

Exploratory closure of item (c) via thermal-field-theory framing: postulates β_KO = KO mod 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^(-1). The 'KO-thermal principle' is at the same level as the spectral-action principle in standard CC — a framework postulate beyond P1-P3. Comp 99 is HELD AS EXPLORATORY, not as a closure: PST's claim is that all structure derives from P1-P3 alone, so a CC-framework postulate weakens that claim. SUPERSEDED by Comp 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
.py.txtcomputation_100.py

Closes the SUBSTRATE-SIDE content of item (c) of Comp 96 from P1-P3 via tensor-product structural argument. 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 (factor as product over independent bits), f must satisfy f(x+y) = f(x) f(y), unique smooth solution f(x) = exp(-x). At matched scaling Λ² = D: (1/2^D) Tr exp(-Δ/D) = ((1 + exp(-2/D))/2)^D → e^(-1). The factor '2' comes from one-bit Clifford structure, NOT KO mod 8 (numerically equal but structurally independent). Comp 99's KO-thermal postulate SUPERSEDED. SCOPE NOTE : the matching identification λ_SM(M_*) = b · Z_H(β_KO) — bridge premise (B) itself — remains the single open foundational step; standard Wilsonian quartic matching is additive in loop corrections, not multiplicative. Conditional on (B): λ_SM(M_*) = b · e^(-1) = 0.0920 matches observed 0.0927 at 0.8% Buttazzo-RGE precision.

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

Attacks the concern that standard Wilsonian threshold matching of a quartic is additive in loop corrections, not multiplicative. Reinterprets λ_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 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 exp(-β_KO · H_Higgs(C))? (2) why is the vertex renormalisation trivial (Z_Γ4 = 1)? Both are tractable extensions of Comp 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.

Computation 102 — Item 1.2 resolution: the two '2's are structurally independent
.py.txtcomputation_102.py

Addresses the concern that the identification c = −(KO_total mod 8) = −2 in the Bernoulli MGF is asserted, not motivated. Comp 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 mod 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 Comp 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; remaining open content is external-validation / literature-search task on whether partition-function-level CC has standing in the spectral-action literature (not new research).

Computation 103 — Bridge premise (B) attack: explicit path-integral Jacobian (negative result)
.py.txtcomputation_103.py

Attempts Comp 101's sub-question (1): is the matched-scaling Π's path-integral Jacobian equal to exp(-β_KO · H_Higgs(C))? Two natural interpretations both fail: (A) linear-map determinant of Π : V_0 ⊕ V_1 → SM_low is configuration-independent (constant); (B) Cramér-Bernoulli LDP near the vacuum X̄ = 1/2 gives Gaussian weight in (X̄ − 1/2), not exponential in X̄ (numerically verified I(x) = x ln(2x) + (1-x) ln(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 the with a concrete negative result. Sub-question (2) [Z_Γ4 = 1] is satisfied at tree level (standard CC inner fluctuations). Net: item 1.1 is RE-CHARACTERISED, not closed. Bridge premise (B) is a structural postulate of the partition-function-level CC framework at the same level as the spectral-action principle in standard CC. HONEST CLOSURE STATUS: PST has P1-P3 + partition-function-level CC postulate; substrate-side derivation of e^(-1) is real (Comp 100), SM-side matching is the structural input.

Computation 104 — Bridge premise (B) factored: substrate side derived, SM side inherited from standard CC
.py.txtcomputation_104.py

Attempts to extend Mosco convergence from Dirichlet forms to spectral actions, with goal of deriving (B) at same level Mosco delivers SM kinetic term. Structural analysis shows (B) FACTORS into two parts: (B-substrate) (1/2^|D|) Tr f(Δ/|D|) → e⁻¹, CLOSED at proof-detail level from P1-P3 by Comp 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 (Chamseddine-Connes 1996, Connes-Marcolli 2008, follow-up work), 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); it does not invent the SM-side framework. The standing procedural question is addressed: YES, the framework PST uses has literature standing. PST's foundational claim refines to 'P1-P3 + standard CC spectral-action framework' — substrate-side contribution is the genuinely new PST result (Comp 100), SM-side framework inherited from established CC literature.

Computation 105 — Two-layer one-loop diagram: no double-counting of M_* zero modes
.py.txtcomputation_105.py

Verifies non-double-counting between the substrate χ_∅ singleton-cylinder loop (Layer A, Λ = √D → M_*) and the emergent SM Higgs C_H loop (Layer B, M_* → IR) via three independent criteria: (i) disjoint momentum support — the two integration shells meet only at the single matching scale M_* (measure zero in both integrals); (ii) distinct dynamical variables — χ_∅ on P(D) vs C_H on M = ℝ × S³; (iii) distinct counterterm structures — the matching is a boundary condition between the two layers, not a loop integration. Numerical verification: the substrate full-trace partition function (Comp 100) decomposes EXACTLY as Z_sub^full = Z_sub^IR-shell + Z_sub^UV-shell across the M_*-shell partition with machine-zero residual at D ∈ {10, 20, 50, 100}. The substrate-vs-emergent distinction of Comps 93/97 is therefore robust at the diagram level.

Computation 106 — SM gauge block-diagonality on ⊕_k G_k (no tree-level FCNC)
.py.txtcomputation_106.py

Explicit proof that SM gauge generators preserve each generation sector in the Dixon-synthesis decomposition V_Dixon = V_core ⊕ G_1 ⊕ G_2 ⊕ G_3. Three-part proof: (1) SU(3)_c is the M_3(ℂ) factor of A_F acting on the colour triplet WITHIN one Fano-plane sub-algebra 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, 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 O(1) off-block magnitudes. The absence of tree-level flavour-changing neutral currents is therefore 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 is sector-blind (closes 1.3 (i)+(ii))
.py.txtcomputation_107.py

Closes sub-tasks (i) configuration-to-particle map and (ii) cross-sector coupling universality at M_* of open-research item 1.3 (full-SM extension of the boson-fermion cancellation) via a single observation: the substrate-level binomial cancellation N_B - N_F = 0 is SECTOR-BLIND. The projection chain Π : χ_S → field on M = ℝ × S³ → SM particle via A_F inner fluctuations is surjective onto SM at M_* by construction (paper eq:eft; A_F via Comps 3, 19, 106; three generations via Comps 46-48). Substrate count is invariant under ANY partition of the 2^D Walsh modes into hypothetical SM sectors — per-sector totals vary widely (IR Yukawa hierarchy expressed at substrate level) but the grand total stays zero by binomial identity. Numerical verification across D ∈ {4, 6, 8, ..., 20} confirms exact cancellation; random sector partitions at D = 12 give grand total = 0 across all trials. Empirical SM Veltman residual at M_* ≈ 1.573 TeV is ≈ -4 (IR Yukawa hierarchy expressed at M_*, not a substrate contradiction). The explicit χ_S → SM-particle map per particle is downstream proof-detail work, not required for the cancellation argument.

Computation 108 — Inner-fluctuation gauge bosons are bosonic (closes 1.3 (iii))
.py.txtcomputation_108.py

Closes sub-task (iii) of open-research item 1.3: inner-fluctuation gauge bosons carry BOSONIC statistics by three independent readings. (R1) Algebraic: Connes 1-forms A = Σᵢ aᵢ [D, bᵢ] for aᵢ, bᵢ ∈ 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 (automatically bosonic). (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 dim Cl_even(P(D)) = 2^(D-1) accommodates these at D ≥ 6 with excess capacity (lifts to higher-derivative corrections at M_*⁻¹ suppression). Combined with Comp 107, closes open-research item 1.3 at the structural level.

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.

v26.15