Heterotic line bundle models on elliptically fibered Calabi-Yau three-folds

We analyze heterotic line bundle models on elliptically fibered Calabi-Yau three-folds over weak Fano bases. In order to facilitate Wilson line breaking to the standard model group, we focus on elliptically fibered three-folds with a second section and a freely-acting involution. Specifically, we consider toric weak Fano surfaces as base manifolds and identify six such manifolds with the required properties. The requisite mathematical tools for the construction of line bundle models on these spaces, including the calculation of line bundle cohomology, are developed. A computer scan leads to more than 400 line bundle models with the right number of families and an SU(5) GUT group which could descend to standard-like models after taking the ℤ2 quotient. A common and surprising feature of these models is the presence of a large number of vector-like states.


JHEP04(2018)087
1 Introduction Heterotic string model building on Calabi-Yau (CY) manifolds requires an understanding of the gauge bundle on the compactification space and this constitutes one of the main technical challenges in constructing realistic models, particularly for gauge bundles with a non-Abelian structure group. Indeed, only a relatively small number of quasi-realistic models based on such gauge bundles are known in the literature [1][2][3][4]. However, it has been realised [5,6] that quasi-realistic models can also be constructed using gauge bundles with Abelian structure groups which are considerably easier to handle than their non-Abelian counterparts. For these models, the vector bundle is a direct sum of line bundles which are classified and can be analyzed systematically. By scanning this space of heterotic line bundle models, a large number of examples which lead to the standard model spectrum has been found [5,6]. For technical reasons, these models have been obtained using the arguably simplest type of CY manifolds -complete intersections in product of projective spaces -and, to a lesser extent, CY hypersurfaces in toric four-folds. The main purpose of the present paper is to study heterotic line bundle models for another class of CY manifolds, namely elliptically fibered CY three-folds.
There are a number of motivations for addressing this problem. Firstly, we would like to develop the necessary tools for constructing line bundle models on elliptically fibered CY three-folds, including the calculation of line bundle cohomology on these manifolds. Secondly, by studying line bundle models on another class of CY manifolds, we would like to gain some insight into which of their properties are generic and which are related to the particular type of underlying CY manifold. Finally, we are motivated by heterotic F-theory duality which is formulated for elliptically fibered CY manifolds. More specifically, this duality normally relies on spectral cover bundles [7], usually with a non-Abelian structure group, on elliptically fibered CY manifolds. It would be interesting to understand in detail how the duality works for line bundle models. In the present paper we focus on the construction of models on the heterotic side as a first step in this direction while the discussion of heterotic F-theory duality for these models will be the subject of future work.
There are strong indications [8] that the construction of phenomenologically interesting heterotic CY models requires a CY manifold with a non-trivial first fundamental group which facilitates the presence of a Wilson line. In such constructions, the part of the bundle with a non-flat connection is used to break to a GUT group, such as SU (5), while the Wilson line breaks further to the standard model group. The correct chiral asymmetry is easily engineered at the GUT level by imposing a single condition on the index of the bundle. The Wilson line, while breaking the GUT multiplets up into standard model multiplets, does not change this chiral asymmetry. In this way three chiral families are easily obtained by imposing one condition on the bundle. On the other hand, the direct breaking to the standard model group without the presence of a Wilson line requires imposing one index condition per standard model multiplet. It can be shown [8] that the combination of these conditions leads to serious phenomenological problems.
For this reason we would like to follow the standard two-step construction with an intermediate GUT group, which we choose to be SU (5), and subsequent Wilson line break-

JHEP04(2018)087
ing. The required CY manifolds with a non-trivial first fundamental group are usually constructed by starting with a simply-connected "upstairs" CY manifold X with a freelyacting discrete symmetry Γ and then taking the quotient X/Γ. We would like to follow this approach and should, therefore, construct elliptically fibered CY three-folds with freelyacting symmetries. For our purpose, we consider the simplest case where Γ = Z 2 , that is, we consider freely-acting involutions. The construction and analysis of such elliptically fibered CY three-folds with a freely-acting involution constitutes another technical complication, partly addressed in earlier work [9,10], which we review and further develop for our purposes. For the vector bundle to descend under the quotient map, it must be 'equivariant' under the symmetry action. In this work, we restrict ourselves to a necessary check for equivariance, as this is better suited to a systematic model scan.
The main idea for constructing such freely-acting involutions is as follows. Consider an elliptically fibered CY three-fold X with projection π : X → B, two-fold base B and with a section σ : B → X. We can describe the typical elliptic fiber, E b = π −1 (b) where b ∈ B, by its Jacobian, that is, by a complex w plane with identifications w ∼ w + 1 and w ∼ w + τ b , where τ b is the complex structure of E b . An obvious starting point is to consider an involution ι X of X which acts on the fiber as a half-shift, that is, as ι E : w → w + 1/2. This looks promising since this action of ι E is already fixed point free. However, there are two complications which can be inferred from the following argument. The presence of the half-shift on the fibers means that the elliptic fibration has to have a second section, ζ = ι E • σ, in addition to σ. Such elliptic fibrations with two sections are known to have a special structure and, in particular, the discriminant locus (the locus on the base B over which the torus fiber degenerates) splits into two components. It turns out that on one of these components the fiber degenerates such that the entire CY manifold becomes singular. To produce a smooth manifold, the singular loci need to be blown-up and this leads to an additional (second and fourth) cohomology class. The other component of the discriminant locus does not lead to singularities of X but the action of ι E on the corresponding degenerate fibers is not fixed point free. To remedy this problem we construct ι X by combining the half-shift, ι E , on the fiber with a simultaneous action ι B on the base. The latter does not have to be fixed point free on B. In order to ensure that ι X is fixed point free it is sufficient to require that the fixed point locus of ι B should not intersect the above mentioned second component of the discriminant locus where ι E has fixed points on the fiber. This is generically the case if ι B has at most fixed points (rather than fixed curves) on B and this is what we will require. This construction has implications for the possible choices of base spaces B. Smoothness of a generic Weierstrass model over B implies that B should be weak Fano [11,12]. Further, finding a freely-acting involution ι X along the lines described above requires the existence of an involution ι B with at most fixed points. We will focus on two-dimensional toric Fano base spaces B and find that there are six such spaces which lead to smooth CY three-folds and have a suitable involution ι B .
For such elliptically fibered three-folds X with a freely-acting involution and two sections we systematically develop the required model building tools. This includes the construction of a suitable integral basis of the second and fourth homology, calculation of the JHEP04(2018)087 intersection ring, the construction of Kähler and Mori cones and the analysis of line bundles L → X and their properties. In particular, we show how to calculate the cohomology of line bundles by combining the Leray spectral sequence with known methods for calculating line bundle cohomology on toric manifolds [13,14].
Based on these results, we scan rank five line bundle models for the six suitable base spaces and find more than 400 SU(5) GUT models with the correct chiral asymmetry. Upon taking the quotient by the involution and introducing a Z 2 Wilson line in the hypercharge direction these will lead to models with the standard model gauge group and three chiral families. A common feature of all models is the presence of a large number of vectorlike states. In SU(5) language all models have at least one 10-10 vector-like pair and at least 20 (!) 5-5 vector-like pairs. Particularly the latter number is surprising compared with the results obtained for heterotic line bundle models on complete intersection CY manifolds [5,6,15]. For complete intersection CY manifolds, imposing the right chiral asymmetry frequently meant the absence of 10-10 vector-like pairs and a small number of 5-5 vector-like pairs.
The plan of the paper is as follows. In the next section, we will introduce elliptically fibered CY manifolds and discuss their properties including line bundles on these spaces. In section 3, heterotic line bundle models are reviewed and in section 4 we present the results of our model scan. We conclude in section 5. In the main part of the paper, we will keep the discussion informal and as non-technical as the subject allows; technical details can be found in the appendices.

Elliptically fibered CY three-folds and involutions
In this section, we discuss the required background on elliptically fibered CY three-folds, specifically in the presence of a freely-acting involution. Much of the material is not new and can be inferred by combining results in the literature [9,10] but we would like to present a coherent exposition as required for systematic model building. This will be kept somewhat informal with focus on the main results while technical details can be found in appendices A and B.

Elliptically fibered CY three-folds with a single section
As a warm-up, we consider elliptically fibered Calabi-Yau three-folds X over a two-fold base B with projection π : X → B and a single section σ : B → X. Points on the base are denoted by b ∈ B and K B is the canonical bundle of the base. A Weierstrass model for X is given by the equation which, for each base point b ∈ B, describes an elliptic curve, E b = π −1 (b), embedded in P 2 with homogeneous coordinates x, y, z. Here, f and g are sections of the line bundles K −4 B and K −6 B , respectively, which encode the variation of the elliptic curve over the base. In these P 2 coordinates the section can be explicitly written as σ(b) = (b, (0, 1, 0)), that is, it is "located" at the point (x, y, z) = (0, 1, 0) ∈ E b in each fiber.

JHEP04(2018)087
A typical torus fiber E b can also be described by its Jacobian, that is, by a complex w plane with identifications w ∼ w + 1 and w ∼ w + τ b , where τ b is the complex structure of E b . For each point b ∈ B, this complex structure is related to the sections f and g by the equation (2. 2) The denominator on the right-hand side of this expression, that is, and is called the discriminant. The discriminant locus, defined by {b ∈ B | ∆(b) = 0}, is a curve in the base over which the fiber E b degenerates. This will be discussed in more detail in the next sub-section for the case of two sections.
A basis of curves and divisors of X can be obtained from a basis of curves on B by using the maps π −1 and σ. Other relevant properties of X, such as the intersection ring, the Mori cone and the Kähler cone, can also be obtained in terms of properties of the base. Since we are primarily interested in the two-section case, we will defer the details of this to the following sub-sections.

Elliptically fibered CY three-folds with two sections
As indicated earlier, we would like to construct elliptically fibered CY three-folds X with a freely-acting involution ι X , starting with an involution ι E on each fiber which acts as a half-shift w → w + 1/2. The presence of ι E implies the presence of a second section ζ = ι E • σ : B → X of the fibration, located at the "two-torsion point" in each fiber, in addition to the section σ : B → X. We should, therefore, discuss the structure of elliptically fibered CY three-folds with two sections.
Starting from the general Weierstrass model (2.1), a Weierstrass model with a second section can be found by an appropriate tuning of parameters, the general form of which has been found in ref. [16]. As explained in appendix A.2, a second section located at the two-torsion point of the elliptic fiber emerges by choosing the specific, factorisable form (see also ref. [9]) of the Weierstrass model, where α and β are sections of K −2 B and K −4 B , respectively. Comparison with eq. (2.1) shows that this corresponds to choosing the sections f and g in the general Weierstrass model as The two sections of such a fibration are given by so they are located at (x, y, z) = (0, 1, 0) and (x, y, z) = (α, 0, 1), respectively.

JHEP04(2018)087
As can be inferred from the factorised form of the Weierstrass model, this model has singularities which we need to resolve. This can be explicitly seen be working out the discriminant (2.3) which becomes (2.7) Using this notation and (shifted) affine coordinates X = x/z − α, Y = y/z on the patch of P 2 where z = 0, it is shown in appendix A.2 that the half-shift ι E on regular fibers can be written as From the above discriminant, there are singularities, corresponding to fibers of Kodaira type I 2 , over the curve {b ∈ B | ∆ 2 (b) = 0} in the base. In order to obtain a smooth CY three-fold, these singularities need to be resolved, after which the fibers over the locus ∆ 2 (b) = 0 become a pair of P 1 's touching in a point. In addition, the action of ι E has to be extended to these blown-up fibers. Appendix A.2 provides the details of the blow-up procedure and shows explicitly that ι E can indeed be extended to the blow-up and has no fixed points on the blown-up fibers. More specifically, it turns out that ι E exchanges the two P 1 curves of the blown-up fibers. For these reducible fibers, the section σ takes values in one of the P 1 curves and the second section ζ takes values in the other P 1 . From eq. (2.7), there is another component of the discriminant locus, given by the curve {b ∈ B | ∆ 1 (b) = 0} in the base. The degenerate fibers over this locus do not lead to singularities of the entire three-fold so that there is no need for a further blow-up. However, as discussed in appendix A.2 , the map ι E does have fixed points on the fibers over this locus, while it is fixed point free on all other fibers. This means that the involution ι X , if it is to be fixed point free, cannot simply be taken to be ι E with a trivial action on the base. However, we can obtain a fixed point free involution ι X by combining ι E with an involution ι B on the base whose fixed point locus does not intersect the curve {b ∈ B | ∆ 1 (b) = 0}. This is generically the case if ι B has at most fixed points (rather than fixed curves) on B and this is what we will require. The presence of such an involution ι B places constraints on the allowed base manifolds which we will discuss below.
Alternatively, one may realise elliptic Calabi-Yau three-folds with a second section at the two-torsion point by embedding the elliptic fiber into the Hirzebruch surface F 2 , see [16][17][18]. We will review this approach in appendix A.3. This approach has the advantage that, for a toric base space B, the corresponding elliptic Calabi-Yau three-fold can be presented as a hypersurface in a toric variety. We will use this alternative method to realise the relevant CY manifolds to check some of our results, particularly in relation to line bundle cohomology.

Curves and divisors
We would now like to construct a number of objects, as required for our model building purposes, including a basis of curve and divisor classes, on elliptically fibered three-folds X with a freely-acting involution and two sections, in terms of the corresponding objects on the base B. Technical details can be found in appendix A.4.

JHEP04(2018)087
We begin by introducing an integral basis {C i } of the second homology of the base B, as well as a dual basis, {C i }, of curve classes such that Here and in the following we use indices i, j, . . . to indicate the index range 1, . . . , h 1,1 (B). The intersection forms for these basis curves are denoted by and it is easy to see that g ij is the inverse of g ij . Further, these "metrics" can be used to raise and lower indices, that is, C i = g ij C j and C i = g ij C j . For later purposes, it is also useful to introduce the quantities where K B = −c 1 (B) is the canonical bundle of the base and λ i = g ij λ j . We can use either one of the two sections σ and ζ to raise these curve classes on the base to curve classes on the entire three-fold. The three-fold has two more classes which cannot be obtained in this way. These are the class F of a generic fiber and the new class N introduced by the blow-up, the latter chosen such that the component of the reducible fibers over {b ∈ B | ∆ 2 (b) = 0} which is met by the zero section σ has class F − N . For a general curve class, C, on the base we have the relation which shows that lifts of base curves with the two sections σ and ζ are linearly related. It is, therefore, sufficient to consider lifts by one of the sections and we will use the zero section σ for this purpose. Accordingly, we introduce a basis {C I }, where I = (0,0, i), of curve classes on X by Divisor classes on X can be obtained from curve classes, C, on the base by the inverse image π −1 (C). There are two further classes, the images σ(B) and ζ(B) of the base under the two sections, which cannot be obtained in this way. Hence, we introduce a basis {D I } of divisor classes, where I = (0,0, i), by (2.14) From the intersections in table 4 in appendix A. 4 we can see that this basis of divisor classes is dual to the above basis of curve classes, that is with all other components either fixed by symmetry from the above or else vanishing. For our model building purposes we also require the Mori and Kähler cones of X. We begin with the Mori cone, M X of X which, following appendix A.5, can be written as where M B is the Mori cone of the base. If we write Kähler forms as J = t I D I with the Kähler moduli t = (t I ) relative to the divisor basis {D I }, the Kähler cone is the dual of the Mori cone, that is, (2.20) and can, hence, be explicitly worked out once the Mori cone is known. We should also discuss how the above basis of curve and divisor classes relates to the involution ι X on X and its action ι B on the base. We find where C is a curve class on the base. In particular, as is evident from the second line, ι X exchanges D 0 = σ(B) with D0 = ζ(B). Finally, the second Chern class and Euler number of X are given by (see refs. [9,10]) We have thus succeeded in expressing all basic properties of elliptically fibered CY threefolds with two sections in terms of corresponding properties of the base.

Line bundles
We would now like to collect properties of line bundles L → X on elliptically fibered CY three-folds with a freely-acting involution and two sections. As usual, we denote by O X (D) a line bundle with first Chern class or character where c 2I (X) are the components of the second Chern class of X, defined by c 2 (X) = c 2I (X)C I , given in eq. (2.24). We also require an expression for the slope of a line bundle.
With the Kähler form written as J = t I D I the slope of a line bundle L = O X (k I D I ) is defined by Later, we will be interested in line bundles L whose slope vanishes somewhere in the (interior,K X , of the) Kähler cone (2.19), so we have to solve the quadratic equation To build heterotic line bundle models, we also need to know which line bundles L admit an equivariant structure under the involution ι X . In fact, in this paper, we will be content checking invariance of L, a necessary but not sufficient condition for equivariance which is better suited for a systematic model search.
From eqs. (2.22), the last condition can be worked out more explicitly as where I B is a matrix which describes the action of the involution ι B on the basis {C i } of curve classes on the base such that Stated differently, invariant line bundles must be of the form where n ∈ Z and L = O B (k i C i ) is a line bundle on the base which satisfies I i Bj k j = k i . To determine the full spectrum of line bundle models we need to compute line bundle cohomology, rather than merely line bundle indices. Later on, we will be interested in line bundles L → X whose slope µ X (L) vanishes somewhere in the interior of Kähler cone. It is useful to note that, from a general vanishing theorem [19], such line bundles satisfy with the exception of the trivial line bundle. Even for such line bundles, the index does not provide the full information and at least one more cohomology needs to be computed. The cohomology of line bundles L → X on an elliptically fibered CY three-fold X can be expressed in terms of line bundle cohomology on the base B, using the direct image π * L and the associated higher direct images, R q π * L, together with the Leray spectral JHEP04(2018)087 sequence [7,20,21]. Here we outline the structure of this calculation and its main results. Further details can be found in appendix C, including details on the Leray spectral sequence and computations of the direct images and higher direct images.
The starting point of the computation is the exact sequence where which follows from the Leray spectral sequence. For our model building effort we require invariant line bundles, that is line bundles of the form (2.32), and we will restrict the cohomology calculation to such cases. To work out the higher direct images of such invariant line bundles the formula R i π * (O X (nΣ) ⊗ π * (L)) = R i π * (O X (nΣ)) ⊗ L is of some help and shows that all we require is the (higher) direct images of the line bundles O X (nΣ). These are explicitly given by (2.37) The above results imply that E 2,0 2 = 0 for n < 0 and E 0,1 2 = 0 for n > 0 and, hence, that the sequence (2.34) splits for all n = 0 . For n = 0 this cannot be inferred in general. However, it turns out that, for our choices of base spaces and line bundles with n = 0, one of E 2,0 2 and E 0,1 2 is always zero. This means that the sequence (2.34) splits for all cases of interest. In conclusion, for line bundles invariant under the involution ι X , that is line bundles of the form L = O X (nΣ) ⊗ π * (L) (with L a line bundle on the base), we have Hence, we can compute the relevant line bundle cohomologies on the CY three-fold X from line bundle cohomologies on the base B.

Base space choices and involutions
We have now expressed all relevant properties of elliptically fibered CY three-folds X with two sections in terms of properties of the base manifold B. The final step in our model building set-up is to find a suitable explicit class of base manifolds and to develop all their required characteristics. Here we will present the main results with details relegated to appendix B.
We begin by discussing the general constraints imposed on the base manifolds. Recall from section 2.2 that we require the existence of a freely-acting involution ι X of X which is constructed by combining a half-shift ι E on the elliptic fibers with an involution ι B on the base. The presence of this involution implies the existence of two sections, σ and ζ, of the fibration, as well as the factorisation of the discriminant, that is, ∆ = ∆ 1 ∆ 2 2 . This leads to two requirements on the base manifold B. Firstly, the three-folds X → B described by eq. (2.4) should be smooth after resolving the singularity over the locus ∆ 2 = 0. Secondly, in order for ι X to be fixed point free, the fixed point locus of ι B should not intersect the locus ∆ 1 = 0, where ι E is not freely-acting on the fibers. To ensure this we require that ι B has at most fixed points on B.
The first requirement means in particular that a generic Weierstrass model (2.1) over B should be smooth, which happens only when the base space B is a weak Fano two-fold [11,12]. If the base space B is not weak Fano, it is still possible to construct a corresponding smooth elliptic CY three-fold X by resolving the singular Weierstrass model. After such a resolution, X will have additional divisors beyond the two sections and divisors inherited from the base, which need to be taken into account when constructing and analysing models. The properties of such additional divisors depend on the choice of the base space which complicates the model building significantly. For this reason, we focus on the case of weak Fano base spaces in the present paper. More specifically, we will only consider toric weak Fano two-folds. From the 61,359 toric base spaces giving rise to flat elliptic CY three-folds [11], 16 are weak Fano and lead to generically smooth fibrations [11,12]. These correspond to the 16 inequivalent reflexive polygons in two dimensions. (A list of these 16 reflexive polytopes can be found in figure 1 of ref. [22].) To implement the second requirement -the existence of an involution ι B with at most fixed points -we have searched these 16 surfaces for involutions which can be realised linearly on the homogeneous coordinates. We find that 6 out of these 16 manifolds allow for an involution of this kind which has at most fixed points. The polygons for these 6 possible base manifolds B, on which we will focus for the construction of line bundle models, are shown in figure 1. Some basic properties of these six base spaces are listed in table 1 with more details available in appendix B.
From eq. (2.17), the intersection numbers of the CY three-fold X can be expressed in terms of intersections on the base, that is, in terms of the quantities λ i and λ given in table 1 and in terms of the intersection forms g ij = C i · C j given in appendix B. The second Chern class of the CY three-fold X, which will be important to check that models are anomaly-free, can be computed from the data in table 1 inserted into eq. (2.23). This JHEP04(2018)087 Figure 1. The 6 reflexive polytopes corresponding to toric base spaces B which lead to smooth elliptically fibered CY three-folds and allow for an involution with at most fixed points. The homogeneous coordinates x i associated to the rays are indicated and will be used throughout the paper. Table 1. Some basic properties of the six toric bases spaces B which lead to smooth elliptically fibered CY three-folds and allow for an involution with at most fixed points. The coordinates x i have been defined in figure 1 and are also used to denote the curve classes defined by x i = 0. The row denoted h 1,1 inv (B) gives the dimension of the second homology invariant under the involution ι B . The first Chern class of B is obtained from c 1 (B) = −K B and the dual basis, {C i }, of curve classes can be obtained from the intersection forms g ij given in appendix B, via the relation C i = g ij C j . Also recall the definitions λ i = K B · C i and λ = λ i λ i .

JHEP04(2018)087
leads to The Mori cone of X can be obtained from the Mori cone of the base B via eq. (2.18) and the generators of M B for the six base spaces are given in appendix B.
We will also need to know the explicit action of the involution ι B on the homogeneous coordinates x i . For each of the base space choices B, there are multiple distinct actions on the homogeneous coordinates that give involutions with only fixed points. Most of these differ only by signs, and many of these result in the same action, I B , on the curve classes in B (see eq. (2.31)). For the purpose of this paper, it is only the action I B on the curve classes which enters the discussion, while the action on the coordinates is not explicitly used. For this reason, it is sufficient to consider one representative action ι X for each I B . These are explicitly given by Note that there are two inequivalent choices, referred to as cases (a) and (b), for the base space F 15 . These two cases do indeed lead to different actions I B on the curve classes of B and should, hence, both be taken into account. The explicit matrices I B in each case are provided in appendix B.
We recall from eq. (2.30) that line bundles on the CY three-fold which are invariant under the involution ι X involve curve classes on B which are ι B invariant. It is, therefore, important to have explicit expressions for all such ι B invariant curve classes k i C i , that is, classes which satisfy I i Bj k j = k i . Using the explicit matrices I B given in appendix B these JHEP04(2018)087 invariant curves can be characterised as follows.
Finally, we need to be able to compute the cohomology for line bundles on B. Methods to compute line bundle cohomology on toric spaces have been developed in refs. [13,14] and we will use the accompanying code cohomCalg for our calculations.

Line bundle models
All mathematical ingredients for the construction of heterotic line bundle models on elliptically fibered CY three-folds with two sections are now in place. In this section, we review the general construction of line bundle models and the structure of their low-energy spectrum (see refs. [6] for a comprehensive account) as well as some particular features of line bundle models for elliptically fibered CY three-folds with two sections. A systematic line bundle model search on these CY manifolds will be presented in the next section.

Construction of line bundle models
A heterotic E 8 × E 8 line bundle model is defined by two ingredients: a CY three-fold X and a line bundle sum with rank r = rk(V ), where L a are line bundles on X. We think of V as the bundle in the "observable" E 8 sector and will comment on the fate of the hidden sector below. In practice, it is useful to represent line bundles relative to an integral basis {D I } of divisor classes on X and we write L = O X (k) for a line bundle L with first Chern class c 1 (L) = k I D I . Using this notation, the line bundles can each be represented by an integer vector k a and the entire line bundle sum V by a h 1,1 (X) × r integer matrix K = (k I a ). The data (X, V ) is subject to three conditions which are required for the consistency of the model. Firstly, we need to be able to embed the structure group of the bundle V JHEP04(2018)087 into E 8 . Apart from the obvious rank constraint , r ≤ 8, this requires the vanishing for the first Chern class of V , that is, Provided this is satisfied the structure group of V is generically S(U(1) r ) which allows for an embedding into E 8 via the sub-group chain The second constraint originates from the requirement that the bundle V preserve supersymmetry. It can be formulated in terms of the slope of a line bundle where d IJK = D I · D J · D K are the triple intersection numbers of X. Supersymmetry of V requires that the slope of all line bundles L a vanishes simultaneously somewhere in the interior,K X , of the Kähler cone of X. This means that the equations for a = 1, . . . , r should have a common solution t ∈K X . The third requirement is anomaly cancellation which demands the existence of a hidden bundleṼ (which should also be supersymmetric) and a holomorphic curve W ⊂ X (around which a five-brane wraps) such that (3.6) A practical way to ensure that this condition can be satisfied is to demand that where M X is the Mori cone (that is, the cone of effective curve classes) on X. Indeed, provided eq. (3.7) holds, we can always choose a holomorphic curve W with [W ] = c 2 (T X)− c 2 (V ), so that the anomaly condition (3.6) is satisfied for a trivial hidden bundle (although different choices for [W ] combined with a non-trivial hidden bundleṼ are usually possible as well). To check the condition (3.7) explicitly, we introduce a basis {C I } of curve classes on X, dual to our basis {D I } of divisor classes such that C I · D J = δ I J . Then, the second Chern class of V can be written as and can be compared with the second Chern class of the tangent bundle expressed in the same basis as c 2 (T X) = c 2I (T X)C I . Provided the above three conditions are satisfied we have a consistent line bundle model with defining data (X, V ). For phenomenological purposes, we would like to quotient this model by a freely-acting symmetry, Γ, of the CY three-fold X and obtain a model on JHEP04(2018)087 the quotient CY manifoldX = X/Γ. The symmetry Γ should lift to the bundle V (in mathematical terminology, the bundle V should have a Γ-equivariant structure) so that it descends to a bundleV →X on the quotient. The final step is to add a Wilson line bundle W so that the complete "downstairs" bundle isV ⊕ W.
For the purpose of this paper, we would like to carry out a necessary (but not quite sufficient) check, adapted to our algorithmic model search, for equivariance. We will check that the line bundle sum V in eq.
amounts to a permutations of the various line bundles. As a further simplification, we focus on cases where these permutations are trivial, so that every line bundle L a is invariant by itself. Given that the symmetry group we consider is Γ = Z 2 , generated by an involution ι X , we therefore demand that for all a = 1, . . . , r.

Spectrum of line bundle models
The (observable) low-energy gauge group is the commutant of S(U(1) r ), the structure group of the line bundle sum V , within E 8 . From a phenomenological point of view, the most attractive choice is r = 5 and in this case the low-energy gauge group is given by Embedding G ⊂ SU(5) × SU(5) the representations on the right-hand side branch further into the G-representations (3.12) Here, r q denotes an SU(5) representation r with S(U(1) 5 ) charge q, a five-dimensional integer vector, defined up to multiples of (1, 1, 1, 1, 1). Further, e a , where a = 1, . . . , 5, denote the five-dimensional standard unit vector. This means that, for example, the multiplet 10 ea carries charge +1 under the a th U(1) symmetry in S(U(1) 5 ) and is uncharged under the others. For simplicity, we will frequently replace the sub-scipt e a by a, so that, for example, 10 ea is written as 10 a .
The G-representations on the right-hand sides of eq. (3.12) provide a list of possible multiplets which can arise in the effective theory. If we formally assign the charge e a to the line bundle L a then every tensor product of these line bundles and their duals acquires an induced charge, in the obvious way. For example, the line bundles L a ⊗ L b then carry JHEP04(2018)087 Table 2. A list of SU(5) × S(U(1) 5 ) multiplets in the low-energy theory and their associated line bundles. The multiplicity of each multiplet is computed by the cohomology h 1 (X, L) of the associated line bundle L.
charges e a + e b . Then, every multiplet r q in (3.12) can be associated to the line bundle, The first cohomology, h 1 (X, L), of this associated line bundle counts the multiplicity of a multiplet. The details of this correspondence are provided in table 2. We recall that all line bundles L a (and, hence, their duals and all their tensor products) need to have a vanishing slope which implies that h 0 (X, L) = h 3 (X, L) = 0 for all line bundles L which appear in table 2 (with the exception of L a ⊗ L * a = O X ). This means that the chiral asymmetry of multiplets is computed by the index, that is, What should be required of this spectrum for a physically promising model? First, we remark that the additional gauge bosons associated to the S(U(1) 5 ) symmetry do not constitute a phenomenological problem. They are either heavy due to the Green-Schwarz effect or can acquire mass due to spontaneous symmetry breaking induced by the SU(5) singlet fields in the last row of table 2. For the correct total chiral asymmetry in the 10-10 sector we should require that This condition guarantees three chiral 10 families after dividing by the symmetry Γ with order |Γ|. Fortunately, for SU(5) bundles V , we have ind(V ) = ind(∧ 2 V ) so that the above condition also guarantees the correct chiral asymmetry in the5-5 sector. For models where each line bundle L a is equivariant by itself we should also ensure that the chiral asymmetry in every S(U(1) 5 ) charge sector has the "correct" sign, that is, that there is no excess of 10 a over 10 a multiplets for any a and no excess of 5 a,b over5 a,b multiplets for any a, b. To avoid such "wrong" chiral asymmetries we impose that for all a, b = 1, . . . , 5 with a < b.

JHEP04(2018)087
Further, for a "clean" spectrum we should demand the absence of 10-10 vector-like pairs and the presence of precisely one5-5 vector-like pair to account for the Higgs doublets. These requirements can be expressed as If all these conditions are satisfied, a model with the precise MSSM spectrum can usually be obtained after dividing by the symmetry Γ. The only exotic states left in the upstairs theory are the Higgs triplets contained in the5-5 vector-like pair. They can normally be projected out for a suitable choice of equivariant structure and Wilson line, while the Higgs doublets can be kept. While the conditions (3.17) give rise to a clean standard model spectrum in this way they are by no means necessary. Additional vector-like states as they arise when the conditions (3.17) are violated can receive a mass from superpotential terms of the form 1 10 10 and 1 55 (or even higher-dimensional operators with multiple singlet insertions), when the SU(5) singlet fields acquire a vacuum expectation value. It is worth noting that these singlet fields, which appear in the last row of table 2, carry a S(U(1) 5 ) charge so that the presence of such terms is constrained by the S(U(1) 5 ) symmetry. Switching on singlet vacuum expectation values (which, in the language of the effective theory, has to be done preserving D-and F-flatness) corresponds to moving in the bundle moduli space and deforming the line bundle sum V to a bundle with non-Abelian structure group. It is perfectly possible and likely to happen in many cases, that unwanted vector-like states can be removed in this way, but checking this is a matter of detailed analysis within each model.

Line bundle models and elliptic fibrations
We would like to discuss a number of specific features which arise for line bundle models on elliptically fibered CY three-folds with a freely-acting involution of the kind considered in this paper.
First, recall that on such CY three-folds, the line bundles L a can be written as divisor classes has been defined in eq. (2.14). Hence, every line bundle can be represented by an integer vector Alternatively, we will also represent the entire line bundle sum V by an h 1,1 (X) × 5 integer matrix

JHEP04(2018)087
It is interesting to consider the total chiral asymmetry of such models which, from eqs. (2.28) and (3.3), is given by Since the pure base intersection numbers d ijk vanish from eq. (2.17) this means not all of the integers k 0 a (and k0 a ) can be zero if we want to obtain a chiral model. In other words, some of the line bundles L a need to have a first Chern class with a non-zero coefficient in the direction of σ(B) (and ζ(B)), leading to non-trivial bundle upon restriction to the fibers, in order to generate a non-vanishing chiral asymmetry. Heterotic F-theory duality in its current formulation requires bundles which are flat on the fibers [7]. Hence, at present, there is no obvious F-theory dual for chiral line bundle models. It would be interesting to study if heterotic F-theory duality can be extended to such cases.
Finally, to calculate line bundle cohomology, we use the results of section 2.4 to express the cohomology of line bundles on X in terms of line bundle cohomology on the toric twofold base B. The latter is then computed using the code for line bundle cohomology on toric spaces developed in refs. [13,14].
The CY manifolds used in this paper can also be realised as hypersurfaces in toric four-folds, as described in appendix A.3 and this provides us with an alternative method to calculate line bundle cohomology. More specifically, a line bundle L on the toric ambient four-fold A and its restriction L to the CY three-fold X are related by the Koszul sequence 22) where N is the anti-canonical bundle of A. The cohomologies of L and N * ⊗L can again be computed using the code in refs. [13,14], this time applied to the ambient toric four-fold A.
The cohomology of L can then be inferred from the long exact sequence associated to the above Koszul sequence. We have carried this out for ι X invariant line bundles and we find that the cohomology of L can be determined easily from the long exact sequence without the need to compute ranks of maps (that is, the sequence always splits). Moreover, the results always agree with the previous method, based on the Leray spectral sequence.

Systematic model search
We have now collected all ingredients required for the construction of heterotic line bundle models on elliptically fibered CY three-folds with a freely-acting involution. In this section, we describe the results of a systematic scan, searching for physically promising models, which covers the six possible base manifolds in table 1 and rank five line bundles sums.
We also illustrate our results by explicitly presenting a specific model found in this scan.

Scan results
Our search has been carried out for all six base spaces in  Table 3. The number of phenomenologically interesting models found for each of the six base manifolds. The scan was carried out over all line bundle models with |k I a | ≤ k max and k mod gives the largest value of |k I a | which arises in a physically interesting model. constituent line bundle is ι X invariant) and whose entries are bounded by |k I a | ≤ k max . The quantity k max has been maximised in view of computational limitations on a desktop machine and its values are listed in table 3.
From the line bundle sums generated in this way, we have selected the physically promising ones by the following set of criteria.
• There is a locus in Kähler moduli space where the slopes of all line bundles L a vanish, that is, eq. (3.5) is satisfied for all L a . This means the line bundle sum V preserves supersymmetry.
• Following eq. (3.15), the index of the line bundle sum satisfies ind(V ) = −6. This guarantees three chiral families of quarks and leptons after taking the quotient by the involution ι X .
• The indices of L a and L a ⊗ L b are constrained by eqs. (3.16) to avoid a chiral asymmetry with the wrong sign in any S(U(1) 5 ) charge sector.
The number of models satisfying these conditions is given in the last column in table 3.
The list of all integer matrices K for those physically promising models can be downloaded from ref.
[23]. The number of models found for each base manifold B can be qualitatively understood by considering the number h 1,1 inv (B)+1 (see table 1 for the values of h 1,1 inv (B)) of independent integers which specify a ι X invariant line bundle. For B = F 2 , F 4 we have h 1,1 inv (B) = 2 and this does evidently not provide enough freedom to allow for interesting models. For B = F 7 , F 9 we have h 1,1 inv (B) = 3 and in these case we are able to find all physically promising models by extending the scan to a sufficiently large k max .   and we have only carried out a partial scan for k max = 3. For those cases, the number of interesting models exceeds the numbers given in table (3). In total, we find 442 models for all six base manifolds.
For those 442 models, we have also determined the complete spectrum by computing all relevant line bundle cohomologies, as explained in section 3.2. The results of these computations are summarised in figures 2, 3 and 4 which provide frequency plots for the number of 10 multiplets, 5 multiplets and singlets, respectively. As these plots show, there is unfortunately no model without additional vector-like pairs. From  exists always at least one 10-10 vector-like pair and frequently many more. Figure 3 shows that the situation is worse for 5-5 vector-like pairs, where the minimal number is 20. This large number of vector-like pairs comes as a surprise, given the experience with line bundle models on complete intersection CY three-folds [5,6,15] where imposing the correct chiral asymmetry frequently resulted in the absence of additional vector-like states.
As we have argued above, these vector-like states do not necessarily render the models unphysical since they can be given a mass via couplings to singlet fields with non-vanishing vacuum expectation values. Figure 4 shows that our models do indeed have a significant number of such singlet fields and it is likely that they can be used to remove unwanted vector-like pairs in many cases. Analysing this is a matter of more detailed model building which is beyond the scope of the present paper.

An example model
As an illustration, we will now present one of the physically interesting models from the previous sub-section in detail. Our example is for the base space B = F 7 = dP 3 which, following appendix B.3, has a basis of curve classes (4.1) Writing L a = O X (k I a D I ) as before, the integer matrix K = (k I a ) which defines our example is given by

JHEP04(2018)087
with every column representing one of the line bundles L a . First, we note that the columns sum up to zero so that the constraint (3.3), c 1 (V ) ∼ a k a = 0, is indeed satisfied. Further, we see that the matrix is consistent with ι X invariance of each line bundle. Specifically, the first and second rows are identical, so that k 0 a = k0 a , and the last four rows satisfy k 1 a = k 2 a + k 3 a + k 4 a , in accordance with eq. (3.19). Following the list of required properties in section 4.1, we should next check that the slope of all line bundles vanishes somewhere in the Kähler cone. Given the structure of the matrix (4.2) it is sufficient to do this for the line bundles L 1 and L 2 whose slopes are explicitly given by It can be verified that µ X (L 1 ) = µ X (L 2 ) = 0 for Comparison with eq. (B.29) shows that this point it indeed in the interior of the Kähler cone of X. Next, we should verify the anomaly condition for this model. The second Chern class of the bundle V is given by and comparing this with the second Chern class of the tangent bundle (2.40) gives In summary, we have seen that this line bundle sum is invariant under the involution and provides a consistent model with the correct chiral asymmetry.
To determine the spectrum in more detail we consider line bundle cohomologies. For the line bundle sum V we find h • (X, V ) = (0, 7, 1, 0) .  8, 8, 0) . (4.14) Combining these results the chiral spectrum of the model is  This illustrates the aforementioned proliferation of vector-like pairs. In addition, across all S(U(1) 5 ) charge sectors, the model has h 1 (V ⊗ V * ) = 94 singlet fields whose S(U(1) 5 ) charges we do not list explicitly.

Conclusions
In this paper, we have studied heterotic line bundle models on elliptically fibered CY threefolds. Standard heterotic model building requires a non-trivial first fundamental group of the CY three-fold which is normally realised by starting with a simply-connected CY threefold X with a freely-acting symmetry Γ and then taking the quotient X/Γ. This has led us to study elliptically fibered CY three-folds with the simplest type of symmetry, a freelyacting involution ι X . We have realised this involution by combining a half-shift ι E on the elliptic fibers with an involution ι B on the base B of the fibration. Such elliptic fibrations necessarily have two sections which are exchanged by ι E . We have systematically developed the tools required for heterotic line bundle model building on such manifolds, including the calculation of line bundle cohomology, by expressing all relevant properties in terms of properties of the base B.
The choice of base spaces is restricted by requiring a generically smooth Weierstrass elliptic CY three-fold X and a fixed point free involution ι X . The latter can be realised if the involution ι B on the base has at most fixed points. In this paper, we have focused on toric two-fold base spaces B and, from this class, the above requirements single out six spaces, represented by the reflexive polygons shown in figure 1.
For those six base spaces, we have systematically searched for physically promising models, that is models based on rank five line bundles sums V , with each constituent line bundle being ι X invariant, and a total chiral asymmetry of six. After taking the quotient by ι X and including a Wilson line these give rise to theories with the standard model group and three chiral families. Over the six possible base manifolds, we have found a total of 442 models of this kind and we have computed the complete spectrum for all these models.

JHEP04(2018)087
A generic feature is the presence of vector-like states, specifically at least one 10-10 pair and at least 20 5-5 pairs in each case. In this respect, the results are very different from the ones obtained for line bundle models on complete intersection CY three-folds [5,6,15], where vector-like pairs were absent for most models with the correct chiral asymmetry. The underlying reason for this difference seems to be the different structure of line bundle cohomology for the two constructions. Complete intersection CY three-folds are defined in ambient spaces consisting of products of projective spaces. They inherit to some degree the relatively simply structure of line bundle cohomology on projective spaces, where at most one cohomology of a line bundle can be non-vanishing. On the other hand, for the elliptically fibered CY three-folds considered in this paper, line bundle cohomology is determined from line bundle cohomology on the toric two-fold base (or, alternatively, from line bundle cohomology on the ambient toric four-fold), which tends to be more complicated as compared to projective spaces.
The presence of vector-like states does not mean that our models are phenomenologically ruled out. All models contain standard model singlet fields (which should be interpreted as bundle moduli) and, depending on details, may allow for superpotential couplings of the form 1 10 10 or 1 55 (or similar, higher-dimensional operators with multiple singlet insertions). Non-trivial singlet vacuum expectation values which correspond to deformations of the bundle away from a line bundle sum and to a bundle with non-Abelian structure group can then generate masses for the vector-like pairs and remove these states from the low-energy spectrum. A detailed study of this is beyond the scope of the present paper and a subject of future work. Specifically, it would be interesting to check if all unwanted vector-like states can be removed in this way while the desired pair of Higgs doublets can be kept light at the same time.
We have also seen that chiral line bundle models on elliptically fibered CY three-folds necessarily have a bundle V which restricts non-trivially to the fibers. Since heterotic Ftheory duality is formulated in terms of spectral cover bundles which, by construction, are flat on the fibers it is not clear whether the models found in this paper have an F-theory dual. It would be interesting to investigate the possible relation to F-theory further and work along these lines is in progress.
It is likely that our model building strategy can be generalised to base spaces which are not weak Fano. Besides finding appropriate involutions on those base spaces, this would require us to find an extension of the half-shift involution ι X to the reducible fibers present in such models. As such models have more divisors than those associated with sections and divisors of the base, the set of possible line bundles and the Mori cone are larger and more complicated. In practice, the realisation of elliptically fibered CY three-folds as hypersurfaces in toric four-folds, as described in appendix A.3, may provide the appropriate framework for such a generalisation to base spaces which are not weak Fano. Pursuing this is an interesting direction for future work.

A Elliptic fibrations and half-shifts
In this appendix we review the required mathematics for elliptically fibered CY three-folds with a freely-acting involution and two sections. Much of the material is not new but can be found in various places in the literature [9,10].

A.1 The group law on elliptic curves and half-shifts
An elliptic curve E can be described by the polynomial equation, where x, y, z are homogeneous coordinates on P 2 and f , g are complex constants. This does indeed describe a smooth curve as long as the discriminant, defined by is non-vanishing. If we denote the nowhere vanishing (1, 0)-form on E by Ω and introduce a basis of one-cycles A, B with A · B = 1 on E the periods can be defined by For a suitable normalisation of Ω we have τ B = 1. The other period, τ = τ A , is called the modular parameter of the torus. We can introduce the lattice Λ ⊂ C generated by 1 and τ and the Jacobian C/Λ of the elliptic curve. Then, the Abel-Jacobi map E → C/Λ is defined by where γ is a path linking p with a given point p 0 ∈ E (conveniently taken as p 0 = [x : y : z] = [0 : 1 : 0]). From now on we work in the patch of P 2 where we can set z = 1. In this patch, the inverse of the Abel-Jacobi map can be expressed in terms of the Weierstrass ℘-function as where w ∈ C/Λ [24]. Via this inverse map, the obvious addition of points in C/Λ can be translated into an addition law on the elliptic curve E. Specifically, if we consider three points w 1 , w 2 , w 3 ∈ JHEP04(2018)087 C/Λ with w 3 = w 1 +w 2 , and introduce the corresponding points (x i , y i ) = (℘(w i ), 1 2 ℘ (w i )), where i = 1, 2, 3, on the elliptic curve, the addition law on E takes the form For fixed (x 2 , y 2 ), this defines a map from (x 1 , y 1 ) to (x 3 , y 3 ) which, provided the discriminant is non-zero, is an automorphism of the elliptic curve.
Since we are interested in constructing involutions of elliptic CY three-folds it is natural to start with involutions on elliptic curves. In the Jacobian description, C/Λ, obvious involutions are half-shifts of the form w → w + ω i , with the three possible choices ω 1 = 1 2 , ω 2 = τ 2 or ω 3 = 1 2 + τ 2 for the shift. Note that these half-shifts are clearly freely-acting and, due to the quotient by the lattice Λ, they do indeed square to the identity. Via the map (A.5), the half-shift points ω i ∈ C/Λ can be translated to certain points (x, y) of the elliptic curve E. It turns out that ℘ (ω i ) = 0 and, hence, that the resulting points are of the form (x, y) = (e i , 0) ∈ E, where e i := ℘(ω i ). Since the y-coordinates of these three points vanish the x-coordinates e i must be zeros of the right-hand side of the Weierstrass equation (A.1).
To be specific, we focus on the first half-shift by ω 1 with corresponding half-shift point (x, y) = (α, 0), where α := e 1 . Since α must be a zero of the right-hand side of eq. (A.1), the Weierstrass equation (for z = 1) factorises as where β is another complex number. Comparison with the standard form (A.1) of the Weierstrass equation shows that f = β − α 2 and g = −αβ. The discriminant (A.2) factorises as well and can be written as In terms of shifted coordinates, defined as X := x − α, Y := y, the above Weierstrass equation takes the form Setting x 2 = α, y 2 = 0, x 1 = x, y 1 = y, x 3 = x and y 3 = y in eqs. (A.6), (A.7) we can translate the half-shift w → w + ω 1 on the Jacobian into a map on the elliptic curve which we denote by ι E . In the affine (x, y) coordinates this map reads explicitly

JHEP04(2018)087
In terms of the shifted X, Y coordinates and their primed counterparts, these transformations can be re-written as Finally, using the Weierstrass equation (A.10) to replace Y 2 in the equation for X , the map ι E can be cast into the simple form Clearly, since the half-shift w → w + ω 1 is freely-acting on the Jacobian, the action of ι E is also free, as long as the elliptic curve is smooth. For later purposes, we note that, in the smooth case, ι E maps the point [ To see this we have to convert the two points into affine coordinates on the patch where z = 1. The point [α : 0 : 1] ∈ P 2 then corresponds to (x, y) = (α, 0) while the point [0 : 1 : 0] ∈ P 2 is mapped to infinity, x, y → ∞, with y 2 = x 3 . By taking limits carefully, it can indeed be shown that (A.14) exchanges these two affine points.
Let us discuss what happens when the elliptic curve becomes singular, that is, when ∆ 1 = 0 or ∆ 2 = 0. If ∆ 2 = 0, the Weierstrass equation (A.8) acquires a double root at x = α (that is, the second factor on the right-hand side of eq. (A.8) develops a root at α). In this case, the transformation (A.14) clearly degenerates badly as all points (X, Y ) are mapped to the origin (X , Y ) = (0, 0). As we will explain, in the context of elliptic fibrations, such fibers will be blown-up and the degeneracy of the transformation will be removed in this way.

A.2 Elliptically fibered CY three-folds with involutions
We would now like to apply the discussion of the previous sub-section to the construction of elliptically fibered CY three-folds, X, with an involution ι X . The idea is to construct ι X by combining the above half-shift ι E on the elliptic fibers with a suitable involution ι B on the base B of the fibration.
We start with an elliptically fibered CY three-fold X with two-fold base B, projection π : X → B and section σ : B → X. Each elliptic fiber E b = π −1 (b) over b ∈ B is described by a Weierstrass equation

JHEP04(2018)087
We would now like to consider a situation where a half-shift ι E acts on every (smooth) elliptic fiber as in eq. (A.14). From the previous sub-section we know that this leads to a factorisation of the Weierstrass equation, namely where α and β should now be seen as sections of K −2 B and K −4 B , respectively. This relates to the standard Weierstrass form (A.16) via As we have seen before, the discriminant factors as As a consequence, the discriminant locus on the base B consists of the two components We should now discuss the singular fibers which arise over the discriminant locus and the fate of the involution ι E over these fibers. We begin with the component of the discriminant locus defined by ∆ 2 = 0. As can be seen from eq. (A.10), the three-fold has an A 1 singularity at X = Y = 0 over this locus which needs to be resolved in order to arrive at a smooth CY three-fold.
A convenient technique to find crepant resolutions of such singularities is to first promote ∆ 2 to a coordinate of the ambient space, which then has to satisfy the equation The blow-up is realised by introducing the new coordinatesX,Ŷ ,∆ 2 and ξ, excising the locus {X =Ŷ =∆ 2 } and modding out by a C * action with weightŝ These new coordinates relate to the old ones as After the proper transform, the Weierstrass equation (A.10) and eq. (A.21), turn intô

JHEP04(2018)087
which describe a smooth CY three-fold. Note that due to the special form of eq. (A.10), the proper transform removes two powers of ξ, so that the blow-up is crepant. The previously singular fibers over ∆ 2 = 0 have now been replaced by two irreducible components,∆ 2 = 0 and ξ = 0, both of which are P 1 s. These two P 1 's are explicitly given bŷ and they evidently touch in the two pointŝ What remains to be discussed is the action of the half-shift (A.14) on the blown-up space. First note that as long as ∆ 2 = 0, we can set ξ = 1 and our equations return to the ones before the blow-up, so that our earlier discussion applies. This means that the action of ι E remains free for fibers which are away from the discriminant locus and that it has fixed points on fibers over ∆ 1 = 0.
In order to discuss the situation for fibers over the locus ∆ 2 = 0, we first need to properly define the half-shift in terms of the new coordinatesX,Ŷ ,∆ 2 and ξ. Of course, we wish to preserve the action of the map on X = ξX and Y = ξŶ , which should still be given by eq. (A.14). Translating these equations into the new coordinates implies We also know that, away from the discriminant locus, the half-shift exchanges the two sections and that, for the reducible fibers over the locus ∆ 2 = 0, each P 1 component is met by one of the sections. Since the two P 1 's are given by∆ 2 = 0 and ξ = 0, respectively, this suggests that ι E should exchange the coordinates∆ 2 and ξ. Hence the appropriate generalisation of the half-shift involution to the blown-up space is This transformation respects eq. (A.28) and it exchanges the two P 1 's, defined by∆ 2 = 0 and ξ = 0, in the reducible fibers. What about the fixed points for this refined version of the involution ι E ? By construction, eq. (A.29) reduces to the previous action for all fibers away from the locus ∆ 2 = 0. Hence, it is fixed point free on fibers away from the discriminant locus and it has fixed points on fibers over the locus ∆ 1 = 0 (certainly as long as ∆ 2 = 0). For fibers over the locus ∆ 2 = 0, the map (A.29) exchanges the two fiber components ξ = 0 and∆ 2 = 0. These two P 1 components touch in the two points (A.27) which are also swapped by the JHEP04(2018)087 map (A.29). Hence, as long as these two points are different the action of ι E on fibers over the locus ∆ 2 = 0 is fixed point free. The two points coincide precisely when ∆ 2 = ∆ 1 = 0 and the resulting fixed point is what we expect on fibers over the locus ∆ 1 = 0. In summary, the involution ι E , as defined by eq. (A.29), is fixed point free over all fibers with ∆ 1 = 0 and it has fix points over fibers with ∆ 1 = 0.
It remains to discuss the two sections σ and ζ for the blown-up version of the CY three-fold. The zero section σ is specified by mapping a point in the base to the point whereX andŶ are taken to infinity such thatŶ 2 =X 3ξ is satisfied along the path. (This is required for the first equation in (A.24) to be satisfied in the limit.) The second section, ζ, maps a point in the base tô whereX andŶ are taken to zero such thatŶ 2 =X∆ 2 is satisfied along the path. (This is required for the first equation in (A.24) to be satisfied in the limit.) It is straightforward to check that the involution (A.29) swaps the two section points (A.30) and (A.31) and, hence, swaps the two sections σ and ζ.

A.3 Alternative realisation of elliptic Calabi-Yau three-folds with half-shifts
There exists another presentation of elliptic Calabi-Yau three-folds with a second section at the two-torsion point which uses a different embedding of the elliptic fiber. Embedding the elliptic fiber in the Hirzebruch surface F 2 (which can also be found as a blow-up of P 2 112 ), an elliptic fibration with two sections has the realisation Here, the x i are the homogeneous coordinates of a toric variety with weight system and Stanley-Reisner ideal generated by {x 1 , x 2 }, {x 3 , x 4 }. Note that we may always find this form by appropriately redefining the coordinates x 1 and x 3 of a general hypersurface in F 2 . In order for a fibration with the above elliptic curve to form a CY space, the b i must be sections of the line bundles −K ⊗i B . The two sections of such a model are located at x 4 = x 1 ± x 2 x 2 3 = 0. Choosing one of the two as the zero section, the other section is located at the two-torsion point if we set b 3 ≡ 0. The involution corresponding to the half-shift is then simply given by and the discriminant becomes the same as (A.19), when using the replacements b 2 = −6iα and b 4 = ∆ 1 = 4β − α 2 . This action is free as long as b 4 = ∆ 1 = 0.

JHEP04(2018)087
In contrast to the realisation via a specialised Weierstrass model, no further blow-ups are required in order to find a smooth Calabi-Yau three-fold X. Furthermore, for a toric base space B, this presentation allows a description of X as a toric hypersurface as follows. Let us denote the rays of the fan of the toric base B by w i ∈ Z 2 . If the convex hull of the vectors    forms a reflexive polytope ∆ • , we can realise (A.32) as a generic toric hypersurface which sits in a fibration of F 2 over B. For the base spaces we consider, which are toric weak Fano two-folds, this condition is always met. In order to find candidate free involutions, we may then restrict to cases where b 3 = 0. In this form, our examples should also be found using the algorithm of ref. [25]. A further upshot of this alternative realisation is that it allows for an easy alternative method to compute line bundle cohomologies. Instead of exploiting the elliptic fibration, as described in appendix C.5, we may simply compute cohomologies on the toric ambient four-fold, using the code cohomCalg [13,14], and then apply the Koszul sequence to find the cohomologies on X. We have used this alternative method to check our results for line bundle cohomology.

A.4 Intersection theory on CY three-folds with involutions
The main purpose of this section is to introduce a suitable integral basis of curve and divisor classes on the CY three-fold X in terms of a basis of curve classes on the base B and to study the resulting intersection properties.
We begin by introducing the required objects on the base B. These consist of a basis, {C i }, where i, j, . . . = 1, . . . , h 1,1 (B), of curve classes, its dual basis {C i } such that and the various intersection numbers where K B is the canonical bundle of B. Note that g ij is the inverse of g ij and that these two metrics can be used to raise and lower indices. In particular, we have λ i = g ij λ j . Next, we should construct a basis of curve classes on the CY three-fold X which can be done by "lifting" the above curve classes C i on the base, using the sections σ or ζ. There are two further curve classes on X which cannot be obtained in this way, namely the class of the generic fiber and the new class introduced by the blow-up. We denote the class of the generic fiber by F and the new class, represented by the ξ = 0 component of the reducible fibers, by N . The other component, defined by∆ 2 = 0, of the reducible fibers then has the class F − N . It turns out that, for a general curve class, C, on the base we have the relation [9] ζ(C) = σ(C) Table 4. Intersection properties of curve classes and divisor classes. Here, C and C are two curve classes on the base. This means that the lifts of a base curve by the two sections are linearly related and that we can focus on one of the sections for the purpose of constructing a basis of curve classes on X. With this in mind we introduce a basis {C I }, where I = (0,0, i) and i = 1, . . . , h 1,1 (B), of curve classes on X by setting The last term in the definition of C i has been included for later convenience.) Then, an arbitrary second homology class in H 2 (X, Z) can then be written as C = n I C I , where n I ∈ Z.
To obtain an integral basis of divisor classes on X, we proceed as follows. We can find divisors of X by lifting curves of the base using the inverse projection map, π −1 , of the fibration. There are only two further divisor classes, the images σ(B) and ζ(B) of the base, which cannot be obtained in this way. Hence, we introduce the basis {D I }, where I = (0,0, i) and i = 1, . . . , h 1,1 (B), by As table 4 shows, this basis is dual to the above basis of curve classes, that is We also require the intersections of two divisors which are summarised in table 5.  Combining the information from table 4 and table 5, we can work out the triple intersection numbers which are explicitly given by with all other components fixed by symmetry or else vanishing.

JHEP04(2018)087
Finally, it can be shown [9,10] that on an elliptically fibered Calabi-Yau X with two sections, the second Chern class c 2 (X) and the Euler number χ(X) can be expressed in terms of properties of the base as

A.5 Mori cone and Kähler cone
For the purpose of checking the heterotic anomaly condition we need to know the Mori cone of the CY three-fold, after the blow-up. As we will see, the Mori cone M X of X can be expressed in terms of the Mori cone M B of the base. First we note that for two effective curve classes C, C ∈ M B in the base, the curve class where n, m ∈ Z ≥0 , must be in the Mori cone of X. We would now like to argue that all effective curve classes on X are, in fact, obtained in this way. It is sufficient to show that all irreducible effective classes can be represented in the form (A.47) (since an extremal ray of the cone of effective classes is an irreducible effective class). We write a general classĈ ∈ H 2 (X, Z) as, with C an element of H 2 (B, Z) andm,n ∈ Z. (Note that, by virtue of the relation (A.38), we do not need to include a term ζ(C ) in this expression.) Let us now require thatĈ is effective and irreducible. Effectiveness means that the projection ofĈ onto the base, which is simply C, is an effective class on the base. Irreducibility ofĈ implies, in particular, that a curve in this class intersects non-negatively any effective divisor in which it is not contained, or more precisely whenever the generic element of the curve class does not sit inside the generic element of the divisor class. IfĈ is not contained in σ(B) or ζ(B), then exploiting this fact leads to the conditions We have seen that the curve C is effective and, for the base spaces we consider, we have −K B · C ≥ 0 for all effective curves. Hence, we conclude thatm ≥ 0 andn ≥ 0. In conclusion, if the effective (and irreducible) classĈ is not contained in σ(B) or ζ(B) then it is of the form (A.47). Now supposeĈ is contained in σ(B). Then it follows thatĈ = σ(C) and, since C is effective on the base, this is already of the form (A.47). A similar argument applies ifĈ is contained in ζ(B). In this case,Ĉ = ζ(C ) with an effective C , which is again of the form (A.47). In summary, the Mori cone of the elliptically fibered CY three fold X consists of classes of the form (A.47) and can, hence, be written as where M B is the Mori cone of the base B.

JHEP04(2018)087
In the context of line bundle model building, we will need to check if the slope of line bundles vanishes somewhere in the Kähler cone, K X , of the CY three-fold X and, hence, we need an explicit description of K X . For the manifolds in question the Kähler cone is the dual of the Mori cone (see for example Theorem 1.4.9 in [26]). Writing a general Kähler form as J = t I D I , where t = (t I ) are the Kähler moduli relative to a basis {D I } of divisor classes, and writing a general curve class as n I C I relative to a dual basis {C I } (which satisfies C I · D J = δ J I ), the Kähler cone can be expressed as and can, hence, be obtained from the Mori cone of X.

B Base spaces
So far, we have expressed all relevant properties of elliptically fibered CY three-folds with two sections in terms of corresponding properties of the base. In this appendix, we will focus on the six suitable toric base spaces which we have identified. More specifically, from the 61, 359 possible toric two-folds classified in ref. [11] only the 16 cases identified in refs. [11,12]  For this basis, we provide the intersection numbers and we note that the intersection matrix g ij := C i · C j is the inverse of g ij . The metric g ij and its inverse g ij can be used to lower and raise indices so that, for example, λ i = g ij λ j . We will also provide a set of generators of the Mori cone M B . We note that for any of these bases B these properties can be read off from the corresponding polytope in figure 1. First note that for any vector v in the plane of the polytope there is a linear relation between divisors,

JHEP04(2018)087
where n is the number of external points on the polygon, and where x i is the vector of the point on the edge of the polytope. Note that there are two independent such relations, and hence h 1,1 (B) = n − 2. To compute intersection numbers, first note that two [x i ] corresponding to adjacent points have mutual intersection number 1. Self-intersections [x i ] · [x i ] can then be computed by combining this fact with the linear relations between divisors. Additionally we can read off the Chern classes, Finally, the Mori cone is generated by the collection of the [x i ].
We have searched for involutions which can be linearly realised on the homegeneous coordinates and which have at most fixed points. We will explicitly specify the action of these involution ι B on the homogeneous coordinates. In addition, we also provide the matrix I B which describes the action of ι B on the curve classes in line with the relation Finally, we indicate the most general form of curve classes C = k i C i which are invariant under ι B as a list of constraints on the integers k i and the resulting dimension, h 1,1 inv (B), of the sub-space of invariant classes.
As for the resulting CY three-folds over these base spaces B, recall from eqs. (2.23), (2.25) and (2.18) that we can express the second Chern class c 2 (X), the Euler characteristic χ(X) and the Mori cone M X in terms of properties of the base B. We also recall that the Kähler cone of the three-fold X is obtained as the dual cone to the Mori cone. Concretely this means one can obtain the inequalities constraining the Kähler moduli by dotting the vector (t I ) of Kähler moduli into each of the Mori cone generators (seeing the latter as vectors relative to the basis C I ). Finally, by noting that h 1,1 (X) = h 1,1 (B) + 2 we can also obtain the Hodge numbers of the three-fold. Using these results, we will list explicitly, for each base space B, the properties of the associated elliptically fibered CY three-fold X with two sections over that base. The intersections of l 1 and l 2 are clearly given by l 1 · l 1 = l 2 · l 2 = 0 and l 1 · l 2 = 1 and, hence, we can introduce a basis of curve classes and its dual by (C i ) = (l 1 , l 2 ) , (C i ) = (l 2 , l 1 ) . (B.7)

JHEP04(2018)087
The relevant intersection numbers are then given by The Mori cone can be written as where the angle brackets are used to denote the list of generators. With the homogenous coordinates (x 1 , x 2 , x 3 , x 4 ) as in figure 1 (where (x 1 , x 3 ) correspond to one P 1 factor and (x 2 , x 4 ) to the other), the action of the involution takes the form This base has the basic properties and the Mori cone is given by For the action of the involution ι B on the homogeneous coordinates we have The resulting transformation on the curve classes is trivial so that I B = 1 2 . Consequently, every curve class on B is ι B invariant and, hence, h 1,1 inv (B) = 2. For the associated CY three-fold this leads to This space equals the del Pezzo surface dP 3 which can be seen as P 2 blown-up in three distinct points. Its basic topological properties are Here, l is the hyperplace class and E a , where a = 1, 2, 3, are the classes of the three blow-ups. In terms of the coordinates in figure 1, these classes can also be expressed as The intersections are given by the well-known formulae so that a suitable choice for the basis of curve classes and its dual is given by In terms of this basis choice the intersection numbers read The Mori cone is generated by The action of ι B on the homogeneous coordinates leads to the action on the curve classes specified by The most general ι B invariant curve class k i C i is then characterised by the constraint so that h 1,1 inv (B) = 3. The associated CY three-fold has the properites The basic properties of this space are given by where x 1 , . . . , x 6 are the homogeneous coordinates as indicated in figure 1. The standard basis of curve classes and its dual can be introduced as This leads to the intersections The Mori cone is generated by  and, in terms of the coordinates in figure 1, a suitable basis of curves is defined by From the intersection numbers we can infer the form of the dual basis in terms of classes [x i ] by using the relation C i = g ij C j . The Mori cone of this space is generated by The involution acts on the homogeneous coordinates as leading to an action on the curves classes specified by The basic topological characteristics of this space are A suitable choice for the basis of curve classes is given by and the dual basis can be inferred from C i = g ij C j using the intersection numbers The Mori cone is generated by For the involution, we find there are two inequivalent choices, which we refer to as cases (a) and (b). Their actions on the homogeneous coordinates are given by which means that h 1,1 inv (B) = 4 in either case. For the associated CY three-fold we find

JHEP04(2018)087
C Line bundles on elliptically fibered CY manifolds In this section we collect useful information on line bundles on elliptically fibered CY manifolds, in particular on how to calculate their cohomology. The material of this appendix can also be found in the literature [7,20,21] and we collect and summarise all the properties relevant to heterotic line bundle model building.

C.1 Basic properties of line bundles
We consider CY three-folds X with a basis {D I } of divisor classes and a dual basis {C I } of curve classes (such that D I · C J = δ J I ) and denote the triple intersection numbers by d IJK = D I · D J · D K . A general Kähler form is written as J = t I D I , where t I are the Kähler moduli relative to the chosen basis of divisor classes.
The line bundle L on X associated to the divisor D = k I D I is denoted by L = O X (D) or sometimes, when it is clear which divisor basis is being referred to, simply by L = O X (k). The first Chern class of this line bundle is given by For the Chern character one finds The Todd class of a CY three-fold is given by Td(X) = 1 + c 2 (X)/12, where c 2 (X) = c 2I (X)C I is its second Chern class. Then, using the index theorem ind(L) = ch(L)Td(X), the index of a line bundle can be written as ind(L) = ch 3 (L) + 1 12 Another quantity we will require is the slope of the line bundle L, defined by These above results for single line bundles can be easily generalised to line bundle sums with rank r = rk(V ). The first Chern class for such a line bundle sum is We will normally be interested in line bundle sums V with c 1 (V ) = 0 and in this case the second Chern class and the index of V are given by (C.9)

C.2 Preliminaries for cohomology calculations
In this appendix, we collect a number of mathematical results which will enter the calculation of line bundle cohomology. We will make use of the Serre duality theorem. Let X be a compact complex ndimensional manifold with canonical bundle K X , and let V → X be a vector bundle with dual bundle V * . Then, Serre duality states that (C.10) In particular, for a CY three-fold, the canonical bundle K X is trivial, so that H q (X, V ) ∼ = H 3−q (X, V * ).
In order to exploit the fibration structure of the CY three-folds we will require (higher) direct images. Let X be an elliptically fibered CY three-fold, with base B and projection map π : X → B. Given a short exact sequence of vector bundles, V i → X, one can ask if such a sequence implies the existence of a short exact sequence for the direct images (or push-forwards), π * V i , which are bundles on the base B. This question finds a general answer in the language of derived functors and higher direct images, as follows. Suppose we have the short exact sequence of bundles The direct image functor π * is left exact, so applying it to the above sequence leads to the exact sequence 0 → π * V 1 → π * V 2 → π * V 3 . (C.12) How can this sequence be continued on the right, so as to keep it exact? The answer is (see, for example, chapter III, theorem 1.1A of ref. [27]) where R q π * are called the 'right derived functors' of the direct image functor, and R q π * V i are called the 'higher direct images' of the bundles V i . It can be shown (see, for example, chapter III, proposition 8.1 of [27]) that the higher direct image R q π * V (q ≥ 1) is the sheaf associated to the pre-sheaf There are two very useful identities for higher direct images that will be helpful below. The first of these identities is the projection formula, R q π * (π * L ⊗ V ) = L ⊗ R q π * V , (C.18)

JHEP04(2018)087
where L is a bundle on the base B and V a bundle on the total space (see, for example, Chapter III.8 in ref. [27]). The second identity, referred to as relative duality, is given by (see, for example, III.12 in ref. [28]) This result is frequently applied in the context of spectral cover bundles, as for example in ref. [7]. Finally, we will make extensive use of the Leray spectral sequence (see, for example, Chapter 5 of ref. [29]). This sequence relates cohomologies of bundles on the total space X to cohomologies of bundles on the base B. For our purposes the following will suffice. For a bundle V , write H p ≡ H p (X, V ) and define E p,q 2 ≡ H p (B, R q π * V ). From the Leray spectral sequence one can show the existence of the exact sequence In our case, E 3,0 2 is a third cohomology on the two-dimensional base B and, hence, vanishes. Similarly, E 0,2 2 vanishes, since it involves R 2 π * V , and this is associated to second cohomologies on the one-dimensional fiber. As a result we have Ker H 2 → E 0,2 2 = H 2 and the above exact sequence becomes If this sequence splits (as it will in all cases we are interested in) it determines H 1 and H 2 in terms of cohomology on the base B. The other two cohomologies can also be expressed in terms of cohomology on the base, via the relations H 0 = E 0,0 2 and H 3 = E 2,1 2 .

C.3 Line bundle cohomology on the base
For most of the base spaces considered in this paper, no analytic expressions for line bundle cohomology are available, and we use the code cohomCalg [13,14] to compute cohomologies on the base. However the cohomology of line bundles on del Pezzo surfaces has been computed, in appendix B of ref. [30]. Since one of the base spaces we consider is a del Pezzo surface, we briefly present these results 1 for reference. Consider a del Pezzo surface dP r with r ≤ 8, given by a P 2 blown-up at r distinct points p a , where a = 1, . . . , r, with hyperplane class l and exceptional divisors E a = ρ −1 (p a ) and the blow-down map ρ. A line bundle on dP r can then be written as where n ∈ Z and b a , c a ∈ Z ≥0 . For n ≥ −2, the cohomology of the line bundle L is given by while the cohomology for n < −2 can be obtained by Serre duality, h q (B, L) = h 2−q (B, L * ⊗ K dPr ), using the expression for the canonical bundle. Here A capa (n) is the dimension of the space of homogeneous polynomials of degree n on P 2 that have order c a zeroes at the points p a .
C.4 Line bundle cohomology for elliptic fibrations with a single section As a preparation, we first consider line bundle cohomology on an elliptically fibered CY three-fold X over base B, with projection π : X → B and a single section σ : B → X.
To make some formulae easier to read, in this subsection and the next we will sometimes abuse notation and write σ for the image of the section σ(B), and similarly ζ for ζ(B), and it should be clear from context which one is meant. Let V → X be a vector bundle over X and, as before, we write H p = H p (X, V ) and E p,q 2 ≡ H p (B, R q π * V ). From appendix C.2 we know that the Leray spectral sequence implies the exact sequence and, in addition, the results H 0 = E 0,0 2 and H 3 = E 2,1 2 . A line bundle L on X can be written as L ≡ O X (nσ) ⊗ π * L , (C. 26) where L is a line bundle on B and n is an integer. Using the projection formula, R q π * (V ⊗ π * L) = R q π * V ⊗ L, we can write the exact sequence explicitly as To work this out further we require the (higher) direct images of the line bundles O X (nσ).
These have been worked out in appendix C of ref. [20] and the result is for n ≤ −2 . (C.29) From these results, we should distinguish the following cases.

C.5 Line bundle cohomology for elliptic fibrations with two sections
As we have seen, elliptically fibered CY three-folds X with two sections, σ : B → X and ζ : B → X, have to be blown-up and, as a result, contain a further curve and divisor class, relative to the single section case. The most general line bundle L → X can now be written as where m, n ∈ Z and L → B is a line bundle on the base. In fact, for our model building purposes we only require line bundles in a sub-class which consists of line bundles of the form L = O X (nΣ) ⊗ π * L , (C.34) where n ∈ Z and Σ = σ(B) + ζ(B).
In order to compute the cohomology of such line bundles, we follow appendix A of ref. [21]. Note that the subsequent derivation only holds if H 1 (B, K −n B ) = 0 for all n ≥ 0, a condition which is satisfied for toric weak Fano base spaces B and, hence, for all base spaces considered in this paper. That this condition holds for toric weak Fano bases can be seen as follows (for more information see for example ref. [31]). The toric Kawamata-Viehweg theorem states that on a compact toric variety B, if D is a nef and big divisor then If B is weak Fano, then K −1 B is of the form O B (D) where D is nef and big, and from this it follows that this is true also of K −n B for n > 1. Hence for weak Fano bases B, using K −n−1 B in the theorem we have that H 1 (B, K −n B ) = 0 for n ≥ 0. We begin by recalling relative duality which states that R 1 π * L = (π * L * ) * ⊗ K B . (C.36) In particular, this implies that R 1 π * O X (nΣ) = 0 for n > 0. Now consider the short exact sequence,

JHEP04(2018)087
Now we consider the associated long exact sequence of higher direct images, but, since R 1 π * O X (σ) = 0, this leads to the short exact sequence for the direct images. Here, we have used that π * O X (σ) = O B and π * (O X (Σ) ⊗ O ζ ) = K B . The space of extensions defined by the sequence (C.39) is isomorphic to H 1 (B, K * B ) and, from our assumption about the base space B, this cohomology vanishes. As a result, the sequence (C.39) splits and we have (C.40) Starting from this result, we will now prove by induction that for all n > 1 . (C.41) The starting point for the proof is the short exact sequence along with the relation O X (nΣ) ⊗ O X | Σ = O X (nΣ)| Σ . Since R 1 π * O X (mΣ) = 0 for all m > 0, we have R 1 π * O X ((n − 1)Σ) = 0 for n > 1, so that the associated long exact sequence of higher direct images truncates to the short exact sequence We note that π * O X (nΣ)| Σ = K n B ⊕K n B . Recall that we are working with a base B for which H 1 (B, K −m B ) = 0 for all m ≥ 0. Then, by the induction hypothesis for π * O X ((n − 1)Σ), we see that this exact sequence splits, and that the result (C.41) follows. An expression for R 1 π * O X (nσ) can then be obtained from relative duality.
In summary, the results for the direct images and higher direct images of the line bundles O X (nΣ) for all n ∈ Z, where Σ ≡ σ(B) + ζ(B), are: for n ≤ −2 . (C.45) Deriving the cohomology for line bundles L of the form (C.34) then proceeds in complete analogy with the single section case discussed in the previous sub-section. All we have to do is replace σ by Σ and, instead of using the results for the (higher) direct images of O X (nσ), we use the above results for O X (nΣ).