Chapter 14: Protocol B (Mass Spectrum Computation)

The derivation of the lepton mass hierarchy in Chapter 8 identifies the muon and tau as the 11th and 17th harmonics of the electron's dodecahedral cage. To audit these results against the engine ledger, we employ Protocol B: a deterministic verification of the projected lattice geometry. This protocol computes the charged-lepton mass ratios from the Tier 2 harmonic-ladder mechanism plus the Tier 3 integer anchors $N=11$ and $N=17$ per Parameter Ledger §0.1; first-principles unique-eigenvalue extraction remains open under O.5/O.14/O.15.

---

14.1 Methodology

14.1.1 Algorithm: Projected Lattice Dynamics

The analysis models the vacuum as an Ammann-Kramer-Neri (AKN) tiling, generated via the Cut-and-Project method from the 6-dimensional hyper-lattice $\mathbb{Z}^6$. We construct the Geometric Resonance Matrix ($\mathbf{G}$), which represents the coupled harmonic structure of the $N=144$ dodecahedral cage.

Unlike standard molecular dynamics, which operate in 3D Euclidean space, the GCT algorithm operates on Projected Degrees of Freedom. Each node $i$ in the 6D parent lattice possesses a 6D displacement vector $\mathbf{U}_i$. Upon projection into the physical manifold, this vector is partitioned into physical ($\mathbf{u}_\parallel$) and internal ($\mathbf{u}_\perp$) components. The interaction energy of a bond between nodes $i$ and $j$ is defined by the stiffness hierarchy derived in Volume 2: $V_{ij} = \frac{1}{2} K_\parallel (\Delta \mathbf{u}_\parallel \cdot \hat{\mathbf{n}}_\parallel)^2 + \frac{1}{2} K_\perp (\Delta \mathbf{u}_\perp \cdot \hat{\mathbf{n}}_\perp)^2$ The analysis solves for the Node-Centric Eigenvalues (eigenstates of the center-of-mass motion), distinguishing these from the interface stress energies (surface area mismatches) utilized for the hadron spectrum.

14.1.2 Parameters: The Phason Stiffness Ansatz ($D=18$) [Tier 2 postulate + Tier 3 specific exponent — Parameter Ledger §0.1 P3; Open Problem O.15]

The analysis requires exactly one material input: the ratio of phason stiffness to phonon stiffness.

Theorem (normalization convention): The Gram eigenvalue ratio of the 6D$\\to$3D icosahedral projection gives a tree-level stiffness ratio $\\eta_0 = \\phi^{-2}$ once the parallel-sector kinetic term is normalized to the lattice-unit convention used in App K. Proof: Before normalization, the parallel Gram factor contains the golden-ratio contribution $1+\phi^{-1}=\phi$, not $2$. The Protocol B convention rescales the parallel kinetic term to the unit phonon stiffness used for $K_\parallel$, so the displayed tree-level phason-to-phonon stiffness ratio is the normalized lattice-unit ratio $\eta_0=\phi^{-2}$. The equality $1+\phi^{-1}=\phi$ is the algebraic identity; the factor $2$ belongs only to the chosen kinetic-term normalization.

In 3D icosahedral quasicrystals, phonon-phason coupling drives K_\\perp under RG flow toward zero (Lubensky et al. 1985). The stiffness is suppressed from $\\eta_0 = \\phi^{-2}$ by the phonon-phason coupling at the relevant physical scale.

We set $n=18$ as a Tier 2 postulate + Tier 3 specific exponent (Parameter Ledger §0.1 P3 canonical disposition; first-principles RG closure tracked under Open Problem O.15). This equals $(\\eta_0)^9 = (\\phi^{-2})^9$, consistent with 9 phonon-phason coupling channels (3 perp dimensions $\\times$ 3 phonon polarizations) — one of three heuristic motivations for the exponent (cube-of-Gram-determinant ratio per App K §K.3, $H_3$-Coxeter rank × Galois × dim per App K §K.4, and the 9-channel-squared count here). Empirical support: measured phason stiffness in i-AlMnSi and i-AlPdMn is $\\sim 10^{-2}$ to $10^{-1}$ in raw $K_{\\text{phason}}/K_{\\text{phonon}}$ ratio [Tier 3 — quasicrystal materials measurements], 100-1000× larger than the bare GCT prediction $\\phi^{-18} \\approx 1.73 \\times 10^{-4}$; the Penrose-Toner chemical-bonding sweep in App H O.15 / App K §K.4b quantifies the metallic-alloy floor at $K_{PT}/K_{phonon}\in[9.3\\times10^{-3},1.09\\times10^{-1}]$, while the first-principles RG derivation of the correction remains open.

The input parameter for the phason stiffness is set as a rigid structural property of the icosahedral vertex: \\frac{K_\\perp}{K_\\parallel} = \\phi^{-18} \\approx 1.733 \\times 10^{-4}

14.1.3 Computational Setup: Harmonic Ansatz Evaluation

The $6N \times 6N$ geometric projection matrix is analyzed to extract the spectrum of eigenfrequencies $\{\omega_k\}$. Because the density of states in a quasicrystal is singular continuous, eigenvalues are clustered in a complex hierarchy. Protocol B utilizes a Harmonic Ansatz Evaluation to target the specific pseudo-gaps ($\phi^{11}$ and $\phi^{17}$) without the algorithm becoming lost in the dense background of uncoupled phason modes.

The evaluation identifies the fundamental breathing mode resonance ($\omega_0$), which represents the electron's mass-energy.

14.1.4 Numerical Convergence and the Local Isomorphism Shield

Stability is probed by increasing the cluster size from the $N=144$ core toward an $N=2000$ node neighborhood. Because the dodecahedral cage satisfies the local matching rules of the icosahedral projection, it is protected by the Local Isomorphism Theorem. The aperiodic background acts as a Topological Shield; while every individual node is unique, the $N=144$ symmetry acts as a resonant cavity that filters out local fluctuations, yielding a monochromatic mass eigenvalue protected against local environmental disorder [Tier 2 mechanism — Local Isomorphism Theorem]. A specific numerical convergence rate is not asserted here: the finite-$N$ discretization artifact present in the direct graph extraction (§14.2; protocol_lepton_spectrum.py) is why first-principles non-perturbative eigenvalue closure remains open (O.5/O.14/O.15).

---

14.2 Validation and Results

[!CAUTION] The Linear Boundary Constraint — read this before the precision tables below. Extensive HPC analysis of the full sparse Hessian (see Appendix Q) confirms that the bare linear dynamical matrix yields normalized eigen-ratios of [1.0, 1.67, 2.1, 2.52...]. Neither $\phi^{11}$ nor $\phi^{17}$ emerge spontaneously from linear diagonalization.

The leptons are non-linear phason solitons whose exact mass scaling is governed by the $6N \pm 1$ symmetry saturation, not linear oscillatory modes. The $\phi^{11}$ and $\phi^{17}$ gap alignments are therefore Tier 2 Geometric Postulates (the harmonic-ladder structure committed to in Parameter Ledger §0.1 P4/P5) plus the Tier 3 specific integer choices $N=11$ (muon) and $N=17$ (tau). The non-perturbative computational extraction of these solitonic eigenvalues from first principles remains an open goal for the GCT framework. The ppm-level precision figures in §14.2.3 below are scored against pole-mass values that already include SM 2-loop EW + 3-loop QED corrections (see App R §R.2.1 Loop-Order Discipline); the bare Tier 2 geometric precision is ~0.25%, with the sub-percent figures requiring the GCT geometric form combined with the SM-equivalent radiative-correction provenance.

14.2.1 The Harmonic Ansatz and the Linear/Non-Linear Boundary [Tier 2 Postulate]

The primary non-linear pseudo-gaps of the dodecahedral cage correspond to the 11th and 17th harmonics of the inflation operator. Under the Geometric Harmonic Ansatz:

• Muon Target: $\phi^{11} \approx 199.005$ [Tier 2 Postulate]
• Tau Target: $\phi^{17} \approx 3571.34$ [Tier 2 Postulate]

14.2.2 Precision Mass Audit

We apply the phason drag corrections derived in Chapter 8 to the geometric eigenvalues. For the Muon, we include the second-order Electroweak Mixing Correction ($\phi^8\alpha^2$ [Tier 2], verified computationally; see Appendix Q). The residual error of 21 ppm ($0.0021\%$) for the Muon (lepton-spectrum protocol; App R §R.1) sits within the $\sim$10–40 ppm higher-loop theory floor disclosed for the SM-loop-equivalent comparison discipline (App R §R.2.1). The next geometric correction enters at third order carrying a $\phi$-power coefficient comparable to the $\phi^8$ second-order term — $\phi^8\alpha^3 \sim 1.8 \times 10^{-5}$ [Tier 3 — order-of-magnitude estimate by analogy to the second-order coefficient], of the same scale as the residual — whereas the bare uncoefficiented Third-Order Radiative Correction $\alpha^3 \sim 3.9 \times 10^{-7}$ [Tier 1 — QED perturbation theory] alone is two orders smaller. Native first-principles loop closure remains open under O.5/O.14/O.15.

14.2.3 Final Lepton Mass Table vs. PDG

14.2.4 The Tau Pre-diction

It is critical to note that for the Tau mass, GCT is currently more precise than experimental measurement. While the Particle Data Group (PDG) uncertainty is $\pm 0.12$ MeV ($\sim$0.007%), GCT predicts the mass to five decimal places. Protocol B thus serves as a Targeted Pre-diction for future high-luminosity lepton colliders.

Conclusion: Protocol B treats the charged lepton mass hierarchy as a deterministic spectral fingerprint candidate of the 6D hyper-lattice. The precision of the match across multiple generations using the icosahedral geometric framework (parameterised by $\phi$ as the structural constant, plus the integer-valued harmonic exponents $N = 11, 17$ which enter as Tier 3 specific values per Parameter Ledger §0.1 P4/P5 — not as additional fitted parameters but as structural commitments of the framework) provides conditional support for the harmonic-ladder mechanism while first-principles unique-eigenvalue extraction remains open under O.5/O.14/O.15. References to $\phi$ as the structural constant are structural-constants-only language; the full 5-postulate-plus-1-anchor bare gauge+lepton sub-sector, expanding to 5-postulate-plus-3-anchor when native-RGE endpoint and measurement-anchored precision-comparison rows are included (Ledger §0.1), is the load-bearing tier discipline.