The Geometry of the Standard Model:
Deriving the Higgs Mass, Lagrangians, and
Gravity Echoes from Lattice Saturation
Raghu Kulkarni
Independent Researcher, Calabasas, CA
raghu@idrive.com
February 21, 2026
Abstract
We propose that the Standard Model Lagrangian is the continuum limit of a dis-
crete vacuum geometry. Using the framework of a saturated Cuboctahedral Lattice
(K = 12), we demonstrate how the fundamental sectors of particle physics emerge
from geometric first principles:
1. The Lagrangian Sector: We review how the standard Klein-Gordon (Scalar),
Dirac (Spinor), and Yang-Mills (Gauge) Lagrangians naturally emerge as the
continuum elastic limits of Face-Centered Cubic (FCC) lattice tension, topo-
logical braiding, and plaquette preservation. Crucially, we formalize how the
Non-Bipartite topology of the simplicial lattice resolves the Nielsen-Ninomiya
“Fermion Doubling Problem,” allowing chiral fermions to exist without spuri-
ous mirrors.
2. The Higgs Sector: Using the holographic ratio of the punctured topological
boundary surface (108) to the saturated bulk volume (1728), we derive a ge-
ometric coupling λ ≈ 0.125, predicting a bare Higgs mass of 123.11 GeV. We
explicitly derive the necessary factor of 2 from the bi-directional topological
flux across shared lattice interfaces.
3. The Gravity Sector: We interpret the Event Horizon as a physical phase
boundary where the lattice melts from a K = 12 saturated crystal to a dis-
ordered fluid. This topological melting slows gravitational wave propagation
to v ≈ c/4. Applying this to the recent LIGO detection GW250114, we pre-
dict a prompt scattering echo delay of ∆t ≈ 2.7 ms, seamlessly matching the
post-merger overtone window.
Contents
1 The Vacuum as a Crystal 3
2 Review: The Scalar Sector and Inertia 3
2.1 The Klein-Gordon Continuum Limit . . . . . . . . . . . . . . . . . . . . . 3
2.2 Geometric Origin of Mass . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1
3 Review: The Spinor Sector and Matter 4
3.1 The Dirac Equation on a Lattice . . . . . . . . . . . . . . . . . . . . . . . 4
3.2 Resolution of the Fermion Doubling Problem . . . . . . . . . . . . . . . . . 4
3.2.1 The Hypercubic Trap (Bipartite Symmetry) . . . . . . . . . . . . . 4
3.2.2 The FCC Solution (Non-Bipartite Topology) . . . . . . . . . . . . . 5
3.2.3 Dispersion Analysis at the Zone Boundary . . . . . . . . . . . . . . 5
4 Review: The Gauge Sector and Plaquettes 5
4.1 Yang-Mills Action from Wilson Loops . . . . . . . . . . . . . . . . . . . . . 6
4.2 The Geometric Origin of the Plaquette . . . . . . . . . . . . . . . . . . . . 6
5 The Higgs: Freezing the Lattice 6
5.1 The Geometric Origin of λ . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
5.2 The Mass Prediction and RG Flow . . . . . . . . . . . . . . . . . . . . . . 7
6 Falsifiability: Gravitational Echoes 7
6.1 The Vacuum Phase Transition . . . . . . . . . . . . . . . . . . . . . . . . . 7
6.2 The Echo Prediction (2.7 ms) . . . . . . . . . . . . . . . . . . . . . . . . . 8
6.3 Comparison with Observation . . . . . . . . . . . . . . . . . . . . . . . . . 8
7 Conclusion 8
2
1 The Vacuum as a Crystal
The quest for a unified theory faces a persistent obstacle: the incompatibility between
the smooth manifold of General Relativity and the quantized nature of particles. Modern
Quantum Field Theory (QFT) treats space-time as a continuous background, leading to
singularities (UV divergences) that require renormalization. Furthermore, the Standard
Model cannot internally explain the origin of its own continuous parameters (m
h
≈ 125
GeV, v ≈ 246 GeV).
We propose a geometric paradigm: Geometric Saturation, formalized in the Selection-
Stitch Model (SSM) [1]. We posit that the vacuum is a discrete physical mesh of tetra-
hedra. This entanglement process naturally saturates at a strict geometric coordination
number of K = 12 (the Kepler limit for sphere packing), creating a rigid face-centered
cubic (FCC) structure. If the vacuum is a discrete crystal, then “particles” are simply its
macroscopic topological and vibrational modes.
2 Review: The Scalar Sector and Inertia
While the derivation of continuous Lagrangians from discrete lattice Hamiltonians is a
well-established discipline (e.g., Wilson [5], Creutz [6], Rothe [7]), we briefly review these
mechanics within the specific context of the FCC lattice to establish notation for the novel
Higgs and Gravity derivations in Sections 5 and 6.
We model the vacuum as a coupled harmonic system. Let ϕ
n
(t) represent the geometric
deviation of a node n from its equilibrium position.
2.1 The Klein-Gordon Continuum Limit
The Hamiltonian is defined by the Kinetic Energy (time evolution) and Stitch Tension
(gradient energy):
H
scalar
=
X
n
"
1
2
m
0
˙
ϕ
2
n
+
1
2
κ
X
m∈neigh
(ϕ
n
− ϕ
m
)
2
+ V
geo
(ϕ
n
)
#
(1)
To recover standard physics, we take the continuum limit where lattice spacing a → 0.
The discrete difference operator transforms into the Laplacian via Taylor expansion:
ϕ
n+1
− ϕ
n
≈ a · ∇ϕ +
1
2
a
2
∇
2
ϕ (2)
Substituting this into the potential term, we recover the continuum action. Setting the
wave speed c
2
= κa
2
/m
0
= 1 (natural units), we obtain the exact Klein-Gordon La-
grangian:
L
scalar
=
1
2
(∂
µ
ϕ)(∂
µ
ϕ) −V
eff
(ϕ) (3)
2.2 Geometric Origin of Mass
In standard QFT, bare mass is an arbitrary coupling parameter. In a discrete framework,
mass is a geometric constraint derived strictly from the lattice topology [2]. The lattice
potential V
geo
enforces the integer topology of the mesh (K = 12). A node cannot exist
3
“between” integer sites. This creates a periodic potential well. Thus, mass is identified
strictly as the resonance frequency of a defect vibrating within the geometric trap of the
vacuum.
3 Review: The Spinor Sector and Matter
Standard scalar fields describe bosons (force carriers). To describe Fermions (Matter), we
must model defects as Topological Braids (Twists) that possess orientation and cannot
be continuously deformed to zero.
3.1 The Dirac Equation on a Lattice
Lattice hopping for a directed twist requires a mechanism to preserve directional infor-
mation (Spin). A standard Laplacian ∇
2
loses this information. We effectively “take the
square root” of the Laplacian using the Clifford algebra of the tetrahedral basis vectors.
We define the Directed Stitch Operator D:
D = −iγ
µ
∂
µ
(4)
where γ
µ
satisfy the anticommutation relation {γ
µ
, γ
ν
} = 2g
µν
. This ensures that D
2
=
−∇
2
, recovering the correct energy-momentum relation (E
2
= p
2
). The equation of
motion for a braid with topological tension M is the Dirac Equation:
(iγ
µ
∂
µ
− M)ψ = 0 (5)
3.2 Resolution of the Fermion Doubling Problem
A major theoretical obstacle for discrete field theories is the Nielsen-Ninomiya “No-Go”
Theorem, which dictates that discretizing the Dirac equation on a standard hypercubic
lattice unavoidably produces spurious “doubler” fermions (mirror particles) at the corners
of the Brillouin zone [8].
It is a known phenomenon in condensed matter physics that non-bipartite (triangu-
lar/simplicial) lattices can evade this symmetry trapping. We formalize this mechanism
here for the K = 12 vacuum geometry, based on the rigorous derivations in [3].
3.2.1 The Hypercubic Trap (Bipartite Symmetry)
On a standard square/cubic lattice, the discrete Dirac operator in momentum space is:
D
cubic
(k) = iγ
µ
sin(k
µ
a)
a
(6)
The energy spectrum contains zeros whenever sin(k
µ
a) = 0. This occurs at the physical
origin (k = 0) but also at the zone boundary (k = π/a). The k → k + π symmetry is a
direct mathematical consequence of the strictly Bipartite Nature of the hypercubic graph
(checkerboard coloring).
4
3.2.2 The FCC Solution (Non-Bipartite Topology)
The vacuum is Face-Centered Cubic (FCC), defined by 12 nearest neighbors. The basis
vectors are the permutations of n =
a
√
2
(±1, ±1, 0). The discrete Dirac operator sums
over these 12 non-orthogonal links:
D
F CC
(k) =
12
X
j=1
Γ
j
e
ik·n
j
(7)
where Γ
j
are the spin projections along the bond directions.
Crucially, the FCC lattice is formed of Tetrahedra (Simplexes). A tetrahedron contains
triangular faces (odd cycles of length 3).
• Topological Frustration: It is mathematically impossible to two-color a triangle.
Therefore, the lattice is strictly Non-Bipartite.
• Breaking the Symmetry: The shift k → k + π is no longer a symmetry of the
Hamiltonian because the odd loops introduce phase factors that do not cancel.
3.2.3 Dispersion Analysis at the Zone Boundary
We analyze the dispersion relation E(k) at the boundary of the FCC Brillouin Zone (a
Truncated Octahedron). The energy squared is proportional to:
E(k)
2
∝
X
j
(1 −cos(k · n
j
)) + Cross Terms (8)
Unlike the cubic case where terms decouple, the FCC geometry intrinsically mixes mo-
mentum components.
• At Origin (k = 0): E = 0 (The physical massless chiral fermion).
• At the L-Point (Zone Face Center): The specific geometry of the truncated oc-
tahedron ensures that the dispersion is maximally massive. The spurious “doubler”
mode is lifted entirely to the UV cutoff scale (E ≈ 1/a).
Thus, the mirror fermion is naturally decoupled from the low-energy spectrum. The
K = 12 lattice naturally supports a single species of Chiral Fermion without requiring
artificial Wilson mass terms [3].
4 Review: The Gauge Sector and Plaquettes
Forces arise from the requirement to preserve lattice integrity under local rotation (Gauge
Invariance). If we rotate a node ψ
n
by a phase e
iθ
, the connecting stitch to neighbor m is
strained. To preserve the stitch energy, we introduce a Compensating Link Variable U
nm
(the Gauge Field).
5
4.1 Yang-Mills Action from Wilson Loops
As originally formulated by Wilson [5], the energy of the gauge field is the cost of “curving”
the lattice. We measure this curvature by transporting a vector around a closed lattice
loop (plaquette). The operator is the Wilson Loop:
U
loop
= P exp
ig
I
A
µ
dx
µ
(9)
Using the Baker-Campbell-Hausdorff formula for small lattice spacing a, the deviation
from unity is:
U
loop
≈ 1 + iga
2
F
µν
−
1
2
g
2
a
4
F
2
µν
(10)
The real part of the energy density is therefore proportional to F
2
µν
. This yields the
standard Yang-Mills Lagrangian:
L
Gauge
= −
1
4
F
µν
F
µν
(11)
This confirms that Electromagnetism and the Strong Force are simply the Elastic Potential
Energy of the discrete vacuum resisting topological curvature.
4.2 The Geometric Origin of the Plaquette
In the Cuboctahedral vacuum (K = 12), these Wilson Loops are not abstract math-
ematical constructs; they are physically realized as the Square Bases of the stabilizing
continuous pyramids.
• Matter (Fermions): Generated by the Tetrahedral (triangular) sectors via non-
bipartite topological twisting.
• Forces (Bosons): Generated by the Square Pyramidal (bipartite) sectors via con-
tinuous flux screening.
5 The Higgs: Freezing the Lattice
In the Standard Model, the Higgs potential V (Φ) = −µ
2
|Φ|
2
+ λ|Φ|
4
is assumed to break
electroweak symmetry. We derive the self-coupling parameter λ directly from discrete
combinatorics.
5.1 The Geometric Origin of λ
We interpret λ as the geometric stiffness of the vacuum voxel. It is defined as the ratio
of Topological Boundary Constraints to Volumetric Configurations in the saturated unit
cell.
• Configuration Space (N
V
= 1728): The configuration space volume of a discrete
dynamical system is defined by Ω = (States per Site)
Dimension
. Since the local con-
nectivity is K = 12 and the vacuum manifold has Hausdorff dimension D = 3, the
total macroscopic phase space volume per saturated defect is Ω = 12
3
= 1728 [2].
6
• Holographic Boundary Constraint (N
S
= 108): The boundary constraint is
defined by the interaction surface of the stable proton defect (the Trefoil knot, c = 3).
The absolute phase space of an undisturbed 2D holographic lattice membrane is
K
2
= 144. However, the continuous 1D core of the knot punctures this membrane
at its c = 3 topological crossings. Each crossing acts as a 0D constraint that
rigidly locks one full local coordination shell (K = 12), subtracting exactly c ×K =
36 degrees of freedom. The effective locking constraint is therefore the residual
boundary surface tension: N
S
= K
2
− cK = 144 − 36 = 108.
• Derivation of the Factor 2 (Bi-Directional Flux): Crucially, because the
vacuum is a tessellation, every face is a shared interface between two adjacent unit
cells. The topological flux is therefore strictly bi-directional (Cell A ↔ Cell B),
introducing a necessary topological factor of 2 to the effective surface interface.
The geometric bare coupling is precisely:
λ
geo
= 2 ×
N
S
N
V
= 2 ×
108
1728
=
2
16
= 0.125 (12)
5.2 The Mass Prediction and RG Flow
Using the standard continuous relation m
h
= v
√
2λ and the measured vacuum expectation
value v = 246.22 GeV, the bare lattice prediction is:
m
bare
= 246.22 ×
p
2(0.125) ≈ 123.11 GeV (13)
This represents the Bare Coupling at the discrete lattice cutoff (the Planck Scale, M
P
). A
perturbative continuous shift from λ
bare
= 0.125 to λ
obs
≈ 0.129 corresponds to an approx-
imately 3% renormalization correction, which perfectly mirrors the expected magnitude
for radiative corrections running from M
P
down to the Electroweak scale [9].
6 Falsifiability: Gravitational Echoes
To rigorously test the geometric vacuum hypothesis, we look to the extreme regime where
lattice saturation fails: the Black Hole Event Horizon. The concept of post-merger gravita-
tional echoes arising from Planck-scale horizon modifications is an active and contentious
area of research [11]. We ground our prediction in the specific mechanics of lattice melting.
6.1 The Vacuum Phase Transition
Standard General Relativity models the horizon as a coordinate singularity on a smooth
manifold. In a discrete framework, it is a physical phase boundary where the lattice
“melts” from the saturated Crystal Phase (K = 12) to a disordered Fluid Phase.
As the rigid FCC crystal melts into an amorphous fluid, the network loses the highly
coordinated transverse shear planes required to sustain v = c. Relying instead on un-
coordinated longitudinal percolation, wave propagation is heavily retarded. As derived
computationally in our geometric renormalization framework [4], the propagation speed
in this completely melted boundary layer drops to the fundamental longitudinal limit:
v
fluid
≈
1
4
c (14)
7
6.2 The Echo Prediction (2.7 ms)
If the horizon is a physical boundary, infalling gravitational waves will partially scatter.
The “echo” is a reflection from this slow zone. Assuming the melted phase boundary
extends roughly one Schwarzschild radius (R
s
) outward from the classical horizon surface,
the one-way signal delay ∆t through this layer is calculated as the thickness divided by
the local velocity:
∆t =
R
s
c/4
=
4R
s
c
=
4(2GM/c
2
)
c
=
8GM
c
3
(15)
Case Study: GW250114
For the recently detected, exceptionally high-SNR merger (M
f
≈ 68.5M
⊙
) [12]:
∆t ≈
8 ×(6.674 × 10
−11
) ×(68.5 × 1.989 × 10
30
)
(2.998 ×10
8
)
3
≈ 2.7 ms (16)
6.3 Comparison with Observation
Standard continuous analyses of GW250114 typically search for quasi-normal mode over-
tones (such as the “221 mode”) containing significant residual energy in the window
t < 3.0 ms post-merger.
• Standard Model: Interprets this strictly as a coherent damped sinusoid (Ring-
down) governed by a smooth continuous manifold.
• Discrete Lattice: Interprets this as a prompt, incoherent lattice echo (Scattering)
emerging from the boundary transit delay.
The predicted geometric delay time (2.7 ms) falls seamlessly into the observed signal
duration window. We predict that a Spectral Flatness Test on this residual energy will
reveal a flat noise spectrum (Ξ ≈ 1), characteristic of disordered boundary scattering,
rather than a pure coherent continuous overtone.
7 Conclusion
We have demonstrated how the structural forms of the Standard Model Lagrangians natu-
rally emerge from the discrete principles of a saturated K = 12 geometry. By interpreting
the vacuum as a rigid FCC mesh, Scalars manifest as lattice vibrations, Fermions as topo-
logical braids, Forces as curvature corrections, and the Higgs field as the lattice freezing
transition.
The model offers two precise, falsifiable numerical forecasts derived purely from dis-
crete topological limits: 1. Bare Higgs Mass: 123.11 GeV (Verified to within ≈ 1.7% of
the observed renormalized mass). 2. Gravitational Echoes: 2.7 ms delay for GW250114
(Matching the established post-merger overtone window).
These results, combined with the rigorous geometric resolution of the Fermion Dou-
bling Problem via non-bipartite topology, provide strong theoretical evidence that the
continuous parameters of the Standard Model are not arbitrary phenomenological con-
stants, but the inevitable mathematical consequences of discrete boundary geometry.
8
References
[1] Kulkarni, R. “The Selection-Stitch Model (SSM): Space-Time Emergence via Evo-
lutionary Nucleation in a Polycrystalline Tensor Network.” Zenodo (2026). https:
//doi.org/10.5281/zenodo.18138227
[2] Kulkarni, R. “The Geometric Origin of Mass: Holographic Mass of the Proton from
K = 12 Lattice Geometry.” Zenodo (Submitted concurrently, 2026). https://doi.
org/10.5281/zenodo.18253326
[3] Kulkarni, R. “Fermion Chirality from Non-Bipartite Topology: Resolving the Dou-
bling Problem via Lattice Saturation.” Zenodo (2026). https://doi.org/10.5281/
zenodo.18410364
[4] Kulkarni, R. “Geometric Renormalization of the Speed of Light and the Origin of the
Planck Scale in a Saturation-Stitch Vacuum.” Zenodo (2026). https://doi.org/10.
5281/zenodo.18447672
[5] Wilson, K. G. “Confinement of Quarks.” Physical Review D 10.8 (1974): 2445.
[6] Creutz, M. Quarks, Gluons and Lattices. Cambridge University Press (1983).
[7] Rothe, H. J. Lattice Gauge Theories: An Introduction. World Scientific (2005).
[8] Nielsen, H. B., & Ninomiya, M. “Absence of neutrinos on a lattice: (I). Proof by
homotopy theory.” Nuclear Physics B 185.1 (1981): 20-40.
[9] Workman, R. L., et al. (Particle Data Group). “Review of Particle Physics.” PTEP
2022, 083C01 (2022).
[10] Rolfsen, D. Knots and Links. Publish or Perish (1976).
[11] Abedi, J., Dykaar, H., & Afshordi, N. “Echoes from the abyss: Tentative evidence
for Planck-scale structure at black hole horizons.” Physical Review D 96.8 (2017):
082004.
[12] Abac, A. G., et al. (LIGO Scientific Collaboration, Virgo Collaboration, and KA-
GRA Collaboration), “Black Hole Spectroscopy and Tests of General Relativity with
GW250114.” arXiv:2509.08099 [gr-qc] (2025). https://arxiv.org/abs/2509.08099
9
