The 1.37 Billion Year Big Bang:
Deriving a Universal Age Gradient and
Co-Aligned Structure Dipoles
from a Single-Origin Vacuum Crystallization
Raghu Kulkarni
SSMTheory Group, IDrive Inc., Calabasas, CA 91302, USA
raghu@idrive.com
March 14, 2026
Abstract
The Standard Cosmological Model (ΛCDM) postulates that primordial reheating oc-
curred simultaneously across a spacelike hypersurface, predicting a strictly isotropic
universe. However, this perfectly symmetric prediction (A = 0) is in > 3σ tension
with the intrinsic hemispherical power asymmetry (A = 0.066 ± 0.021) observed
by the Planck satellite [1]. In this paper, we demonstrate that this anomaly is the
geometric fingerprint of a single-origin discrete vacuum phase transition. In the
Selection-Stitch Model (SSM) [4], the “Big Bang” is a singular nucleation event that
propagates outward as a volumetric crystallization front. We derive the temporal
lag of this front using kinematic interface drag, showing that the tunneling proba-
bility per lattice sweep (p = e
−3
≈ 0.0498) establishes a permanent, scale-invariant
fractional age gradient of ∆t/t ≈ 0.0996 across the observable horizon. Evaluated
at the present epoch, one hemisphere of our observable universe is approximately
1.37 billion years older than the antipodal hemisphere. Because this gradient origi-
nates from a single, unidirectional wavefront, the CMB dipole, high-redshift quasar
evolutionary states [10], and large-scale structure maturity must all be perfectly
co-aligned along a single preferred axis.
1 Introduction: The Symmetry Paradox
Standard ΛCDM cosmology assumes that cosmic inflation was driven by a continuous
scalar field decaying on a spacelike hypersurface of constant field value [2]. By construc-
tion, this mandates that reheating happened everywhere simultaneously. The predicted
age difference between any two antipodal points on the Cosmic Microwave Background
(CMB) horizon is exactly zero. Yet, independent analyses of WMAP and Planck satellite
data have consistently confirmed a significant, intrinsic hemispherical power asymmetry—
a dipole modulation in the temperature fluctuation amplitude across the sky [1,3]. Planck
1
2018 explicitly bounds this dipole at A = 0.066 ± 0.021, excluding the ΛCDM null hy-
pothesis (A = 0) at greater than 3σ significance.
In this paper, we model the Big Bang as a physical phase transition originating from a
singular nucleation event within the discrete spacetime framework of the Selection-Stitch
Model (SSM) [4, 5]. We derive the finite temporal lag of the propagating crystallization
front from first principles, connect it to a macroscopic age gradient, and show that this
gradient quantitatively matches the observed CMB hemispherical asymmetry.
Interactive 3D visualizations. Readers can explore the geometric phase transi-
tions and macroscopic wave propagation discussed in this paper through two interac-
tive WebGL applications:
1. Spacetime Crystallization: The K = 6 → K = 4 → K = 12 topological
relaxation, visually illustrating how the volumetric lock-in fundamentally trails the
2D causal expansion:
https://raghu91302.github.io/ssmtheory/ssm_regge_deficit.html
2. Big Bang Wavefront: A macroscopic simulation of the single-origin crystal-
lization front sweeping across a comoving volume, demonstrating the origin of the
temporal lag and the resulting structure dipoles:
https://raghu91302.github.io/ssmtheory/ssm_bigbang_wavefront.html
2 The Single-Origin Crystallization Front
In the SSM, the early universe undergoes a geometric phase transition from a hot, topo-
logically frustrated tetrahedral foam (K = 4) into a cold, saturated Face-Centered Cubic
continuum (K = 12) [4]. This is not a simultaneous global event. There is exactly one
singular nucleation event. The entire observable universe traces back to this one origin
point propagating outward.
2.1 Two propagation mechanisms
The geometry of the expansion is strictly stratified by the causal limits of the discrete
network:
The 2D holographic boundary. The causal thermal conduit—the boundary sheet of
the tensor network—races outward from the singular origin at the primary lattice veloc-
ity v
front
. This speed is set by the translational structure tensor: v
front
= S
trans
× v
lattice
=
4v
lattice
= c [6]. It provides instantaneous thermal equilibration across the expanding vol-
ume relative to the bulk lock-in, naturally resolving the horizon problem without invoking
a separate continuous inflaton field.
The 3D volumetric lock-in. The actual K = 4 → K = 12 geometric relaxation—the
crystallization that releases the latent heat of reheating—trails behind the causal front.
It cannot proceed instantly; it is limited by the microscopic nucleation probability.
2
2.2 From activation barriers to macroscopic lag
The connection between the microscopic phase transition and the macroscopic cosmolog-
ical age gradient requires an explicit kinematic derivation.
2.2.1 The topological activation barrier (p = e
−3
)
The K = 4 → K = 12 transition at a single lattice site requires rearranging the local
coordination of a tetrahedral void. The minimum topological operation to resolve a K = 4
frustrated state into the saturated K = 12 background requires simultaneously severing
and reconnecting the 3 skew-edge pairs of the tetrahedron. This corresponds strictly to
the geometric crossing number c = 3 [7].
If the fundamental entanglement bond energy of the lattice is ε, the activation barrier
to execute this topological rearrangement is E
a
= 3ε. To evaluate the transition prob-
ability without introducing arbitrary free parameters, we introduce the thermodynamic
ansatz of the SSM: at the critical boundary of the volumetric phase transition (the Big
Bang quench), the ambient thermal energy of the local heat bath equals the latent bond
energy of the lattice (k
B
T
reheat
= ε).
Therefore, the tunneling probability per causal sweep (τ
sweep
∼ t
P
) follows the standard
Arrhenius transition rate:
p = exp
−
E
a
k
B
T
reheat
= exp
−
3ε
ε
= e
−3
≈ 0.04979 (1)
This is mathematically derived from the tetrahedral geometry (c = 3) evaluated exactly
at the critical phase boundary. Because this geometric crystallization is a first-order phase
transition, it proceeds isothermally. The macroscopic release of latent bond energy at the
3D lock-in boundary precisely maintains the local front temperature at T
reheat
as it sweeps
through the unrelaxed bulk (analogous to water remaining at exactly 0
◦
C while freezing).
Consequently, the tunneling probability p = e
−3
remains strictly constant throughout the
entire macroscopic duration of the transition.
2.2.2 The structural interface drag equation
Unlike classical bulk nucleation, this is a single-origin interface. We model the tethered
expansion of the universe as a kinematic interface drag. By the definition of a discrete
kinematic lattice, the absolute maximum speed of information transfer (c) is exactly one
lattice spacing a per fundamental timestep (τ
sweep
). Thus, the 2D causal front acts as
a strict kinematic ceiling, advancing its radial position by exactly x
c
(t) = ct. The 3D
structural lock-in front, x
s
(t), is topologically tethered to this causal front but is retarded
by the tunneling barrier.
We derive this delay rigorously as a Bernoulli process. During each lattice sweep of
duration τ
sweep
= a/c, the structural front attempts to advance by one lattice constant a.
The probability of remaining frustrated is p; the probability of successful lock-in is (1−p).
The expected spatial advancement per sweep is exactly ⟨∆x⟩ = a(1 −p) + 0(p) = a(1 − p).
The macroscopic effective velocity of the structural front is therefore the expected rate of
this progression: v
s
= ⟨∆x⟩/τ
sweep
= c(1 − p).
The physical distance between the causal front and the structural lock-in front grows
linearly over time: x
c
(t) − x
s
(t) = cpt. The accumulated temporal lag ∆t at any given
3
point is simply this spatial separation divided by the causal speed c:
∆t =
x
c
(t) − x
s
(t)
c
= p · t (2)
This Bernoulli derivation explicitly proves why the temporal lag is strictly linear with
respect to the crossing time.
2.3 Deriving the age gradient
Consider an observer situated within the fully crystallized bulk. Let R
0
be the comoving
radius of the observable universe, corresponding to a causal look-back time of t
0
(the
present cosmic age). The comoving diameter of the observable universe is therefore D
obs
=
2R
0
, representing a maximum causal crossing time of:
t
cross
= 2t
0
(3)
Because the wavefront swept past our local comoving volume from one specific direc-
tion, the hemisphere facing the origin crystallized first, while the antipodal hemisphere
crystallized last. Substituting the maximum crossing time into our linear drag formula
(Equation 2), the total accumulated temporal difference between the two poles of the
observable universe is:
∆t = p(2t
0
) (4)
Dividing by the observer’s isotropic baseline age (t
0
) yields the scale-free fractional age
gradient across the sky:
∆t
t
0
= 2p = 2e
−3
≈ 0.0996. (5)
This establishes a rigid mathematical chain:
Crossing number c = 3 → Activation barrier e
−3
→ Kinematic Drag p · t
cross
→ Age
gradient ∆t/t
0
= 2p.
Given the single thermodynamic ansatz (k
B
T
reheat
= ε), every subsequent step is
determined strictly by the tetrahedral geometry and interface kinematics. There are no
additional fitted phenomenological parameters.
3 Baryogenesis as Incomplete Crystallization
Standard cosmology relies on unobserved high-energy mechanisms satisfying the Sakharov
conditions to explain the origin of matter. Within the SSM, the K = 4 → K = 12
phase transition provides a direct geometric mechanism. Baryonic matter consists of
frozen remnants of the pre-crystallization vacuum—isolated K = 4 tetrahedral voids
permanently trapped within the K = 12 bulk [8]. The observed cosmic baryon-to-photon
ratio (η ∼ 6 × 10
−10
) is not an arbitrary free parameter; it is a direct physical measure of
the crystallization efficiency—the volumetric fraction of K = 4 domains that topologically
locked and failed to transition.
4
4 The Duration of the Big Bang
The fractional age gradient (Equation 5) is scale-free and epoch-independent. Evaluating
at the accepted cosmic age (t
0
= 13.80 Gyr) [1]:
t
older
≈ t
0
(1 + p) = 13.80 × 1.04979 ≈ 14.49 Gyr, (6)
t
younger
≈ t
0
(1 − p) = 13.80 × 0.95021 ≈ 13.11 Gyr. (7)
The absolute chronological difference between the two poles of our observable universe
is:
∆t = 2e
−3
× t
0
≈ 1.374 billion years. (8)
In plain language: The 2D causal front (traveling at c) crosses the comoving diameter
of the observable universe in 2t
0
≈ 27.6 Gyr. Because the 3D volumetric crystallization
travels at an effective velocity v
s
= c(1−p), the interface lags the causal front. The differ-
ential lag accumulated strictly during the time the front traverses our specific observable
patch is exactly p × 2t
0
≈ 1.37 Gyr.
This is the literal physical duration of the Big Bang phase transition as experienced
across our comoving volume.
5 Universal Amplitude, Observer-Dependent Direction
The single-origin model yields a profound philosophical and observational consequence.
What standard cosmology terms the “Big Bang” was not a simultaneous global explosion,
but rather a localized crossing event. The crystallization front that ignited our observable
volume is merely a small patch of a vastly larger phase transition. This macroscopic
wavefront originated long before our local “time zero” and continues to sweep through
the unobservable bulk far beyond our cosmic horizon. Our position within this network
is not strictly special; we reside at an arbitrary location that the front happened to pass
through ∼ 13.8 Gyr ago.
Because the fractional gradient (2p ≈ 10%) is a scale-free property of the phase tran-
sition’s trailing temporal lag, it is a universal constant. Every observer, anywhere in
the crystallized universe, will measure a dipole of the exact same fractional magnitude
(A ≈ 0.049). However, the direction of the dipole is observer-dependent: it must always
point radially away from the singular absolute cosmic origin. The dual-front propagation
also resolves the CMB smoothness paradox. Because the 2D causal conduit propagates
instantaneously relative to the 3D lock-in lag, thermal information crosses the universe
dozens of times before crystallization completely fixes the structure. The observed ∼ 10
−5
temperature uniformity is natural; the dipole asymmetry is purely a structural byproduct
of varying lock-in times.
6 Observational Comparison: The Power Asymmetry
It is crucial to distinguish between a kinematic temperature dipole and an intrinsic
power asymmetry. The Planck satellite measures a hemispherical power asymmetry
(A = 0.066 ± 0.021)—a modulation in the actual amplitude of the primordial density
fluctuations (∆T/T ) across the sky [1].
Standard ΛCDM cannot produce this because simultaneous reheating dictates that
perturbations everywhere evolved for the exact same amount of time. However, in the
5
SSM, the fractional age gradient (∆t/t
0
= 2p ≈ 0.0996) means that one hemisphere of the
observable universe is older, meaning its localized density perturbations have had more
time to structurally evolve and grow.
In standard cosmological perturbation theory, super-horizon density perturbations
(δ = δρ/ρ) during the radiation-dominated era grow linearly with cosmic time (δ ∝ t).
Via the cosmological Poisson equation in Fourier space (Φ ∝ δ/k
2
), the gravitational po-
tential Φ inherits this identical linear time dependence (Φ ∝ t). Assuming scale-invariant
initial conditions and a hemisphere-symmetric transfer function, this structural maturity
couples directly to the observable sky via the Sachs-Wolfe effect (∆T /T = −Φ/3), which
linearly links the temperature fluctuations to the underlying gravitational potential.
The Planck collaboration defines the hemispherical asymmetry phenomenologically as
a linear dipole modulation of the temperature field: ∆T (ˆn) = ∆T
iso
(ˆn)[1+Aˆp·ˆn]. Because
the linear amplitude of the temperature fluctuation scales directly with the perturbation
growth (Φ ∝ t), the modulation amplitude A (measuring the maximum deviation from
the mean) is exactly half the total fractional age gradient across the sky:
A =
∆Φ
Φ
=
1
2
∆t
t
0
=
1
2
(2p) = p = e
−3
≈ 0.049 (9)
Because the SSM structurally predicts A ≈ 0.049, we can cleanly compare its performance
against the strict isotropic demands of ΛCDM (Table 1). While 0.049 sits slightly below
the Planck central value of 0.066, it falls comfortably within the 1σ error margin (±0.021).
Metric ΛCDM SSM
Reheating mechanism Simultaneous hypersurface Propagating front
Origin geometry N/A (global) Single nucleation event
Predicted power asymmetry A 0.000 0.049
Observed asymmetry A 0.066 ± 0.021 [1] 0.066 ± 0.021
Tension with observation > 3σ excluded Within 1σ (0.8σ)
Present-day age difference 0 Gyr 1.37 Gyr
Table 1: Comparison of hemispherical asymmetry predictions. ΛCDM predicts zero asym-
metry, ruled out by Planck. The SSM single-origin wavefront naturally accommodates
the observed dipole amplitude.
7 Falsifiable Predictions: Co-Aligned Structure Dipoles
A literal 1.37 billion-year age gap across the modern sky must leave macroscopic imprints
on the evolution of matter. Because the gradient originates from a singular wavefront, all
structural anomalies must be perfectly co-aligned along the same preferred axis.
Quasar density and evolution. Supermassive black holes in the older hemisphere had
significantly more time to grow. We predict a statistically significant dipole in the evo-
lutionary states of high-redshift AGN co-aligned with the CMB asymmetry axis. Cru-
cially, Sinha et al. (2023) [10] measured a quasar spectral index dipole pointing toward
(l, b) ≈ (201.5
◦
, −29.4
◦
), closely co-aligned with the CMB hemispherical power asymme-
try axis at (l, b) ≈ (221
◦
, −27
◦
) [1]. The SSM predicts this is not coincidence but the
direct macroscopic imprint of the ∼ 1.37 Gyr chronological gradient.
6
Galaxy cluster maturity. Massive clusters in the older hemisphere had an additional
∼ 1.37 Gyr to accrete mass and virialize, predicting statistically higher mean cluster
masses and more relaxed morphologies at equivalent redshifts.
BAO scale asymmetry. A ∼ 5% age difference at the CMB epoch subtly alters the
physical sound horizon r
d
between hemispheres, translating to a shift of roughly ∼ 7 Mpc
in the standard BAO ruler. The DESI survey [11] measures r
d
to ∼ 0.5% global precision;
a hemispheric split would yield ∼ 7σ detection of this geometric axis.
8 Conclusions
By treating the early universe as a discrete tensor network undergoing a volumetric phase
transition from a single nucleation event, we derive a scale-free age gradient ∆t/t
0
≈
2e
−3
≈ 0.0996 through an explicit chain: the tetrahedral crossing number c = 3 dic-
tates the topological activation barrier, setting the tunneling probability p = e
−3
, which
determines the fractional kinematic lag of the 3D lock-in phase. Evaluated today, one
hemisphere is ∼ 1.37 billion years older than its antipode. The derived CMB dipole am-
plitude A = p ≈ 0.049 falls within 1σ of the Planck measurement (0.066 ± 0.021), while
the ΛCDM prediction of A = 0 is excluded at > 3σ. The framework predicts co-aligned
structure formation dipoles in quasar evolution, cluster maturity, and BAO scale that are
testable with DESI, Euclid, and the Vera C. Rubin Observatory.
Data Availability
No new observational data were generated. The interactive 3D visualizations of the SSM
phase transition sequence and the macroscopic Big Bang wavefront are openly available
at https://raghu91302.github.io/ssmtheory/ssm_regge_deficit.html and https:
//raghu91302.github.io/ssmtheory/ssm_bigbang_wavefront.html.
References
[1] Planck Collaboration, “Planck 2018 results. VII. Isotropy and statistics of the CMB,”
Astron. Astrophys. 641, A7 (2020).
[2] A. H. Guth, “Inflationary universe: A possible solution to the horizon and flatness
problems,” Phys. Rev. D 23, 347 (1981).
[3] H. K. Eriksen et al., “Asymmetries in the Cosmic Microwave Background anisotropy
field,” Astrophys. J. 605, 14-20 (2004).
[4] R. Kulkarni, “Geometric Phase Transitions in a Discrete Vacuum: Deriving Cosmic
Flatness, Inflation, and Reheating from Tensor Network Topology,” Preprint (In
Review), Zenodo: 10.5281/zenodo.18727238 (2026).
[5] R. Kulkarni, “Constructive Verification of K=12 Lattice Saturation: Exploring Kine-
matic Consistency in the Selection-Stitch Model,” Preprint (In Review), Zenodo:
10.5281/zenodo.18294925 (2026).
[6] R. Kulkarni, “Geometric Emergence of Spacetime Scales,” Preprint (In Review), Zen-
odo: 10.5281/zenodo.18752809 (2026).
7
[7] R. Kulkarni, “A Topological Ansatz for the Proton-to-Electron Mass Ratio,” Preprint
(In Review), Zenodo: 10.5281/zenodo.18253326 (2026).
[8] R. Kulkarni, “Matter as Frozen Phase Boundaries: Quark Structure, Fractional
Charges, and Color Confinement from Tetrahedral Defects in a K = 12 Vacuum
Lattice,” Preprint (In Review), Zenodo: 10.5281/zenodo.18917946 (2026).
[9] N. J. Secrest et al., “A Test of the Cosmological Principle with Quasars,” Astrophys.
J. Lett. 908, L51 (2021).
[10] K. Sinha et al., “Probing Cosmological Principle using the spectral index of quasar
flux distribution,” arXiv:2312.07596 (2023).
[11] DESI Collaboration, “DESI 2024 VI: Cosmological Constraints from the Measure-
ments of Baryon Acoustic Oscillations,” arXiv:2404.03002 (2024).
8
