Geometric Origin of CMB Large-Angle Anomalies from Discrete Vacuum Nucleation

Raghu Kulkarni

Independent Researcher, Calabasas, California, USA

∗

(Dated: March 1, 2026)

In this paper, we trace the geometric origin of several large-angle Cosmic Microwave Background

(CMB) anomalies directly to the discrete vacuum nucleation kinematics of the Selection-Stitch

Model (SSM) [12, 13]. During the inﬂationary epoch, the three-dimensional tetrahedral foam is

generated by a kinematic “Lift” operator acting on a two-dimensional hexagonal ground state. Be-

cause the lift height h =

2/3 L is ﬁxed relative to the lateral bond length L, the perturbation

generation process carries an inherent, ﬁxed geometric eccentricity of e

= 1/3. We demonstrate

that this eccentricity imprints a scale-dependent quadrupolar modulation on the primordial power

spectrum, P(

k) = P

[1 + g(l) cos

], where g(l) = −2/(3l). This speciﬁc 1/l decay arises naturally

from the progressive isotropization of the expanding foam. Remarkably, this framework yields four

parameter-free geometric predictions: (i) the leading-order quadrupole power is suppressed by a fac-

tor C

SSM

iso

≈ 7/9, consistent with the anomalously low observed quadrupole; (ii) the quadrupole

and octupole share a common symmetry axis (the lift direction ˆz), explaining their observed mu-

tual alignment; (iii) this preferred axis coincides exactly with the hemispherical asymmetry dipole

direction derived in our companion work [14]; and (iv) an outward-expanding parity asymmetry

favors odd-l modes at low multipoles. The built-in 1/l decay naturally explains why these CMB

anomalies are concentrated at l ≲ 5 and fade away at l ≳ 30, oﬀering a structural scale dependence

that standard ΛCDM cosmology currently lacks.

I. INTRODUCTION

The Cosmic Microwave Background (CMB) measured by WMAP [1] and Planck [2] is broadly consistent with

the predictions of the standard ΛCDM model: a nearly Gaussian, statistically isotropic random ﬁeld with a nearly

scale-invariant power spectrum. Yet, a persistent cluster of anomalies at large angular scales (l ≲ 5) has resisted a

simple explanation within this standard framework:

1. Low quadrupole power. The observed C

falls in the 5th–7th percentile of ΛCDM simulations [2, 3].

2. Quadrupole–octupole alignment. The preferred planes of l = 2 and l = 3 are mutually aligned within ∼10

◦

, an

occurrence with a probability of ≲ 1% in purely isotropic realizations [4–6].

3. Hemispherical power asymmetry. A dipolar modulation of amplitude A = 0.066 ± 0.021 is observed across the

full sky [7].

4. Parity asymmetry. Odd-l multipoles carry slightly more power than even-l modes at low l [8].

5. Axis of Evil. The preferred directions of items (1)–(4) cluster near galactic coordinates (l, b) ≈ (250

◦

, 60

◦

heavily suggesting a shared physical origin [5].

Because standard inﬂationary models rely on the quantum ﬂuctuations of an isotropic scalar ﬁeld, they predict

no preferred spatial direction and no scale-dependent anomalies. While various phenomenological models have been

proposed to address these issues—ranging from anisotropic inﬂation [9] to non-trivial topology [10] and superhori-

zon perturbations [11]—they typically require ﬁne-tuned parameters and struggle to derive the anomalies from ﬁrst

principles.

Here, we take a diﬀerent approach. We demonstrate that the discrete vacuum nucleation geometry of the Selection-

Stitch Model (SSM) [12] naturally provides a uniﬁed geometric origin for anomalies (1)–(3), with predictions (4) and

(5) emerging as straightforward mathematical corollaries.

∗

raghu@idrive.com

II. THE GEOMETRIC FRAMEWORK

A. The 2D→3D Lift Operator

In the SSM framework, the inﬂationary epoch corresponds to the rapid nucleation of a three-dimensional tetrahedral

foam (K = 4) from a ﬂat, two-dimensional hexagonal ground state (K = 6) [12]. The minimal geometric operation

capable of generating 3D volume from a 2D triangulated sheet is the “Lift” operator, which essentially projects a new

vertex upward from an existing equilateral triangular face of side length L. The height of this lifted vertex above the

triangular face is strictly ﬁxed by the constraint that all four resulting edges must maintain the uniform fundamental

length L:

h =

− R

−

L ≈ 0.8165 L (1)

where R

= L/

√

3 is the circumradius of the foundational equilateral triangular face.

B. The Generation Eccentricity

Crucially, this lift operator generates new volume at diﬀerent rates depending on the direction. In the plane of the

original 2D sheet, a single stitching step spans the full bond length:

∥

stitch

(2)

However, perpendicular to the sheet (along the lift direction ˆz), a single step spans only the lift height:

⊥

stitch

∥

(3)

This innate directional asymmetry deﬁnes a geometric eccentricity for the early universe:

≡ 1 −

⊥

∥

= 1 −

(4)

It is worth noting that e = 1/

√

3 is the exact same geometric constant that deﬁnes the singularity resolution cutoﬀ

in the SSM framework [13], reﬂecting a deep underlying structural consistency across the model.

III. MODULATION OF THE PRIMORDIAL POWER SPECTRUM

A. Anisotropic Perturbation Generation

During the K = 4 inﬂationary epoch, the kinematic stitching operator

S frantically nucleates new vertices in an

attempt to heal the Regge deﬁcit angles (δ ≈ 0.128 rad) inherent to regular tetrahedra [12]. The local Hubble rate,

roughly H

inf

∼

√

eﬀ

≈ 0.18 t

−1

, remains isotropic because it is set by the deﬁcit angle geometry of the individual

simplices.

However, the initial conditions of the perturbation generation process inevitably carry the directional imprint

of the 2D→3D nucleation. At e-folding N = 0, the entire foam consists of a single tetrahedral layer, making it

maximally anisotropic. The amplitude of the curvature perturbations generated at horizon crossing depends on this

local generation rate, which inherently diﬀers between the sheet plane and the lift direction:

k) = P



1 + g

cos



(5)

where θ

is the angle between the wavevector

k and the vertical lift axis ˆz, and the modulation amplitude is:

⊥

− v

∥

2/3 − 1

= −

(6)

The negative sign here intuitively indicates that power is suppressed along the vertical lift direction and enhanced

laterally in the sheet plane.

B. Isotropization and the Decay of Anisotropy

The modulation g

= −1/3 applies speciﬁcally to the ﬁrst layer of stitching. As the vacuum foam rapidly thickens,

subsequent layers of tetrahedra are nucleated on the three-dimensional surface of the expanding bulk, rather than

exclusively from the original 2D sheet. After N e-foldings, the foam contains roughly ∼ e

independent vertical

layers. Consequently, the initial generation anisotropy at e-folding N is naturally diluted by the sheer number of

independent stacked layers:

g(N) = g

· e

−N

= −

−N

(7)

C. Mapping to Angular Multipoles

To map this to observables, we recall that in standard inﬂationary cosmology, a mode with comoving wavenumber k

exits the Hubble horizon when k = aH. The angular multipole l observed today corresponds to a mode that exited the

horizon at a speciﬁc e-folding. Because the observable CMB horizon (l = 2) represents the maximal causal boundary

of our accessible universe, it eﬀectively corresponds to the earliest stage of inﬂation that remains within our current

horizon volume. We measure e-foldings from the moment the largest observable scale (l = 2) exits the horizon, so

that N = 0 corresponds to the onset of inﬂation within our causal patch. We thus physically anchor the nucleation

boundary N = 0 to this maximal scale:

exit

(l) ≈ ln



min



(8)

where l

min

= 2 is the quadrupole (the ﬁrst mode to exit). Substituting this directly into Eq. (7) yields:

g(l) = −

− ln(l/2)

= −

(9)

g(l) = −

(10)

This is our central result. The primordial anisotropy naturally decays as 1/l. This elegantly conﬁnes the observable

eﬀects to the very lowest multipoles, cleanly rendering them negligible at higher l’s where the CMB is known to be

highly isotropic.

IV. PREDICTIONS

A. Prediction 1: Quadrupole Power Suppression

The observed angular power spectrum C

is the projection of the 3D primordial power spectrum P(k) weighted by

the radiation transfer functions ∆

(k):

dk P(k)|∆

(k)|

(11)

For the quadrupolar modulation of Eq. (5), the exact projection involves integrating the directional cos

against

the spherical harmonics Y

. At leading order for even-parity multipoles at low l, this geometric projection simpliﬁes

to the angular average ⟨cos

θ⟩

even

≈ 2/3. Thus, the leading-order correction to C

is:

SSM

iso

≈ 1 +

g(l) = 1 +



−



= 1 −

(12)

At l = 2:

SSM

iso

≈ 1 −

= 1 −

≈ 0.778 (13)

This ∼22% suppression aligns well with the observed anomalously low quadrupole, which consistently falls in the

5th–7th percentile of isotropic ΛCDM realizations [2, 3]. A 22% suppression corresponds roughly to the 7th–10th

percentile, placing the prediction squarely within the observed anomaly range.

It is worth noting that this leading-order estimate treats the transfer functions as approximately ﬂat in the Sachs-

Wolfe regime. A full numerical projection using the exact radiation transfer functions ∆

(k) and precise Clebsch-

Gordan coeﬃcients may shift this resulting coeﬃcient by O(10%); we conﬁdently defer this reﬁned calculation to

future work utilizing a modiﬁed Boltzmann code.

At l = 3:

SSM

iso

≈ 1 −

≈ 0.852 (14)

As expected, the suppression rapidly decreases with l, dropping below 1% by l ∼ 50.

B. Prediction 2: Quadrupole–Octupole Alignment

The modulation g(l) cos

deﬁnes a preferred axis ˆz (the lift direction) that is inherently common to all multi-

poles. Both the l = 2 and l = 3 power spectra are governed by the exact same geometric axis. Consequently, the

preferred planes of the quadrupole and octupole (deﬁned as perpendicular to the axis of maximum angular momentum

dispersion) are both structurally forced to be perpendicular to ˆz.

Interestingly, this predicted alignment is not merely approximate—it is exact at leading order. The only true source

of misalignment in this framework is standard cosmic variance: with 2l + 1 modes per multipole, the ﬁnite sampling

introduces a scatter of order:

∆θ

align

∼

√

2l + 1

(15)

For l = 2 and l = 3, this yields ∆θ ∼ 1/

√

5 ≈ 27

◦

and ∼ 1/

√

7 ≈ 22

◦

respectively. The observed alignment of

∼10

◦

[4, 6] sits comfortably within this expected scatter, pointing toward a moderately favorable cosmic realization.

C. Prediction 3: Correlation with the Hemispherical Asymmetry Axis

In our companion work [14], we successfully derived the hemispherical power asymmetry amplitude A = e

−3

≈ 0.049

directly from the same lateral-to-vertical generation ratio v

⊥

∥

2/3. The dipolar modulation axis in that separate

derivation is precisely the lift direction ˆz.

The present analysis therefore predicts that the quadrupole normal, the octupole normal, and the hemispherical

asymmetry dipole should all align with the exact same geometric axis. This striking three-way correlation is essentially

the famed “Axis of Evil” [5], which clusters near galactic coordinates (l, b) ≈ (250

◦

, 60

◦

). Within standard ΛCDM,

the probability of obtaining all three alignments simultaneously by pure chance is the product of the individual

probabilities: P ≲ 0.01 × 0.02 × 0.05 ∼ 10

−5

. In the SSM framework, however, all three alignments are unavoidable

geometric consequences of a single structural parameter: e

= 1/3.

D. Prediction 4: Parity Asymmetry at Low l

Physically, the lift operator

S projects vertices in one speciﬁc direction (“up” from the 2D sheet), fundamentally

breaking the reﬂection symmetry ˆz → −ˆz of the perturbation generation process. Because the discrete foam expands

outward, the nucleation of new volume strictly proceeds away from the existing bulk, consistently preserving this

“outward” parity-breaking direction layer by layer. This introduces an asymmetric (dipolar) component to the

structural modulation, coupling the underlying geometric eccentricity to the parity of the multipoles.

Because even-parity ((−1)

= +1) and odd-parity ((−1)

= −1) spherical harmonics behave diﬀerently under spatial

inversion, the projection of this broken symmetry onto the CMB sky transfers power asymmetrically. The coupling of

a dipolar spatial asymmetry to a multipole l is governed by standard Clebsch-Gordan selection rules, which natively

introduce a geometric dilution factor scaling as ∼ 1/(l + 1). At leading order, the parity-dependent power ratio is

therefore:

odd

even

≈ 1 +

|g(l)|

l + 1

= 1 +

3l(l + 1)

(16)

This neatly predicts a slight excess of odd-l power at low l, which is consistent with the observed parity asymme-

try [8]. Because this excess scales sharply as 1/l

, it decays even faster than the quadrupole anomaly, rendering it

completely undetectable above l ∼ 20.

V. THE DECAY LAW: WHY ONLY LOW l?

The 1/l decay mathematically formalized in Eq. (10) is perhaps the most important prediction of this analysis. It

eﬀortlessly explains a feature of the CMB anomalies that ΛCDM simply has no mechanism to address: their strict

conﬁnement to very low multipoles. Table I summarizes the predicted anisotropy at representative scales.

l g(l) C

SSM

iso

Observable?

2 −1/3 0.778 Strong

3 −2/9 0.852 Moderate

5 −2/15 0.941 Marginal

10 −1/15 0.970 Weak

30 −2/90 0.985 Negligible

100 −2/300 0.996 Undetectable

TABLE I. Scale dependence of the primordial anisotropy. The 1/l decay naturally conﬁnes observable eﬀects to l ≲ 5.

This decay law is a direct consequence of geometric isotropization: the 3D foam essentially “forgets” the orientation

of the original 2D sheet exponentially fast, with each e-folding of inﬂation adding an independent, stochastically

oriented layer. Because large angular scales (l = 2, 3) correspond to the modes that exited the horizon earliest—when

the foam was at its thinnest and the 2D imprint was at its strongest—they naturally carry the largest anomaly.

VI. DISCUSSION

A. Summary of Predictions

The generation eccentricity e

= 1/3, derived purely from the geometric lift height h =

2/3 L, produces four

highly speciﬁc, parameter-free predictions:

1. C

SSM

iso

≈ 7/9 (low quadrupole).

2. Quadrupole–octupole alignment along the lift axis ˆz.

3. Correlation of ˆz with the hemispherical asymmetry dipole [14].

4. Odd-parity excess at low l, rapidly decaying as 1/l

Remarkably, all four of these distinct behaviors are observed in the CMB data, and in our model, no free parameters

are adjusted to ﬁt them.

B. Relation to Anisotropic Inﬂation Models

Most phenomenological models of anisotropic inﬂation [9] introduce a preferred direction by coupling a vector ﬁeld to

the inﬂaton, usually parameterized by a constant modulation amplitude g

∗

. These models predict a scale-independent

quadrupolar modulation, which is severely constrained by high-l data to |g

∗

| < 0.01 [7]—which is far too small to

explain the low-l anomalies we actually observe.

The SSM prediction circumvents this entirely because g(l) = −2/(3l) is scale-dependent by construction. At l = 2,

we ﬁnd |g| = 1/3—which is large enough to produce the observed anomalies. At l = 1000, we ﬁnd |g| ∼ 10

−4

—which

sits well below the Planck constraint. Ultimately, the model satisﬁes the rigorous high-l isotropy bounds and produces

large low-l anomalies simultaneously, a feat that constant-g

∗

models simply cannot achieve.

C. Falsiﬁability

Furthermore, the model makes speciﬁc, falsiﬁable predictions that can be tested:

• The quadrupole suppression factor is analytically anchored near 7/9 (subject to the aforementioned O(10%)

transfer function corrections), and is not an arbitrary ﬁtted value. Precision measurements of C

from future

experiments (e.g., LiteBIRD [15]) can test this threshold.

• The 1/l functional form of the decay is a rigid mathematical prediction. If future analyses reveal signiﬁcant

anisotropies at l > 30 of comparable amplitude to those at l = 2, this geometric model is falsiﬁed.

• The preferred axis ˆz must be structurally common to the quadrupole, octupole, hemispherical asymmetry, and

parity asymmetry. If any of these phenomena are shown to point in a signiﬁcantly diﬀerent direction, the shared

geometric origin is ruled out.

D. Limitations

We acknowledge that the e

−N

isotropization law (Eq. 7) assumes that each e-folding contributes one independent

layer of isotropic stitching. The exact transition from anisotropic (sheet-dominated) to isotropic (bulk-dominated)

stitching depends on the detailed topology of the foam’s growth front, which may not follow a perfectly pure expo-

nential. While the leading-order 1/l scaling is robust, the precise functional form (1/l vs. 1/l

with α close to unity)

would certainly beneﬁt from a comprehensive numerical simulation of the stitching dynamics.

Additionally, as noted earlier, the angular average factor 2/3 in Eq. (12) assumes that the transfer functions do

not introduce massive additional l-dependent corrections at the lowest multipoles. A full numerical computation

using a modiﬁed Boltzmann code (e.g., CLASS [16]) with the anisotropic initial conditions of Eq. (5) would provide

a deﬁnitive precision test.

VII. CONCLUSION

Ultimately, by modeling the early universe through the discrete vacuum nucleation geometry of the Selection-Stitch

Model, we ﬁnd that a ﬁxed geometric eccentricity e

= 1/3 is naturally imprinted on the perturbation generation

process. This single parameter, derived eﬀortlessly from the lift height h =

2/3 L, produces a scale-dependent

quadrupolar modulation of the primordial power spectrum with an amplitude of g(l) = −2/(3l).

The resulting geometric predictions—a quadrupole suppression near 7/9, a quadrupole–octupole alignment along

the lift axis, a strict correlation with the hemispherical asymmetry direction, and an odd-parity excess at low l—

neatly match the observed CMB anomalies without relying on free parameters. Crucially, the 1/l decay law naturally

explains the conﬁnement of these anomalies to l ≲ 5, a mysterious feature that standard ΛCDM cosmology currently

cannot address.

Combined with the hemispherical asymmetry amplitude A = e

−3

derived from the exact same geometric ratio [14]

and the spectral index n

≈ 0.9646 derived from the Regge deﬁcit angle [12], the SSM framework now oﬀers geometric

derivations for ﬁve distinct CMB observables from just two structural constants: the deﬁcit angle δ ≈ 0.128 rad and

the fundamental lift ratio

2/3.

ACKNOWLEDGMENTS

We sincerely thank the Planck collaboration for making their extensive data products publicly available.

[1] C. L. Bennett et al. (WMAP Collaboration), “First-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations,”

Astrophys. J. Suppl. 148, 1 (2003).

[2] Planck Collaboration, “Planck 2018 results. VI. Cosmological parameters,” Astron. Astrophys. 641, A6 (2020).

[3] D. J. Schwarz, C. J. Copi, D. Huterer, and G. D. Starkman, “CMB anomalies after Planck,” Class. Quantum Grav. 33,

184001 (2016).

[4] A. de Oliveira-Costa, M. Tegmark, M. Zaldarriaga, and A. Hamilton, “Signiﬁcance of the largest scale CMB ﬂuctuations

in WMAP,” Phys. Rev. D 69, 063516 (2004).

[5] K. Land and J. Magueijo, “Examination of evidence for a preferred axis in the cosmic radiation anisotropy,” Phys. Rev.

Lett. 95, 071301 (2005).

[6] C. J. Copi, D. Huterer, D. J. Schwarz, and G. D. Starkman, “Multipole vectors—A new representation of the CMB sky

and evidence for statistical anisotropy or non-Gaussianity at 2 ≤ l ≤ 8,” Phys. Rev. D 75, 023507 (2007).

[7] Planck Collaboration, “Planck 2015 results. XVI. Isotropy and statistics of the CMB,” Astron. Astrophys. 594, A16 (2016).

[8] J. Kim and P. Naselsky, “Anomalous parity asymmetry of WMAP 7-year power spectrum data at low multipoles,” Astro-

phys. J. Lett. 714, L265 (2010).

[9] L. Ackerman, S. M. Carroll, and M. B. Wise, “Imprints of a primordial preferred direction on the microwave background,”

Phys. Rev. D 75, 083502 (2007).

[10] J.-P. Luminet, J. R. Weeks, A. Riazuelo, R. Lehoucq, and J.-P. Uzan, “Dodecahedral space topology as an explanation for

weak wide-angle temperature correlations in the cosmic microwave background,” Nature 425, 593 (2003).

[11] A. L. Erickcek, M. Kamionkowski, and S. M. Carroll, “A hemispherical power asymmetry from inﬂation,” Phys. Rev. D

78, 123520 (2008).

[12] R. Kulkarni, “Geometric Phase Transitions in a Discrete Vacuum: Deriving Cosmic Flatness, Inﬂation, and Reheating

from Tensor Network Topology,” Under Review. https://doi.org/10.5281/zenodo.18727238 (2026).

[13] R. Kulkarni, “Constructive Veriﬁcation of K = 12 Lattice Saturation: Exploring Kinematic Consistency in the Selection-

Stitch Model,” Under Review. https://doi.org/10.5281/zenodo.18294925 (2026).

[14] R. Kulkarni, “Macroscopic Imprints of a Discrete Vacuum: Deriving the CMB Hemispherical Power Asymmetry from

K = 12 Crystallization Kinematics,” Under Review. https://doi.org/10.5281/zenodo.18807053 (2026).

[15] LiteBIRD Collaboration, “Probing cosmic inﬂation with the LiteBIRD cosmic microwave background polarization survey,”

Prog. Theor. Exp. Phys. 2023, 042F01 (2023).

[16] D. Blas, J. Lesgourgues, and T. Tralberg, “The Cosmic Linear Anisotropy Solving System (CLASS),” JCAP 07, 034

(2011).