The Cosserat Vacuum
Emergent Gauge Fields and Chiral Fermions from the Cosserat
Mechanics of a Saturated FCC Vacuum
Raghu Kulkarni
∗
Independent Researcher, Calabasas, CA
February 2026
Abstract
We present a first-principles derivation of the Standard Model scalar, spinor, gauge, and
gravitational kinematic sectors from the discrete mechanics of a Face-Centered Cubic (FCC)
vacuum lattice saturated at the Kepler packing limit (K = 12). By modeling the pre-
geometric vacuum as a Chiral Cosserat Solid—a discrete medium where nodes possess in-
dependent translational (u) and microrotational (ω) degrees of freedom—we extract the
continuum Effective Field Theory (a → 0) without assuming a pre-existing continuous man-
ifold. We rigorously demonstrate four primary results: (i) the discrete translational kinetic
energy natively converges to the isotropic Klein-Gordon Lagrangian via a strict geometric
tensor identity; (ii) the non-bipartite triangular topology dynamically generates a geometric
Wilson mass that lifts all fermion doublers, yielding a unique massless chiral fermion while
naturally partitioning the Brillouin zone into a 1 ⊕ 3 flavor-suggestive topology; (iii) the
structural face topologies of the unit cell natively break the A
3
lattice root system down
to the Standard Model SU (3) × SU(2) × U(1) gauge group, analytically fixing the bare
Weak Mixing Angle at sin
2
θ
W
= 3/13 ≈ 0.2308 (experiment: 0.2312); and (iv) the Cosserat
elastic displacement strictly encodes metric deformation, with the discrete Regge calculus
of the FCC simplicial complex converging exactly to the Einstein-Hilbert action.
∗
Corresp ondence: raghu@idrive.com
1
1 Introduction
A persistent and foundational challenge in modern theoretical physics is the tension between the
discrete nature of quantum mechanics and the continuous manifold assumed by both General
Relativity and the Standard Model of particle physics. While continuous formulations—ranging
from standard Quantum Field Theory (QFT) to Grand Unified Theories (GUTs) [1] and String
Theory—have achieved remarkable phenomenological success, they universally require the ad
hoc insertion of arbitrary free parameters. The number of fermion generations, the absolute
scale of gauge couplings, and the exact matrix symmetries must be tuned to match empirical
observation [2], precisely because continuous manifolds impose few rigid structural constraints.
Consequently, significant theoretical effort has been directed toward discrete spacetime mod-
els. Causal Set Theory [3] successfully reproduces a discrete Lorentzian signature but struggles
to natively recover the smooth local metric and Standard Model particle content. Loop Quan-
tum Gravity (LQG) [4] robustly quantizes continuous space into discrete spin networks, yet
coupling chiral fermions to these networks without generating anomalies remains an open prob-
lem. Standard Lattice Gauge Theory [5] excels in computing non-perturbative QCD effects,
but its reliance on simple hypercubic grids famously generates spurious, non-physical fermion
doublers—a phenomenon strictly codified by the Nielsen-Ninomiya No-Go Theorem [7].
In this paper, we approach the discretization problem not from abstract quantum geom-
etry, but from the rigorous classical mechanics of saturated elastic media. The fundamental
mathematical constraint on any regular 3D spatial network is geometric saturation. The Kepler
Conjecture, rigorously proven by Hales [8], dictates that the maximum theoretical density for
packing identical spheres (or uniform information nodes) in three dimensions is achieved by
the Face-Centered Cubic (FCC) lattice, possessing a strict coordination number of K = 12.
The physical selection mechanism for this specific lattice geometry is developed in a companion
work [9], where we show that a pre-geometric network undergoing thermodynamic crystal-
lization under entropy maximization naturally condenses into the K = 12 saturated config-
uration. The coherence properties of this lattice—including the primary acoustic harmonic
λ = ⌈K
√
2⌉ = 17—have been independently verified through direct phonon dispersion simula-
tion [10]. In the present work, we take the saturated FCC geometry as a given structural input
and focus exclusively on extracting its continuum field-theoretic content.
To map this geometry to physical fields, we rely on the framework of Cosserat (or mi-
cropolar) continuum mechanics [11]. Unlike classical Cauchy elasticity, which treats points as
structureless masses, a Cosserat solid assigns both translational (u) and microrotational (ω)
degrees of freedom to every point. As highlighted by Hehl and Obukhov [12], the inclusion of an
independent rotational kinematic field is strictly necessary to support spinor fields in an elastic
spacetime framework.
By treating the K = 12 saturated vacuum as a discrete Chiral Cosserat Solid, we derive
its long-wavelength continuum limit. We do not attempt to construct a complete “Theory of
Everything”; rather, we demonstrate four specific, highly constrained geometric results. First,
we show that the K = 12 geometry naturally yields an isotropic scalar Lagrangian. Second, we
prove that the non-bipartite triangular faces of the FCC lattice dynamically generate a Wilson
mass term [5] that circumvents the fermion doubling problem, yielding a unique massless chiral
spinor. Third, we map the local symmetries of the lattice bonds (the A
3
root system) to the
Standard Model gauge group, demonstrating that the geometric segregation of flux channels
provides a zero-parameter derivation of the bare Weak Mixing Angle. Finally, we demonstrate
that the elastic Cosserat displacement natively encodes spacetime curvature, reproducing the
Einstein-Hilbert action.
2
2 The Lattice Geometry
2.1 The Basis Vectors
We define the vacuum lattice Λ as the set of discrete points generated by the primitive vectors of
the FCC lattice. The physical lattice spacing is a. To preserve clear algebraic integers, we scale
the momentum space vectors such that the 12 nearest-neighbor vectors n
j
are permutations of
(±1, ±1, 0):
n
1
= (1, 1, 0), n
2
= (1, −1, 0), n
3
= (−1, 1, 0), n
4
= (−1, −1, 0),
n
5
= (1, 0, 1), n
6
= (1, 0, −1), n
7
= (−1, 0, 1), n
8
= (−1, 0, −1),
n
9
= (0, 1, 1), n
10
= (0, 1, −1), n
11
= (0, −1, 1), n
12
= (0, −1, −1).
(1)
These 12 vectors precisely define the root system of the Lie algebra A
3
(corresponding to
the group SU(4)) [13], providing the underlying algebraic foundation for the allowed gauge
symmetries.
2.2 Topology: Triangles and Squares
The fundamental coordination envelope of the K = 12 FCC lattice is the Cuboctahedron, which
possesses exactly 14 faces. A critical property of this geometry is that it strictly segregates these
faces based on loop connectivity distance:
• 8 Triangular Faces: Formed exclusively by nearest-neighbor bonds. A minimal closed
loop requires three vectors (e.g., n
1
+ n
7
+ n
12
= 0). In graph theory, odd-cycle loops are
strictly non-bipartite.
• 6 Square Faces: Bounded by the exact same nearest-neighbor bonds as the triangles,
but defining orthogonal planar loops. Even-cycle square loops are strictly bipartite.
This structural dichotomy (non-bipartite vs. bipartite) governs the propagation of internal
phase fields across the lattice boundaries.
3 The Cosserat Action
In a discrete Cosserat lattice [11], the kinematic state of the vacuum is defined by two indepen-
dent geometric fields defined at each node x ∈ Λ:
1. Displacement u
µ
(x): A continuous vector field representing spatial translation.
2. Microrotation ω
µ
(x): A pseudo-vector field representing internal orientation (spin).
The classical Action of this medium is given by:
S =
Z
dt
X
x∈Λ
1
2
m
˙
u
2
+
1
2
I
˙
ω
2
− V(u, ω)
(2)
The discrete potential energy density V is the harmonic sum over the 12 nearest neighbors j:
V =
1
2
12
X
j=1
α(∆
j
u)
2
+ β(∆
j
ω)
2
+ γ(∆
j
u × ω)
(3)
where ∆
j
ϕ = ϕ(x + n
j
) − ϕ(x) is the finite difference operator along the lattice bonds.
3
4 The Scalar Sector: Emergent Isotropy from the Tensor Iden-
tity
We naturally identify the translational displacement magnitude |u| with the fundamental scalar
field ϕ. To extract the long-wavelength continuum Lagrangian, we expand the finite difference
operator in a Taylor series for a small lattice spacing a:
ϕ(x + n
j
) − ϕ(x) ≈ a(n
j
· ∇)ϕ +
a
2
2
(n
j
· ∇)
2
ϕ (4)
Substituting this into the discrete potential term
P
(∆
j
ϕ)
2
yields the spatial kinetic energy:
12
X
j=1
[a(n
µ
j
∂
µ
ϕ)]
2
= a
2
∂
µ
ϕ∂
ν
ϕ
12
X
j=1
n
µ
j
n
ν
j
(5)
We evaluate the geometric tensor sum T
µν
=
P
12
j=1
n
µ
j
n
ν
j
explicitly over the K = 12 basis. For
diagonal elements (µ = ν = x), exactly 8 vectors have non-zero components (±1)
2
= 1, yielding
a sum of 8. For off-diagonal elements (µ = ν), the terms cancel in exact pairs due to the lattice’s
inversion symmetry, yielding 0.
Thus, T
µν
= 8δ
µν
. This geometric tensor identity demonstrates that the discrete kinetic en-
ergy density natively converges to the continuous, directionless isotropic scalar (Klein-Gordon)
Lagrangian:
L
kin
∝ 4a
2
(∂
µ
ϕ)
2
(6)
Despite the discrete lattice possessing only cubic symmetry, the macroscopic continuous scalar
field propagates with full Euclidean isotropy.
5 The Spinor Sector: Doubler Removal
We identify the fundamental spinor field with the Cosserat microrotation ω.
5.1 The Discrete Dirac Operator and the Wilson Mass
Naively discretizing the Dirac equation on a simple hypercubic lattice generates 2
D
massive
copies of the fermion at the Brillouin zone boundaries [7]. In standard Lattice Gauge Theory,
an ad hoc Wilson mass term (1 − cos(k)) is manually added to the action to suppress these
doublers [5]. In the FCC geometry, we find this term is dynamically generated by the local
lattice topology.
We construct the full discrete Dirac operator on the FCC lattice in momentum space:
D(k) =
12
X
j=1
[i(γ · n
j
) sin(k ·n
j
) + I(1 − cos(k ·n
j
))] (7)
The (1 −cos) mass term arises intrinsically from second-order hopping processes along the non-
bipartite triangular loops. On an FCC lattice, the sum of two nearest-neighbor vectors naturally
equates to a third nearest-neighbor vector (e.g., (1, 0, 1) + (−1, 1, 0) = (0, 1, 1)). The effective
second-order hopping amplitude through the intermediate site at x + n
a
is proportional to
t
2
/∆E, which in momentum space yields a term proportional to cos(k ·n
c
) where n
c
= n
a
+ n
b
.
Because this vector sum n
c
is itself a nearest-neighbor vector in the FCC lattice, this geometric
correction enters at the exact same O(1) magnitude as the bare hopping term.
4
5.2 The Uniqueness Theorem
To prove that the discrete FCC vacuum hosts exactly one physical, massless fermion species,
we evaluate the effective mass gap M
eff
(k) =
P
12
j=1
(1 − cos(k · n
j
)) across the Brillouin Zone.
Theorem 1. The effective mass function M (k) on an FCC lattice vanishes if and only if
k = (0, 0, 0) (modulo reciprocal lattice vectors).
Proof. For a strictly massless mode (M(k) = 0), we require cos(k · n
j
) = 1 for all 12 vectors
simultaneously. Consider the specific subset: n
1
= (1, 1, 0) and n
2
= (1, −1, 0).
k ·n
1
= 0 =⇒ k
x
+ k
y
= 0
k ·n
2
= 0 =⇒ k
x
− k
y
= 0
Adding these relations implies k
x
= 0 and k
y
= 0. Substituting this result into n
5
= (1, 0, 1)
yields k
x
+ k
z
= 0, which strictly enforces k
z
= 0. Thus, the origin k = (0, 0, 0) is the unique
massless solution in the fundamental Brillouin zone.
By natively generating the Wilson mass through its triangular face topology, the K = 12
vacuum naturally evades the fermion doubling problem.
5.3 Topological Implications for Generational Flavor
The reciprocal space of the FCC lattice forms a Truncated Octahedron Brillouin Zone. The
kinematic zero-modes of the bare Dirac operator separate into distinct topological classes: the
fundamental mode sits at the zone center (Γ), while the primary doublers are isolated at the
centers of the 8 hexagonal faces (the L-points). Because k and −k are physically identified,
there are exactly 4 independent L-points forming an isolated topological quartet deep in the
bulk.
These 4 L-points rigidly transform under the discrete tetrahedral group S
4
. In representation
theory, this 4-dimensional representation is reducible into a singlet and a triplet:
4 → 1 ⊕ 3 (8)
While a full mapping of these specific modes to physical fermion generations—including precise
mass hierarchies and mixing matrices—remains outside the scope of this present structural
paper, the Brillouin zone topology naturally admits a geometric decomposition that is highly
suggestive of the 3 + 1 generational structure observed in nature (three active generations plus
a sterile state).
6 The Gauge Sector: Weak Mixing Angle
6.1 Root System Symmetry Breaking
The fundamental set of 12 hopping vectors is isomorphic to the root system of SU(4) (the A
3
algebra). The unit cell’s heterogeneous face topology inherently breaks this parent symmetry.
Using the standard Dynkin Diagram decomposition for A
3
(◦ − ◦ − ◦), removing an end node
structurally yields A
2
⊕ A
1
⊕ U (1) [13].
1. A
2
(8 Roots): This subsystem maps identically to the 8 non-bipartite Triangular Faces.
Because these frustrated loops forbid perfect dipole screening, they natively support a
confining phase analogous to the Strong Force SU (3).
2. A
1
⊕U(1) (4 Generators): This subsystem maps to the 6 bipartite Square Faces, which
permit unfrustrated charge alternation and macroscopic flux screening. This structural
symmetry aligns with the Electroweak Force SU(2) × U(1).
5
6.2 Derivation of the Bare Weak Mixing Angle
The Weinberg mixing angle (θ
W
) measures the proportional mixing between the spatial gauge
projection (SU(2)) and the internal scalar phase (U(1)). In standard Lattice Gauge Theory,
the squared coupling constant is inversely proportional to the number of available discrete flux
channels (N): g
2
∝ 1/N .
A fundamental localized geometric particle state within the lattice requires a complete topo-
logical vertex star. In standard Lattice Gauge Theory, the smallest localized, gauge-invariant
structural unit (such as the simplest closed Wilson loop enveloping a volume) must close on
the vertex star—the central node plus its surrounding coordination shell [6]. For the K = 12
FCC lattice, this minimal bounded 3D bulk cluster strictly consists of the 1 central reference
node plus its 12 immediate nearest neighbors. This yields a mathematically rigid topological
reservoir of N
bulk
= 1 + 12 = 13 discrete degrees of freedom. We analytically partition this
holistic geometric reservoir:
• Weak Sector (SU(2)
L
): The weak generators correspond to the fundamental orthogonal
rotational planes of 3D space. The spatial weak coupling g is therefore distributed over
exactly N
spatial
= 3 channels.
• Hypercharge Sector (U(1)
Y
): The scalar hypercharge must couple to the remaining,
non-spatially oriented internal degrees of freedom within the 13-node bulk: N
hypercharge
=
13 − 3 = 10.
Substituting these geometric flux capacities into the standard gauge mixing formula yields:
sin
2
θ
W
=
g
′2
g
2
+ g
′2
=
1/10
1/3 + 1/10
=
3
13
≈ 0.2308 (9)
This parameter-free topological limit precisely aligns with the empirical Z-pole measurement
(≈ 0.2312) [2].
7 The Gravitational Sector
7.1 Metric Deformation and Analogue Gravity
The Cosserat action defined in Equations (2) and (3) intrinsically contains the gravitational de-
grees of freedom. In the long-wavelength continuum limit, the displacement field u
µ
(x) encodes
the geometric deformation of the effective metric, where macroscopic strain directly corresponds
to spacetime curvature. This formal identification of elastic strain with the gravitational field
is a well-established equivalence in analogue gravity models [14].
7.2 Regge Calculus on the FCC Simplicial Complex
The precise geometric mapping to General Relativity is achieved through the native simplicial
structure of the lattice. The Delaunay triangulation of the Face-Centered Cubic lattice naturally
decomposes the vacuum into a close-packed simplicial complex of tetrahedra and octahedra.
Applying Regge calculus [15] to this simplicial complex, the spatial curvature is fundamentally
concentrated at the discrete structural hinges (edges) h with length L
h
and angular deficit δ
h
.
The discrete Regge action over these hinges takes the form:
S
Regge
=
1
8πG
X
h
L
h
δ
h
(10)
Crucially, as rigorously proven by Cheeger, M¨uller, and Schrader [16], this discrete simplicial
action strictly converges to the continuum Einstein-Hilbert action
R
d
4
x
√
−g
R
16πG
in the a →
0 limit.
6
7.3 Gravitational Coupling and Lorentz Invariance
In this saturated elastic framework, the phenomenological gravitational coupling is strictly
determined by the discrete mechanics. The elastic moduli α and β from Equation (3) define the
bulk modulus of the lattice. Newton’s gravitational constant G emerges inversely proportional
to the stiffness of the vacuum, governed by the relation G ∼ 1/(ρv
2
a
2
), where ρ is the discrete
node density and v is the acoustic propagation speed. A precise determination of G from the
lattice parameters requires specifying the absolute energy scale of the elastic moduli, which lies
beyond the scope of the present kinematic analysis.
Furthermore, the exact O
h
point group symmetry of the FCC lattice ensures that the emer-
gent gravitational sector and macroscopic wave propagation remain isotropic to order O(a
2
/L
2
),
where L is the macroscopic observation scale. Because the lattice spacing a is fixed near the
Planck length, this geometric isotropy safely satisfies the stringent bounds on Lorentz invariance
violation established by modern high-energy astrophysical observations, such as those from the
Fermi-LAT gamma-ray telescope [17].
8 Conclusion
We have presented a rigorous, mathematical extraction of Standard Model kinematics from the
classical elasticity of a saturated K = 12 discrete lattice. By relying strictly on the topology
of the Face-Centered Cubic structure, we demonstrated that an isotropic scalar field natively
emerges via exact geometric tensor identities. Furthermore, we proved that the non-bipartite
triangular faces of the unit cell dynamically generate a Wilson mass that resolves the fermion
doubling problem, yielding a unique massless chiral fermion and a Brillouin zone topology highly
suggestive of a three-generation flavor structure. We then mapped the A
3
root system to the
lattice face geometries, analytically partitioned the SU(3) ×SU(2) ×U(1) gauge fields, and de-
rived the bare Weak Mixing Angle as sin
2
θ
W
= 3/13. Finally, we demonstrated that the elastic
Cosserat displacement natively encodes spacetime curvature, reproducing the Einstein-Hilbert
action through the rigorous convergence of the FCC Regge complex. These mathematical re-
sults strongly suggest that the core phenomenological parameters of continuous field theory and
gravity are the exact structural footprints of a discrete, maximally saturated quantum vacuum.
References
[1] Georgi, H., and Glashow, S. L. “Unity of All Elementary-Particle Forces,” Phys. Rev. Lett.
32, 438 (1974).
[2] Workman, R. L. et al. (Particle Data Group). “Review of Particle Physics,” Prog. Theor.
Exp. Phys. 2022, 083C01 (2022).
[3] Bombelli, L., Lee, J., Meyer, D., and Sorkin, R. D. “Space-time as a causal set,” Phys.
Rev. Lett. 59, 521 (1987).
[4] Rovelli, C. Quantum Gravity, Cambridge University Press (2004).
[5] Wilson, K. G. “Confinement of quarks,” Phys. Rev. D 10, 2445 (1974).
[6] Rothe, H. J. Lattice Gauge Theories: An Introduction, World Scientific (2012).
[7] Nielsen, H. B., and Ninomiya, M. “No-go theorem for regularizing chiral fermions,” Nucl.
Phys. B 185, 20 (1981).
[8] Hales, T. C. “A proof of the Kepler conjecture,” Annals of Mathematics 162, 1065-1185
(2005).
7
[9] Kulkarni, R., “Geometric Phase Transitions in a Discrete Vacuum: Deriving Cosmic Flat-
ness, Inflation, and Reheating from Tensor Network Topology.” Submitted to Physics Let-
ters B, 2026. https://doi.org/10.5281/zenodo.18727238
[10] Kulkarni, R., “Constructive Verification of K = 12 Lattice Saturation: Exploring Kinematic
Consistency in the Selection-Stitch Model.” Submitted to Physics Letters B, 2026. https:
//doi.org/10.5281/zenodo.18294925
[11] Cosserat, E., and Cosserat, F. Theorie des Corps Deformables, Hermann et Fils, Paris
(1909).
[12] Hehl, F. W., and Obukhov, Y. N. Foundations of Classical Electrodynamics: Charge, Flux,
and Metric, Birkh¨auser (2003).
[13] Fulton, W., and Harris, J. Representation Theory: A First Course, Springer (1991).
[14] Barcel´o, C., Liberati, S., and Visser, M. “Analogue Gravity,” Living Reviews in Relativity
14, 3 (2011).
[15] Regge, T. “General relativity without coordinates,” Il Nuovo Cimento 19.3, 558-571 (1961).
[16] Cheeger, J., M¨uller, W., and Schrader, R. “On the curvature of piecewise flat spaces,”
Communications in Mathematical Physics 92, 405-454 (1984).
[17] Abdo, A. A., et al. (Fermi LAT Collaboration). “A limit on the variation of the speed of
light arising from quantum gravity effects,” Nature 462, 331-334 (2009).
8
