Defect-Bound Modes of the Naive Dirac Operator on the FCC
Lattice:
Spectral Multiplet Structure at Tetrahedral and Octahedral Voids
Raghu Kulkarni
∗
SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA
May 2026
Abstract
The face-centered cubic (FCC) lattice has two classes of interstitial void: tetrahedral,
with coordination 4, and octahedral, with coordination 6. A trapped extra node at either
void is a localized point defect coupled to its surrounding FCC sites through nearest-neighbor
bonds. We compute the bound-mode multiplet of the naive Dirac operator at each defect
class and decompose it under the residual point-group symmetry. At a tetrahedral defect
the multiplet has dimension 4 and decomposes under T
d
∼
=
S
4
as A
1
⊕T
1
. At an octahedral
defect the multiplet has dimension 6 and decomposes under O
h
as A
1g
⊕ E
g
⊕ T
1u
. The
3-dimensional T
1
piece of the tetrahedral multiplet has the dimension and group-theoretic
properties needed to carry an internal SU(3) color triplet; the matter paper of Kulkarni
[Phys. Open 27, 100423 (2026)] builds such an action on the three perfect matchings of K
4
,
the bonded graph of the tetrahedral defect’s surrounding sites. The 3-dimensional T
1u
piece
of the octahedral multiplet is the spatial-vector representation of O
h
, not an internal triplet.
The bonded graph of the octahedron’s six surrounding sites is the complete tripartite K
2,2,2
rather than K
6
, because the three antipodal vertex pairs sit at second-nearest-neighbor dis-
tance in the FCC lattice and are not physical bonds; K
2,2,2
has 8 perfect matchings and 30
skew-edge pairs, neither giving the 3-color structure needed for the matter-paper construc-
tion. The same obstruction is reached on independent combinatorial grounds in a companion
dark-matter paper [15], which identifies octahedral defects with a dark-sector matter species;
the spectral analysis here gives independent representation-theoretic evidence for that iden-
tification. Tetrahedral defects admit the matter-paper color algebra; octahedral defects do
not. As supporting background we give a complete analytic classification of the free-field
FCC zero-mode spectrum from an exact factorization of the kinetic vector field, with three
topological classes (a singlet at Γ, a quartet at L-points, and three boundary nodal loops)
summing to zero index by Nielsen–Ninomiya. The defect mechanism is unlike the minimally
doubled and supersymmetric-lattice constructions, which act on the bulk Dirac operator;
the analysis here leaves the bulk operator unchanged. The local result depends on the
cuboctahedral K = 12 first-shell structure of FCC, which is not present in HCP (anticuboc-
tahedral D
3h
, with a preferred c-axis at rank 4 in the bulk), BCC (K = 8 with broken K
4
at tetrahedral voids), or hypercubic lattices. Among the two 3D close-packed lattices only
FCC carries the full O
h
first-shell symmetry needed for an isotropic substrate.
1 Introduction
The naive Dirac operator on a hypercubic lattice in D dimensions has 2
D
geometrically equiv-
alent doublers [1, 2]. Wilson fermions [3] give them a mass term that breaks chiral symmetry.
∗
raghu@idrive.com
1
Staggered fermions [4, 5] cut the count to 2
⌈D/2⌉
tastes at the cost of a residual flavor symmetry.
Minimally doubled formulations [6, 7, 8, 9, 10] keep only two species in the continuum limit by
breaking hypercubic symmetry. Supersymmetric lattice constructions [11, 12] use root lattices
such as A
∗
4
to preserve part of the supersymmetry algebra exactly. Each of these acts on the
bulk operator.
We take a different starting point. Consider the standard naive Dirac operator on the face-
centered cubic (FCC) lattice, and look at the spectrum it produces when a localized point defect
is inserted at an interstitial void. FCC has two void classes per conventional cubic cell: eight
tetrahedral voids (coordination 4) and four octahedral voids (coordination 6). A trapped extra
node at either is a natural physical defect. The bound modes at the defect form a multiplet
under the residual point-group symmetry, and the structure of that multiplet decides whether
the defect can carry an internal symmetry algebra of a given type.
The setting is the matter paper [14] of the SSMTheory program, in which physical matter
is identified with localized lattice defects rather than with a bulk fermion field. There, color
quantum numbers are assigned by combinatorial counting: the complete graph K
4
on the
four valence bonds of a tetrahedral defect has exactly three perfect matchings, giving a three-
dimensional internal color Hilbert space on which an su(3) algebra acts. That construction
is group-theoretic and uses the local bond connectivity only, with no input from the Dirac
spectrum. The question this paper takes up is whether the spectral content of the Dirac operator
at the same defect is consistent with the assigned color structure — that is, whether the bound-
mode multiplet contains a representation of the right dimension and the right transformation
properties — and what the corresponding analysis gives for the octahedral defect.
Results. At a tetrahedral defect the bound-mode multiplet has dimension 4 and decomposes
under the tetrahedral point group T
d
∼
=
S
4
as
4
tet
= A
1
⊕ T
1
, dim 1 + 3. (1)
The 3-dimensional T
1
piece has the dimension required to carry an internal SU(3) color triplet.
The matter paper’s K
4
perfect-matching construction defines such an action on it (Section 4).
At an octahedral defect the multiplet has dimension 6 and decomposes under the octahedral
point group O
h
as
6
oct
= A
1g
⊕ E
g
⊕ T
1u
, dim 1 + 2 + 3. (2)
The 3-dimensional T
1u
piece is the vector representation of O
h
, not an internal symmetry.
The bonded graph of the six surrounding sites is the complete tripartite graph K
2,2,2
, not K
6
:
the three antipodal pairs of the octahedron sit at second-nearest-neighbor distance in the FCC
lattice, not nearest-neighbor, so they are not physical bonds. K
2,2,2
has eight perfect matchings,
not the three required for the matter-paper algebra. No internal SU(3) structure of that form
is supported on the octahedral bound-mode multiplet (Section 5). The companion dark-matter
paper [15] reaches the same conclusion combinatorially and identifies octahedral defects with a
dark-sector matter species. The spectral analysis here is independent representation-theoretic
evidence for that identification. Tetrahedral defects support the matter-paper color algebra;
octahedral defects do not.
A note on what this result is and is not. The Dirac analysis does not derive SU(3) color
from the spectrum; the matter-paper construction is logically prior and uses different input (K
4
perfect matchings, not Dirac eigenmodes). What we add is a spectral consistency check at each
defect class, which the tetrahedral defect passes and the octahedral defect fails. Failure of the
spectral check at the octahedral defect, together with the combinatorial K
2,2,2
obstruction, is
the substantive content distinguishing the two FCC defect classes as candidate matter species.
2
Bulk background. Sections 2 and 3 give the free-field bulk spectrum on FCC. The ki-
netic vector field factors exactly (Theorem 1), so the full zero-mode structure can be classi-
fied by cases. Three topological classes appear: a singlet at Γ, a quartet at the four L-points
(±π/2, ±π/2, ±π/2), and three boundary nodal loops along the X–W segments of the truncated-
octahedron BZ. The chirality sum is +1 − 4 + 3 = 0, as required by Nielsen–Ninomiya. The
three classes carry distinct Wilson masses (0, 12/a, 16/a) and are not equivalent to the eight
uniformly-massed doublers of the naive simple-cubic operator. This bulk analysis is supporting
context for the defect-bound-mode analysis in Sections 4 and 5, which is the main new content.
Relation to other reduced-fermion constructions. The defect mechanism we use is un-
like minimally doubled fermions [6, 7, 8], supersymmetric lattices [11, 12], twisted-mass, and
domain-wall constructions: those modify the bulk operator. Here the bulk operator is un-
changed; the physical content sits in localized topological excitations whose internal symmetries
are set by the defect’s local geometry. Section 6 discusses the comparison in detail and works
out where the cuboctahedral first-shell structure of FCC fits among the alternatives. The clos-
est 3D candidate is HCP, with the same coordination number K = 12, but HCP’s first shell
is the anticuboctahedron and gives a preferred-axis anisotropy at the rank-4 bond-tensor level;
FCC is the only 3D close-packed lattice with a fully cubic-symmetric rank-4 bond tensor and
the corresponding full T
d
/ O
h
symmetry at its interstitial voids.
2 The FCC Lattice and the Naive Dirac Operator
The FCC lattice has twelve nearest-neighbor vectors (in units of a/2):
N = {(±1, ±1, 0), (±1, 0, ±1), (0, ±1, ±1)}. (3)
The convex hull of these twelve points is the cuboctahedron [17], with symmetry group O
h
of
order 48. The reciprocal lattice is body-centered cubic; the first Brillouin zone is a truncated
octahedron [18] with high-symmetry points
Γ = (0, 0, 0), L = (±
π
2
, ±
π
2
, ±
π
2
), X = (π, 0, 0), W = (π,
π
2
, 0) (4)
(and cubic permutations of X, W ).
The naive Dirac operator in three spatial dimensions uses two-component spinors with no
γ
5
. Chiralities below are local degrees of the kinetic map f : T
3
→ R
3
, not eigenvalues of
γ
5
[16]. The kinetic vector field is
f
µ
(k) =
1
a
X
n∈N
n
µ
sin(k · n), µ = 1, 2, 3, (5)
and the Dirac operator is D(k) = iσ ·
f(k). Zero modes occur wherever |
f(k)| = 0.
Theorem 1 (Factorization). The kinetic vector field of Eq. (5) factors as
f
µ
(k) =
4
a
sin(k
µ
)
cos(k
ν
) + cos(k
ρ
)
, (6)
where (µ, ν, ρ) is any permutation of (1, 2, 3).
Proof. The twelve FCC neighbors partition into three sheets S
µν
containing the four vectors
with zero ρ-component. Only sheets S
µν
and S
µρ
contribute to f
µ
(sheet S
νρ
has n
µ
= 0). From
S
µν
, the four vectors (n
µ
, n
ν
) ∈ {(±1, ±1)} contribute
sin(k
µ
+ k
ν
) + sin(k
µ
− k
ν
) −sin(−k
µ
+ k
ν
) −sin(−k
µ
− k
ν
) = 4 sin(k
µ
) cos(k
ν
),
using sin A + sin B = 2 sin
A+B
2
cos
A−B
2
. Sheet S
µρ
contributes 4 sin(k
µ
) cos(k
ρ
) identically.
Summing gives Eq. (6).
3
3 Bulk Zero Modes and Their Topological Indices
Theorem 2 (Complete classification). The zero set f
−1
(0) on the FCC Brillouin zone consists
of exactly:
1. Γ = (0, 0, 0): one isolated zero.
2. L = (±π/2, ±π/2, ±π/2): four isolated zeros after antipodal identification.
3. Nodal lines on the BZ boundary connecting X- to W -points.
No other zeros exist.
Proof. By Theorem 1, f
µ
= 0 requires sin(k
µ
) = 0 (call this S
µ
) or cos(k
ν
) + cos(k
ρ
) = 0 (C
µ
).
Exhausting all 2
3
= 8 sign combinations:
Case SSS (S
1
∧S
2
∧S
3
): each k
µ
∈ {0, π}. BZ constraints admit Γ and the three X-points.
Case CCC: subtracting pairs gives cos k
1
= cos k
2
= cos k
3
; substituting back gives
2 cos k
1
= 0, so k
µ
= ±π/2. These are the four L-points (inside the BZ).
Cases SSC, SCS, CSS: two momenta in {0, π} with their cosines summing to zero; the third
is free. These produce nodal lines: (0, π, t), (π, 0, t), and permutations — the X–W boundary
segments.
Cases SCC, CSC, CCS: require sin k
µ
= 0 and cos k
µ
= −cos k
ν
= −cos k
ρ
. If k
µ
= 0:
k
ν
= k
ρ
= π, but |k
ν
| + |k
ρ
| = 2π > π (outside BZ). If k
µ
= π: k
ν
= k
ρ
= 0 (X-point, already
in case SSS).
The eight cases are exhaustive.
3.1 Topological indices and chiralities
The topological index at an isolated zero k
i
is χ
i
= sgn(det J), where J
µν
(k) = ∂f
µ
/∂k
ν
=
1
a
P
n
n
µ
n
ν
cos(k · n).
Proposition 1 (Γ-point). At Γ, J
µν
= (8/a)δ
µν
, det J = 512/a
3
> 0, χ
Γ
= +1.
Proof. At k = 0, cos(k ·n) = 1 for all n. The sum
P
n
n
µ
n
ν
= 8δ
µν
: each diagonal entry receives
a contribution of +1 from the eight FCC vectors with |n
µ
| = 1, and off-diagonal entries cancel
by inversion symmetry.
Proposition 2 (L-points). At L = (π/2, π/2, π/2), J(L) = −(4/a)(E − I), where E is the
all-ones matrix and I the identity. Eigenvalues are {+4, +4, −8}/a, det J = −128/a
3
< 0,
χ
L
= −1.
Proof. At L, k ·n ∈ {0, ±π}. Six neighbors have cos = +1; six have cos = −1. Diagonal entries
P
n
n
2
µ
cos(k ·n) = 0 (equal numbers of +1 and −1 along each axis). For the off-diagonal entries
the four neighbors with n
µ
n
ν
= 0 all give cos(k · n) = −1 and n
µ
n
ν
= ±1, contributing −4/a.
So J(L) = −(4/a)(E − I), eigenvalues {+4/a, +4/a, −8/a}.
With four independent L-points, χ
total
L
= 4 ×(−1) = −4.
3.2 Boundary chiralities via winding number
Along a nodal line in k
2
(for example k = (π, t, 0)), the transverse field (f
1
, f
3
) winds around
the origin. We computed the winding number on a small circle of radius ϵ = 0.01 in the (k
1
, k
3
)
plane, sampled with 2000 angular points, obtaining w = +1.000±0.003 independent of position
t along the line. Topological stability follows from the factorization in Eq. (6): at k
1
= π,
sin(k
1
) = 0 forces f
1
≡ 0, so any zero of |f | on this face must satisfy f
2
= f
3
= 0, confining
the locus to the nodal line. The three independent boundary loops (one per spatial axis) each
contribute χ = +1, so χ
boundary
= +3.
4
The Nielsen–Ninomiya sum is
χ
total
= +1
|{z}
Γ
+ −4
|{z}
L
+ +3
|{z}
boundary
= 0. ✓ (7)
Figure 1 shows the three classes on the truncated-octahedron Brillouin zone.
0
k
1
[ /
a
]
0
k
2
[ /
a
]
0
k
3
[ /
a
]
( = +1)
L
( = 1)
X
--
W
nodal loop ( = +1)
Figure 1: FCC first Brillouin zone (truncated octahedron) with the three zero-mode classes of
the naive Dirac operator. Green circle: Γ-point singlet, χ = +1. Red squares: four L-points at
(±π/2, ±π/2, ±π/2), χ = −1 each, giving χ
total
L
= −4. Blue outlines: the six square boundary
faces, paired by inversion into three independent X–W nodal loops with χ = +1 per loop. The
Nielsen–Ninomiya sum +1 − 4 + 3 = 0 is realized by these three classes.
3.3 Wilson masses
A standard isotropic Wilson term gives the bulk zero-mode classes mass M
W
= 0 at Γ, 12/a
at each L-point, and 16/a on the boundary loops. The three classes are therefore stratified by
Wilson mass at distinct scales [3], in contrast with the uniform 6/a produced by the standard
Wilson term on the simple cubic lattice. An anisotropic Wilson term with enhanced coupling
along (1, 1, 1) splits the L-quartet under S
4
⊃ S
3
as 1 ⊕ 3, with ∆M = 4/a.
4 Tetrahedral Defects
This section gives the main new content of the paper. The bulk classification of Theorem 2
sets the stage; we now turn to the localized bound-mode multiplet at a tetrahedral interstitial
defect (Figure 2(a)).
4.1 Defect geometry
A tetrahedral void of the FCC lattice is bounded by four FCC sites that form a regular tetra-
hedron. An interstitial node placed at the void centroid sits at equal distance from each of
5
these four sites. Label the surrounding sites {1, 2, 3, 4} and the bonds from the defect to them
{b
1
, b
2
, b
3
, b
4
}. The defect site and the four bonds form a star graph whose point group is
T
d
∼
=
S
4
, the full tetrahedral group of order 24.
The naive Dirac operator restricted to this star has matrix elements between the defect
spinor and each of the four surrounding-site spinors. The Hilbert space localized on the star
factors as C
2
⊗ (C ⊕ C
4
): the spinor index times (defect site ⊕ four surrounding sites). What
sets the bound-mode multiplet structure is the four-dimensional permutation representation of
T
d
on the surrounding sites, with the spinor index decorating each multiplet member.
4.2 Multiplet structure under T
d
The permutation representation of S
4
on four sites has character χ = (4, 2, 0, 1, 0) on the
conjugacy classes (e, (12), (12)(34), (123), (1234)). Projecting onto S
4
irreps gives the following
decomposition.
Proposition 3. The permutation representation of T
d
∼
=
S
4
on the four surrounding sites of a
tetrahedral defect decomposes as
4
tet
= A
1
⊕ T
1
, (8)
with dimensions 1 + 3 = 4.
Proof. Direct projection against the S
4
character table. For the trivial representation A
1
with
character (1, 1, 1, 1, 1) the multiplicity is
⟨A
1
, χ⟩ =
1
24
1 ·1 · 4 + 6 · 1 ·2 + 3 ·1 · 0 + 8 · 1 · 1 + 6 · 1 · 0
= 1.
For the standard rep T
1
with character (3, 1, −1, 0, −1) the multiplicity is
⟨T
1
, χ⟩ =
1
24
1 ·3 · 4 + 6 · 1 ·2 + 3 ·(−1) · 0 + 8 · 0 · 1 + 6 · (−1) · 0
= 1.
The other irreps A
2
, E, T
2
have multiplicity 0.
The A
1
singlet is the symmetric combination |s⟩ = (|1⟩+ |2⟩+ |3⟩+ |4⟩)/2. The T
1
triplet is
spanned by three orthogonal antisymmetric combinations. With spinor decoration each of the
four multiplet members carries a two-component spinor index, giving eight bound-mode states
in total.
4.3 Compatibility with the K
4
color construction
The matter paper [14] builds a three-color quark Hilbert space at a tetrahedral defect from the
perfect-matching combinatorics of K
4
, the complete graph on the four valence bonds. K
4
has
six edges; they partition into three perfect matchings (pairs of disjoint edges covering all four
vertices):
M
1
= {b
12
, b
34
}, M
2
= {b
13
, b
24
}, M
3
= {b
14
, b
23
}, (9)
where b
ij
joins surrounding sites i and j. The Weyl group of the resulting su(3) algebra is
S
3
= S
4
/V
4
, where V
4
= {e, (12)(34), (13)(24), (14)(23)} is the Klein four-group; V
4
acts trivially
on the three matchings, so the quotient acts faithfully.
The matter-paper color Hilbert space H
C
= span{|M
1
⟩, |M
2
⟩, |M
3
⟩} is logically independent
of the Dirac bound-mode space, since it lives on matchings rather than on sites. What needs to
be checked is whether the bound-mode space of Proposition 3 can carry an SU(3) action that
respects the assigned color structure.
Proposition 4. The three-dimensional T
1
sub-multiplet of the tetrahedral bound-mode space
has the dimension required to carry an SU(3) fundamental representation. The matter paper’s
K
4
-derived su(3) algebra can be defined on this sub-multiplet by identifying T
1
basis vectors
with linear combinations of matching states {|M
1
⟩, |M
2
⟩, |M
3
⟩} that respect the S
3
= S
4
/V
4
Weyl-group action.
6
Dimensions match exactly: dim T
1
= 3 = dim H
C
. The group action matches up to a
linear identification of bases. The standard T
1
of S
4
restricted to S
3
⊂ S
4
decomposes as
A ⊕E (dimensions 1 + 2), and the matching representation A
1
⊕ E of S
3
on H
C
has the same
decomposition. The two 3-dimensional spaces are therefore isomorphic as S
3
-modules. The
color representation of the matter paper [14] therefore acts on the bound-mode T
1
via this
isomorphism.
4.4 Numerical confirmation
The 1 + 3 site-permutation decomposition is confirmed by direct diagonalization of the naive
Dirac operator on a finite cluster with a defect. The verification script (verify_defects.py,
included with the submission) builds the 55-site FCC cluster (x + y + z even, |r|
2
≤ 9.5) and
models the trapped-node defect at the tetrahedral void (1/2, 1/2, 1/2) as a strong attractive on-
site potential at the four surrounding FCC sites — the Schur-complement reduction of the extra-
site model after integrating out the defect site. The resulting Hermitian 110 ×110 Hamiltonian
is diagonalized exactly. The eight most-localized eigenstates (= 1 + 3 spatial × 2 spinor) carry
the A
1
and T
1
site representations with total weights converging to 2 and 6 respectively as
the potential strength increases, with ratio T
1
/A
1
→ 3 — the expected ratio for the A
1
⊕ T
1
decomposition.
5 Octahedral Defects
5.1 Defect geometry
An octahedral void of the FCC lattice is bounded by six FCC sites at the vertices of a regular
octahedron. The surrounding sites and the defect site form a local star graph with point group
O
h
of order 48. The defect couples to each of the six surrounding sites through one bond.
5.2 Multiplet structure under O
h
The permutation representation of O
h
on the six octahedron vertices decomposes into O
h
irreps
as follows.
Proposition 5. The permutation representation of O
h
on the six surrounding sites of an octa-
hedral defect decomposes as
6
oct
= A
1g
⊕ E
g
⊕ T
1u
, (10)
with dimensions 1 + 2 + 3 = 6.
This is the standard O
h
decomposition on six octahedron vertices, derivable from the O
h
character table [18, 17]. The A
1g
singlet is the symmetric sum. The E
g
doublet captures two
diagonal differences between antipodal-pair amplitudes. The T
1u
triplet is the vector represen-
tation, with basis vectors transforming as (x, y, z) under O
h
.
5.3 Why T
1u
is not an internal SU (3) triplet
The 3-dimensional T
1u
has the same dimension as an SU(3) fundamental, but its transformation
properties identify it as a spatial-vector representation of O
h
, not an internal-symmetry triplet.
Its defining property is that it transforms under the proper rotations and inversions of O
h
the
same way (x, y, z) do. An internal SU(3) triplet, by contrast, must transform under an SU(3)
action that commutes with O
h
, since SU(3) is an internal symmetry rather than a spatial one.
Identifying T
1u
with an internal triplet would break that commutativity. The two structures
are different objects.
7
defect
1
2
3
4
(a) Tetrahedral defect:
K
4
, 3 perfect matchings
K
4
perfect matchings
M
1
= {12, 34}
M
2
= {13, 24}
M
3
= {14, 23}
defect
1
1'
2
2'
3
3'
(b) Octahedral defect:
K
2, 2, 2
, 8 perfect matchings
K
2, 2, 2
edges (12)
antipodal pairs (not bonds)
one perfect matching
(8 total)
Figure 2: Bond graphs for the two interstitial defect classes of the FCC lattice. (a) Tetrahedral
defect: the four surrounding FCC sites form a regular tetrahedron, all six pairs at nearest-
neighbor distance, giving the complete graph K
4
. The three perfect matchings M
1
, M
2
, M
3
are
color-coded; the matter paper [14] builds the SU(3) color triplet on this 3-element matching
set. (b) Octahedral defect: the six surrounding FCC sites form a regular octahedron, but only
the twelve non-antipodal pairs sit at nearest-neighbor distance (gray edges); the three antipodal
pairs (red dotted) lie at second-nearest-neighbor distance
√
2 L and are not physical bonds. The
bonded graph is therefore K
2,2,2
, with eight perfect matchings rather than three (one matching
shown in green). The matter-paper color algebra does not extend to this graph, as shown in
Proposition 6 and independently in the dark-matter paper [15].
8
5.4 Combinatorial obstruction from K
2,2,2
The matter paper’s su(3) algebra uses perfect-matching combinatorics on the complete graph
of the four mutually nearest-neighbor valence bonds of a tetrahedral defect. For the octahedral
defect, the analogous graph is built from the six surrounding sites and the bonds joining those
pairs at nearest-neighbor distance in the FCC lattice. The six octahedron vertices fall into
three antipodal pairs at second-nearest-neighbor distance
√
2 L (with L the nearest-neighbor
distance), and these antipodal pairs are not bonded. Only the twelve non-antipodal pairs are
physical bonds. The bonded graph is therefore the complete tripartite graph
K
2,2,2
: 6 vertices in three antipodal pairs, 12 edges between non-antipodal pairs, (11)
not the complete graph K
6
on six vertices (compare Figure 2(a) and (b)). Direct enumeration
on K
2,2,2
(verified by the accompanying script) gives:
• Number of edges: 12.
• Number of perfect matchings (triples of disjoint edges covering all six vertices): 8.
• Number of skew-edge pairs (pairs of disjoint edges): 30.
The same counts appear in the companion dark-matter paper [15], where the skew-edge pair
count c
(O)
skew
= 30 enters as combinatorial input to the mass derivation.
Proposition 6. The matter-paper su(3) algebra constructed from the three perfect matchings
of K
4
does not extend to a structurally analogous algebra on the bond graph K
2,2,2
. K
2,2,2
has
eight perfect matchings, not three, and a skew-pair count of 30; neither factors as a three-color
representation [15]. No internal SU(3) algebra of the matter-paper form acts on the bound-mode
multiplet of an octahedral defect.
The bound-mode multiplet at an octahedral defect therefore carries the A
1g
⊕ E
g
⊕ T
1u
decomposition derived above, but no embedded internal triplet of the matter-paper type. The
same obstruction is reached combinatorially in [15]; the spectral analysis here is independent
representation-theoretic evidence.
5.5 Implication: a structurally distinct second sector
Tetrahedral defects support the matter-paper color construction. Octahedral defects do not.
The two defect classes are not equivalent matter sectors of the FCC substrate. Only the tetra-
hedral class, with coordination 4, is consistent with the SU(3) color algebra of the SSMTheory
program. The octahedral class is identified in the companion dark-matter paper [15] as a dark-
sector matter species; the O
h
multiplet decomposition of Proposition 5 and the bond-graph
identification of Eq. (11) are independent inputs consistent with that identification.
6 Comparison with Other Lattices and Reduced-Fermion Con-
structions
The defect mechanism of Sections 4–5 is unlike the established reduced-fermion constructions
in lattice gauge theory. We retain the standard naive Dirac operator on FCC; physical fermion
content is carried by localized excitations at substrate defects, not by changes to the bulk
operator. The point group of the defect, together with the bonded graph it generates, fixes the
multiplet structure. Table 1 summarizes the comparison across the lattices considered below.
9
HCP: same K = 12, broken isotropy along the c-axis. Hexagonal close-packed (HCP)
and FCC are the two close-packed lattices in 3D, both with coordination K = 12 and the
same nearest-neighbor distance. The first shells differ structurally. FCC’s first shell is the
cuboctahedron, with point group O
h
(order 48) including inversion symmetry. HCP’s first shell
is the anticuboctahedron (the Johnson solid J
27
, also called the triangular orthobicupola), with
point group D
3h
(order 12) and a distinguished c-axis. The two solids are related by rotating
one triangular face of the cuboctahedron by 60
◦
relative to its opposite face. That single twist
breaks O
h
to D
3h
and singles out the stacking direction.
The consequence at the bond-tensor level is sharp. Both lattices give T
(2)
µν
= K
2
δ
µν
, so
ω ∝ |
k| at long wavelengths in both cases. The rank-3 tensor T
(3)
µνρ
vanishes for both: by full
inversion symmetry for FCC, and by the horizontal mirror σ
h
acting on the anticuboctahedron
for HCP. At rank 4 the two diverge. For FCC, direct computation gives
T
(4), FCC
xxxx
= T
(4), FCC
yyyy
= T
(4), FCC
zzzz
,
T
(4), FCC
xxxx
T
(4), FCC
xxyy
= 2, (12)
which is anisotropic at rank 4 relative to the fully isotropic value 3, but with all three coordinate
axes equivalent. The unbroken symmetry of the rank-4 corrections is the full cubic group, so
the anisotropy is the standard discrete cubic anisotropy of any 3D lattice. For HCP at ideal
c/a =
p
8/3, direct computation gives
T
(4), HCP
xxxx
= T
(4), HCP
yyyy
=
5
2
= T
(4), HCP
zzzz
=
8
3
,
T
(4), HCP
xxxx
T
(4), HCP
xxyy
= 3,
T
(4), HCP
zzzz
T
(4), HCP
xxzz
= 4. (13)
The in-plane ratio 3 is what fully-isotropic rank-4 tensors give, but the c-axis breaks this: T
zzzz
differs from T
xxxx
, and the c-axis ratio differs from the in-plane ratio. HCP’s rank-4 anisotropy
contains a preferred-axis term forbidden by full SO(3), not merely a discrete cubic anisotropy.
A substrate intended to support an emergent rotationally invariant continuum theory cannot
single out a spatial direction at the O(a
2
) correction level of its bulk dispersion. HCP fails this
requirement, FCC passes it (subject to the standard discrete cubic anisotropy expected of any
3D lattice). Of the two close-packed 3D lattices, only FCC is therefore a candidate for a
substrate of isotropic physical spacetime. In four dimensions the analogous question singles out
the D
4
root lattice, whose 24-cell first shell gives an exactly isotropic rank-4 bond tensor.
At the local defect-bound-mode level relevant to the rest of this paper, HCP’s reduced
symmetry has a different consequence. The combinatorial bond graphs K
4
at a tetrahedral void
and K
2,2,2
at an octahedral void are preserved (they are determined by which pairs of bounding
vertices sit at nearest-neighbor distance, which depends only on the local coordination), so the
matter-paper K
4
perfect-matching construction still defines a three-element color basis at an
HCP tetrahedral defect. The local point group at the voids, however, is reduced: from T
d
to
a C
3v
subgroup at a tetrahedral void, and from O
h
to D
3d
at an octahedral void. The bound-
mode multiplet decompositions split more finely under the reduced symmetry. The A
1
⊕ T
1
structure of Section 4 refines under T
d
→ C
3v
as A
1
⊕A
2
⊕E (dimensions 1 + 1 + 2), with the
T
1
triplet of FCC splitting into a C
3v
singlet plus doublet. The T
d
-degeneracy among the three
colors is lifted by the broken symmetry. Whether this HCP refinement is physically interesting
in its own right is a question we do not pursue here; HCP is not a candidate substrate for the
SSMTheory framework on the rank-4 isotropy grounds described above.
BCC: different coordination, no analog. Body-centered cubic has K = 8 with cubic first
shell. Tetrahedral voids in BCC exist but are distorted; the four bounding vertices are not at
equal nearest-neighbor distance from one another (the bonded graph of a tetrahedral void in
BCC is not K
4
), and the T
d
local symmetry is broken. The matter-paper K
4
perfect-matching
construction does not directly apply. Gauge fields on BCC were studied by Celmaster [13];
10
the related fermion analysis on BCC, including whether a different local symmetry-based color
construction is available, is outside the scope of the present paper. The local structure of
FCC and HCP — not BCC or hypercubic Z
3
— is what supports the defect-bound-mode color
algebra.
Minimally doubled fermions. Karsten [6], Wilczek [7], and Boriçi–Creutz [8, 9] introduced
operators that produce exactly two continuum doublers, breaking hypercubic symmetry to
obstruct the doubling theorem. Bedaque, Buchoff, Tiburzi, and Walker-Loud [10] analyzed the
symmetry structure systematically. The modification is at the bulk-operator level. The defect
approach here is orthogonal: it leaves the bulk operator alone and looks at localized structure.
A future minimally-doubled construction on FCC could be combined with the defect analysis;
we do not undertake that here.
Supersymmetric lattices. Kaplan, Katz, and Ünsal [11] and Catterall [12] construct super-
symmetric lattice gauge theories on root lattices such as A
∗
4
, exploiting algebraic structure of
the root system to preserve part of the continuum supersymmetry. Defect-bound modes are not
in scope. The FCC lattice is the A
3
root lattice; a supersymmetric extension of the construction
here is consistent with these methods but is not pursued.
Wilson, staggered, twisted-mass. The standard hypercubic constructions [3, 4, 5] remove
or relabel doublers through mass terms, phase factors, or chirally rotated masses. All act on
the bulk operator on Z
D
. Staggered fermions on a 3D hypercubic lattice produce 4 tastes at the
eight cubic-BZ corners, all carrying Wilson mass 6/a. The FCC spectrum derived here has zeros
at the BZ interior (Γ, four L-points at (±π/2, ±π/2, ±π/2)) and on boundary nodal lines, with
three topological classes at three Wilson masses (0, 12/a, 16/a). The two Brillouin zones have
different shapes (truncated octahedron versus cube), and no linear change of lattice variables
maps one to the other while preserving the O
h
symmetry and nearest-neighbor structure. The
FCC operator is not a reparametrization of staggered fermions.
Lattice K First shell Point group Defect-mechanism status
Simple cubic (naive) 6 cube O
h
no analog (distorted tet void)
Simple cubic (staggered) 6 cube O
h
no analog
BCC (naive) 8 cube O
h
distorted tet void, K
4
broken [13]
HCP 12 anticuboctahedron D
3h
K
4
, K
2,2,2
preserved; reduced void symmetry; c-axis anisotropy
Minimally doubled varies varies varies bulk modification, no defect content
FCC (this work) 12 cuboctahedron O
h
K
4
, K
2,2,2
with full T
d
, O
h
at voids
Table 1: Comparison across common 3D lattices. Among the close-packed lattices, only FCC
has the cuboctahedral first shell that gives full O
h
symmetry without a preferred axis at the
rank-4 bond-tensor level. HCP shares the combinatorial bond graphs at its voids but breaks
isotropy along its c-axis at O(a
2
) in the bulk dispersion and reduces the local point group at
each void to a C
3v
or D
3d
subgroup. BCC has K
4
-broken tetrahedral voids and does not support
the matter-paper color construction directly.
FCC is uniquely selected by three independent criteria. A substrate compatible with
the matter-paper defect-and-color picture, with rotationally invariant continuum physics, and
with maximum-density packing must satisfy three requirements simultaneously: (i) full O
h
first-
shell symmetry, so the rank-4 bond tensor has no preferred axis beyond the standard discrete
cubic anisotropy of any 3D lattice; (ii) the K
4
bond graph at tetrahedral interstitial voids
and the K
2,2,2
bond graph at octahedral voids, so the matter-paper su(3) algebra on K
4
perfect
matchings is defined at the tetrahedral class; (iii) the maximum kissing number K = 12 in three
dimensions and the maximum Bravais sphere-packing density π/(3
√
2) ≈ 0.7405 established by
11
Kepler’s theorem (proved by Hales [20]), so the substrate is the densest 3D arrangement allowed
by sphere-packing geometry [19]. Each of the alternatives in Table 1 fails at least one of these.
Simple cubic (K = 6) and hypercubic Z
3
have full O
h
but no tetrahedral void with K
4
bond
structure, and K = 6 is not the maximum kissing number. BCC (K = 8) has full O
h
but its
tetrahedral voids are distorted, with bounding vertices at unequal distances and a K
4
-broken
bond graph; K = 8 is also below the maximum. HCP (K = 12) attains the maximum kissing
number and packing density and has the same combinatorial bond graphs at its voids as FCC,
but its anticuboctahedral first shell has only D
3h
symmetry and develops a preferred c-axis
at the rank-4 bond-tensor level, failing requirement (i). FCC is the unique 3D Bravais lattice
satisfying all three criteria. The construction in this paper is therefore not lattice-agnostic. It
picks out FCC specifically on independent geometric grounds (isotropy of the bulk dispersion),
combinatorial grounds (the K
4
structure at the tetrahedral void), and packing-theoretic grounds
(maximum kissing number and density). The same selection extends to four dimensions through
the D
4
root lattice, where the 24-cell first shell gives exact rank-4 isotropy and the FCC structure
is recovered on each constant-time slice.
7 Discussion
What this paper establishes. Two specific structural results. First, the bound-mode mul-
tiplet at a tetrahedral defect of the FCC lattice has dimension 4 and decomposes under T
d
as A
1
⊕ T
1
(Proposition 3), with the three-dimensional T
1
sub-multiplet dimensionally and
group-theoretically compatible with the SU(3) color triplet constructed in the matter paper
from K
4
perfect matchings (Proposition 4). Second, the bound-mode multiplet at an octahe-
dral defect has dimension 6 and decomposes under O
h
as A
1g
⊕E
g
⊕T
1u
(Proposition 5), with
the three-dimensional T
1u
sub-multiplet carrying the spatial vector representation rather than
an internal triplet; the corresponding bonded graph K
2,2,2
admits 8 perfect matchings and 30
skew-edge pairs, neither factoring as a three-color structure (Proposition 6). The matter-paper
color construction applies to tetrahedral defects and fails for octahedral defects on both spectral
and combinatorial grounds. This identifies tetrahedral defects with the quark-class matter of
the SSMTheory program and octahedral defects with the dark-sector species of the companion
dark-matter paper [15], whose analysis uses the same K
2,2,2
graph and skew-pair counts derived
here. The free-field bulk classification (Theorem 2) of three topological classes summing to
vanishing Nielsen–Ninomiya index supports this analysis as the spectral background.
What this paper does not establish. The Dirac analysis does not derive SU(3) color
from the spectrum; the matter paper’s construction is logically prior, using K
4
combinatorics
independently of the Dirac operator. The Dirac analysis confirms spectral compatibility for
tetrahedral defects and incompatibility for octahedral defects — a non-trivial check, but a
check rather than a derivation. The physical interpretation of the octahedral sector, including
its mass scale, cosmological abundance, and indirect-detection signatures, is the subject of the
companion dark-matter paper [15] and is not duplicated here.
Open questions. Two natural extensions. First is the dynamical analysis: turning on the
SU(3) gauge field of the SSMTheory program on the FCC bonds and computing how the bound-
mode multiplet responds to gauge fluctuations would close the dynamical picture. Second is
the 4D extension: an FCC slice × discrete time gives a D
4
root-lattice structure with exact
SO(4) isotropy in the long-wavelength limit; the analogous defect-bound-mode analysis on that
4D lattice would address the propagating-fermion content of the substrate.
12
Acknowledgments
I thank a referee for a careful reading of an earlier version of this work; that feedback prompted
the present reorganization around the defect-bound-mode analysis.
Data availability
A Python script reproducing all numerical and algebraic claims of this paper is included
with the submission as verify_defects.py; it runs 55 explicit assertions including the finite-
cluster Dirac diagonalization of Section 4, and exits with status 0 on success. The same script
is also available online at https://github.com/raghu91302/ssmtheory/blob/main/verify_
defects.py. The full LaTeX source, figure-generating Python scripts, and verification script
constitute the complete reproducibility package; no separate data files are required.
AI contributions statement
The author used Anthropic’s Claude (model: Claude Opus 4.7) as a writing and verification
assistant during preparation of this manuscript. Specifically, Claude was used to: refine prose
and improve clarity of explanations; check LaTeX builds and identify formatting issues; cross-
check internal consistency of citations, dimensions, and notation; and draft passages of the
supplementary documentation (README, comments in the verification script). All scientific
claims, mathematical results, and the verification-script logic — including the finite-cluster
Dirac diagonalization and its representation-theoretic interpretation — were derived, verified,
and accepted by the human author, who takes full responsibility for the manuscript’s content.
References
[1] H. B. Nielsen and M. Ninomiya, Nucl. Phys. B 185, 20 (1981); Erratum: Nucl. Phys. B
195, 541 (1982). doi:10.1016/0550-3213(81)90361-8.
[2] H. B. Nielsen and M. Ninomiya, Phys. Lett. B 105, 219 (1981). doi:10.1016/0370-
2693(81)91026-1.
[3] K. G. Wilson, Phys. Rev. D 10, 2445 (1974). doi:10.1103/PhysRevD.10.2445.
[4] L. Susskind, Phys. Rev. D 16, 3031 (1977). doi:10.1103/PhysRevD.16.3031.
[5] J. Kogut and L. Susskind, Phys. Rev. D 11, 395 (1975). doi:10.1103/PhysRevD.11.395.
[6] L. H. Karsten, Phys. Lett. B 104, 315 (1981). doi:10.1016/0370-2693(81)90133-7.
[7] F. Wilczek, Phys. Rev. Lett. 59, 2397 (1987). doi:10.1103/PhysRevLett.59.2397.
[8] A. Boriçi, Phys. Rev. D 78, 074504 (2008). doi:10.1103/PhysRevD.78.074504.
[9] M. Creutz, JHEP 04, 017 (2008). doi:10.1088/1126-6708/2008/04/017.
[10] P. F. Bedaque, M. I. Buchoff, B. C. Tiburzi, and A. Walker-Loud, Phys. Lett. B 662, 449
(2008). doi:10.1016/j.physletb.2008.03.020.
[11] D. B. Kaplan, E. Katz, and M. Ünsal, JHEP 05, 037 (2003). doi:10.1088/1126-
6708/2003/05/037.
[12] S. Catterall, JHEP 06, 027 (2008). doi:10.1088/1126-6708/2008/06/027.
13
[13] W. Celmaster, Phys. Rev. D 26, 2955 (1982). doi:10.1103/PhysRevD.26.2955.
[14] R. Kulkarni, Matter as Incomplete Crystallization: Quark Charges, Color Confinement,
and the Proton Mass from a Single Extra Node in the Vacuum Lattice, Phys. Open 27,
100423 (2026). doi:10.1016/j.physo.2026.100423.
[15] R. Kulkarni, Dark Matter as a Trapped K = 6 Remnant in the Octahedral Voids of the
FCC Vacuum Lattice, submitted to Phys. Open (2026). doi:10.5281/zenodo.20047901.
[16] L. H. Karsten and J. Smit, Nucl. Phys. B 183, 103 (1981). doi:10.1016/0550-3213(81)90549-
6.
[17] H. S. M. Coxeter, Regular Polytopes, 3rd ed., Dover, New York (1973).
[18] N. W. Ashcroft and N. D. Mermin, Solid State Physics, Saunders College, Philadelphia
(1976).
[19] J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, 3rd ed., Springer,
New York (1999). doi:10.1007/978-1-4757-6568-7.
[20] T. C. Hales, A proof of the Kepler conjecture, Annals of Mathematics 162, 1065 (2005).
doi:10.4007/annals.2005.162.1065.
14
