Kreuzer–Skarke Axiverse Overview

• Kreuzer–Skarke Axiverse is the ensemble of axion effective theories derived from toric Calabi–Yau threefolds via reflexive polytope triangulations.
• The framework utilizes detailed toric data—such as Hodge numbers, intersection tensors, and cone constraints—to infer axion multiplicities and low-energy couplings.
• Universal statistical features from large ensembles constrain axion masses, decay constants, and photon couplings, impacting observable string phenomenology.

Searching arXiv for relevant papers on the Kreuzer–Skarke axiverse and foundational database. The Kreuzer–Skarke axiverse is the ensemble of axion effective theories obtained from string compactifications on Calabi–Yau threefold hypersurfaces associated with the Kreuzer–Skarke classification of four-dimensional reflexive polytopes. In the underlying toric construction, a fine, regular, star triangulation of a reflexive polytope defines a toric ambient space, and a generic anticanonical hypersurface in that space is a Calabi–Yau threefold. The resulting geometric data—Hodge numbers, intersection tensors, Chern classes, and Mori and Kähler cone information—provide the inputs from which axion multiplicities, kinetic terms, instanton actions, and low-energy couplings are inferred (Altman et al., 2014). Subsequent work turned this toric dataset into an axiverse in the technical sense: at large $h^{1,1}$, the geometry generically produces many Ramond–Ramond four-form axions with highly suppressed nonperturbative masses, together with distinctive statistical regularities across enormous ensembles of compactifications (Demirtas et al., 2018).

1. Database origin and toric construction

The foundational object is the Kreuzer–Skarke list of all $473{,}800{,}776$ reflexive polytopes in four dimensions, together with $30{,}108$ distinct Hodge pairs $(h^{1,1},h^{2,1})$. In the toric hypersurface construction, a four-dimensional reflexive polytope $\Delta$ and its dual $\Delta^\ast$ encode singular ambient toric varieties in which Calabi–Yau threefold hypersurfaces can appear. A triangulation of $\Delta^\ast$ yields a simplicial toric fourfold $A$, and the Calabi–Yau threefold $X\subset A$ is realized as an anticanonical divisor (Altman et al., 2014).

The companion database to the toric Calabi–Yau construction reports, for each threefold obtained from such triangulations, the Hodge numbers $h^{1,1}(X)$, $473{,}800{,}776$0, the Euler number $473{,}800{,}776$1, the triple intersection numbers $473{,}800{,}776$2, the second Chern class $473{,}800{,}776$3, Mori cone generators, Kähler cone inequalities, toric divisor classes, the Stanley–Reisner ideal, ambient Chern classes, and triangulation-specific cone data. It also distinguishes between different triangulations that yield genuinely distinct Calabi–Yau geometries and triangulations that represent different Kähler phases of the same topology, gluing the latter when topological equivalence is established via matching of invariants such as $473{,}800{,}776$4 and $473{,}800{,}776$5 (Altman et al., 2014).

The explicitly computed subset is complete for $473{,}800{,}776$6 and partial for $473{,}800{,}776$7, comprising $473{,}800{,}776$8 polytopes, $473{,}800{,}776$9 threefolds, and about $30{,}108$0 favorable geometries. The dependence on triangulation is not a minor technicality. For the $30{,}108$1 polytope $30{,}108$2, the database finds two distinct glued geometries with identical Hodge numbers $30{,}108$3, $30{,}108$4, but different $30{,}108$5, $30{,}108$6, and Mori/Kähler cones; one geometry is assembled from two triangulations, the other from one. This directly shows that Hodge data alone do not determine axion physics (Altman et al., 2014).

At the classification level, the Kreuzer–Skarke program relies on canonical normal forms for lattice polytopes. The PALP normal form was used by Kreuzer and Skarke to classify $30{,}108$7D and $30{,}108$8D reflexive polytopes, preventing double-counting and organizing the database by isomorphism classes. This infrastructure does not itself address axion physics, but it is essential for the uniqueness and scalability of the datasets on which KS axiverse studies depend (Grinis et al., 2013).

2. Topological data, mirror symmetry, and cone geometry

For toric hypersurfaces, the Hodge numbers are determined combinatorially by Batyrev’s formulas. If $30{,}108$9 denotes the number of lattice points and $(h^{1,1},h^{2,1})$0 the number of interior lattice points of a face $(h^{1,1},h^{2,1})$1, then for reflexive $(h^{1,1},h^{2,1})$2 and $(h^{1,1},h^{2,1})$3,

$(h^{1,1},h^{2,1})$4

$(h^{1,1},h^{2,1})$5

Mirror symmetry exchanges $(h^{1,1},h^{2,1})$6, so $(h^{1,1},h^{2,1})$7 and $(h^{1,1},h^{2,1})$8 (Altman et al., 2014).

For favorable polytopes, inherited toric divisors generate $(h^{1,1},h^{2,1})$9 and the effective cone studied in large-scale axiverse analyses. This favorable sector is especially important because it provides explicit divisor representatives and intersection data suitable for systematic scans (Demirtas et al., 2018).

A major structural result of the KS axiverse is that the Kähler cone becomes very narrow at large $\Delta$0. In the large-$\Delta$1 study of $\Delta$2 Calabi–Yau hypersurfaces, including at least one triangulation for every polytope with $\Delta$3, the ambient-space Mori cone $\Delta$4 and its dual Kähler cone $\Delta$5 were computed for all models, while triple intersection numbers were computed for $\Delta$6. An outer approximation $\Delta$7 to the Calabi–Yau Kähler cone $\Delta$8 is defined by requiring positivity of the volumes of $\Delta$9, all inherited divisors, and all toric intersection curves, yielding

$\Delta^\ast$0

The corresponding stretched cones $\Delta^\ast$1 impose the condition that relevant cycle volumes are at least $\Delta^\ast$2 in string units (Demirtas et al., 2018).

Quantitatively, the narrowness is expressed through the growth of the minimal Euclidean distance from the origin to the stretched cones: $\Delta^\ast$3 Under the condition that every effective curve has volume $\Delta^\ast$4, the typical divisor and threefold volumes grow as powers of $\Delta^\ast$5: $\Delta^\ast$6

$\Delta^\ast$7

Because the number of nonvanishing $\Delta^\ast$8 grows only linearly with $\Delta^\ast$9, and the RMS size of nonzero intersections is roughly independent of $\Delta^\ast$0, this scaling is attributed primarily to cone geometry rather than to explosive growth in intersection combinatorics (Demirtas et al., 2018).

3. From toric geometry to axion effective theories

The Kähler form is expanded in a divisor basis as $\Delta^\ast$1, and the toric data determine the standard geometric quantities

$\Delta^\ast$2

The Calabi–Yau condition gives $\Delta^\ast$3, while the second Chern class is obtained from the ambient toric data by adjunction, with $\Delta^\ast$4 in the formulation used in the database (Altman et al., 2014).

In type IIB orientifolds, the principal axions of the KS axiverse arise from the Ramond–Ramond four-form $\Delta^\ast$5 reduced on four-cycles. Writing the complexified Kähler moduli as $\Delta^\ast$6 or $\Delta^\ast$7, the axions are the imaginary parts of these moduli. The tree-level Kähler potential is

$\Delta^\ast$8

and the axion kinetic matrix is derived from the corresponding Kähler metric. The nonperturbative superpotential and potential are taken schematically as

$\Delta^\ast$9

$A$0

with instanton actions scaling as $A$1. In the hierarchical approximation one uses $A$2 (Cheng et al., 16 Jul 2025).

The same toric data control the allowed large-volume regions. The Mori cone generators constrain the Kähler parameters through linear inequalities $A$3, so $A$4 and the Kähler cone together determine the reachable $A$5, and hence the instanton actions and decay constants. The database quantity $A$6 can be combined with $A$7 to compute

$A$8

which is used as a diagnostic for divisor rigidity in instanton analyses (Altman et al., 2014).

Axion counting depends on compactification details. In a simplifying type IIB setup often used in KS ensemble studies, the orientifold is assumed not to project out Kähler moduli, so the number of axions is taken to be $A$9. More generally, the orientifold splits $X\subset A$0 into $X\subset A$1 and $X\subset A$2, affecting multiplicities, kinetic terms, and couplings, and the KS hypersurface database by itself does not supply the involution needed to determine that splitting. In type IIA orientifolds, axions predominantly arise from $X\subset A$3 reduced on three-cycles, with count $X\subset A$4 up to orientifold parities (Cheng et al., 16 Jul 2025).

4. Universal statistics in the KS axiverse

A later development was the identification of universal divisor-volume statistics across the KS landscape. In the large-volume regime, restricting to favorable toric hypersurfaces and prime toric divisors inherited from the ambient toric variety, one defines normalized divisor volumes by first rescaling along a Kähler-cone ray so that $X\subset A$5, and then setting

$X\subset A$6

Sampling favorable polytopes with $X\subset A$7, with $X\subset A$8 polytopes per $X\subset A$9 for $h^{1,1}(X)$0 and all favorable polytopes above that threshold, the resulting normalized distributions were found to be universal across geometries, triangulations, and random positions in the extended Kähler cone. The evidence was quantified using $h^{1,1}(X)$1-Wasserstein distances and adjusted $h^{1,1}(X)$2, with pairwise percentage errors in $h^{1,1}(X)$3 typically below $h^{1,1}(X)$4 and decreasing with $h^{1,1}(X)$5 (Cheng et al., 16 Jul 2025).

The total normalized distribution develops a second peak at moderate $h^{1,1}(X)$6, explained by superposition of sub-distributions associated with three toric divisor types: vertex, edge, and face. These sub-distributions are themselves universal, while the prevalence of face divisors grows with $h^{1,1}(X)$7 (Cheng et al., 16 Jul 2025).

To model this statistically, the paper proposes a geometry-independent single-function fit using the Gaussian Orthogonal Ensemble level-spacing distribution,

$h^{1,1}(X)$8

together with an $h^{1,1}(X)$9-dependent mean

$473{,}800{,}776$00

Within the approximations adopted—large volume, single-instanton dominance, near-diagonal dominance in the Kähler metric, a simplified divisor basis, and volume domination by the largest $473{,}800{,}776$01-cycle parameter—this model reproduces distributions of decay constants, masses, and the stringy PQ-quality estimator with good accuracy; for $473{,}800{,}776$02, the fitted and exact ensemble distributions match with $473{,}800{,}776$03 $473{,}800{,}776$04 percentage error (Cheng et al., 16 Jul 2025).

A plausible implication is that much of the KS axiverse is governed statistically by divisor-volume hierarchies rather than by detailed, geometry-by-geometry inversion of the full Kähler metric. The universality claim, however, is explicitly tied to the large-volume regime and to generic saxion vevs in the extended Kähler cone, not to finely tuned stabilization scenarios such as KKLT- or LVS-like regions where many $473{,}800{,}776$05 are forced to be comparable (Cheng et al., 16 Jul 2025).

5. Phenomenology and observational windows

The KS axiverse has been used to study axion couplings and masses in explicit large ensembles. In type IIB toy models constructed from the KS database, axion–photon couplings were found to be systematically suppressed relative to the inverse axion periodicity for two distinct reasons. First, couplings to the QED $473{,}800{,}776$06-angle are suppressed for axion mass eigenstates that are light compared to the scale set by stringy instantons on the cycle supporting QED. Second, at large $473{,}800{,}776$07 the intersection matrix is sparse, so kinetic mixing is weak. In an ensemble of approximately $473{,}800{,}776$08 toy models constructed from $473{,}800{,}776$09 polytopes at $473{,}800{,}776$10, this yields many effectively invisible axions and only a handful with potentially observable photon couplings. The resulting phenomenology includes freeze-in production and decay of thermal axions, cosmic microwave background birefringence, X-ray spectrum oscillations, and supernova constraints on the QCD axion (Gendler et al., 2023).

The same paper finds that the number of axions per model with $473{,}800{,}776$11, where $473{,}800{,}776$12, remains small even at large $473{,}800{,}776$13. The maximum $473{,}800{,}776$14 increases with $473{,}800{,}776$15 because decay constants decrease, but the fraction of visible axions shrinks because kinetic isolation strengthens with $473{,}800{,}776$16. Under the assumptions used there, the QCD axion mass often reaches $473{,}800{,}776$17 at $473{,}800{,}776$18, near current SN1987A bounds (Gendler et al., 2023).

A separate KS scan imposed grand-unification and proton-decay constraints on type IIB O3/O7 orientifolds. Sampling up to $473{,}800{,}776$19 favorable polytopes and $473{,}800{,}776$20 triangulations per polytope for each $473{,}800{,}776$21, producing approximately $473{,}800{,}776$22 compactifications, the study fixed a random prime toric divisor to host QCD, rescaled to $473{,}800{,}776$23, discarded geometries with any divisor volume $473{,}800{,}776$24, and evaluated EFT quantities at the tip of the stretched Kähler cone. Imposing $473{,}800{,}776$25, the fraction of compactifications satisfying this condition drops to zero for $473{,}800{,}776$26, leading to the conclusion that the KS axiverse compatible with these requirements contains at most $473{,}800{,}776$27 closed-string axions. In this restricted ensemble the QCD axion mass lies within

$473{,}800{,}776$28

and the majority of compactifications satisfying $473{,}800{,}776$29 also satisfy the stated PQ-quality criterion, except at very small $473{,}800{,}776$30 (Benabou et al., 21 May 2025).

That study also reports that no KS candidates were found with $473{,}800{,}776$31 and $473{,}800{,}776$32, so it finds no axions suitable for the DESI-evolving-dark-energy scenario within the scanned ensemble (Benabou et al., 21 May 2025). This suggests that apparently different phenomenological conclusions in the KS literature often reflect different ensemble cuts: generic large-$473{,}800{,}776$33 toy-model analyses, GUT-filtered scans, and fixed visible-sector assumptions probe different slices of the same geometric landscape.

6. Statistical inference on moduli space and current limitations

The KS axiverse has also been recast as a problem in Bayesian inference on Calabi–Yau moduli space. In that framework, one samples the invariant Weil–Petersson measure on Kähler moduli space using MCMC and normalizing flows, with explicit implementations up to $473{,}800{,}776$34. The real-parameter-space metric is

$473{,}800{,}776$35

and the corresponding measure is

$473{,}800{,}776$36

which behaves as $473{,}800{,}776$37, i.e. log-flat in the Calabi–Yau volume. The sampled region is the stretched Kähler cone, with all effective curve volumes $473{,}800{,}776$38 and all prime toric divisor volumes $473{,}800{,}776$39, together with explicit volume cutoffs $473{,}800{,}776$40 and $473{,}800{,}776$41 for $473{,}800{,}776$42 and $473{,}800{,}776$43 (Jain et al., 28 Nov 2025).

Using this prior, axion mass and decay-constant distributions were constructed for all inequivalent KS Calabi–Yau threefolds with $473{,}800{,}776$44. Under the restricted WP measure, Calabi–Yau volume histograms are log-flat between the cutoff scales, while divisor volume spectra depend strongly on divisor topology: K3 fibre divisors and point blowups tend to be smaller and less correlated with $473{,}800{,}776$45, and a trimodality organized by vertex, edge, and face divisors persists under WP sampling (Jain et al., 28 Nov 2025).

The same framework was used to ask inverse questions. For mock haloscope detections at $473{,}800{,}776$46, a relatively heavy QCD axion at $473{,}800{,}776$47 with $473{,}800{,}776$48 statistically favors rigid divisors, blowup divisors, elliptic fibrations, and divisors shrinking on higher-dimensional boundaries of the cone. By contrast, DMRadio-like neV detections are described as typical at small $473{,}800{,}776$49 and therefore less constraining on topology (Jain et al., 28 Nov 2025).

A full forward model incorporating Planck 2018 temperature and polarization data together with eBOSS Ly$473{,}800{,}776$50 compressed likelihoods was then constructed for a specific $473{,}800{,}776$51 geometry. In that example, the CY-informed posterior localizes the moduli-space region supporting an ultralight axion with $473{,}800{,}776$52, $473{,}800{,}776$53, and $473{,}800{,}776$54 peaked near $473{,}800{,}776$55, narrowing the EFT-like posterior from below in mass (Jain et al., 28 Nov 2025).

Across the literature, several limitations recur. The KS hypersurface database does not specify orientifold involutions, $473{,}800{,}776$56 splittings, O-plane data, or D-brane configurations; most ensemble studies work in the large-volume regime, often in the orientifold-even sector, and do not perform explicit moduli stabilization. Additional effects—ED3 zero-mode constraints, Freed–Witten consistency conditions, loop and higher-derivative corrections, backreaction, alignment, and swampland constraints—are repeatedly identified as relevant but not fully incorporated (Altman et al., 2014). The present state of the subject is therefore best understood as a highly explicit geometric framework for statistical axion phenomenology, rather than a complete moduli-stabilized census of all string-realized axion sectors.