The Kreuzer-Skarke Axiverse

We study the topological properties of Calabi-Yau threefold hypersurfaces at large $h^{1,1}$. We obtain two million threefolds $X$ by triangulating polytopes from the Kreuzer-Skarke list, including all polytopes with $240 \le h^{1,1}\le 491$. We show that the K\"ahler cone of $X$ is very narrow at large $h^{1,1}$, and as a consequence, control of the $\alpha^{\prime}$ expansion in string compactifications on $X$ is correlated with the presence of ultralight axions. If every effective curve has volume $\ge 1$ in string units, then the typical volumes of irreducible effective curves and divisors, and of $X$ itself, scale as $(h^{1,1})^p$, with $3\lesssim p \lesssim 7$ depending on the type of cycle in question. Instantons from branes wrapping these cycles are thus highly suppressed.


Introduction
As a step toward understanding general properties of quantum gravity in four spacetime dimensions, one can study compactifications of weakly-coupled string theories on sixmanifolds whose curvatures are small in string units. Understanding what is possible in such compactifications can shed light on what is possible in quantum gravity.
Bounds on topological or geometric properties of a class of compactifications can imply interesting statements about the corresponding effective theories, such as bounds on the number of fields, the rank of the gauge group, or the diameter of moduli space. Moreover, an understanding of generic properties of compactifications can inform lowenergy model building, e.g. the fact that typical Calabi-Yau threefolds have scores or hundreds of moduli and axions suggests considering inflationary and dark sectors that might be discarded as excessively complicated from a bottom-up perspective.
The Kreuzer-Skarke database of four-dimensional reflexive polytopes [1] is a fount of data on Calabi-Yau compactifications. A fine, regular, star triangulation (FRST) of any of the 473,800,776 polytopes ∆ • in the list determines a toric variety V , in which a generic anticanonical hypersurface is a smooth Calabi-Yau threefold X. However, only the most elementary data, such as the Hodge numbers of X, can be obtained directly from the database without computation. To more fully characterize a compactification on X, one needs to compute and manipulate an FRST of ∆ • .
A key measure of the difficulty of this computation is the number N of relevant lattice points in ∆ • , or equivalently the Picard number h 1,1 of X: when X is favorable, N = h 1,1 + 4. The number of possible triangulations of a given ∆ • , and the complexity of each triangulation, both grow rapidly with N . Most publicly-available software, such as Sage, is effective only for N 10. Through a major computational effort, Altman et al. obtained the data of all compactifications with N ≤ 10 [2]. However, there have been few studies, none systematic, of compactifications on Calabi-Yau threefold hypersurfaces with h 1, 1 1. This is a critical gap in our understanding: the Kreuzer-Skarke list contains threefolds with h 1,1 as large as 491, and vast numbers of triangulations -corresponding to potentially-distinct threefolds -are possible at h 1, 1 1. To the best of our knowledge, most Calabi-Yau threefold hypersurfaces have not yet been examined.
In this work we initiate a study of Calabi-Yau threefold hypersurfaces with large Picard number. We obtain FRSTs of 2,031,335 reflexive polytopes with 2 ≤ h 1,1 ≤ 491, including one triangulation for each polytope with 240 ≤ h 1,1 ≤ 491. 1 We compute the Mori cones of the associated toric varieties V , and for 2 ≤ h 1,1 ≤ 100 we compute the intersection numbers of Calabi-Yau hypersurfaces X ⊂ V .
We first use these data to bound the Kähler cone K X of X. We know of no efficient algorithm to compute K X directly in a hypersurface with h 1,1 10, so we instead place upper and lower bounds by computing cones containing K X , and contained in K X . The Kähler cone K V of V obeys K V ⊂ K X , while a cone K ∩ associated to the intersections of divisors D ⊂ V (see §2 for a precise definition) obeys K X ⊂ K ∩ . Equipped with bounds on the Kähler cone of X, we examine the α expansion in a compactification of string theory 2 on X. For the α expansion to be under control, we will require that every holomorphic curve Σ on X has volume obeying with c a dimensionless constant, so that worldsheet instantons wrapping Σ give corrections to the effective action e −2πc . Although we will suppose that c is of order unity, one can immediately extend our findings to any desired numerical value of c. We henceforth set 2 s = 1 and c = 1, and so the constraint (1.1) reads We argue in §2 that (1.2) is a useful proxy for control of perturbative and nonperturbative corrections in the α expansion. The requirement that (1.2) holds for all Σ typically implies that some irreducible holomorphic curves have volumes 1. Moreover, some irreducible effective divisors have even larger volumes, and the total threefold volume V is larger still. At first glance these trends appear unsurprising: the number of nonvanishing triple intersection numbers κ ijk must grow with h 1,1 , and so too should V = 1 6 κ ijk t i t j t k , where t i , i = 1, . . . , h 1,1 , are the volumes of a basis of H 2 (X, Z). However, obtaining the intersection numbers κ ijk for a hypersurface with h 1, 1 1 is computationally expensive, so prior studies of this point have been very limited. In this work we precisely quantify the growth of curve, divisor, and threefold volumes with h 1,1 : our computation of K ∩ leads to lower bounds on these quantities. The volumes grow far more quickly with h 1,1 than can be accounted for by the growth of the intersection numbers alone. We find that there are only O(h 1,1 ) nonvanishing intersection numbers in each geometry, with mean size independent of h 1,1 , yet we find that V grows as (h 1,1 ) p with p ∼ 7: see §6.
The primary cause of this rapid growth of volumes is the narrowness of the Kähler cone. The Kähler cone conditions enforce inequalities relating the various curve volumes, and with increasing h 1,1 , this effect becomes more pronounced, because the number of inequalities grows. Intuitively, the Kähler cone becomes very narrow for h 1, 1 1, so to be well-separated from every wall one must be very far from the origin of the cone.
One physical consequence of this finding is that requiring control of the α ex-pansion, in the sense of (1.2), typically leads to ultralight axions, unless h 1,1 is small. As an example, in a compactification of type IIB string theory on an orientifold of a hypersurface X, the Ramond-Ramond four-form C 4 gives rise to axion fields that are massless to all orders in perturbation theory, and acquire mass from Euclidean D3-branes. Suitable holomorphic four-cycles (i.e., suitable effective divisors) support superpotential contributions [4], which are well-understood, while non-holomorphic four-cycles can support contributions to the Kähler potential. We find that for typical geometries in our ensemble, every basis constructed from generators of the cone of effective divisors (cf. §4) contains elements with volume (h 1,1 ) 3 : see Figure 6. Thus, superpotential couplings 3 give extremely small masses to some of the axions. In every geometry in our ensemble with h 1,1 > 22, the lightest axion is essentially massless, with m < 10 −33 eV. An important caveat is that our finding of rapid growth of volumes with h 1,1 is a consequence of the requirement (1.2). It is possible that α corrections to the fourdimensional action are small in some cases even if some effective curves have volumes violating (1.2). Constraining this possibility would be worthwhile, but would likely require advances in computing perturbative corrections in the α expansion. Moreover, our qualitative results would be unaffected unless most curves can be made small in string units.
We also study the radius of the axion fundamental domain for each geometry in our ensemble. Understanding whether super-Planckian displacements of an inflaton field can occur in well-controlled compactifications is a pressing problem, and one way forward is to search for geometries in which the axion field space has radius R M pl . Prior work in [5] has shown that R O(1) in every Calabi-Yau hypersurface with h 1,1 ≤ 4. Here we extend the analysis of [5] to 2×10 6 hypersurfaces with 5 ≤ h 1,1 ≤ 100. We show that R M pl for most of the geometries in our ensemble. However, in a small fraction of cases we cannot exclude the possibility of radii R M pl in the parameter regime where (1.2) holds and the α expansion is well-controlled. Obtaining definitive results in these intriguing cases would require advances in computing the Kähler cones of Calabi-Yau hypersurfaces per se, rather than just the Kähler cones of the corresponding ambient toric varieties.
The organization of this note is as follows. In §2 we review basic facts about the Kähler cone of a Calabi-Yau hypersurface in a toric variety. In §3 we introduce the notion of a stretched Kähler cone, and in §4 we explain how upper bounds on axion masses can be obtained by computing cycle volumes in an appropriate stretched Kähler cone. In §5 we describe our algorithm for computing the Kähler cones, and approximations to the Kähler cones, in an ensemble of Calabi-Yau threefold hypersurfaces constructed from the Kreuzer-Skarke database. We present our results in §6.
In §7 we explore the implications of our findings for the axion mass spectrum in type IIB compactifications. We conclude in §8. Although our findings directly involve the volumes of holomorphic cycles, in Appendix §A we discuss how instantons wrapping non-holomorphic volume-minimizing chains could be governed to good approximation by the growth of volume that we establish in the holomorphic case.

The Effective, Kähler, and Mori Cones
In this section we recall the definitions and basic properties of the effective cone, the Kähler cone, and the Mori cone of a projective algebraic variety X, and we explain how to compute approximations to these cones when X is a Calabi-Yau threefold hypersurface in a toric variety. From the data of these convex cones one can read off properties of the effective theory arising in a string compactification on X.

The effective cone
Let X be a projective algebraic variety of complex dimension n. A Weil divisor D on X is a finite formal sum of irreducible codimension-one subvarieties D A , The divisor D is called effective if the n A are all nonnegative. We define the effective cone Eff(X) to be the convex cone in H 2n−2 (X, R) spanned by the classes of effective divisors. The relevance of the effective cone is that a Euclidean D3-brane wrapping a divisor D in an orientifold of a Calabi-Yau threefold X can contribute to the superpotential only if D is effective. Intuitively, effective divisors consist of finite collections of irreducible holomorphic hypersurfaces, each of which can support BPS D-branes.

Effective divisors of a Calabi-Yau hypersurface
Let ∆ • be a four-dimensional reflexive polytope. An FRST of ∆ • defines a fan that corresponds to a simplicial toric fourfold V . The generic anticanonical hypersurface X ⊂ V is a smooth Calabi-Yau threefold [6].
Each lattice point v I on the boundary of ∆ • corresponds to a homogeneous toric coordinate x I , whose vanishing defines a prime toric divisor D I . The prime toric divisors are irreducible effective divisors on V . A subset {v A } ⊂ {v I } of the points on the boundary of ∆ • are not interior to 3-faces (facets) of ∆ • , but instead lie in faces of dimension ≤ 2. Each such lattice point v A not interior to a facet corresponds to a prime toric divisor that intersects X transversely. The restriction to X then defines a divisor D A ⊂ X, that is effective on X. Points interior to facets, on the other hand, define divisors of V that do not intersect a generic Calabi-Yau hypersurface X. In triangulating ∆ • we may therefore ignore lattice points interior to facets; such a triangulation corresponds to a maximal projective crepant partial (MPCP) desingularization, in the sense of [7]. We will restrict ourselves to such partial desingularizations.
In general, D A may be a reducible divisor on X, even though D A is irreducible on V . This occurs if and only if v A corresponds to a point in the strict interior of a 2-face f ⊂ ∆ • , and * (f • ) > 0, where * (f • ) is the number of lattice points in the strict interior of the dual face f • ⊂ ∆. The condition that all of the prime toric divisors D A on V that intersect X in fact restrict to irreducible divisors on X is thus where the sum is over all 2-faces f ⊂ ∆ • . A polytope obeying (2.3) is called favorable, and by extension we refer to the associated V and X as being favorable.
For simplicity we will confine our attention to the case where X is favorable, though we expect the results of our analysis to extend into the non-favorable regime. For X favorable, there are exactly h 1,1 (X) + 4 prime toric divisors D A . We call the inherited prime toric divisors on X.

Inherited and autochthonous divisors
The inherited prime toric divisors D A of a Calabi-Yau threefold hypersurface X ⊂ V are effective divisors on X that are inherited from effective divisors on V . We call the cone in H 4 (X, R) generated by the classes of the {D A } the inherited effective cone Eff ι (X). Clearly, Eff ι (X) ⊆ Eff(X). Because V is specified by combinatorial data, it is straightforward to compute Eff ι (X). However, in many cases Eff ι (X) Eff(X): that is, there are effective divisors on X that are not inherited from any effective divisor on V . We call such a non-inherited divisor an autochthonous divisor.
In this work, we approximate Eff(X) by Eff ι (X). In particular, in computing axion masses in compactifications of type IIB string theory on X, we here consider only Euclidean D3-branes wrapping inherited effective divisors. Autochthonous divisors of Calabi-Yau hypersurfaces are studied in our forthcoming work [8]. Among other things, we show there that Euclidean D3-branes wrapping autochthonous divisors do not sig-nificantly affect the axion mass hierarchies found here, and so for present purposes it suffices to study the conceptually and computationally simpler inherited effective cone.

The Kähler cone and the Mori cone
Let X be a projective algebraic variety, and let J ∈ H 1,1 (X, R) be a closed (1,1)-form on X. For a k-dimensional subvariety U ⊂ X, we define We define the Kähler cone of X, K X , as the subset of H 1,1 (X, R) consisting of cohomology classes of Kähler forms J on X, i.e. J such that Vol J (U ) > 0 for all subvarieties U . The Kähler cone K X , also called the ample cone, is an open convex cone whose closure K X is the cone of nef (1,1) classes. 4 We next define the Mori cone of X, M X , to be the cone in H 2 (X, R) generated by irreducible algebraic curves C a on X. (The Mori cone of X is often denoted NE(X) in other parts of the literature.) The Kähler cone and the Mori cone are related by i.e. the dual of the Mori cone is the closure of the Kähler cone. When X is a Calabi-Yau threefold hypersurface, the subvarieties of interest are the curves C a , the divisors D A , and the threefold itself. The volumes of these subvarieties are t a :=Vol J (C a ) = Ca J , It is convenient to expand J in terms of the Poincaré duals [D i ] of the divisors D i , (2.10) The h 1,1 Kähler parameters t i , which are not necessarily positive when J is inside the Kähler cone, should not be confused with the curve volumes t a , which are positive for J ∈ K X .

The Stretched Kähler Cone
One of the aims of this work is to determine the volumes of holomorphic submanifolds in X, when every effective curve in X is constrained to have volume > 1, as in (1.2). We therefore need to determine the cone of effective curves, i.e. the Mori cone M X . To our knowledge there is no established algorithm for computing M X , even for the well-studied ensemble of Calabi-Yau threefold hypersurfaces. However, we will identify two cones M in and M out that bound M X on the inside and the outside, and it is these bounding cones that we will study. The duals of these cones will then provide cones that bound K X on the outside and the inside, respectively: defining K in := M ∨ out and K out := M ∨ in , and writing K in for the interior of K in , and K out for the interior of K out , we have As we shall see, the Kähler cone K V of the ambient toric variety V can play the role of K in , while a new cone, K ∩ , provides the outer bound K out [10].
Although computing M X is challenging, the Mori cone M V of the toric variety V can be computed efficiently from the fan using an algorithm due to Berglund, Katz, and Klemm [11], which is equivalent to the classical algorithm of Oda and Park [12]. By (2.6), the dual of M V is the closure K V of the Kähler cone K V of V . Restricting the Kähler parameters t i so that t i [D i ] ∈ K V ensures that all holomorphic submanifolds of V have positive volume, and therefore this restriction also guarantees that all holomorphic submanifolds of X have positive volume. We therefore have We remark that subvarieties of V that correspond to simplices interior to facets do not intersect a generic X, and therefore any triangulations of ∆ • that differ only by simplices interior to facets define isomorphic Calabi-Yau hypersurfaces, but with different toric ambient spaces V α . It is then natural to glue the Kähler cones K Vα together and define K ∪ [7]: However, such a process appears prohibitively complicated at large h 1,1 , and will not play a role in our analysis.
Consider the following set of surfaces in V : The intersection of any of the S AB with a generic anticanonical hypersurface X, if nonempty, is transverse and defines a corresponding curve in X, The curve C AB lies in M X , but in general not every element of M X can be written in the form (3.6). Because the {C AB } are the curves inherited from intersections of distinct prime toric divisors, we call the {C AB } toric intersection curves. The volumes of the toric intersection curves are We define the intersection cone K ∩ as the space of Kähler parameters t i for which the volumes V, τ A and t AB are all positive: As these conditions are necessary, but in general not sufficient, 5 to ensure that t i [D i ] ∈ K X , we have the inclusions The stretched Kähler cone: In order to study the effect of demanding that all cycles satisfy the minimal volume constraint (1.2), we introduce the notion of a stretched Kähler cone. Let X be a projective algebraic variety, let J ∈ H 1,1 (X, R) be a closed (1,1) form on X, and let W = {W } be a set of subvarieties W ⊂ X. Given a number c > 0, we define the (c, W)-stretched Kähler cone of X, (3.10) The first stretched Kähler cone we consider is the stretched Kähler cone of X, We next define the stretched intersection cone as the region in which all intersection curves C AB , all inherited prime toric divisors D A , as well as X itself, have volume ≥ 1. In all cases we have K X ⊆ K ∩ , but because the curves C AB typically do not generate M X , we typically have K X K ∩ . Finally, noting that for favorable X, H 1,1 (V, R) can be naturally identified with H 1,1 (X, R), we define the stretched Kähler cone of V, In a complete toric variety, any curve is rationally (and thus numerically) equivalent to an effective sum of toric curves [13]. A curve C ⊂ X ⊂ V is also a curve in V , and so in homology C can be expressed as a non-negative integral linear combination of toric curves. It follows that K V ⊂ K X .
We have therefore bounded the stretched Kähler cone: (3.14)

Axion Couplings
Consider a compactification of type IIB string theory on an orientifold 7 of a Calabi-Yau threefold hypersurface X. The four-dimensional theory contains h 1,1 axions from reduction of the Ramond-Ramond four-form C 4 . In this section we explain how the kinetic and potential couplings of the axion fields are computed from geometric data.

Kinetic term
In terms of a basis {D i }, i = 1, . . . , h 1,1 for H 4 (X, Z), we define to be the corresponding dimensionless axions. The Kähler coordinates on Kähler moduli space are the complexified divisor volumes . The axion kinetic term is then 8 where the Kähler metric K ij is obtained from the Kähler potential K = −2 log V.

Nonperturbative superpotential
The axions are perturbatively massless and receive mass only nonperturbatively, from instantons: specifically, from Euclidean D3-branes wrapping four-cycles. 9 The leading-order bosonic action S for a Euclidean D3-brane is given by the Dirac-Born-Infeld action plus an imaginary Chern-Simons term that provides the coupling to the axion (see e.g. [14]). Consider Euclidean D3-branes wrapping the four-cycles for some n i α ∈ Z, α = 1, . . . , N , and for some N > 0. The action S α of the Euclidean D3-brane wrapping Σ α is then Although one can in principle consider Euclidean D3-branes wrapping any four-cycle Σ α ∈ H 4 (X, Z), the situation is best-understood when Σ α is an effective divisor, i.e. when [Σ α ] ∈ Eff(X): precisely in that case, Σ α is calibrated by the Kähler form J, and so obeys , determining the volume of the minimum-volume representative of the class [Σ α ] is in general very difficult, as we explain in Appendix A. Moreover, Euclidean D3-branes wrapping a representative Σ α of a class [Σ α ] ∈ Eff(X) cannot contribute to the superpotential. They may contribute to the Kähler potential, but such effects are not well understood.
For now we will focus on effective divisors, and we suppose that superpotential terms arise from Euclidean D3-branes wrapping the divisors for some q i α ∈ Z, α = 1, . . . , p, and for some p > 0. The superpotential then takes the form [4,15] W where W 0 is the classical flux superpotential [16]. The Pfaffians A α depend on the complex structure moduli, and will be set to unity in our analysis. The axion potential can then be written as We will make the conservative choice W 0 ∼ 1: a smaller value of the flux superpotential would make our upper bounds on axion masses more stringent. Performing a The Hessian of the canonically-normalized axions is and we denote its eigenvalues by h 2 1 ≤ . . . ≤ h 2 h 1,1 . The potential (4.9) has a rich structure of minima and critical points, cf. e.g. [17][18][19], and finding the global minimum numerically is expensive when h 1,1 1 and p h 1,1 (for p slightly larger than h 1,1 , which does not hold here, the methods of [19] could be used). In the remainder, by axion masses-squared we mean the Hessian eigenvalues h 2 i , evaluated at the origin 0, i.e. at θ 1 = θ 2 = · · · = θ h 1,1 = 0. By minimum axion mass-squared we mean One should bear in mind that these quantities could change slightly upon shifting the axion vev to a minimum, but we have found no evidence for changes large enough to invalidate our parametric results.

Geometric field ranges
The effective Lagrangian for the axions is usefully rewritten as where the mass scales Λ a are determined by the instanton actions S α , and the charge matrix Q has the entries where a = 1, . . . p(p + 1)/2 ≡ P . The rows involving q i β − q i γ arise from cross terms in the F-term potential, see [5,20].
Because the potential is periodic it is natural to define the axion fundamental domain F [5,20], given by the hyperplane constraints: The fundamental domain is compact when Q has rank h 1,1 . A quantity of key interest for axion inflation is the geometric field range, i.e. the maximum distance R from the origin to the boundary of F, measured with respect to K ij . That is, where d ρ is the matrix of the vertices of F, and K is the Kähler metric. The walls and vertices of F are determined by the integers q i α , i.e. by the set of effective divisors {D α } in (4.7) that support superpotential terms. The problem of identifying those effective divisors of a Calabi-Yau threefold hypersurface that support nonvanishing superpotential terms has not been fully solved, cf. [21]. For the purposes of the present work we will assume that every prime toric divisor D A supports a Euclidean D3-brane superpotential term, cf. [5].
Computing R directly from (4.16) is prohibitively expensive at large h 1,1 , since the number of vertices that must be checked is at least 2 h 1,1 −1 . We will instead consider an upper bound on R. By performing a basis transformation where Q is a rank h 1,1 subblock of Q, we can trivialize 2h 1,1 of the hyperplane constraints. The metric in the ϑ basis is then with eigenvalues ξ 2 1 ≤ · · · ≤ ξ 2 h 1,1 . An upper bound for R is then given by where ξ 2 h 1,1 is the largest eigenvalue of Ξ. When Q is not square, R bound depends on the choice of Q, but each choice does provide an upper bound on R. Because we have assumed that each of the D A supports a Euclidean D3-brane superpotential term, we can choose h 1,1 of the toric coordinates for Q, in such a way that Q is the h 1,1 × h 1,1 identity, and Ξ = K.

Masses
Suppose that, for some specified Kähler form J, and for some positive number L, every basis of H 4 (X, Z) contains at least k ≥ 1 members Σ α with Vol(Σ α ) > L in string units. Then at least k axions must have mass e −2πL . One can therefore place upper bounds on the masses of the lightest axions by placing lower bounds on the volumes of four-cycles furnishing bases for H 4 (X, Z).
Let us first consider placing upper bounds on superpotential contributions to axion masses, by placing lower bounds on the volumes of effective divisors. As explained in §2.1.2, in this work we approximate Eff(X) by Eff ι (X); corrections to this approximation will be described in [8]. The inherited prime toric divisors {D A }, A = 1, . . . , h 1,1 + 4, provide a set of generators of Eff ι (X), and also, in the above approximation, of Eff(X). For any J ∈ H 1,1 (X, R), not necessarily inside K X , we can compute the volumes τ A := 1 2 D A J ∧ J. There are at most h 1,1 +4 sets {D i } of h 1,1 prime toric divisors that furnish bases for H 4 (X, Q), and for each such basis B we can compute the volumes τ B 1 ≤ · · · ≤ τ B h 1,1 of the basis generators. Define B min to be the basis choice that minimizes τ B h 1,1 . Roughly speaking, B min is a minimum-volume basis of generators of the effective cone. We write denoting explicitly the dependence on the choice of J. We can now give an upper bound on the magnitude of the leading superpotential term involving the lightest axion, for the given J: Furthermore, given any region R ⊂ H 1,1 (X, R), not necessarily inside K X , we can compute τ R last := min J∈R τ last (J) . We then write Using (3.14), we have τ ∩ last ≤ τ X last ≤ τ V last . Thus, when the condition (1.2) for control of the α expansion is imposed, the superpotential for the lightest axion is bounded above by For h 1,1 1 the exponentials in (4.9) are parametrically dominant, and in evaluating the dependence of (4.9) on the lightest axion we can omit factors that are only polynomial in the volumes, including the effect of canonical normalization. 10 We then arrive at an upper bound on the mass-squared m 2 min of the lightest axion from (4.26), What about axion mass terms from instanton contributions to the Kähler potential, resulting from Euclidean D3-branes wrapping classes [Σ α ] ∈ H 4 (X, Z) that are outside Eff(X), and admit no holomorphic representative? Could such instantons give masses m min ? We discuss this question in Appendix A, and find that present knowledge of minimum-volume representatives of classes outside Eff(X) is not sufficient to give a definite answer, but at the same time there is no evidence of such a parametric enhancement in known threefolds. We find it plausible that masses from K are least parametrically comparable to those from W , and so are approximately given by (4.27).

Computation
We obtained the topological data of Calabi-Yau threefold hypersurfaces as follows. For each value of h 1,1 that we studied, we drew a number N (h 1,1 ) of polytopes at random from the Kreuzer-Skarke database: see Table 1. 11 We manipulated the polytopes using Sage [22]. For each polytope that was favorable, we used TOPCOM to obtain a fine and regular (but not star) triangulation. We removed the lines in the strict interior of the polytope and included a line from the origin to each point in the polytope, thus producing an FRST T [23]. Such a triangulation defines a fan, and in turn defines a toric variety V . As explained in §2.1.1, to study a generic Calabi-Yau threefold hypersurface, one can omit simplices of T that pass through facets of ∆ • . We denote the set of remaining simplices by T , and abuse language slightly in calling T an FRST as well.
Because we have restricted to favorable polytopes, there are h 1,1 + 4 prime toric divisors D A ⊂ V , each corresponding to a ray of the fan determined by T . We picked a basis for H 4 (X, Q) by selecting a set of h 1,1 of the inherited prime toric divisors D A := D A ∩ X that are linearly independent. Using Sage, we computed the triple intersection numbers κ ijk in the chosen basis. Finally, we computed the Mori cone M V of the toric variety in Mathematica using the algorithm described in [11], which is equivalent to that of [12], but easier to implement.
With this data in hand, we turned to analyzing the resulting cones. For each geometry the stretched cones K V and K ∩ were constructed as described in §3. We minimized the volumes τ A and V inside K V and K ∩ using IPOPT, a software package for largescale nonlinear optimization, which is included in version 11 of Mathematica. Because IPOPT uses an interior point algorithm that finds a local solution to the optimization problem, we performed the minimization multiple times, from different starting points, in an attempt to find the global minimum. Finding even one feasible starting point for the optimization algorithm is challenging at large h 1,1 , as the cones K V and K ∩ become very narrow. We made use of IBM's optimization software CPLEX as well as the LinearProgramming function of Mathematica to find such points.
Note that we computed one FRST for each favorable polytope studied. With our methods it takes of order a day to obtain the topological data of all FRSTs of all threefolds with h 1,1 ≤ 6, but for larger h 1,1 it quickly becomes infeasible to compute all triangulations. In order to provide a better point of comparison for the data we can obtain at h 1, 1 1, we limited ourselves to one FRST per polytope even for small h 1,1 . The values of h 1,1 that we studied, and the numbers N (h 1,1 ), were chosen to balance the computational expense at h 1,1 1 against the potential for illuminating scaling laws. Obtaining more extensive data at large h 1,1 is an obvious next step [8]. In fact, the present work has established the feasibility of obtaining the topological data of at least one threefold (i.e., one FRST) for each polytope in the Kreuzer-Skarke database. A very rough estimate is that such a computation could require a few CPUcenturies, absent any improvements to the algorithms.

Results
The primary topological data produced by our analysis are the generators of the Mori cones M V of toric varieties V , and the intersection numbers κ ABC of inherited prime toric divisors D A of Calabi-Yau hypersurfaces X ⊂ V . Taking these data and imposing the condition (1.2), we can compute the stretched Kähler cones K V and K ∩ , which bound the stretched Kähler cone K X of X from the inside and the outside, respectively, cf. (3.14). Then, for any holomorphic 2k-cycle Σ 2k (1 ≤ k ≤ 3) in X, K ∩ determines a lower bound on Vol(Σ 2k ) ≡ 1 k! Σ 2k ∧ k J. In this section we report salient features of the intersection numbers, Mori cones, volumes of holomorphic cycles, geometric field ranges, and masses in our ensemble.

Topological Data:
The volumes V, τ A and t AB depend on the intersection numbers κ ijk , as given in (2.7). Since κ ijk depends on a choice of basis of H 4 (X, Z), we instead report the statistical properties of κ ABC , which is basis-independent.
We first examine the sparseness of κ ABC . The number of nonvanishing intersection numbers per geometry increases approximately linearly with h 1,1 , as shown in Figure 1. As a result, κ ABC becomes very sparse at large h 1,1 . In Figure 2 we show the root mean square (RMS) size of the nonvanishing intersection numbers for each geometry.
The cone K V is given by the intersection of the half-spaces defined by the linear inequalities M ai t i > 0. (6.1) As h 1,1 increases, the number of inequalities grows and K V becomes very narrow. A conceptually straightforward way to quantify the narrowness of the cone K V would be to analyze the behavior of the solid angle subtended by K V as a function of h 1,1 . However, this becomes computationally expensive when h 1,1 15. Instead, we characterize the narrowness of K V by computing the cosine of the smallest angle between two hyperplanes, denoted M a and M b : As the angle θ min between two hyperplanes approaches zero, the cone becomes infinitely narrow. This can also be understood from the perspective of the dual cone M V . When M V has two generating rays M a and M b that are almost antiparallel (such that K V has facets whose normals are almost antiparallel), it is difficult to find a Kähler form J such that both of the associated curves have positive volumes simultaneously. Figure  3 shows cos(θ min ) as a function of h 1,1 .     As K V becomes more narrow, the stretched cone K V , defined in (3.13), gets pushed further away from the origin. Another measure of the size of K V is therefore the shortest distance d V min between the origin and any point of and we denote the minimum-distance point by t V d . See Figure 4. Although K V is computationally accessible (even for h 1,1 = 491), and the size of K V is generally correlated with the size of K X , K V can in principle be much more narrow than K X . Analysis of K V alone can therefore provide only estimates of the volumes of holomorphic cycles in X, for t i [D i ] ∈ K X , rather than definite bounds. To obtain lower bounds on cycle volumes, we instead examine K ∩ , which contains K X . The tradeoff is that K ∩ is a complicated cone defined by linear, quadratic and cubic constraints, and defining a quantity analogous to θ min is difficult. We can, however, compute d ∩ min , the shortest distance between the origin and any point of K ∩ , and we denote the minimum-distance point by t ∩ d . See Figure 5.

Volumes:
To compute lower bounds on τ last and V, for each prime toric divisor D A we numerically minimize the divisor volume τ A in K V and in K ∩ . We then calculate τ V last and τ ∩ last as described in §4 and §5. The resulting bounds are shown in Figures 6-9. 12 (a) log 10 (τ V last ) vs. log 10 (h 1,1 ). The fit is log 10 (τ V last ) = −1.9 + 4.3 log 10 (h 1,1 ).

Geometric field ranges:
As explained in §4, we estimate the radius R (4.16) of the axion fundamental domain F by assuming that all prime toric divisors D A contribute to the superpotential. The radius depends on the Kähler parameters t i , and we report upper bounds on R at two locations. We define t V V and t ∩ V to be the points in K V and K ∩ , respectively, where the threefold volume V is minimized, and we define We first compute the Kähler metric K ij at t V V . We next trivialize 2h 1,1 of the hyperplane constraints, as in (4.17), taking Q to be the h 1,1 × h 1,1 identity subblock of the charge matrix Q corresponding to a choice of h 1,1 of the toric coordinates. This yields an upper bound R V bound ≥ R V , shown in Figure 10. Computing K ij instead at t ∩ V and repeating the trivialization, we obtain the upper bound R ∩ bound ≥ R ∩ , shown in Figure  10.

Axion masses:
Now consider type IIB string theory compactified on an orientifold of a hypersurface X from our ensemble. The large divisor volumes lead to powerful suppression of superpotential contributions to the potential for C 4 axions θ i . We find that in every geometry in our ensemble 13 with h 1,1 > 22, the lightest axion is essentially massless, with the canonically-normalized field having mass m < 10 −33 eV .  Let us also give a heuristic estimate of the expected mass of the lightest axion. By the definition (4.22), every basis for H 4 (X, Q) consisting of elements of Eff(X) has members with volume ≥ τ ∩ last . As seen from the fit in Figure 6, τ ∩ last ∼ 0.02(h 1,1 ) p with p ∼ 3. Hence, one or more of the h 1,1 axions θ i receives no superpotential contributions larger than The exponent p changes, within the range 3 p 6, depending on whether one examines τ ∩ last -which is the most direct and conservative -or instead a more computable proxy such as τ V last or (d V min ) 2 . However, such changes do not alter our central finding that one or more axions are extremely light when h 1,1 1 and J ∈ K ∩ .

Summary:
A root cause of our findings is that the Kähler cones of Calabi-Yau threefold hypersurfaces are very narrow for h 1,1 1, as shown in Figure 3. The condition (1.2) that every effective curve has volume ≥ 1, which we have used as a proxy for control of the α expansion, then implies that the Kähler form J ∈ H 1,1 (X, R) is far from the origin in H 1,1 (X, R), in the sense of (6. Furthermore, the volume V of X itself is large (Figures 8-9), the geometric field range is generally small (Figure 10), and the eigenvalues of the axion kinetic matrix are small ( Figure 11). The minimum axion mass is small, and strongly correlated with τ ∩ last ( Figure 12). 14

Implications for Axion Cosmology
The overall picture that emerges from our analysis is that in a compactification of type IIB string theory on an orientifold of a Calabi-Yau threefold hypersurface X with h 1,1 1, in the regime of control of the α expansion, X and most of its subvarieties have very large volumes in string units. The resulting effective theory has many axions, some of which are essentially massless, 15 with m 10 −33 eV. The axion kinetic matrix has small eigenvalues, and the radius of the axion fundamental domain is sub-Planckian. 16 In summary, we find an axiverse with hundreds of axions, some of them 14 In fact, m min (t ∩ d ) is almost perfectly correlated with τ last (t ∩ d ). Note that by (4.20), τ ∩ last = τ last (t ∩ d ). 15 Many authors use the term 'ultralight axion' for axions with m 10 −33 eV that could make up part of the dark matter, as in [24,25]. We avoid the term 'ultralight' when speaking of the far lighter axions found here, with m 10 −33 eV; these we instead call 'massless', even though strictly speaking their masses are negligibly small, not zero. 16 As explained in §6, in a small fraction of cases we cannot exclude the possibility of super-Planckian radii, but neither can we prove that all curves in X have positive volume in these cases. For the present discussion we consider only the better-established examples with J ∈ K V , for which the radii are sub-Planckian. massless, and all with small periodicities. In this section we will mention a few of the implications of these findings for the cosmology of compactifications with h 1, 1 1. Axions and axion-like particles with appreciable couplings to the Standard Model are strongly constrained by a wealth of data from diverse channels, including terrestrial appearance experiments such as helioscopes and haloscopes; red giant evolution; supernovae; CMB spectral distortions; and X-ray production in galactic or cosmological magnetic fields. See [26] for a review. To apply these constraints to the large-h 1,1 axiverse that we have described, it would be necessary to make specific assumptions about the realization of the Standard Model, and its couplings to the axion sector. While very interesting, such an analysis is extremely model-dependent.
Cosmological effects of the gravitational couplings of axions present a more direct, but still somewhat model-dependent, set of constraints. Sufficiently light axions, with m 10 −33 eV, are indistinguishable from vacuum energy unless excitations of the axion field (i.e., particles) are produced as dark radiation, for example through the decay of an associated modulus. The limits on dark radiation are rather stringent, cf. [27][28][29], but again depend on the details of post-inflationary evolution. For example, if the lightest modulus decays only to a single axion, as well as to Standard Model particles, the dark radiation constraints are insensitive to the existence of other axions and moduli [29], but can be severe nonetheless [27,28]. Axions with m ∼ 10 −33 eV can be quintessence fields [30], and in special cases could even alleviate the "why now" problem [31]. Axions with m ∼ 10 −22 eV could constitute a portion of the dark matter, and could give rise to small-scale structure in better agreement with observations than that predicted by cold dark matter models [24] (for recent work, see e.g. [25,[32][33][34]). Overproduction of axion dark matter -and in some mass ranges, isocurvature perturbations in the CMB -provide serious constraints [35][36][37], especially in models with many axions [38].
Perhaps the most interesting constraints on the large-h 1,1 axiverse come from black hole superradiance [39]. Axions in the mass range 10 −10 eV − 10 −20 eV, even if not present as a cosmologically abundant population, can trigger instabilities of black holes. Detailed modeling of moduli stabilization would be necessary to make precise statements, but as a rough estimate, we find that approximately half the geometries in our ensemble have an axion in the mass range 10 −10 eV ≤ m ≤ 10 −20 eV. Superradiance limits on many-axion theories have been obtained in [40]. However, the analysis of [40] is only directly applicable to theories with relatively large decay constants, f 10 14 GeV. Axions with smaller periodicities suffer from nonlinear interactions, potentially changing the limits of [40]. A dedicated study of superradiance constraints on the Kreuzer-Skarke axiverse would be a worthwhile topic for the future.

Conclusions
We have initiated a survey of compactifications on Calabi-Yau threefold hypersurfaces with arbitrary h 1,1 , i.e. of the entire Kreuzer-Skarke list.
This work extends and complements the complete enumeration carried out by Altman et al. [2] for h 1,1 ≤ 6. The large h 1,1 regime presents evident computational challenges, a few of which we have overcome. Publicly-available software such as Sage [22] generally produces FRSTs only for h 1,1 10, and the improved triangulation algorithms that have been implemented on a large scale in the past are expensive, and function only for h 1,1 30 [2,41]. Moreover, the sheaf cohomology computations needed for studying divisors D ⊂ X likewise explode in difficulty for h 1,1 10. These limitations have led to a perception that systematic enumeration and study of hypersurfaces with h 1, 1 10 -corresponding to the bulk of the Kreuzer-Skarke database -is not possible at present. In this work we have demonstrated, on the contrary, that large-scale studies are feasible for any range of h 1,1 arising in the Kreuzer-Skarke list, given only modest computational resources.
A key step was implementing the triangulation algorithm described in [23], which allowed us to obtain fine regular star triangulations ∼ 5 · 10 3 times faster (per CPU) than was possible in [2,41]. With our methods, finding one FRST takes just seconds even for h 1,1 = 491. 17 However, although we can efficiently generate large numbers of compactifications at any desired h 1,1 , several challenges remain. In this work, we used Sage to obtain the intersection numbers of hypersurfaces with h 1,1 ≤ 100, at a computational cost of order half a CPU-hour per hypersurface at h 1,1 = 100. Significant further gains are possible in this area, and allow efficient computation of intersection numbers for any h 1,1 ≤ 491, as we will show in [8]. Even so, one thing that remains out of reach is a complete enumeration of hypersurfaces at large h 1,1 , simply because the number of such hypersurfaces -corresponding to the number of inequivalent triangulations of reflexive polytopes with many lattice points -appears to be vast.
The principal raw data produced by our analysis are FRSTs of four-dimensional reflexive polytopes; the Kähler cones of the corresponding toric varieties V ; and the intersection numbers of generic Calabi-Yau threefold hypersurfaces X ⊂ V . These data provide a wealth of information about four-dimensional effective theories arising from string compactifications on such X. In this paper we studied two of the most salient physical findings, axion mass hierarchies and axion field ranges, leaving a complete characterization of the physical implications of our topological and geometric data as a significant task for future work.
The first observable we studied is the set of relations among cycle volumes enforced by the Mori cone conditions, which control the structure of the potential generated by instantons. We found that enforcing that every effective curve has volume at least one in string units, as a proxy for ensuring control of the α expansion, has -for h 1,1 1 -a striking consequence: the volumes of many irreducible curves and divisors on X, and of X itself, become extremely large. We found that these volumes scale roughly as (h 1,1 ) p , with the exponent 3 p 7 depending on the type of cycle considered.
One consequence is that in a compactification of type IIB string theory on an orientifold of a typical Calabi-Yau threefold hypersurface with h 1, 1 1, one of the following holds: 1. One or more axions are effectively massless.
3. The axion mass terms produced by Euclidean D3-branes wrapping non-holomorphic four-cycles are parametrically larger than those from holomorphic four-cycles.
When condition (3) holds, the breakdown of the α expansion cannot be detected by computing the volumes of calibrated cycles, while condition (2) suggests but does not guarantee the existence of large perturbative and nonperturbative corrections in the α expansion. Thus, we have established a tension between nonvanishing masses for all axions, and manifest control of the α expansion. The second observable we studied is the metric on Kähler moduli space, which is relevant for understanding quantum gravity constraints on large-field inflation. We found that the eigenvalues of the axion kinetic matrix are typically small at large h 1,1 , primarily because of the large volume of X. In each geometry we computed an approximation to the radius of the axion fundamental domain. The radius depends strongly on how restrictive a condition one imposes on the Kähler form J. For Kähler forms such that every curve in the ambient toric variety V has volume ≥ 1, we found field ranges M pl in every example. In the less restrictive case of Kähler forms in the region defined by (3.12), corresponding to the outer approximant to the stretched Kähler cone of X, we found super-Planckian axion field ranges in a small fraction of geometries, at each h 1,1 . While intriguing, this finding cannot be taken as evidence for large field ranges in the regime of control of the α expansion, because without a direct computation of K X we cannot exclude the possibility that in each example giving an apparent large field range, one or more effective curves C ∈ M X has volume < 1, or indeed < 0. 18 Overcoming this limitation is an important task for the future.
Because our results are drawn from a statistical study of an ensemble of geometries, they should be taken as statements about typical compactifications, and there can be special Calabi-Yau threefolds that violate the behavior observed in this work.
For instance, the Mori cone constraints in the geometry studied in [42] are much milder than those of a typical hypersurface, possibly because the Kähler moduli space in [42] has a symmetry corresponding to exchanging the Kähler moduli that correspond to blowups. It would be interesting to understand the incidence of such geometries.
Our results give a sharper picture of the spectrum of axion masses and decay constants arising in geometric compactifications of string theory. Reasonable a priori estimates of these spectra, as well as studies in families of examples, have been made and used in the study of the string axiverse [29,39,[43][44][45], and our ensemble provides a foundation for refining these estimates.
There are several directions for future work. By applying computational resources on a larger scale, one could more finely sample the Kreuzer-Skarke database. A rough estimate is that in under a few hundred CPU-years one could find one FRST for every one of the 473,800,776 polytopes in the list. Intersection numbers could be obtained at comparable cost given the improved methods of [8]. As we have stressed, however, the number of distinct triangulations is plausibly vastly larger than the number of polytopes, and so it remains to be seen whether it is possible even to store the topological data of every compactification arising from the Kreuzer-Skarke database.
The geometric data obtained here can serve to answer questions about which sorts of effective theories are possible in compactifications on Calabi-Yau hypersurfaces. To answer such questions, it would be natural to use machine learning [46][47][48][49][50][51][52][53], among other tools, given the scale and complexity of the data.

A Non-holomorphic Instantons
We noted in §2 that computing contributions to the superpotential from Euclidean D3branes wrapping holomorphic four-cycles (i.e., effective divisors) is much simpler than computing contributions to the Kähler potential from Euclidean D3-branes wrapping non-holomorphic four-cycles. On the physics side, one reason for the disparity is that superpotential terms are constrained by holomorphy. Geometrically, the difference between the two computations is that a holomorphic four-cycle Σ 4 is calibrated by the Kähler form J, and enjoys the relation so that once J is given, Vol(Σ 4 ) is determined by topological data. Similarly, an antiholomorphic four-cycle Σ 4 has orientation opposite to that of a holomorphic cycle, and obeys However, it is much more difficult to compute the volume of a cycle that is neither holomorphic nor antiholomorphic, as we now explain.

A.1 Volume-minimizing currents
Suppose that X is a compact Kähler manifold of dimension n, 19 i.e. MinVol([α]) is the volume of the smallest-volume representative of the class α. However, it is not obvious that the variational problem implied by (A.3) is wellposed: does one search over all representatives of α, or just representatives obeying an appropriate smoothness condition? It is also not clear a priori how smooth the volumeminimizing configuration will be: in fact, one can easily find examples in which the volume-minimizing configuration has singularities, at least at complex codimension one.
Fortunately, the problem of finding the minimum-volume representative of a given homology class is one of the central questions of geometric measure theory, and was put on sound footing in the 1960s by Federer and Fleming. They defined objects called integral p-currents, which roughly correspond to formal sums of p-dimensional submanifolds, except for sets of p-dimensional Hausdorff measure zero. Federer and Fleming showed that the class of integral p-currents has a compactness property that is very useful in formulating variational problems: in fact, they proved that for each class 20 [α] ∈ H 2n−2 (X, Z), there exists an integral current of least volume [54]. In other words, (A.3) actually does define the solution of a well-posed variational problem, provided that α is understood to vary over integral currents, not just over smooth submanifolds.

A.2 Non-holomorphic instantons and volume reduction
Consider a Euclidean D3-brane in a homology class [α] ∈ Eff(X), which necessarily cannot contribute to the superpotential, but may contribute to the Kähler potential. The real part of the action of such a Euclidean D3-brane is plausibly proportional to MinVol([α]), which is well-defined thanks to geometric measure theory. Even so, computing MinVol([α]) is nontrivial.
As a toy example, suppose that X is such that four-cycles α 1 and α 2 are a basis for H 4 (X, Z) = Z 2 , and α 1 and α 2 also generate Eff(X). For a given Kähler form J, define τ i := 1 2 α i J ∧ J and θ i := α i C 4 , i = 1, 2. If J is such that τ 1 , τ 2 L for some L 1, then Euclidean D3-brane terms in the superpotential are no larger than e −2πL .
In this situation, one should ask about contributions to the Kähler potential from Euclidean D3-branes wrapping a representative γ of a non-effective class such as [γ] := [α 1 − α 2 ] ∈ Eff(X). Because α 1 and α 2 are calibrated by J, we have MinVol([α i ]) = τ i L. The action of instantons on γ is determined by τ γ := MinVol([γ]). However, we cannot conclude that τ γ ≥ τ 1 + τ 2 . (A.4) If α 1 and α 2 are disjoint, then (A.4) actually does hold, but more generally the intersection locus of the minimum-volume representatives of [α 1 ] and [α 2 ] can be deformed to produce a representative of [γ] with volume < τ 1 + τ 2 . When τ γ = τ 1 + τ 2 − ∆τ for ∆τ > 0, we will say that recombination has led to volume reduction by an amount ∆τ > 0. The question of volume reduction is best-understood for two-dimensional currents. Building on work of Almgren [55], Chang proved that in any Riemannian manifold, the singular set of a volume-minimizing two-dimensional current consists of isolated branch points [56]. It is therefore tempting to conjecture that in a Kähler manifold, a volume-minimizing two-dimensional current consists of a union of holomorphic and antiholomorphic curves, intersecting at points; one consequence would be that there is no volume reduction for two-dimensional currents in a Kähler manifold. However, in [57], for X a K3 surface, Micallef and Wolfson gave an explicit example of a noneffective class [α 1 − α 2 ] ∈ H 2 (K3, Z) whose minimum-volume representative is not a union of holomorphic and antiholomorphic curves, 21 and for which MinVol([α 1 − α 2 ]) < MinVol([α 1 ]) + MinVol([α 2 ]). The volume reduction is proportional to the small parameter ε measuring the deviation from the orbifold limit of K3.
The issue, returning to four-cycles, is then the following. If a significant volume reduction ∆τ ∼ O(τ 1 , τ 2 ) could occur in some setting, so that τ γ τ i , then ensuring τ 1 , τ 2 L would not place any upper bound on the size of Euclidean D3-brane terms in the Kähler potential. The axion masses from non-holomorphic instantons in K would be parametrically larger than those from holomorphic instantons in W .
Although the Micallef-Wolfson construction proves that nonzero volume reduction can occur in a Calabi-Yau compactification, we are not aware of any example of parametrically large volume reduction in a comparable setting. Moreover, the cycle volume determines only the leading semiclassical action of a Euclidean D-brane, and one should compute corrections to this action, such as the fluctuation determinant, before drawing conclusions about the relative sizes of physical effects. 22 In summary, determining whether Euclidean D3-branes wrapping non-holomorphic cycles can contribute axion masses that are parametrically larger than those arising from holomorphic cycles is an open problem. The available evidence does not exclude this possibility, but also does not, in our view, strongly support it. Our results on axion masses rely on our computation of the volumes of holomorphic cycles, and could be affected if large volume reduction occurs and causes non-holomorphic instantons to dominate in the potential. This proviso should be kept in mind when interpreting our findings.