Microstate counting of AdS4 hyperbolic black hole entropy via the topologically twisted index

We compute the topologically twisted index for general N=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}=2 $$\end{document} supersymmetric field theories on ℍ2×S1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{\mathbb{H}}}_2\times {S}^1 $$\end{document}. We also discuss asymptotically AdS4 magnetically charged black holes with hyperbolic horizon, in four-dimensional N=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}=2 $$\end{document} gauged supergravity. With certain assumptions, put forward by Benini, Hristov and Zaffaroni, we find precise agreement between the black hole entropy and the topologically twisted index, for ABJ M theories.


Introduction
Black holes have an entropy that fits neatly in a thermodynamics framework as originally established in the works of Bekenstein and Hawking in the early 1970's. The microscopic origin, that is, the nature of the degrees of freedom that this entropy counts, has been an outstanding challenge for many decades. Any candidate to a theory of quantum gravity must provide an answer to this fundamental question. String theory, in the works of Strominger and Vafa, has successfully passed this test for a particular type of black holes [1]. In the context of the AdS/CF T correspondence, the original work of Strominger and Vafa can be interpreted as an instance of AdS 3 /CF T 2 . A natural question pertains higher dimensional versions of the AdS/CF T correspondence. Recent work by Benini, Hristov and Zaffaroni addresses the microscopic counting of the entropy of certain black holes from the point of view of AdS 4 /CF T 3 [2].
In this manuscript we explore the topologically twisted index, originally introduced by Benini and Zaffaroni in the framework of N = 2 supersymmetric three-dimensional field theories in S 2 × S 1 [3] (see also [4][5][6][7]), for the case of supersymmetric theories in H 2 × S 1 , where H 2 is the hyperbolic plane. Although we provide the ingredients for arbitrary N = 2 supersymmetric theories, we will particularize our results for a specific deformation of ABJM theory. The holographic dual of such deformation is thought to be a hyperbolic black hole. In this work, our main motivation comes from the prospect of understanding the D = 3 SCFT representation of the appropriate AdS 4 black hole microstates. With this aim we are driven to explore four dimensional N = 2 gauged supergravity and find black hole solutions with H 2 horizon. Hyperbolic black holes have been discussed in the context of AdS/CF T in, for example, [8].
Asymptotically AdS 4 black holes in N = 2 gauged supergravity, which are sourced by magnetic fluxes, have been widely studied [9][10][11][12][13]. Roughly speaking, from the bulk perspective, the presence of fluxes allows to define the black hole as interpolating from the UV AdS 4 to the near horizon AdS 2 × S 2 . As a result of our study we are able to identify the role of such fluxes from the dual SCFT perspective. These flavor fluxes, together with a continuous of color fluxes, generate a one-parameter hierarchy of Landau levels on H 2 , that determines the value of the ABJM index. What we are set to explore in this paper, is whether the leading behavior in the large N limit of the topologically twisted index of a specific deformation of ABJM , evaluated on the Hilbert space composed by the aforementioned Landau levels, coincides with the Bekenstein-Hawking expression for the semiclassical entropy of the black holes in question. We will find that indeed both results coincide.

JHEP08(2017)023
Another important motivation for our work, is the intrinsically interesting field theory problem of localization of supersymmetric field theories in non-compact spaces. This problem naturally appears in the context of localization of supergravity theories, for an understanding of exact black hole entropy counting [14,15]. The same problem appears in holographic approaches to Wilson loops where the world volume of the classical configuration contains an AdS 2 factor. For example, the excitations on a D3 brane which is dual to a Wilson loop in the totally symmetric rank k representation [16] were identified to correspond, to an N = 4 vector multiplet in H 2 × S 2 [17]. Localization in non-compact spaces has recently been addressed in [18] and [19], our work constitutes an extension to the topologically twisted case.
The manuscript is organized as follows. In section 2 we discuss the preliminary ingredients we need, for example, our guidance principle on the field theory side: supersymmetric localization [20], the background metric, spin connection, and supersymmetric structure of the actions needed to compute BP S observables in a generic three-dimensional N = 2 Chern-Simons-Matter theory on H 2 × S 1 . To complete section 2, we discuss the boundary conditions to be used in the manuscript. In section 3 we present the space of square and delta-normalizable functions that will be used to integrate upon, and their respective discrete and continuous spectrum. In section 4 we compute the one loop super-determinants. In section 5 we assemble our results to write down the ABJM index on H 2 × S 1 , and then move on to compute its leading contribution in the large N expansion, by following the procedure pioneered in [2]. In section 6 we find what we believe to be the dual AdS 4 black holes and compare its Bekenstein-Hawking entropy to the leading contribution in the large N expansion of the ABJM index on H 2 × S 1 . In section 7 we conclude with a short summary of our results and comment on interesting open and related problems. In a series of appendices we discuss more technical aspects such as, for instance, the construction of square integrable modes in appendix A.
2 Towards the index on H 2 × S 1 In this section we summarize the building blocks that will be needed in order to compute the topologically twisted index of a generic N = 2 Chern-Simons-Matter theory on H 2 ×S 1 . The zero locus will be parametrized by a continuous of color fluxes and holonomies. On H 2 , these flux BP S configurations are non-normalizable but they are part of the zero locus: the localizing term Q V , which is constructed to be semi-positive definite, will vanish at them.
First, we review the SUSY localization method to compute the partition function of 3d Chern-Simons-Matter defined over a Euclidean space M with off-shell supersymmetry charge Q. The space M is usually taken to be compact. The localization principle is well known and has been elegantly summarized and interpreted in various reviews, for example, [21,22]. However, given some of the intricacies we face, for the case we discuss, we review it here, with the goal of setting up our guiding principle, notation and to highlight some of the points on which we will make a particular emphasis.
To close this section we elaborate on the specific set of boundary conditions that we shall use for background and fluctuations. JHEP08(2017)023

SUSY localization principle
The SUSY localization method is summarized in the following steps • Select a "middle dimensional" section Γ in the space of complex fields, such as a 3d vector multiplet {A µ , σ, D, λ,λ} of your theory. The path integral defining the SUSY partition function of a classical action S cl , Z[Γ] is to be performed over Γ. The path Γ must be a consistent path of integration of S cl .
• The contour Γ intersects a set of Q -BP S configurations that will be denoted as BP S [Γ] and that is better known as: the localization locus.
• For each Γ there should exist a Q V local functional of fields whose bosonic part is semi-positive definite at Γ and vanishes at BP S[Γ].
• Given the previous conditions, the strict limit τ → ∞ can be taken in such a way that the final result for the partition function is guaranteed to be where δ (2) Q V, X (0) is the quadratic expansion of Q V about X (0) . We have omitted the integration over the space M to ease the reading, but remember it is there. Let us review the semiclassical reduction (2.1). The starting point, is to notice that the partition function Z[Γ] does not change if the initial classical action S cl is deformed by an arbitrary Q -exact deformation τ Q V provided the measure of integration in field configuration space is Q invariant and that there are not contributions from the boundary of the latter. Under the aforementioned conditions, we can choose a deformation term Q V with semi-positive definite bosonic part and thereafter perform a field redefinition X → X (0) + 1 √ τ X (1) .

Background geometry and supersymmetry
In this subsection we introduce the basic elements needed for the evaluation of the localization formula for the topologically twisted index of a generic N = 2 Chern-Simons-Matter theory. Specifically, we are interested in U (N ) k × U (N ) −k Chern-Simons theory coupled to matter in the bi-fundamental representation: ABJM [23], living in the non-compact space M = H 2 × S 1 whose metric we will represent as We shall use in this paper the following signature (−, −, −) on the 3d boundary theory. The flat space metric is η = diag(−1, −1, −1).
In the conventions used in this section, the non trivial spin connection component is The 2d space H 2 has infinite volume. When dealing with extensive quantities on H 2 we will use a cut-off at large θ and drop out the dependence on such cut-off in the very end. More precisely, this recipe has been used in the context of black hole entropy in [14,24] and, in the context of holographic computations for Wilson loops it was discussed in [25]; it amounts to defining the volume of H 2 as: As general principle, we will consider background configurations that grow asymptotically as the volume element, or slower. As for extensive quantities constructed out of such non normalizable backgrounds, we shall apply the previous regularization recipe. 1 The results of these sections allow to compute the topologically twisted index of any N = 2 Chern-Simons theory coupled to matter. As mentioned before we are interested in the particular case of ABJM . The ABJM theory is composed by two vector multiplets and four matter multiplets in the bi-fundamental of the gauge group. Specifically Chern-Simons ± k : A µ , σ, D, λ q=1 ,λ q=1 ±k , where 2 q is the charge of the corresponding field under the R−symmetry flux (2.14). 1 For example, to work with boundary objects -like the boundary action (2.52)-with finite limit in the cut-off θ0 → ∞, we follow [26]. The idea is to use coordinates θ := θ0 − θ,φ = 1 2 e θ 0 ϕ , in such a way the metric dθ 2 + sinh θ 2 dϕ 2 , whereφ is a periodic coordinate with periodβ = πe θ 0 and 0 <θ < θ0. 2 The supersymmetry transformation rules are defined over the complex conjugated of (φ,ψ,F ), which are denoted as (φ † ,ψ † ,F † ).

JHEP08(2017)023
We can represent ABJM theories by the following standard quiver diagram: ABJM theories have N = 8 superconformal symmetry for level k = 1, 2 and N = 6 for level k ≥ 3, the global symmetry that is manifest in the N = 2 notation is SU where each SU(2) acts upon the doublet whose components are denoted by the corresponding labels [23]. The usual charge assignment is are the Cartans of SU(2) 1,2/3,4 and U(1) R is our selection of R-symmetry. We will turn on fluxesñ a and fugacities∆ a for the latter global symmetry generators: J a . The action of the flux 4 b=1ñ b J b -appearing in covariant derivatives -upon the Φ a 's, will be written -for convenience -as -there is not summation over the index a in the r.h.s. -where n 1 2 =ñ 1 +ñ 2 +ñ The action of the holonomy operator e i 4 b=1∆ b J b upon the Φ a 's can be written as follows where We are interested in a specific deformation of ABJM . Part of such deformation is a classical background for the R-symmetry potential (2.14)

JHEP08(2017)023
The background (2.14) is non normalizable. However, V goes like the volume element of H 2 for large θ. The deformation (2.14) has non trivial consequences in the final result of the localization procedure. From (2.11), (2.13) and (2.14) the following relations come out 4 a=1 n a = 2, The R−symmetry background allows for the presence of a Killing spinor with Rcharge q = 1. The Killing spinor equation (KSE) being and σ a , σ b being Pauli matrices. As we are using negative signature it is important to keep in mind that σ a = −σ a . (2.18) In fact, the algebrae and actions that will be defined later on, are obtained out of the results in [7] by the appropriate change of signature, and (2.18). The most general normalized solution to the KSE (2.17), is proportional to Out of we can construct an off-shell supercharge Q . Before dealing with the construction of Q , it is convenient to perform the following field redefinition In terms of the new variables, the offshell algebra takes the following form for the vector and matter multiplet 3

JHEP08(2017)023
The action of the gauge covariant derivative beinĝ (2.23) It can be shown, that the Chern-Simons theory is annihilated by (2.21), up to a total derivative There can also be a mixed CS term whenever we have several Abelian factors: Where k ij is symmetric and i = j, in this case one similarly gets boundary pieces The discussion of the topological current in [3] is valid for any Chern-Simons theory. In the case of ABJM , the topological U(1) T global symmetry is generated by the conserved The supersymmetry completion of such term, is a particularization of the action (2.26). Such particularization, is given by picking k ij = k ji = 1 and regarding just a couple of indices (i = 1, 2). The index "1" labels a background Q -spurion vector multiplet, and the index "2" labels a U(1) dynamical vector multiplet. In such a way, we obtain the corresponding mixed supersymmetric Chern-Simons action, out of (2.26). For instance, in the case of gauge group U(N ), there is a unique dynamical U(1), and the bosonic term of the latter action is In the very end, we will fix the v.e.v of the spurion vector supermultiplet to specific Q -BP S values. 4 At this point, we must select a "middle dimensional" contour of integration in field space. Let us introduce a contour Γ consistent with the one of [3] Γ vector : Bosonic F ields = (Bosonic F ields) * , e. g. D = −(D) * . (2.29) The contour Γ vector will cross a specific family of Q -BP S configurations.

JHEP08(2017)023
where m and u are Cartan valued arbitrary constants. The u are the Coulomb moduli and parametrize the Coulomb branch of the theory. Expression (2.30) is the most general solution -single valued at the S 1 factor and without fermionic zero modes -to the BP S equation along the contour (2.29).
As for the matter multiplet we define In our case, the zero locus of matter is Finally, we define Q exact terms. The Q exact terms must be semi-positive definite along Γ, as already stressed. In the case of the vector multiplet and the choice of Γ (2.29), such a term is The bosonic and fermionic part of (2.34) are where λ 2 is the lower component of the gaugino λ = λ 1 λ 2 .
For the matter multiplet, and given the choice of Γ in (2.29) and (2.32), such a term is The bosonic and fermionic part of (2.37) are where P ∓ := 1∓σ 3 2 and V µν , W µν are the field strengths of R and flavor symmetry backgrounds, respectively.

JHEP08(2017)023
The term Q V matter B is semi-positive definite when expanded around BP S[Γ vector ] and over Γ matter . As shall be shown in due time, this last statement is implied by the requirement of square integrability over H 2 . Square integrability over H 2 , imposes bounds on the spectrum of eigenvalues of the relevant magnetic Laplacian. The aforementioned bounds imply the convergence of the Gaussian path integral X (1) e −δ (2) Chern-Simons, being a gauge theory, requires gauge fixing, which we choose to be the axial conditionÂ 3 = const. (2.40) In contradistinction to 3d pure Yang-Mills theory, in 3d Chern-Simons coupled to Yang-Mills and/or matter, the constraint (2.40), fixes the gauge degeneracy completely. In the latter theory there are 3 − 1 = 2 physical off-shell vector degrees of freedom (DoF) meanwhile in 3d pure Yang-Mills there is 3 − 2 = 1 massless vector offshell DoF. For a nice review on the canonical quantization of 3d Chern-Simons theory, see for instance [27].
To implement the gauge fixing, we use BRST method [20,28] and enlarge the vector multiplet, by adding the ghost fields (c,c,b). We enlarge the algebra (2.21), by the following transformation rules  Our BRST construction is conceptually that of Pestun [20] and it has been previously presented in the 3d case by Kapustin, Willet and Yaakov [28].

JHEP08(2017)023
Before moving on, let us comment on the potential addition of Super Yang-Mills (SYM) action to the analysis. This issue was treated in details in [7] -see section 2.3 -, so we do not repeat the analysis there given, but just comment upon the conclusions there reached. First, we need to construct a SYM action, supersymmetric with respect to the localizing supercharge. In fact, such a SYM action is a Q-exact term by itself. Consequently, applicability of localization method -with the supercharge Q -implies independence on g Y M . 5

Boundary conditions
In non-compact manifolds like H 2 × S 1 or manifolds with boundary, appropriate boundary conditions must be imposed in order to have a well defined variational -Lagrangianproblem. Once a proper classical theory has been defined, quantization is in order. Let and X (1) = δA µ , δσ, δD, δλ, δλ, δc, δc, δb, δφ, δφ, δψ, δψ, δF, δF , (2.46) be the non trivial zero locus background fields and offshell fluctuations respectively. As for the X (0) we define the following boundary condition with κ ≥ 1 2 . The value of κ defines important features of the spectrum of the associated S 1 quantum mechanics, if a sort of dimensional reduction is possible to perform in this case.
The following table sketches the relation between the boundary conditions (2.48) and the results reported in the next section: The total derivative part of the offshell variation of the Chern-Simons Lagrangian (2.24), multiplied by the volume element (2.50)
To have a well-defined variational principle, we redefine the classical action from Chern-Simons to (2.52) As already stated, provided δÂ 3 is constrained to go to zero faster enough -as θ goes to infinity and zero -, the boundary term S bdy will compensate the off shell -and particularly the gauge -variation of S CS on H 2 × S 1 . Upon the latter conditions, the action (2.52) has no boundary gauge anomaly and hence we do not need to -and shall not -consider other boundary contributions -that are frequently used as anomaly compensators -, such as WZW terms.
It is immediate to check that the supersymmetric transformation of L cl is trivial by construction: the supersymmetry variation of S bdy cancels the integration of the total derivative term (2.25), as it should.
The classical action evaluated on the zero locus is Where u · m := u i m i = 1 2 T r(u m). In our conventions h i and h j are Cartan generators in the Chevalley basis, and consequently T r[h i h j ] = 2δ ij .
The contributions proportional to cosh θ M ax cancel out, θ M ax being the large cut off in θ. The divergent terms coming from the integral over H 2 and the boundary term (2.52) cancel each other. Whenever we have contributions which diverge like the volume, we regulate them as we regulate the volume in (2.8), and boundary terms are regulated as explained in footnote in page 4.
It is convenient to write down the exponential of −(2.53)

JHEP08(2017)023
The total derivative part of the variation with respect toφ † of the bosonic localizing action of matter is Under off-shell boundary conditions (2.48) and (2.49), the integration of (2.55) gives (2.56) The term (2.56) vanishes, when evaluated in the functional space we are going to integrate over. The explanation of the latter fact, shall be given in the beginning of subsection 4.1.
Notice, that the ghosts vanish at the boundary, due to (2.48). Consequently, BRST gauge transformations do not affect the boundary.
Having established the localization locus, the next step is to compute one loop determinant contributions Z X (0) . In order to do that, we need to define an appropriate functional space to integrate upon. That, will be the scope of the next section. Thereafter, we can compute Z X (0) and use equation (2.1) to evaluate our final result for the topologically twisted index.

The spectrum on H 2 with flux
The spectrum of the Laplace operator on H 2 has a long history ( [29] and references therein). Even though this section might seem just a technical remark, the result of its analysis is very relevant in reaching our conclusions. We therefore choose to include it in the main body of the text and provide more details in an appendix.
The eigenvalue problem solved in appendix A, is related to the propagation of a scalar particle in the presence of a flux on H 2 . The hierarchy of modes to be reported in this section, can be interpreted as a series of Landau levels on H 2 , that emerge due to the presence of a flux s [30][31][32]. These alternative viewpoints deserve further attention and we are keen to pay them so in forthcoming works.
In this section, we will present the outcome of the analysis that shall be reported in appendix A. We encourage the reader looking for a detailed understanding, to go through that appendix.
The Laplacian in the presence of a flux s, coming from a potential is given by The equation that defines the functional space upon which we will compute determinants is

The discrete spectrum
First, we parametrize the eigenvalue as: The quantization conditions together with equation (3.3), define a finite (resp. infinite) dimensional space of square integrable functions on H 2 , that we denote as The explicit form of the eigenfunctions f ∆,j 3 , is defined in appendix A.1 and A.2, respectively, for the case s > 1 2 . The case s < − 1 2 can be worked out analogously. The range of j 3 is given by the relations (3.10) There are additional constraints to the value of j. Indeed, the eigenfunctions f (1,2) ∆,j 3 and henceforth the spectrum ∆ = j(j + 1) − s 2 , are invariant under the transformation j → − (j + 1) .
Thenceforth we must restrict j to be either In order to have square integrable functions For the choice j > − 1 2 , restrictions (3.9) and (3.10), imply an upper bound for j Interestingly, for there are no square integrable modes. A particular conclusion of this last statement, is the known fact that in H 2 there are no square integrable scalar modes. In the presence of flux s : |s| > 1 2 , square integrable modes emerge. In appendix A.7, it is proven that in the case |s| = 1, our square integrable eigenmodes match those well known discrete modes, of the vector Laplace-Beltrami operator in H 2 , with helicity s = ±1. This last statement, suggests to explore the possibility that our spectrum encodes the full tower of higher spin square integrable eigenmodes of the Laplace-Beltrami operator on H 2 . We hope to come back to this point in the future.
The relevant scalar product is As already mentioned, and proven in appendix A, square integrability of f is interrelated to the specific bounds on j 3 and j that were previously written.
Different states f Summarizing, the space of square integrable modes for a given s is In the next section, we will refer to the following spaces The spaces (3.17) are subspaces of (3.16).

The continuous spectrum
The continuous spectrum is a direct generalization of the spectrum reported by Higuchi and Camporesi [29] to the case where there is a constant flux on H 2 . The corresponding eigenmodes solve the defining equation

JHEP08(2017)023
and boundary conditions Conditions (3.19) and (3.20) are given in coordinates x, but they are equivalent to the particularization κ = 1 2 of (2.48). The final solution to the boundary problem just presented, is obtained by imposing (3.20) on the most general solution (A.3). The result is The norm of f ∆ λ j 3 under the scalar product (3.15) is infinite. By choosing appropriately the remaining integration constant one can set (3.23) Comments. Our thermal cycle is not the S 1 inside the H 2 , but the trivially fibered one. The latter fact, is related to an important conceptual difference between the physical framework of our approach and the one of, for instance, [33]. In physical terms, our H 2 modes are not probing the near horizon limit of a black hole in the presence of electric flux [34][35][36][37][38][39][40], but the boundary dynamics of a magnetically charged hyperbolic AdS 4 black hole.
That said, if we interpret our ϕ-cycle as the thermal one, our hierarchy of square integrable modes is certainly probing Euclidean AdS 2 , in the presence of an electric flux deformation. That is closely related to the problem addressed in [33]. To have a selfconsistent approach, coming from supersymmetric localization, to the problem studied in [33], one should try to localize an appropriate off-shell supercharge on the quantum gravity side. In that spirit, in [41], AdS 2 × S 2 was shown to be the unique ungauged BP S localizing solution, to 4d N = 2 Super-Conformal gravity in a convenient gauge fixing. It is plausible, that other electric or magnetic AdS 2 × Σ localizing solutions can be found, by relaxing some of the conditions used in [41], as suggested by the results in [13]. If that is the case, it would be quite interesting to explore what can supersymmetric localization say about the problems addressed in [33].

JHEP08(2017)023
term is where the result in the second line follows from the asymptotic behavior (2.48) and (2.49). This boundary term vanishes because, as proven in appendix A, the only modes of φ that do not vanish at the contractible cycle θ = 0, have the following angular dependence and they are annihilated byD Due to a careful choice of boundary conditions, total derivatives are irrelevant for the current discussion, as they do not contribute to the 1-loop determinant.
After integration by parts, the quadratic expansion of the bosonic part of the Lagrangian density of the matter localizing term (2.38), takes the form Notice that the operator O B is positive definite on "representations" Ξ j (s) labeled by j running at step 1 down from |s| − 1 but larger than − 1 2 . For "representations" labeled by j : , the operator ( s − s) has eigenvalues that obey the following inequality −j(j + 1) + s(s − 1) ≥ 0 and hence is semi-positive definite. Consequently, O B is positive definite if (ρ(u) + i∂ t ) 2 > 0. This last condition is guaranteed, provided we avoid points in the Coulomb branch such that ρ(u) ∈ Z.
Having in mind the particular case at some stages, we will denote the aforementioned set of j s as follows The union of the aforementioned Ξ j (s), is the maximal space of square integrable modes (3.16). For later convenience, let us define JHEP08 (2017)023 where k = i∂ t . The result for s < − 1 2 is obtained analogously. As φ is a complex scalar, the functional integration

Fermionic localizing operator
To compute the fermionic determinant, we used the square of the kinetic operator that appears in the quadratic expansion of the fermionic part (2.39) of the localizing term, we also specify the space of functions on which each component acts .
While reproducing the computations that will be reported in section 4.5, it will be convenient to use the following identity We use ζ-function regularization to compute the determinants of O B upon the functional space . We stress that in the case s > 1 2 and after cohomological cancellations, these zero modes are the only ones that contribute to the one loop super-determinant.
In order to compute the heat kernel K(0, 0) associated to the eigenspaces in question, we need to analyze the relevant case j 3 = s. In the latter case and after particularizing to j = s − 1 or −s, the square integrable modes f (1) and f (2) drastically simplify to The constant χ is determined from the normalization condition

JHEP08(2017)023
The 2 in the l.h.s. of (4.8) is the line element in coordinates x. From the value of |χ| 2 , we obtain the heat kernel at origin From (4.9) we obtain the zeta function After using ζ-function method, we obtain the desired determinant In appendix B.1 we obtain the same result, by using an alternative procedure. Notice that in (4.13) we could have also written (4.14) It is possible that such a change in the election of sign, changes the value of the partition function. From now on, we will ignore this second choice, except for specific steps where having it in mind will be useful. In the case s < 1 2 , the contribution to the super-determinant is coming from the following set of eigenfunctions which is the space of zero modes of In the case s < 1 2 , and after cohomological cancellations, these zero modes are the only modes that contribute to the one loop super-determinant.
Following the very same steps described in the previous section, we focus on the solutions obtained for j 3 = s − 1 and j = −s (or s − 1). In this case, the zero mode solutions reduce to f from where it is straigthforward to compute the ζ-function.

JHEP08(2017)023
The constant χ is determined from the normalization condition From the value of |χ| 2 , we obtain the heat kernel at origin From (4.18), we obtain the zeta function Finally, we obtain the desired determinant

Super-determinant
In the computation presented in this section, we will assume j, j 3 , s ∈ Z, the other cases can be worked out in complete analogy. The super-determinant in the case s > 1 2 is Notice that in the r.h.s. of the second line, we have a quotient of two identical infinite products. This cancellation, occurs due to the supersymmetric pairing of eigenmodes: cohomological cancellations.
Let us comment about the particular case s = 1. In that case, there are not cohomological cancellations. The reason is, that when s = 1, only the space of zero modes j = s−1 JHEP08(2017)023 (or j = −s), for the scalar φ, and the "chiral" spinor ψ + exist. In the case s = 1, there is not "anti-chiral" square integrable mode ψ − on H 2 , because for such spinors the effective flux is s − 1 = 0.
Next, let us compute the super-determinant in the case s < 1 (4.23) The final expression coincides with the one of the case s > 1 2 . However, in the case of s < 1 2 the unpaired modes are the "anti-chiral" square integrable modes ψ − , labeled by a radial number j = −s, and perceiving a flux s − 1 on H 2 .
Let us treat separately the case s = 0. In that case, there are not square integrable "chiral" modes ψ + , neither scalar ones φ. However, there exist "anti-chiral" square integrable modes ψ − , perceiving an effective magnetic flux of −1 and labeled by radial number j = 0. In this case the regularized super-determinant is given by Collecting partial results, not only for the case j, j 3 , s ∈ Z, but for the most general spectrum j, j 3 and s = − ρ(m)−q R 2 obeying "discreteness" conditions (3.6), is straigthforward to obtain the final result for the one loop super-determinant part of the index on H 2 The restriction to s = 1 2 , is a necessary condition to have square integrable modes. However, |(ρ(u) + k)| s− 1 2 = 1 and consequently s = 1 2 can very well be ignored and the result for Z matter 1−loop (H 2 , m, q R ) will be the same. Interestingly, under GNO conditions coincides with the square root of the analog result on S 2 × S 1 , under the identification of s H 2 with s S 2 , for each mode (ρ, k). Namely, under GNO quantization conditions Let us remark that we are not forced to impose GNO conditions on H 2 . Consequently, we shall not impose GNO quantization conditions The one loop determinant of a matter multiplet in the adjoint representation of the gauge group, with R-charge q R = 2 (which coincides with the vector multiplet superdeterminant, see appendix B.2) in the presence of flux is Notice that the result for the vector multiplet (4.25) -which is independent on m -, matches -up to a constant that can be absorbed in the regularization procedure -with the result of [18] in their flat space limit L → ∞ and under the transformation u → iu. Notice that their transformation of A µ under , δ A µ matches with ours, if only if A µ is substituted by iA µ .
We tried to obtain the cohomological cancellations from ζ-regularisation of the non zero modes. However, the heat kernel method for spinors of [29], definitely, does not respect supersymmetry (in the case of the discrete spectrum when the absolute value of the total flux felt by the mode must be larger than 1 2 ) and breaks the cancellations between fermions and bosons, unless the normalization associated to the heat kernel of antichiral modes ψ − is modified in a non elegant way.

What about the continuous spectrum?
So far, we have focused on the contribution of the discrete spectrum to the index. It is time to find out what is the contribution of the continuous spectrum. The eigenvalues of the relevant differential operators when acting upon the eigenfunctions of the continuous JHEP08 (2017)023 spectrum are (4.27) The super-determinant to compute is . (4.28) The determinants in (4.28) are computed by the method of Heat Kernel. Once the f ∆ (λ,s) ,j 3 are normalized as in (3.23) the Heat Kernel of an operator O with eigenvalues E[λ, s] when acting upon f ∆ (λ,s) ,j 3 is defined as where p = {θ, ϕ, t} and p labels the set of coordinates of a given point. We do not need the full heat kernel, since for the ζ-function all we need is its value at the origin p = p = 0. As in the case of square integrable modes, all the eigenmodes f ∆ (λ,s) ,j 3 vanish at the origin, except for those with j 3 = s. After some work, the spectral function µ (s) (λ) is found to be 1 vol . (4.30) Having the spectral function in (4.30) we are ready to compute the ζ-functions by using the following definition Notice that, the spectral function obeys the following property 33) and the eigenvalues of the operators O B and ((ρ(u) + k) 2 + + s−1 ) upon the respective eigenfunctions f ∆ (λ,s) ,j 3 and f ∆ (λ,s−1) ,j 3 are the same, see equations (4.27). From the latter facts and after the definitions (4.31) and (4.32) it follows the triviality of (4.28) In conclusion, the continuous spectrum provides a trivial contribution to the topologically twisted index.

GNO condition?
The gauge potential representative (3.1) is singular at the contractible cycle θ = 0. There are ways to solve this issue. One of them is to impose the holonomy of the gauge potential θ = 0 to be in the centre of the group, we shall not resort to this way. A second way is to simply perform a non trivial gauge transformation with parameter The new potential is regular at θ = 0 and has the same behavior at the boundary θ → ∞ as (3.1). In fact (4.36) is the analytic continuation of the section on the north chart of the magnetic monopole bundle on S 2 to the single chart that covers H 2 .
To appreciate the consequence of the non triviality of (3.1) at θ = 0, let us comment on its effect on matter. Out of the hierarchy of eigenmodes, the only modes that do not vanish at the contractible cycle, are those with j 3 = s, see equation (4.2). When cycling around the contractible cycle, these modes exhibit an Aharonov-Bohm phase of 2πs, (4.37) due to the non triviality of (3.1) at θ = 0. Should we impose the scalars to be periodic at θ = 0, the phase (4.37) must be an integer multiplet of 2π. In consequence, periodicity of scalars is not implying GNO quantization conditions. The GNO conditions [42] are where α is any element of the root lattice of the gauge group G. Condition (4.38) states that m is in the co-root lattice of G.
If we particularize (4.2) to ρ = α and q R = 0 then (4.38) implies As already explained, (4.39) is consistent with square integrability. However, if s ∈ Z + 1 2 , the Aharonov-Bohm phase is an odd multiple of π and the scalar modes are multivalued at θ = 0. Notice that with the smooth representative (4.36) such issue is no longer present. In this representation the only modes that do not vanish at θ = 0 are those with j 3 = 0. This is because by performing a gauge transformation to a regular gauge potential the j 3 gets substituted by j 3 − s. In fact, in this smooth representation, the Aharonov-Bohm phase at any ϕ-cycle becomes an integer multiple of 2π in virtue of our quantization conditions (3.6).
Conditions (3.6) are less restrictive than GNO conditions as they include not just (4.39), but a continuous family of flux configurations. Consequently, we must not restrict our zero locus BP S[Γ] by the GNO quantization conditions.
Notice that in the case of S 2 , the monopole bundle consists of two charts. In that case, the GNO condition comes from imposing single-valuedness of the structure group JHEP08(2017)023 transformation that relates the sections at north and south [43]. In the case of H 2 , there is not such a feature.
We must say, however, that in our case, relaxing GNO conditions has consequences on global gauge invariance -see for instance section 2.1 of [44], to appreciate a related discussion -. Let us analyze the case of ABJM . In that case, indeed, under a large gauge transformation u i → u i + 2π andũ i →ũ i + 2π, the Chern-Simons terms x ikm i /2 i andx −ikm i /2 i , change by a phase of e πikm i and e −πikm i , respectively. Those phases can be absorbed by a couple of topological U(1) T holonomies 7 of the form after the change of labels It is time to analyze the ABJM index on H 2 × S 1 . In this section, we borrow notation and strategy from section 2.1 of [2]. The final scope is to obtain the leading large-N behavior of the index, in terms of flavor fluxes and holonomies. In order to do that, we will show that the corresponding large-N Bethe Ansatz equations (BAE), are equivalent to the ones defined in [2]. In fact, the leading large-N solution presented in [2] will be a solution to our BAE, and consequently, can be used to evaluate the leading large-N result for the ABJM index on H 2 × S 1 . To obtain (4.41) we have considered the following non trivial v.e.v'sÂ T L3 = πp and A T R3 = −πp for the L and R spurion multiplets. Allover our discussion, we will fix the spurion fluxes to zero tL =tR = 0.

JHEP08(2017)023
Let us start by writing down the localization formula for the ABJM index. The Chern-Simons plus boundary term contribution is After collecting classical and 1-loop contributions, we can write down the expression of the ABJM index on H 2 × S 1 where c := 1 p,p∈Z 1 , n a are flavor fluxes and y a = e i∆a flavor holonomies. It is useful to remember that n a and ∆ a obey the relations (2.15) and (2.16), respectively. Notice that we have written back the sign degeneracy, mentioned below equation (4.13). If we suppose the former signs and the absolute values between parenthesis, become spurious in the contour of integration to be defined below, and consequently we drop them from now on. As we are interested in the large-N limit, (5.2) is enough to our purposes. The integration over fluxes and eigenvalues is dictated by the localization principle: they are the zero locus BP S[Γ] associated to our contour of field-integration Γ and their values are not fixed by boundary conditions. The color holonomies x i = e iu i ,x j = e iũ j are integrated along S 1 as follows from our reality conditions on Γ: Im[u i ] = 0 and 2π−periodicity of the integrand dependence on u andũ.
The general idea is to pick up certain residues in the large-N limit. In order to find out the position of the relevant simple poles, we will need to compute the large-N solution to the very same Bethe ansatz equation of [2]. In [2] it was suggested that in the large-N limit, the set of Bethe ansatz eigenvalues, which is the set of simple poles enclosed by our contour, condense to a support included in the region 8 if a = 3, 4 .

(5.3)
8 Even though the evidence presented in [2] is quite convincing, it would be nice to have a proof of the absence of extra eigenvalues outside the region X(∆a). We believe this is a point that deserves further understanding, but we will leave that analysis for future work.

JHEP08(2017)023
The region X(∆ a ) is the union of the regions covered by angles of the (i, j)-complex planes with coordinatex ji := e i(ũ j −u i ) . Each (i, j)-angle is defined as (5.4) To compute our residues we have two possibilities. Either we deform the (ij)-S 1 : |x ji | = 1, to the perimeter of the inner region of the (i, j)-angle including the origin, or we deform it to the perimeter of the outer region, and close the contour at infinityx ji = ∞. 9 Both choices are equivalent, provided we include the respective "boundary contribution" atx ji = 0 or x ji = ∞, depending on the case. The "inner" choice corresponds to selecting poles with Im(ũ j − u i ) > 0. The "outer" choice corresponds to poles with Im(ũ j − u i ) < 0. Notice that for the large N solution of our interest, an outer (i, j)-pole implies the presence of an inner (j, i)-pole and vice versa (see equation (2.39) of [2]). In fact in the limit N → ∞ all of such poles will condense either at 0 or ∞, except for the case i = j. In the particular case i = j, the poles condense to an arc of (i, i)-S 1 . For each (i, j), we will choose to deform the (i, j)-S 1 contour in the way that encloses the poles at the "bulk" (these, are the poles whose positions are the eigenvalues in equation (2.39) of [2]).
The projection of the domain X(∆ a ) upon our integration contour, which is obtained by demanding Im[u] = Im[ũ] = 0 on the former, is a parametrization of the maximally connected region without Coulomb branch singularities. 10 Actually, the intersection of the boundary of the complex region X(∆ a ), ∂X(∆ a ) with our contour of integration S 1 , is a parametrization of the domain of such singularities. By a singularity of the Coulomb branch we mean a point (u i ,ũ j ) such that the quantity that defines the one loop contributions becomes 0 or ∞. Notice that, as we are not imposing the GNO conditions, we have to integrate over the values of fluxes m i andm i along the Cartan directions. As we have already stated, the fluxes m andm are non-normalizable modes, even though they are in BP S[Γ]. In that respect, our approach is reminiscent of the one advocated in [45]. In [45], integration over non-normalizable modes belonging to the zero locus of the relevant supercharge, was suggested for the localization on H 2 × S 1 . Although the localization performed there, was on the branched sphere, the integration over the Coulomb branch parameter completes the nice picture suggested by (4.15) of [45]. 11 9 By inner (resp. outer) region we mean the region of the (i, j)-angle that is in (resp. out) of the (i, j)-S1. 10 In fact in integrating xi over S1 we should avoid colliding with such singularities, either by slightly deforming the contour, or by turning on an infinitesimal mass regulator in the matter one loop determinant. 11 A second possibility we will not explore in this work, is to fix the values of non normalizable modes in Q [Γ] to specific values. However, in order to match the final result to the supergravity dual, an extremization procedure should be engineered for those values. In some sense, integration over these non normalizable color -not flavor-modes is such sort of extremization. In this second, more open, line of thought, perhaps one could relax our reality condition on ui and simply define the latter extremization as

JHEP08(2017)023
We use a couple of very large cut-offs M ± ,M ± > 0, because the volume of the moduli space of fluxes (m,m) is infinite. After computing the integral over the holonomies for fixed values of M ± andM ± , we are free to send one and only one, of either M + (resp.M + ) or M − (resp.M − ), to infinity. The other one, remains as a regulator that we will redefine as M (resp.M ). Thereafter, we pick up the residues of the analytical continuation of the regulated integrand, the final result will be independent on M andM , and we are free to take M,M → ∞ on such residues. 12 Next, we shall show independence of the regulated expression on the cutoffs M and M . For clarity, it is convenient to re-organize the r.h.s. of (5.1) as follows integration over the complex and more abstract Jeffrey-Kirwan (JK) contour. It would be quite remarkable if with such an alternative approach, one obtains the same result for the index. In that case, the approach followed in this manuscript, would provide a less abstract viewpoint of the JK contour. We suspect this is indeed the case, but as it is not the final goal of our study, we shall not check so in this manuscript. 12 The consequences of using this regularization procedure are somehow reminiscent of the consequences of applying the Jeffrey-Kirwan recipe in [2].
The next step is to evaluate the residues of the analytical continuation of the regulated integrand in (5.10), at the simple poles enclosed by our contour. Such poles are located at positions x * andx * , defined by the eigenvalues of the BAE e iB i (x * , x * ) = 1, e i B i (x * , x * ) = 1. (5.11) Solutions of (5.11) are not onto our integration contour, 13 but are enclosed by it. In this way we solve the remaining integrals in (5.10). Notice that we are naively focusing on the contribution coming from simple poles generated by the BAE (5.11). Next, we shall see how the solution to (5.11) is independent on the topological holonomies p andp. The integrand of (5.10) is also independent on p andp. Consequently, the average over p and p will be trivial. The final result for the index is Notice that the cut-off dependence disappeared in (5.12) due to imposition of the BAE. We are interested in the large-N limit of Z H 2 ×S 1 . In that limit, there is a unique BAE solution, thence the summation over the label I becomes spurious.
It is convenient to define (5.14)

JHEP08(2017)023
In terms of D(z), the l.h.s. of the BAE are In terms of the quantity the matrix B takes the form The next step is to write down the BAE in "angular" coordinates u i ,ũ i , and ∆ a , which are defined from In these coordinates, the constraint a y a = 1 looks like a ∆ a = 0 (mod 2π).
To fix the values of n i andñ j , we use the identity together with the assumption of absence of "long range interaction" [2]. The latter condition, implies

JHEP08(2017)023
It is key to observe that (5.17) and (5.18), imply that solutions to the BAE are independent of the topological holonomies p,p. Additionally, (5.15) is also independent of p andp and we conclude that the average on p andp can be substituted by one. The solution of (5.11), obeying (5.19), is precisely the one used in [2], to obtain their quite nice result, (2.89), in the large-N limit. In the next section, for completeness, we evaluate the aforementioned solution. We have tried other choices, such as a ∆ a = 0. However, as pointed out in [2], there are always some issues with the potential solutions.

Large-N behavior of the index
From now on, we take the limit N → ∞, assume the Chern-Simons level k = 1, introduce the density ρ(t) := 1 N di dt and the quantity δv(t), as precisely done in section 2.3 of [2]. In this continuous limit, the BAE arise from the variations of the auxiliary Lagrangian It is easy to follow the steps in [2]. Indeed for a ∆ a = 2π and under the assumptions

From (5.22) it follows the ordering of transition times
From the normalization condition To obtain the leading free energy in the limit N → ∞, one evaluates (5.12) at the BAE solution (5.24). The final result can be easily inferred, given the fact that our BAE solution is the same one found in [2], for the case k = 1, a ∆ a = 2π. In the latter case, and from the fact that the summand in (5.12) is the square root of the one in eq. (2.24) of [2], it results that Re log Z S 2 ×S 1 | Large N BAE solution + sub. terms Finally, one can arrive to the result After extremizing with respect to ∆ 1 , ∆ 2 and ∆ 3 , we obtain the following relation between fluxes and holonomies

JHEP08(2017)023
that will prove to be useful later on, when comparing with the conjectured AdS/CF T dual quantity.

Comments on the index
The index computed in this work, which follows closely [3], has the canonical interpretation of a Witten index counting ground states according to Z(n a , ∆ a ) = Tr(−1) F e −βH e iJa∆a . From the 3d perspective, we are simply counting operators with the corresponding relation among quantum numbers. Now, assuming that the deformed ABJM theory flows to an effective quantum mechanics in the IR, then the index computes the degeneracy of ground states in the quantum mechanics. Since the index is an invariant of the flow, we connect directly the 3d and 1d perspective.
On the gravity side, we have, similarly, the possibility of viewing the counting from the 4d or 2d perspectives. The better formulated one, at the moment, turns out to be the 2d perspective, which Sen has developed in the framework of AdS 2 /CF T 1 [14,26]. In this context, the ground state degeneracy in the quantum mechanics is computed by an AdS 2 partition function with specific boundary conditions, which leads precisely to the quantum black hole entropy. There are, of course, some open issues with the application of Sen's proposal, in the context of asymptotically AdS black holes, but it certainly provides a solid starting point.

The hyperbolic AdS 4 black hole
In this section we construct what we believe are the holographic dual to the ABJM configuration discussed thus far. Namely, we construct magnetically charged, asymptotically AdS 4 black holes with non-compact H 2 horizon that are embedded in M-theory. Our construction follows similar spherical solutions in N = 2 gauged supergravity, see, for example, [9][10][11][12][13]. We will comment on some similarities and differences with the solutions with spherical horizon in subsection 6.3. We shall focus on the case of n V = 3 vector multiplets. In this way the n V + 1 vector fields -counting also the one in the graviphoton vector multiplet -are set to be identified as dual to the global charges of ABJM .
6.1 A brief summary of 4d N = 2 gauged SUGRA with n V = 3 For completeness, let us briefly introduce the concepts that we shall use. The central object is the pre-potential which is a holomorphic function of the holomorphic sections X ∆ (z i ), ∆ = 1, 2, 3 and 4. The symplectic sections are functions of the physical scalars z i with i = 1, 2, 3. Another important object is the Kähler potential We will also need to use the period matrix

JHEP08(2017)023
and the following auxiliary variables where the covariant derivative D i is defined as D i := ∂ z i + K i . In our case we will be interested in real holomorphic sections To construct black holes, we shall set the fermions, and fermionic variations, to zero. The supersymmetry variation of the gravitino is where the covariant derivative of the Killing spinor ε A is The supersymmetry variation of the gaugino will be used too. In equations (6.4) and (6.5) we are discarding higher order terms in fermions. These terms are not relevant to our discussion. To reproduce our results it will be useful to have the following definitions [46] F − Λ :=

Hyperbolic black holes
To avoid confusion, the index Λ = {1, 2, 3, 4} is equivalent to the index a = {1, 2, 3, 4} that will be introduced in the next subsection. We are interested in the STU model. Thence, we fix the Fayet-Iliopoulus parameters in an isotropic manner The relevant pre-potential will be We consider real sections, with the following parametrization JHEP08(2017)023 and propose the following static, spherically symmetric ansatz, for the metric and sections ds 2 = −U −2 (r)dr 2 − h 2 (r)dθ 2 − h 2 (r) sinh 2 θ dϕ 2 + U 2 (r)dt 2 , (6.8) The non trivial components of the spin connection are 14 In this section 4123 = 1, η ab = (−1, −1, −1, 1), We use the following Ansätze for the functions U (r) and h(r) where g, c and d are constants. The corresponding black holes, are sourced by magnetic fluxes p Λ : The non trivial components of the anti-selfdual field strength are The chiral and anti-chiral Killing spinors A and B , have to obey the following relation -these conditions are obtained from the vanishing of the gravitino supersymmetric transformation- 14) The most general solution to (6.14) is (6.15) 14 In this section we used different conventions than in section 2. We have used standard conventions on four-dimensional N = 2 gauged supergravity, which are the ones given in [46] (see also [11]). For instance, the definition of spin connection is minus the one used in section 2. Consequently, in the covariant derivatives there is a relative minus sign in front of the term proportional to the spin-connection between this section and section 2.

JHEP08(2017)023
Solving the BP S conditions, leads to relations Notice that the constant c is positive. The relation between fluxes and the parameters β a is also obtained from the BP S conditions The warping of the Killing spinor is also fixed by the BP S conditions We have, finally, completely solved the BP S conditions and constructed our hyperbolic AdS 4 black holes.

Spherical black holes
A prevalent intuition in the context of supergravity states that changing the horizon from spherical to hyperbolic, leads from black holes to naked singularities and vice versa [47,48].
In this section we explore the details of this intuition in the context of the magnetically charged black holes. Let us first solve the BP S equations for the spherical black hole ansatz ds 2 = −U −2 (r)dr 2 − h 2 (r)dθ 2 − h 2 (r) sin 2 θ dϕ 2 + U 2 (r)dt 2 , (6.21) with U (r) and h(r) defined in (6.10) and (6.11). The non vanishing components of the spin connection are For technical convenience, let us parametrize the gauge potential as follows 15 The definition of field strength used in this section differs from the one used in section 2.

JHEP08(2017)023
The non trivial components of the antiselfdual potential are . (6.24) After solving the BP S equations, we arrive to and The warping of Killing spinor is also fixed by the BP S conditions In contradistinction to the hyperbolic solution, in the spherical case, c can be negative, see equation (6.26). The position of the curvature singularity is If β a = 0, the position of the curvature singularity is r s = 0. In the case β a = 0, one encounters a hyperbolic black hole since c > 0 and there is a horizon. However for the spherical solutions, c < 0, and we encounter a naked singularity. It is straightforward to check that, when β a = 0, the change (r, θ, t, p a , c) → (ir, iθ, it, − p a , − c) (6.30) transforms the hyperbolic black hole BP S solutions of the previous subsection, into the spherical BP S solution of this subsection, particularized to β a = 0. The latter has a naked singularity. Actually, the exchange p a → −p a can be cancelled by an exchange of Killing spinor -only because β a = 0 -. By an exchange of Killing spinor, we mean to change the choice of sign in the constraint (6.14). Such a change, has physical meaning, as it leads to a configuration that is BP S with respect to a different supercharge. We emphasize that the intuition emanating from [47,48], is restricted to the case of β a = 0, that is, the case of constant sections; it is, therefore, relaxed in the case β a = 0. In this subsection, we compare the entropy of hyperbolic and spherical black holes. We first consider the case of isotropic fluxes. In the end, we will find that the entropy density of the hyperbolic solution coincides with the entropy density of the spherical one which were discussed in [2].
From now on, we particularize our hyperbolic solutions to the following case A sufficient -not necessary -condition for the existence of hyperbolic AdS 4 black holes is where the domain of the radial coordinate r is r > r s . The constant r s is the radial position of the singularity, which is covered by the horizon at r h = √ c. A particular solution to these conditions is In that case, the fluxes are The regularized area density of the horizon is Next, we compare the entropy density (6.34) to the one of spherical black holes used in [2]. The isotropic solution presented in section 4.1 of [2] is From the following quantities [2]

JHEP08(2017)023
one arrives to the following expression for the area density, of the spherical horizon in terms of the flux p We can do the comparison with the result obtained from scratch, with our spherical solutions, however, in order to match our results with the ones in [2], we report the comparison by using theirs. Notice that the large flux limits (6.36) and (6.39), do coincide. In fact, it can be checked that for any value of u the entropy density as a function of flux p and p coincide, as Equation (6.40) can be checked by comparing (6.34) and (6.38), order by order in Taylor expansions in u about 0 and ∞, or by simply working out the expressions.

Matching results
In this final section, we compare the AdS/CF T dual results. On one side, we have the result for the ABJM index on H 2 × S 1 (5.25). On the other side, we have the entropy of the hyperbolic magnetic AdS 4 black holes (6.8). The first thing to do, is to compute the Bekenstein-Hawking entropy of (6.8). Thereafter, we check the relation between the classical entropy and the value of the holomorphic sections X a -or as were denoted in the previous subsection X Λ -at the horizon r h . The aforementioned relation is identical, up to a relabelling of variables, to the relation between the logarithm of the ABJM index and the holonomies ∆ a (5.25). Finally, we prove that under the appropriate relabeling of variables and extremization of the logarithm of the ABJM index (5.25) with respect to the holonomies ∆ a , the bulk and SCFT results coincide, as was the case in [2]. For generic p 1 , p 2 and p 3 , we have checked that the classical entropy where Ψ(β 1 , β 2 , β 3 ) := 8d 2 β 1 2 + β 1 β 2 + β 1 β 3 + β 2 β 3 + β 2 2 + β 3 2 + 1, (6.42) coincides with the expression The value of β 4 is determined by (6.18). The p a 's as function of β a 's have been given in equation (6.20) which follows from the BP S equations. We will prove next, that equation (6.43) -that comes from the analysis in the bulkis equal to the extremal value of the SCFT topologically twisted index (5.25), under the specific dictionary √ 2p a ↔ n a , (6.45) −2πX a (r h ) ↔∆ a , a = 1, 2, 3, 4, (6.46) where the∆ a are the solutions to the variables ∆ a that come out of the inversion of equation (5.27).
There are many ways to prove that (6.43) is equivalent to the extreme value of (5.25), the simplest one is to evaluate (5.27) on ∆ a = ∆ a (β 1 , β 2 , β 3 ) = −2πX a (r h ), (6.47) to obtain Notice that Equation (6.48) = √ 2 Equation (6.20), (6.49) under (6.31) and identification (6.45). The relation (6.49) implies that the positions∆ a of the saddle points of (5.25) coincide with the values of the sections X a at the horizon (6.44), under identification (6.45). As the logarithm of the ABJM index (5.25) and the Bekenstein-Hawking entropy (6.43) are the same under identification (6.45) and we have thence proven that under the aforementioned identifications, boundary "degeneracy of states" and bulk black hole entropy coincide. As a final comment, we notice that the identifications (6.45) are not directly obtained from the AdS/CF T dictionary. The AdS/CF T dictionary is naturally formulated in the UV, the UV value of the holomorphic sections X a is − 1 4 . To obtain agreement, the use of extremization principle of the result of the SCFT side is crucial [2]. We believe there is a proper way to clarify some of these ad hoc issues but we leave the discussion for future work.

JHEP08(2017)023 7 Conclusions
In this manuscript we have first studied topologically twisted localization of N = 2 supersymmetric field theories in H 2 × S 1 . Our work differs in various important points from recent work on localization on this space by [18]. In particular, we have crucially considered topologically twisted theories and extended the type of theories under consideration beyond vector multiplets, to include, for example, matter multiplets.
At a technical level, we have also discussed explicitly subtle aspects of the eigenvalue problem corresponding to the Laplacian in the presence of a background magnetic field and we expect that such results could have wide application in the general context of localization. Quite interestingly, we have found a hierarchy of normalizable modes and its corresponding discrete spectrum. A particular sub-family of the aforementioned hierarchy, corresponds to the vector zero modes of the Laplace-Beltrami operator on H 2 , that were introduced in [29] and figure prominently in, for example, [49] and more recently in [18]. The full hierarchy of normalizable modes exists due to the presence of magnetic fluxes s over a specific threshold: |s| > 1 2 . We strongly suspect, that the discrete spectrum is encoding the full hierarchy of higher spin normalizable modes of the Laplace-Beltrami operator on H 2 . If this is indeed the case, it would be very interesting to pursue a study of the potential traces of 2d higher spin symmetry, on the set of black holes microstates [50]. As a first step toward formulating such a problem, it would be useful to start by identifying the square integrable modes in the language of [50].
We have also studied N = 2 gauged supergravity and found magnetically charged supersymmetric solutions with hyperbolic horizon. We have shown that under assumptions similar to those advanced in [2] the entropy of these solutions coincides with the real part of the logarithm of the topologically twisted index of the dual field theory. In conclusion, we have provided evidence in favor of identifying the set of square integrable modes in the presence of a constant flux on H 2 ×S 1 , precisely speaking a very restricted set of zero modes out of the maximal set, as the boundary microstates responsible for the Bekenstein-Hawking entropy of the AdS 4 hyperbolic black holes presented in section 6. One important further test for this identification, would be to compute quantum corrections to the Bekenstein-Hawking entropy, on both sides of the duality.
On the 3d SCFT side we have made crucial use of the extremization approach advocated in [2]. The result of this approach is consistent with the constraints from the BP S equations on the gravity side and has been argued to be equivalent to the attractor mechanism. Under these conditions we have found precise agreement between the leading large-N results on the two sides. However, it would be important to elucidate the role of extremization intrinsically in the field theory but also from the gravity perspective. This is particularly important because in some cases the attractor mechanism has been shown to apply away from the strictly supersymmetric context.
Another natural generalization of this work, following [51], is to extend the analysis to dyonic black holes. More generally, it would be interesting to consider mapping the full space of deformations on both sides of the correspondence and, in particular, its modifications on the free energy and the entropy. Another interesting direction concerns potential JHEP08(2017)023 factorizations of the index on S 2 × S 1 , introduced in [3], in terms of blocks given by the partition functions in H 2 × S 1 . A similar factorization principle has been uncovered in various theories and in different dimensions, see, for example, [7,[52][53][54][55]. In this manuscript we have found a particular relation but it should be pointed out that we have set all fermionic zero modes to zero and have integrated over a particular set of modes. Clearly, to achieve a bona fide factorization formula we will need to consider more general boundary conditions and contemplate retracing some of the steps suggested in [45]. Indeed, such an approach with general boundary conditions has been implemented for GLSM's in [54].
Finally, it would be interesting to discuss the microstate counting of magnetically charged strings in asymptotically AdS 5 spacetimes. Such magnetically charged solutions have a long history in supergravity dating back to explorations in [47]. It is logical to expect that the microscopic explanation should be found within 4d topologically twisted field theories on S 2 × T 2 or possibly H 2 × T 2 . Indeed, as a natural starting point along these lines, the topologically twisted index introduced by Benini and Zaffaroni in [3] for supersymmetric field theories on S 2 × S 1 , was briefly discussed for 4d theories in S 2 × T 2 in their original work, and was also addressed in [56]. It has recently been shown that, in the high temperature limit, the index produces a central charge that matches the supergravity answer [57] therefore providing a strong argument in favor of the identification. We hope to report on some of these interesting directions soon.

A Comments on the discrete spectrum
In this appendix we report details on the construction of the square integrable modes defined in section 3.
The general solution to the defining equation (3.3) is

JHEP08(2017)023
where χ 1,2 are arbitrary integration constants. It is very useful to work in the following coordinates In x coordinate, (A.1) takes the form where the parameters a 1,2 , b 1,2 and c 1,2 are Let us particularize ∆ to The choice (A.5) completes perfect square inside the square roots in (A.4). The asymptotic behavior of the solutions proportional to χ 1 and χ 2 -from now on f 1 and f 2 -, is and Regularity at the contractible cycle x = 0, conditions to pick up We demand C ∞ -differentiability at the contractible cycle x = 0. In the vicinity of Notice that from (A.12), it follows that: j 3 ∈ Z implies s ∈ Z.. However, we are not forced to impose integrality of j 3 , j or s. For the time being, let us assume s ≥ 0. In due time, we extend the analysis to the case of generic s.
It is important to stress that square integrability condition is equivalent to impose the conditions For pedagogical reasons, let us assume for the time being We shall see in due time that assumption (A.15) is not necessary. As for j, let us not assume nothing else at this point. On the track, we will comment on the restrictions that arise for j. Condition (A.14) selects the solution f 1 16) Notice that this solution is invariant under the transformation j → − (j + 1) and in consequence we have to restrict j to be either as preferred. Notice that j = − 1 2 is left invariant by the transformation above. For j = − 1 2 both independent solutions have the same asymptotic behavior x j and x −j−1 , and they are not square integrable.
In order to have square integrability it is necessary to have We exclude the particular case j = − 1 2 .

JHEP08(2017)023
Before writing down the quantization conditions, let us comment on the strategy. It turns out that the quantization conditions are the conditions for which the hypergeometric factor in (A.16) truncates to a specific polynomial. The sum of the degree of such polynomial with the degree of the leading power of the prefactor in (A. 16) in the limit x → −∞ must equate to The conditions to obtain the previously mentioned goal are Together with (A.14) these conditions are compactly written in the following one It is straightforward to check that, if we assume (A.15) together with (A.22) the desired truncation holds: The degree of the polynomial d (1) being At this point is easy to check that indeed the aforementioned asymptotic behavior of (A. 16) about x = 0 and x = −∞ holds.
As the previous case, this solution is invariant under the transformation j → − (j + 1) and in consequence at some stage we shall be forced to assume (A.17). Let us, however, not assume the latter restriction on j yet. Let us just assume (A.15). The quantization conditions are Together with (A.14), these conditions are compacted in the following one It is straightforward to check that, if we assume (A.15) together with (A.27) the desired truncation holds: The degree of the polynomial d (2) being ∆(s), j 3 is square integrable. Notice that, there are not square integrable modes for s = 0, as well known.
A.3 The case of negative flux s < 0 So far, we have focused on the case of positive magnetic flux s > 0, or being more specific on the case s > 1 2 . However, there are square integrable modes when s < 0 too -as parity preservation dictates -. To find those, it is convenient to use the identity upon the previously written solutions f ∆(s), j 3 , to obtain ∆(s), j 3 = χ 2 x −s (x (x − 1)) j 3 +s 2 2 F 1 (j 3 − j, 1 + j 3 + j, 1 + j 3 − s; x) .

JHEP08(2017)023
Again, these eigenfunctions are invariant under the change (A.17) and in consequence The hypergeometric factors written above, truncate to polynomials -and in consequence f

A.4 Generalized conditions
So far we have been assuming However, the aforementioned GNO conditions -see subsection 4.7 -can be relaxed. As already stated, to achieve regularity at the contractible cycle x = 0, the following necessary condition j 3 − s ∈ Z, must hold. To have discrete spectrum there are necessary conditions too: Notice that the conditions in the right (resp. left) side, follow from a linear combination of the condition of regularity at the contractible cycle, and the respective conditions in the left (resp. right) side. Hence, we can write down the more compact and equivalent statement In the table below, we write down the explicit form of the spectrum. For simplicity of presentation but without lack of generality, let us take j > − 1 2 . In that case, the relevant spectrum is ∀s such that s >

JHEP08(2017)023
A particular case is when j, j 3 , s ∈ Z + 1 2 . In that case, the table above reduces to ∀s such that s > . The corresponding eigenfunctions can be recovered from the summary that shall be presented next, and the results in previous sections.

A.5 Collecting the eigenfunctions
The maximal functional space of square integrable modes is Ξ(s) := Of special interest will be the following limiting spaces These spaces are the ones that contribute to the super-determinant that concerns us, when s > 1 2 and cohomological cancellations are performed. It will be useful to keep in mind that for every j, in the direct sum space Ξ j (s) := Ξ  j (s) the angular number j 3 will range at step 1, departing from the lower (resp. upper) bound given below max(|j|, |j + 1|) < j 3 < ∞ if s > 1 2 , −∞ < j 3 < − max(|j|, |j + 1|) if s < − 1 2 .

JHEP08(2017)023
Comment. Notice that upon square integrable "representations" Ξ j (s), which are labeled by j running at step 1 down from |s| − 1 and greater than − 1 2 , namely j : − 1 2 < j ≤ |s| − 1, the bosonic operator O B (4.4), is positive definite if (ρ(u) + k) 2 > 0. Indeed, that operator needs to be positive definite in order to have convergence of the functional integral of the exponential of the quadratic expansion of the bosonic localizing term.

A.6 Normalizable modes from asymptotics
One can also find the discrete spectrum by looking at the asymptotic expansion of the general solutions (A.16) (one can repeat the procedure for the other solution). We choose to focus on f 1 and for values s > 1 2 , from regularity and smoothness at x = 0 it follows that j 3 ≤ s and the difference s − j 3 ∈ Z + .
As before, we define ∆ := j(j + 1) − s 2 . At x = −∞ f 1 is The coefficients above are: Suppose j > − 1 2 , to cancel out the x j behavior of f 1 while preserving the x −1−j , we need to have that χ + 1 vanishes, this is achieved when either of the arguments of the Γ's in the denominator is 0 or a negative integer. Then 1 − j 3 + j = −n or s + j = −n with n ∈ Z + , (A. 35) the second choice is out of order given our assumptions (s > 1 2 , j > − 1 2 ) therefore 1−j 3 +j = −n when replacing this value in χ − 1 one has to be careful, since the arguments of Γ's in the denominator might be also a negative integer, We have not replaced the value of j in terms of n and j 3 . The Γ functions in the denominator, given our assumptions, will be a number different from zero and finite and will not play a role. Looking at the denominator in the last expression one naively conclude that there are values of j 3 and n for which the argument is negative integer (since j 3 > n + 1), and therefore χ − 1 is also zero, but it is not the case since, for each of this values, the argument JHEP08(2017)023 In words, − → X 0 is an eigen-tensor of rank one of the Laplace-Beltrami operator ∇ µ ∇ µ with eigenvalue 1.
Notice that the operator in the l.h.s. equation above coincides with our s=1 if and only if j 3 > 0.
In fact, the equation above for X 0θ implies that X 0θ obeys our defining equation if and only if 0 = j (or − 1 = j) < s = 1 ≤ j 3 ∈ Z.

B.2 Vector multiplet
In this appendix, we prove that the index of a vector multiplet coincides with the index of matter multiplet with R-charge q R = 2.
The quadratic actions coming out of the localizing terms (2.35) and (2.36), along the complex path (2.29) and after imposing gauge fixing condition (2.40) The integration of δλ 2 and δλ † 2 , multiplied by the integration of δc † and δc, following from the BRST action (2.44), gives The result (B.9) is a divergent quantity. To regularize these objects we use the zetaregularization procedure but only after co-homological cancellations are performed.

JHEP08(2017)023
If the reader is interested in the file, please write an email to us. If there is interest, we are more than happy to share it. The 4d gamma matrices