Black hole entropy function for toric theories via Bethe Ansatz

We evaluate the large-$N$ behavior of the superconformal indices of toric quiver gauge theories, and use it to find the entropy functions of the dual electrically charged rotating $\mathrm{AdS}_5$ black holes. To this end, we employ the recently proposed Bethe Ansatz method, and find a certain set of solutions to the Bethe Ansatz Equations of toric theories. This, in turn, allows us to compute the large-$N$ behavior of the index for these theories, including the infinite families $Y^{pq}$, $X^{pq}$ and $L^{pqr}$ of quiver gauge theories. Our results are in perfect agreement with the predictions made recently using the Cardy-like limit of the superconformal index. We also explore the index structure in the space of chemical potentials and describe the pattern of Stokes lines arising in the conifold theory case.


Introduction
The pursuit of a theory of quantum gravity is one of the greatest challenges of modern theoretical physics. The description of the microscopic origin of black hole entropy constitutes a highly important step towards such a theory. One of the tremendous successes in this direction was achieved in string theory which provides a microscopic explanation for the entropy of a class of asymptotically flat black holes [1].
Another source that might shed light on this question could be the study of the entropy of asymptotically AdS black holes. This is by virtue of the AdS/CFT correspondence which establishes a connection between a quantum theory of gravity and a conformal field theory (CFT) living on the boundary of the AdS-space. This connection provides us with the natural playground for the study of black hole microstates. In this setting the problem of microstate counting of asymptotically AdS black holes translates into the counting of BPS states in the dual CFT. For a long time, many approaches employed in the attempt to match these dual descriptions had fallen short. However, recent years have seen substantial progress in this direction.
While considerable progress was made in the context of static magnetically charged black holes (in the AdS/CFT setting), the situation as far as rotating electrically charged black holes are concerned had remained unclear. The simplest and canonical example of this kind is that of BPS black holes arising in type IIB supergravity on AdS 5 × S 5 [21][22][23][24]. On the dual side, microstate counting of these black holes should correspond to the counting of 1 16 -BPS states in 4d N = 4 SU(N ) Super Yang-Mills (SYM) theory on S 3 × R. The latter are counted by the superconformal index [25][26][27]. Yet, an attempt to match these two quantities has failed already at the leading order in N , since the number of BPS states grew as N 0 , which is much slower than the growth of N 2 corresponding to the black hole entropy. This slow growth was explained by the huge cancellation between fermionic and bosonic states counted by the index. Later, indices of other theories including the conifold theory [28], the Y pq family [29] and finally toric theories in general [30], were also evaluated at the large-N limit. In each case the situation repeated the one of N = 4 SYM case: the indices appeared to be of order O(1) at large N and even the matching with a supergravity index not involving any black hole microstates took place. This puzzle has been recently resolved in several steps. The first inspiration came from the extremization principle used to resolve the problem in the case of magnetically charged black holes. In particular, in [31] the entropy function of the rotating AdS 5 black hole with two angular momenta J 1,2 and three charges Q 1,2,3 was proposed. The entropy function is the Legendre transform of the black hole entropy. It is a function of the chemical potentials for the electric charges ∆ 1,2,3 and angular momenta ω 1,2 , correspondingly. The chemical potentials, in turn, appear to be constrained by the relation ∆ 1 +∆ 2 +∆ 3 = ω 1 +ω 2 −1. Later on, similar constructions were proposed for electrically-charged black holes in other dimensions [32,33]. Furthermore, it was also realized that the black hole entropy function may be reproduced by the evaluation of the complexified supergravity on-shell action in the bulk [34].
On the dual CFT side, the entropy function has a clear formulation in terms of the corresponding field theory parameters, but a precise relation to the counting of BPS states had remained obscure. An answer to this question appeared in [35] where the authors considered the superconformal index with complex fugacities. The latter allowed them to avoid large cancellations between different states and to obtain an N 2 scaling for the growth in the number of BPS states. More importantly, considering the Cardy-like limit [36], also referred to as the high-temperature limit, the authors were able to identify the superconformal index of N = 4 SYM given by with the corresponding entropy function of the electrically charged AdS black hole. Here ∆ a and ω a are chemical potentials conjugated to three R-charges of SO(6) R-symmetry and two angular momenta J 1,2 on S 3 , correspondingly. Just as in the entropy function, these chemical potentials satisfy the constraint mentioned above.
This important result was subsequently extended to various 4d N = 1 SCFTs [37][38][39][40][41]. For the purposes of the present paper, the most relevant is the observation made in [41], where the authors considered a Cardy-like limit of the superconformal index for the class of toric quiver gauge theories. The following expression was obtained at the large N limit where ∆ I are chemical potentials for (d − 1) U(1) I global symmetries satisfying the constraint d−1 I=1 ∆ I − ω 1 − ω 2 = 1. The coefficients C IJK are associated with the triangle anomalies of the corresponding U(1) global symmetries that can be obtained directly from the toric data of the gauge theory quiver. Interestingly, expression (1.2) was already proposed for the entropy function of AdS 5 black holes with d electric charges in Appendix A of [32].
Despite the great success in computing the superconformal index at the Cardy-like limit, the importance of that particular limit and the reason why it precisely reporduced the black hole entropy, even at finite temperatures, were unclear. The proper way to treat the problem is to calculate the large N limit of the index with complex fugacities, relying upon no further restrictions or limits. This kind of computation was performed in [42] for the canonical case of N = 4 SYM, for which the authors exploited the Bethe Ansatz technique [43,44]. Its main idea is to reformulate the integral representation of the superconformal index as a sum of residues of the integrand. The positions of the integrand's poles, in turn, are defined by the solutions to the Bethe Ansatz Equations (BAEs) system. Generally, it is very hard to solve this system of equations. However, in [42] it was shown that using a particular class of relatively simple solutions to the BAEs, it is possible to reproduce the entropy function (1.1) at large N , but finite temperature limit, within a certain range of values for the chemical potentials ∆ I .
Another intriguing observation made in [42] is the highly complicated structure of the large N index in the space of complex chemical potentials. In particular, there are many competing exponential contributions to the index. The full space of chemical potentials divides into separate regions with one of the exponential contributions dominating in each of them. These regions are separated by codimension one surfaces, Stokes lines, on top of which different exponents contribute the same and hence there is no dominant contribution. This also explains why the calculation in [26] resulted in an O(N ) index since it has been performed precisely on such Stokes lines where O(N 2 ) saddles cancel each other.
The BAE approach is considerably more involved than Cardy limit calculations, but it provides us with a large N answer that is precise in temperature. It further grants access to the Stokes phenomenon and the full large N index structure. It is therefore natural to apply this approach to theories beyond N = 4 SYM. Here, similarly to [41] we consider the superconformal index of various toric theories. In particular, we show that the basic solutions to the BAEs for N = 4 SYM used in [42] also solve the BAEs for all toric quiver theories. This allows us to readily evaluate the leading order contribution of the basic solution to the superconformal index of these theories. Note that in principle there might be other solutions contributing at the same order of N which are difficult to find due to the complicated structure of the BAEs. However, just as in the N = 4 SYM case, we find that the Cardy limit prediction of [41] presented in (1.2) is exactly reproduced by the basic solution contribution, and is precise in temperature for a certain range of values for the chamical potentials. We thus provide a more definite support in favor of the conjectured form of the multi-charge black hole entropy function proposed in [32].
Unfortunately, entropy functions (1.2) appear to be quite complicated for performing extremization. In particular, in [41] the only cases considered in detail were those of Y pp and L aba . However, the BAE approach also grants us access to the Stokes phenomenon. Since reproducing the full structure of the large-N index in the space of chemical potentials seems too complicated a task, we consider certain limits. In particular, as was done in [42] we consider the special case of equal chemical potentials ∆ I and also take ω 1 = ω 2 = τ as required by the BAE technique. Consequently, the only parameter on which the large-N index depends is τ , and we can define the positions of Stokes lines and regions of dominance of certain exponents. Interestingly, for the conifold theory we obtain the structure of the index in the complex plane of τ exactly reproducing the one presented in [42] for the case of N = 4 SYM.
The paper is organized as follows. In Section 2 we briefly review the BAE method for the calculation of the index and present the basic solutions to the BAEs for toric theories. We also present a general formula that can be used to evaluate the contribution of these solutions to the index at the large-N limit. Next, in Section 3 we study in detail the simplest possible case of all toric theories, namely the conifold theory. We exploit the relative simplicity that characterizes this case, to clearly illustrate how most of the arguments presented in Section 2 work. In Section 4, omitting details we directly apply the results of Section 2 to many different cases of toric quiver theories, including the infinite families Y pq , X pq , L pqr . Finally, in the last Section 5 we consider the large-N index of the conifold theory derived in Section 3 and study its structure in the case where all the chemical potentials ∆ I are equal. In particular, we find the positions of Stokes lines and regions of dominance for various contributions to the index.
Note added: When we were finalizing our paper the preprint [45], which has a significant overlap with our work, has been posted on arXiv. The authors' conclusions disagree with ours at various points.

Large N index for toric theories
In this section we will consider some generalities of index computations that we will later use in order to obtain large-N limits for the indices of toric theories. Let us consider general N = 1 quiver gauge theories with gauge group G, flavour symmetry G F and non-anomalous U (1) R -symmetry. The matter consists of n f chiral multiplets in the representation R a of the gauge group G. Since we are interested in toric theories, we will concentrate only on the cases of adjoint and bifundamental matter leaving the discussion of matter in the fundamental and other representations for the future. The integral representation of the superconformal index is then given by [27,46] Here p and q are complex fugacities for the angular momentum, v a are flavor fugacities with a = 1, . . . , rk(G F ). The integration variables z i parametrize the maximal torus of G with the integration contour given by rk(G) unit circles. The roots of G are parametrized by vectors α while ρ a are the weights of the representations R a in which the chiral multiplets transform. Also, |W G | is the order of the Weyl group. Finally, Γ (z; p, q) and (a; q) ∞ are standard elliptic Gamma function [47] and q-Pochhammer symbol: Throughout the paper, we will also be using the chemical potentials τ, σ, ξ and complexified holonomies u i as follows p = e 2πiτ , q = e 2πiσ , v a = e 2πiξa , z i = e 2πiu i . (2.4) All the functions in the index (2.1), and hence index itself, are well-defined for |p| < 1 and |q| < 1. We will be interested in the large N behavior of the index. In order to estimate it we will employ the technique of Bethe Ansatz Equations (BAE) [43,44], which has been recently used to find the large N behavior of the superconformal index of N = 4 SYM [42]. The obtained results in this case perfectly matched the entropy of the corresponding AdS 5 black hole. Below, we briefly describe this method closely following [42,43].
The main observation behind the method is that for τ and σ satisfying we can recast the integral representation (2.1) of the index as the sum over poles located at the solutions to certain transcendental equations called BAEs. Notice that as was shown in [43], the set of τ and σ satisfying (2.5) is dense in the domain {|p| < 1, |q| < 1} so the method is, in principle, applicable for any fugacities p and q. However, even when τ and σ are generic parameters satisfying (2.5), it is hard to perform BAE calculations. Therefore, for simplicity and without loss of generality we will consider the simplest possible case of τ = σ. Under this condition, as was shown in [43,44], the poles of the integral (2.1) are located at the solutions to the following equations written in terms of BA operators Q i : where the function P (u) is expressed through the usual θ 0 (u; τ ) ≡ e 2πiu ; e 2πiτ ∞ e 2πi(τ −u) ; e 2πiτ ∞ as follows: After finding all the solution to the equations (2.6), we can evaluate the residues of the integrand (2.1) and represent the supreconformal index in the following form: where we introduce the following notations. Z (û; ξ, τ ) is the integrand of the index (2.1) and can be generally written in the following form where we introduce a newΓ-function, which is just the usual Γ-function defined as the function of chemical potentials,Γ (u; τ, σ) ≡ Γ e 2πiu ; e 2πiτ , e 2πiσ . (2.10) The κ G in (2.8) denotes the constant prefactor of the index integral representation in (2.1) and is defined as The final ingredient to be explained in (2.8) is the Jacobian H (u; ξ, τ ) for the change of variables u i → Q i (u): Notice that as promised we have concentrated on the particular case of p = q for the index. However, all the equations above may be written for the more general case of τ /σ ∈ Q + . All the corresponding expressions can be found in [43].

Solutions to the BAE
Let us start by very briefly reviewing the solution to the BAE for N = 4 SU(N ) SYM that was used in [42] in order to match the black hole entropy with the large N behavior of the index. The BAE (2.6) in this case take the following form: where u ij ≡ u i − u j . Notice that due to the double periodicity of Q i , if u i is the solution to the BAE then so are u i + 1 and u i + τ . Therefore, there is in fact an infinite number of solutions to the BAE, but we can group them into equivalence classes. In other words, our holonomies u i live on a torus with modular parameter τ . Following [43] we also introduce the convenient parametrization ∆ for the flavor fugacities defined as follows: where r a are the R-charges of the corresponding multiplets. The for the case of N = 4 SYM, this parametrization just takes the form ∆ a = ξ a + τ +σ 3 . Finally, λ in (2.13) is the Langrange multiplier required here since we are considering SU(N ) gauge group instead of U(N ). For the same reason, the holonomies u i satisfy the following constraint: The exact same equations appeared in [9] in the context of the twisted index of N = 4 SU(N ) SYM on T 2 ×S 2 with the topological twist on S 2 . The authors also found the following solutions to these equation at the high-temperature limit τ → 0: 1 In [42] it was shown that the solution above is actually not limited to the high-temperature limit and, in fact, solves the BAE (2.13) exactly for both finite N and τ .
It was also shown in [9] that a similar solution, namely is valid at the high temperature limit of the BAE for Klebanov-Witten (KW) (or conifold) theory [48] and in general for any toric quiver. The index (a) in the expression above refers to the node of the quiver. We will demonstrate here that, just as in the case of N = 4 SYM, this high temperature solution is actually exact in N and τ . Finally, just as in the case of N = 4 SYM we added the constant u to satisfy the SU(N ) constraint (2.15) at each of the nodes. We will refer to the solution (2.17) as basic solution.
To understand how this works, notice that in toric theories each node has an equal number of bifundamentals and anti-bifundamentals. The contribution of each such pair to the BAE is given by where ∆ 1,2 are fugacities corresponding to each of the chiral multiplet and u j . Again, all upper indices denote quiver nodes. In particular, the node (a) is the one that we write the equation for, while (b) and (c) are the ones that connect to (a) through bifundamental chirals. The same expression describes the contribution of one adjoint multiplet provided b = c = a. In this case, ∆ aa is just the chemical potential associated with this multiplet.
Substituting the basic solution (2.17) in the expression above we can use the following chain of manipulations to drastically simplify it N j=1 θ 0 u (ba) Here, between the first line and the second we introduce new summation indices k = i − j in the numerator and k = j − i in the denominator. Between the second line and the third we performed a shift of the index k → k + N in the product i−N k=−1 in the numerator and in the product 1−i k=−1 in the denominator. Now all the terms which do not depend on l, including the ratio of θ-functions product, can be omitted since we can always absorb them into the Lagrange multiplier Λ (a) = e 2πiλ (a) .
Combining the l-dependent part with the one coming from the B(u; τ ) parts we get: , m(n − r) gcd(n, r) , where the last entry of the triplet {m, n, r} is always understood modulo n. We will always be considering only the contribution of the basic solution (2.17) and of the more general family (2.21). In principle, it is not obvious that there are no other solutions to the BAE contributing to the index at the leading N order. This is already true even for the simplest case of N = 4 SYM, let alone for complicated quiver theories, where the situation becomes more involved. Nevertheless, from [42] we know that the class of solutions (2.21) gives the leading contribution to the index which matches the supergravity result. Motivated by this result, we will also proceed to check the contribution of (2.21) against AdS black hole entropy and previous high-temperature index computations [41].

Large-N superconformal index
In this section we will describe the contributions to the index (2.1) at large-N of the basic solution (2.17) and of the more general solutions (2.21) to the BAE. For this, we should just employ Eq. (2.8) which is straightforward though quite a tedious task. However, as we have shown in the previous section, the solutions we exploit are basically the same solutions that have been previously used in [42] in order to find the large-N behavior of the N = 4 SYM index. Hence, we can exploit the derivations used in the latter case. Below, we briefly summarize the crucial expressions and derivations that are required for the calculation and explain, where relevant, the generalization to quiver theories. For all the details of the derivations we refer the reader to the original paper [42].
Contribution of the integrand. The first and most important thing to estimate is the contribution of eachΓ-function (2.10) in the integrand Z (û; ξ, τ ) given in (2.10) in the most general form. In [42] it has been shown that upon substitution of the basic solution (2.17), the leading term of the correspondingΓ-function in the large-N limit is given by: where we have introduced the periodic and discontinuous function There is a number of useful properties of this function: i.e. the function [∆] τ shifts ∆ along the real line and brings it to the fundamental domain For convenience we show this strip in Fig. 1. We also introduce the notation γ for the right boundary of this domain, i.e. the line going between the points 0 and τ . Notice that the function [∆] τ is discontinuous along the lines γ + Z. The latter ones are Stokes lines that divide the complex ∆-plane into regions corresponding to different dominating large-N contributions. Along these lines the different contributions compete with each other and a more accurate calculation including an evaluation of the subleading terms is required.
As a result, expression (2.23) is not defined in the case ∆ = 0 which sits exactly on the Stokes line. At the same time, this value of ∆ is relevant for the contribution of the vector multiplets in (2.10). In this case, as has been shown in [42], the leading-order contribution of theΓ function is given by These two expressions constitute all the building blocks we need for the evaluation of the integrand contribution to the superconformal index.
Contribution of the Jacobian. The second possible contribution to the index comes from the Jacobian H of the transformation between the u i and Q i variables given in (2.12). However, since we consider a quiver theory with SU(N ) groups it is convenient to write the definition of the Jacobian in a clearer form as In the expression above, we have first of all precisely written out the quiver indices 1, . . . , p, where p is the number of nodes. The resulting expression can be seen as a division of the full N p × N p matrix into p × p blocks of the size N × N . Second, notice that due to the SU(N ) gauge groups in the quiver, each set of holonomies u first. As we said, since our quiver is balanced each node (a) has an equal number of outgoing (bifundamental chiral) arrows (ab) and incoming (antibifundamental chiral) arrows (ca). Then, we can write the corresponding operator Q (a) i using the general expression (2.6) in the following form: where we denote by (a, b) the set of all arrows going from the node (a) to the node (b) of the quiver. Therefore, when we take the product over all b ∈ (a, b) for example, the product is over all the nodes that have incoming arrows originating from the node (a). As discussed before, the number of terms of the products in the numerator and the denominator of (2.29) is the same since each node of the quiver has an equal number of incoming and outgoing arrows. Moreover, as before: u Now, we can easily find the logarithmic derivative of the BAE operator (2.29): where for convenience we introduced the following function: Now we can estimate the contributions of each of the logarithmic derivatives (2.30) to the Jacobian (2.28) upon the substitution of the basic solutions (2.17) to the BAE. In particular, we should use Notice that the last two terms appear because we consider only N − 1 holonomies u (a) i as independent variables, while due to the SU(N ) constraint at each node u (a) i . Then, we can rewrite the matrix elements of the Jacobian matrix H as follows: where in the first equality we have used the definition of the BAEs Q (a) i = 1. In addition, we have introduced δ b∈(a,b) , which is one if b is in the set of arrows going from node (a) to node (b) and zero otherwise.
Finally, we also need the elements of the Jacobian (2.28) containing the derivatives with respect to the Lagrange multipliers λ (a) : . Indeed, upon substituting the basic solution (2.17) in the function F (u), we see that its argument is τ 1 At the same time, potential divergences can arise only at zeros of the θ 0 (u) function, which take place at u = m + nτ , m, n ∈ Z. As we see from the above, the latter condition can only be satisfied if the chemical potential ∆ ab ∈ m + nτ . In this case, on the boundaries of the Bethe roots the distribution F u would contribute at the leading order in N [42]. In particular, one can notice that both the integrand Z and the Jacobian H are invariant under τ → τ + r. Hence, once again only the integrand Z contributes to the index at the leading order and the value can be computed from (2.23) and (2.27) by simply substituting τ → τ + r. At the same time, in [42] that S-transformed basic solutions contribute to the index only at order N and are thus subleading. Based on this, the authors argued that among all the solutions (2.21) labeled by {m, n, r}, only those obtained from the basic one (2.16) contribute at leading order. Since our basic solutions and all the terms in the index (2.8) have exactly the same form as in the case of N = 4 SYM, the arguments above are also valid in our case and we will assume that only the T -transformed solutions of (2.17) contribute at leading order to the index.
Summary. Gathering the above derivations, we can write down the leading N 2 term in the superconformal index (2.1). Since we are now convinced that for all toric theories only the basic solution (2.17) and its T -transform contribute at the leading order, we may now directly substitute these solutions into the index (2.8). Hence each time we have an adjoint or a chiral bifundamental multiplet with corresponding chemical potential ∆, we should add a contribution (2.23) to log I. For each vector multiplet, it i.e. for each node of the quiver, we should add a contribution of (2.27) to the log I. Since, as we saw, the Jacobian H always contributes at the subleading order, these are all the relevant expressions we need to calculate the total contribution of the basic solution to the index, namely where p is the total number of the gauge group nodes and in the last term of the Θ({∆}; τ ) function, we summed over all oriented arrows in the quiver graph. This is done to avoid counting twice the contribution of each bifundamental chiral multiplet. In addition, if a = b in the (a, b) pair, then the adjoint multiplet is pictured as one arrow on the graph and should be counted only once.
As discussed above, we should finally include the contributions of the T -transformed solutions (2.36). This is done simply by shifting τ → τ + r in Θ({∆}; τ ). This corresponds to the contribution of one of the T -transformed solutions. In the final step we should find a T -transformation (or equivalently such r ∈ Z) that maximaizes the real part of −πiΘ: This would give us the final expression for the dominant large-N contribution to the index. Still, this expression is not highly informative on its own. In the following sections of the paper, we will thus be considering various toric theories case by case, and applying (2.39) we will be deriving a compact form for the final answer such that the relation between the U (1) triangular anomalies and large-N superconformal index would be clear.
The function I ∞ has a complicated structure. The space of complex parameters ({∆}, τ ) is divided into regions separated by Stokes lines, which are codimension-one surfaces. In each of the regions one solution determined by the value of r that maximizes ImΘ ({∆}, τ + r) dominates over the other solutions. On the Stokes lines themselves, however, two or more values of r 1 , r 2 , · · · ∈ Z yield the same contribution, it i.e.
In this case we have two or more equivalent exponents contributing to I ∞ . We should sum these exponents, but in this case their relative sign should be known in order to compute the final answer for I ∞ . That is because the contributions could potentially cancel each other, forcing us to consider the subleading orders. The particular structure of the parameter space and Stokes lines depends on the particular theory considered. We will consider some details pertaining to the Stokes phenomenon for various toric theories in Section 5.

The conifold theory
We are now ready to apply the derivations of the previous section to particular examples of toric theories. In principle, as discussed above we can directly apply the final expression (2.39) to any toric theory. However, in this section we will study in detail the BAEs, solutions to it and the large-N index of the simplest possible toric theory beyond N = 4 SYM, namely the conifold theory. This theory was originally proposed in [50] as the worldvolume theory of a stack of N D3-branes probing the tip of the conical singularity xy − zt = 0. In the near-horizon limit this gives the AdS 5 × T 1,1 dual of the theory. Here T 1,1 is the conifold which can be seen as the coset SU(2) × SU(2)/U(1). The conifold is in fact a toric manifold with a toric diagram identified with the vectors: The theory is N = 1 supersymmetric with SU(N ) × SU(N ) gauge group and two pairs of bifundamental chiral multiplets A 1,2 and B 1,2 transforming in the N, N and N, N representations of the gauge groups, correspondingly. The corresponding quiver is shown on the Fig.3. The superpotential term is given by The global symmetries of the theory include U(1) R R-symmetry factor, two SU(2) flavor symmetries rotating the A and B doublets, correspondingly, and a U(1) B baryonic symmetry.
These symmetries can also be identified with the isometries of T 1,1 . It is further convenient to use instead of SU(2) flavor symmetries Cartans of their combinations, which we denote by U(1) 1 and U(1) 2 . The matter fields' charges under these global symmetries are summarized in the table below.
Also notice that conifold theory is a particular case of the infinite Y pq class of theories, corresponding to p = 1 and q = 0. We will be considering this class of theories in section 4.1 with less detail.
We now have all the ingredients at our disposal to write down the integral representation (2.1) of the superconformal index for the conifold theory.
where we turned on fugacities v F1,2 for the flavor symmetries U(1) 1,2 , and v B for the baryonic symmetry U(1) B . As usual, we introduce the corresponding chemical potentials according to The expression for κ G is given by (2.11) Throughout this paper we focus on the case where p = q and use the chemical potential τ defined in (2.4). From expression (3.3), it is obvious that the following parameterization of the chemical potentials would be convenient to work with: In particular, with this parametrization the BAEs (2.6),(2.29) for the conifold case may then be written in the following compact form Notice that upon substituting the functions B(u) given in (2.7) into the BAE operators (2.29) we will get certain terms depending on ∆'s but independent of the holonomies u (1,2) i in the exponents. However, we ignored these terms in the expressions written above since they can always be adsorbed into a redefinition of the Lagrange multipliers Λ (i) ≡ e 2πiλ (i) .
It is noteworthy that the parametrization (3.6) can be obtained in two natural steps. First, one should perform a mixing of the flavor and baryonic symmetries. This leads to the new basis of U(1) 1,2,3 global symmetries with the corresponding chemical potentialsξ 1,2,3 turned on for each of them. The charges of the conifold theory matter fields in this new basis are summarized in the table below.
Actually, the proper choice of mixed U(1)-symmetries can be seen directly from the geometric data of the toric manifold. To this end, one should apply algorithms defining the matter content charges from the toric diagrams as proposed in [51]. The same parametrization of the global symmetries has been used in the calculation of the superconformal indices of toric theories at the high-temperature limit in [41]. In Section 4 we will be considering many examples of toric theories, and we will always be using a choice of global symmetries that is determined by the corresponding toric data. The relevant algorithm will be briefly reviewed in that section.
Finally, after the flavor symmetry mixing we should perform the shift (2.14) ∆ a =ξ a + r a τ , (3.8) leading exactly to expressions (3.6). It will often prove important to remember in future derivations that the chemical potentials ∆ a are not independent, rather they satisfy the constraint Let us now show how the basic solution (2.17) satisfies the BAEs (3.7). In particular, let us substitute the basic solution into the first BAE and apply a chain of relations similar to the one used in (2.19). We first wish to concentrate on the transformation of θ-functions: where for the notational brevity we introduced Using (3.10) in the first BAE (3.7), and similar relations in the second, it turns out that all l-dependences drop out and we are left simply with The equations above merely define the values of the Lagrange multipliers Λ (1,2) ; these are irrelevant for the large-N contribution of the basic solution, as discussed in Section 2. However, it is an interesting observation that the values of the Lagrange multipliers obtained from (3.12) precisely reproduce the results of [9] at the high-temperature limit τ → 0. 4 After showing how the basic solution (2.17) solves the BAEs (3.7) of the conifold theory, we are now ready to directly evaluate the index (3.3) using our BAE formula (2.8). In principle, calculations of all possible contributions were performed in Section 2. Nevertheless, to illustrate our estimates for the contribution of the Jacobian (2.12) to the index, let us take a look at the particular explicit form the Jacobian assumes for the conifold theory case, that is Here, each of the diagonal blocks J i has the form (3.14) where the elements A (i) j are all of order N and are given in (2.33). Notice that J i coincides with the Jacobian appearing in the calculations for the case of N = 4 SYM. It can be clearly seen from the expression above that indeed, as expected, the logarithm of the Jacobian log H contributes at order N log N at most. Now using the results of Section 2, we can directly obtain expression (2.39) with the function Θ given by where the function [∆] τ is defined in (2.24). However, this form may be drastically simplified if one remembers that our chemical potentials ∆ a satisfy the constraint (3.9), leading to To simplify the function Θ written above, we then need to understand which of the stripes shown in Fig. 1 Case III :  Fig. 1. The three different relations are Upon substituting these expressions back into the chemical potentials constraint (3.16) and into (3.15), we obtain the final expression for the Θ-function: where [∆ 4 ] τ is defined in (3.16) and the function F ([∆] τ ; τ ) is defined as follows: (3.20) and (2.39) completely define the leading term in the superconformal index of the conifold theory. To find the final expression for the index, it is still necessary to perform an r-extremization as prescribed by (2.39). However, it is quite a technically involved task due to the Stokes phenomenon and at the moment we postpone this question to the future. In Section 5 we will consider the extremization procedure for the particular case of equal chemical potentials ∆ a . Notice that in case I our large-N index behavior (3.19) is in line with the Cardy-like limit computations (1.2) performed in [41].
Finally, as we have mentioned in the beginning of this section the conifold theory is a particular example of an infinite Y pq class of toric theories with the specification p = 1, q = 0. As we will see further in Section 4.1, the general result (4.15) for the large-N index of Y pq theories indeed reduces to the T 1,1 result (3.19) upon this identification.

Other toric models
In this section we extend our analysis to more toric theories and find the large-N asymptotic forms of their indices using the general prescription of section 2, thereby obtaining the entropy functions of the dual black holes. We consider the models discussed in [41], and for most of them (where it is not cumbersome) present the expressions for the entropy function in all the different domains, as in the previous section.
Before going into particular examples of toric quiver theories, let us consider some of their general properties which will be of particular use for us in this section. Toric quiver gauge theories arise as the world-volume theories of D3 branes probing the tip of a toric cone over various five-dimensional Sasaki-Einstein manifolds. In the near-horizon limit, these constructions provide us with the holographic dual description of the corresponding gauge theories. A nice and useful property of the toric models is that all their data can be connected to the geometry of the cones [52,53]. In particular, all this data can be encoded in the toric diagrams of the corresponding cones, which are convex polytopes P built from vectors of the form V i = (·, ·, 1). We summarize below the algorithm of [51] used to extract the assignment of U (1) charges for the chiral fields, as well as other gauge theory data. To illustrate this algorithm, we show all the steps in the example of the conifold theory whose toric diagram is shown in Fig. 2 and the corresponding vectors V I are given in (3.1).
As a first step, and without any extra constructions, we can extract the total number of gauge group nodes F from the area of the polytope P : Another important fact is that the total number of global symmetries (including R-symmetry) d is equal to the number of external points in the diagram. In the case of the T 1,1 diagram, we can see that the total number of gauge group nodes is 2 × 1 and the total number of global symmetries is 4, as we expect from the gauge theory content.
Next, we continue with defining two-component vectors that connect nodes of the toric diagram: v i ≡ V i+1 − V i , and assign numbers for each of them as follows: a i = {a R i , a I i } where I = 1, . . . , (d − 1). In the case of the conifold theory, these vectors are given by To each pair of vectors (ij) = (v i , v j ) such that v j can be rotated into v i clockwise with less than a 180 • turn, we associate a chiral field with a multiplicity of det v T i v T j and charges a i+1 + a i+2 + · · · + a j . In the conifold example, this gives four pairs in total: pair (12) (23) (34) (41) multiplicity In the last line of the table we identify a specific bifundamental multiplet of the conifold theory for each of the pairs. The algorithm for this identification requires an introduction for dimer constructions [52] and we omit it here for simplicity 5 Now, in order to fix the charges of the fields we are free to make any choice of the numbers a i that satisfies the constraints: Notice that for an index computation the precise identification of pairs with particular fields of the quiver theory is not even required. (4.5) Then, a R i will set the assignment of R-charges according to the table of pairs and a I i will assign the charges under the (d − 1) global symmetries U(1) I . The choice presented above precisely reproduces the charges used in Section 3.
The final ingredients we will need to extract from the toric data are the triangular anomalies of the U (1) I global symmetries C IJK . They are simply given by the area of the triangles built from the vertices (V I , V J , V K ) of the toric diagram [54]. Now let us discuss the form of the large-N index we expect to obtain. As we have mentioned in the introduction, the Cardy-like limit computations of [41] has led to the general result (1.2) for the indices of the toric theories. This result, in turn, supported the conjectured form of the entropy function of multi-charged AdS 5 black holes proposed in [31]. However, our computations are slightly different since we do not limit ourselves to any chamber of the chemical potentials. Still, we expect to obtain the same form for this expression, with ∆ a substituted by [∆ a ] τ , at least in some regions of the ∆-space, similarly to the N = 4 SYM case in [42]. In terms of the function Θ, defining the index according to (2.38) we expect the following contribution of the basic solution in one of the chambers: It can be easily checked that in the case of the conifold theory, (3.19) precisely reproduces this expression in case I. As we will see in this section, the same situation takes place for every toric theory we consider.

The family Y pq
Let us start our discussion with the infinite family of models Y pq , given by quiver gauge theories with 2p gauge groups and bifundamental chiral fields [55]. In Fig. 5 we show an example of the Y 32 quiver gauge theory. The notation Y pq denotes a toric manifold with the corresponding diagram shown in Fig. 4, and is parametrized by the vectors: Using the algorithm summarized above, the charges of the fields under the three global U(1) symmetries and the U(1) R R-symmetry can be identified from the form of the toric diagram shown on Fig. 4 as follows, As before, we turn on the chemical potentials ξ i for each of the three U(1) i global symmetries and τ for the R symmetry, and perform the shift (2.14): where we have also defined the usual auxiliary chemical potential ∆ 4 . Now we have everything we need in order to evaluate the superconformal index using the expression (2.39), with the function Θ ({∆ i }; τ ) (2.38) given for the case of the Y pq matter content by the following lengthy expression: Just as in the case of the conifold theory, this expression can be simplified and written in a compact form. In the expression above we have the bracket functions of various sums of chemical potentials, like [∆ 2 + ∆ 3 ] τ for example. In order to simplify them we need to understand how these bracket functions are expressed in terms of the bracket functions of the summands, like [∆ 2 ] τ + [∆ 3 ] τ in the example of the previous sentence. To do it, we first distinguish between different cases similarly to the previous section, and summarize the results in the table below. In this table we fill in the number of the strip in the complex plane to which the particular sum of [∆ i ] τ belongs, with the n-th strip being defined as follows: More explicitly, the green stripe shown in Fig. 1 corresponds to n = 1, while the strips to its left (right) correspond to n = 2, 3, . . . (n = 0, −1, −2, . . . ). As a simple example that illustrates the notations used in this type of table, which we will also employ below when analyzing other theories, let us consider the first line of the table corresponding to case I. As we see, all the relevant sums of [∆ i ] τ are in the first strip, hence we get that case I is given by: Next, after we distinguished between the various cases, we find the following relations for the bracket functions: -24 -Substituting these relations into (4.9), we get the final form of the function Θ (∆ i ; τ ) in the different cases: Case III (4.15) where, defining for convenience the functions (a = {1, 2, 3, 4}): Θ K (∆ i ; τ ) are given by: Notice that Θ I (∆ i ; τ ), which can be written more explicitly as has exactly the same structure as the expected result (4.6). In addition to that, it can be shown that the corresponding coefficients are exactly the triangular anomalies which can be obtained from the toric diagram in Fig. 4. Note also that exactly the same result has been obtained in [41] using the Cardy limit of the superconformal index, and that by specifying to p = 1, q = 0 the expressions in (4.15) reduce to (3.15) obtained in the conifold case.

The family L pqr (p = r)
We next turn to the family of models L pqr , built from quiver gauge theories with p + q gauge groups and bifundamental chiral fields [56][57][58]. The notation L pqr denotes a toric manifold with the corresponding diagram parametrized by the vectors: where we have P q = r − ks, s = p + q − r and we assume that p = r. Using our previous discussion, the charges of the fields under the global symmetries can be read from the toric diagram, and are given in the table below.
As in (4.8), we turn on the chemical potentials ξ i and τ for the U(1) i and U(1) R symmetries and perform the shift (2.14): where we have defined ∆ 4 as before. Then, the index is given at large N by (2.39), where To simplify (4.22) as we did above, we notice that the different cases here are exactly the same as for the Y pq family. Therefore, we can use the same classification along with the relations (4.12), (4.13) and (4.14), and get We notice that here as well, Θ I (∆ i ; τ ) (which corresponds to Case I) has the expected form (4.6) with the coefficients given by the triangular anomalies, and matches the result obtained in [41] using the Cardy limit of the index.

The family X pq
This family is built from quiver gauge theories with 2p + 1 gauge groups and bifundamental chiral fields [59]. The toric diagram of the manifold X pq is parametrized by the vectors: and as discussed above, the corresponding charges of the fields under the global symmetries are as follows, A 234 q − 1 0 1 1 1 6 5 We turn on the chemical potentials ξ i and τ for the U(1) i and U(1) R symmetries and perform the shift (2.14): . . . , 4) , where we have defined the auxiliary ∆ 5 in the usual way. Then, the index is given at large N by (2.39), where As before, to simplify this expression we need to consider the different domains of the bracket functions. However, for compactness we only present here the result for Case Ia, in which Case Ia : (note that all the sub-sums are also in the same strip). For this case, we have which is of the same form as (4.6) with the coefficients given by the triangular anomalies, and matches the result obtained in [41] using the Cardy limit of the index.

The model SPP
We consider the SPP (short for "suspended pinch point") model, the theory living on a stack of N D3-branes probing the SPP conical singularity x 2 y = wz. This gauge theory is described by the quiver given in Fig. 6 below. The toric manifold in this case is parametrized by the vectors: and the corresponding charges of the fields under the global symmetries are given in the following table: We turn on the chemical potentials ξ i and τ for the U(1) i and U(1) R symmetries and perform the shift (2.14): . . . , 4) , where we have defined the auxiliary ∆ 5 in the usual way. Then, the index is given at large As usual, to simplify (4.28) we need to consider the different domains of the bracket functions.
Substituting these relations into (4.28), we get the expression for Θ (∆ i ; τ ) in the different -30 -cases: Here, Θ Ia (∆ i ; τ ) (which corresponds to Case Ia) is of the same form as (4.6) with the coefficients given by the triangular anomalies, and matches the result obtained in [41] using the Cardy limit of the index.

4.5
The models F 0 , dP 1 and dP 2 The large-N indices of these models can be easily obtained from those corresponding to some of the theories we considered above, and will therefore not be analyzed separately. In particular, the large-N index of the model F 0 is just twice that of the conifold theory, while the large-N indices of the models dP 1 and dP 2 are the same as those of the theories Y 21 and X 21 , respectively.

The model dP 3
We turn to the dP 3 model, the theory living on a stack of N D3-branes probing the tip of the Calabi-Yau cone whose base is dP 3 (i.e. the third del Pezzo surface). This gauge theory is described by the quiver given in Fig. 7 below.  Figure 7: Quiver diagram corresponding to the dP 3 model. Each node represents an SU (N ) gauge group and each X ij field transforms in the bifundamental representation, i.e. in the N of the i-th node and in the N of the j-th node.
The toric manifold in this case is parametrized by the vectors: (1, 0, 1) , and the corresponding charges of the fields under the global symmetries are given in the following table: We turn on the chemical potentials ξ i and τ for the U(1) i and U(1) R symmetries and perform the shift (2.14): where we have defined the auxiliary ∆ 6 in the usual way. Then, the index is given at large N by (2.39), where To simplify this expression, we need to consider the different cases for the bracket functions. However, for compactness we only present here the result for Case Ia, in which Case Ia : (note that all the sub-sums are also in the same strip). For this case, we have which is of the same form as (4.6) with the coefficients given by the triangular anomalies, and matches the result obtained in [41] using the Cardy limit of the index.

The model (P)dP 4
We consider the (P)dP 4 model, defined using the fourth del Pezzo surface. The quiver corresponding to this theory is given in Fig. 8 below.
We turn on the chemical potentials ξ i and τ for the U(1) i and U(1) R symmetries and perform the shift (2.14): where we have defined the auxiliary ∆ 7 in the usual way. Then, the index is given at large N by (2.39), where To simplify this expression, we need to consider the different cases for the bracket functions. However, for compactness we only present here the result for Case Ia, in which Case Ia : (note that all the sub-sums are also in the same strip). For this case, we have which is of the same form as (4.6) with the coefficients given by the triangular anomalies, and matches the result obtained in [41] using the Cardy limit of the index.

Extremization and Black Hole Entropy
As discussed in the introduction, black hole entropy can be obtained from the Legendre transform of the entropy function [31,34]. The latter may be formulated in terms of the dual CFT parameters, and has been shown to coincide with the Cardy-like limit of the index for the case of N = 4 SYM dual to the AdS 5 × S 5 black hole with three electric charges Q 1,2,3 and two angular momenta J 1,2 [21,22]. Its entropy function in particular was shown to be equal to where ω 1,2 are chemical potentials conjugate to two angular momenta and X 1,2,3 are the ones conjugate to three charges. In [42] the large-N behavior of the index was derived using a BAE analysis. In particular, it has been shown that the basic contribution to the index is dominant for the case of the dual single-centered black hole 6 . Moreover, the authors have shown that (5.1) is reproduced by their large-N index in this case upon the identification where ∆ a are the same chemical potentials (2.14) we use in the present work and d = 3 for the SYM case. Their observations hold under an assumption that the ∆s belong to the region analogous to the one we call "case I" through Sections 3 and 4. In this case, [∆ a ] τ satisfy the constraint which is precisely the constraint that should be satisfied by X a . Another argument in favor of the identification (5.2), as shown in [42] for the case of SYM, is that X a should satisfy the same strip inequality as [∆ a ] τ do, namely Similar arguments work also for black hole with d > 3 charges, which on the dual side correspond to theories with d global U(1) symmetries. We thus propose, based on the results of Sections 3 and 4, the following entropy functions for black holes dual to toric theories: Here, the coefficients C IJK are defined by the triangular anomalies of the corresponding U(1) global symmetries, and the chemical potentials X a satisfy the constraint (5.3), Notice that this is exactly the proposal of [32] confirmed at the Cardy-like limit in [41].
Let us now briefly revise the procedure of extracting the entropy from the entropy functions S, giving a concrete illustration of the simplest case of N = 4 SYM. The black hole entropy is computed as the Legendre transform of the entropy function S, given by the following function: computed at the critical point Here, as usual Q I are the black hole d electric charges and J i are its two angular momenta. We additionally introduce a Lagrange multiplier Λ to impose the constraint (5.6). At this point one should, in principle, perform the technically challenging task of solving (5.8) and substituting the solution into (5.7). However, Ref. [31] proposed a nice trick to circumvent this complication and obtain the black hole entropy rather easily. Specifically, it was noticed that the entropy function (5.1) is a monomial of X so it can be rewritten in the following form In more general cases of toric theories, the entropy function (5.5) is a homogenous polynomial of degree three so the same argument applies. This relation together with the equations of motion (5.8) allows us to express the black hole entropy (5.7) in terms the Lagrange multiplier Λ alone, as follows: Now, the only thing left to do is to express Λ in terms of the charges Q I and angular momenta J i . For example, in the N = 4 SYM case the equations of motion (5.8) take the following form: Then it can be seen that the charges and angular momenta enjoy the following cubic relation: where we introduced the coefficients p i defined as From the requirement that the entropy be real and positive, one can easily deduce that the charges should satisfy the constraint p 0 = p 1 p 2 and that the entropy of the dual black hole is given by which precisely reproduces the supergravity result [21]. A similar analysis has also been performed for the black hole duals of Y pp and L aba theories in [41]. However, the algorithm described above is technically complicated for other examples of toric theories. As we will be seeing momentarily for the conifold theory case, one of the problems is that the cubic relation (5.12) does not always exist. At the same time, solving the equations of motion (5.8) directly is technically too complicated. It would be interesting to understand if these problems can somehow be worked around in the future. Taking another route, however, we will try to analyze the structure of the large-N index of the conifold theory in the space of chemical potentials ∆ I , for the specific case dual to a black hole with equal charges. A similar analysis was performed in [42]. In particular, the authors identified Stokes lines and regions, and found that for certain values of charges and angular momenta single-centered black hole solutions become unstable, as the index corresponding to these values ends up sitting on a Stokes line.

The case of equal charges
Let us now consider the large-N index of the conifold theory given in (3.19) and (3.20) and concentrate on the following simple choice of chemical potentials: where the last inequality comes from the constraint (5.3). Notice that we are now left with the single chemical potential τ . Also notice that the index I 2τ −1 4 ; τ is periodic under the shift τ → τ + 2, from which it follows that we may restrict to 0 ≤ Reτ < 2. As we will later see, the choice (5.15) has a nice interpretation in terms of a dual black hole with equal charges Q I .
The contribution of the basic solution and its T -transforms to the logarithm of the superconformal index reduces in this case to The definition of the different cases may be found in (3.18). To work with this expression, we must understand what are the values of [∆] τ +r for various r ∈ Z and which cases in (5.16) they correspond to. In particular, it can be checked that we should distinguish between cases of even and odd r as follows: It can also be checked that for even r, one should use "case I" in (5.16), while for odd r "case III" is the suitable one. Now we finally would like to understand the structure of the large-N index in the complex plane of the single parameter we are left with, namely τ . For this we need to understand whichr are maximizing ImΘ (∆; τ + r) in the different regions. If in some region there is no r maximizing this function, one then needs to sum over infinitely many exponents contributing at the same order. This would require extra information and probably the inclusion of subleading terms to evaluate index. That situation precisely corresponds to being on the Stokes line. In our case, the imaginary part of the Θ function (5.16) is given by ImΘ (∆; τ + r) = This is remarkably similar to the SYM result [42], with only two differences between the two cases. One difference is that in N = 4 SYM there are three possibilities with r = 0, 1, 2 mod 3. For r = 1 mod 3, [∆] τ is not defined since this value is precisely on the Stokes line. The second is that the expressions for r = 0 and r = 2 mod 3 exactly reproduce our expressions for even and odd r, correspondingly, up to an overall factor. However, this is irrelevant since this coefficient only defines the value of the index in the particular region but not the region itself.
As illustrated in Fig. 9, there exists a maximum of ImΘ (∆; τ + r) given in (5.18) with an even r if τ resides in the interior of the semicircle defined by Reτ +r even > 3 |τ +r even | 2 . This pair of semicircles appearing in the fundamental range 0 ≤ Reτ < 2 periodically features in all the images of the fundamental range under τ → τ + 2. 7 The region external to the semicircles admits no dominant contribution; there, the supremal value for ImΘ as a function of r is approached asymptotically in the limit where r → ±∞ and is equal to 1 2 Imτ , as can be readily seen from (5.18). Now, when we know the structure of the large-N index for the conifold theory in the complex τ -plane, we can also try to explore the position of the critical points of the entropy function in the same plane. In particular, according to our proposal (5.5) the entropy function of the black hole dual to the conifold theory should be given by with the chemical potentials X a satisfying the constraint (5.6). The equations of motion (5.8) 7 The same picture emerges for SYM, with the difference that there the periodicity is τ → τ + 3 and, correspondingly, the fundamental range is 0 ≤ Reτ < 3. in this case take the following form: In general, it is very hard to solve these equations and a relation of the form (5.12) is also impossible to find. However, if we require equal chemical potentials ∆ (5.15) in the index, then the identification (5.2) immediately implies: Substituting the above expression into the equations of motion (5.22), and assuming ω 1 = ω 2 = τ we obtain that all electric charges are equal as well as two angular momenta are equal to each other: We are now able to use the usual trick and construct a simple cubic relation between the charges and the angular momenta: (Q + Λ) 3 + c (J − Λ) 2 = Λ 3 + p 2 Λ 2 + p 1 Λ + p 0 , where p 0 = Q 3 + cJ 2 , p 1 = 3Q 2 − cJ , p 2 = 3Q + c , c ≡ 27 16 ν . (5.27) It is convenient to measure the charges in units of cν, so we introduce Q ≡ cνQ , J ≡ cνJ , Λ ≡ cνΛ. In this case all our expressions, except for the equations of motion (5.25) appear to be the same as in the case of N = 4 SYM in [42]. In particular, we obtain the following constraint for the charges and entropy 0 = p 1 p 2 − p 0 = c 3 ν 3 8Q 3 + 3Q 2 − 2 3Q + 1 J −J 2 , S = 2π √ p 1 = 2πcν 3Q 2 −J . It is interesting to note that the solution for τ is exactly the same as for the case of N = 4 SYM [42]. Moreover, even the parametrization we use in (5.29) is the same as in the SYM case, so the expression appears to be exactly the same. The form of the solution for the chemical potential X is also identical to the one found in the case of N = 4 SYM, up to an overall coefficient only.
It is interesting to see how these results look like in Fig. 9, where we show the structure of the large-N index dual to the entropy of a black hole with equal charges. Using the parametrization (5.29) for the solution (5.30), the line of critical points is shown in orange in Fig. 9. It is not surprising that the picture we obtain is identical to the one obtained for N = 4 SYM in [42]. In particular, we see that for small values of µ below a certain critical value µ < µ * , the solution is in the region with no dominant contribution to the index. Hence, in this region we should take into account other solutions beyond the basic one to evaluate the index, and probably the leading order contribution will no longer be of order N 2 . In [42] it was proposed that such behavior signals the instability of small (i.e. with charge below some critical value, Q < Q * ) single-centered black hole towards other supergravity solutions. On the other hand, when we increase µ at some point it crosses the semicircle (5.19) and the solution we found becomes dominant. Therefore, large single-centered black holes (with charge Q above the threshold Q * ) appear to be stable solutions of supergravity. The critical value of µ * can be found from the equation Reτ = 3|τ | 2 defining the boundary of the semicircle (5.19). Using the parametrization (5.29) together with the solution (5.30) we can define the following values of critical charges and entropy: µ * = 2 3 , τ * = 1 + i 6 , Q * = 3 2 ν, J * = 13 8 ν, S * = 9π 4 ν . (5.31) An identical analysis can be performed for the interior of the semicircle (5.20), where the large-N index and the corresponding entropy function are dominated by the r = 1 contribution and take the following form log I ∞ (∆; τ ) = −πiN 2 Θ (∆; τ + 1) == −πiN 2 (2τ + 1) 3 16(τ + 1) 2 − 2 . (5.32) In Fig.9 we show the line of critical points in green. As we can see, the form of the curve just mirrors the line of critical points for the case of r = 0.
As we have mentioned repeatedly in this section, the picture we obtain for the case of equal black hole charges Q I is exactly identical to the one obtained in [42]. It would be interesting to understand if the same pattern takes place for other toric theories. However, note that conifold theory is special because in general the requirement (5.23) of equal chemical potentials does not translate into a condition of equal electric charges Q I of the dual black hole. This is due to more a complicated structure of the equations of motion (5.8).

Discussion and Outlook
In this paper we used a BAE analysis to obtain the large-N behavior of the superconformal index for a variety of toric gauge theories. We demonstrated that the basic solution (2.16) to the BAEs of N = 4 SYM can be directly generalized to (2.21), which solves the BAEs for any toric quiver theory. Using these solutions, we derived the relatively simple expressions (2.38) and (2.39) for their contribution to the large-N limit of the superconformal index for arbitrary toric quivers. Studying these expressions for a number of theories, including the infinite families Y pq , X pq and L pqr , we found a way to rewrite them in a form that coincides with (1.2) for certain ranges of parameters. The latter is the result of the index computation in the Cardy-like limit [41]. Remarkably, this also reproduces a conjecture made in [32] for the entropy function of a multi-charged rotating AdS 5 black hole.
Finally, we also studied the large-N index structure for the conifold theory in the specific case corresponding to a dual black hole with four equal charges and two equal angular momenta. In this case, we identified Stokes lines and regions in the space of the single complex parameter τ , which is just the chemical potential dual to the angular momentum of the black hole. We further performed an extremization for this case. Curiously, our picture exactly reproduces the one found in the N = 4 case presented in [42].
Many related open questions are still waiting to be answered. First, our results are based only on the basic solutions (2.17) and their SL (2, Z) transforms. However, one should study closely the BAEs for particular theories and try to understand if there are other solutions, at which order of N they contribute to the index, and what are the corresponding observables in the bulk of AdS.
Another open question is the extremization procedure of the large-N index. In this paper, we presented results for the very particular case of T 1,1 with all the chemical potentials being equal. The authors of [41] managed to perform an extremization in the cases of the Y pp and L aba theories. However, for all the other toric cases we have failed to perform a Legendre transform since the algorithms of [31] do not work for them and a straightforward solution is too complicated. It would be interesting to resolve these problems in one way or another.
Finally, it would be interesting to understand if the BAE approach can be used in order to study superconformal indices of larger classes of theories. Of particular interest are theories with matter in the fundamental or anti-fundamental representations. It looks like the basic solution does not work in this case and that one needs to perform a more complicated analysis of the corresponding BAE equations.