A Self-Correcting Discrete Vacuum: Baryon Conservation and Color, from an Emergent Qutrit to Temporal su(3) Closure, on One FCC Code Complex

A Self-Correcting Discrete Vacuum:
Baryon Conservation and Color, from an Emergent Qutrit
to Temporal su(3) Closure, on One FCC Code Complex
Raghu Kulkarni
SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA
raghu@idrive.com
June 2026
Abstract
A single FCC code complex on a face-centered-cubic vacuum, with a quark a node trapped
in a tetrahedral void, protects both the conserved quantities of baryonic matter and its color,
and we develop both sectors at one standard of rigor, every quantitative claim reproduced by
an accompanying script. The [[192, 130, 3]] binary CSS code with weight-twelve vertex and
octahedral-void stabilizers, together with its Z
3
triality refinement on the same chain complex,
is the common origin. In the external sector it makes the defect’s winding charge a Gauss-law
super-selection invariant, yielding the baryon ladder
, n, p,
++
with charges Q = 1 + W ;
the void orientation, by contrast, is not protected (a negative result we record); the migrating
defect obeys E = γE
rest
; and migration is a code-corrected string. In the internal sector
the same complex protects color. A color qutrit emerges as the symmetric sector of the cage
kinetic term T = (J
3
I
3
) J
2
; the three skew-pair matchings carry an S
3
Weyl symmetry
and a rank-two torus, and among the rank-two simple Lie algebras both the Weyl order and the
existence of a three-dimensional representation independently select A
2
=
su(3); anchor-induced
color/dark leakage is a weight-one correctable error. The FCC complex is then gauge-ready—
its minimal plaquettes are triangles, each carrying one bond of every color—and supports a
genuine Z
3
triality code, whose vertex Gauss law enforces q
0
+ q
1
+ q
2
0 mod 3 and which
subsumes the octahedral baryon operator; the non-abelian SU(3) Gauss law at a vertex has a 513-
dimensional gauge-invariant space gated by triality, and the triangular plaquette Hamiltonian is
gauge-invariant, so a Kogut-Susskind SU(3) Hamiltonian is well-posed on the FCC plaquettes.
Building the color connection from the vacuum’s processes gives a layered dynamical result.
The spatial channels transport color only diagonally, an N(T ) = T S
3
connection short of
SU(3) by the root generators because the colors are bond-disjoint. We then test whether the
stabilizer correction cycle removes this obstruction—whether a syndrome-conditioned recovery,
dressing the bare transport, induces an off-diagonal logical map—and find on the real code that
it does not: the minimum-weight recovery of a leakage event is single-bond and same-color, so
the QEC-dressed transporter stays in N(T ). The off-diagonal direction is instead supplied by
the temporal link, the cage time evolution with off-diagonal generator J
3
I
3
, which with the
torus closes the full eight-dimensional su(3). Leading strong-coupling confinement holds on the
triangular plaquettes with a positive string tension, and the single-plaquette spectrum has the
fundamental Casimir gap C
2
(3) = 4/3. The color algebra, both protected sectors, and the well-
posed Hamiltonian are thus attained on one code complex, with the full color algebra appearing
only when temporal qutrit evolution is included; the continuum dynamics—confinement in the
continuum limit, the running coupling, Λ
QCD
—are not addressed and remain open.
1
Contents
1 Introduction 3
2 The vacuum code 3
2.1 Lattice, qubits, and stabilizers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Commutation and the CSS property . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.3 Logical count and distance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.4 The [4, 4, 4] vertex split . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3 External sector I: orientation is not protected 5
4 External sector II: baryon charge as a Gauss-law super-selection rule 5
4.1 Charge as enclosed flux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
4.2 The winding and the baryon ladder . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
5 External sector III: relativistic kinematics 6
6 External sector IV: migration as code-corrected transport 7
7 Internal sector I: the emergent color qutrit 8
7.1 Three colors from K
4
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
7.2 The kinetic-term factorization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
8 Internal sector II: Weyl group, torus, and the selection of A
2
9
8.1 The Weyl group from tetrahedral symmetry . . . . . . . . . . . . . . . . . . . . . . . 9
8.2 The maximal torus from bond phases . . . . . . . . . . . . . . . . . . . . . . . . . . 9
8.3 Selection of A
2
among the rank-two simple algebras . . . . . . . . . . . . . . . . . . 9
9 Internal sector III: code stabilization of color 9
10 Internal sector IV: the gauge-ready geometry 10
10.1 Triangular plaquettes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
10.2 Color-completeness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
11 Internal sector V: the Z
3
triality code 11
11.1 The chain complex and the CSS condition . . . . . . . . . . . . . . . . . . . . . . . . 11
11.2 Logical count and homology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
11.3 Triality, correctability, and subsumption . . . . . . . . . . . . . . . . . . . . . . . . . 11
12 Internal sector VI: a well-posed Kogut-Susskind SU(3) Hamiltonian 12
12.1 The non-abelian Gauss law and the 513-dimensional vertex space . . . . . . . . . . . 12
12.2 The magnetic term is gauge-invariant . . . . . . . . . . . . . . . . . . . . . . . . . . 12
13 Internal sector VII: the dynamical connection—spatial obstruction 12
13.1 The migration channel realizes the Weyl group . . . . . . . . . . . . . . . . . . . . . 12
13.2 Spatial channels reach only N(T ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
13.3 The other spatial channels are color-diagonal . . . . . . . . . . . . . . . . . . . . . . 13
14 Internal sector VIII: QEC-dressed logical transport is tested and remains diag-
2
onal 13
14.1 The dressed transporter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
14.2 The test on the real code . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
14.3 Consequence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
15 Internal sector IX: temporal resolution—the root direction from cage evolution 14
15.1 The temporal link is color-mixing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
15.2 The temporal generator closes su(3) . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
16 Internal sector X: dynamical footholds 15
16.1 Strong-coupling confinement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
16.2 The single-plaquette spectrum . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
17 Unification: one code complex, two sectors 17
18 Scope and status 17
1 Introduction
The strong sector of the Standard Model is specified by ingredients taken as input: the gauge group
SU(3), the color number N
c
= 3, baryon-number conservation as an accidental global symmetry,
and the confining dynamics. The Selection-Stitch Model proposes a common origin for these in a
discrete vacuum—an FCC entanglement network carrying a stabilizer code, in which a quark is a
node trapped at the centroid of a tetrahedral void [1]. This paper assembles, in one place and at
one standard of rigor, what that proposal establishes about the strong sector, organized around
a single claim: that one FCC code complex—the same chain complex in its binary distance-three
realization and its Z
3
triality refinement—protects two sectors.
The thesis is that the [[192, 130, 3]] FCC stabilizer code [2] is simultaneously the origin of baryon-
number conservation (the external sector: charge, kinematics, migration) and of color (the internal
sector: the qutrit, the algebra, the gauge connection). Figure 1 shows the architecture. The external
sector occupies Sections 26; the internal sector occupies Sections 716, ascending through three
levels—representation, kinematics, and dynamics. Section 17 states the unification, and Section 18
the boundary: the color algebra, both sectors, and the well-posed theory are reached; the confining
continuum dynamics are not.
The discipline throughout is to state, at each step, exactly what is computed, and to mark negative
results and open questions as plainly as positive ones. Every quantitative claim is reproduced by a
script in the data-availability archive.
2 The vacuum code
2.1 Lattice, qubits, and stabilizers
The vacuum is the FCC lattice [12], with the four-atom conventional basis (0, 0, 0), (
1
2
,
1
2
, 0), (
1
2
, 0,
1
2
), (0,
1
2
,
1
2
).
We work on a periodic 2 × 2 × 2 block, which contains N = 32 sites. The nearest neighbors of a
site are the twelve 110/
2 directions at distance a/
2 (the FCC kissing number 12), and on the
block there are E = 192 distinct nearest-neighbor bonds. A qubit is placed on each bond, giving a
192-qubit code space.
3
FCC code complex
binary [[192, 130, 3]] code + Z
3
triality refinement
External sector (baryon number) Internal sector (color)
charge Q = 1 + W (Gauss law)
baryon ladder
, n, p,
++
E = γE
rest
; migration string
qutrit, S
3
, torus, A
2
=su(3)
Z
3
triality; well-posed KS Hamiltonian
N(T ) spatially su(3) with time
Figure 1: One FCC code complex, two protected sectors. The same chain complex—in its binary
[[192, 130, 3]] realization and its Z
3
triality refinement—makes the defect’s winding charge a Gauss-
law invariant (external) and protects the emergent color qutrit (internal); the internal sector ascends
from representation through kinematics to dynamics.
Two stabilizer families are defined. To each vertex v associate the star operator
A
v
=
Y
estar(v)
Z
e
, |star(v)| = 12, (1)
the product of Z over the twelve bonds incident to v. The octahedral voids of the FCC lattice
(the octahedral interstitial sites) number 32 on the block; the six atoms nearest a void o form an
octahedron whose twelve edges are nearest-neighbor bonds, and to each void associate the cage
operator
B
o
=
Y
eoct(o)
X
e
, |oct(o)| = 12, (2)
the product of X over those twelve bonds. Both families are weight-twelve.
2.2 Commutation and the CSS property
A
v
is a product of Z and B
o
a product of X, so they commute if and only if their supports overlap
in an even number of bonds. A vertex v lying on the octahedron of void o contributes the bonds
incident to v within that octahedron; a direct enumeration shows this overlap is either 0 (when v is
not a vertex of the octahedron) or exactly 4 (when it is), both even. Hence [A
v
, B
o
] = 0 for all v, o,
and the code is a Calderbank-Shor-Steane (CSS) code [5, 6], with H
Z
the vertex–bond incidence
and H
X
the void–bond incidence over F
2
.
2.3 Logical count and distance
Each bond is incident to exactly two vertices, so every bond appears in exactly two star operators;
therefore
Q
v
A
v
= (each Z
e
appears twice and Z
2
e
= ), one relation among the 32 Z-stabilizers,
leaving 31 independent. Each bond borders exactly two octahedral voids, so likewise
Q
o
B
o
= ,
leaving 31 independent X-stabilizers. The independent stabilizer count is 31 + 31 = 62, and the
logical-qubit count is
k = E 62 = 192 62 = 130, (3)
the code [[192, 130, 3]]. For the distance: a single-bond Z error on bond e anticommutes with
exactly the cage operators of the two voids bordering e, a weight-two X-syndrome that is distinct
for each bond, so single-bond Z errors are detectable and correctable; symmetrically a single-bond
4
X error flags the two vertices of e. Weight-one errors are thus correctable, consistent with distance
three.
2.4 The [4, 4, 4] vertex split
The twelve bonds incident to a vertex lie in the three coordinate planes—four in the xy plane, four
in xz, four in yz—giving a [4, 4, 4] partition of the star. This partition is invisible in the external
sector but is the seed of color: in Section 8 the three classes are the three colors, and in Section 11
the Z
3
vertex Gauss law factorizes through it into the triality condition. The two stabilizer types
are the two halves of a discrete gauge structure, in the manner of the toric code [7]: the vertex
Z-stars are Gauss-law generators, the octahedral X-cages are magnetic operators, and the same
FCC chain complex below carries both sectors, through its binary distance-three realization and
its Z
3
triality refinement.
3 External sector I: orientation is not protected
Before the positive external results we record a natural candidate invariant that fails; the failure
sharpens what the code protects. A tetrahedral void has an orientation—the sense in which its
trapped node binds its four cage vertices—and one might expect it to be topologically protected. It
is not. Consider the two ways a relocation can treat the orientation. Preserving it keeps the node
bound through the six second-neighbor bonds of length a/
2 that maintain the original binding
sense; flipping it rebonds the node through the four nearest cage bonds of length a/2. The elastic
cost scales with total bond length, so
cost(flip)
cost(preserve)
=
4 · (a/2)
6 · (a/
2)
=
2a
6a/
2
=
2
3
0.667. (4)
The orientation-flipping relocation is therefore cheaper, by a factor 2/3, than the orientation-
preserving one. Orientation is not a protected quantity. The lesson is that protection in this
vacuum is not conferred by local geometric labels, which can be relaxed by cheaper rebondings, but
by global Gauss-law fluxes, which cannot; the protected external invariant is the winding charge of
the next section.
4 External sector II: baryon charge as a Gauss-law super-selection
rule
4.1 Charge as enclosed flux
The protected external quantity is the defect’s winding charge, and it is protected because it is
a Gauss-law flux. The vertex stars A
v
generate a discrete Gauss law: on a state with a trapped
defect, the product of A
v
over a closed surface Σ enclosing the defect equals the net Z-flux through
Σ, and this is a stabilizer-measured, conserved quantity. It is conserved because the boundary
map’s rows sum to zero:
Q
v
A
v
= globally (Section 2), so the total enclosed flux cannot change
under any local operation that does not cross Σ. The charge is therefore a super-selection label,
invariant under every local rebonding interior to Σ—unlike the orientation, which is a local label
and is not invariant.
5
4.2 The winding and the baryon ladder
A trapped node bonds to four cage vertices; three carry the valence structure and one is the anchor
fixing the baseline. Each valence bond carries a binary winding w
i
{0, 1} (whether its Bell-pair
phase has wound), and the net winding is W =
P
3
i=1
w
i
{0, 1, 2, 3}. The enclosed flux is the
anchor baseline plus the valence winding,
Q = 1 + W, (5)
the 1 from the anchor and +W from the valence bonds. The per-quark charge is correspondingly
q
i
=
1
3
+ w
i
, fractional because the anchor baseline 1 is shared democratically across the three
colors. Summing,
P
i
q
i
= 1 + W = Q, consistent with Eq. (5). The four winding values give
the baryon ladder of Figure 2 and Table 1: W = 0 gives
with all quarks at
1
3
and Q = 1;
W = 1 the neutron, one quark at +
2
3
and two at
1
3
, Q = 0; W = 2 the proton, Q = +1; W = 3
the
++
, all quarks at +
2
3
, Q = +2. The integer quark charges and the baryon ladder are thus
super-selection consequences of Eq. (5), not separate inputs.
baryon W quark charges (
1
3
+ w
i
) Q = 1 + W
0 (
1
3
,
1
3
,
1
3
) 1
n 1 (+
2
3
,
1
3
,
1
3
) 0
p 2 (+
2
3
, +
2
3
,
1
3
) +1
++
3 (+
2
3
, +
2
3
, +
2
3
) +2
Table 1: The baryon ladder generated by the winding W . Each unit of winding flips one valence
quark from
1
3
to +
2
3
and raises Q by one.
0 1 2 3
winding W =
X
i
w
i
−1
0
1
2
electric charge Q = ¡1 + W
¢
¡
(¡
1
3
; ¡
1
3
; ¡
1
3
)
n
(+
2
3
; ¡
1
3
; ¡
1
3
)
p
(+
2
3
; +
2
3
; ¡
1
3
)
¢
+ +
(+
2
3
; +
2
3
; +
2
3
)
Figure 2: The baryon spectrum as a Gauss-law super-selection ladder; Q = 1+W , with per-quark
charges shown.
5 External sector III: relativistic kinematics
A defect set in motion is a soliton of the bond dynamics. In the long-wavelength limit the bond-
phase field obeys a sine-Gordon-type equation, whose static kink is the defect; in lattice units its
6
rest energy is E
rest
= 8. The sine-Gordon kink is a Bogomolnyi (topological) soliton [11], for which
the boosted energy is exactly the Lorentz-transformed rest energy,
E(v) = γ E
rest
, γ =
1
p
1 v
2
/c
2
, (6)
with c the lattice signal speed. We verify this on the lattice: boosting the discrete kink profile
by velocity v and evaluating the lattice energy functional reproduces γE
rest
to six decimal places
across v/c [0, 0.95] (Figure 3). The relation is exact for the one-dimensional Bogomolnyi soliton;
a fully three-dimensional defect inherits E = γE
rest
by Lorentz covariance of the long-wavelength
effective theory, which is the level at which the claim is made, and which we do not strengthen to
an exact three-dimensional statement.
0.0 0.2 0.4 0.6 0.8
defect velocity v=c
1.0
1.5
2.0
2.5
3.0
E=E
rest
= °
boosted lattice soliton
matches °E
rest
to 6 d.p.
Figure 3: The migrating defect’s energy follows E = γE
rest
with E
rest
= 8, verified against the
lattice soliton boost to six decimal places.
6 External sector IV: migration as code-corrected transport
Baryon transport is defect migration, and migration is realized as a string of single-bond distur-
bances that the code corrects en route. Translating the defect by one lattice step toggles the bonds
along its path. By Section 2 each toggled bond is a weight-one error with a definite two-void (or
two-vertex) syndrome, hence correctable; the correction repairs the local disturbance while leaving
the enclosed flux—the charge of Eq. (5)—unchanged, since the migration string does not cross the
surface defining Q. Migration is therefore code-assisted string transport: the external charge is
conserved not by a continuous symmetry imposed by hand but by the error-correcting structure of
the vacuum, which detects and undoes every local step that would otherwise leak it. The string
transports the charge as a logical invariant, in direct analogy with the transport of a logical op-
erator along a corrected path in a stabilizer code. This is the external half of the “one code, two
sectors” thesis; the internal half follows.
7
7 Internal sector I: the emergent color qutrit
7.1 Three colors from K
4
The same defect carries color. It is bounded by four sites through four valence bonds, and the
complete graph K
4
on the bonds admits exactly three perfect matchings into skew pairs,
M
1
= {AB, CD}, M
2
= {AC, BD}, M
3
= {AD, BC}, (7)
so the color number is the rigid combinatorial count
N
c
=
1
2
4
2
= 3, (8)
with no fourth pairing, and none for an odd cage (K
3
has no perfect matching). The dimension of
the color Hilbert space is fixed before any algebra is written (Figure 4).
A
B
C
D
node
M
1
= {AB, CD}
M
2
= {AC, BD}
··· M
3
= {AD, BC}
three skew pairings
= three colors
Figure 4: The trapped node bonds to the four cage vertices A, B, C, D; K
4
has exactly three perfect
matchings into skew pairs (the three colors), and the pairs are bond-disjoint.
7.2 The kinetic-term factorization
Color is not merely a label; it emerges as a dynamical sector. The cage degrees of freedom factor as
the three-matching space Z
3
tensored with the two-state space of each skew pair (the two members
of a pair). On this 3 × 2 space the nearest-neighbor kinetic term, which hops amplitude between
bonds sharing a vertex, takes the factorized form
T = (J
3
I
3
) J
2
, (9)
where J
3
is the 3×3 all-ones matrix (hopping between distinct matchings), I
3
the identity (removing
the diagonal), and J
2
=
1 1
1 1
the all-ones matrix on the two members of a skew pair. The
eigenstates of J
2
are the symmetric and antisymmetric combinations
|+
i
=
1
2
|e
i1
+ |e
i2
, |−
i
=
1
2
|e
i1
|e
i2
, (10)
with J
2
|+
i
= +2|+
i
and J
2
|−
i
= 0, so the antisymmetric states are annihilated by the hopping.
Since T = (J
3
I
3
) J
2
, the symmetric sector and antisymmetric sector decouple, and on the
symmetric (color) sector T restricts (absorbing the factor 2 into the hopping rate Γ) to
T
color
= J
3
I
3
=
0 1 1
1 0 1
1 1 0
, (11)
a Hamiltonian on the qutrit H
C
= span
C
{|+
1
, |+
2
, |+
3
⟩} with six color transitions |+
i
|+
j
.
The antisymmetric states |−
i
are dark: they carry eigenvalue 0 of J
2
and decouple from the color
dynamics entirely. The qutrit is therefore the symmetric sector of the cage hopping, derived from
Eq. (9) rather than assumed. The structural fact used repeatedly below is that the three matchings
are bond-disjoint: each cage bond belongs to exactly one skew pair, hence to exactly one color.
8
8 Internal sector II: Weyl group, torus, and the selection of A
2
8.1 The Weyl group from tetrahedral symmetry
The symmetric group S
4
acts on the four bond labels. The Klein four-group V
4
S
4
of double
transpositions {e, (AB)(CD), (AC)(BD), (AD)(BC)} fixes each matching as a set, so the induced
action on the three matchings is through the quotient
S
4
/V
4
=
S
3
, (12)
which permutes {M
1
, M
2
, M
3
} as the full symmetric group on three objects. This S
3
is the Weyl
group of A
2
. Direct enumeration of the 24 bond permutations confirms exactly the six color
permutations are realized.
8.2 The maximal torus from bond phases
The three valence bonds carry Bell-pair phases ϕ
1
, ϕ
2
, ϕ
3
. The diagonal combination ϕ
1
+ ϕ
2
+ ϕ
3
is the overall U(1) phase, which is not color and is identified with electromagnetism elsewhere in
the program; the two independent traceless combinations form the rank-two maximal torus T
2
of
the color algebra, with the Cartan generators λ
3
, λ
8
as their conjugate charges.
8.3 Selection of A
2
among the rank-two simple algebras
The data—a three-dimensional representation with Weyl group S
3
and a rank-two torus—select
A
2
=
su(3), and the selection is sharp. There are exactly three rank-two simple Lie algebras [9],
and the geometry supplies two independent data, each of which excludes the other two (Table 2).
The Weyl-group orders are 6, 8, 12, so the order-six S
3
is the Weyl group of A
2
alone. The smallest
nontrivial representation dimensions are 3, 4, 7, so a three-color (three-dimensional) representation
exists only for A
2
. Either datum excludes B
2
and G
2
.
algebra Weyl-group order smallest nontrivial rep selected by
A
2
= su(3) 6 (S
3
) 3 both
B
2
=
so(5) 8 4 neither
G
2
12 7 neither
Table 2: Among the rank-two simple Lie algebras, both the S
3
Weyl symmetry (order 6) and the
existence of a three-dimensional representation independently select A
2
=
su(3).
This is a selection of the rank-two representation and Weyl structure, not a derivation of the
continuous gauge dynamics. Once a genuine qutrit exists, the traceless operators on it close on
su(3) automatically; that closure is not independent evidence and is not offered as such.
9 Internal sector III: code stabilization of color
The color qutrit is a protected logical subspace of the code of Section 2. The anchor, needed to fix
the charge baseline and the Weyl structure, breaks the cage symmetry between the two members
of a skew pair and so couples the color qutrit to the dark sector. Modeling the anchor as an on-site
asymmetry between the bonds of a pair, the leakage operator connects |+
i
to |−
i
at a rate set by
9
the anchor strength,
Γ
leak
=
Γ
A
Γ
2
, (13)
with Γ
A
the anchored and Γ the symmetric hopping rate; the leakage vanishes in the symmetric
limit Γ
A
= Γ. The color/dark exchange is the exchange parity of a single skew pair, which is
a single-bond Z operator, hence a weight-one error of the code—exactly the correctable class at
distance three, flagged by the bordering-void syndrome of Section 2.
The protection operates in two regimes. Without active correction, the steady-state dark population
is perturbatively bounded by
leak
/Γ)
2
: below one percent for a weak anchor
A
/Γ 1.2, giving
0.7%), rising to about fifteen percent only at the strong value Γ
A
/Γ = 2. With periodic stabilizer
projection, the leakage is suppressed by the quantum Zeno effect when the correction rate exceeds
the leakage rate, γ
corr
Γ, the dark population then falling as γ
2
corr
. Color is, in this precise and
quantified sense, a protected logical degree of freedom, on the same code that protects the external
charge—but the protection is conditional on the anchor being weak or the correction fast, a caveat
we state rather than suppress.
10 Internal sector IV: the gauge-ready geometry
10.1 Triangular plaquettes
A lattice gauge theory is built on plaquettes, the minimal closed loops of the bond graph. The
FCC nearest-neighbor graph is close-packed, and its minimal loops are triangles: three mutually
adjacent sites. On the block there are 256 such triangular plaquettes, and each nearest-neighbor
bond lies in exactly four of them, so the triangles tile the two-skeleton uniformly. A gauge action
on the FCC vacuum is therefore intrinsically triangular rather than square; this is unusual but not
pathological.
10.2 Color-completeness
The triangles are color-complete: every one of the 256 carries exactly one bond of each of the three
coordinate-plane color classes. To see this, take a triangle on sites 0, (
1
2
,
1
2
, 0), (
1
2
, 0,
1
2
); its three edges
have displacement vectors (
1
2
,
1
2
, 0), (
1
2
, 0,
1
2
), and (0,
1
2
,
1
2
), which lie in the xy, xz, and yz planes
respectively—one of each class. The same holds for all 256 triangles. Consequently a plaquette
term Re Tr(U
ab
U
bc
U
ca
) on an FCC triangle couples one bond of each color, and the elementary
loops are not color-blind (Figure 5)—the geometric prerequisite for color-mixing dynamics.
R
GB
each of the 256 triangular plaquettes
carries one bond of every color
Figure 5: A color-complete triangular plaquette: every FCC triangle has one bond from each
coordinate-plane color class (R, G, B), so a Wilson plaquette term mixes the three colors by con-
struction.
10
11 Internal sector V: the Z
3
triality code
11.1 The chain complex and the CSS condition
The binary code of Section 2 protects only the Z
2
reduction of triality. The triangular two-cells lift
it to genuine Z
3
. We build the cellular chain complex over Z
3
with 0-cells the vertices, 1-cells the
bonds (now carrying qutrits with the Z
3
clock and shift operators), and 2-cells the 256 triangles.
The boundary maps are
1
(bond its two endpoints, with +1 at the head and 1 at the tail)
and
2
(triangle its three bonds, with cyclic orientation). The vertex Z-stabilizers are the rows
of
1
and the triangular X-stabilizers the rows of
2
, so the CSS condition is
H
X
H
Z
=
2
1
= (
1
2
)
= 0 mod 3, (14)
which holds automatically, being the topological identity
1
2
= 0. (By contrast the octahedral-
void operators do not survive the Z
3
lift: the four-bond vertex–void overlap is 4 1 mod 3, and
the naive qutrit generalization fails to commute, with 152 violations on the block.)
11.2 Logical count and homology
A rank computation over GF(3) gives rank H
Z
= 31 (the boundary rank of a connected graph,
N 1) and rank H
X
= 158, so
k = 192 31 158 = 3, (15)
equal to the first Betti number b
1
(T
3
) = 3 of the three-torus: the three logical qutrits are the three
independent noncontractible cycles, the correct homological count rather than an artifact. The
integer complex closes as well,
1
2
= 0 over Z, so an integer-charge code exists with triality as its
mod-three reduction; comparing dim H
1
over Q, F
2
, and F
3
gives 3 in every case, so H
1
= Z
3
is
free with no torsion at any prime. Triality here is therefore the mod-three reduction of the integer
charge, not a topologically forced class—a negative we record rather than overstate.
11.3 Triality, correctability, and subsumption
At each vertex the twelve incident bonds split [4, 4, 4] by color (Section 2), and the Z
3
vertex Gauss
law factorizes through this split into the three class charges q
0
, q
1
, q
2
, imposing
q
0
+ q
1
+ q
2
0 mod 3, (16)
genuine triality—the center Z
3
SU(3) governing color superselection—rather than the Z
2
shadow
of the binary code. Single-bond Z errors have weight-four triangular syndromes (the four triangles
containing the bond), all 192 distinct, hence correctable; the total triality charge is conserved since
the boundary rows sum to zero. Finally, the weight-twelve octahedral baryon operator, suitably
Z
3
-oriented, is a boundary in the triangular complex—it lies in the image of
2
, hence is a product
of triangular stabilizers and is already implied by them. The triangular Z
3
code therefore subsumes
the octahedral structure: one Z
3
code carries both the triality (internal) and the baryon (external)
sectors.
11
12 Internal sector VI: a well-posed Kogut-Susskind SU(3) Hamil-
tonian
12.1 The non-abelian Gauss law and the 513-dimensional vertex space
Triality is the abelian center; the non-abelian question is whether SU(3) link variables admit a
consistent gauge-invariant Hamiltonian in the Kogut-Susskind formulation [4] of Wilson’s lattice
gauge theory [3]. Place a fundamental on each bond, oriented so a vertex sees 3 on outgoing and
¯
3
on incoming links. A twelve-link vertex with six of each orientation has gauge-invariant (physical,
uncharged) states given by the singlets in 3
6
¯
3
6
. Decomposing 3
6
into irreps by iterated
fusion and summing the squared multiplicities (equivalently, integrating |χ
3
|
12
over SU(3)) gives a
singlet multiplicity of
dim
singlets in 3
6
¯
3
6
= 513. (17)
The local gauge-invariant Hilbert space at a vertex is 513-dimensional. It is nonzero precisely
because the net triality is 6 6 0 mod 3; an unbalanced vertex such as 3
5
(triality 2) has no
singlet, and gauge invariance is forbidden. The solvability of the non-abelian Gauss law is thus
controlled, at the center, by the same triality condition (16) that the abelian code enforces: the
qutrit code is the center-reduction of the non-abelian Gauss law.
12.2 The magnetic term is gauge-invariant
A Kogut-Susskind Hamiltonian adds the magnetic (plaquette) operator Re Tr(U
ab
U
bc
U
ca
) on each
triangular plaquette, and consistency requires it to commute with the vertex Gauss-law generators
G
a
v
. Because a triangular plaquette is a closed oriented loop on which each of its three vertices
touches exactly two links, an infinitesimal gauge rotation g
v
at a vertex acts on the two incident
plaquette links by left- and right-multiplication and cancels inside the trace. Explicitly, applying
the eight SU(3) generators as infinitesimal vertex rotations to a random plaquette holonomy gives
|d Tr(W )| 10
16
, zero to machine precision. A gauge-invariant Kogut-Susskind SU(3) Hamilto-
nian is therefore constructible on the triangular FCC plaquettes: the Hamiltonian construction is
well-posed (a consistent operator on the gauge-invariant space, not a derivation of QCD dynamics).
13 Internal sector VII: the dynamical connection—spatial obstruc-
tion
13.1 The migration channel realizes the Weyl group
A well-posed Hamiltonian is kinematics; a dynamical gauge structure needs a transporter U
xy
between neighboring color fibers, built from the vacuum’s processes. The elementary migration
channel, the orientation-flipping hop, identifies three cage vertices and swaps one, inducing a per-
mutation of the four valence-bond labels and hence a permutation of the three matchings. Enu-
merating the 24 bond permutations realizes exactly the six color permutations, the Weyl group
S
3
. The hop transports color by relabeling it—a Weyl element—not by rotating it through the
continuous group.
13.2 Spatial channels reach only N(T )
Combining the Weyl permutations with the Cartan torus T
2
of Section 8, the reachable transporters
are exactly the monomial matrices (one nonzero entry per row and column), which form the
12
normalizer of the torus,
T, S
3
= N(T ) = T S
3
SU(3). (18)
A scan of 2×10
4
random torus–Weyl products confirms every one is monomial. The links are N(T )-
valued: a genuine real-space color connection, with U
xy
7→ G
x
U
xy
G
1
y
under N (T) frame changes
and gauge-invariant closed-loop holonomy Tr U
γ
. It is not full SU(3): the monomial matrices
are a measure-zero, two-dimensional (plus six-sheet) subset of the eight-dimensional group, and a
continuous root-direction element e
itλ
1
is not monomial.
13.3 The other spatial channels are color-diagonal
Neither of the remaining spatial channels supplies the missing off-diagonal direction. Single-bond
error correction: the recoveries are single-bond Pauli operators, and projected onto the qutrit they
are diagonal—a phase correction Z
b
acts as the identity when b is not in color i, and maps |+
i
to
its dark partner (projected away) when it is; generating the algebra of all single-bond recoveries
yields only diagonal maps (maximum off-diagonal element zero). This is examined in full, on the
real code, in Section 14, where the conjecture that the correction cycle promotes the transporter to
full rank is stated and tested. The anchor: H
anchor
= I
3
σ
z
is off-diagonal in the color/dark axis,
but the induced color-to-color map obtained by leaking to the dark sector and recovering back,
(P
C
H
anchor
P
D
)(P
D
H
anchor
P
C
), is the identity—off-diagonal color content zero. The common cause
is bond-disjointness: every operator built from the spatial cage bonds touches at most one color
and is therefore color-diagonal on the qutrit; the off-diagonal root generators E
α
cannot be built
from bond operators at all (supplying them through enlarged link Hilbert spaces is the route taken
by quantum-link models [8]). The obstruction is genuine, but it constrains operators built from the
spacelike bonds and says nothing about the temporal direction.
14 Internal sector VIII: QEC-dressed logical transport is tested
and remains diagonal
The spatial obstruction reaches only N(T ) for bare transport. A natural conjecture is that the
obstruction is an artifact of bareness—that the physical transport is not bare stitching but stitch-
ing dressed by the stabilizer correction cycle, and that the syndrome-conditioned recovery could
promote the rank-deficient bare transporter to a full-rank one carrying the root generators. We
state the conjecture precisely and test it on the real code; it fails.
14.1 The dressed transporter
Let H
C,x
= span{| +
1
(x), | +
2
(x), | +
3
(x)⟩} be the protected color qutrit at defect site x,
with projector P
C,x
. Bare cage-to-cage transport is M
bare
xy
= P
C,x
S
xy
P
C,y
, with S
xy
the physical
stitching/migration map; by Section 13 it is monomial. The QEC-dressed transporter inserts the
recovery R
x
that the code applies after detecting the stitching disturbance,
M
QEC
xy
= P
C,x
R
x
S
xy
P
C,y
, (19)
and the conjecture is that M
QEC
xy
can be full rank, its determinant-one polar unitary U
xy
=
polar(M
QEC
xy
)/ det(···)
1/3
then serving as a candidate SU(3) link.
13
14.2 The test on the real code
The recovery R
x
is not free: it is fixed by the code as the minimum-weight operator carrying the
leakage syndrome. We compute it on the actual [[192, 130, 3]] code. A color/dark leakage is the
relative-sign flip within a matching, a single-bond Z on one cage edge; its real syndrome is the two
octahedral voids bordering that edge (weight two). Exhaustive search over all Z-operators returns a
unique minimum-weight recovery of weight one—the same single bond—touching exactly one color.
There is no competing lower- or equal-weight recovery, and no cross-color syndrome degeneracy:
all six cage edges have distinct syndromes (the distance-three property), so the decoder can never
be steered to a wrong-color recovery. When a leakage hits two colors at once, the minimum-
weight recovery factors into two independent single-bond fixes, one per color, whose product is
block-diagonal in color. Consequently
M
QEC
xy
N(T ), (20)
not full SU(3): R
x
is color-diagonal, so P
C,x
R
x
adds only a diagonal phase to the already-monomial
M
bare
xy
, and the polar unitary remains monomial. The search did find weight-three syndrome-free
cage operators spanning all three colors (one bond from each), but these are products of three
commuting single-bond Z’s—simultaneous diagonal phases, i.e. Cartan/torus elements already in
T , not off-diagonal root generators. The correction cycle does not remove the bond-disjointness
obstruction; it reproduces it.
14.3 Consequence
The QEC recovery cycle stabilizes color—it is exactly the weight-one correctable structure of Sec-
tion 9—but it does not supply the off-diagonal direction. This is a tested-negative result, and it
sharpens the picture: the search space of code-internal spatial operations has been checked, bare
and dressed alike, and all of it lands in N (T). The root generators must come from elsewhere. The
remaining channel, the temporal one, is the subject of the next section.
15 Internal sector IX: temporal resolution—the root direction
from cage evolution
15.1 The temporal link is color-mixing
A four-dimensional Wilson formulation carries, in addition to the spatial links, temporal links
connecting a defect to itself at the next time step, and a temporal link is not a bond operator: it
is the time-evolution transporter of the color qutrit,
U
t
= exp
iH
cage
dt
H
C
. (21)
Its generator is the cage Hamiltonian projected onto the color sector, which by the factorization (9)
(with J
2
acting as +2 on the symmetric sector, absorbed into Γ) is
H
cage
H
C
= J
3
I
3
=
0 1 1
1 0 1
1 1 0
, (22)
off-diagonal in the color index. The temporal link is therefore color-mixing—non-monomial for
any dt = 0 (off-diagonal weight 0.29 at dt = 0.3, 0.665 at dt = 1.0), connecting color i to color j
14
directly. This is exactly the operation the spatial channels could not perform, and it is available
here because temporal transport is generated by the cage kinetic term, which mixes the matchings,
rather than by a single-bond operator, which cannot. Bond-disjointness constrains the spacelike
bonds; it does not constrain time evolution.
15.2 The temporal generator closes su(3)
A single off-diagonal generator suffices. The operator J
3
I
3
is the S
3
-symmetric combination of
the three root directions, and together with the two Cartan generators λ
3
, λ
8
supplied by the bond
phases it closes under commutation to the full algebra. Computing the Lie closure,
λ
3
, λ
8
= 2, λ
3
, λ
8
, J
3
I
3
Lie
= 8 = dim su(3), (23)
(a full Gell-Mann basis returns 8 by the same routine, confirming the count). The commutators
[λ
3
, J
3
I
3
] and [λ
8
, J
3
I
3
] generate the remaining root directions that no spatial channel reached.
The four-dimensional construction therefore attains the entire color algebra, and the off-diagonal
generator is not introduced by hand: it is the cage’s own temporal evolution, the same operator
(J
3
I
3
) whose symmetric sector defined the qutrit. Table 3 collects the logical color action of
every channel.
Channel Logical action on the qutrit Reached
Spatial migration color permutation S
3
(Weyl)
Bond phases (Cartan) diagonal phase T (torus)
Single-bond recovery diagonal T (torus)
Anchor (round trip) identity {1}
Spatial combined monomial N(T ) = T S
3
Temporal link off-diagonal J
3
I
3
root direction
Temporal + torus Lie closure full su(3)
Table 3: Logical color action of each vacuum channel. The spatial channels are color-diagonal and
combine to N(T ); the temporal link is off-diagonal and closes su(3) with the torus.
16 Internal sector X: dynamical footholds
With the algebra and the well-posed Hamiltonian in hand, two bounded dynamical results follow
on the triangular plaquettes; both are leading-order strong-coupling statements, not continuum
results.
16.1 Strong-coupling confinement
The leading strong-coupling term of a Wilson loop [3] tiles the minimal surface the loop bounds
with plaquettes, each contributing one power of the fundamental character coefficient; the character-
expansion technique and the resulting area law are standard [10]. A side-L triangular Wilson loop
in an FCC {111} plane (a triangular lattice of FCC bonds) bounds a minimal surface of exactly
L
2
triangles, with perimeter 3L bonds, so
W f
A
min
= f
L
2
, (24)
15
which decays with the enclosed area, not the perimeter: an area law, hence confinement. The
per-plaquette factor is f = c
f
/N, with c
f
(β) the SU(3) fundamental character coefficient and the
extra 1/N from the link integrations in the tiling. Computing c
f
(β) by Haar Monte Carlo gives
c
f
(β) β/6 at small β (not β/18, which is the per-plaquette factor f after the link integration),
so the string tension is
σa
2
= ln f ln
2N
2
β
= ln
18
β
, β 0, (25)
a finite, positive tension per triangle (Figure 6). The leading strong-coupling expansion confines
on the triangular geometry with no obstruction. It is a small-β statement: by β 6, the regime
that would correspond to coarse-grained continuum QCD, the per-plaquette factor exceeds unity
and the expansion breaks down, so the weak-coupling/coarse-grained regime is not accessible this
way.
0 1 2 3 4 5 6 7 8
¯
0
1
2
3
4
string tension ¾a
2
expansion
breaks down
small-¯: ln(18)
¾a
2
= ¡ln(c
f
=N) (Haar MC)
Figure 6: Leading strong-coupling string tension on the triangular plaquettes (Haar Monte Carlo),
with the small-β form ln(18). The expansion breaks down near β 6; confinement here is the
strong-coupling statement only.
16.2 The single-plaquette spectrum
The smallest gauge-invariant sector is a single plaquette, whose Kogut-Susskind Hamiltonian [4]
on the gauge-invariant (class-function) sector is
H = H
E
λH
B
, H
E
|r = C
2
(r)|r, H
B
=
1
2
(ˆχ
3
+ ˆχ
3
), (26)
with H
E
the quadratic Casimir (diagonal in irreps) and H
B
multiplication by the fundamental
character (off-diagonal via the SU(3) fusion rules), truncated by a Casimir cutoff. The spectrum is
real (Hermitian H
B
) and converges fast: the gap stabilizes by cutoff C
2
10 (dimension 15). In
the strong-coupling limit λ 0 the ground state is the trivial irrep (C
2
= 0, the confined vacuum)
and the gap is exactly
= C
2
(3) =
4
3
, (27)
the single-plaquette fundamental Casimir gap, with no tuning; as λ grows the ground state spreads
over higher irreps and the levels reorganize in the known crossover (Figure 7). The full H
E
+ H
B
is thus a sensible operator with this gap. This is the universal single-plaquette problem—at one
16
plaquette the holonomy is a single SU(3) element regardless of plaquette shape—so it confirms the
Hamiltonian but does not by itself probe triangular-specific physics or continuum confinement, and
the value 4/3 is the fundamental Casimir, not a derivation of the physical QCD mass gap.
0 1 2 3 4 5 6 7 8
magnetic coupling ¸ (weaker gauge coupling ! )
−14
−12
−10
−8
−6
−4
−2
0
2
energy
gap = C
2
(3) =
4
3
E
0
E
1
E
2
Figure 7: Single-plaquette SU(3) Kogut-Susskind spectrum versus magnetic coupling λ. The strong-
coupling gap is C
2
(3) = 4/3; the ground state spreads over higher irreps as λ grows.
17 Unification: one code complex, two sectors
The two sectors are protected by one FCC code complex through one mechanism. The external
charge of Eq. (5) is the vertex Gauss-law flux; the internal triality of Eq. (16) is the Z
3
refinement
of the same vertex condition; and the 513-dimensional non-abelian Gauss law has that triality as
its center-reduction. The binary distance-three code and the Z
3
triality code are not literally the
same stabilizer code—the octahedral void operators do not survive the Z
3
lift (Section 11)—but
they are two realizations on the same FCC chain complex, the binary one carrying the external
Z
2
structure and the Z
3
one the internal triality. The chain is a single Gauss law refined in three
steps,
Z (integer charge) Z
3
(triality) SU(3) (color), (28)
each the next refinement of the vertex condition, all on the same complex. Both sectors share
the distance-three protection: a single-bond disturbance is the correctable error that protects the
migrating charge (external) and the color/dark leakage (internal) alike. And the octahedral baryon
operator, which carries the external sector at the Z
2
level, is subsumed as a boundary by the
triangular triality code that carries the internal sector (Section 11). One code complex, one error-
correcting mechanism, one Gauss law refined from Z to Z
3
to SU(3): charge and color are two
readings of a single protected structure.
18 Scope and status
The results reached are definite and should be stated with their boundaries, and the boundaries
are sharp enough to list in three classes.
17
Reached. In the external sector: the winding charge is a Gauss-law super-selection invariant giving
the baryon ladder; the orientation is explicitly not protected; the migrating defect obeys E = γE
rest
;
migration is code-corrected transport. In the internal sector: the color qutrit, the S
3
Weyl group,
the rank-two torus, and the selection of A
2
=
su(3) are derived; the qutrit is code-protected
(conditionally on a weak anchor or fast correction); the FCC complex supports a genuine Z
3
triality code and a well-posed SU(3) Hamiltonian; the bare spatial connection is N (T )-valued by
bond-disjointness; and the temporal cage evolution supplies the off-diagonal root generators and
closes the full su(3). Leading strong-coupling confinement holds on the triangular plaquettes, with
single-plaquette fundamental Casimir gap C
2
(3) = 4/3.
Tested and negative. The QEC correction cycle was tested as a possible source of the off-diagonal
SU(3) direction, M
QEC
xy
= P
C,x
R
x
S
xy
P
C,y
. On the real code the minimum-weight recovery R
x
is
single-bond and same-color, hence color-diagonal, so M
QEC
xy
N(T ) and the dressed transporter
does not generate root directions (Section 14). The QEC cycle stabilizes color but does not remove
the spatial obstruction; the root direction comes from the temporal link, not from recovery.
Not reached. The continuum dynamics. That the algebra closes and the Hamiltonian is well-
posed does not establish continuum confinement, asymptotic freedom, or the value of Λ
QCD
; the
strong-coupling area law and the single-plaquette gap are footholds, and the expansion delivering
them breaks down at the couplings corresponding to the coarse-grained regime. The correction-
cycle generator L
corr
, the loop stiffness κ
γ
, and the effective coupling g
eff
are not derived from the
microscopic vacuum dynamics. The E = γE
rest
relation is exact for the one-dimensional soliton and
inherited by the three-dimensional defect by covariance, not proven directly in three dimensions.
The closure of su(3) is at the level of the color algebra; whether the temporal-link construction
is fully consistent with the baryon and triality protection of the spatial code, and whether the
four-dimensional theory reproduces QCD under coarse-graining over the physical FCC vacuum,
are open.
The central claim, stated with these boundaries, is therefore the following: one FCC code complex
protects baryon charge and color, the geometry selects the qutrit and triality, the bare spatial
transport is obstructed to N(T ), the QEC recovery cycle does not remove that obstruction, and
the temporal cage Hamiltonian supplies the missing root directions and closes su(3). The full
color algebra appears only when temporal qutrit evolution is included. The boundary this paper
draws is between that structure—two protected sectors, a closed color algebra, and a well-posed
kinematics—which is reached, and the confining continuum dynamics, which are not.
Data availability
The Python scripts reproducing every result are available as a single archive, ssm strong scripts.zip,
downloadable at https://github.com/raghu91302/ssmtheory/raw/main/ssm _strong_script
s.zip. Requiring NumPy (some also SciPy/Matplotlib), they cover: the orientation cost and code
build of Sections 23; the baryon spectrum and kinematics of Sections 45; the migration string of
Section 6; the qutrit emergence, Weyl/rank-two selection, and QEC stabilization of Sections 79;
the gauge-ready geometry, Z
3
triality code, and SU(3) Gauss law and plaquette of Sections 10
12; the spatial N (T ) obstruction, the tested-negative QEC-dressed recovery (the minimum-weight
recovery computation of Section 14), and the temporal su(3) closure of Sections 1315; and the
strong-coupling and single-plaquette computations of Section 16. No other data were generated or
analyzed in this study.
18
References
[1] 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 (2026) 100423.
doi:10.1016/j.physo.2026.100423.
[2] R. Kulkarni, A face-centered-cubic stabilizer code for the vacuum lattice, arXiv:2603.20294
[quant-ph] (2026). arXiv:2603.20294.
[3] K. G. Wilson, Confinement of quarks, Phys. Rev. D 10 (1974) 2445.
doi:10.1103/PhysRevD.10.2445.
[4] J. B. Kogut and L. Susskind, Hamiltonian formulation of Wilson’s lattice gauge theories, Phys.
Rev. D 11 (1975) 395. doi:10.1103/PhysRevD.11.395.
[5] A. R. Calderbank and P. W. Shor, Good quantum error-correcting codes exist, Phys. Rev. A
54 (1996) 1098. doi:10.1103/PhysRevA.54.1098.
[6] A. M. Steane, Error correcting codes in quantum theory, Phys. Rev. Lett. 77 (1996) 793.
doi:10.1103/PhysRevLett.77.793.
[7] A. Yu. Kitaev, Fault-tolerant quantum computation by anyons, Ann. Phys. 303 (2003) 2.
doi:10.1016/S0003-4916(02)00018-0; arXiv:quant-ph/9707021.
[8] S. Chandrasekharan and U.-J. Wiese, Quantum link models: a discrete approach to gauge
theories, Nucl. Phys. B 492 (1997) 455. doi:10.1016/S0550-3213(97)80041-7; arXiv:hep-
lat/9609042.
[9] H. Georgi, Lie Algebras in Particle Physics, 2nd ed., Westview Press, Boulder, 1999.
[10] I. Montvay and G. unster, Quantum Fields on a Lattice, Cambridge University Press, Cam-
bridge, 1994.
[11] R. Rajaraman, Solitons and Instantons, North-Holland, Amsterdam, 1982.
[12] J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, 3rd ed., Springer,
New York, 1999.
19