Branes wrapped on orbifolds and their gravitational blocks

We construct new supersymmetric $\mathrm{AdS}_2\times \mathbb{M}_4$ solutions of $D=6$ gauged supergravity, where $\mathbb{M}_4$ are certain four-dimensional orbifolds. After uplifting to massive type IIA supergravity these correspond to the near-horizon limit of a system of $N$ D4-branes and $N_f$ D8-branes wrapped on $\mathbb{M}_4$. In one class of solutions $\mathbb{M}_4 = \Sigma_{\mathrm{g}}\ltimes\Sigma$ is a spindle fibred over a smooth Riemann surface of genus $\mathrm{g}>1$, while in another class $\mathbb{M}_4 = \Sigma\ltimes\Sigma$ is a spindle fibred over another spindle. Both classes can be thought of as orbifold generalizations of Hirzebruch surfaces and, in the second case, we describe the solutions in terms of toric geometry. We show that the entropy associated with these solutions is reproduced by extremizing an entropy function obtained by gluing gravitational blocks, using a general recipe for orbifolds that we propose. We also discuss how our prescription can be used to define an off-shell central charge whose extremization reproduces the gravitational central charge of analogous $\mathrm{AdS}_3\times \mathbb{M}_4$ solutions of $D=7$ gauged supergravity, arising from wrapping M5-branes on $\mathbb{M}_4$.


Introduction
For two decades several examples of the AdS/CFT correspondence have been engineered considering setups involving branes wrapped on compact cycles, preserving supersymmetry via a topological twist on the brane world-volumes [1]. On the supergravity side, these give rise to supersymmetric solutions involving a warped product of an anti-de Sitter (AdS) spacetime, together with an internal compact manifold. The corresponding dual field theories are obtained from suitably twisted compactifications of the parent theories on smooth Riemann surfaces embedded holomorphically in Calabi-Yau manifolds, or different higher-dimensional calibrated submanifolds in special holonomy manifolds.
Solutions with an AdS 2 factor are particularly interesting because on general grounds these are expected to correspond to the near-horizon limit of supersymmetric black holes arising in string or M-theory. In this case, the AdS/CFT correspondence gives access to the microscopic entropy of these black holes, obtained via counting the microstates of the dual supersymmetric quantum mechanics, which becomes particularly useful when the latter arises from a twisted compactification of a higher-dimensional gauge theory. The paradigmatic example of this successful strategy is that of the four-dimensional static black holes of [2], for which a microscopic account of the entropy has been obtained in [3].
The construction presented in [4] changed somewhat the rules of the game, showing that one can obtain examples of holographic dual pairs with genuine novel features by considering D3-branes wrapped on the orbifold Σ = WCP 1 [n − ,n + ] , where n − , n + are coprime integers, sometimes referred to as spindle. Perhaps the most striking aspect of this construction is that supersymmetry is preserved via a novel mechanism different from the topological twist, that was later called anti-twist. Further examples of various branes wrapping the spindle have been discussed in [5][6][7][8][9][10][11][12][13][14][15][16] and closely related constructions appeared in [17][18][19][20][21][22][23][24]. In general, the only two possible ways to preserve supersymmetry on a spindle are by the anti-twist of [4] or the twist introduced in [9], which is a global version of the standard topological twist [12].
Once it has been established that it makes sense to wrap branes on spaces with orbifold singularities, it is natural to extend the exploration of the landscape of holographic dualities to branes wrapping orbifolds of dimension higher than two. Simple examples of branes wrapped on four-dimensional orbifolds of the type M 4 = Σ g × Σ were discussed in [7] for M5-branes and in [11,13] for D4-branes (see also [25]). In this paper we will focus mainly on D4-branes, extending the results of [11] in various ways. In particular, we will present two classes of new supersymmetric solutions of gauged D = 6 supergravity of the type AdS 2 × M 4 , where M 4 are four-dimensional orbifolds. After uplifting to massive type IIA supergravity these describe a system of N D4-branes and N f D8-branes wrapped on M 4 , and are expected to correspond to the near-horizon limit of novel supersymmetric AdS 6 black holes. Previously, supersymmetric AdS 2 × M 4 solutions of the same theory, where M 4 is either a negatively curved Kähler-Einstein manifold or the product of two Riemann surfaces Σ g 1 × Σ g 2 (with g 1 , g 2 > 1), were discussed in [26] and [27], with the latter reference also presenting numerical solutions for supersymmetric AdS 6 black holes with Σ g 1 ×Σ g 2 horizons. Below we will construct two new classes of solutions. In one class M 4 = Σ g ⋉ Σ is a spindle fibred over a smooth Riemann surface of genus g > 1, while in the other class M 4 = Σ ⋉ Σ is a spindle fibred over another spindle. Both classes can be thought of as generalizations of Hirzebruch surfaces and analogous AdS 3 × M 4 solutions of D = 7 gauged supergravity, arising from wrapping M5-branes on M 4 , appeared in [28].
One of the motivations for this paper was to begin investigating entropy functions whose extremization reproduces the entropy associated with the corresponding supergravity solutions, comprising orbifolds of dimension higher than two. In the two-dimensional case, entropy functions associated with compactifications on the spindle were discussed in [6,9,29] and a conjectural extension to more general off-shell free energies was put forward in [11]. In particular, for all cases, precise formulas for the off-shell free energies and for the associated constraints obeyed by the magnetic fluxes and fugacities were given in terms of data at the north and south poles of the spindle. These are the fixed points of the Killing vector generating its U(1) azimuthal symmetry and the contributions to the off-shell free energies arising at these fixed points are referred to as gravitational blocks [30]. Here we will extend the conjecture of [11] presenting an orbifold entropy function whose contributions (blocks) arise from fixed points of a U(1) 2 action on M 4 , together with a recipe for determining the appropriate constraints.
The key for formulating this proposal is to provide a toric description of our solutions. Despite the fact that our orbifolds are not symplectic and therefore there are no moment maps with associated polytope, assuming a U(1) 2 action, we can still define a set of toric data from the co-dimension two loci where one U(1) ⊂ U(1) 2 degenerates. In particular, these can be determined by studying the Killing vectors of the explicit metrics on M 4 . We can then associate a compact polytope defined using these data as normal vectors to the edges and identify the vertices as fixed points of the torus action, as in symplectic toric geometry. Our proposal for the entropy function reproduces the results of the explicit supergravity solutions presented here, but we formulate it in a way that is in principle applicable to more generic toric orbifolds, including toric manifolds, for which explicit solutions have not yet been found.
Finally, we discuss how our prescription can be extended to describe analogous supersymmetric AdS 3 × M 4 solutions of D = 7 gauged supergravity arising from wrapping M5-branes on M 4 . In these cases one can employ the method of integrating the M5-branes anomaly polynomial on M 4 to obtain an off-shell central charge of the corresponding d = 2, N = (0, 2) SCFTs and our proposal reproduces the results available in the literature. In particular, when M 4 is actually a toric manifold, our recipe reduces to the construction of [31], whereas for the specific case of M 4 = Σ ⋉ Σ we recover the results of [28].
The rest of the paper is organized as follows. In section 2 we present the new solutions in D = 6 and demonstrate their supersymmetry by solving the relevant Killing spinor equations. In section 3 we uplift our solutions to massive type IIA supergravity and calculate their associated geometric entropies. In section 4 we work out the toric data of the relevant solutions. In section 5 we present our conjectural entropy functions and use these to recover the entropies of the above solutions, as well as the central charge associated with the solutions in [28]. Appendix A contains a discussion about the different quantization conditions on the fluxes and their resolution. In appendix B we present a special limit in which a Riemann surface can be retrieved from a spindle. In appendix C we derive toric data from different metrics, including the solutions of [28].
Note added: as this work was being finalized we became aware that there would be significant overlap with the results of [32], which appeared on the arXiv on the same day.

Solutions in D = 6 gauged supergravity
In this paper we discuss solutions of a D = 6 gauged supergravity with gauge group U(1) 2 , comprising two gauge fields A 1 , A 2 , a two-form B and two real scalar fields ϕ = (ϕ 1 , ϕ 2 ). This model can also be obtained as a sub-sector of an extension of Romans F (4) gauged supergravity [33], coupled to three vector multiplets [34,35]. The bosonic part of the action reads 1 where F i = dA i , H = dB and the scalar fields ϕ are parameterized as The scalar potential is with g the gauge coupling and m the mass parameter, and where for later convenience we defined X 0 = (X 1 X 2 ) −3/2 . A solution to the equations of motion of the model is supersymmetric if and only if it satisfies also the following set of Killing spinor equations [34]: These follow from setting to zero the supersymmetry variations of the fermionic fields of the theory with three vector multiplets [34], that do not vanish automatically in the sub-truncation that we are considering.
The SU(2) indices A, B are raised and lowered as ǫ A = ε AB ǫ B and ǫ A = ǫ B ε BA , where ε AB = −ε BA and its inverse matrix ε AB is defined such that ε AB ε AC = δ B C . The supersymmetry parameter ǫ A is an eight-component symplectic-Majorana spinor, hence it satisfies the condition where B 6 is related to the six-dimensional charge conjugation matrix C 6 by B 6 = −i C 6 Γ 0 . 1 Here and in what follows we define, for any p-form ω, |ω| 2 = 1 p! ω µ1...µp ω µ1...µp .

AdS
The first local solutions to the supergravity action (2.1) that we consider are of the type AdS 2 × Σ g ⋉ Σ 2 , which consist of a two-dimensional spindle Σ 2 non-trivially fibred over a Riemann surface Σ g , with genus g > 1. When g = 0 a solution to the equations of motion still exists, but it is not supersymmetric. Specifically, the backgrounds are described by the following set of fields solving the equations of motion where ds 2 AdS 2 denotes the unit radius metric on AdS 2 , ds 2 Σg the unit radius metric on Σ g and ω Σg is such that dω Σg = vol(Σ g ). The functions h i and F are given by with q 1 , q 2 two real parameters. A curvature singularity lies at y = 0, hence without loss of generality, in what follows we will restrict to y > 0. The first step we take in the analysis of these solutions is to prove that they are supersymmetric by constructing the local form of the Killing spinors solving the equations (2.4)-(2.6). We employ the following orthonormal frame whereêâ,â = 0, 1, is the zweibein on AdS 2 , whose coordinates are denoted as xμ, anděǎ, a = 2, 3, is the zweibein on Σ g , whose coordinates are denoted as xμ. Equation (2.4) then splits into the following system where we definedh ≡ log(h 1 h 2 ) and we applied the gauge transformation A i → A i +α i dz. Equations (2.5) and (2.6) yield the same constraint We note that the equations along y and z in (2.12) and equation (2.13) are the same as in the AdS 4 × Σ system (cf. equations (3.3) and (3.4) of [11]), a fact that will play an important role later on. From the fourth equation in (2.12) we see immediately that setting 14) leads to Killing spinors independent of z and we will adopt this choice in the reminder of this section. Before tackling the Killing spinor equations, we note the following fact about the AdS 4 ×Σ solution of [11]: choosing a suitable decomposition of the six-dimensional gamma matrices, it is possible to write the six-dimensional Killing spinors as tensor products of a single AdS 4 spinor and a single spinor on Σ. Specifically, we consider the decomposition where γã,ã = 0, . . . , 3, are the (Lorentzian) gamma matrices in D = 4, ρî,î = 1, 2, are the (Euclidean) gamma matrices in D = 2 and ρ * = −i ρ 1 ρ 2 is the related chiral matrix. For ρî we choose the following representation and we take B ρ 2 = −σ 2 . For consistency, the six-dimensional matrix B 6 decomposes as where γ 5 = i γ 0 γ 1 γ 2 γ 3 is the four-dimensional chiral matrix. Adopting this decomposition, the symplectic-Majorana condition (2.8) and the Killing spinor equations (3.3)-(3.4) of [11] can be solved to give where ϑ (4) = ϑ (4) (xμ) is a Majorana Killing spinor on AdS 4 and Here, ξ η is a complex constant and we have defined which satisfy F (y) = f 1 (y) f 2 (y). Notice that η 1 and η 2 are related by η 2 = i σ 1 (η 1 ) * , which descends from the symplectic-Majorana condition (2.8).
Going back to the AdS 2 × Σ g ⋉ Σ 2 background, in order to solve its Killing spinor equations we adopt the same decomposition of the six-dimensional gamma matrices as in (2.15). Moreover, we decompose the four-dimensional gamma matrices γã as where βâ are the (Lorentzian) gamma matrices in D = 2, τî,î = 1, 2, are the (Euclidean) gamma matrices in D = 2 and τ * = −i τ 1 τ 2 is the related chiral matrix. Explicitly, we adopt the same representation as in (2.16) and we take B τ 2 = σ 1 . The four-dimensional matrices B 4 and γ 5 decompose as where β * = −β 0 β 1 is the (Lorentzian) chiral matrix in D = 2. The ansatz for the sixdimensional Killing spinors is where ϑ = ϑ(xμ) is a Majorana Killing spinor on AdS 2 , hence∇μϑ = 1 2 βμϑ and ϑ * = B β 2 ϑ. χ A = χ A (xμ) are two two-component spinors defined on the Riemann surface Σ g , while η 1 and η 2 are given in (2.19). Putting all the ingredients together, the symplectic-Majorana condition (2.8) and the Killing spinor equations (2.12)-(2.13) can be solved to give where ξ η is a complex constant. Notice that we can express χ 2 as We can now count the number of supersymmetries preserved by our AdS 2 × Σ g ⋉ Σ 2 solution. ϑ is a Majorana spinor, hence it has two real degrees of freedom, while the tensor product χ A ⊗ η A is fully determined by the product ξ χ ξ η , which has two real degrees of freedom. Therefore, there are four real independent Killing spinors, that is one quarter of the number of supersymmetries of the six-dimensional N = (1, 1) theory, hence the solution is 1/4-BPS.
We now proceed to the global analysis of the local solution (2.9), for which we set m = 2g/3 without loss of generality 2 . Fixed a generic point on the Riemann surface Σ g , the metric on the spindle Σ 2 reads In order to have a well-defined metric and positive scalars X i we need to take F > 0, h 1 > 0, h 2 > 0 in a closed interval not containing the curvature singularity in y = 0, thus without loss of generality we restrict to y > 0. Moreover, denoting [y N , y S ] the range of the coordinate y, Σ 2 is a proper spindle given that [11] gF ′ (y N ) where n ± are two co-prime integers and ∆z is the periodicity of the z coordinate. These two relations ensure that at the poles y = y N,S there are Z n ∓ orbifold singularities, respectively. An additional constraint comes from the quantization of the magnetic fluxes across Σ 2 , which arises from the requirement that gA i be well-defined connection oneforms on U(1) bundles over Σ 2 . Given any point of Σ g , we must have It can be shown that [11] p 1 + p 2 = n + + n − , (2.29) meaning that the R-symmetry gauge field A R ≡ g(A 1 + A 2 ) is a connection one-form on the line bundle O(n + + n − ) over Σ 2 , and therefore the integers p i can be conveniently parameterized as where z is an appropriate rational number. The value of the sum p 1 + p 2 implies that A R realizes on the bundle the type of twist that was dubbed "twist", as for the D4 [11] and M5-branes [8] wrapped on a spindle. All the conditions we presented were studied in [11], to which we refer for a detailed analysis, and it was proven that they can be satisfied if n − < n + , p 1 < 0, p 2 > 0, the expressions of q i and y N,S are determined as in equations (3.39) and (3.31) of [11] and z has periodicity Here, χ 2 , the Euler characteristic of Σ 2 , and µ are defined as and x is the only solution of the quartic equation lying inside the range 0 < x < 1.
As we already mentioned, the internal space has the structure of a spindle Σ 2 fibred over a Riemann surface Σ g . This fibration is well-defined in the orbifold sense if the one-form η describing the fibration is globally defined, i.e.
This requirement yields a quantization condition relating the spindle data (n ± , z) and the genus g, namely We refer to appendix A for the analysis of this constraint.
The last conditions follow from the quantization of the fluxes F i through two-cycles in Σ g ⋉ Σ 2 . Their integration through Σ 2 were computed in (2.28), hence they are already appropriately quantized. Next, we define the two-cycles S − ≡ {y = y N } and S + ≡ {y = y S }, corresponding to the two poles of Σ 2 , and compute, e.g., where we substituted t as in (2.36). Similarly to what happened in (2.34), the quantization condition requires (2.37) to be an integer. In appendix A we present some examples in which this constraint is satisfied. The quantization through S + can be addressed easily since an explicit computation shows that therefore the related condition is automatically satisfied due to the fact that n ± divide t, as explained in appendix A. We also notice that (2.38) agrees with the homology relation S + − S − = t Σ 2 . The fluxes of F 2 can be obtained flipping the sign of z in the previous formulas and, as a consequence, the quantization of its fluxes automatically holds exchanging p 1 with p 2 . Focussing on the R-symmetry gauge field, its fluxes read and, as expected, they are correctly quantized because n ± divide t.

AdS
A second class of solutions to the model described by action (2.1) is the AdS 2 × Σ 1 ⋉ Σ 2 background, where a two-dimensional spindle Σ 2 is non-trivially fibred over a different two-dimensional spindle Σ 1 . Remarkably, there exists a specific limit of the coordinates and the parameters in which the AdS 2 ×Σ g ⋉Σ 2 solutions are retrieved, both their bosonic content and the Killing spinors. All the details and computations are left to appendix B. Going back to the systems under exam, they read where ds 2 AdS 2 denotes the unit radius metric on AdS 2 and a is a real parameter. The functions h i and F are the same as in (2.10), while q is given by (2.41) Notice that formally taking a = 0 we retrieve the AdS 4 × Σ background studied in [11]. In order to construct the local form of the Killing spinors and to demonstrate that this solution is supersymmetric, we must first specialize the Killing spinor equations (2.4)-(2.6) to our system employing the orthonormal frame whereêâ,â = 0, 1, is the zweibein on AdS 2 , whose coordinates are denoted as xμ. Equation (2.4) then splits into the following system whereh ≡ log(h 1 h 2 ) is defined as before and we employed the gauge transformation A i → A i + α i dz. Equations (2.5) and (2.6) yield, again, the same constraint Also in this case, the Killing spinor equations along y and z and the last constraint are identical to those obtained in the AdS 4 × Σ background. Once again we fix the gauge such that α 1 + α 2 = 2m/g, in order to have a z-independent Killing spinors. The explicit construction of the Killing spinors proceeds in a way similar to the AdS 2 × Σ g ⋉ Σ 2 case. First of all, we employ the following decomposition of the six-dimensional gamma matrices (cf. equations (2.15) and (2.21)): which implies, where βâ are the (Lorentzian) gamma matrices in D = 2, ρî and τî,î = 1, 2, are the (Euclidean) gamma matrices in D = 2 and ρ * and τ * are the related chiral matrices. In the Euclidean sectors we adopt the representation (2.16) and we take B ρ 2 = −σ 2 and B τ 2 = σ 1 . The six-dimensional Killing spinors are assumed to consist of the tensor product of spinors living on AdS 2 , Σ 1 and Σ 2 . Namely, where ξ χ is a complex constant and which satisfy q(x) = Q 1 (x) Q 2 (x). Notice that the relation χ 2 = −i σ 2 (χ 1 ) * holds also in this case. The count of the real degrees of freedom is identical to the AdS 2 × Σ g ⋉ Σ 2 case and shows that the solution is 1/4-BPS, since it preserves four real supercharges. The first part of the global analysis of the AdS 2 ×Σ 1 ⋉Σ 2 solution (2.40) proceeds along the path traced in the previous section. Chosen a generic point on the base spindle Σ 1 , the metric of Σ 2 is the same as in (2.26), therefore are identical also the conditions for it to describe a proper spindle with Z n ∓ orbifold singularities. In particular, ∆z must be as in (2.31). Moreover, at a generic point of Σ 1 (different from the north and south poles) the magnetic fluxes are quantized as in (2.28): with p 1 + p 2 = n + + n − . The base spindle Σ 1 , described by the metric was thoroughly studied in [5], thus we will take advantage of the results presented therein, setting j = 0. Σ 1 is indeed a spindle, characterized by the two co-prime integers m ± , with m − < m + 3 , if the parameter a and the periodicity of the coordinate ψ are The coordinate x is restricted between the two middle roots of the quartic polynomial q(x), At the north and south poles of Σ 1 , namely x = x N,S , are present Z m ∓ orbifold singularities, respectively. Moreover, defined the vector field its magnetic flux through the spindle Σ 1 is In the four-dimensional theory this implies that 2A 4d is a connection one-form on the line bundle O(m + − m − ) over Σ 1 and, as a consequence, that we have anti-twist.
Similarly to the Σ g ⋉ Σ 2 system, the fibration of Σ 2 over Σ 1 is well-defined given that where, in this case, We derive the following constraint on the quantum numbers n ± , m ± and z which is solved in appendix A. Like we did in the previous section, we study the quantization of the fluxes of F i , now through two-cycles in Σ 1 ⋉ Σ 2 . The magnetic fluxes across Σ 2 were computed in (2.28) and are, thus, already quantized. Again we define the two-cycles S − ≡ {y = y N } and S + ≡ {y = y S }, which corresponds to two copies of the base spindle Σ 1 , and compute (2.60) In analogy to (2.37), the flux (2.60) must be an integer divided by (m + m − ), and we refer to appendix A for the analysis of this constraint. The same quantization condition, but with F 1 integrated across S + , is automatically satisfied due to the relation which can be both computed explicitly and obtained from the homology relation On the other hand, the fluxes of F 2 can be obtained switching the sign of z, i.e. exchanging p 1 with p 2 , in the previous formulas and, as a consequence, the quantization follows directly. The fluxes of the R-symmetry gauge potential read and are well quantized, i.e. they are integers divided by m + m − , because n ± divide t.

Truncation to D = 4 minimal gauged supergravity
An accurate inspection of the AdS 2 × Σ g ⋉ Σ 2 and AdS 2 × Σ 1 ⋉ Σ 2 solutions given in (2.9) and (2.40), respectively, hints at the existence of an underlying consistent truncation of matter-coupled D = 6, U(1) 2 gauged supergravity on a spindle Σ 2 , down to D = 4, N = 2 minimal gauged supergravity. Even though we do not demonstrate the consistency of this truncation 4 , we conjecture its validity in analogy with the results of [28] and we will provide supporting evidence. For completeness, we report here the action describing the four-dimensional theory along with the corresponding Killing spinor equations where F = dA , {γ a , γ b } = 2η ab and the AdS 4 radius has been set to 1. Given a solution of D = 4 minimal gauged supergravity comprising a metric ds 2 (4) and a gauge potential A , we propose the following ansatz for the six-dimensional solution: Here, ⋆ 4 is the Hodge dual with respect to the metric ds 2 (4) and (cf. (2.10)) This ansatz may be justified comparing the two six-dimensional solutions already mentioned with the following AdS 2 × Σ g [36] and AdS 2 × Σ 1 [5] backgrounds In both cases, plugging the four-dimensional solutions into the ansatz (2.65) gives the corresponding ones in D = 6. An additional example is given by the AdS 4 × Σ 2 solution of [11], which can be obtained starting from the AdS 4 vacuum.
In addition to the truncation of the bosonic sector, we also propose an ansatz connecting the Killing spinors of the two supergravity theories in D = 4 and D = 6. This conjecture is driven both by the construction presented in [28] and by the fact, already observed, that some of the Killing spinor equations are the same for the two background considered in this paper, as well as for the AdS 4 × Σ 2 solution.
The truncation ansatz for the Killing spinors heavily depends on the choice of the decomposition of the six-dimensional gamma matrices. In what follows we will consider (2.15) Γã = γã ⊗ σ 3 , where γã,ã = 0, . . . , 3, are the (Lorentzian) gamma matrices in D = 4. For consistency, B 6 decomposes as where γ 5 = i γ 0 γ 1 γ 2 γ 3 is the four-dimensional chiral matrix and B 2 = −σ 2 . We conjecture that the first Killing spinor ǫ 1 takes the expression where ζ is a Killing spinor of D = 4 minimal gauged supergravity and η is given by 5 where ξ is a complex constant and satisfying F (y) = f 1 (y) f 2 (y). The second Killing spinor ǫ 2 is determined by the symplectic-Majorana condition (2.8), explicitly where we defined the charge conjugate spinor ζ c = B −1 4 ζ * . As for the bosonic part of the truncation, we can justify this ansatz making contact with the Killing spinors of D = 4 minimal gauged supergravity. Assuming, additionally, the decomposition of the four-dimensional gamma matrices in (2.21), the Killing spinors of the AdS 2 × Σ g [36] and AdS 2 × Σ 1 [9] (see also [5]) backgrounds are respectively 6 (2.76) 5 Here we relabelled η 1 with η and ξ η with ξ with respect to (2.19). 6 Here ζ is a superposition of the two Killing spinors ǫ 1 and ǫ 2 of [9].
where ϑ is a Majorana Killing spinor on AdS 2 and we defined . Uplifting these two spinors by means of the proposed recipe we obtain exactly the Killing spinors of the corresponding six-dimensional solution. For what concerns the AdS 4 × Σ 2 background, we must take ζ = γ 5 ϑ (4) , with ϑ (4) a Killing spinor on AdS 4 , due to a slightly different convention on the four-dimensional Killing spinor equations. The existence and the consistency of the truncation for the Killing spinors would prove that any supersymmetric solution of D = 4 minimal gauged supergravity gives rise to a supersymmetric solution of D = 6, U(1) 2 gauged supergravity after uplifting on Σ 2 via (2.65).

Uplift to massive type IIA supergravity
The local solutions presented in section 2 can be embedded in massive type IIA supergravity by means of an appropriate ansatz. A consistent truncation for the field content of interest is unfortunately not available. However, if we set m = 2g/3 we can consider an enlarged version of the ansatz presented in [11]. This proposal 7 comes as a superposition of the formulas presented in [37] and [38], which we recover in the respective limits. In the former only one scalar is considered, but the six-dimensional two-form B is non-zero, while the latter encompasses four scalar fields, but vanishing B field. The ansatz for the metric in the string frame and the dilaton is the same as in [11] but now we have also a non-vanishing two-form field Here ds 2 6 stands for the line element of the six-dimensional space M 6 and the one-forms σ i ≡ dφ i − gA i are built up from the six-dimensional gauge fields. The angular coordinates φ 1 , φ 2 have canonical 2π periodicities, whereas the coordinates µ a , with a = 0, 1, 2, satisfy the constraint µ 2 a = 1. Lastly, the warp factor reads A possible parameterization of the coordinates µ a is µ 0 = sin ξ , µ 1 = cos ξ sin η , µ 2 = cos ξ cos η , (3.5) with η ∈ [0, π/2] and ξ ∈ (0, π/2], where the range of ξ is fixed by the necessity of having µ 0 > 0. In this way the ten-dimensional metric (3.1) parameterizes, at each point of M 6 , a four-dimensional hemisphere, denoted as S 4 from now on. The metric on S 4 is in general squashed, and becomes the metric on "half the round four-sphere" when X 1 = X 2 = 1. The RR sector of massive type IIA supergravity comprises the ten-dimensional Romans mass and the two-form flux along with the four-form flux, conveniently written in terms of its Hodge dual as 8 where we defined The Hodge star operators ⋆ 10 and ⋆ 6 are computed using the string frame metric (3.1) and ds 2 6 , respectively. Even though we did not prove that this truncation is consistent, we tested it successfully against our AdS 2 × Σ g ⋉ Σ 2 (2.9) and Provided the equations of motion of the six-dimensional gauged supergravity are satisfied, the above ansatz should solve the equations of motion of massive type IIA supergravity, which can be derived from the string frame action 9 , (3.9) where, in these conventions, the field strengths are (3.10) As noticed in [11], the ten-dimensional equations of motion are invariant under a scaling symmetry, which defines an "improved uplift" and whose action is where λ is any strictly positive constant. This additional parameter λ enters in the regularity analysis of the uplifted ten-dimensional solutions. In particular, its presence will be necessary for the fluxes to be correctly quantized [11].
In the next subsections we will specify our proposed truncation ansatz to the solutions presented in (2.9) and (2.40). After having applied the scaling symmetry (3.11), we will quantize the fluxes, thus ensuring the global regularity of the ten-dimensional solutions.

AdS
We begin this section by presenting the key ingredients for the flux quantization and the computation of the entropy. The quantization conditions are where M 4 = Σ g ⋉ Σ 2 . At ξ = 0, where µ 0 = 0 and the warp factor is singular, there are an O8-plane and N f = 8 − n 0 coincident D8-branes, while N D4-branes are wrapped over Σ g ⋉ Σ 2 . These integers get contributions from the following fluxes while the metric and the dilaton enter in the computation of the entropy For convenience, we defined the functions Additionally, we report for completeness also the other ten-dimensional fields, which however do not play an important role in our analysis Plugging our solution in the first two conditions of (3.12) gives It is worth noting that, being N and n 0 integers, the second equation would be inconsistent with λ = 1. Such a problem arises from the fact that, without introducing λ, the first two constraints in (3.12) would have to be imposed on an unique dimensionless parameter, namely (gℓ s ). This makes the scaling symmetry (3.11) crucial to establish the regularity of the ten-dimensional uplifted solution. The integration in the third condition of (3.12) is performed along a representative of M 4 , which we take to be at the pole of the hemisphere S 4 , i.e. in ξ = π/2. The result can be implicitly expressed in terms of the integer parameters (g, n ± , z, N) by means of equations (3.31) and (3.38) of [11], and reads The analysis of this quantization condition is left to appendix A. We now move to the computation of the entropy. The first step is to write the tendimensional metric in the form There is a U(1) 2 symmetry acting on S 4 , with gA i connections on the associated circle bundles, and the correct quantization of the corresponding magnetic fluxes t i and s ± i , as discussed in section 2.1, ensures that the ten-dimensional solution is a well-defined orbifold. The entropy of our solution can be read from the two-dimensional effective Newton constant G (2) as Applying this to metric (3.15), we obtain where F S 3 ×Σ 2 is the free energy of d = 3, N = 2 SCFTs that arise from a system of N D4-branes and N f D8-branes wrapped on a spindle [11] Similarly to the previous section we present the explicit expressions of the fluxes, which are then quantized. Apart from the Romans mass, which is the same as in (3.13), the complete uplifted solution includes the four-flux 25) and the metric and dilaton where ∆ h and U h were defined in (3.17). Again, even if they are not fundamental for the computations, we also have The quantization conditions are the same as in (3.12), where now M 4 = Σ 1 ⋉ Σ 2 , and the first two give again the relations (3.19). The last one, expressed in terms of (m ± , n ± , z, N), reads We refer to appendix A for a detailed analysis. The ten-dimensional geometry is similar to the solution of the previous subsection, with the internal space M 8 having a fibration structure S 4 ֒→ M 8 → M 4 . Again, there is a U(1) 2 symmetry acting on S 4 , with gA i connections on the associated circle bundles. Therefore, the uplifted solution can be considered a well-defined orbifold thanks to the correct quantization of the magnetic fluxes of gA i , as discussed in section 2.2.
The entropy of this solution can be computed following the same steps of the previous case and the final result is where F S 3 ×Σ 2 is given in (3.24) and A Σ 1 is the area of the horizon of the four-dimensional supersymmetric black hole solutions with AdS 2 × Σ 1 near-horizon studied in [5], namely (3.31)

Toric geometry
In this section we will provide a description of the orbifold Σ 1 ⋉ Σ 2 in terms of certain combinatorial data, that will be used in the next section to construct the entropy function. Since Σ 1 ⋉ Σ 2 admits a U(1) 2 action, corresponding to the U(1) 2 isometry of the metric (2.40), it is natural to seek a description in terms of toric geometry. Below we will discuss a method for extracting a set of toric data that is sufficient for the purposes of this paper, leaving a more detailed study of the geometry for future work. As a warm-up, we will begin with an S 2 ⋉ Σ 2 geometry, namely a spindle fibred over a smooth twosphere. This can be obtained formally as an analytic continuation of the Σ g ⋉ Σ 2 family and although it is not an actual supersymmetric background, it solves the equations of motion. This case can be analysed in the framework of symplectic toric geometry and we will derive a labelled polytope [39] from the image of an associated moment map. From the algebraic geometry point of view, the corresponding fan of labelled normals to the polytope does not lead to a toric variety, but to a generalization known as a "stack". The reason is that in ordinary toric varieties orbifold singularities along divisors cannot appear, as all the singularities are in co-dimension higher than one. In the broader context of stacks, the key combinatorial gadget is called a stacky fan, which roughly speaking is a simplicial fan with a distinguished lattice point on each ray of the fan (see e.g. [40][41][42][43]). For the more general Σ 1 ⋉ Σ 2 family we will not be able to find a symplectic structure and therefore we will have to use an indirect method to extract the toric data. These will again take the form of a stacky fan and will be taken as inputs for constructing an entropy function in terms of gravitational blocks in the next section. As further illustration of our approach, in appendix C we present the derivation of the Delzant polytope of the first Hirzebruch surface F 1 , starting from the explicit metric found in [44], as well as the toric data relevant for the seven-dimensional solutions of [28].

AdS 2 × S 2 ⋉ Σ 2 solutions
Before studying the Σ 1 ⋉ Σ 2 orbifold, in which case the analysis is more involved, we begin with a geometry characterized by a smooth two-dimensional base. Since a Riemann surface with g > 0 is not toric, we need to consider its analytic continuation to a twosphere. Although non-supersymmetric, the AdS 2 × S 2 ⋉ Σ 2 backgrounds form a family of solutions to the equations of motion very similar to (2.9), but with a different value of the radius of the Riemann surface, now a sphere. Specifically, the space we consider is the four-dimensional toric orbifold S 2 ⋉ Σ 2 with metric In this case, the condition on t (2.36) in order to have a well-defined fibration reads where t results to be negative. For fixed values of y, metric (4.1) describes an S 3 /Z −t , where S 3 is a squashed three-sphere written as a Hopf fibration, therefore a natural basis of an effective two-torus action is [44] 10 where we introduced the 2π-periodic coordinate ν 2 = 2π ∆z z. In total we have four fixed points v ℓ (ℓ = 1, . . . , 4) under the action of the two-torus, which are all the possible combinations obtained by pairing the poles of the sphere (θ = 0, π) and the poles of the spindle (y = y N,S ) We now consider the conformally rescaled metric ds 2 = Γ(y) ds 2 S 2 ⋉Σ 2 , with Γ(y) > 0, with compatible symplectic two-form 11 , the two-form (4.5) is closed and can be written as From this expression we can derive the moment maps with respect to the basis (4.3) and compute the image of the fixed points (4.4) which are vertices of the moment polytope. Vertices are connected by facets and we define the facet D ℓ to be the one that joins v ℓ−1 and v ℓ , where ℓ is modulo 4 12 . In order to correctly draw the polytope we notice that, since Γ ′ (y) < 0, we have Γ(y N ) > Γ(y S ). Due to the orbifold nature of the spindle, we can construct an associated labelled polytope following [39], which extends the classic Delzant construction. In particular, if the preimage under moment maps of a facet has an orbifold singularity locally modelled on C/Z m ℓ , we label the facet with the corresponding integer m ℓ , otherwise there is no label (formally, we have m ℓ = 1). Specifically, the S 2 ⋉ Σ 2 orbifold is characterized by The resulting labelled polytope is depicted in figure 1, and if the labels are stripped off this is exactly the polytope corresponding to the Hirzebruch surfaces F −t . This polytope can also be described uniquely by a set of inequalities µ, m ℓ n ℓ ≤ λ ℓ , (4.10) where ·, · is the standard inner product between vectors in R 2 , µ is a generic point in the (µ 1 , µ 2 )-plane, λ ℓ is a real number and each n ℓ is a primitive element of Z 2 , which represents the outward-pointing normal vector to the facet D ℓ . In particular, we have The method we employed to construct the labelled polytope above is rigorous and the result fits in the theory of symplectic toric orbifolds [39]. However, we will now study this toric orbifold from a different point of view, in order to derive a method that can be applied also to orbifolds which do not admit a conformally rescaled symplectic form. A direct computation using the metric (4.1) shows that there are four distinct loci D ℓ where a linear combination of the two Killing vectors ∂ ψ and ∂ z degenerates. These are precisely the toric divisors which in the symplectic context correspond to the preimages of the four facets of the polytope and intersect at the four fixed points. Explicitly, denoting with ξ (ℓ) the Killing vector that degenerates at D ℓ , they are normalized so to have unitary surface gravity (4.14) We notice that ξ (ℓ) = m ℓ n ℓ · (e 1 , e 2 ) , where we packed the basis vectors (4.3) into a two-dimensional vector (e 1 , e 2 ) and ( n ℓ , m ℓ ) are the toric data of the corresponding labelled polytope. We can use relation (4.15) to infer the toric data starting from the information about the degenerate Killing vectors. In particular, given a toric divisor D ℓ where the Killing vector ξ (ℓ) is degenerate, studying the (possible) orbifold singularity on D ℓ we can obtain the corresponding integer m ℓ . Expanding ξ (ℓ) on a suitable basis-in which the twotorus act effectively-according to (4.15) we can deduce the normal vectors n ℓ , up to a global sign depending on the conventions. This sign can be chosen arbitrarily for one of the vectors and for the others it follows requiring the vectors n ℓ to be ordered counterclockwise, in our conventions.
For an un-labelled, convex, rational polytope (not necessarily a Delzant polytope) one can pass to the algebraic point of view, where the primitive normals to the polytope represent the fan of a toric variety. In this context, the one-dimensional cones (rays) of the fan are in a one-to-one correspondence with the toric divisors and the order of the orbifold singularities at the intersection between two divisors can be read off from the determinant of two adjacent generators, namely det( n ℓ , n ℓ+1 ). Orbifold singularities along divisors cannot arise in toric varieties, so the information of a labelled polytope cannot be accommodated in an ordinary fan. An equivalent way of thinking about the labels m ℓ is by dropping the primitive requirement and considering n ℓ = m ℓ n ℓ instead of n ℓ [45]. We will refer to the n ℓ as "short" vectors and to the n ℓ as "long" vectors. In algebraic geometry the fan obtained from such "long" (non-primitive) normals can be understood in the broader context of stacks and it is referred to as a stacky fan. This turns out to be the correct combinatorial data needed in the construction of the entropy function in the next section. The "long" vectors characterizing the S 2 × Σ 2 orbifold then read Let us make some comments on the orbifold Σ g ⋉ Σ 2 , with g > 1. This cannot be toric in any sense, therefore none of the previous considerations can be applied. Nevertheless, for our purposes, we can think about extending the results obtained for the S 2 ⋉ Σ 2 toric orbifold to this case, and adopt as "toric data" the same vectors n ℓ (4.11) and the same integers m ℓ (4.9). The information on the genus of the Riemann surface does not enter in any way in the "toric data", but it will appear in the twisting procedure of section 5, as discussed in the examples. Clearly, for Σ g ⋉ Σ 2 the divisors D 2 and D 4 have no physical meaning, and likewise the four fixed points v ℓ , but they are all artefacts needed for the application of the conjectural entropy function we will present in section 5.

AdS
After the warm-up with the S 2 ⋉ Σ 2 case, we are now ready to study the Σ 1 ⋉ Σ 2 orbifold. In this case we can not construct a conformal symplectic structure, therefore we do not have a moment map to derive the labelled polytope. Nevertheless, we are still able to obtain the toric data that we need by analysing where the Killing vectors of Σ 1 ⋉ Σ 2 degenerate. Recall that the metric on Σ 1 ⋉ Σ 2 is (4.17) A simple direct calculation shows that it has four Killing vectors degenerating at the poles of the two spindles, defining the four loci intersecting at the the four fixed points normalizing the degenerating Killing vectors so that they have unitary surface gravity, these read where, in addition to ν 2 = 2π ∆z z, we introduced also the 2π-periodic coordinate ν 1 = 2π ∆ψ ψ. Reading off the labels associated with the D ℓ turns out to be a bit tricky. We will first present a naive argument leading to slightly wrong conclusions, and then we will explain how to amend this. Each of the Killing vectors in (4.20) is in a one-to-one correspondence with a specific pole of either Σ 1 or Σ 2 , whose rotation around the axis is generated by the 2π-periodic Killing vector ∂ ν 1 or ∂ ν 2 , respectively. This suggests that the integer coefficient in front of this rotational Killing vector corresponds to a C/Z m ℓ quotient singularity, with m ℓ = (n − , m − , n + , m + ). While for D 1 and D 3 this is correct, for D 2 and D 4 this conclusion is not precise, as we will show below.
In order to derive the normal vectors n ℓ we need to find an appropriate basis in which to expand the ξ (ℓ) . Inspired by (4.3), where ∂ ν 2 is the first basis element and ξ (4) plays the role of e 2 , we consider the following basis in which the degenerating Killing vectors are decomposed as where we made use of the quantization condition on t in (2.59) in the form Note that, even though x N,S and a are real numbers, all the coefficients in (4.22) are rational. Taking the labels as above, for the "short" vectors we obtain where we picked the sign of n 1 as in (4.11) and the others followed accordingly. However, since these vectors are not primitive elements of Z 2 , they cannot be regarded as the normal vectors to the facets of a polytope, in the sense of [39]. Even adopting the alternative description in which the labels are included in the normal vectors, the "long" vectors associated with (4.24) are not valued in Z 2 : Although in the mathematical literature there are examples where the normal vectors are allowed to take rational [45,46] as well as real values [47], the fact that the Killing vectors in (4.22) are not linear combinations of the basis {e 1 , e 2 } with integer coefficients, indicates that the two-torus is not acting effectively. A better basis can be obtained rotating the vectors n ℓ through the following SL(2, Q) matrix where r − is an integer prime to m + . Defining ŵ ℓ ≡ S n ℓ , we have where r ± ∈ Z are chosen such that t = r + m + + r − m − , and exist by Bézout's lemma 13 , for co-prime m ± . The vectors ŵ ℓ now take values in Z 2 , with primitive ŵ 2 , ŵ 4 , hence m 2 = m 4 = 1, and non-primitive ŵ 1 , ŵ 3 , corresponding to the labels m 1 = n − and m 3 = n + . We can now define a new basis {E 1 , E 2 }, obtained as E i = S −1 ji e j , such that and the degenerating Killing vectors take the form All the coefficients of this decomposition are now integers, indicating that in this basis the torus action is effective. For completeness we also write the basis vectors in terms of the Killing vectors ∂ ν 1 and ∂ ν 2 The combinatorial data can be presented in the form of a stacky fan, as depicted in figure 3a. Note that for n + = n − = 1 this reduces to the "Hirzebruch orbifold" discussed in [48] in the context of toric stacks. We can then read off the order of the quotient singularities at the intersections of the four divisors, corresponding to the four cones τ ℓ , which are given by C 2 /Z det( n ℓ , n ℓ+1 ) and thus read C 2 /Z n − m − , C 2 /Z n + m − , C 2 /Z n + m + , C 2 /Z n − m + , for ℓ = 1, . . . , 4 respectively. From the toric description we see that the total space Σ 1 ⋉Σ 2 is smooth not only away from the poles of Σ 2 , but also at the poles of Σ 1 , in agreement with [28]. Despite the fact that we do not have a moment map, it is of course possible to formally define a labelled polytope, encapsulating the same combinatorial data. This is obtained from the stacky fan in figure 3a, considering the vectors generating the rays of the fan and interpreting them as normal vectors to each facet of the polytope. The labels are read off from the "multiplicities" of the non-primitive vectors. The result is the rational convex polytope depicted in figure 3b, described by the labels m ℓ = (n − , 1, n + , 1) and the primitive normal vectors w ℓ = ŵ ℓ /m ℓ , which together are packed into what we shall refer to as toric data ( w ℓ , m ℓ ).

Toric divisors and their intersection
In this section we will discuss toric divisors with the aim of computing their intersection matrix, that will be a key ingredient in the construction of the entropy function. We will begin recalling some basic facts about divisors in the context of toric varieties (thus, ignoring the labels) and then we will indicate the modifications needed to pass to the broader setting in which our toric orbifolds fit. Recall that all the toric varieties defined by fans are normal and therefore all the (orbifold) singularities are in complex co-dimension higher than one. For toric surfaces M 4 (i.e. in real dimension 4), this implies that only points can be singular, while there cannot be singularities along divisors. Divisors that are torus-invariant are called toric divisors and are in a one-to-one correspondence with the primitive generators of the rays of the fan. These are in general Weil divisors, but for toric surfaces, since all the fans are necessarily simplicial, all Weil divisors D have an integer multiple lD that is a Cartier divisor-namely they are Q-Cartier divisors 15 . Below we will consider only such Q-Cartier toric divisors, associated with the primitive generators of a fan, v ℓ , that will be denoted as D ℓ . As in the previous sections, ℓ = 1, . . . , n runs over the generators of the fan in cyclic order. Furthermore, associated with any toric divisor D ℓ there is an equivariant line bundle L ℓ , whose first Chern class c 1 (L ℓ ) ∈ H 2 (M 4 , Q) is Poincaré dual to D ℓ . While the second co-homology group of M 4 , H 2 (M 4 , Q), is (n − 2)dimensional, its equivariant extension, generated by the c 1 (L ℓ ), has clearly dimension n.
The intersection matrix of the toric divisors is defined as and by a standard calculation in toric geometry, it reads where we used the notation D ℓ,ℓ ′ [ v ℓ ] to emphasize that it is the intersection matrix of the toric divisors associated with a fan with generators v ℓ . It can be easily proven that this satisfies the relation implying that D ℓ,ℓ ′ , viewed as an n×n matrix, has two null eigenvalues, and it is therefore not invertible. This of course is consistent with the fact that dim(H 2 (M 4 , Q)) = n − 2 implies that only n − 2 toric divisors D ℓ are independent in homology. For the special case of toric manifolds, all det( v ℓ , v ℓ+1 ) = 1, and the above discussion reduces to that in [31]. On the other hand, for our purposes we need to incorporate the effect of the labels. We will refrain from attempting a rigorous treatment, indicating a natural way to take into account this additional information. If we have a labelled polytope/stacky fan described by "long" vectors v ℓ = m ℓ v ℓ , it is natural to associate a divisorD ℓ with v ℓ , which is related to D ℓ simply bŷ 15 Varieties with such a property are called Q-factorial.
In particular, this implies that for any two-form Λ the integrals over the two sets of divisors are related as m ℓ (4.35) and the intersection matrix of the divisorsD ℓ is given bŷ Correspondingly, the first Chern classes of the dual line bundles are related as c 1 (L ℓ ) = m ℓ c 1 (L ℓ ). Since we can use either c 1 (L ℓ ) or c 1 (L ℓ ) as basis for H 2 (M 4 , Q), these two descriptions are essentially interchangeable. We note thatD ℓ,ℓ ′ is invariant under SL(2, Z) transformations of the basis and more generally for any transformation in SL(2, R). This follows immediately from the fact that given any two vectors v 1 , v 2 ∈ R 2 , det( v 1 , v 2 ) is invariant under SL(2, R) transformations of the two vectors. In particular, taking a generic matrix S ∈ SL(2, R) and transforming the vectors as v ′ a = S v a , we have that S acts on the matrix ( v 1 , v 2 ) simply by matrix This property implies that formally the intersection matrix of the divisors computed in the two bases ŵ ℓ and n ℓ are identical, being related by the SL(2, Q) transformation (4.26).
In order to make contact with the computations in supergravity, we need to describe the relation between the toric divisors defined above and the representative of the cycles, that we used to integrate fluxes in our main example of Σ 1 ⋉ Σ 2 . The S ± introduced in section 2.2 were precisely two copies of the base Σ 1 at the south pole (section at infinity) and north pole (zero section) of the fibre Σ 2 and we can identify these with the toric divisors D 1 and D 3 as On the other hand, D 2 and D 4 correspond to copies of the fibre Σ 2 at the poles of the base Σ 1 . In the smooth case (m + = m − = 1), these would be simply homologous to Σ 2 . However, the orbifold singularities of the base imply that instead where X ± ≡ Σ 2 /Z m ± , respectively. Therefore, in order to compare the integral of a twoform Λ over X ± with the integral over Σ 2 at a generic point of Σ 1 we should use the relation m ± Since in section 2.2 the fluxes were computed integrating on Σ 2 (see e.g. (2.51)), the integrals computed using X ± will have to be multiplied by a factor of m ± to be compared with the fluxes computed in section 2.2.

The entropy function
We conjecture that the solutions presented in section 2.2 are holographically dual to onedimensional SCQMs obtained compactifying on a Σ 1 ⋉Σ 2 the five-dimensional SCFTs dual to the solution of [49]. In this section we provide supporting evidence for this conjecture, proposing an entropy function, whose extremization reproduces the entropy (3.30) of the (putative) black hole near-horizon AdS 2 × Σ 1 ⋉ Σ 2 . Starting from first principles, this function should be derived from the localized partition function of the d = 5 SCFT, placed on the background of S 1 × Σ 1 ⋉ Σ 2 , and then taking the large N limit. Following the idea of "gravitational blocks" [30], we will propose a large N entropy function on S 1 × Σ 1 ⋉ Σ 2 obtained by suitably gluing the S 5 free energy of the d = 5 theories. With some further assumptions, this prescription can be applied also to the AdS 2 × Σ g ⋉ Σ 2 solutions of section 2.1. The prescription that we propose extends the results of [11] to a broader class of configurations, including solutions arising in D = 7 supergravity.

Twisting data
In this section we consider the compactification of five-and six-dimensional SCFTs on a four-dimensional toric orbifold M 4 . Specifically, we place the theory on M 4 and couple it to two background gauge fields A i , with appropriately quantized magnetic fluxes. All the information about M 4 is encoded in the toric data ( v ℓ , m ℓ ), which can be thought of as the data defining a labelled polytope. The vectors v ℓ must be primitive, Z 2 -valued and, in our conventions, ordered counter-clockwise. The recipe we will formulate is in principle applicable to toric orbifolds (including smooth toric manifolds) with an arbitrary number n ≥ 3 of fixed points. Consider the two line bundles E i on which the one-forms −A i are the connections. Since the set of c 1 (L ℓ ) form a basis for the equivariant extension of H 2 (M 4 , Q), we can decompose the first Chern class of E i as The "physical fluxes" are defined as 16 and using (5.1) we obtain the relation where D ℓ,ℓ ′ is the intersection matrix defined in (4.32). We recall that H 2 (M 4 , Q) has dimension n − 2, thus the whole set of p (ℓ) i is a redundant parameterization. In particular, for fixed i, only n − 2 of the p (ℓ) i are linearly independent. The two additional degrees of freedom are the expression of a gauge symmetry that we will discuss in the next subsection.
A further consequence of relation (4.33) is that, combined with (5.3), it gives This is a vector equation and, for fixed i, we have n fluxes q (ℓ) i constrained by two equations, one for each component of the two-dimensional vectors v ℓ . Therefore, only n − 2 physical fluxes are linearly independent. In what follows, given a generic toric orbifold with arbitrary fluxes, we shall require them to satisfy the constraint (5.4). We can now go back to relation (5.3). As we explained, it represents a linear system of rank(D ℓ,ℓ ′ ) = n−2 independent equations for n − 2 independent unknowns p where m ℓ are the labels of the labelled polytope associated with M 4 and, a priori, σ (ℓ) are n arbitrary coefficients. Denoting as E R the R-symmetry (orbifold) line bundle, equation (5.5) can be rewritten as 6) or, following standard practice, as c 1 (E R ) = − ℓ σ (ℓ)D ℓ . Note that the standard topological twist corresponds to identifying E R with the orbifold canonical line bundle K orb M 4 = − ℓD ℓ and hence σ (ℓ) top-twist = (+, . . . , +). See e.g. [50] for a related discussion of orbifold line bundles.
In analogy with the case of the spindle [12], we conjecture that the only possible values of σ (ℓ) are ±1. This hypothesis is supported by the different examples analysed, as we shall see. In other words, we expect that it should be possible to compactify a SCFT on a toric orbifold M 4 , turning on a background R-symmetry gauge field with magnetic fluxes given in (5.5), with σ (ℓ) parameterizing the different supersymmetry-preserving twists. Of course, we have not proven that all these different twists preserve supersymmetry, nor that there cannot exist more general twists. It would be interesting to carry out such an analysis extending the results of [12], where it was demonstrated that, in whole generality, on a spindle the only two supersymmetry-preserving twists are the twist and the anti-twist.

The recipe
We conjecture that given a general class of SCFTs in d = 5, 6 compactified on a fourdimensional toric orbifold M 4 , with an arbitrary twist parameterized by a set of signs σ (ℓ) as in (5.5), the corresponding entropy/central charge, respectively, is determined extremizing the following off-shell free energy (5.7) The sum runs over the n fixed points v ℓ of the toric orbifold M 4 , k d is a numerical constant, which, a posteriori, takes the values The variables Φ with ǫ = (ǫ 1 , ǫ 2 ) 17 . Here, ǫ 1 and ǫ 2 may be interpreted as the fugacities associated with the two U(1) rotational symmetries and parameterize their mixing with the R-symmetry. η (ℓ) d = ± are signs that at this stage must be tuned by hand: in particular we believe that in general 18 in d = 6 they are all +, while for d = 5 they are related to the type of twist and we speculate that η (ℓ) 5 = σ (ℓ) σ (ℓ+1) . In the next section we will see some explicit examples applying this prescription to the AdS 2 × Σ 1 ⋉ Σ 2 family. In this construction ϕ i and ǫ i are variables with respect to which one has to extremize F (ϕ i , ǫ i ), v ℓ and m ℓ are data describing the toric orbifold M 4 and p The variables ϕ i and ǫ i are subject to the constraint where W is a two-dimensional constant vector parameterizing a "gauge invariance" of the problem, that we discuss below. This condition is inherited from the R-symmetry constraint ∆ 1 + ∆ 2 = 2, where ∆ i are the fugacities parameterizing the R-symmetry within the Cartan subgroup of the global symmetries of the d-dimensional parent theory. The vector W can be determined imposing 13) 17 The two variables ǫ 1,2 are different from the ones adopted in [11]. Specifically, ǫ here i = 2ǫ there i . 18 Due to different conventions they could also be all −, as explained at the end of subsection 5.4. which is a system of n equations -one for each value of ℓ -for the two components of W , and it is thus overdetermined. Yet, a solution always exists. By means of equations (5.3) and (5.13) we can write the following chain of equalities 14) where, in the last step, we made use of (4.33). First of all, this relation proves the consistency of (5.13), since, as it should be, the physical quantities q (ℓ) R do not depend on the auxiliary vector W . Moreover, it shows that n − 2 independent linear combinations of equations (5.13), obtained contracting them with D ℓ,ℓ ′ , are already satisfied once (5.5) is imposed. The result is a linear system of two independent equations with two unknowns, that can always be solved.
As we already mentioned, the quantity det( v 1 , v 2 ) is invariant under SL(2, R) transformations of the two vectors v 1 , v 2 ∈ R 2 . This property is fundamental for the consistency of our construction. A given toric orbifold can be described by an infinite number of polytopes equivalent under SL(2, Z) transformations of the vectors v ℓ . These transformations are generated by rotations of the basis vectors (e 1 , e 2 ) that give rise to effective torus actions. On the other hand, we expect the free energy not to be affected by these transformations and, indeed, this is the case. If we perform an SL(2, Z) rotation of the vectors v ℓ , the expression of the intersection matrix D ℓ,ℓ ′ is retained, and the explicit form of the off-shell free energy is not modified provided we apply exactly the same transformation to ǫ and W .
The off-shell entropy function we constructed enjoys also another symmetry, the gauge symmetry we mentioned before. Thanks to (4.33), relation (5.3) is left invariant under the "gauge transformation"p for any two-dimensional constant vector λ i . This transformation affects the value of the vector W appearing in the constraint. Specifically, since σ (ℓ) and m ℓ are fixed, we must imposep which returns W = W + λ 1 + λ 2 . Introducing the new variablesφ i = ϕ i + det( λ i , ǫ) it is possible to keep the functional expression of Φ p). Remarkably, the constraint changes accordinglỹ As a result, both the off-shell free energy and the constraint retain the same functional expression under the gauge transformation (5.15), thus the value of the former at its critical points will be the same. The proposed prescription takes inspiration from [31] and extends their construction to compactifications on generic toric orbifolds. Moreover, it also applies to d = 5, generalizing the formulas presented in [11] for compactifications on Σ g 1 × Σ g 2 and Σ × Σ g , to the case of fibred spaces. The constrained extremization of F (ϕ i , ǫ i ) can be performed defining the function and extremizing it with respect to ϕ i , ǫ i and the Lagrangian multiplier Λ. In order to solve this problem it might be helpful to notice that (5.7) is homogeneous of degree h in ϕ i and ǫ i (h = 1 for d = 5, h = 2 for d = 6). As a consequence, by Euler's theorem, we

Application to the AdS
The first example to which we can apply the prescription presented in the previous section is the AdS 2 × Σ 1 ⋉ Σ 2 solution (2.40). We now briefly collect some of the results needed in the construction. The toric orbifold M 4 = Σ 1 ⋉ Σ 2 can be described by the toric data (see section 4.2) with t = r + m + + r − m − . The vectors w ℓ are integers, primitive and ordered counterclockwise, as required by the recipe. Keeping in mind all the observations about divisors made in section 4.3, the physical fluxes defined in (5.2) read where s + i , s − i and t i are given, respectively, in (2.61), (2.60) and (2.28). The intersection matrix D ℓ,ℓ ′ , computed straightforwardly using (4.32), reads i /m − = p (4) i /m + . Therefore, the relation (5.3) yields the identifications In particular s * 1 explicitly reads and s * 2 can be obtained flipping the sign of z. With this definition, s ± i can be written as System (5.13) can be solved to get W , which gives the constraint (5.12) The ingredients ǫ (ℓ) 1,2 and Φ (ℓ) i can be constructed by means of (5.11) and (5.10), whereas the off-shell free energy (5.7) reads where the relative signs have been fixed using the prescription η (ℓ) . We can now proceed to the extremization of the off-shell free energy and, after some work, we get the critical values where the sign ambiguity η = ±1 arises by solving the equations over the complex numbers 19 . Inserting these values back into (5.27) we obtain where F S 3 ×Σ 2 is given in (3.24). In order to get a positive entropy we need to pick η = −1 and in this case the result agrees with the gravitational entropy (3.30). We propose that the procedure we formulated can be extended also to non-toric fourdimensional orbifolds, e.g. in presence of Riemann surfaces as submanifold. In this case we can adopt the same toric data obtained replacing Σ g (g > 1) with S 2 , setting g = 0 whenever it occurs in (5.13) and (5.5). The genus is preserved in the explicit expression of the fluxes. Specifically, we will test our conjecture against the AdS 2 ×Σ g ⋉Σ 2 system (2.9). The "toric data" are given in (4.11) and (4.9), the physical fluxes are the same as in the previous example, with m ± = 1, and now s ± i can be read in (2.38) and (2.37). The fluxes p (ℓ) i are identified as before and s * i retains the same form, but, again, with m ± = 1 and t as in (2.36). The constraint on the variables ϕ i and ǫ i is whereas the off-shell free energy reads where again the relative signs were chosen to be η (ℓ) 5 = σ (ℓ) σ (ℓ+1) . These different gluing signs are due to the different types of twist that we have. In particular, on both Σ g and Σ 2 the twist is realized, in contrast to the previous case, in which we had anti-twist on Σ 1 . The off-shell free energy is extremized by 32) to which corresponds the critical value Also in this case, the result returned by the extremization agrees with the entropy computed from the ten-dimensional supergravity solution (3.23).

Application to the AdS
The conjectured prescription presented at the beginning of this section can also be applied to another interesting family, the AdS 3 × Σ 1 ⋉ Σ 2 solutions constructed in [28] 20 . The toric data of the Σ 1 ⋉ Σ 2 toric orbifold are derived in appendix C.3 and we recall them here for the reader's convenience where p 1 and p 2 are two integers related to the magnetic fluxes as in (2.11) of [28]. The physical fluxes can be computed as in the previous example and, compactly, they read where s − 1 corresponds to the flux computed in (3.15) of [28], s − 2 and s + i can be derived from the former and t i = p i n + n − . The gauge symmetry allows us to impose p i /m + , which yields, through (5.3), the identifications For later convenience we defined s * 1 ≡ , whose expression is The fluxes s ± i are related to s * i through The constraint can be derived following the path traced in the previous case and reads Notice that this expression is identical to (5.26). As we discussed in section 5.2, in d = 6 we set η (ℓ) 6 = + for all ℓ. Remarkably, this implies that F is quadratic in the ϕ i , as expected for the off-shell central charge of a two-dimensional SCFT [51]. Explicitly, this reads The extremization procedure, realized by means of Lagrangian multipliers, gives the critical values 21 ϕ which, plugged into the off-shell central charge, gives where a 4d is the central charge of d = 4, N = 1 SCFTs that arise from N M5-branes wrapped on a spindle [8] a 4d = 3p 2 1 p 2 2 (s + p 1 + p 2 ) 8n + n − (n + − p 1 )(p 2 − n + )[s + 2(p 1 + p 2 )] 2 N 3 . (5.44) This matches exactly the central charge of the AdS 3 ×Σ 1 ⋉Σ 2 system computed in [28]. In general, we expect that the anomaly polynomial computation of [28] and the prescription presented in this paper should be equivalent and connected by a suitable gauge choice and a possible redefinition of ǫ 1,2 , mixing the two related U(1) isometries. As explained in the previous section, our recipe can be applied also to non-toric fourdimensional orbifolds. In this class we consider AdS 3 × Σ g ⋉ Σ 2 , whose "toric data" are deduced in appendix C.2 for AdS 3 × S 2 ⋉ Σ 2 and read Global regularity leads to , (5.46) with s given in equation (5.35). We take the physical fluxes to be as in (5.36), where t i = p i n + n − is unchanged and s ± i can be read from (4.7) of [28]. The gauge symmetry allows us to identify the p (ℓ) i as in equation (5.37), with s ± i = s * i ± t 2 t i and We can then solve equation (5.5) to obtain σ (ℓ) = (+, −, +, −), use it to identify the auxiliary vector W through (5.13) and write down the constraint (5.12) This is the same as in (5.30), with t → −t and ǫ 1,2 → −ǫ 1,2 , due to the different conventions used in [28]. Taking all the η as well as the off-shell central charge at the extremum with a 4d given in (5.44). This result agrees with the central charge of the AdS 3 × Σ g ⋉ Σ 2 system studied in [28]. The origin of the choice of the η (ℓ) 6 can be traced to the sign of the gauge potential of the five-dimensional solution that the authors of [28] uplifted to obtain the AdS 3 × Σ g ⋉ Σ 2 backgrounds. The component along AdS 3 of their Killing spinors satisfy the corresponding Killing spinor equations, but with a minus, namely∇μϑ = − 1 2 αμϑ, which implies that in the two-dimensional dual SCFTs (2, 0) supersymmetries are preserved, instead of (0, 2) as in our case. This fact leads to a different sign in the extraction of the central charge from the anomaly polynomial.

Discussion
In this paper we have presented two new families of supersymmetric AdS 2 solutions of massive type IIA supergravity, associated with D4-branes wrapped on four-dimensional orbifolds M 4 , which may be viewed as two different generalizations of the Hirzebruch surfaces. In one case M 4 = Σ g ⋉Σ is a spindle Σ fibred over a smooth Riemann surface Σ g of genus g > 1, while in the other case M 4 = Σ⋉Σ is a spindle Σ fibred over another spindle Σ. We have argued that these are dual to N = 2 SCQMs arising from different twisted compactifications of the d = 5, N = 1 USp(2N) supersymmetric gauge theories on M 4 . The structure of these solutions is analogous to that of the AdS 3 solutions constructed in [28] and indeed for the M 4 = Σ⋉Σ family we have provided a toric geometry description and formulated a conjectural extremal problem that applies to both classes.
As in [28], the solutions that we found have less parameters than the data specifying the orbifolds M 4 . In particular, the twisting parameter t is not an arbitrary integer, but it is unnaturally related to the spindle data n ± , m ± . We therefore expect more general solutions to exist, although the ansatz describing them is likely to be quite different from those employed here and in [28]. Furthermore, the toric description of the M 4 = Σ ⋉ Σ solutions strongly suggests that there exist solutions corresponding to toric orbifolds with an arbitrary number of fixed points, and could include the case of toric manifolds.
While in this paper we provided a basic description of the orbifolds, which was adequate for characterizing the present solutions, it would be desirable to study more systematically their underlying structure. For example, it would be useful to perform a global analysis of the gauge fields and Killing spinors on M 4 , determining the bundles of which these are sections. Relatedly, given a generic compact four-dimensional toric orbifold M 4 , it would be nice to determine what are the allowed supersymmetry-preserving twists. It could also be useful to cast our solutions in terms of classifications of supersymmetric solutions of massive type IIA supergravity, or to characterize them directly in the context of D = 6 and D = 7 gauged supergravities. For example, we expect that generic M 4 's will be endowed with a canonical complex structure 22 . Another natural extension of our results is to investigate whether supersymmetric AdS 6 black holes with M 4 horizons can be constructed.
In this paper we have proposed an entropy function constructed assembling gravitational blocks, extending the proposal of [11] from the spindle to a very general class of four-dimensional orbifolds. We believe that our approach will be applicable to orbifolds more general than those for which we constructed explicit solutions here. However, finding additional explicit examples may lead to a more refined formulation of our recipe. Similarly to the spindle off-shell free energies [11], it is striking that the structure we introduced applies to both D = 6 and D = 7 supergravities, with minor tweaks. In the context of D = 7 supergravity, our proposal reproduces the results obtained by integrating the M5-branes anomaly polynomial on compact toric four-manifolds [31] and on the M 4 = Σ ⋉ Σ orbifolds [28]. On the other hand, in the context of D = 6 supergravity, at present we do not have access to any field theory computation and therefore our proposal gives a prediction for the large N limit of the localized partition function of the d = 5,

A Some solutions to the Diophantine equation
In this appendix we address the quantization conditions we stated through the paper. Although we were not able to prove that these constraints can be satisfied appropriately tuning the parameters of the solutions, nevertheless we could find some examples in which the quantization can be realized.
A priori, there are no constraints on x apart from being solution of the quartic equation (2.33) and lying in the range (0, 1). Nevertheless, in order to construct the integer quantities we are searching for, it is easier to work with rational values of x, if they exist. Indeed they exist, and in table 1 we present some of them, along with the corresponding values of the integers n + , n − and p 1 (recall that p 2 = n + + n − − p 1 ). Inspired by [28] we propose an algorithm to construct a family of solution to the quantization conditions. For the combinations of parameters presented in table 1, x is rational and positive, thus we can write it as x = p/q, with p, q ∈ N. Then, we define t as with u ∈ N, which makes t integer by construction and, moreover, multiple of both n + and n − . We now focus on the AdS 2 × Σ 1 ⋉ Σ 2 system, the case with Riemann surface proceeds similarly. Substituting the expression of t in the constraint (2.59) we obtain which can be solved to find m ± in terms of n ± , z and u. The quantization condition that follows from (2.60) is automatically satisfied given t as in (A.1), indeed it becomes which is evidently an integer. The last requirement (3.29) can be met considering that 3µ(x 2 + 1) − x(x 2 + 5) > 0 (cf. [11]) and that all the quantities appearing are rational, hence N can be tuned appropriately so that K ∈ N.

B Riemann surfaces from spindles
In presence of twist, it was explicitly shown in [12] that the spindle factor of the known AdS D−2 × Σ solutions, with D = 4, 5, 7, can be turned into a round sphere performing a particular limit of the coordinates and parameters of the system. In the same spirit, as we will see in this appendix, a similar mechanism applies to the AdS D−2 × Σ solutions with anti-twist, so far constructed only in D = 4 and D = 5. In this case the spindle becomes the two-dimensional hyperbolic space, which can later be quotiented to give a constant curvature Riemann surface with genus g > 1.

B.1 AdS 2 × M 4 solutions
Let us begin considering the local spindle metric (2.52) where, we recall, q(x) = x 4 − 4x 2 + 4ax − a 2 . Performing the change of coordinates and computing the limit ǫ → 0 + , the spindle metric becomes which is the local metric of the hyperbolic space H 2 , with radius squared equal to 1/2. This space can then be quotiented to obtain a Riemann surface with g > 1. In order for metric (B.3) to genuinely describe H 2 , the coordinates ξ and φ must have the correct range of definition, i.e. |ξ| ∈ (1, +∞) and φ ∈ [0, 2π]. To address this issue, we begin considering that the coordinate x is restricted between x N and x S , where [5] x By means of (B.2), these bounds are mapped to in the limit ǫ → 0 + , thus correctly giving ξ ∈ (−∞, −1). The periodicity of φ can be studied analysing the conditions under which Σ 1 is proper spindle with m ± two co-prime integers such that m − < m + . By means of the change of coordinates (B.2), these conditions yield hence we have ∆φ = 2π if and only if, at first order, we see that it diverges when m − → 0 and m + = 1. The reason for this behaviour can be found in the definition of Σ g as a quotient of the hyperbolic space, whose Ricci curvature, integrated over the whole space, diverges, since it is not compact. Therefore, we need a way to "regularize" this limit. To this end, we notice that, adopting the substitution (B.8), that is, we obtain a simple pole at ǫ = 0 plus the Euler characteristic of a Riemann surface. We propose the following prescription, that we will test against different examples: given a system comprising a Σ g factor that can be obtained from Σ 1 through limit (B.8), the global quantities of the former can be computed multiplying by 2(1−g) the corresponding quantities of the latter and taking the finite term in the Laurent series.
We can now transfer this technology to higher-dimensional spacetimes. Applying the change of coordinates (B.2)-driven by the underlying limit (B.8)-to the AdS 2 × Σ 1 system described by (2.68) we obtain the AdS 2 × Σ g background in (2.67). Likewise are related the respective Killing spinors. Indeed, in this limit Q 1 vanishes and Q 2 becomes constant, thus (2.76) turns into (2.75), apart from an irrelevant overall factor. For what concerns the total charge of the two systems, applying the recipe illustrated for the Euler characteristic χ 1 we obtain (cf. (2.56)) where a −1 is the residue. The change of coordinates (B.2) also maps the AdS 2 × Σ 1 ⋉ Σ 2 solution and Killing spinors to the AdS 2 × Σ g ⋉ Σ 2 corresponding quantities. This relation reflects in a formal connection between the respective entropies by means of (B.8): where S Σ 1 ⋉Σ 2 is given in (3.30) and S Σg⋉Σ 2 in (3.23).

B.2 AdS 3 × M 4 solutions
Moving to the local spindle metric presented and studied in [4] 23 where q(x) = 4x 3 − 9x 2 + 6ax − a 2 , we perform the change of coordinates (B.14) In the ǫ → 0 + limit, the spindle metric becomes Once again, we can quotient the hyperbolic space to give a constant curvature Riemann surface with genus g > 1. It can be shown that the metric obtained in this limit describes a genuine hyperbolic space if Applying this recipe to the AdS 3 × Σ 1 solution of [4] we have where ω H 2 is such that dω H 2 = vol(H 2 ). This is the AdS 3 × H 2 near-horizon of the black string solution of D = 5, N = 2 minimal gauged supergravity constructed in [53] 24 . Like in the six-dimensional case, applying the change of coordinates (B.14), followed by the ǫ → 0 limit, to the AdS 3 × Σ 1 ⋉ Σ 2 solution of [28] gives the AdS 3 × Σ g ⋉ Σ 2 background therein. Moreover, it is possible to retrieve the central charge of the latter, presented in (5.50), multiplying the central charge of the former, given in (5.43), by 2(1 − g), substituting m ± as in (B.16) and taking the finite part of the Laurent series. Explicitly,

C More toric data from explicit metrics
In this appendix we will provide further examples of how to extract toric data from explicit metrics, as discussed in the main body of the paper. The first example we consider is the Hirzebruch surface F 1 , that is a well known symplectic toric manifold. Its associated polytope is a Delzant polytope (see e.g. [54]) and below we show how to obtain it starting from an explicit metric on F 1 [44], using two methods. Namely, computing the image of the moment map associated with a symplectic structure compatible with a conformally rescaled metric and studying the Killing vectors generating the U(1) 2 isometry of the metric. We will then employ these methods to extract the toric data associated with the metrics presented in [28]. These are used in section 5 to show that our recipe reproduces the central charge computed in [28].

C.1 Hirzebruch F 1
We now illustrate an example of Delzant polytope and its construction studying the first Hirzebruch surface F 1 , whose metric is given by (cf. equation (7.11) of [44]) 25 ds 2 F 1 = 1 − y 6 dθ 2 + sin 2 θ dφ 2 + dy 2 w q + l 2 w q 36F (dΩ + cos θ dφ) 2 , (C.1) where the expression of the functions w(y), q(y) and F (y) can be found in [44]. The angular coordinates lie in the ranges θ ∈ [0, π], φ ∈ [0, 2π], Ω ∈ [0, 4π], while y spans the range [y 1 , y 2 ], where y 1 and y 2 are, respectively, the negative and the smallest positive roots of q(y). To ensure the correct signature of the metric we require 1 − y > 0. After having rescaled the metric as ds 2 = Γ(y) ds 2 F 1 , we consider the symplectic two-form ω = Γ(y) 1 − y 6 sin θ dθ ∧ dφ + l 6F 1/2 dy ∧ (dΩ + cos θ dφ) . (C.2) When Γ(y) is such that Γ ′ (y) = F 1/2 −l (1−y)F 1/2 Γ(y), the two-form (C.2) is closed and can be written as For fixed values of y, metric (C.1) describes a squashed, smooth three-sphere, corresponding to the system studied in section 4 of [44] for m = 1. Following their arguments, we take as a basis of an effective torus action 26 where V 1 generates the rotations of the two-sphere about the equator with weight one, hence V 1 = ∂ φ , and V 2 the rotations of the S 1 fibre, also with weight one, hence V 2 = ∂ ν , where we introduced the new radial coordinate ν = Ω/2, with period 2π. Once we have chosen a basis for the torus action, we can derive the related moment maps µ = 1 − y 3 Γ(y) 1, cos θ + 1 2 . (C.5) In total we have four fixed points, corresponding to all the possible combinations pairing the poles of the base two-sphere (θ = 0, π) and the poles of the fibre two-sphere (y = y 1,2 ) where we defined, for convenience, the positive constants Ξ i = 1−y i 3 Γ(y i ). A numerical analysis shows that Ξ 1 > Ξ 2 , therefore the moment polytope of F 1 can be drawn as in figure 4a. A quick investigation proves that it is a Delzant polytope in R 2 .
The polytope is uniquely defined by the relations µ, n 1 ≤ Ξ 1 , µ, n 2 ≤ 0 , µ, n 3 ≤ −Ξ 2 , µ, n 4 ≤ 0 , (C.8) where the Z 2 primitive normal vectors n ℓ are n 1 = (1, 0) , n 2 = (−1, 1) , n 3 = (−1, 0) , n 4 = (0, −1) . (C.9) The fan generated by these vectors is depicted in figure 4b and is dual to the Delzant polytope of figure 4a.   A different approach to the construction of the toric data of a given toric manifold goes through the analysis of the degenerating Killing vectors of the corresponding metric. In total we have four Killing vectors ξ (ℓ) degenerating at the loci D ℓ , namely We notice that the four loci D ℓ are in a one-to-one correspondence with the facets of the polytope in figure 4a and that the Killing vectors ξ (ℓ) can be written in terms of the related normal vectors n ℓ as ξ (ℓ) = n ℓ · (e 1 , e 2 ).
We can now construct the polytope, considering the four-dimensional conformallyrescaled metric ds 2 = Γ(y) ds 2 S 2 ⋉Σ 2 and the related symplectic two-form ω = Γ(y) 1 3 sin θ dθ ∧ dψ + y 1/2 2P 1/2 dy ∧ dz − (C.25) Since we have chosen the same basis as in (4.3), the resultant polytope is similar to the one constructed in section 4.1. Indeed, it can be read off from figure 1 exchanging v 1 ↔ v 3 , v 2 ↔ v 4 and n + ↔ n − . The normal vectors n ℓ change accordingly with respect to equation (4.11), i.e. n 1 ↔ n 3 and n 2 ↔ n 4 , together with t → −t. We can quickly verify that the quantities we derived satisfy relation (4.15).

C.3 AdS 3 × Σ 1 ⋉ Σ 2 solutions
We now move to the AdS 3 × Σ 1 ⋉ Σ 2 solutions presented in section 3 of [28] 29 . The metric on the toric orbifold is given by where functions h i , P and Q are the same as in (C.13) and The coordinate x ranges between x 1 and x 2 , the two smallest roots of f (x) with x 1 < x 2 , and y lies between y 2 and y 3 as before. The fibration is globally well-defined if which gives the relation t m + m − = − a(x 2 − x 1 ) 3x 1 x 2 ∆ψ ∆z .