Thermodynamic Emergence:
Deriving the Cuboctahedral Vacuum from
Geometric Saturation and Topological
Ground States
Raghu Kulkarni
∗
Independent Researcher, Calabasas, CA
March 2, 2026
Abstract
We explore the thermodynamic and topological origins of the foundational space-
time geometry in the Selection-Stitch Model (SSM) [1, 2]. Rather than assuming a
background manifold, we posit that the continuous vacuum emerges from a discrete
quantum tensor network evolving to minimize its free energy while resolving the
geometric frustration of the early universe. In this framework, the minimal-energy
ground state of the Hamiltonian (T → 0) is identified as the two-dimensional K = 6
hexagonal sheet. As the tetrahedral foam (K = 4) of Cosmic Inflation expands, it
thermodynamically seeks to saturate its open deficit angles. We derive a topologi-
cal free energy functional directly from a microscopic partition function, utilizing a
controlled cluster expansion bounded by the network’s hard-sphere exclusion limit.
This yields a physically motivated kinetic rate equation via non-conserved order
parameter dynamics. Drawing on the Kepler Conjecture [6], the absolute maxi-
mum density for 3D Euclidean packing caps the coordination at K = 12. While
this density is shared by Face-Centered Cubic (FCC) and Hexagonal Close-Packed
(HCP) lattices, the FCC lattice is selected during the Alder transition. Under an
adiabatic cosmological quench (Γ ≫ H ), its isotropic point-group symmetry maxi-
mizes the phononic density of states, making its vibrational entropy advantage [7, 8]
deterministic in the thermodynamic limit. Finally, we demonstrate that the FCC
Cuboctahedron is the unique low-energy variational minimum preserving the 2D
K = 6 zero-stress Euler state across its 3D bulk.
Contents
1 Introduction: The Geometry Selection Problem 3
2 The Thermodynamic Free Energy Functional 3
2.1 The Euler Ground State (⟨k⟩ = 6) and Pre-Geometric Thermodynamics . . 3
2.2 Resolving 3D Geometric Frustration . . . . . . . . . . . . . . . . . . . . . . 4
2.3 Microscopic Derivation of the Free Energy . . . . . . . . . . . . . . . . . . 5
2.4 Mean-Field Fluctuations and Stability . . . . . . . . . . . . . . . . . . . . 6
∗
raghu@idrive.com
1
3 The Avoidance of Topological Glass and Kibble-Zurek 6
4 Breaking the Degeneracy: The Perfect Geometric Extension 7
4.1 Vibrational Entropy in the Thermodynamic Limit . . . . . . . . . . . . . . 8
4.2 A Discrete Variational Principle for Geometric Extension . . . . . . . . . . 8
5 The Universality of Geometric Ground States 9
5.1 Falsifiable Consequences: The Silver Ratio Sequence . . . . . . . . . . . . . 9
6 Conclusion 10
2
1 Introduction: The Geometry Selection Problem
A foundational critique of any discrete model of quantum gravity is the problem of ge-
ometric selection: if spacetime is not a smooth, pre-existing continuous manifold, but
an emergent network of discrete quantum information, what mechanism determines its
ultimate macroscopic structure?
While approaches like Causal Dynamical Triangulations (CDT) [3] and Spin Foam
models explore the quantum superposition of geometries, and tensor networks (like MERA)
[4] establish bulk-boundary entanglement, the Selection-Stitch Model (SSM) focuses on
the strict thermodynamic selection of a singular, deterministic crystalline ground state.
Within the geometric framework of the SSM, the present-day vacuum is mathematically
modeled as a continuous Face-Centered Cubic (FCC) lattice, defined by the Cuboctahe-
dron unit cell (K = 12). A natural question arises: Why this specific topology? Why
did the universe not freeze into an amorphous topological glass, a simple cubic grid, or
remain as a fractured Tetrahedral Foam (K = 4)?
We propose that the selection of the K = 12 FCC vacuum is a natural consequence
of Thermodynamic Efficiency and Topological Saturation. In the pre-geometric
phase of the early universe, quantum geometric fluctuations were stochastically driven by
extreme topological frustration. As the universe expanded, these fundamental nodes of
information evolved to minimize their free energy, structurally packing together to form
the most energetically stable spacetime metric permissible by strict topological boundary
conditions.
Throughout this work, we establish a consistent notational convention: we denote
the target macroscopic lattice coordination state as K (e.g., the K = 12 geometry), the
discrete per-vertex coordination variable as k
i
, and its statistical mean-field expectation
value as ⟨k⟩.
2 The Thermodynamic Free Energy Functional
2.1 The Euler Ground State (⟨k⟩ = 6) and Pre-Geometric Ther-
modynamics
To understand the 3D vacuum, we must first establish its absolute minimal-energy topo-
logical ground state. In the extreme low-energy limit (T → 0), the network seeks to close
minimal entanglement loops (triangles) without inducing stress. The global topology of
a 2D triangulated surface is dictated by the Euler characteristic χ = V − E + F. For
a closed manifold constructed entirely of triangles, 3F = 2E, yielding the macroscopic
curvature constraint:
χ = V
1 −
⟨k⟩
6
(1)
For an infinite planar triangulation where boundary terms vanish, the global topology
dictates χ = 0. Transverse (out-of-plane) dimensionality carries a volumetric entropic
penalty, requiring additional energetic action to maintain orthogonal degrees of freedom.
Each transverse degree of freedom requires a minimum excitation energy of order k
B
T to
sustain, penalizing configurations that extend beyond the minimal dimensionality at low
temperatures.
Applying thermodynamic concepts to a pre-geometric network requires a formal def-
inition of the thermal ensemble. In the Euclidean path integral formulation of discrete
3
quantum gravity, the partition function over all valid graph topologies is evaluated as
Z =
P
exp(−S
R
), where S
R
is the Euclidean Regge action. While a standard Boltz-
mann ensemble samples state configurations within a fixed spatial geometry, here the
graph topology itself is the dynamical variable. In the mean-field approximation, sum-
ming over these local graph topologies maps mathematically to a lattice gas of structural
bonds [5], where the formal network “temperature” k
B
T is rigorously defined as being
inversely proportional to the gravitational coupling. Physically, this temperature acts as
the magnitude of topological quantum fluctuations (the edge-rewiring rate of the tensor
network).
During the onset of inflation, extreme vacuum energy drives rapid edge-swapping (T →
∞). As the universe expands and the inflationary potential dilutes, the effective coupling
strength rises, acting as the inherent cooling mechanism that quenches the network toward
T → 0. As the system cools, the Hamiltonian minimizes when the lattice achieves its
densest 2D packing: the flat hexagonal plane (⟨k⟩ = 6, forcing χ = 0). Therefore, the
foundational thermodynamic ground state of the universe naturally presents as a 2D
boundary sheet.
2.2 Resolving 3D Geometric Frustration
As established in our companion derivation of the Early Universe [1], the epoch of Cosmic
Inflation was mechanically driven by the volumetric projection of this 2D state into a
3D Tetrahedral Foam (K = 4). Regular tetrahedra cannot smoothly tile 3D Euclidean
space; when five regular tetrahedra are packed around a single shared edge, they consume
5 × 70.528
◦
= 352.64
◦
of the available radial volume. This leaves open an unavoidable
Regge deficit angle of δ ≈ 7.36
◦
(0.128 rad), as illustrated in Figure 1.
Figure 1: The geometric origin of Cosmic Inflation. Packing five regular tetrahedra
(K = 4) around a shared central edge leaves a strict topological void (the Regge deficit
angle of δ ≈ 0.128 rad). The vacuum thermodynamically seeks to close this energetic gap
by transitioning to the fully saturated K = 12 continuous geometry.
4
This unclosed metric gap acts as a compressed geometric spring. To exit inflation and
reach a stable continuum state, the tensor network seeks to perfectly close these deficit
angles by maximizing the local coordination per unit volume.
2.3 Microscopic Derivation of the Free Energy
Rather than asserting a phenomenological functional, we derive the topological free energy
F (⟨k⟩) directly from the microscopic partition function of the network. We define an
effective mean-field Hamiltonian driven by pairwise bonding:
H
eff
= −
α
2
X
i
k
i
(2)
A complete quantum gravity Hamiltonian includes an infinite series of complex, non-
linear multi-body curvature penalties. However, we can rigorously bound the leading
neglected term to justify this truncation. In a Regge framework, a three-body curvature
penalty between adjacent hinges scales with the square of the deficit angle, O(δ
2
). With
δ ≈ 0.128 rad, the leading non-linear perturbation per interaction loop is strongly sup-
pressed (δ
2
≈ 0.016). Because the geometric exclusion limit strictly caps the number of
adjacent triangles sharing any single bond (e.g., exactly 4 in the FCC lattice), the total
non-linear curvature correction per bond scales as ∼ 4δ
2
≈ 0.064. This bounded 6.4%
perturbative correction shifts the exact value of the critical transition temperature, but
is mathematically insufficient to destabilize the deep global minimum of the topological
ground state, safely justifying the linearized approximation.
The fundamental energy scale α is derived from the standard Regge action, which
assigns a spatial hinge h an energy proportional to
1
G
δ
h
A
h
. At the Planck scale (A
h
∼ l
2
P
),
the area and Newton’s constant cancel in natural units, yielding α ≈ δ · E
P
≈ 0.128E
P
.
We treat this specific magnitude as a strict self-consistency condition: the fundamental
quantum of activation energy driving the phase transition must identically match the
specific quantum of geometric frustration generated by the K = 4 tetrahedral foam. The
factor of 1/2 corrects for pairwise double-counting.
The canonical partition function is Z =
P
{k
i
}
Ω(k
i
)e
−H
ef f
/k
B
T
. Drawing on Thomas
Hales’ formal proof of the Kepler Conjecture [6], the absolute maximum density for hard
spheres rigidly caps the geometric slots at K = 12. The combinatorial multiplicity (phase
space volume Ω) of assigning k identical bonds into 12 available slots is the binomial
coefficient Ω(k) =
12
k
.
In the mean-field limit, applying Stirling’s approximation, the macroscopic entropy of
this topological mixing is S ≈ k
B
[12 ln 12 − ⟨k⟩ln⟨k⟩ − (12 − ⟨k⟩) ln(12 − ⟨k⟩)]. Evaluat-
ing F = ⟨H
eff
⟩ − T S yields the macroscopic functional:
F (⟨k⟩) = −α
⟨k⟩
2
+ k
B
T [⟨k⟩ln⟨k⟩ + (12 − ⟨k⟩) ln(12 − ⟨k⟩)] (3)
(dropping the constant 12 ln 12 term). This establishes that the excluded-volume entropy
penalty is a rigorous combinatorial consequence of the finite phase space of a bounded
topological network.
5
2.4 Mean-Field Fluctuations and Stability
Minimizing this functional (∂F/∂⟨k⟩ = 0) reveals the equilibrium state of the mean-field
coordination:
−
α
2
+ k
B
T ln
⟨k⟩
12 − ⟨k⟩
= 0 (4)
Solving for ⟨k⟩ yields the classic sigmoidal saturation curve:
⟨k⟩
eq
=
12
1 + e
−α/2k
B
T
(5)
To confirm this stationary point is a stable global minimum, we evaluate the second
derivative:
∂
2
F
∂⟨k⟩
2
= k
B
T
1
⟨k⟩
+
1
12 − ⟨k⟩
=
12k
B
T
⟨k⟩(12 − ⟨k⟩)
> 0 (6)
Because this is strictly positive for all 0 < ⟨k⟩ < 12, the equilibrium is definitively stable.
The variance of the topological fluctuations is given by:
⟨(δk)
2
⟩ =
k
B
T
∂
2
F
∂⟨k⟩
2
=
⟨k⟩(12 − ⟨k⟩)
12
(7)
This exactly recovers the variance of a binomial distribution Np(1 − p) where N = 12
and p = ⟨k⟩/12. It explicitly demonstrates that at absolute zero (T = 0), the topological
variance vanishes, e
−∞
→ 0, and the vacuum strictly saturates at the Kepler limit of
K = 12.
3 The Avoidance of Topological Glass and Kibble-
Zurek
A related objection considers Random Close Packing (RCP), a state where discrete nodes
form a disordered topological glass with a localized coordination of K ≈ 12, but a signif-
icantly lower macroscopic density (ρ ≈ 0.64). Why did the vacuum not freeze into this
amorphous, glassy state?
In condensed matter physics, systems fall into glassy local minima when quenched too
rapidly. However, the vacuum transition is governed by the Alder Transition [9]. In
statistical mechanics, a gas of hard spheres will spontaneously undergo a first-order phase
transition into a crystal, not to minimize interaction energy, but purely to maximize its
global entropy once a critical packing fraction is reached.
How does macroscopic cosmic expansion drive a transition dependent on compression?
We must distinguish between metric volume and topological density. In a graph, the
packing fraction ϕ = ⟨k⟩/12 represents the ratio of realized edges to maximum possible
edges. As inflation injects K = 4 tetrahedra into a causal patch, the 0.128 rad deficit
acts as an energetic penalty. To minimize this action, the network relentlessly stitches
new edges to close the angular gaps. This internal topological knitting actively drives the
local graph density ϕ upward, even as the macroscopic boundary expands.
We derive the kinetic relaxation of this local patch using non-conserved order parame-
ter dynamics (Model A), where the rate of change is proportional to the negative gradient
of the free energy:
d⟨k⟩
dt
= −Γ
∂F
∂⟨k⟩
= Γ
α
2
− k
B
T ln
⟨k⟩
12 − ⟨k⟩
(8)
6
where Γ is the transition mobility. The Alder transition famously occurs at a critical
packing fraction of ϕ ≈ 0.494. As the geometric “spring” of the K = 4 foam drives the
local packing fraction toward the K = 12 Kepler limit (ϕ = π/
√
18 ≈ 0.740), it inevitably
sweeps across this critical density threshold.
To avoid freezing into a network of topological defects via the Kibble-Zurek mechanism
[10, 11], the local topological relaxation time (1/Γ) must be significantly faster than the
macroscopic Hubble expansion rate (1/H). We provide an order-of-magnitude plausibility
bound for this adiabatic condition. While deriving Γ from a full quantum kinetic the-
ory remains an open problem, dimensional necessity bounds the fundamental topological
relaxation rate near the Planck frequency (Γ ∼ t
−1
P
∼ 10
43
Hz). The Hubble expansion
rate during inflation is bounded by the GUT scale, typically H ∼ 10
37
Hz. Therefore, the
ratio Γ/H ∼ 10
6
≫ 1. This massive timescale separation establishes the plausibility of
an adiabatic quench, allowing crystalline domains to anneal over super-horizon scales.
4 Breaking the Degeneracy: The Perfect Geometric
Extension
Critically, two distinct lattice configurations share this exact maximum density and K =
12 coordination limit: the Face-Centered Cubic (FCC) and the Hexagonal Close-Packed
(HCP) lattices. A purely density-driven model cannot distinguish between them. How-
ever, thermodynamics offers a selection mechanism via Entropy Maximization and Sym-
metry Breaking.
Figure 2: Breaking the K = 12 topological degeneracy. While both geometries satisfy
the Kepler limit for maximum 3D density, the FCC configuration (left) preserves central
isotropy (O
h
). The HCP configuration (right) forces a preferred vertical axis (D
6h
). The
vacuum selects FCC to maximize vibrational entropy.
7
4.1 Vibrational Entropy in the Thermodynamic Limit
As illustrated in Figure 2, the topological difference lies entirely in the stacking sequence:
the FCC lattice is structurally isotropic due to its cyclic (ABCABC) sequence, while the
HCP lattice possesses a preferred vertical axis due to alternating (ABAB) stacking.
In the statistical mechanics of hard spheres, the FCC lattice possesses a higher vi-
brational (phonon) entropy than the HCP lattice [7, 8]. Applying this classical result to
a quantum tensor network is structurally robust because the phononic density of states
is a direct group-theoretic consequence of the local potential well’s symmetry. Because
the FCC geometry possesses a higher point-group symmetry (O
h
) than HCP (D
6h
), its
perfectly isotropic potential natively accommodates a wider, more symmetric spectrum
of transverse zero-point fluctuations. This geometrically enforces a higher vibrational
entropy regardless of the specific microscopic bond interactions.
While this entropic advantage is famously microscopic locally (on the order of 10
−3
k
B
per particle), phase transitions in the early universe are collective, global phenomena.
Because Γ ≫ H, the crystallization front propagates adiabatically, allowing the causal
patch to maintain thermal equilibrium where the total free energy difference scales as
∆F = −N(T ∆S
vib
). In the thermodynamic limit (N → ∞), the probability of the
universe falling into the HCP state over the FCC state is given by the Boltzmann weight:
P
HCP
P
F CC
= exp
−
N∆S
vib
k
B
→ 0 (9)
The massive scale of the universe transforms a microscopic local advantage into a deter-
ministic outcome in the thermodynamic limit. Random stacking faults are exponentially
suppressed, strictly driving the selection of the isotropic FCC state.
4.2 A Discrete Variational Principle for Geometric Extension
Beyond simple density and entropy, the selection of the FCC lattice satisfies a discrete
topological variational principle. When transitioning from the 2D boundary to the 3D
bulk, the Hamiltonian seeks to minimize spatial shear by matching the zero-stress 2D
Euler ground state (⟨k⟩ = 6) as closely as possible along any continuous cross-section.
To establish uniqueness, we restrict our search to the set of K = 12 close-packed
Euclidean vertex figures. Distorted or non-regular geometries inherently require unequal
bond lengths, invoking a severe harmonic elastic energy penalty that scales as ∆E ∼
κ(∆l/l)
2
per bond. In the high-tension environment of the expanding early universe,
this immediately removes distorted geometries from the low-energy ground state. The
minimum-energy candidates are thus strictly restricted to the three maximally symmetric
contact graphs for 12 uniform spheres: the Icosahedron, the HCP Twisted Cuboctahedron,
and the FCC Cuboctahedron.
The Icosahedron entirely lacks regular planar equators, geometrically forcing k
i
= 5
across all skewed intersecting planes. The HCP figure contains exactly one regular k
i
= 6
equatorial plane, breaking spatial regularity along its vertical axis. The FCC Cuboctahe-
dron is mathematically unique in possessing exactly four regular k
i
= 6 equatorial cross-
sections (the {111} crystallographic planes), perfectly satisfying the discrete ⟨k
2D
⟩ = 6
minimum-energy condition across all spatial axes simultaneously.
8
5 The Universality of Geometric Ground States
A frequent critique of discrete spacetime models is the application of classical sphere-
packing logic to purely quantum objects (vacuum nodes). However, the mapping of ran-
dom tensor networks to lattice models (such as the Ising or Potts models) guarantees that
geometric phase transitions in graphs fall into established universality classes [5]. This
ensures that the collective topological phase transition is independent of the underlying
microscopic quantum gravity action.
We observe distinct physical analogies in condensed matter physics where purely quan-
tum mechanical constituents spontaneously self-organize into classical geometric lattices
to minimize their free energy:
• Abrikosov Lattices: Quantized magnetic flux tubes naturally self-assemble into
a perfect 2D Hexagonal Lattice [12], achieving the absolute mathematical limit of
2D circle packing (ρ = π/
√
12).
• Coulomb Crystals: Laser-cooled ions geometrically arrange themselves into con-
centric shells that mimic classical close-packing FCC structures [13].
• Wigner Crystallization: A purely quantum electron gas “freezes” into a rigid,
classical Body-Centered Cubic (BCC) lattice to minimize long-range Coulomb re-
pulsive potential [14].
Unlike electrons optimizing a long-range Coulomb field, the emergent vacuum tensor
network is governed by strictly local topological saturation. This optimization criterion
points directly to the Kepler limit, suggesting that the K = 12 continuous geometry is
the natural topological ground state of an emergent spacetime.
5.1 Falsifiable Consequences: The Silver Ratio Sequence
A crucial requirement of this thermodynamic derivation is falsifiability. If the K = 12
vacuum is an arbitrary phenomenological assumption, its geometry carries no mandated
physical signatures. However, if the vacuum was thermodynamically selected specifically
as the FCC lattice to maximize vibrational phonon entropy, the acoustic overtones of this
specific lattice must be imprinted on the universe.
Because the fundamental FCC geometric constant is
√
2, the dynamical matrix of
the isolated K = 12 cuboctahedral cluster produces conjugate vibrational eigenvalues of
ω
2
= 1 ± 1/
√
2. In this discrete framework, the primordial baryon-photon fluid does not
oscillate in a continuous vacuum, but rather couples directly to the discrete density per-
turbations (phonons) of the underlying lattice. The acoustic normal modes of the spatial
cuboctahedron thus physically dictate the overtones of the primordial sound horizon.
As detailed in our companion work regarding acoustic peaks, this exact algebraic
structure simultaneously predicts the locations of the first three acoustic peaks, mandating
a rigid sequence: l
1
= 221, l
2
= 221(1 +
√
2) ≈ 534, and l
3
= 221 × 2/
q
1 − 1/
√
2 ≈
817. Standard ΛCDM extracts this exact same spectrum by continuously tuning multiple
physical parameters, including the spatial curvature and the baryon and dark matter
densities (Ω
k
, Ω
b
, Ω
c
). Conversely, the SSM dictates this specific sequence as a zero-
parameter geometric prior.
While one matching ratio might be dismissed as coincidence, the full sequence locks
into the Planck 2018 (PR3) baseline observations [15]: the predicted l
1
= 221 matches
9
the observed l
1
= 220.9 ± 0.5 (0.05% error); the predicted l
2
≈ 534 matches the observed
l
2
≈ 536 (0.4% error); and the predicted l
3
≈ 817 matches the observed l
3
≈ 810 (0.9%
error). To definitively distinguish the discrete SSM geometry from fitted continuous fluid
models, this framework predicts that this rigid algebraic spacing will persist unfettered
into the higher-order E-mode polarization peaks, providing a rigid, testable diagnostic
template free from the parameter degeneracy of ΛCDM.
6 Conclusion
The Cuboctahedral geometry (K = 12) of the Selection-Stitch Model is not an arbitrary
phenomenological assumption. It emerges naturally as the geometric result of fundamental
thermodynamic drives resolving the instability of the early universe:
1. Topological Saturation: The frustrated vacuum seeks maximum density to close
the Regge deficit angles of inflation, structurally bounded by the exact mathematical
limit of 3D Euclidean space (π/
√
18). We derive this rigorously from a truncated
microscopic partition function, formalized via a free energy functional governed by
non-conserved kinetic dynamics.
2. Entropy Maximization: Driven upward to the critical density packing fraction
of 0.494 by the internal topological knitting of inflation, the vacuum avoids a glassy
state via the Alder transition. Assuming an adiabatic Kibble-Zurek quench (Γ ≫
H), the higher point-group symmetry of the FCC stacking generates a broader
phononic density of states, deterministically breaking the structural degeneracy with
the HCP lattice in the thermodynamic limit.
3. Variational Extension: Verified by geometric exhaustion of the regular close-
packed vertex figures, the selected FCC Cuboctahedron acts as a strict topological
minimum, uniquely preserving the 2D ⟨k⟩ = 6 zero-stress Euler ground state across
the 3D volumetric bulk.
By framing the origin of the continuous vacuum as the thermodynamic emergence
of a maximally dense geometric ground state, we eliminate the need for arbitrary fine-
tuning. The Face-Centered Cubic vacuum appears to be the most efficient topological
configuration that three-dimensional quantum information can occupy.
References
[1] R. Kulkarni, “Geometric Phase Transitions in a Discrete Vacuum: Deriving Cosmic
Flatness, Inflation, and Reheating from Tensor Network Topology.” Under Review.
https://doi.org/10.5281/zenodo.18727238 (2026).
[2] R. Kulkarni, “Constructive Verification of K = 12 Lattice Saturation: Exploring
Kinematic Consistency in the Selection-Stitch Model.” Under Review. https://doi.
org/10.5281/zenodo.18294925 (2026).
[3] J. Ambjørn, J. Jurkiewicz, and R. Loll, “Emergence of a 4D World from Causal
Dynamical Triangulations,” Physical Review Letters, 93(13), 131301 (2004).
10
[4] B. Swingle, “Entanglement renormalization and holography,” Physical Review D,
86(6), 065007 (2012).
[5] P. Hayden, S. Nezami, X. L. Qi, N. Thomas, M. Walter, and Z. Yang, “Holographic
duality from random tensor networks,” JHEP, 2016(11), 9 (2016).
[6] T. C. Hales, “A proof of the Kepler conjecture,” Annals of Mathematics, 162(3),
1065-1185 (2005).
[7] D. Frenkel and A. J. C. Ladd, “New Monte Carlo method to compute the free energy
of arbitrary solid phases,” Journal of Chemical Physics, 81(7), 3188-3193 (1984).
[8] L. V. Woodcock, “Entropy difference between the face-centred cubic and hexagonal
close-packed crystal structures,” Nature, 385, 141-143 (1997).
[9] B. J. Alder and T. E. Wainwright, “Phase Transition for a Hard Sphere System,”
Journal of Chemical Physics, 27(5), 1208-1209 (1957).
[10] T. W. B. Kibble, “Topology of cosmic domains and strings,” Journal of Physics A:
Mathematical and General, 9(8), 1387 (1976).
[11] W. H. Zurek, “Cosmological experiments in superfluid helium?,” Nature, 317, 505-
508 (1985).
[12] A. A. Abrikosov, “On the Magnetic Properties of Superconductors of the Second
Group,” Soviet Physics JETP, 5, 1174 (1957).
[13] G. Birkl, S. Kassner, and H. Walther, “Multiple-Shell Structures of Laser-Cooled
24
Mg
+
Ions in a Quadrupole Storage Ring,” Nature, 357, 310-313 (1992).
[14] E. Wigner, “On the Interaction of Electrons in Metals,” Physical Review, 46(11),
1002 (1934).
[15] Planck Collaboration, “Planck 2018 results. VI. Cosmological parameters,” Astron.
Astrophys., 641, A6 (2020).
11
