Discrete Symmetries of Calabi–Yau Hypersurfaces in Toric Four-Folds

We analyze freely-acting discrete symmetries of Calabi–Yau three-folds defined as hypersurfaces in ambient toric four-folds. An algorithm that allows the systematic classification of such symmetries which are linearly realised on the toric ambient space is devised. This algorithm is applied to all Calabi–Yau manifolds with h1,1(X)≤3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h^{1,1}(X)\leq 3}$$\end{document} obtained by triangulation from the Kreuzer–Skarke list, a list of some 350 manifolds. All previously known freely-acting symmetries on these manifolds are correctly reproduced and we find five manifolds with freely-acting symmetries. These include a single new example, a manifold with a Z2×Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb{Z}_2\times\mathbb{Z}_2}$$\end{document} symmetry where only one of the Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb{Z}_2}$$\end{document} factors was previously known. In addition, a new freely-acting Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb{Z}_2}$$\end{document} symmetry is constructed for a manifold with h1,1(X)=6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h^{1,1}(X)=6}$$\end{document}. While our results show that there are more freely-acting symmetries within the Kreuzer–Skarke set than previously known, it appears that such symmetries are relatively rare.


Introduction
Discrete symmetries of Calabi-Yau (CY) manifolds, both freely and non-freely acting, are important for a number of reasons. As an immediate application, freely-acting symmetries can be divided out, thereby leading to new CY manifolds with smaller Hodge numbers. Such quotient CY manifolds have been used to fill out the previous fairly sparse tip of the Hodge number plot [1] and, due to their relatively small moduli spaces, they are useful for string compactifications. Further, most standard constructions lead to CY manifolds with a trivial fundamental group. Quotients of CY manifolds by freely-acting symmetries, on the other hand, have a non-trivial fundamental group. In the context of string compactifications on CY manifolds, the presence of a non-trivial fundamental group is required whenever gauge field Wilson lines need to be included. For CY compactifications of the heterotic string in particular the inclusion of Wilson lines appears to be the only viable way [2] to arrive at phenomenologically promising models. In fact, the limited knowledge about such freely-acting discrete symmetries is one of the current "bottlenecks" in the attempt to systematically construct heterotic CY vacua [4,5].
Freely (and non-freely) acting discrete symmetries of CY manifolds can be important for string compactifications in yet another way. Provided they are not divided out, such symmetries can translate into discrete symmetries of the resulting lower-dimensional theory. Such discrete symmetries can carry important information about the structure of the effective theory, for example forbid certain unwanted operators, and they have been extensively used in particle phenomenology. The possible origin of such "phenomenological" discrete symmetries within string theory is an important question which provides a further motivation for the present paper.
Let us summarise the systematic knowledge on freely-acting symmetries of CY threefolds to date. The oldest and conceptually simplest set of CY three-folds consists of complete intersection CY three-folds in products of projective spaces (CICYs), some 7890 manifolds which have been classified in Ref. [3]. Considerable progress has been made over the past years in finding freely-acting symmetries of CICYs [6,7], culminating in Ref. [8] which provides a classification of freely-acting symmetries with a linear realisation on the projective ambient space for the entire CICY dataset. It was found that only 195 of the 7890 CICY manifolds have freely-acting symmetries of this kind, although many of these 195 manifolds allow for multiple symmetries. These results have been used to construct new CICY quotients with small Hodge numbers [6,9] and to systematically search for physically viable heterotic models on CICYs [4,5]. Especially the latter work, which has led to the largest set of heterotic standard models to date, would have been impossible to complete without these results on discrete symmetries for CICYs.
Much less is known about discrete symmetries for the largest known set of CY threefolds, hypersurfaces in four-dimensional toric varieties, which were classified by Kreuzer and Skarke in Ref. [10]. Somewhat imprecisely but for ease of terminology we will refer to these manifolds as toric CY (TCY) manifolds. The Kreuzer-Skarke list consists of about half a billion reflexive polytopes leading to an even larger number of associated TCY manifolds obtained by triangulation. In Ref. [11], this process of triangulation has been carried out for all cases with Picard number h 1,1 (X ) ≤ 6 and we will be relying on this dataset. Batyrev and Kreuzer [12] have classified all TCY manifolds X with freely-acting symmetries such that the quotient X/ can again be described as a TCY manifold. Amazingly, among the half a billion reflexive polytopes, there are only 16 freely-acting symmetries of this kind, associated to 16 of these reflexive polytopes.
In general, the quotient manifold associated to a freely-acting symmetry of a TCY manifold does not have to be a TCY manifold itself. For this reason, there is no expectation that these 16 cases provide a complete list of freely-acting symmetries for TCY manifolds. In fact, this is already clear from the overlap between the two sets of CICY and TCY three-folds. For example, it is well-known that the quintic in P 4 , which appears in both data sets, has a freely-acting Z 5 × Z 5 symmetry. The quintic is also among the 16 cases identified by Batyrev and Kreuzer, who find a freely-acting Z 5 symmetry. In other words, only a single Z 5 symmetry is "toric" and appears in the classification of Batyrev and Kreuzer, while the full Z 5 × Z 5 symmetry is not toric and is, hence, not obtained by their method. Further examples of this kind are provided by the bi-cubic in P 2 × P 2 (with a toric Z 3 symmetry among the 16 cases and a non-toric Z 3 × Z 3 symmetry group) and the tetra-quadric in P 1 × P 1 × P 1 × P 1 (with a toric Z 2 symmetry among the 16 cases and various larger, non-toric symmetry groups). These examples certainly show that the 16 spaces identified by Batyrev and Kreuzer can have larger symmetry groups than the toric ones identified in their paper. It is also expected that there are more TCY manifolds with freely-acting (but non-toric) symmetries and the present paper provides an example that shows this is indeed the case.
The main purpose of this paper is to take first steps towards a classification of discrete symmetries for TCY three-folds. While in itself a fairly technical undertaking, this is highly relevant for the above mentioned tasks, including the construction of CY quotient manifolds with small Hodge numbers, finding phenomenologically promising CY vacua of the heterotic string and understanding the string origin of discrete symmetries in particle physics. Given the size and complexity of the Kreuzer-Skarke dataset, a complete classification is also a formidable task well beyond the present scope of available algorithms and computing power. The required triangulations of the reflexive polytopes in the Kreuzer-Skarke list have only been carried out for h 1,1 (X ) ≤ 6 [11] and this sets an upper limit on what we can currently hope to achieve. In this paper we will, in fact, be more modest and, for our systematic search, restrict to TCY manifolds with h 1,1 (X ) ≤ 3. This amounts to a class of some 350 spaces. A further restriction concerns the type of symmetries we take into consideration. For a TCY three-fold X embedded into a toric ambient four-fold A, we first determine the group G A of symmetries of A acting linearly on the homogeneous coordinates of A. Then we classify symmetries groups of X which can be embedded into G A , that is, symmetries which descend from linear actions on the ambient space.
Broadly, the algorithm for finding such symmetries of TCY three-folds X , defined as the zero loci of families of polynomials p in a toric four-fold A, proceeds as follows.
• Find the ambient space symmetry group G A from the toric data for A.
• For a given finite group , study all group monomorphisms → G A .
• For each such group monomorphism, find the family of invariant polynomials p.
• Check if the hypersurfaces defined by such polynomials p are generically smooth.
• For freely-acting symmetries, check if the symmetry action is fixed point free on X .
For the case of freely-acting symmetries, the Euler number, η(X ), of the CY manifolds needs to be divisible by the group order | |. A given TCY manifold X (with η(X ) = 0), therefore, has a finite list of candidate symmetry groups . There are further indices [13], depending on the embedding X ⊂ A, which need to be divisible by the group order | | and which can be used to constrain the possible freely-acting groups further (or, in cases where η(X ) = 0, provide some non-trivial constraints). The net result is a finite list of group orders and, hence, a finite list of possible freely-acting symmetry groups for a given TCY three-fold, each of which can be analysed using the algorithm outlined above.
For non freely-acting groups there is no simple a priori constraint on the group order. However, we can still apply the above algorithm (dropping the last step of checking fixed points) for the finite number of groups at any given group order, starting with group order two and successively increasing up to a desired maximum. In this way, we can find all symmetries (embedded into G A ), freely-acting or not, up to a given maximal group order. In practice, since symmetry groups of CY manifolds tend to be small, this can amount to a complete classification provided the group order can be pushed sufficiently high.
By applying this procedure, we find all freely-acting symmetries, linearly realised on the ambient space, of the ∼ 350 TCY three-folds with h 1,1 (X ) ≤ 3. We recover all symmetries among these manifolds known from the Batyrev-Kreuzer classification, as well as those known from the overlap with CICY manifolds. We find a single new, nontoric symmetry group, Z 2 × Z 2 , on one of the 16 Batyrev-Kreuzer manifolds, where one of the Z 2 sub-groups is the toric symmetry identified by Batyrev and Kreuzer. In total, there are five TCY three-folds with h 1,1 (X ) ≤ 3 and with freely-acting symmetries. Our results are summarised in Sect. 4.
While these systematic results show that there are new freely-acting symmetries, beyond those identified by Batyrev-Kreuzer and those known from CICY manifolds, the search has not led to new manifolds with freely-acting symmetries. To show that such manifolds exist, we also use our algorithm to construct a freely-acting Z 2 symmetry for a TCY manifold with h 1,1 (X ) = 6, which is neither among the 16 manifolds of Batyrev-Kreuzer nor a CICY manifold.
The outline of the paper is as follows. In the next section, we describe the mathematical background for our algorithm at a relatively informal level, with technical details and proofs relegated to the Appendices. Section 2.7 sets out the classification algorithm which is applied to two specific examples in Sect. 3. Our results from the systematic search of all TCY manifolds with h 1,1 (X ) ≤ 3 are given in Sect. 4 and we conclude in Sect. 5.

Construction of Discrete Symmetries
In this section we review some background material on toric geometry and describe our algorithm for classifying freely acting symmetries of TCY manifolds. Our presentation will be informal and we will focus on the main ideas and results relevant for the classification algorithm. Technical details and proofs can be found in the Appendices.

Toric
Calabi-Yau hypersurfaces. The central objects of interest are Calabi-Yau three-folds defined as hypersurfaces in toric four-folds, TCY manifolds for short. This section provides a short review of TCY manifolds, mainly to set the scene and introduce the required notation. For a more thorough discussion the reader may, for example, consult Refs. [11,14,15]. We begin with a broad outline of the construction.
Following Ref. [14], TCY manifolds (along with their mirrors) can be constructed elegantly from a pair of reflexive polytopes and • . Underlying the construction are two four-dimensional lattices M ∼ = Z 4 and N ∼ = Z 4 with pairing · , · , such that ∈ M R and • ∈ N R satisfy A fan can be associated to the polytope • in the following way. First note that reflexivity of • means that the origin of N is the unique interior lattice point of • . All other lattice points v i of • are primitive generators of the rays of the fan . The cones of are given by a triangulation of • , that is by forming subsets of the v i , each of which contains the generators of a cone. We will only consider triangulations which are star, regular and fine. 1 The fan gives rise to a toric four-fold A = P by a construction which we will review below. Within this toric four-fold, TCY three-folds are defined as zero loci of polynomials p, whose constituent monomials are encoded by , the Newton polytope of the hypersurface equation, as explained in more detail below. The so-defined TCY hypersurfaces are denoted by X . Although the ambient space P is not necessarily smooth, for CY three-folds every fine triangulation leads to a smooth hypersurface X [14].
It can be computationally expensive procedure to find all triangulations of a polytope • which give rise to a smooth Calabi-Yau hypersurface. However, not all of the triangulation data is crucial for the geometry of X . In particular, we may ignore everything which happens inside of faces of codimension one of • . Furthermore, it is sufficient to specify a triangulation of • in terms of its cones of maximal dimension, which are denoted by tr( • ). The task of constructing all triangulations for sufficiently simple polytopes, that is, those with few lattice points, has been completed in Ref. [11] and this work represents the natural input data for our algorithm.
Let us now be more specific about the construction. The data in Ref. [11] provides a list v i = (v ir ) r =1,...,4 ∈ • ∩ N , where i = 1, . . . , n, of the lattice points of the polytope • , along with a triangulation tr( • ), that is a list of cones which constitute the fan ⊂ N R .
There are several ways to construct the toric ambient space A = P from the data encoded in the fan . Here we will focus on the "global construction" in terms of homogeneous coordinates which is best suited for our discussion of symmetries. In this approach, P is found as a quotient of an appropriate "upstairs" space by a group G. The starting point for this construction is the space C n with coordinates x i which are in one-to-one correspondence with the ray generators v i of . From this space we need to remove the exceptional set Z ( ) which is obtained as follows. Consider index sets I = {i 1 , . . . , i p } ⊂ {1, . . . , n} for which the corresponding generators {v i | i ∈ I } do not share a common cone in . For any such I , this property is preserved if we enlarge I by adding further indices, so that there is a minimal choice of such index sets generating all of them. This choice is related to the generators of the Stanley-Reisner-ideal and we denote the set of such minimal index sets by I. With the corresponding zero loci we can construct the "upstairs" space of the quotient we want to describe by removing the total zero locus Z ( ) from C n , that is, by To obtain the toric variety A, this upstairs space needs to be divided by the toric group G ⊂ (C * ) n whose elements t = (t 1 , . . . , t n ) ∈ (C * ) n act on C n multiplicatively as and, in general, has a continuous part isomorphic to (C * ) n−4 and a finite part. Note that G only depends on the lattice vectors v i but not on how they are grouped into cones. Consequently, G is the same for any triangulation of a given polytope. The continuous part of G can be explicitly obtained from the linear relations between the generators v i . The coefficients q ri in these relations form a (n − 4) × n matrix Q, the so-called charge matrix, whose columns we denote by q i . The continuous part of the toric group can then be written as where s = (s 1 , . . . , s n−4 ) are the group parameters and we have used the short-hand notation s q i = n−4 r =1 s q ri r . The toric variety A = P associated to the fan can now be constructed as the quotient With the ambient space in hand, the family of polynomials defining the TCY manifolds X ⊂ A is finally given by for a choice of complex constants c m . The sum in this expression runs over all lattice point of the Newton polytope which can be obtained from • using Eq. (2.1). The Hodge numbers h 1,1 (X ) and h 2,1 (X ) of X can also be computed from the data encoded by the polytopes and • using the well-known results from Ref. [14]. The following is an example of a global TCY construction, following the route outlined above. This example will also be useful in order to illustrate our procedure for classifying symmetries and we will repeatedly come back to it as we go along.
Example. Let us consider a polytope • with five vertices v 1 , . . . , v 5 given by the columns of the matrix Furthermore, this polytope also contains the integral points v 6 = (−1, 2, −1, 0) and v 7 = (−1, 0, 1, 0), which are interior to one-dimensional faces, the origin (0, 0, 0, 0), which is the unique interior point, as well as the integral point v 8 = (−1, 1, 0, 0), which is interior to a three-dimensional face and can, hence, be neglected. The triangulation is unique and reads explicitly It follows that the zero set associated to this triangulation is described by the index set and is explicitly given by From the above vertices and Eq. (2.4) the charge matrix (with the columns ordered as v 1 , v 2 , . . . , v 7 ) can be determined as The polar dual polytope can be obtained from Eq. (2.1) and it has the five vertices given by the columns of the matrix ⎛ ⎜ ⎝ The defining polynomial, p, of the TCY manifold can then be computed from Eq. (2.7) by summing over all the lattice points of . For example, the five monomials in p which arise from the above vertices of are explicitly given by (2.14) However, the most general polynomial p has many more terms from the other lattice points of . From the standard formulae, the Hodge numbers of the so-defined TCY manifold X = {p = 0} ⊂ A = P are given by h 1,1 (X ) = 3 and h 2,1 (X ) = 83.

Constraints on the order of symmetry groups.
Freely-acting symmetries of a CY manifold X are severely constrained since certain topological indices need to be divisible by the group order | |. The simplest example for such an index is the Euler number, η(X ), of the manifold. The Euler numbers of X and its quotient X/ by a freely-acting symmetry are related by η(X/ ) = η(X )/| | and this shows that η(X ) must indeed be divisible by the group order | |. For manifolds with η(X ) = 0 this already implies a finite number of possible group orders, and, hence, a finite number of possible freely-acting groups. When the manifold X is embedded into an ambient space A with normal bundle N , as is the case for TCY manifolds, there are further indices which need to be divisible by | |. From Ref. [13], one class of such indices is given by the indices of the bundles N k ⊗ T X l , where k, l ∈ Z ≥0 . Explicitly, they are given by where td(T X) = 1 + 1 12 c 2 (T X) is the Todd class of the CY manifold X . These quantities need to be divisible by | | for all (k, l) but, as has been shown in Ref. [13], this is the case if and only if | | is a divisor of the four indices χ k,l with (k, l) ∈ Another relevant set of indices [13] are the modified signatures of the bundles N k ⊗ T X l , defined by (2.16) where L(T X) = 1 − 1 3 c 2 (T X) is the L-genus of the CY manifold and ch is the Chern character with the curvature R replaced by 2R. It has been argued in Ref. [13] that the only new condition from these signatures arises from σ 1,1 .
In summary, in order to strengthen the constraint on | | from the Euler number (or, in the case η(X ) = 0, to obtain constraints in the first place), we should consider divisibility by | | of the bundle indices χ 0,1 , χ 1,0 , χ 2,0 and χ 3,0 in Eq. (2.15) and divisibility of the signature σ 1,1 in Eq. (2.16). This leads to a finite list of possible group orders | | and, therefore, to a finite list of possible freely-acting groups which can be systematically analysed based on the algorithm described below. For our example this works as follows.
Example. Let us work out the above indices for our standard example. For the Euler number it follows that η(X ) = 2(h 1,1 (X ) − h 2,1 (X )) = −160. Using the explicit formulae in Ref. [13] we find for the other indices The greatest common divisor of these indices is 8 and, hence, the possible orders of freely-acting groups for the manifold are | | ∈ {2, 4, 8}.
As mentioned before, we are not aware of similar constraints on the order of nonfreely acting symmetries. In this case, the best we can do is to apply our algorithm to all finite groups up to a given maximal order.

Symmetry group of the ambient space.
With a description of the ambient toric space in hand, we are now in a position to find the subgroup, G A , of the ambient space automorphism group which is linearly realized on the homogeneous coordinates x i . For simplicity, we will simply refer to G A as the toric symmetry group. More details on the full automorphism group of a toric variety and the toric symmetry group, including proofs for the statements in this section, can be found in Appendix A.
As mentioned above, we will only consider linear transformations on the homogeneous coordinates x i , that is, elements of Gl(n, C). Our first task is to find the sub-group G B ⊂ Gl(n, C) of linear automorphisms of the "upstairs" space B, defined in Eq. (2.3). Then, we will study how such actions descend to the quotient A (2.6), which we will refer to as the "downstairs" space.
We recall that the structure of B is encoded in the set I which contains index sets I ⊂ {1, . . . , n}, one for each zero locus Z (I ) which is removed from C n in order to obtain B. While the index sets I ∈ I can overlap, that is corresponding zero sets Z (I ) can have a non-trivial intersection, it is straightforward to define a refinement of I, denoted J , whose index sets form a partition of {1, . . . , n}. Specifically, we can define an equivalence relation ∼ on {1, . . . , n} by i ∼ j : ⇐⇒ i, j are contained in the same sets I ∈ I. (2.18) The refinement J then consists of the equivalence classes under this relation.
Example. The index set I for our example (2.10) is a clean partition of {1, . . . , 7} and for this reason we have The equality of the two sets is a special feature of this particular example. In general, the index sets I and J can be different.
There are a number of obvious sub-groups of G B which can be written down immediately. These include the discrete permutation groups S = {σ ∈ S n | σ (I ) = I, ∀ I ∈ I} ⊂ S n , R = {σ ∈ S n | σ (I ) ∈ I, ∀I ∈ I} ⊂ S n , (2.20) which, via the obvious embedding S n → Gl(n, C), can also be thought off as sub-groups of Gl(n, C). The first of these groups is the stabiliser group of the zero sets labelled by I ∈ I, that is, it leaves each component, Z (I ), of the zero set unchanged. The second group, R, maps the zero sets into each other. Clearly, both S and R leave the total zero set Z ( ) invariant and are, hence, sub-groups of G B .
It can be shown (see Appendix A) that where S(J ) is the group of permutations on the set J , and P = R/S. The group P can be worked out explicitly by finding all permutation in S n which map the index sets J ∈ J into each other and which preserve the "natural" order of indices within each set J . In other words, P should be thought of as the group which permutes zero sets. Another obvious sub-group of G B is given by where Gl(J, C) denotes the general linear group acting on the coordinates x j , where j ∈ J . It turns out that the full upstairs symmetry group G B can be expressed in terms of the above groups and is, in fact, given by In summary, the upstairs symmetry group G B consists of two parts, the group H B which consists of blocks of general linear groups whose action leaves the components of the zero sets Z (I ) unchanged and a permutation part, P, which permutes components of the zero set.
Example. With the set J for our example in Eqs. (2.19) and (2.21) it follows that the permutation groups S and P are given by where the permutation which generates P (given above in cycle notation) swaps the two blocks In accordance with this structure, let us order the coordinates as ( (1,4), (2,5), (3,6,7)). The elements h ∈ H B can be represented as 7×7 matrices which, relative to this ordering of coordinates, have the form Here h 2 ,h 2 are general linear 2 × 2 matrices, acting on the coordinates (1, 4) and (2, 5), respectively, while h 3 is a general linear 3 × 3 matrix acting on the coordinates (3, 6, 7). For the same ordering of the coordinates, the matrix version of the generator of P in Eq. (2.26) can be written as The upstairs symmetry group G B = P H B can then be viewed as all matrices generated by the matrices in Eqs. (2.28) and (2.29).
The symmetry group G A of the toric variety A is the part of the upstairs symmetry groups G B which descends to the quotient A = B/G. On general grounds, this part is given by the formula (2.30) where N G B (G) denotes the normalizer of G within G B . The task is to evaluate this formula for the upstairs symmetry group G B as given in Eq. (2.24) and this is explicitly carried out in Appendix A. Here, we summarise the result which can be expressed in terms of two further index sets K and L. The index set K is defined as containing the sets K ⊂ {1, . . . , n} which label coordinates with the same charge q i . This means that the toric group G can also be written as where G q K (K ) denotes the group which acts by multiplying the coordinates x k for k ∈ K by s q K . The second index set, L, is simply the refinement of the index sets J and K, that is, it consists of all non-empty intersections J ∩ K , where J ∈ J and K ∈ K. In general, this will be a new index set, distinct from the ones introduced previously. However, for reflexive polytopes, which is our case of interest, it can been shown (see Appendix A) that L = K. The downstairs versions of the groups R and S are defined as where C Gl(n,C n ) (G) is the centralizer of G in Gl(n, C n ). As before, we can write R A as the semi-direct product In other words, the group R A consists of all the elements of r ∈ R which normalize G, that is, all r ∈ R for which the equations have solutionss ∈ (C * ) n−4 for all s ∈ (C * ) n−4 . Then P A contains permutations in R A which preserve the order of indices within each block K ∈ K but permute the blocks with one another. Hence, P A can also be viewed as a sub-group of S(K), the permutation group of the blocks K ∈ K. With this notation the symmetry group of the toric space A can be written as As usual, the above semi-direct product is defined by the multiplication rule for p,p ∈ P A and h,h ∈ H A /G. Evidently, this symmetry group has the same structure as its upstairs counterpart. There is a continuous part, H A /G, with blocks of general linear groups acting on the homogeneous coordinates in a way that is consistent with the structure of A, and a discrete part, P A , which permutes these blocks. For our example, following the above steps leads to the following.
Example. The charge matrix (2.12) for our example shows that the index set K is given by This is indeed merely a refinement of the index set J in Eq. (2.19) so that the intersection of J and K does indeed not lead to a new index set, in line with our earlier claim. If we write the charge vectors (2.12) as (2.38) then, from Eq. (2.31), the toric group can be written as In practice, adopting the coordinate ordering ((1, 4), (2, 5), (3), (6), (7)) as before, the elements of G are given by the 7 × 7 matrices the three parameters of the toric group. From Eq. (2.35), we also find where h 2 ,h 2 are general linear 2 × 2 matrices and h 1 ,h 1 , h 1 ∈ C * . The permutation part P A can be determined by solving the Eqs. (2.34) for this case, using the charge matrix (2.12). The result is that P A is generated by the matrix which simultaneously permutes the blocks (1, 4) with (2,5) and (6) with (7). Hence P A ∼ = S 2 .

Construction of representations.
Having determined the symmetry group G A of the toric ambient space, our next step is to study group monomorphisms R : → G A from a finite group . Given the structure of G A in Eq. (2.35), for γ ∈ , we can always write The multiplication rule (2.36) in G A means that the maps π and r must satisfy and, hence, π : → P A is a group homomorphism (indeed a permutation representation given that P A can be viewed as a sub-group of the permutation group S(K) and r : → H A /G is a π -homomorphism. Therefore, we proceed by first generating the permutation representations π and then, for each π , the corresponding π -representations r . Permutation representations can be obtained by a standard method available in the Mathematics literature, see for example Ref. [16], particularly Section 2, for details. Broadly, this method works as follows. For a permutation representation π : → S(K), where S(K) is the permutation group of the blocks K ∈ K, the set K splits into subsets K i on each of which π acts transitively. For an x ∈ K i we denote by H i ⊂ the stabiliser sub-group of x. Then the coset /H i and K i can be identified via the map H i γ → γ x for γ ∈ . The action of the permutation representation π on K i can now be described, via this identification, by a right-multiplication on the cosets in /H i . Hence, in essence, all that is required to construct the permutation representations π : → S(K) is knowledge of all the sub-groups of , something that can be worked out by standard group theory methods. From those representations π we then have to select the ones whose image is contained in P A ⊂ S(K). Note that, while we are asking for the full representation R to be injective, this does not necessarily have to be the case for π .
The second step is to construct the maps r : → H A /G for each π and there are two technical complications we have to resolve in this context. First, we have to deal with the fact that r is, in general, a π -homomorphisms rather than a regular group homomorphism and, secondly, that r corresponds not to a linear but a (multi-)projective representation, since its target space is not simply H A but the quotient H A /G.
To disentangle these two problems, let us first consider the case where π is the trivial representation so that is actually a group homomorphism. There is a standard method to classify projective representations ρ : → Gl(n, C)/C * of a group using the group's Schur cover,ˆ . The Schur cover is a certain central extension of the original group and can be computed purely from group theory. Then, all projective representations ρ of can be obtained from linear representationsρ :ˆ → Gl(n, C) of the Schur cover by simply projecting to the quotient (see, for example, Ref. [16] or the summary in Appendix B.1.1 for details). This method can be trivially generalised to the multi-projective case, that is, to homomorphisms of the form Specifically, we can study all linear representationsρ :ˆ → K ∈K Gl(K , C) of the Schur cover group and then project to the quotient in order to obtain all such multiprojective representations ρ. These are not quite yet the representations r we are interested in, since the toric group G in Eq. (2.46) is usually a genuine sub-group of (C * ) |K| . However, we can deal with this complication by simply selecting from the representationsρ those representationsr which lead to well-defined homomorphisms upon taking the quotient by G.
Example. Let us illustrate the discussion so far with or standard example, using the group = Z 2 = {1, −1} for simplicity. Note that, from our previous index analysis, | | = 2 is indeed a possible order of a freely-acting group for this TCY manifold.
The permutation part, P A , of the ambient space symmetry group is isomorphic to S 2 so there are two possible choices for the representation π . In line with our simplifying assumption above we will first discuss the case where π is trivial. (The other case will be exemplified below.) The Schur cover of Z 2 is Z 2 so we can simply study linear representationsr : Z 2 → H A , for H A as in Eq. (2.41). There are many possible choices and one illustrative example, for the coordinate ordering ((1, 4), (2, 5), (3), (6), (7)), is Clearly,r (−1) chosen in this way normalises the toric group G in Eq. (2.39) and, hence, this representation descends to a well-defined homomorphism r : We are now ready to describe the algorithm for constructing all π -representations r in the general case. We start with one of the representations π : → P A constructed as described above and recall that P A can be viewed as a subgroup of S(K), the permutation group of the blocks K ∈ K. In general, there are various orbits under the action of π but in order to simplify the discussion let us assume that there is only a single orbit. The multi-orbit case is a straightforward generalization. Let us denote the blocks in this single orbit by K i ∈ K, where i = 1, . . . , b, and their stabiliser sub-groups under the representation π by i ⊂ . We can think of the representation R as being given in the following form where the r i (γ ) act on the blocks K i . The maps r i can now be restricted to i so that they induce projective representationsr i : 1 → Gl(K i , C)/C * . Focusing onr 1 , we can construct all such representations from linear representationsr 1 :ˆ 1 → Gl(K 1 , C) of the Schur coverˆ 1 , as discussed above. The point is now thatr 1 already determines the full π homomorphism r , as is explained in detail in Appendix B. Broadly, this works as follows. First choose γ i ∈ which map the blocks K i to the block K 1 under the representation π , that is, which satisfy π(γ i )(i) = 1. It can then be shown that, for every block i, γ ∈ can be written uniquely in the form where h ∈ 1 . Then it follows that so that all r i and, hence, r is fixed in terms ofr 1 .
Example. Coming back to our standard example, with group = Z 2 = {1, −1}, we would now like to discuss a case with a non-trivial choice of π . The only non-trivial choice is, in fact, π(1) = 1 7 and π(−1) =p, with the matrixp given in Eq. (2.43). There are three orbits under the action of π defined in this way, namely Let us start with the orbit K 1 . The stabilizer group 1 of the block {1, 4} is the trivial group and for the group elements γ i in Eq. (2.50) we can choose γ {1,4} = 1 and γ {2,5} = −1. This means that the correspondence between γ ∈ and h ∈ 1 established by Eq. (2.50) is given by 1 → 1 and −1 → 1. As a result, from Eq. (2.51), r must be trivial on the orbit K 1 . A similar argument shows that r must also be trivial on the orbit K 2 .
Finally, for the orbit K 3 the stabiliser group for the only block {3} in this orbit is the full group 1 = and we can, for example, choose −1 ∈ to act as −1 in this direction. This means the full homomorphism R for the ordering ((1, 4), (2, 5), (3), (6), (7)) is specified by with the matrixp in Eq. (2.43). Note that this action is not freely acting.
To conclude our discussion, we summarise the algorithm to construct the monomorphisms R : → G A .
• Find all representations π : G → P A , not necessarily faithful.
• Focus on one of the representations π and find all orbits of the blocks under its action.
• For each orbit pick one block, K 1 , and consider its stabiliser group 1 .
• Study all projective representationsr 1 : 1 → Gl(K , C)/C * by using the Schur coverˆ 1 of 1 . • For each such representationr 1 , re-construct r on each orbit and, hence, the full representation r . • Check if the combination of π and r does indeed define a monomorphism R : → G A .

Smoothness of hypersurfaces.
A group monomorphism R : → G A , constructed as described in the previous sub-section, provides an automorphism of the ambient space A and it defines an action on the homogeneous coordinates x i of A. For R to give rise to an automorphism of the TCY manifold X ⊂ A the defining polynomial p of X needs to be invariant under the action of R. The general invariant polynomial p R( ) under the action of R can be found from a generic defining polynomial p of X by applying the Reynolds operator to p, that is, by In general, the invariant polynomials p R( ) have fewer independent coefficients than a general p and it is not guaranteed that the generic TCY space X R( ) defined by p R( ) is smooth. In the following, we devise a method to check smoothness for a TCY space X defined as the zero locus of a polynomial p. So far, we have used the global description of a toric variety in terms of homogeneous coordinates. Alternatively, a toric variety may be glued together from patches, one for each cone of maximal dimension (see Appendix C for details). Our smoothness check will be carried in this latter description of a toric variety, proceeding patch by patch.
Concretely, to each four-dimensional cone σ in the fan we can associate a patch V (σ ) with affine coordinates x σ i , which correspond to the rays v i of σ . The defining polynomial, p σ , on this patch can be obtained from p by setting the coordinates x i = x σ i for v i ∈ σ and setting all other coordinates to one. For the toric varieties we are considering, the patches V (σ ) are either of the form C 4 or C 4 /G for a finite Abelian group G. In the latter case, the coordinates x σ i should be interpreted as "upstairs" coordinates on C 4 . In the case V (σ ) = C 4 the patch is clearly non-singular. For V (σ ) = C 4 /G and for fans constructed from fine triangulations of reflexive four-dimensional polytopes, the action of G is such that it has fixed points only at the origin, x σ i = 0, of the patch. As a practical matter, a patch V (σ ) is smooth precisely if the matrix formed by v i ∈ σ has determinant ±1. For all non-smooth patches identified in this way, we have to check whether or not the origin lies on the TCY hypersurface, that is, if p σ (0) = 0. If it is for at least one patch the TCY hypersurface is singular and we can stop.
Otherwise, we can proceed and check if the defining equation p gives rise to singularities. This is done patch by patch, considering the Jacobi ideals Example. We would like to check smoothness for our standard example with symmetry = Z 2 and the symmetry action which contains only the monomials (2.14) which correspond to the vertices of the polytope . We first note that there are four singular patches on the ambient space which correspond to the cones (2.57) As none of the associated loci meets our hypersurface (2.56), these singularities do not induce singularities on X R( ) . To see that there are no singularities induced by the hypersurface equation (2.56), note that this equation and all its partial derivatives can only vanish in one of the patches when x 3 = 0. As the whole situation is furthermore symmetric under swapping any of x 1 ↔ x 4 , x 2 ↔ x 5 or x 6 ↔ x 7 , we only need to check a single four-dimensional cone, which we can take to be the one generated by In the corresponding patch, the hypersurface equation becomes It can be checked that the Jacobi ideal, I σ , of this equation has dimension −1 so that X R( ) is indeed smooth.
Our method for checking smoothness can, hence, be summarised as follows: • Determine the most general invariant family of R( )-invariant hypersurfaces X R( ) by finding the corresponding polynomials p R( ) from Eq. (2.53). • For each four-dimensional cone, σ ∈ , with rays v i introduce affine coordinates x σ i and the affine version, p σ , of the polynomial p R( ) on the patch V (σ ).
• For each four-dimensional cone σ of carry out the following. − To set this up, we denote by q : B → A the projection to the quotient space. We also require liftsˆ andĜ A of the groups and G A to the upstairs space, together with projections :ˆ → and ν :Ĝ A → G A and the upstairs versionR :ˆ →Ĝ A of the homomorphism R such that R • = ν •R. (In practice, G A andĜ A are generated by the same matrices, where the projective relations are taken into account for G A and ignored forĜ A .) Writing γ = (γ ), we have (2.59) For γ ∈ we define the downstairs fixed point set as We would like to compare this set with its upstairs counterpart forγ ∈ˆ , where that is, we obtain the downstairs fixed point set simply by projecting its upstairs counterpart to the quotient. Of course we are not primarily interested in the fixed point set of R( ) in A but in its intersection with the TCY manifold X ⊂ A. We denote this intersection by F γ := F γ ∩ X . ByX ⊂ B we denote a lift of the TCY manifold such that q(X ) = X . Then, the upstairs version of the fixed point set intersected withX given byFγ =Fγ ∩X projects down to F γ , that is (2.63) What we want is for the Calabi-Yau manifold not to intersect the fixed point set, that is, we require F γ to be empty for all γ ∈ . From the previous equation this is the same as saying thatFγ is empty for allγ ∈ˆ . Therefore, from Eq. (2.61), we should, in a first instance, solve the equationsR is an element of the toric group G, parametrized s = (s 1 , . . . , s n−4 ). We can think of these equations as defining an ideal I in the ring C[z, s].
(It may be necessary to multiply some components of Eq. (2.64) with powers of some parameters s r if the corresponding charges q i r are negative, in order to achieve polynomial form.) To this ideal we have to add the defining TCY equation in order to get the ideal which describedFγ . Then, in order for the action to be fixed point free we need the dimension of this ideal to be −1 for allγ ∈ˆ .
Example. Finally, let us check for our standard example that the action of Z 2 symmetry given in Eq. (2.55) is indeed fixed point free.
For this, we need to find all fixed loci on the ambient space and intersect with the hypersurface equation p R( ) = 0, with p R( ) given in Eq. (2.56). The fixed loci of the ambient space are x 3 = x 4 = x 5 = x 6 = 0, together with their images under any permutation x 1 ↔ x 4 , x 2 ↔ x 5 or x 6 ↔ x 7 . Since the situation is symmetric under these permutations, it again suffices to check one such locus only. Solving p R( ) = x 3 = x 4 = x 5 = x 6 = 0 we find that x 1 = 0 or x 7 = 0, a locus which is contained within in the exceptional set Z ( ). In conclusion, the action of the symmetry (2.55) is indeed fixed point free. An implementation of the above algorithm requires the data of a fan obtained from a triangulation of • for each model. We are using the database compiled in [11] as input data containing such models. Furthermore, we are using GAP [17] for the required group-theoretical data. Our computation system is composed of Mathematica modules, and the computations in commutative algebra are performed using Singular [18]. The Calabi-Yau hypersurfaces are generated with random coefficients c m for each monomial ranging from 0 to 100. If the computation in Singular is too time-consuming over the field of rational numbers, we alternatively use finite fields and check over a number of primes.

Examples
We will now carry out our algorithm for two specific examples. The first example is for a freely-acting Z 4 symmetry on the tetra-quadric hypersurface in P 1 ×P 1 ×P 1 ×P 1 . Since the tetra-quadric is a TCY as well as a CICY manifold, this Z 4 symmetry is, in fact, well-known in the literature [8]. This example serves as an illustration of our method in an established context. The second example is for a TCY manifold with h 1,1 (X ) = 6 and a freely-acting Z 2 symmetry which was not previously known.

A Z 4 symmetry of the tetra-quadric.
The tetra-quadric with Hodge numbers h 1,1(X ) = 4 and h 1,2 (X ) = 68 is described by a polytope • with eight vertices given by The zero set is then characterised by The linear relations between the above vertices leads to the charge matrix The greatest common divisor of these numbers and the Euler number is 16 so for freelyacting symmetry groups we should consider the group orders | | ∈ {2, 4, 8, 16}.
Step II: In order to compute the ambient space symmetry group we need to determine the index sets J and K from the above toric data. In this particular case it turn out they are both equal to I so that Following the general procedure explained in Sect. 2.3 it is then straightforward to determine the ambient space symmetry group G A = P A (H A /G). One finds that where P A acts by permuting the four blocks {1, 8}, {2, 7}, {3, 6}, {4, 5} ∈ K.
First, we need to construct the permutation representations π : Z 4 → P A ∼ = S 4 . For illustration purpose we choose π as . We can focus on K 1 (the other orbit K 2 will work in exactly the same way) and note that it has two blocks K 1 = {1, 8} and K 2 = {2, 7}. The stabilizer groups 1 of K 1 is given by 1 We choose the representationr 1 of 1 on the block K 1 asr 1 (γ 2 ) = σ 3 = diag(1, −1). (3.11) The full map r : → H A /G can then be re-constructed fromr 1 and is it straightforward to show that r (γ ) = diag(1 2 , σ 3 , 1 2 , σ 3 ). (3.12) Combining this with the permutation representation (3.10) we finally get The general invariant polynomial under this symmetry action is given in Appendix D.9.
Step IV: To check fixed points we have to write down Eq. (2.64) for all Z 4 group elements using the group action defined by Eq. (3.13) and the elements (3.5) of the toric group. This leads to It can be shown in each case that the solutions to these equations are either contained in the exceptional set (3.3) or are point-like in the ambient space and, hence, do not lie on a generic tetra-quadric hypersurface. The most general defining polynomial p consistent with the Z 4 symmetry contains 21 independent terms and coefficients, including, for example A calculation, for example using Singular [18], shows that a sufficiently general combination of these terms does indeed produce a smooth hypersurface.

3.2.
A Z 2 symmetry of the twisted tetra-quadric. We now proceed to a more complicated example which leads to a freely-acting Z 2 symmetry which was not previously known. Moreover, the underlying manifold is not among the 16 manifolds identified by Batyrev and Kreuzer and, hence, this example shows that the number of manifolds with freelyacting symmetries within the Kreuzer Skarke list is definitely larger than 16. The space is defined by a polytope • with the 10 vertices (1, 1, 1, 1). Step The greatest common divisor of these numbers and the Euler number is 10 so for freelyacting symmetry groups we should consider the group orders | | ∈ {2, 5, 10}.
Step III: Here, we need to study all group monomorphisms R : → G A for all finite groups | | ∈ {2, 5, 10}. For the purpose of this example, we focus on the group = Z 2 and we denote its generator by γ . A simple choice of representation is to take r to be trivial and to set the permutation part to π(γ ) = ((1, 10), (2,9), (3,8), (4,7), (5, 6)). (3.23) Step IV: As described earlier, a generic polynomial symmetric under the Z 2 (3.23) may easily be found by summing a generic polynomial with its image, the result is found in Appendix D. 10. Such a polynomial gives rise to a smooth hypersurface X R(Z 2 ) . To see the fixed point freedom, we have to find solutions to (2.64). For the action above, fixed points are found as solutions of To find the corresponding loci on the toric variety, we may consider the elimination ideal with respect to the variables s i . It is generated by all polynomials of the form x 2 i x 2 j − x 2 10−i x 2 10− j for i, j = 1 . . . 5 and has dimension 5, so that there are only fixed points under this group action. As none of these points meet a generic smooth symmetric polynomial, it follows that (3.24) has no simultaneous solutions with a generic smooth symmetric polynomial and there are not fixed points under the group action (3.23).

Systematic Search for Symmetries
The classification of reflexive four-dimensional polytopes by Kreuzer and Skarke [10] forms a natural starting point to apply our algorithm and find new examples of Calabi-Yau threefolds with freely acting symmetries. A scan over all triangulations off all reflexive polytopes, however, is presently beyond reach. In particular, it is already a daunting task to find all triangulations, as the number of triangulations of each reflexive polytope increases dramatically with the Picard number. In Ref. [11] all triangulations for h 1,1 (X ) ≤ 6 were found. 2 For the purpose of our symmetry classification we will focus on the subset of those triangulations with h 1,1 (X ) ≤ 3, some 350 manifolds to which we apply our classification algorithm.

Scan results for freely-acting symmetries.
In this sub-section we present the result of our classification, a list of all freely-acting symmetries, linearly realised in the ambient space, for TCY manifolds with h 1,1 (X ) ≤ 3. We find that only five of these TCY manifolds allow for freely-acting symmetries of this kind and we briefly list those five cases, together with the most important data, in the following. The invariant polynomials are recorded in Appendix D. Case #1, quintic in P 4 : • Hodge numbers: h 1,1 (X ) = 1, h 2,1 (X ) = 101 • Possible symmetry orders from indices: | | ∈ {5, 25} • Ray generators of fan : • Charge matrix: in all four cases. • Remarks: This is polytope ID# 2 of the database of [11]. The quintic is of course a CICY manifold and all the above freely-acting symmetries are well-known. The Z 5 symmetry is toric and among the 16 cases found by Batyrev and Kreuzer, while the four Z 5 × Z 5 actions are non-toric.
Case #2, bi-cubic in P 2 × P 2 : • Symmetry = Z 3 = γ : in all four cases. • Remarks: This is polytope ID# 10 of the database of [11]. The bi-cubic is of course a CICY manifold and all the above freely-acting symmetries have been found before [8]. The Z 3 symmetry is toric and among the 16 cases found by Batyrev and Kreuzer, while the four Z 3 × Z 3 actions are non-toric.
• Charge matrix: • Symmetry = Z 2 = γ : • Remarks: This is polytope ID# 245 of the database of [11]. This manifold with its freely-acting Z 2 symmetry is among the 16 cases identified by Batyrev and Kreuzer. The polytope has two triangulations realizing different regions of the Kähler cone. They both allow the same fixed-point free actions as given above. However, the Z 2 Z 2 is non-toric and has not been found before. We will take a closer look at this symmetry in the following sub-section. Case #5: • Hodge numbers: h 1,1 (X ) = 3, h 2,1 (X ) = 115 • Possible symmetry orders from indices: | | ∈ {2, 4, 8} • Rays of fan : • Charge matrix: • Symmetry = Z 2 = γ : • Remarks: This is polytope ID# 249 of the database of [11]. This manifold with its freely-acting Z 2 symmetry is among the 16 cases identified by Batyrev and Kreuzer. The polytope has three different triangulations realizing different regions of the Kähler cone. All of them allow the same fixed-point free action given above.

4.2.
A closer look at the Z 2 Z 2 symmetry in case #4. In the previous classification, we have found one new case, not previously known from the classification of symmetries for CICYs or the Batyrev Kreuzer classification. We would now like to discuss this case in more detail. We recall that the vertices of • for this TCY manifold are given by Eq. (4.14), the associated charge matrix by (4.15) and that the Hodge numbers are h 1,1 (X ) = 3 and h 2,1 (X ) = 115. We consider the following triangulation 3 which consists of twelve four-dimensional cones and leads to the zero set The charge matrix (4.15) implies the following toric group: Step The greatest common divisors of the numbers and the Euler number is 8, so the possible orders of freely-acting groups are ∈ {2, 4, 8}.
Step II: From the index set I in Eq. while P A is trivial.
Step III: In general, we should consider all groups of order | | ∈ {2, 4, 8} but for the present discussion we focus on the order four group = Z 2 Z 2 = γ 1 , γ 2 . There are numerous possible group monomorphisms R : → G A and we consider the specific example with Note that, while these two matrices generate a group of order eight, this is reduced to order four once projective identifications are taken into account.
Step IV: The most general TCY manifold of this type is defined by linear combination of 153 monomials (corresponding to the 153 lattice points on ) which, in line with Eq. (2.7), correspond to the points in the polar polytope . These monomials have the form (4.29) Requiring invariance under the action of the symmetry (4.28) reduces this list to 40 independent invariant polynomials of the form It can be checked that a generic linear combination p R( ) of these 40 invariant building blocks leads to a smooth TCY manifold X (R( ) . In order to check fixed point freedom, we focus on the action of R(γ 2 ). With the toric group (4.24), Eq. (2.64) leads to the following fixed point equations for R(γ 2 ). After excluding the zero set (4.23), it can be checked that no component of the so-defined fixed point set intersects the generic TCY manifold X R( ) . The same conclusion can be reached for the action of the other group elements.

Conclusion
In this paper, we have analysed symmetry groups of CY three-folds defined as hypersurfaces in toric ambient spaces (TCY manifolds) and we have developed methods to classify those symmetries which descend from a linear action on the toric ambient spaces. Based on these results, we have set up an algorithm that allows for the classification of symmetries for TCY manifolds obtained by triangulation from the Kreuzer-Skarke list.
We have carried out this classification for a small initial sample of all TCY manifolds with h 1,1 (X ) ≤ 3, some 350 manifolds. We find that only 5 out of these 350 TCY manifolds admit (linearly realised) freelyacting symmetries. Our algorithm reproduces all freely-acting symmetries known previously from either the overlap with the list of CICY manifolds [8] or the classification of toric symmetries by Batyrev and Kreuzer [12]. We find one new, non-toric Z 2 Z 2 symmetry, where one of the Z 2 factors corresponds to one of the 16 symmetries identified by Batyrev and Kreuzer. We have also used our methods to construct a new, freely-acting Z 2 symmetry on a TCY manifold with h 1,1 (X ) = 6 which was not previously known to have freely-acting symmetries. This example shows that the subset of TCY manifolds with freely-acting symmetries is larger than the 16 cases identified by Batyrev and Kreuzer. From our preliminary evidence it appears that the fraction of TCY manifolds with freely-acting (linearly realised) symmetries is relatively small, namely five manifolds out of 350. This figure is broadly in line with what was found in the context of CICY manifolds [8], where about 2.5% of the 7890 types carry linearly realised, freely-acting symmetries.
Naturally, our algorithm can also be used to find non-freely acting symmetries. There is no a-priori constraint on the group order anymore but we can simply search all finite groups up to a maximal order and, of course, omit the fixed point check from the algorithm.
It would be desirable to apply our algorithm to a larger class of TCY manifolds. As it stands this is possible for TCY manifolds with h 1,1 (X ) ≤ 6 since the corresponding triangulations have been worked out in Ref. [11]. To extend our analysis to the entire Kreuzer-Skarke list of about half a billion reflexive polytopes is a formidable challenge, not least because working out the triangulations for the entire list is computationally difficult. Dealing with this will certainly require some modifications and refinements of our algorithm. For example, it is conceivable that necessary conditions for the existence of freely-acting symmetries can be used to select a subset of polytopes for triangulation. These and related issues are the subject of future work.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

A The Automorphism Group of a Toric Variety
A.1 Demazure's structure theorem. The automorphism group of a toric variety is elegantly captured in terms of Demazure's structure theorem [19,20]. It states that the Lie algebra of the automorphism group of a toric variety is generated by the action of the open dense algebraic torus giving it its name, as well as maps where the α iν are lattice points in M for which α iν , v i = −1 and α iν , v k > −1 for all k = i. The α iν are called Demazure roots. We denote the group generated by (A.1) and the action of the algebraic torus by Aut 0 (P ).
The full automorphism group is recovered as follows [19,20]. For any automorphism of the fan, i.e. maps in Gl(N ) which preserve the set of cones in , there is an associated permutation of homogeneous coordinates which gives rise to an automorphism of P . Let us denote the group of all such maps by Aut (N , ). This automorphism is trivial if it is a Weyl reflection, defined by for any pair of semi-simple roots α (a root of v i (α)) and −α (a root of v j (−α). Then, the following relation holds: Note that the Demazure roots, and hence the continuous part of the automorphism group, is completely independent of the triangulation, as it only depends on the one-dimensional cones. We are going to show in section A.6.2 that this remains true when we restrict to the Lie algebra of linearly realized automorphisms. In the following, we will prove a statement analogous to Demazure's structure theorem for automorphisms linearly realized on the homogeneous coordinates of a toric variety which is suited for our classification of symmetries. This is based on the "global" presentation of a toric variety A as the quotient A.2 Zero set and permutation group. We start with an analysis of the symmetries of the zero set Z ( ) of a toric variety corresponding to the index set I. Its elements are sets I ⊂ {1, . . . , n} which encode the components of the zero set as defined in Eq. (2.2). For this index set we define the following permutations groups: (A.5) The group S leaves all sets I in I invariant individually and we refer to it as the stabiliser group while R permutes the sets I into each other. Clearly, S is a sub-group of R. It is actually a normal sub-group since for s ∈ S and r ∈ R we have rsr −1 (I ) = I and, hence, rsr −1 ∈ S. We can, therefore, define the quotient group which can be thought of as the group of permutations on I, induced from permutations in S n . We can write down the short exact group sequence where ı is the inclusion map and [·] denotes taking the equivalence class. We now show that this sequence splits. We begin with the as a disjoint union. Note that for a J ∈ J and I ∈ I with J ∩ I = ∅ it follows that J ⊂ I . Indeed, if j ∈ J ∩ I and k ∈ J then k ∼ j and, hence, k and j must be contained in the same sets of I. Since j ∈ J it follow that k ∈ J , so that J ⊂ I . This means every I ∈ I can be written as a disjoint union Proof. For the first part consider a set J ∈ J , a permutation s ∈ S and any j ∈ J . Assume that i is contained in precisely the sets I 1 , . . . , I p ∈ I. Hence, s(i) ∈ s(I k ) = I k , so that s( j) is also contained in I 1 , . . . , I p ∈ I. Also, s( j) cannot be contained in any other I ∈ I, which is not among I 1 , . . . , I p ∈ I, since j would have to be in I in this case. This implies that i ∼ s( j), so s( j) ∈ J and, hence s(J ) ⊂ J . Since s is injective and J is a finite set this means s(J ) = J . For the second part, again focus on a J ∈ J and a permutation r ∈ R. Pick a j ∈ J and callJ ∈ J the equivalence class of the image r ( j). I want to show that r (J ) ⊂J . Let k ∈ J , so that j and k are contained in the same sets I 1 , . . . , I p ∈ I. Then, both r ( j) and r (k) are contained in r (I 1 ), . . . , r (I p ) ∈ I and in no further I ∈ I. This means that r (k) ∼ r ( j), so r (k) ∈J and, hence, r (J ) ⊂J . Further, every l ∈J is the image under r of r −1 (l) ∈ J , so that r (J ) =J .
We can now understand the structure of the stabilizer group S. Returning to the general story, we call a permutation r ∈ R order-preserving iff for all J ∈ J and all i, j ∈ J with i < j it follows that r (i) < r ( j). In other words, such permutations preserve the natural ordering of the numbers 1, . . . , n but only within each set J ∈ J . The set of all order-preserving permutations forms a sub-group of R which we denote by P .

Lemma A.3. Each equivalence class p ∈ P contains exactly one order-preserving permutation in P .
Proof. For existence, consider an arbitrary representative s ∈ p of an equivalence class p ∈ P and consider the restriction s| J : J → s(J ) to any J ∈ J . We can find a permutation r J ∈ S(J ) of J such that s| J • r J has the order-preserving property. The maps r J : J → J combine into permutation a r ∈ S n which, from Theorem A.2 is an element of the stabilizer group S. Then, s • r ∈ p and it is order-preserving. For uniqueness, consider two permutations s,s ∈ p which are both order-preserving. This means than r ≡ s •s −1 is in S that is, r | J ∈ S(J ) is a permutation of J . For i, j ∈ J with i < j we have r (i) < r ( j) which means that r | J = id J . This is true for all J ∈ J so that r = id and, hence, s =s.
This result means that we can define a group isomorphism τ : P → P ⊂ R which assigns to a class p ∈ P the unique order-preserving permutation τ ( p) ∈ p. Of course, [τ ( p)] = p, so τ splits the sequence (A.7). This means that R is a semi-direct product with the isomorphism ϕ : P × S → R defined by ϕ(( p, s)) = τ ( p)s and the multiplication in P S given by Hence, we can determine P ∼ = P by finding the permutations in S n which map the sets I into each other (that is, which are elements of R) and are order-preserving. The group R can then be obtained from S and P by forming the semi-direct product (A.11).

A.3 Invariance group of the upstairs space.
We now discuss the structure of G B of the upstairs space B = C n − Z ( ). First note that there is an obvious embedding S n → Gl(n, C) defined by σ (e i ) = e σ (i) , with the standard unit vectors e i . Hence, the above subgroups S, R, P ⊂ S n have isomorphic images in Gl(n, C) which we denote by S B , R B , P B , respectively. In fact, from the definition of S, R, P it is clear that  (I )). Hence, the r I can be combined to a map r ∈ Gl(n, C) with r | Z (I ) = r I which permutes the standard basis vectors e i and maps zero sets Z (I ) into each other, so r ∈ R B . Further, definingh = r −1 • g it follows thath ∈ H B . From the previous section, we know that we can write r = ps, where p ∈ P B and s ∈ S B . Finally, with h ≡ sh ∈ H B we have g = ph, where p ∈ P B and h ∈ H B .
To show that H B is a normal sub-group of G B , consider a η ∈ H B and a g = ph ∈ G B . Then Finally, from the definition of R it is clear that R B ∩ H B = S B . But S B contains precisely one order-preserving element, the identity.
The result of the previous Lemma shows that G B is a semi-direct product G B ∼ = P H B .
If we denote by φ : P → P B the isomorphism between order-preserving permutations in P and corresponding Gl(n, C) matrices, then a pair ( p, h) ∈ P H B is mapped to φ( p)h ∈ G B and the multiplication in P H B is (A.14) We have already understood the permutation group P and how to compute it from the basic data. It remains to consider H B . To this end, we define the group where Gl(J, C) is the general linear group acting in the directions of the coordinates Recall that the sets J were defined around Eq. (A.8). The elements of (A.15) can be viewed as block-diagonal n × n matrices with the blocks corresponding to the coordinates x J . Unfortunately, H B can be larger than G(J ). In general we know that the zero sets Z (I ) and their intersections Z (I 1 ) ∩ · · · ∩ Z (I p ) = Z (I 1 ∪ · · · ∪ I p ) have to be invariant under transformations h = (h i j ) ∈ H B . This implies, for any union I 1 ∪ · · · ∪ I p , that We call a zero set I and associated set J regular if for every J ∈ J and all i ∈ J , j ∈ J c there exists an I ∈ I with i ∈ I but j / ∈ I . Then we have: Proof. Focus on a specific J ∈ J and fix a i ∈ J and a j ∈ J c . From regularity there exists an I ∈ I such that i ∈ I and j ∈ I c . Applying (A.16) for I and I c means that h i j = 0. Since i and j were arbitrary this means that h i j = 0 for all i ∈ J and all j ∈ J c . Now consider an i ∈ J c and a j ∈ J . We know that J c is the disjoint union of the other J , so i ∈J ⊂ J c for a particularJ ∈ J . Regularity, applied toJ , means that there exists an I ∈ I with i ∈ I and j ∈ I c . As before, (A.16) applied to I leads to h i j = 0. Since i, j were arbitrary this means that h i j = 0 for all i ∈ J c and all j ∈ J . In summary, for all J ∈ J we have h i j = 0 if i ∈ J and j ∈ J c or if i ∈ J c and j ∈ J . This means that h ∈ G(J ), so H B ⊂ G(J ). The opposite inclusion has already been shown in the previous Lemma, so G(J ) = H B .
Therefore the complete symmetry group of B = C n − Z ( ) is given by A.4 Normalizer in GL(n, C). In order to find the symmetry group G A of the toric space we need to compute the normalizer of the toric group G within G B , the symmetry group of the upstairs space. This is relatively easy if the group elements of G are proportional to the unit matrix within each block (A.15) of G B . In order to sort this out it is instructive to first deal with a related -but simpler -problem, namely to compute the normalizer of G within Gl(n, C).
To this end it is useful to split the various coordinate directions up into disjoint blocks K = {K | K ⊂ {1, . . . , n}}, collecting the directions with the same charges q i , such that where G q K (K ) = {s q K 1 |K | | s ∈ (C * ) n−d } consists of matrices proportional to the unit matrix and by q K we denote the charge in the directions in K . As usual, we think of the symmetric group S n as being embedded into Gl(n, C) via σ (e i ) = e σ (i) for σ ∈ S n . We would like to work out the normalizer N Gl(n,C) (G) = {g ∈ Gl(n, C) | gG = Gg} = {g ∈ Gl(n, C) | ∀γ ∈ G ∃γ ∈ G : gγ =γ g}.
(A.19) It is also useful to introduce the commutant of G within Gl(n, C) as well as the sub-groups of S n defined by They consist of permutations which normalize or commute G, respectively. We start by characterizing the groups R and S in a different way which is more useful for practical calculations.
Proof. i) Consider a permutation σ ∈ S n which normalizes G so that for all γ ∈ G there exists aγ ∈ G such that γ σ = σγ with γ = diag(s q i ) andγ = diag(s q i ). It follows thats q i = s q σ (i) for all i, so thatγ = diag(s s σ (i) ). Consider all i ∈ K for a given block K , so that q i = q K . Then s q σ (i) =s q K for all i ∈ K . This only works if all q σ (i) = q K are equal, so that σ (i) ∈ K for all i ∈ K . Hence, σ (K ) ⊂ K ∈ K.
Applying the same argument to σ −1 and K leads to σ −1 (K ) ⊂ K , so combined we have σ (K ) = K . ii) Clearly, permutations which only act within blocks are in the commutant. Conversely, permutations in the commutant have to be block-diagonal from Schur's Lemma.
So, S is the set of permutations of directions within each block andR contains permutations of blocks and within blocks. Unfortunately, the group R we are actually interested in can be genuinely smaller thanR. If there are more blocks than variables s it might happen that a permutation of these blocks cannot be compensated by a choice ofs in the normalizer condition. Clearly, S ⊂ R is a sub-group and indeed a normal sub-group so that we can define the quotient group There is a monomorphism P → R which involves assigning to an element of P the element of R which permutes the blocks in K in the same way and preserves the natural order in each block K ∈ K. In this way, We are now ready for the following theorem.
Theorem A.9. The normalizer can be expressed in terms of the commutant as N Gl(n,C) (G) = P C Gl(n,C) (G), C Gl(n,C) (G) is a normal sub-group of N Gl(n,C) (G) and P ∩ C Gl(n,C) (G) = 1.
Proof. We begin by showing that every matrix g =∈ N Gl(n,C) (G) can be written as a product of a permutation and a matrix in the commutant. First, g = (g i j ) ∈ N Gl(n,C) (G) means that for all γ ∈ G there exists aγ ∈ G such that gγ =γ g. Write γ = diag(χ i ) and γ = diag(χ i ) where χ i = s q i andχ i =s q i . Then it follows that g i j (χ i −χ j ) = 0 for all i, j. The matrix g is invertible which means that for every row i there exists a column j with g i j = 0 and, moreover, for each row we can choose a different column. This means there exists a permutation σ ∈ S n so that g iσ −1 (i) = 0 for all i. In order to satisfy the normalized condition we then need thatχ i = χ σ −1 (i) . Now define c = σ −1 g so that g = σ c. The normalizer condition then reads cγ c −1 = σ −1γ σ = diag(χ σ (i) ) = diag(χ i ) = γ and, hence, c ∈ C Gl(n,C) (G) is in the commutant. This shows that every g in the normalizer can be written as g = σ c with σ ∈ S n and c in the commutant. With this decomposition, the normalizer condition, gγ =γ g becomes σ γ =γ σ which shows that, in fact, σ ∈ R. Further we know that every r ∈ R can be written as r = ps, where p ∈ P and s ∈ S, so that g = rc = psc. But both s, c and their product sc are in the commutant. This proves the decomposition N Gl(n,C) (G) = P C Gl(n,C) (G).
We now need to show that C Gl(n,C) (G) is a normal sub-group of N Gl(n,C) (G). Choose a c in the commutant and g = pc in the normalizer. Then gcg −1 = pc p −1 wherẽ c = c cc −1 is in the commutant, so we need to show that pc p −1 is also in the commutant, that is pc p −1 γ = γ pc p −1 for all γ ∈ G. This is equivalent tocγ =γc wherẽ γ = p −1 γ p. However,γ has the same block-structure than elements in G (as the action of p permutes the blocks) so thatc andγ indeed commute. Finally, since the action p −1 γ p of a p ∈ P on γ ∈ G permutes the blocks of γ , it follows that the only p ∈ P for which p −1 γ p = γ is, in fact, p = 1, so that N Gl(n,C) (G) and P ∩ C Gl(n,C) (G) = 1.
So, in summary, we have found that the normalizer can be expressed in terms of the commutant as N Gl(n,C) (G) = P C Gl(n,C) (G) , (A.24) and the latter, by means of Schur's Lemma, can be written as A.5 Invariance group of the toric ambient space. In order to find the symmetry group G A of the toric ambient space A we have to compute the normalizer N G B (G). Clearly, this normalizer is given by The two groups on the right-hand side have been explicitly determined in the previous two sections, so our remaining task is to express their intersection in a useful form. Both groups have a similar structure in that they relate to a block-decomposition of the n coordinates and consist of a semi-direct product of a permutation group which exchanges the blocks times block-diagonal matrices which generate general linear transformations within each block. However, the block-decomposition is in general different for the two cases. The group G B relates to the block-decomposition J which follows from the structure of the zero set and is given by The group N Gl(n,C) (G), on the other hand, relates to the block-decomposition K which follows from the toric group G and is given by The complication is that the block-decompositions J and K are not necessarily identical. It is, therefore, useful to define the refined block-decomposition given as the intersection of J and K defined by 4 We introduce the usual collection of groups associated to this block-decomposition. and we expect this to be the continuous part of N G B (G). We also define the intersection of the permutation groups The two sets of discrete groups relate in the following way.
Proof. i) Consider a r ∈ R ∩ R so that r (J ) ∈ J for all J ∈ J and r (K ) ∈ K for all K ∈ K. Then, for any L = J ∩ K ∈ L we have r (L) = r (J ∩ K ) = r (J ) ∩ r (K ) so that r (L) ∈ L. Hence, R ∩ R ⊂ R L . Proof. "⊃": Consider a g = rh ∈ R A H A , where r ∈ R A = R ∩ R and h ∈ H A = H B ∩ C Gl(n,C) (G). It follows that rh ∈ R H B = G B and rh ∈ RC Gl(n,C) = N Gl(n,C) (G) and, hence, g = rh ∈ N G B (G). "⊂": Start with a g ∈ N G B (G) = G B ∩ N Gl(n,C) (G) which can be written as g = ph = πc with p ∈ P,h ∈ H B , π ∈ P,c ∈ C Gl(n,C) (G). It follows that π −1 p =ch −1 ∈ C Gl(n,C) (G)H B . Since π −1 p is a permutation it follows that π −1 p ∈ S S which is the set of permutations in C Gl(n,C) (G)H B . Write π −1 p = σ s −1 , where σ ∈ S and s ∈ S, so that ps = σ π ≡ γ . Since ps ∈ R and πσ ∈ R it follows that γ ∈ It is clear that we can write N G B (G) as a semi-direct product by dividing out the group S A . To find the invariance group, G A , of the toric space we have to divide this by G which results in This completes the calculation of G A .
A.6 Toric varieties obtained from reflexive polytopes. The results outlined above hold for any toric variety constructed from a simplicial fan. Here, however, we are interested in toric varieties which are associated with triangulations of reflexive polytopes, which gives extra structure. Let us first recall the definition of Demazure roots, (A.1). In the context of reflexive polytopes, it directly follows that only x i corresponding to vertices v i can have Demazure roots, which are in turn given by all interior points 5 of the dual faces [3] i . Restricting to actions which act linearly on the homogeneous coordinates we recover the theorem stated in section A.6.1. The set of all roots α for which −α is a root as well are known as semi-simple roots in the literature.
A.6. 1 The set K. Let us review the set K. Its elements K are sets containing a collection of indices {1, .., n} such that q i = q j for all i, j in a specific K . As we discuss now, this has fairly strong implications for any two i, j in the same K .
Theorem A.12. On a simply connected toric variety, any two lattice points v i and v j on • for which q i = q j are vertices of • . Furthermore, there exist points α i j inside of [3] i as well as α ji inside of [3] j for which α i j = −α ji .
Proof. From the equality of q i and q j it follows that the corresponding divisors D i and D j are linearly equivalent. As all linear equivalences take the form for some vector α i j in M ⊗ R, there must be a vector α i j such that For a simply connected toric variety, the set of all vectors v i generate the lattice N . Hence a first consequence of this is that α i j is contained in the dual lattice M. Furthermore, the above implies that it is a lattice point on , as and v j are contained in a threedimensional hyperplane of N ⊗ R passing through the origin. As • is four-dimensional and the convex hull of all of the v l , it follows that v i and v j must be above and below this hyperplane, so that they must be vertices. In this case (A.38) implies that α i j is inside the dual three-dimensional face of v i , [3] i and α ji = −α i j is inside [3] j .
The above provides an alternative algorithm to determine the set K. All indices l corresponding to non-vertices v l form sets K l = {l}. For any vertex v i , we have to find the interior points of the dual face [3] i . If the negative of such an interior point is also inside of a three-dimensional face [3] j (defining a dual vertex v j of • ) there is a set K ∈ K containing both i and j. This defines an equivalence relation ∼ K and the equivalence classes are precisely the sets K ∈ K.
This way of thinking is particularly useful for polytopes • with many integral points, where many points (all non-vertices) can never occur in any non-trivial set K ∈ K.
A.6. 2 The set J . The data exploited in the last paragraph only depends on the structure of the polytopes • and . Let us now consider the triangulation data. We have described the SR ideal by an index sets I ∈ I and have constructed the refinement J of I as the set of equivalence classes under the relation j ∼ J i ↔ i and j are contained in the same sets I ∈ I. (A.39) In other words, J is obtained by first forming all intersections between sets in I and then discarding sets which are already containing smaller sets.
Lemma A.13. If all cones in a fan are symmetrical with respect to exchanging i ↔ j then i ∼ J j.
Proof. To see this, all we have to do is go through the construction of J from the data of cones. We first find the SR ideal by forming the complement of the set of cones in the power set of {1, · · · , n}. By construction this will be symmetrical under i ↔ j. But then also the set of generators I enjoys the same symmetry so that i and j are in the equivalence class in J as claimed.
Theorem A.14. If D i = D j , so that i ∼ K j, it follows that i ∼ J j.
Proof. Due to the lemma, it is sufficient to show that any triangulation of the polytope • is such that the collection of cones is symmetric under the exchange i ↔ j. First of all, the triangulation tr( • ) is completely fixed in terms of its cones of maximal (i.e. 4) dimension. Every one of these cones is spanned by four linearly independent vectors in N . If such a cone contains both v i and v j among its generators, there is nothing to show, as this is symmetric in i and j. A cone which just contains v i but not v j is hence spanned by v i and three vectors v k 1 , v k 2 , v k 3 in the plane orthogonal to α i j . This means the cone in question has a face (which is a three-dimensional cone) spanned only by v k 1 , v k 2 , v k 3 . For a fan constructed from a triangulation, every three-dimensional cone is the intersection of precisely two four-dimensional cones. Hence there must be another four-dimensional cone spanned by v k 1 , v k 2 , v k 3 and another vector. The only vector which turns As a direct consequence of the above, K is a refinement of J (the converse is wrong, i.e. K = J in general) for any triangulation. This means we do not need to consider the index set L = K ∩ J since it equals K. Also K is ultimately independent of the triangulation chosen for • .

B Toric π -twisted Representations
In this appendix, we collect some results on the representation theory of finite groups acting linearly on the homogeneous coordinates of toric varieties. Due to the decomposition (A.36) for the toric symmetry groups G A , our representations R : → G A can similarly be split into a permutation representation π : → P A and a map r : → H A /G into the continuous part of the ambient space symmetry group. Since this continuous part has the block structure H A /G = K ∈K Gl(K , C)/G the map r can be split up into the corresponding blocks which we denote by r i . The semi-direct product structure of the ambient space symmetry group G A means that r is a π -representation, that is, it satisfies r (γγ ) = π(γ ) −1 r (γ )π(γ )r (γ We will discuss the definition of group cohomology below. The extension sequence (B.4) can be defined by fixing a so-called factor set e : × → K via the multiplication rule By Eq. (B.6), the factor sets are then related by These relations define the group cohomology. Two factor sets f , f are in the same class gives us a linear representation which becomes the projective representationr in Gl(d, C)/C * . The question is now if we can 'undo' the factor set f by enlarging the group toˆ and studying linear representations ofˆ . As discussed in the last section, if we can find a map ϕ ∈ Hom(K , C * ) (that is, ϕ can be used to map e(γ ,γ ) → f (γ ,γ ) and we can think of ϕ(k) as a being proportional to the unit matrix) and h(γ ) such that holds for every factor set f (γ ,γ ) our choice of e(γ ,γ ) is sufficiently general. The central statement is now that the Schur cover obtained by using the Schur multiplier is precisely such that this holds. For the choice K = H 2 ( , C * ), δ is indeed surjective, and any projective representation can be lifted to a linear one. Hence, using the Schur cover we can lift every projective representation to a linear one.

B.1.2 Multiple Schur covers
Let us now discuss the next more complicated cases, which are representations on products of projective spaces First of all, if we are given such a representation and a lift l(γ ), (B.6) implies that f is a map f : × → (C * ) n . We can think of f as a matrix that is, it is block diagonal and each block is proportional to the unit matrix. We can of course rewrite Clearly, equivalent f (γ ,γ ) are again related by (B.8). This equivalence is derived using the action of C * on the homogeneous coordinates and we may consider each of the n actions in turn and denote the corresponding map h i : → C * h i (γ ). This gives rise to a similar decomposition into h(γ ) = h i (γ ) with h i : → C * . With this we find that f (γ ,γ ), f (γ ,γ ) define the same class in H 2 ( , (C * ) n ) if all of the f i and f i define the same class in H 2 ( , C * ) We now use the same Schur cover as before, (that is, we have the same exact sequence (B.4) and e : × → K ), and consider maps ϕ, which now maps ϕ : K → (C * ) n . As is f (γ ,γ ), ϕ(k) will be a block-diagonal matrix with each block proportional to the unit matrix. We can hence also decompose ϕ into homomorphisms ϕ i : K → C * , We are now ready to confront (B.11) for this setup. As we have decomposed ϕ, h, f in the same way, a solution to the generalization of (B.11) can be found if we can solve The existence of such h i , ϕ i for each f (γ ,γ ) is already implied if e(γ ,γ ) originates from the standard Schur cover, so that we are done and the same construction applies here. The linear representation is found aŝ The case of a toric variety may now be discussed in a similar fashion. The difference to having independent products of projective linear groups as in Eq. (B.12) is that we may have more factors of Gl(d i , C) than C * actions in (B.12). Alternatively, we may think of the C * actions on each of the Gl(d i , C) 'blocks' as being correlated. Let us assume we have in total b blocks, with c blocks correlated to the other (b − c) blocks, so that G ∼ = (C * ) b−c . We can find all multi-projective representations in this case by treating all of the blocks as uncorrelated and then checking two more conditions to remove the invalid ones. A sufficient condition for a linear representation to descend to a multi-projective representation in the case of correlated blocks is: This reason for this is that given the same γ , varying k does not change the identification ofr (γ ) insider (k, γ ). Denoting the projection fromˆ to byπ we have thatπ(k, γ ) = π(k, γ ) for all k,k ∈ K . Hence by the definition ofπ we haver (k, γ )r (k, γ ) −1 ∈ G. Finally, we need to identify under which circumstances two different linear representations of give rise to the same multi projective representation. Two projective representationsr 1 andr 2 are equivalent only when there exists a automorphism P of the space of homogeneous coordinates and a function θ : → G, such that for each γ ∈ , we have P(r 1 (γ )) = θ(γ )(r 2 (γ )P). Inserting the definition ofr =π(r (k, γ )), this means that P(π(r (k 1 , γ )) = θ(γ )(π(r 2 (k 2 , γ ))P).

(B.22)
We see that choosing different k 1 and k 2 is equivalent to inserting another factor into θ , namely θ is k 1 , k 2 dependent. However in the sense of equivalent representation, we could just make one convenient choice of k 1 , k 2 for each γ , checking the formula (B.22). This is sufficient to remove the redundancies.

B.2 Twisted representations.
Let us now discuss the case of a non-trivial representation π and the resulting interplay between π and r . Before considering the full problem of finding all (multi-)projective π -twisted representations, we discuss how we can manage to find all π -twisted linear representations r : → P A K ∈K Gl(K , C). To start the construction, let us assume that we have found a, not necessarily injective, group homomorphism π : → P A . (B.23) How can we find all compatible π -twisted representations r satisfying (B.1)? Furthermore, we assume for simplicity that π( ) acts transitively on the 'blocks' K i ∈ K (of course all of these have to have the same size) and leaves everything else unchanged. We will discuss the general case below. Labelling the relevant blocks by K i , where i = 1, . . . , b, we can single out the block K i and consider its stabilizer groups i ⊂ under the action of π . We denote by r i : i → G L(B i , C)/C * the map r restricted to the block K i . Equivalently, we can write the entire representation R as the matrix product R(γ ) = π(γ ) · diag(r 1 (γ ), · · · , r n (γ ). (B.24) Note that the restrictionr i and r i to the stabiliser group i defines a group homomorphism (since, by definition, π acts trivially on the stabiliser). We now explain that it is sufficient to fix a single one of the representationsr i to recover the whole π -twisted representation r . For definiteness, let us hence consider 1 , the stabilizer of the first block, and fix the representationr 1 : 1 → Gl(K 1 , C). To reconstruct the whole action of we pick a set of group elements γ i such that π(γ i )(1) = i. For any block i, we can then write an arbitrary group element γ ∈ uniquely as γ = γ π(γ )(i) hγ −1 i , (B.25) where h ∈ 1 . We can hence think of h to depend on γ and i. To see this, note that we can choose h as h = γ −1 π(γ )(i) γ γ i .

(B.26)
This is indeed an element of 1 since π(h)(1) = π(γ π(γ )(i) ) −1 π(γ )π(γ i )(1) = π(γ π(γ )(i) ) −1 π(γ )(i) = 1. (B.27) We may hence write for any i and γ and using the appropriate h ∈ 1 . Going through the definitions this means that we have This allows us to recover all of the matrices r i (γ ) and hence the entire π -twisted representation fromr 1 . The group homomorphismsr 1 can, in turn, be constructed from the linear representations of the Schur coverˆ 1 , as discussed above. In general π : → P A does not act transitively on the blocks and there may be several orbits. In this case, we can apply the above method separately for each of the orbits and then combine all of the data to find the representation R : → P A H A .
In summary, in order to find the representations R = (π, r ) : → P A H A , we have to carry out the following steps.
(2) Find all orbits of the blocks under the action of π .
(3) For each orbit {K i }, pick out the block K 1 and determine its stabiliser group, 1 , under the action of π . (4) Choose the group elements γ i such that π(γ i )(1) = i. (5) Study all representationsr 1 : 1 → Gl(K 1 , C)/C * by considering the Schur cover, 1 of 1 . (6) Re-construct the entire map r fromr 1 by using Eq. (B.26). (7) Repeat this process for each orbit, assemble the results into the complete representation R and check if R is faithful.

C Patches on Toric Varieties and Smoothness of Hypersurfaces
Let us describe how the coordinate patches corresponding to each four-dimensional cone σ ∈ of a toric variety with fan are found. Let us start with cones corresponding to smooth affine varieties, in which case the corresponding patch is simply C 4 . In this case, the rays of σ are generated by four lattice vectors v σ i such that the matrix (v σ 1 , v σ 2 , v σ 3 , v σ 4 ) has determinant one and the v i generate the N-lattice with integer coefficients. The four rays of the dual cone 6 σ ∨ are then generated by four lattice vectors v ∨ i which generate the M-lattice and satisfy v ∨ j , v σ i = δ i j . (C.1) The affine toric variety V (σ ) = (C * ) 4 associated with each smooth four-dimensional cone has a stratification V (σ ) = (C * ) 4 i=1..4 (C * ) 3 j=1..6 (C * ) 2 k=1..4 (C * ) pt (C. 2) in which each stratum (C * ) n is associated with an 4 − n dimensional cone contained in σ . We can find good coordinates for all of these strata using points on the M-lattice. In the present case, the v ∨ i defined above give rise to coordinates on (C * ) 4 which have a well-defined limit to the other strata in (C.2), that is, we can use them as coordinates for the whole of V (σ ). These coordinates can be written aŝ where we denote the homogeneous coordinates corresponding to the v σ i by x σ i and keep the notation of denoting all homogeneous coordinates simply by x i . In the above expression, the first product ranges over all homogeneous coordinates (ray generators of the fan ) and the second one only over homogeneous coordinates associated with rays not contained in σ .
On the patch V (σ ), we can rewrite the defining polynomial p of a TCY hypersurface X in terms of the local coordinatesx i . Every monomial p m in p is associated with a vector m in the M-lattice and we may write m = i=1.. 4 Not surprisingly, we can rewrite every monomial purely in terms of thex i such that p(x σ i ) = p(x σ i ) when we set all other coordinates = 1. Let us take p σ to be the polynomial which is obtained when we set {x j }\{x σ i } to 1. This is of course the same as 'gauge fixing' all coordinates x i except the local coordinates x σ i . For each patch V (σ ) = C 4 we can define the Jacobi ideal of the local defining polynomial p σ by To check smoothness on V (σ ) we have to verify that the dimension of this ideal is −1.
Let us now discuss the singular patches V (σ ). Naively, singular patches arise because gauge fixing all coordinates except the ones in σ may leave some residual discrete group G ⊂ G which still has to be divided out. Hence, the local patch has the structure V (σ ) = C 4 /G and the local coordinates x σ i really parametrize the covering space C 4 . As before, we may find coordinates on (C * ) 4 ⊂ V (σ ) by considering four vectors generating the M-lattice. What we cannot hope for, however, is for the ray generators of σ ∨ i to generate the M-lattice (over Z). We can however extend these to a set of generators of the M-lattice by considering a larger set {v ∨ μ } ⊃ {v ∨ i } . These can be used to define a number of coordinates on V (σ )x As we can expand any vector of the m-lattice in terms of the v ∨ μ with integer coefficients, we can again rewrite p as a function solely of the v ∨ μ upon dividing out {x j }\{x σ i } x j . (This also follows more abstractly from the fact that C[v ∨ μ ] is the polynomial ring of V (σ ).) However, we have more than four v ∨ μ , so that there are a number of linear relations μ α μ v ∨ μ = 0 ( C . 7 ) with α μ ∈ Z. These turn into relations μx α μ μ = 1. (C.8) Hence we now find a complete intersection as the appropriate expression for p σ . The additional defining equations capture the orbifold singularities on V (σ ) = C 4 /G. As a simple example, let us consider a Calabi-Yau (elliptic curve) hypersurface of P 123 . We take homogeneous coordinates (y, x, z) ∼ (λ 3 y, λ 2 x, λz). The fan of P 123 has rays generated by The hypersurface can be described by y 2 = x 3 + f xz 4 + gz 6 (C.10) Let us discuss the cone σ spanned by x and z. Gauge fixing y = 1 leaves a residual Z 3 acting on x and z so that we expect V (σ ) = C 2 /Z 3 . The generators of σ ∨ are −1 0 →ẑ = z 3 /y −2 3 →x = x 3 /y 2 (C.11) Clearly, we cannot rewrite (C.10) in terms of these coordinates alone, which corresponds to the fact that the above vectors do not generate the M-lattice by only a sublattice M such that M/M = Z 3 . We can cure this by defining a new coordinate 1 −1 →ξ ≡ xz/y (C.12) so that (C.10) becomes (after dividing by y 2 ) 1 =x + fξẑ + gẑ 2 = 0 (C.13) However, we now have the additional relation ξ 3 =ẑx (C.14) which is nothing but the defining equation of an A 3 singularity (which is the same as C 2 /Z 3 ) embedded into C 3 . As we have already discussed, we have to make sure our hypersurface X defined by p = 0 must not meet the singularities of the ambient space as this would lead to singularities on X as well. As the singularities arise through fixed points of the group action of G on C 4 , we can lift X ⊂ V (σ ) to a (different) hypersurfaceX in C 4 which descends to X via a free group action of G. Smoothness of the covering space is then equivalent to smoothness of the quotient. We are hence interested in studying coordinates on the cover of V (σ ). These coordinates can be found by writing thex μ as monomials such that (C.8) is automatically solved. This, however, can be done by simply using (C.6) and setting all {x j }\{x σ i } to unity. We have hence found that we can treat all fourdimensional cones in the very same way provided that the orbifold points of the ambient space A do not meet the invariant hypersurface X .