Finite N corrections to the superconformal index of orbifold quiver gauge theories

We investigate the AdS/CFT correspondence for quiver gauge theories realized on D3-branes put on abelian orbifolds by using the superconformal index. We assume that on the gravity side the finite N corrections of the index are reproduced by D3-branes wrapped on three particular three-cycles in the internal space Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{Y} $$\end{document}, the abelian orbifold of S5. We first establish the relation between baryonic charges on the gauge theory side and the D3-brane wrapping numbers and holonomies on D3-branes. Then we confirm our proposal by comparing the results of localization for gauge theories and the results on the AdS side including the contributions of D3-branes and excitation on them for many examples. We only focus on the leading finite N corrections starting from qN , and leave the sub-leading corrections starting at qkN (k ≥ 2) as a task for the future. We find complete agreement for the leading corrections in all examples.


JHEP10(2019)243
1 Introduction An N = 1 gauge theory is realized on the worldvolume of N D3-branes put at the apex of a Calabi-Yau cone X . If N is sufficiently large the D3-brane system is well described as the classical supergravity solution AdS 5 × Y, where Y is the base of the cone X . The AdS/CFT correspondence [1][2][3] claims that the type IIB string theory in this background is dual to the superconformal gauge theory realized on the D3-branes at IR. This duality has been tested by calculating various quantities on the both sides and confirming their agreement. In this paper we focus on the superconformal index [4]. We define the index as a formal power series of q 1 2 with the coefficient of each term being a Laurent polynomial of other fugacities. A BPS operator with the dimension d and the right-handed spin j contributes ∝ q d+j to the index. See (1.4) for an explicit definition.
In the large N limit the agreement of the superconformal index has been confirmed in a large class of the internal spaces Y and the corresponding superconformal field theories. For N = 4 U(N ) SYM corresponding to Y = S 5 , the large N index was calculated on the both sides in [4] and agreement was confirmed. The AdS/CFT correspondence for orbifolds S 5 /Γ was suggested in [5,6]. The agreement of the index for the orbifolds S 5 /Z n with an A n−1 type fixed locus S 1 ⊂ Y and the corresponding quiver gauge theories was confirmed in [7]. It was found that on the gravity side not only the gravity multiplet but also the tensor multiplets living on the fixed locus contribute to the index. The index of the Kaluza-Klein modes in AdS 5 × T 1,1 was calculated in [8]. The large N index of the corresponding quiver gauge theory, the Klebanov-Witten theory [9], was calculated in [10] and the agreement was confirmed.
If N is finite the correspondence is modified for operators with dimension of order N or larger. For example, in N = 4 SO(2N ) SYM the Pfaffian operators with dimension N correspond not to Kaluza-Klein modes of supergravity fields but to D3-branes wrapped around topologically non-trivial cycles in S 5 /Z 2 [11]. Similar relations hold for baryonic operators in quiver gauge theories whose dual geometries have topologically non-trivial three-cycles [12].
D3-branes play a role in finite N corrections even when the internal space does not have topologically non-trivial three-cycles. In the case of N = 4 U(N ) SYM the one-toone correspondence between BPS operators and Kaluza-Klein modes is broken down for operators with dimension of order N or larger. On the gauge theory side this is because single trace operators with the length L > N are not independent but decomposable into shorter single-trace operators. We can explain this on the gravity side by assuming that the operators correspond not to point-like gravitons but to giant gravitons: D3-branes wrapped on topologically trivial three-cycles in S 5 . 1/2 BPS giant gravitons were constructed in [13]. Their angular momentum J in S 5 has the upper bound J ≤ N , and the absence of giant gravitons with J > N corresponds to the absence of independent single-trace operators with the length L > N . Indeed, the BPS partition function [14] of N = 4 U(N ) SYM with finite N was exactly reproduced in [15] based on the idea of [16] by the geometric quantization of 1/8 BPS giant gravitons constructed in [17]. There is also a complementary way to reproduce the same BPS partition function [18] by using giant gravitons expanded in AdS 5 [19,20].

JHEP10(2019)243
The result in [16] was extended in [21] to S-fold theories including the N = 4 SO(2N ) SYM. The exact BPS partition function was derived by the geometric quantization of BPS configurations of D3-branes in S 5 /Z k , where Z k is the S-fold action transforming the (p, q)-string charges non-trivially. It is natural to attempt a similar derivation for the superconformal index, and indeed in [22] it was shown for the S-fold theories that the leading finite N corrections can be reproduced as the index of fluctuation modes on D3branes wrapped around particular three-cycles. We mean in this paper by "the leading finite N corrections" the corrections starting from O(q N ). In the case of S-fold theories we also have "the sub-leading corrections" starting from O(q 2N ), which was not studied in [22].
The purpose of this paper is to calculate the finite N corrections to the superconformal index for ordinary orbifolds Y = S 5 / Γ, which have much more variety than the S-folds. (We use Γ for the orbifold group rather than Γ and Γ for the dual group Hom( Γ, U(1)) because in the following we use the dual group more frequently than the orbifold group itself.) For an N = 1 supersymmetry to be preserved Γ must be a finite subgroup of SU(3) acting on C 3 coordinates (X, Y, Z). 1 We restrict our attention to the toric case. Namely, we assume that Γ is abelian, and is a subset of the Cartan subgroup H = U(1) 2 ⊂ SU (3).
It is convenient to represent h ∈ H as the 3 × 3 matrix acting on (X, Y, Z): and R Z are the Cartan generators of su(4) ∼ so(6) acting on X, Y , and Z, respectively, and w I (h) (I = X, Y, Z) are complex numbers with the absolute value 1 depending on h ∈ H. For this to be an element of SU(3) w I (h) must satisfy Generically the orbifolding breaks the N = 4 supersymmetry down to N = 1. Let Q a and Qȧ be the left-handed and the right-handed unbroken supercharges, respectively. To define the superconformal index we use Q1 such that ∆ = 2{(Q1) † , Q1} is given by The index is defined by 2 Only BPS operators saturating the BPS bound ∆ ≥ 0 contribute to the index and hence the index is independent of x. H is the dilatation and J and J are the left-and rightangular momenta normalized so that the eigenvalues are quantized with unit 1/2. u I are JHEP10(2019)243 SU(3) fugacities satisfying u X u Y u Z = 1. We also use two independent variables u and v related to u I by (u X , u Y , u Z ) = (u, v u , 1 v ). b and ζ collectively represent baryonic charges and the corresponding fugacities, respectively, which will be discussed later in detail.
From the viewpoint on the gravity side this index is expected to be factorized into two factors: (1.5) I KK is the contribution of Kaluza-Klein modes of massless fields, which has been already studied in the literature. I KK gives the exact index in the large N limit. The purpose of this paper is to investigate the other factor, I D3 , which gives finite N corrections due to D3-branes wrapped on three-cycles in Y. Again we focus only on the leading corrections starting from O(q N ), and do not pay attention to the sub-leading corrections starting from O(q kN ) with k ≥ 2 depending on the sector we consider. The rest of this paper is organized as follows. In the next section we summarize the toric diagrams and the quiver gauge theories for abelian orbifolds Y = S 5 / Γ. In section 3 we review how the superconformal index in the large N limit is calculated on the gravity side as the contribution of Kaluza-Klein modes in the orbifold Y = S 5 / Γ. In section 4 we explain a prescription to obtain the leading finite N corrections from wrapped D3-branes based on the analysis in [22]. We first discuss the relation between the wrapping number of D3-branes and the baryonic charges in the quiver gauge theory, and then we give a prescription to calculate the contribution in each wrapping sector. In section 5 we apply the method to examples with N = 2 and confirm the agreement between the results of the localization and those of the D3-brane analysis. Section 6 is devoted to summary and discussions. The appendix contains results for N = 3 for some of orbifolds in section 5.
To write down indices we use characters associated with the spin and flavor symmetries. They are defined as follows.
The spin characters χ J n are defined by The u(2) characters χ n (a, b) are defined by These are used to give the index of a theory with an SU(2) flavor symmetry. In all examples in this paper the SU(2) acts on Y and Z, and the arguments of the characters are u Y = v u and u Z = 1 v . We use the short-hand notation χ n = χ n ( v u , 1 v ). χ (r 1 ,r 2 ) are the su(3) characters for representations with the Dynkin labels (r 1 , r 2 ). For the fundamental and the anti-fundamental representations these are given by The H plane with the lattice P for the orbifold group Γ = Z 5 generated by diag(ω −2 5 , ω 5 , ω 5 ) is shown. L is expressed as the set of intersections of three sets of parallel lines. The lattice P is shown by dots. (b) A triangle picked up from (a) is shown. This is nothing but the toric diagram of the orbifold C 3 /Z 5 .
For a general representation (r 1 , r 2 ) it is given by (1.9) For later use we define ω n by ω n = exp 2πi n . (1.10) 2 Orbifolds and quiver gauge theories

Toric diagrams
Orbifolds we analyze in this paper are a special class of toric Calabi-Yaus, and their structure can be expressed by using toric diagrams. To define the toric diagram it is convenient to define α I (I = X, Y, Z) by The parameters α I can be regarded as (redundant) coordinates of the covering space H = R 2 of H. These are angular variables with the period 1. Let L be the associated lattice defined by α I ∈ Z. H and H are related by H = H/L. In figure 1 (a) the lattice L is expressed as the set of intersections of three sets of parallel lines α I ∈ Z, which give a tessellation of the plane H by congruent triangles. Γ is a finite subgroup of H, and is expressed as a lattice P in H, which includes L as a sublattice. In other words, P is a refinement of L. Thanks to the periodicity and the Z 2 JHEP10(2019)243 rotational symmetry of the lattice P all triangles contain points in P in the same way, and we can pick up one triangle to represent the orbifold group Γ (see (b) in figure 1). This is nothing but the toric diagram of the orbifold. 3 Note that if we express g ∈ Γ in the matrix form (1.1) each of diagonal components w I ( g) is a one-dimensional representation of Γ. Namely, w I can be regarded as elements of the dual group Γ = Hom( Γ, U(1)).
(2.2) Furthermore, the group Γ is generated by the elements w I . We can specify the orbifold X = C 3 / Γ by giving a set of relations satisfied by w I , which must always include w X w Y w Z = e, where e is the identity element of Γ.

Quiver gauge theories
The quiver gauge theory realized on an orbifold X = C 3 / Γ is obtained by the standard prescription [23,24] as follows. We start from N = 4 SYM with the gauge group U(|Γ|N ), where |Γ| is the order of Γ, which is the same as the order of Γ. The N = 4 vector multiplet consists of an N = 1 vector multiplet V and three N = 1 chiral multiplets Φ I (I = X, Y, Z). Because we consider abelian orbifolds we can discuss the projection on each N = 1 multiplet separately. For the vector multiplet V the action of g ∈ Γ is where u ∈ Hom( Γ, U(|Γ|N )) represents the action associated with a holonomy. In this paper we consider only the case that all gauge groups have the same rank, and it is realized by taking where R( g) is the regular representation of Γ. The regular representation of a finite group is the direct sum of all irreducible representations. Because Γ is abelian irreducible representations are identified with elements of the dual group (2.2) and the regular representation is given by Correspondingly, we divide |Γ|N × |Γ|N matrix V into |Γ| 2 blocks V g 1 g 2 of size N × N labeled by g 1 , g 2 ∈ Γ. The action of g ∈ Γ on each block is

JHEP10(2019)243
Then the projection leaves diagonal blocks V gg . Let U(N ) g = SU(N ) g × U(1) g be the symmetry corresponding to the block V gg . We define G is the gauge group. The gauge fields for G 0 B are decoupled in the IR and G 0 B becomes global symmetry. In general G 0 B is broken by anomalies to a subgroup G B ⊂ G 0 B . For the chiral multiplets the orbifold action is the composition of the C 3 rotation given by (1.1) and the holonomy action like (2.3). Corresponding to the irreducible decomposition (2.5) the chiral multiplets are also decomposed into blocks (Φ I ) g 1 g 2 of size N × N . (Φ I ) g 1 g 2 belongs to the bi-fundamental representation (N, N ) of SU(N ) g 1 × SU(N ) g 2 , and in the following we use the notation Φ g 1 →g 2 I . The g ∈ Γ action on each block is 8) and this block remains after the orbifold projection if and only if the relation holds. 4 This condition determines the matter contents of the quiver gauge theory. Namely, we are left with 3|Γ| bi-fundamental fields Φ g→w I g I labeled by I = X, Y, Z and g ∈ Γ. With this information it is easy to draw the quiver diagram. We first draw |Γ| vertices corresponding to the elements of Γ. Each of them represents an SU(N ) gauge group. A chiral multiplet belonging to the (N, N ) representation of SU(N ) g 1 ×SU(N ) g 2 is represented as an arrow from the vertex g 1 to g 2 . Because Γ is a finite abelian group and can be given as a subgroup of U(1) 2 it is natural and convenient to draw the diagram on the torus, and such a diagram is called a periodic quiver diagram.
Once we obtain the field contents of the quiver gauge theory we can in principle calculate the superconformal index for an arbitrary finite N by the localization formula (2.10) The plethystic exponential Pexp is defined as I g v is the single particle index of the SU(N ) g vector multiplet:

JHEP10(2019)243
where χ g adj is the character of the adjoint representation of SU(N ) g . The single-particle index of the chiral multiplet Φ g 1 →g 2 , (2.13) where χ g N and χ g N are the characters of the fundamental and the anti-fundamental representations of SU(N ) g . The explicit forms of the characters are (2.14) Because the gauge groups are SU(N ) the gauge fugacities z g a (a = 1, . . . , N ) are constrained by N a=1 z g a = 1.
We denote the set of the |Γ| fugacities collectively by ζ and it can be regarded as an element of the classical baryonic symmetry G 0 B . If the baryonic symmetry is anomalous and is broken to its subgroup, ζ takes only values in the anomaly-free subgroup.
dµ is the integration over the gauge fugacities z g a defined by (2.16) Although the single-particle indices (2.13) include fractional powers of ζ g only integral powers remain after the integration over the gauge fugacities.

Baryonic charges
Let Λ 0 B be the charge lattice associated with the classical baryonic symmetry G 0 B . An element of Λ 0 B is specified by the set of |Γ| U(1) g charges B g ∈ Z satisfying g∈Γ B g = 0, (2.17) because all matter fields belong to bi-fundamental representations. 5 Let b g be the formal basis vector of the U(1) g charge. A general element of Λ 0 B is given by

JHEP10(2019)243
Now let us consider the effect of anomalies. If the theory is chiral, the presence of instantons causes violation of the conservation laws of the charges B g . Let N g be the instanton numbers of SU(N ) g . Then the baryonic charges change in this instanton background by where n g→g is the number of the arrows in the quiver diagram from g to g . Arrows in the opposite direction are counted by −1. A single SU(N ) g instanton changes the vector (2.18) by For the orbifold quiver gauge theories this is given by (2.21) The quantum baryonic charge lattice Λ B taking account of this breaking of the conservation laws is given by ζ b in the definition of the index (1.4) is a short-hand notation for g∈Γ ζ Bg g . If we treat fugacities ζ g as independent variables among them we can extract the index I b for each b by the expansion (2.24) However, the physical meaning of I b is not clear because the classical baryonic charges are not conserved due to the anomalies. Instead of defining the index for each b we should treat equivalence classes b ∈ Λ B as conserved charges and define the index for each of such classes by The index for a class b can be directly extracted from the index (1.4) by an expansion similar to (2.24). We impose the conditions

JHEP10(2019)243
on the fugacities ζ. This is equivalent to the requirement that ζ is an element of the anomaly free baryonic symmetry group G B ⊂ G 0 B . Then ζ b depends on b only through the equivalence class to which b belongs and we denote it by ζ b . We can rewrite the expansion (2.24) as (2.27)

Large N limit
In this section we review known results for the superconformal index in the large N limit. As is mentioned in the introduction it is given on the AdS side by I KK , the contribution of Kaluza-Klein modes of massless fields. We can divide this into two factors I grav and I tensor : the index of the gravity multiplet in the bulk and the index of the tensor multiplets localized on the fixed loci.

Gravity multiplet
Let us first consider the gravity multiplet contribution I grav . It is given by where P is the orbifold projection operator and I grav is the single-particle index of the Kaluza-Klein modes in AdS 5 × S 5 . For S 5 without orbifolding the modes belong to the [25,26]. (See [27] for the notation for the superconformal representations.) The corresponding index is given by [ The same index is obtained as the large N limit of the index of N = 4 SYM with the gauge group U(N ). Because we consider quiver gauge theories with the gauge group SU(N ) |Γ| , we define I grav by subtracting the contribution of the N = 1 vector multiplet Namely, we define I grav by We define the projection operator P so that it picks up the Γ invariant terms from the index. Let f be a single-particle index defined by summing up modes in S 5 :

JHEP10(2019)243
We here focus on the dependence on the SU(3) fugacities u I . Then the projected index Pf is given by At the second equality we used g = I=X,Y,Z w I ( g) R I .

Tensor multiplet
The tensor multiplet contribution I tensor is present when the orbifold has fixed loci [7]. The presence of fixed loci and their types are easily read off from the toric diagram. The orbifold Y = S 5 / Γ has fixed points if there exist g ∈ Γ such that g = e and w I ( g) = 1 for one of I = X, Y, Z. In the toric diagram such an element is expressed as a dot on an open edge. (By an edge side we mean an edge with two endpoints excluded.) For example, if there are k − 1 dots on an open edge and the edge is divided by them into k pieces, then the k − 1 elements of Γ associated with the dots together with the identity element form Z k ⊂ Γ acting only on two coordinates, and the fixed locus is the set of A k−1 type singularities.
Let us consider the orbifold X = C × C 2 /Z n with the orbifold group generated by diag(1, ω n , ω −1 n ). In this case an N = 2 supersymmetry is preserved. This has the A n−1type fixed locus S 1 ⊂ Y given by Y = Z = 0. In the six-dimensional space AdS 5 × S 1 the A n−1 type N = (2, 0) theory lives, and n − 1 tensor multiplets in the theory contribute to the index [7]. The Kaluza-Klein modes of a single tensor multiplet belong to the series of N = 2 superconformal representations E m(0,0) (m = 1, 2, 3, . . .), and the corresponding single-particle index is (We can easily calculate this index by using the method in [28].) Similarly to the gravity multiplet contribution, we subtract the contribution of the N = 1 vector multiplet (3.3), and define I tensor In [7] this was extracted from the index of the gauge theory associated with the orbifold.
We can generalize this to other fixed loci appearing along intersections of S 5 and X I planes.

JHEP10(2019)243
We denote the corresponding single-particle index by I tensor I , and it is given by Interestingly, these are the same as terms appearing in I grav in (3.4). If Z n is a proper subgroup of Γ and Γ includes elements that act non-trivially on the fixed locus, we need to perform the corresponding projection. It is realized by the projection P defined in (3.6). In a general situation in which the orbifold has type A n I −1 singular locus on the X I -plane the total contribution of the tensor multiplets is given by (3.10)

Finite N corrections from wrapped D3-branes
The purpose of this section is to give the prescription to calculate the leading finite N corrections to the index from wrapped D3-branes based on some assumptions. The wrapping number of D3-branes in Y is related to the baryonic charges. We first define the wrapping number and then establish the relation to the baryonic charges.

Wrapping number
Let us define S 5 as the subset of C 3 by and represent the worldvolume of a D3-brane in S 5 as the intersection of a surface in C 3 and the S 5 . For the D3-brane to be BPS the function F must be holomorphic [17].
Here we are interested in topological aspects and the holomorphy is not assumed. When we consider D3-branes in the orbifold Y = S 5 / Γ the surface (4.2) must be invariant under the orbifold action. This requires the invariance of function F up to an overall factor: where w I ( g) are diagonal elements in the matrix form (1.1) of Γ and w( g) is a phase factor depending on g. Note that w as well as w I is an element of the dual group (2.2). A continuous deformation of the worldvolume is realized by changing the function F continuously without violating the relation (4.3). Because w takes discrete values, it does not change under continuous deformations. Namely, w is a topological invariant associated with the D3-brane worldvolume in Y. We call w the wrapping number (although it is not an integer but an element of Γ).

JHEP10(2019)243
We will later focus on brane configurations described by monomial holomorphic functions F (X I ). If F (X I ) is such a function, by setting X = Y = Z = 1 in (4.3) we obtain This relation directly gives the wrapping number of the configuration given by the function F . In particular, we obtain w = w I for F = X I . Namely, the diagonal components in the orbifold action (1.1) are nothing but the wrapping numbers of three brane configurations X I = 0 which will play a central role in the following analysis.
If the orbifold has fixed loci there are shrinking two-cycles along the loci. The above definition of the wrapping number does not take these shrinking cycles into account. D3branes wrapped around these cycles are regarded as tensionless strings in the (2, 0) theories. The analysis in [7] showed that the index of the gauge theory in the large N limit is reproduced as the contribution of the gravity multiplet in the bulk and the tensor multiplets localized on the loci, and the tensionless strings do not contribute the index.
To fully specify a D3-brane classical configuration we need to give not only the shape of the worldvolume but also the gauge field on the brane. Let us consider a D3-brane wrapped on X I = 0. The topology of the worldvolume is C = S 3 / Γ, which may have fixed loci. The flat connection on C * is classified by H 1 (C, U(1)) = Γ. 6 For general monomial functions the equation F (X I ) = 0 gives a set of three stacks of D3-branes. A stack wrapped on X I = 0 contains n I D3-branes, and U(n I ) gauge theory is realized on it. We can specify a flat U(1) n I ⊂ U(n I ) connection on each stack by a set of n I elements of Γ: These holonomy variables are shifted by the global 1-form symmetry associated with the background NS-NS two-form field, which we call T -symmetry. This symmetry is parameterized by t ∈ H 1 (Y * , U(1)) = Γ, 7 and it acts on the holonomy variables as h I,i → th I,i .

Wrapped D3-branes and baryonic charges
Based on the results of the analysis of S-fold theories in [22] we assume that the finite N corrections can be obtained from D3-branes wrapped around particular three-cycles X = 0, Y = 0, and Z = 0. The relation between wrapped D3-branes and baryonic operators was first pointed out in [11]. In general a wrapped D3-brane has non-trivial topology, and then we need to take 6 If C contains an An+1 type singular locus it should be treated in this equation as a curve with mod n linking. This is because if we move an endpoint of an open string attached on the worldvolume around the singular locus n times the resulting open string can be continuously deformed back into the original configuration.
7 Y * is defined from Y by removing fixed loci.

JHEP10(2019)243
account of the holonomy of the gauge field on the brane [12]. We can easily generalize the relation between wrapped branes and baryonic operators to a general orbifold. Let us consider a D3-brane wrapped over one of the cycles X I = 0. The flat connection on the worldvolume is specified by one holonomy variable h ∈ Γ. Namely, the classical configuration is specified by I and h. On the gauge theory side we have the corresponding baryonic operator 8 (4.7) This identification is consistent with the fact that on the gauge theory side we can identify the T -symmetry as the shift symmetry of the periodic quiver diagram. (We assume that the marginal deformation parameters such as gauge couplings are appropriately tuned.) An element t ∈ Γ acts on the basis vectors of the baryonic charges as b g → b tg . As a simple check we can easily confirm that the dimension dim B I,h = N is reproduced as the mass of the D3-brane. For the orbifold Y = S 5 / Γ the common radius L of the S 5 and the AdS 5 is given in terms of the D3-brane tension T D3 by The mass of a D3-brane wrapped over one of the cycles X I = 0 is M = 2π 2 L 3 T D3 /|Γ| = N/L, and the dimension of the corresponding operator is LM = N . We can also show that the wrapped brane carries the correct R-charges (R I , R I+1 , R I+2 ) = (N, 0, 0), 9 due to the coupling to the background R-R flux [13]. We can generalize this correspondence to the configuration described by a general monomial function F (X, Y, Z) = X n X Y n Y Z n Z . On the gravity side F (X, Y, Z) = 0 gives n X D3-branes wrapped on X = 0, n Y D3-branes wrapped on Y = 0, and n Z D3-branes wrapped on Z = 0. Let h I = {h I,1 , . . . , h I,n } be the set of n I elements of Γ representing the holonomy on the branes wrapped on X I = 0. A corresponding operator is (4.9) By reading off the baryonic charges from (4.9) we obtain the map from the brane configuration with the holonomy on it to the baryonic charges: should be replaced with B I,h = det(Φ c I h→c I w I h I ), where cI are elements of Γ depending on the choice of the origin. We cannot fix cI only by the consistency with the T -symmetry. For each I the choice of cI is a matter of convention. However, relative values amomg cI are physical. In this paper we focus on configurations with a single wrapped brane and the corresponding indices are not affected by cI . For this reason we simply neglect this ambiguity in the following. 9 We treat I = X, Y, Z as a cyclic variable. For example, if I = Z, (RI , RI+1, RI+2) = (RZ , RX , RY ).

JHEP10(2019)243
Now let us establish the relation between the wrapping number w ∈ Γ and the baryonic charge b ∈ Λ B . In general Λ B is larger than Γ and we cannot simply identify b and w. This difference comes from the holonomy degrees of freedom on wrapped D3-branes. To relate the baryonic charges to the wrapping number we need to eliminate the information associated with the T -symmetry from the baryonic charges. Let us define the reduced lattice Λ red B from Λ B by forgetting about the information associated with T symmetry: where the equivalence relation In fact, we can easily show ∆b g T ∼ 0 and hence Up to the T -symmetry the sum of two basis vectors of the classically conserved baryonic charge lattice Λ 0 B is given by This relation implies that the map b e − b g → g is an isomorphism, and Λ red B ∼ = Γ. We can simply identify Λ red B with the group of the wrapping number. The homomorphism ϕ from Λ 0 B to Λ red B is given by The index for a specific w is given by where the summation is taken over b satisfying ϕ(b) = w. Just like I b and I b we can extract I w from the total index I by the ζ-expansion. For this purpose we impose the constraint corresponding to (4.12). Then ζ b depends on b through the corresponding wrapping number w = ϕ(b) and we denote it by ζ w . We can rewrite (2.24) as Because of (2.17) I is invariant under the overall phase rotation ζ g → e iθ ζ g . With this redundancy we can set ζ e = 1. Then the relation (4.17) means that the map g → ζ g is a homomorphism. Namely, we can regard the fugacities satisfying (4.17) as an element of the orbifold group Γ. ζ w in (4.18) is nothing but the pairing of ζ ∈ Γ and w ∈ Γ.

Index from wrapped D3-branes
In the next section we calculate indices on the both sides of the duality for each sector specified by the wrapping number w ∈ Γ and compare them. For distinction, we denote the index on the gauge theory side by I gauge w and the index on the gravity side by I AdS w . On the gravity side the index I AdS w is given by where I D3 F,h is the index of D3-branes with the worldvolume F = 0 and the holonomy h on it. The summation is taken over monomials F satisfying F (w I ) = w and holonomies h for each wrapped D3-brane configuration F = 0. If F is a monomial of order k the index is O(q kN ). If k = 1 the holonomy varianble is a single element of Γ, while if k ≥ 2 h is a set of k elements of Γ.
Because we are interested in the leading finite N corrections we focus on the four monomials F = 1, X, Y , and Z, which give contribution to the sectors with w = e, w X , w Y , and w Z , respectively. For w = e, w I the index I w is of order O(q 2N ) or higher, and we will not pay attention to them.
Let us first consider the case that e and w I are all different.
The leading contribution to I AdS The leading contribution to I AdS w X is given by F = X, which gives a single D3-brane wrapped over X = 0. The index on the gravity side is given by where k is the lowest order of monomials F = X satisfying F (w I ) = w X . The worldvolume theory on the single D3-brane is a U(1) gauge theory and no fields couple to the holonomy. Therefore, the index I X,h is independent of h and we can replace the summation over h by the factor |Γ|:

JHEP10(2019)243
On the gauge theory side the baryonic charges of sectors giving leading contribution are b g − b w X g . Due to the T -symmetry the index I gauge bg−bw X g does not depend on g, and this degeneracy gives the factor |Γ|: (4.25) By comparing (4.24) and (4.25) we obtain We also have similar relations for w = w Y and w = w Z . The index for a D3-brane wrapped on X I = 0 was calculated in [22] for S-fold theories, and the result can be used for orbifolds, too, by simply replace the S-fold projection operator used in [22] by the orbifold projection defined in (3.6). I D3 X I ,0 are given by where the single-particle indices are the same as what are used in [22]: . (4.28) (We set the U(1) A fugacity η introduced in [22] to realize the S-fold action to be 1.) The relations in (4.22) and (4.26) hold only when e, w X , w Y , w Z are all different. If some of them coincide, the corresponding index I w is given as the sum of their contributions.
Let us consider the case that one of w I , say w X , is e. This is the case for the orbifold C × C 2 /Z n with Z n acting on Y and Z. In this case the two configurations given by F = 1 and F = X contribute to I e : I e = I KK (1 + |Γ|I D3 X,0 + · · · ). (4.29) w X = e means that the wrapping number of the configuration X = 0 is trivial and we can unwrap the brane by a continuous deformation. Correspondingly, the single particle index includes a term proportional to q −1 , which we call the tachyonic term.
When we calculate the plethystic exponential such a tachyonic term is treated as follows.
Although the physical meaning of this is not clear it was found in [22] that this gives the correct leading finite N corrections for N = 4 SYM, and we simply assume that this works for orbifolds, too.

JHEP10(2019)243
If w = w Y = w Z = e, the system has SU(2) flavor symmetry mixing Y and Z. This is the case if all dots in the toric diagram except ones at corners are aligned on a median. The two configurations Y = 0 and Z = 0 contribute to the w sector; In this case, two configurations Y = 0 and Z = 0 belong to the same homology class, and we can continuously deform them to each other via the intermediate configurations given by bY + cZ = 0. (4.33) The coefficients b and c are the homogeneous coordinates of the configuration space CP 1 .
In this case each of the single particle indices PI D3 Y and PI D3 Z includes a term proportional to q 0 . We call such terms "zero-mode terms." The plethystic exponentials of these zero-mode terms give the fractional factors These are combined into the character of the SU(2) flavor symmetry: where χ N is the u(2) character defined in (1.7). If w = w X = w Y = w Z are satisfied and the system has the flavor symmetry SU(3). This is the case for C 3 and C 3 /Z 3 . For C 3 /Z 3 the index of the w-sector is given as the sum of three contributions up to higher order corrections: w O(q 0 ) and O(q N ) higher order Table 1. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C 3 /Z 3 . The shaded row is irrelevant to the leading corrections.

Examples
In this section we confirm that by the prescription given in the previous section we can obtain the leading finite N corrections for N = 2. The results for N = 3 are shown in appendix A.

C 3 /Z 3
If n = 1 the relations among w I are The toric diagram and the quiver diagram are shown in figure 2. There are three sectors labeled by w = e, w X , w 2 X (see table 1). Let us first consider the sector with trivial wrapping w = e, which we call the mesonic sector. The relations (5.2) show that any monomial w n 1 X w n 2 Y w n 3 Z of order n 1 + n 2 + n 3 = 3k (k = 0, 1, 2, . . .) is e, and the corresponding function F (X, Y, Z) = X n 1 Y n 2 Z n 3 gives a brane JHEP10(2019)243 configuration F = 0 contributing to the index of the mesonic sector. Such a monomial with the lowest order is F = 1. The corresponding D3-brane configuration, F = 0, does not contain any D3-brane, and the corresponding index is I D3 1 = 1. We show this monomial in the second column in table 1. There exist infinite number of higher order monomials which give higher order contributions to the index. An example of such higher order monomials with the lowest order, X 3 , is shown in the third column in table 1. This higher order contribution is O(q 3N ), and hence we have the relation I D3  The comparison of (5.4) and (5.5) shows that the relation (5.3) certainly holds for N = 2. Next, let us consider the sectors with w = e. We call such sectors baryonic sectors. In the example of C 3 /Z 3 we have two baryonic sectors with w = w X and w = w 2 X . For the former the corresponding functions F with the lowest order are F = X, F = Y , and F = Z, which are shown in the second column in table 1. For each of D3-brane configurations X = 0, Y = 0, and Z = 0, the holonomy on the brane is specified by an element of Γ = Z 3 . We should sum up contributions of three D3-brane configurations with all possible holonomies. They gives the index of order O(q N ). We also have higher order functions F . Among such functions an example with the lowest order, X 4 , is shown in the third column in table 1. The D3-branes corresponding to such higher order functions give O(q 4N ) contributions to the index. Therefore, the D3-brane contribution for the e = w X sector is I D3 . On the gauge theory side, the operators corresponding to the brane configuration X I = 0 are det(Φ g→w X g I ) (g ∈ Z 3 ), which have the baryonic charges b = b g − b w X g . Three sectors labeled by these baryonic charges give the same index due to the T -symmetry, and we obtain the relation I gauge w X = 3I be−bw X + O(q 4N ). The discrepancy of order O(q 4N ) comes from the existence of other baryonic charges corresponding to w = w X . By combining the relations above we obtain Let us confirm this relation by the numerical calculation for N = 2. The result of the localization for the left hand side of (5.6) is The right hand side of (5.6) is In this equation "· · · (identical terms) · · · " mean that the terms are the same as the corresponding terms in (5.7). Comparison of (5.7) and (5.8) confirm the relation (5.6).
For the other baryonic sector with the wrapping number e = w 2 X the index is O(q 2N ). We do not analyze such sectors in this paper.

C 3 /Z 5
Let us move on to the case with n = 2: X = C 3 /Z 5 . In this case w X = w Y = w Z and the system has a flavor symmetry SU(2). The toric diagram and the quiver diagram are shown in figure 3. There are five sectors labeled by the wrapping numbers w = w k Y (k = 0, 1, 2, 3, 4). The order of indices of these five sectors can be easily read off from table 2. The sector with w = e is the mesonic sector and sectors with w = w Y and w = w 3 Y = w X are baryonic sectors whose indices are O(q N ). The other sectors with JHEP10(2019)243 w O(q 0 ) and O(q N ) higher order Table 2. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C 3 /Z 5 . The shaded rows are irrelevant to the leading corrections. w = w 2 Y and w = w 4 Y , which are shaded in table 2, give indices of O(q 2N ) and we do not discuss them.
The expected relations are Let us confirm these relations for N = 2. For the mesonic sector with w = e the result of localization on the gauge theory side is The formula for the Kaluza-Klein modes (3.1) gives I grav = · · · (identical terms) · · · + 0q 6 + O(q 7 ). (5.11) (5.10) and (5.11) are consistent with the first relation in (5.9). For the baryonic sector with w = w Y the result of localization with N = 2 is and on the gravity side we have For the baryonic sector with w = w 3 Y the result of localization with N = 2 is 2 ), (5.14) w O(q 0 ) and O(q N ) higher order Table 3. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C 3 /Z 7 . The shaded rows are irrelevant to the leading corrections. and on the gravity side we have I D3 X,0 = · · · (identical terms) · · · − u 2 χ 5 + uχ 3 q 7 + O(q

C 3 /Z 7
For n = 3 the orbifold is X = C 3 /Z 7 with the orbifold group generated by diag(ω −2 7 , ω 7 , ω 7 ). The toric diagram and the quiver diagram are shown in figure 4. There are seven sectors specified by the wrapping numbers w = w k Y (k = 0, 1, 2, . . . , 6) (see table 3). We are interested in the mesonic sector with w = e and two baryonic sectors with w = w Y and w = w 5 Y . The other sectors shaded in table 3 give higher order corrections which we do not discuss here. The expected relations are Let us confirm these relations for N = 2.

JHEP10(2019)243
For the mesonic sector with w = e we have on the gauge theory side , (5.17) and on the gravity side These are consistent with the first relation in (5.16). For the baryonic sector with w = w Y the index on the gauge theory side is and on the gravity side we have These are consistent with the second equation in (5.16). For the baryonic sector with w = w 5 Y on the gauge theory side we have and on the gravity side we have I D3 X,0 = · · · (identical terms) · · · + u 2 −χ J 2 + χ J 4 q 8 + O(q 9 ). (5.22) These are consistent with the third equation in (5.16).

The other C 3 /Z 7
There is the other Z 7 orbifold different from the one studied in the last subsection. Its orbifold group Z 7 is generated by diag(ω 7 , ω 2 7 , ω 4 7 ) and the dual group Γ is characterized by the relations The corresponding toric diagram and the quiver diagram are shown in figure 5. Among the seven sectors specified by w = w k X (k = 0, 1, 2, . . . , 6) we are interested in one mesonic sector w = e and three baryonic sectors w = w X , w = w Y = w 2 X , and w = w Z = w 4 X (see w O(q 0 ) and O(q N ) higher order

JHEP10(2019)243
The latter three equations in (5.24) for the baryonic sectors are related by cyclic permutations of X, Y , and Z, and we explicitly check only the second in (5.24). On the gauge theory side the index is calculated as and on the gravity side we obtain These two results are consistent with the second relation in (5.24).
5.3 C × C 2 /Z n (n = 2, 3, 4, . . .) Let us consider the series of orbifolds X = C × C/Z n (n = 2, 3, 4, . . .). We assume that Z n acts on Y and Z as The dual group Γ is characterized by the relations If n = 2, the relation w Y = w Z holds and the system has an SU(2) flavor symmetry while if n ≥ 3 w Y = w Z and there is no continuous flavor symmetry. An important feature of these orbifolds is the existence of the fixed locus Y = Z = 0, and we need to take account of the contribution of the tensor multiplets. Although these systems preserve N = 2 supersymmetry the indices we calculate in this section are not N = 2 indices because of the difference in the treatment of the trace parts in the N = 1 vector multiplets and adjoint chiral multiplets. We include the trace contribution for the adjoint chiral multiplets while we do not for the vector multiplets. w O(q 0 ) and O(q N ) higher order Table 5. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C × C 2 /Z 2 .
We can numerically confirm the relations in (5.31) for N = 2. For the mesonic sector with w = e the calculation on the gauge theory side gives On the gravity side we obtain For the baryonic sector with w = w Y the localization on the gauge theory side gives w O(q 0 ) and O(q N ) higher order Table 6. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C × C 2 /Z 3 . and the analysis on the gravity side gives

C × C 2 /Z 3
Let us consider the n = 3 case with the orbifold X = C×C 2 /Z 3 . The toric diagram and the quiver diagram are shown in figure 7. There are three sectors specified by w = e, w Y , w 2 Y . The sector with w = e is the mesonic sector and the ones with w = w Y and w = w 2 Y are baryonic sectors. See table 6.
The expected relations are

JHEP10(2019)243
We numerically confirm these relations for N = 2. For the mesonic sector the calculation on the gauge theory side gives On the gravity side we obtain These two results are consistent with the first equation in (5.36). Two baryonic sectors with w = w Y and w = w 2 Y are related by the Weyl reflection swapping Y and Z, and hence we check only the relation for w = w Y explicitly. The localization on the gauge theory side gives and on the gravity side we obtain

C 3 /Z 4
Let us consider the orbifold group Z 4 generated by diag(ω −2 4 , ω 4 , ω 4 ). The toric diagram and the quiver diagram are shown in figure 8. w O(q 0 ) and O(q N ) higher order Table 7. Wrapping numbers w ∈ Γ and the corresponding monomials F (X I ) are shown for the orbifold X = C 3 /Z 4 . The shaded row is irrelevant to the leading corrections.
The dual group Γ = Z 4 is characterized by the relations There are four sectors specified by w = w k Y (k = 0, 1, 2, 3) (table 7). The expected relations are Let us confirm these relations for N = 2.
For the mesonic sector with w = e, the calculation on the gauge theory side gives On the gravity side we obtain I grav I tensor X = 1 + 2u 2 q 2 + 5u 4 − χ 2 u + χ 4 q 4 + 0q 5 + O(q 6 ).

JHEP10(2019)243
For the baryonic sector with w = w Y the gauge theory calculation gives 2 ), (5.45) and the analysis on the gravity side gives For the baryonic sector with w = w 2 Y the calculation on the gauge theory side gives 2 ), (5.47) and the analysis on the gravity side gives These results are consistent with the third relation in (5.42).
There are four sectors specified by w = e, w X , w Y , w Z (table 8).

JHEP10(2019)243 6 Summary and discussions
We investigated the finite N corrections to the superconformal index for quiver gauge theories realized on D3-branes put in abelian orbifolds X = C 3 / Γ. We focused on the leading corrections starting from q N and gave a prescription to calculate them on the gravity side as the contribution of wrapped D3-branes. Based on the analysis in [22] for S-fold theories we assumed that the corrections are reproduced as the contributions of D3-branes wrapped over three particular three-cycles, X = 0, Y = 0, and Z = 0 in the internal space Y = S 5 / Γ. The D3-brane worldvolume of such a configuration is expressed by using a monomial function F (X, Y, Z) as F (X, Y, Z) = 0. We established the relation between the D3-brane wrapping number and the baryonic charges. We gave a simple prescription to obtain baryonic charges from the function F supplemented by holonomy variables. We compared indices calculated on the both sides of the AdS/CFT correspondence for each sector specified by the wrapping number of D3-branes.
In this paper we focused on the leading finite N corrections, and did not study subleading corrections starting at q kN with k ≥ 2 depending on the sector we considered. As far as we have checked for examples we found the complete agreement up to higher order terms of sub-leading corrections. This fact strongly suggests that our method is correct at least for configurations with a single D3-brane.
There are many open questions. An interesting and important task is the test of our method for the sub-leading corrections. In particular, analytic expressions for the exact Schur index is known for a few theories with finite N [29,30] and it would be very nice if we could reproduce the complete index on the AdS side. Our assumption claims that they should be reproduced by D3-brane configurations including two or more wrapped D3-branes. Naively, the index for such configurations would be calculated as the index of non-abelian gauge theories realized on the coincident wrapped D3-branes. In addition, if the configuration contains branes wrapped on two different cycles open strings connecting them may contribute to the index. Further investigation is required to establish the method to calculate these contributions.
In the examples X = C × C/Z n we found that D3-branes wrapped on a topologically trivial cycle give non-trivial contribution to the index. These branes are often referred to as giant gravitons, and such non-trivial contribution to the index had been also found in N = 4 SYM in [22]. On such branes a "tachyonic" degree of freedom lives and there exists a q −1 term in the corresponding single-particle index. We formally deal with the plethystic exponential of this term by the relation Pexp(q −1 ) = −q Pexp(q). (6.1) The negative sign appearing in this relation correctly reproduces the finite N correction due to "the absence of giant gravitons." At present, unfortunately, we have no physical understanding of this mysterious agreement. One potential way to justify our method rigorously is to use localization. Let M be the configuration space of D3-branes in the internal space Y. The three configurations, JHEP10(2019)243 X = 0, Y = 0, and Z = 0 are fixed points of a rotation of M induced by an isometry of Y. If we can extend M so that the path integral of the fields on D3-branes are treated as the integration over the extended space, and we can localize the integral to the points corresponding to the branes wrapped on the three cycles by using a localization theorem, then it would give a justification of the prescription we proposed in this paper.
In this paper we studied abelian orbifolds. There seems no obstacles to extend the method to non-abelian orbifolds and toric Sasaki-Einstein manifolds. In particular, for the latter class of manifolds the dual gauge theories are generically strongly coupled and the agreement of the superconformal index is highly non-trivial. There exists systematic prescription [31] to obtain the field contents of quiver gauge theory from the toric data of the dual internal space Y, and it is in principle possible to calculate the index by using the localization formula. On the gravity side, as a naive analogy of the prescription for orbifolds we expect that the finite N correction may be reproduced as the contribution of D3-branes wrapped on three-cycles in Y. Although it may not be straightforward to calculate the single-particle indices for excitations on such three-cycles, once we obtain them we can easily test the proposal. We hope that we can come back to this problem in near future.

JHEP10(2019)243
A.3 C × C 2 /Z 2 The first relation in (5.31) is confirmed for N = 3 by comparing the following results.