Partition functions on 3d circle bundles and their gravity duals

The partition function of a three-dimensional N=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}=2 $$\end{document} theory on the manifold ℳg,p, an S1 bundle of degree p over a closed Riemann surface Σg, was recently computed via supersymmetric localization. In this paper, we compute these partition functions at large N in a class of quiver gauge theories with holographic M-theory duals. We provide the supergravity bulk dual having as conformal boundary such three-dimensional circle bundles. These configurations are solutions to N=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}=2 $$\end{document} minimal gauged supergravity and pertain to the class of Taub-NUT-AdS and Taub-Bolt-AdS preserving 1/4 of the supersymmetries. We discuss the conditions for the uplift of these solutions to M-theory, and compute the on-shell action via holographic renormalization. We show that the uplift condition and on-shell action for the Bolt solutions are correctly reproduced by the large N limit of the partition function of the dual superconformal field theory. In particular, the Σg × S1 = ℳg,0 partition function, which was recently shown to match the entropy of AdS4 black holes, and the S3 ≅ ℳ0,1 free energy, occur as special cases of our formalism, and we comment on relations between them.


JHEP05(2018)116
Recently, there has been much progress in performing exact, nonperturbative computations for superconformal field theories (SCFTs) on curved manifolds via the technique of supersymmetric localization (see the review [1] and references therein). Such methods have greatly developed in the past several years, providing a tool to study a wide variety of SCFTs in various dimensions and backgrounds, leading to non-trivial tests of holography and other known dualities.
In particular, recently these techniques have been successfully applied to the computation of the partition function of three-dimensional superconformal Chern-Simons-matter theories on Σ g ×S 1 in presence of background magnetic flux for the R-and flavor symmetries through the Riemann surface, Σ g [2][3][4][5]. By performing a partial topological twist [6] on Σ g , one obtains the so-called "topologically twisted Witten index" [3]. This was shown in [7] to reproduce in the large N limit 1 the macroscopic entropy of supersymmetric magnetic AdS 4 black holes in theories of 4d FI-gauged supergravity. These black hole configurations, first found in [8], consist of M2-branes wrapped around Σ g , and thus they implement the partial topological twist for the QFT describing the low-energy dynamics of the M2-branes.
This recent success led to several extensions and developments. First of all, the entropy matching was performed on more general supergravity backgrounds, including dyonic black holes [9], black hole configurations arising from massive IIA supergravity truncations [10,11] and solutions with hyperbolic horizon [12]. Moreover, unexpected relations were discovered between the topologically twisted index on S 2 × S 1 and the corresponding partition function on S 3 in the large N limit [13]. Specifically, the computation of the twisted index involves, as an intermediate step, the computation of the twisted superpotential, or Bethe potential, as a function of the flavor fugacities [3]. Then this quantity was shown to coincide, for a suitable mapping of parameters, with the large N limit of the S 3 partition function of the same N = 2 theory [13]. Given that the partition function on S 3 of three-dimensional superconformal theories, such as the ABJM theory [14], has been extensively studied, and has connections to entanglement entropy and the F -theorem [15][16][17], it is natural to ask if such a correspondence has a deeper meaning.
In parallel with these developments, a new class of partition functions for general 3d N = 2 gauge theories was computed in [18] utilizing a three-dimensional uplift of the 2d A-model [19]. These partition functions are defined on the manifold M g,p , a U(1) bundle of Chern degree p ∈ Z over a Riemann surface Σ g , where g ∈ Z ≥0 denotes the genus of the Riemann surface. This set of manifolds includes in particular the three-sphere S 3 and the product spaces Σ g × S 1 M 0,1 S 3 , M g,0 Σ g × S 1 . (1.2)

JHEP05(2018)116
Thus, these partition functions include both the topologically twisted index of [2][3][4][5] and the round S 3 partition function of [20][21][22] as special cases. This then provides a natural framework to address the relation between the topologically twisted index and S 3 partition function. At the same time, constructing explicit supergravity backgrounds whose boundary is such a circle bundle and computing their renormalized on-shell action provides a viable holographic check for these field theory computations. In more detail, the partition function on M g,p can be computed by a sum over supersymmetric "Bethe vacua," [23], where the index I runs over the set S BE of vacua of the theory. Here m i and s i are, respectively, real masses and fluxes for background flavor symmetry gauge fields, and F I , H I , and Π I i are certain functions appearing in the 3d uplift of the A-model, described in section 2 below. We will argue that, for a class of quiver gauge theories with holographically dual M-theory descriptions, in the large N limit this sum can be approximated by a single dominant term, I dom , and we find the result: leading to a very simple dependence on the geometric and flux parameters. We find the partition function exhibits the expected N 3/2 scaling, and reproduce and generalize the results of [4,7] in the case p = 0. However, we find that a large N solution exists only under certain conditions on the mass and flux parameters. In the case of M 0,1 = S 3 , these conditions differ from those under which previous large N computations of the S 3 partition function were carried out, e.g., in [15], and we comment on this discrepancy in section 4 below. We reproduce this result holographically, by providing supergravity backgrounds having boundary M g,p in the framework of minimal N = 2 U(1) gauged supergravity. Such solutions can be embedded locally in 11d on 7-dimensional Sasaki-Einstein manifolds. We construct Euclidean regular solutions which preserve 1/4 of the supersymmetries and have appropriately quantized magnetic flux. Starting from the analysis of [24][25][26], we find that the boundary can be filled with multiple gravity configurations, with different topology. In particular, for the boundary S 3 case we can have regular "NUT" solutions, with topology R 4 , and for S 3 /Z p one finds mildly singular NUT/Z p solutions. On the other hand, for general M g,p we find regular "Bolt" solutions, with topology O(−p) → Σ g . The different topology has non-trivial consequences for the uplift of these solutions. Indeed, while there are no requirements for the NUT solution to lift to M-theory, the Bolt uplifts to eleven dimensions only for certain values of g and p, depending on the geometrical properties of the internal Sasaki-Einstein 7-manifold. Interestingly, the same constraints are recovered in the field theory computation by setting all fluxes equal, thus reproducing the universal twist which corresponds to minimal gauged supergravity. 2

JHEP05(2018)116
The computation of the on-shell action of these two distinct bulk solutions is obtained via standard techniques of holographic renormalization. The resulting on-shell action for the NUT configuration coincides with the free energy of the corresponding theory on S 3 . The renormalized on-shell action for Bolt solutions is instead of the form Vol(S 7 ) Vol(Y 7 ) (4(1 − g) ∓ p) , (1.5) with the additional constraint ± p + 2(g − 1) = 0 (mod I(Y 7 )) , (1.6) where I(Y 7 ) is the Fano index of the internal 7-manifold. In particular, for g = 0 we retrieve the results of [25].
We are able to show that, for this solution in minimal gauged supergravity, the on-shell action in the gravity side matches with the partition function of the corresponding field theory (1.4) − log Z Mg,p = I Bolt , (1.7) as expected. In case of trivial fibration, p = 0, our formulas (1.4) and (1.5) find agreement with those of [27]. In this particular case the on-shell action of the Euclidean solution coincides with the entropy of supersymmetric 1/4 BPS black holes with constant scalars and higher genus horizon. 3 Along with the matching with the Bolt solutions, we study the relation between the S 3 partition function as computed by [15] and the result we obtain for the M g,p partition function, in light of the result of [13]. In particular, we elaborate on how the interesting relation between the extremal value of the twisted superpotential and the large N partition function on S 3 , discovered in [13], fits in our framework by relating these both to the partition function on the lens space S 3 /Z 2 .
The main text of the paper is organized as follows: in section 2 we provide the details of the computation of the large N partition function of a class of N = 2 3d quiver gauge theories on M g,p , focusing on the example of the ABJM theory. In section 3 we describe Euclidean minimal gauged supergravity solution whose boundary is M g,p , and examine their supersymmetry properties, along with their moduli space for regularity. Moreover, we compute the on-shell action via holographic renormalization. In section 4 we show the matching between the renormalized on-shell action of the Bolt solutions and the partition function of the dual field theory for the ABJM theory. In section 5, we consider more general quiver gauge theories, including the V 5,2 theory, and describe the truncation to minimal supergravity for these theories, obtaining a generalization of the universal twist of [27]. In section 6, we discuss the relation between the twisted superpotential and S 3 partition function observed by [13], and relate these to the lens space partition function. Finally in section 7, we discuss some open issues and future directions. Several appendices complete this paper, and they are devoted to the construction of the explicit Killing spinor

JHEP05(2018)116
for the supergravity solutions, to the description of their moduli space, and to the explicit details for the computation of the partition function Z Mg,p .
2 M g,p partition function at large N We start in this section by discussing the computation of the M g,p partition function for 3d N = 2 field theories. We first describe the supersymmetric background on M g,p and review the computation for general, finite N theories, and then turn to the large N computation for a class of U(N ) quiver gauge theories with M-theory duals.

Supersymmetric background on M g,p
Following [18], we consider manifolds which are U(1) bundles over a Riemann surface, Σ g , with p ∈ Z the Chern-degree of the bundle, which we take to be non-zero in this subsection. On this space we take the following metric: 4 ds 2 = β 2 (dψ − C(z,z)) 2 + 2g zz dzdz , (2.1) where z,z are local coordinates on Σ g , with g zz the metric on Σ g , ψ ∼ ψ + 2π is a coordinate along the U(1) fiber, which has length 2πβ, and C is a locally defined 1-form on Σ g , satisfying: We may define a Killing vector K = 1 β ∂ ψ pointing along the U(1) fiber, or equivalently, a 1-form: To preserve supersymmetry on this space, we must turn on additional fields in the background supergravity multiplet [30,31]. These include the R-symmetry gauge field, A R µ , a scalar, H, and a vector, V µ . These lead to the Killing spinor equation: We find a solution on the above geometry, in the local coordinates above, when we take: In [18] the scalar parameter κ was set to zero, but for comparison to the supergravity background below we will take κ = −2pβ; these choices will lead to the same Killing spinor ζ. Here the last term in A R µ corresponds to a contribution from a flat connection, and we will describe this in more detail below.

JHEP05(2018)116
Although in principle we may take an arbitrary smooth metric on Σ g , and an arbitrary connection C subject to (2.2), for concreteness, and to compare to the bulk supergravity solution, we will consider constant curvature metrics and connections: where the connection a(θ, φ) is given by: In all cases, φ ∼ φ + 2π is an angular coordinate. For g = 0, θ ∈ [0, π] and this is the usual round metric on S 2 . For g = 1, we identify θ ∼ θ + 1, obtaining the flat, rectangular metric on the torus. For g > 1, the second and third terms in (2.6) describe the metric on the hyperbolic plane, H 2 , and we form Σ g by taking an appropriate quotient of the hyperbolic plane by a Fuchsian subgroup [32], with a fundamental domain D g . In all cases we have normalized the metric on Σ g so that vol(Σ g ) = π. The connection a has curvature da proportional to the volume form on Σ g , and satisfies: In this constant curvature background, the background supergravity fields are H = −ipβ and V µ = 0, and the R-symmetry gauge field is: Here γ is a flat connection with Chern class generating Z p ⊂ H 2 (M g,p ), so that this gauge field has torsion flux g − 1 (mod p). Then the Killing spinor equation, (2.4), becomes: Background gauge fields. In addition to the fields in the background supergravity multiplet, we may include background gauge multiplets coupled to the global symmetries of the theory. If i = 1, · · · , r H runs over a basis of the Cartan of the flavor symmetry group, H, then we may turn on background gauge multiplets V i in configurations labeled by: where σ i is the real scalar in the background gauge multiplet, and the gauge field A i is given by: where η and a are as in (2.3) and (2.8), and π * is the pullback along the projection map π : M g,p → Σ g . Explicitly, we can write, for p = 0: with γ as in (2.9), so that the gauge field has torsion flux s i (mod p).
The M g,p partition function we compute below will be functions of the parameters m i and s i , depending holomorphically on the former. Note that shifting: does not change the connection A i (modulo gauge transformations), and we will see below this is an invariance of the partition function.

Computation of M g,p partition function
In this subsection we review the computation of the M g,p partition function for a general, finite N , 3d gauge theory, in preparation for computing the partition function at large N in the next subsection. We refer to [18] and references therein for more details.
There are two equivalent methods to compute the M g,p partition function. First, one may couple the UV action of the 3d gauge theory to the supersymmetric background on M g,p discussed above, along the lines of [31], and use localization to reduce the path integral to a finite dimensional integral. As shown in [18], one arrives at the following integral formula for a theory with gauge group G with rank r G and Weyl group W : Here u a , m a , a = 1, · · · , r G are holonomies and fluxes, respectively, for the gauge field, and similarly for m i , s i , i = 1, · · · , r H for background gauge fields coupled to the flavor symmetry group, H. The functions F , H, Π, and H depend on the field content of the theory, and are defined in (2.26) below. The integral is taken over a compact contour C JK , the so-called "Jeffrey-Kirwan" contour [33,34]; we refer to [18] for the precise definition.
In the case p = 0, this contour integral can be deformed into one over a non-compact "Coulomb" contour, C Coul , and one obtains the equivalent formula: where the fluxes m a now takes values in Z p rather than Z. Roughly speaking, 5 C coul = iR r G is the imaginary slice in the complex u plane, and upon making the identification u → iσ, 5 More precisely, this statement is true only when suitable conditions on the R-charges of the chiral multiplets are satisfied; more generally the contour may be deformed to pass around certain poles coming from the contributions of the chirals. See [18] for more details.

JHEP05(2018)116
one recovers the usual integral formula for the S 3 partition function [20][21][22] in the case g = 0, p = 1. We return to the connection to previous computations of the S 3 partition function in sections 4 and 6 below. The second method starts from the observation that the M g,p partition function is computed by a certain 2d topological quantum field theory (TQFT) on the base space, Σ g , of the fiber bundle. Specifically, this TQFT is the "A-twist" of the 2d N = (2, 2) theory obtained by compactifying the 3d gauge theory on a circle and studying the low energy effective action. Then, on general grounds, we expect the partition function to be given by a sum over the supersymmetric vacua of the 3d theory on a circle. Explicitly, one finds: involving the same functions as in (2.15). We will describe this formula in more detail below. These two methods can be shown to be equivalent for an arbitrary 3d gauge theory [18]. For the purpose of taking the large N limit, we will utilize the sum-over-vacua formula, (2.17), for the remainder of this section. However, despite these formulas being equivalent for finite N , there are some subtleties in relating the large N limit obtained by the two methods. We will return to this issue in section 4.
Let us now describe the formula (2.17) in more detail, starting by reviewing the compactification on S 1 and the vacuum structure of the resulting system.
Twisted superpotential and Bethe vacua on R 2 × S 1 β . Given a 3d N = 2 gauge theory, we may place it on R 2 × S 1 β , obtaining at low energies an effective 2d N = (2, 2) description. Then the vacuum structure of the theory is determined by the "Bethe equations" [23]: where W(u a , m i ) is the effective twisted superpotential of this effective 2d N = (2, 2) system, defined in (2.20) below. Here we define: 19) where σ is the real scalar in the dynamical gauge multiplet, and A t the component of the gauge field along S 1 β , and we expand them in a basis of the Cartan subalgebra of G, and similarly for the parameters, m i , i = 1, · · · , r H , for background gauge multiplets coupled to the flavor symmetry, H. Note that u a ∼ u a + 1 and m i ∼ m i + 1 due to large gauge transformations around S 1 β . The twisted superpotential, W(u a , m i ), depends on the matter content and UV Lagrangian of the 3d theory. We consider a general theory with gauge group G of rank r G , and chiral multiplets in some representation of G. We expand the chiral multiplets in weights of G, such that they have charges Q a α , a = 1, · · · , r G , α = 1, · · · , M , where M is the dimension of the space of chiral multiplets. The chiral multiplets also have charges, S i α , i = 1, · · · , r H ,under the global symmetry group, H, which may be restricted by superpotential terms in the Lagrangian. Finally, we allow Chern-Simons terms k ab , k ij , and k ai ,

JHEP05(2018)116
for the gauge, flavor, and mixed CS terms, respectively. Then the twisted superpotential is given by the following function of u a and m i : (2.20) Here W χ (u) is the contribution of a chiral multiplet, regulated with a level − 1 2 CS term to ensure gauge invariance, given by: (2.21) With this definition of the chiral multiplet contribution, the bare CS levels, k ab , k ij , and k ai , must all be integers to ensure gauge invariance. We note that the twisted superpotential in general has branch cuts, and is only defined modulo shifts of the form: for n a , n i , n ∈ Z. However, one can check that (2.18) is invariant under these shifts, and is a polynomial equation in the variables x a = e 2πiua and µ i = e 2πim i . For non-abelian gauge theories, one must discard solutions which are not acted on freely by the Weyl group, W , as supersymmetry is broken at these putative vacua, and we consider the remaining solutions up to Weyl symmetry. Then we define the set of Bethe vacua as: ∂W ∂u a = 1, a = 1, · · · , r G , w ·û a =û a , ∀w ∈ W /W . (2.23) Ingredients in the M g,p partition function. With this background, let us now return to the computation of the M g,p partition function. The coupling of this theory to the curved background of M g,p depends on the choice of a U(1) R symmetry, which is used to perform a partial topological twist along the Σ g directions. Since we will introduce a non-trivial flux for this R-symmetry, we must pick the R-charges, r α , of the chiral multiplets to be integers, so that they live in a well-defined vector bundles over M g,p . Given such a choice of R-symmetry, we define the "effective dilaton," Ω(u a , m i ): where k R a , k R i , k RR ∈ Z correspond to contact terms involving the R-symmetry. Here the second line is the contributions from the W -bosons of the gauge group, and the sum is over the weights of the adjoint representation of G. Note the W -bosons do not contribute to W.

JHEP05(2018)116
Then, as argued in [18], the M g,p partition function is given by the following sum over Bethe vacua: where the "fibering operator," F , "handle-gluing operator," H, and "flux operators," Π i , are defined in terms of the twisted superpotential, W, and effective dilaton, Ω, defined above: Such a formula arises due the topological invariance along Σ g , which implies that the operations of gluing a handle to Σ g , adding a unit of flux for the S 1 fibration, or adding a unit of flavor symmetry flux, are all implemented by local operators, H, F , and Π i , respectively, giving rise to the simple formula (2.25). To take a simple example, the M g,p partition function of a single chiral multiplet is: where m, s, and r are the its mass, flavor symmetry flux, and R-charge, respectively, and: (2.28) Note this depends only on the combination: In other words, a shift of the R-charge, r → r + c, is equivalent to a shift of the flavor symmetry flux, s → s + c(g − 1), and reflects a mixing of the R-symmetry with this flavor symmetry. If we take: (m, ) → (m + 1, + p) , (2.30) and use the difference equation: we see the partition function, (2.27), is invariant. This is consistent with the invariance of the background gauge field, (2.13), under this shift of parameters, and reflects the fact that, for p = 0, the fluxes are torsion, and take values in Z p . For a general gauge theory, we may explicitly write the summand in (2.25) as: 6 . (2.32) 6 Here for simplicity we work in a basis of the flavor symmetry group where mi corresponds to the mass of the ith chiral multiplet, and we take the flavor and R-symmetry CS terms to vanish. We may also treat the contribution of the vector multiplets as that of an R-charge 2 chiral multiplet in the adjoint representation of G.

JHEP05(2018)116
On-shell twisted superpotential. We may conveniently construct the terms in the sum above using the "on-shell" twisted superpotential and effective dilaton, defined by: where the index I runs over S BE . Here there is a branch cut ambiguity in defining W I , but this is partially fixed by imposing: which is a stronger condition than (2.18), and fixes the freedom to shift W by n a u a . Then one has: One can check that the remaining branch cut ambiguities in W I and Ω I drop out of these expressions, and they are well-defined. Then we may construct the partition function as:

Large N computation
We will be interested in computing this partition function for a large N gauge theory. In the case p = 0, this problem was studied in [7]. There they found that, although the number of Bethe vacua, |S BE |, grows with N , in many cases there is a single vacuum, with index I dom , whose contribution is dominant compared to all other terms in (2.36). 7 When this occurs, we expect that (2.36) may be approximated as: Note, in particular, that the partition function has a very simple dependence on the geometric parameters, g and p, and the fluxes, s i . In cases where the theory has a holographic dual, this suggests the holographic free energy has a similar simple dependence on these parameters, which is rather non-trivial. Below we will verify this relation holds quite generally.

U(N ) quiver gauge theories
In this section we will focus on the ABJM model [14], which we describe in more detail below. This is a special case of a more general class of U(N ) quiver gauge theories. The ingredients in the computation of the twisted superpotential and M g,p partition function for these quivers is very similar, so we describe these general ingredients in the next few subsections, returning to a more detailed analysis of these theories in section 5. Specifically, the class of theories we will discuss, following [13,15], consists of N = 2 quiver gauge theories with several U(N ) gauge factors, labeled by an index α = 1, · · · , n. We allow bifundamental chiral multiplets connecting two gauge groups, (anti-)fundamental chiral multiplets in a single gauge group, and Chern-Simons levels, k α , for the αth gauge group. However, we impose the following restrictions: • The sum of all Chern-Simons levels is zero: (2.38) • For each gauge node, α, there is a superpotential constraint which imposes that, for all bifundamental chiral mutiplets with a leg in this node: where Q I is the charge of the Ith such bifundamental chiral multiplet under any flavor symmetry, and r I is its R-charge. Here adjoint chirals are counted twice in the sum.
• The number of bifundamental chiral multiplets entering a node is the same as the number exiting the node.
• The total number of fundamental and anti-fundamental chiral multiplets in the quiver are equal.
These restrictions are to ensure the theory has a well behaved M-theory dual description at large N , with a characteristic N 3/2 scaling of the number of degrees of freedom. For such quivers, following [7], we will take the following large N ansatz for the eigenvalues u α a : u α a = v α a + iN 1/2 t a , a = 1, · · · , N, α = 1, · · · , n . (2.40) In the large N limit the eigenvalues become dense, and we may parameterize them by the continuous variable t, defining: and corresponding functions, v α (t).
Our strategy in the rest of this section is as follows. First, we compute the twisted superpotential at large N , using the above ansatz, and find the eigenvalue distribution JHEP05(2018)116 which extremizes it. Then, as in (2.37), we may assume that the dominant contribution to the Bethe sum computing the M g,p partition function is determined by this extremal distribution. Thus, we evaluate the summand in (2.37) at this extremal distribution to compute the leading behavior of the M g,p partition function.

The twisted superpotential at large N
We start by reviewing the computation of the twisted superpotential at large N , as first computed in [7] for the ABJM theory, and studied for more general quivers of the above type in [13]. 8 For a given set of eigenvalues, u α a , approximated by the distributions ρ(t) and v α (t) above, we may compute the value of the effective twisted superpotential at these eigenvalues as a functional: Let us briefly summarize the various ingredients in the functional W[ρ, v α ], as computed in [7,13]. We review the derivation of these ingredients in appendix C. First, the CS terms, which satisfy α k α = 0, contribute: Next, a bifundamental chiral multiplet connecting the αth and βth groups contributes: where we defined δv = v α −v β , m is the mass of the bifundamental, and we have introduced the following notation for the "fractional part," [u], of a complex number u: [u] = u − n, n ∈ Z, such that 0 < Re([u]) ≤ 1 . (2.46) The function g(u) is given by: Here we have imposed the constraint (2.39), which implies: We also impose that the total number of incoming and outgoing edges at each node in the quiver are equal. Note then that the (v α ) 3 terms contributed by W bif will cancel, so that the functional is in general quadratic in the v α . 8 Let us state the relation between the notations used here and those used in [7,13]. We have: These changes propagate into the large N ansatz, e.g., t them = 2πtus, ρ them = 1 2π ρus, etc..

JHEP05(2018)116
In addition, there are contributions to W from a bifundamental chiral that are subleading in N , but whose derivatives with respect to the v α get large near special points in parameter space, and so they affect the extremization of W. Specifically, these contributions become important when δv I + m I =n I ∈ Z for some bifundamental chiral multiplet, with index I. Then if we write: for some positive function Y I (t), one finds an additional "tail contribution:" Finally, an (anti-)fundamental chiral multiplet contributes: with the + (−) sign for a fundamental (anti-fundamental) chiral.

Extremal value and the ABJM theory
As described above, we will need to find the eigenvalue distribution which extremizes W.
To do this we vary the functional W[ρ, v α ] with respect to ρ(t) and the v α (t). We also include a Lagrange multiplier term, iN 3/2 µ dtρ − 1 , to impose correct normalization of ρ. The solution is in general defined piecewise, bounded by points where δv I + m I becomes an integer, after which δv I becomes locked to this value to leading order, varying at subleading order as in (2.49).
Let us consider as our main example the ABJM theory [14]. This has U(N ) k ×U(N ) −k gauge group, with two bifundamentals in the (N,N ), with masses m 1,2 , and two in the (N , N ) representation, with masses m 3,4 . We assume k > 0; the case with k < 0 can be obtained by exchanging the two gauge groups. This theory includes a quartic superpotential which imposes the following constraints on the masses: The twisted superpotential is periodic under m i → m i + 1 (up to branch jumps), and so depends only on the fractional part of the masses, [m i ]. The functional we obtain for ABJM turns out to depend only on δv = v 1 − v 2 . We look for solutions with: 9

JHEP05(2018)116
Then we simply have [±δv + m i ] = ±δv + [m i ] (where here and below, for "±" we take + for i = 1, 2 and − for i = 3, 4), and then the functional becomes: The extremal distribution was first derived in [7]. Note that (2.52) imposes that where (here we take µ > 0, so these are in ascending order): . (2.57) Then one computes the extremal value of W as: There is a similar solution when i [m i ] = 3, related by m i → 1 − m i and W → −W. However, for i [m i ] = 2, we see the quadratic term in δv vanishes, and we do not find a solution.

M g,p partition function at large N
Next we consider the functional computing the M g,p partition function. Specifically, we compute the contribution to log Z Mg,p from a Bethe vacuum which is, approximately at large N , given by a distribution of eigenvalues u a corresponding to the functions ρ and v α ,

JHEP05(2018)116
as in (2.40). Once we have found the dominant such eigenvalue distribution, as above, we may plug this in to this functional to compute the leading behavior of the partition function.
Here we list the various ingredients, which are derived in appendix C. The Chern-Simons terms contribute: A bifundamental chiral multiplet with mass m, flavor flux s, and R-charge r contributes: where: (2.61) Here n = u − [u] is the integer part of u, and = s + (g − 1)(r − 1), as in (2.29). A vector multiplet (which, recall, does not contribute to W) contributes as above with δv + m → 0 and → (g − 1), giving: Here we have imposed the constraints in (2.39). Once again, since we impose the number of incoming and outgoing edges at each node are equal, the cubic terms in δv will cancel, and this gives an expression quadratic in the v α . We also have contributions from the tail regions, where δv I + m I ≈n I ∈ Z for some I. Here we find: where I is as above, and the sum is over all such tail regions. Finally, for an (anti-)fundamental chiral multiplet, we have: We note that the expressions above are invariant under: where we recall shifting m → m + 1 entails shifting the corresponding integer part, n → n + 1. This reflects the fact that the fluxes are defined modulo p, as in (2.30).
Let us now return to the ABJM example. Then the partition function is a function of the masses, m i = [m i ] + n i , flavor fluxes, s i , and R-charges, r i , where the latter enter in JHEP05(2018)116 the combination i = s i + (g − 1)(r i − 1). Due to the quartic superpotential, these satisfy the constraints: (2.66) Then, one finds the functional computing the M g,p partition function is given by: plus the contribution of the tails. Now let us plug in the extremal solution in (2.55), which, recall, required i [m i ] = 1. Then we must impose: Plugging in the eigenvalue distribution found above and evaluating the integral, one eventually obtains the following simple result: Let us make a few comments about this formula. First, in the case p = 0, where M g,p=0 ∼ = Σ g × S 1 , this reproduces the results of [4,7]. 10 Next, recall this formula only applies when i [m i ] = 1; when i [m i ] = 3, we find another solution related by m i → 1 − m i , explicitly: For i [m i ] = 2, we do not find a solution. We will return to this point in section 4 below. Also, note that the result (2.69) has the expected form (2.37): where here W I dom is given by W ABJM ext in (2.58), and: (2.72) We will see in section 5 that this relation holds also for more general quiver gauge theories. 10 To compare to their results, one makes the identifications in footnote 8, as well as n them

JHEP05(2018)116 3 The supergravity dual
In this section our aim is to find supergravity solutions whose boundary is the manifold M g,p , a circle bundle over a closed Riemann surface Σ g , which can be locally uplifted in M-theory. 11 According to the AdS/CFT dictionary, we expect the on-shell action of these solutions, suitably renormalized, to match with the M g,p partition function of the dual field theory, as computed above. We will return to this comparison in the next section.

Minimal N = 2 gauged supergravity
Our starting point is minimal N = 2 four-dimensional gauged supergravity, whose bosonic action reads [35,36] where G 4 is the four-dimensional Newton's constant and l is the AdS radius, related to the cosmological constant via Λ = −3/l 2 . We work in Euclidean signature. The gravitino supersymmetry variation is where is a Dirac spinor and γ µ are the generators of Cliff(4, 0) and so they satisfy {Γ a , Γ b } = 2g ab . We follow here closely the conventions of [25]. The Einstein's equations coming from (3.1) read and Maxwell's ones are d F = 0 . (3.5) We will restrict our analysis to a set of solutions where Σ g has constant curvature, and to configurations with a real metric. Solutions to the system of equations of motion (3.4)-(3.5) have been obtained in [37] and they have the following form 12 In this case, κ denotes the curvature of Σ g : κ = 1 for S 2 , κ = 0 for R 2 and κ = −1 for H 2 . The 2d area element dΩ 2 κ reads The gauge field has this form In these solutions, M is the mass parameter and r is the radial coordinate. τ parameterizes a circle fibered over a 2-dimensional constant curvature surface Σ g spanned by the coordinates θ and φ. The fibration is due to the presence of the NUT parameter, which we denote by s because of its relation with the squashing of the U(1) fiber relative to the base. 13 Solutions of this kind for κ = 1 were first discovered by Taub [39] and Newman, Unti and Tamburino [40], hence the name Taub-NUT. Their structure and thermodynamics properties were later studied in [38] and [37]. In the latter, (Lorentzian) solutions with NUT charge and planar and hyperbolic horizon were analyzed as well.
In the asymptotic limit r → ∞ the metric approaches where we have defined ψ = τ /(2s). In other words, the boundary is a circle bundle over Σ g , and in the particular case for which f (θ, φ) = cos θdφ and ψ is periodic with period ∆ψ = 4π, the boundary is a squashed 3-sphere with squashing parameterized by 4s 2 /l 2 . We will see in the next subsection that, once one takes into account the appropriate compactifications, the boundary metric in eq. (3.12) coincides with the one considered for the field theory computation, eq. (2.6), up to a rescaling of coordinates. The bulk solutions we have described in this section can be of the type AdS-Taub-NUT and AdS-Taub Bolt, depending on the value of the parameters appearing in the warp factor. These solutions are characterized by different topologies, since one of the Killing vectors has a zero-dimensional fixed point set ("nut") or a two-dimensional one ("bolt"). We discuss the requirements for the regularity for NUTs and Bolts separately below, along with the conditions of periodicity of the coordinates. We focus first on the spherical case, Σ g = S 2 , so κ = 1.

NUTs and Bolts
NUTs and Bolts with S 2 base.
NUTs. For the NUT solution the Killing vector ∂ τ has a fixed point where the S 2 has zero radius: This ensures that the Killing vector has a zero-dimensional fixed point. Moreover, absence of Dirac-Misner [41] strings constrain the period of τ to be ∆τ = 4s∆φ , (3.14) and since ∆φ = 2π, this yields ∆τ = 8sπ, or equivalently ∆ψ = 4π (see formula (3.12)). The coordinate θ goes from 0 to π. The last condition concerns the absence of conical singularities at the location of the nut, r = s and is in the following In particular, for the NUT solutions the warp factor λ(r) has a double root at r = s and the metric is defined for r ≥ s. The point r = s is a NUT-type coordinate singularity and the metric is a smooth metric on R 4 with the origin identified with r = s. For s = 1/2 the squashing vanishes, and the boundary is that of a round S 3 , and one recovers AdS 4 space. The AdS Taub-NUT space is a self-dual Einstein space: the Weyl tensor is self-dual and the gauge field has a self-dual field strength, with Q 2 − P 2 = 0.
Notice that in principle, if we allow for conical singularity, we could also take into consideration quotients of the Taub-NUT space, obtaining the AdS-Taub-NUT/Z p geometry. Quotients of R 4 suffer from conical singularities at the origin, however such backgrounds were studied in [42] and were shown to be sensible backgrounds for holography.
Bolts. For the Bolt solution the Killing vector ∂ τ has a two-dimensional fixed point, so the only condition is that, at a radius r b > s, with r b a single zero of λ(r). The absence of conical singularities at the location of the bolt requires This condition ensures that the metric near the Bolt takes the form [25] where we have defined a new radial coordinater which parameterizes the distance to the bolt in this way:

JHEP05(2018)116
The metric of the Bolt is defined for r ≥ r b and the topology is The boundary is a biaxially squashed lens space S 3 /Z p . At the specific point in parameter space we obtain the quaternionic Eguchi-Hanson [43] solution. Among the class of Bolt solutions, the Eguchi-Hanson solution provides an example of regular self-dual Einstein space.
NUTs and Bolts with Σ g base. The procedure to determine regularity of NUTs and Bolts carries over similarly in the planar and hyperbolic case, along the lines of [37]. We first discuss the case in which the coordinates are not compact.
For the planar case, we look for self-dual solutions and we impose that the warp factor λ(r) admits two coincident roots at r = s. The warp factor then assumes the form In the hyperbolic case, as already noticed in [37], NUT solutions do not exist: once imposing self-duality, the largest roots of the warp factor always lie at a radius r s > s. Bolt solutions, in contrast, are present in both planar and hyperbolic case.
In order to find solutions with M g,p boundary we need to compactify the coordinates by imposing suitable boundary conditions. It is instructive to consider first the metric at the boundary. Consider a 3d metric of the form where ψ is compactified with period ∆ψ: ψ ∼ ψ + ∆ψ. Here θ and φ are coordinates on the compact Riemann surface Σ g , which we obtain by a suitable quotient of R 2 (in the case g = 0) or H 2 (in the case g > 1), and D g is the fundamental domain of this group action, as in section 2.1. Then in order to have a well-defined fiber bundle we must impose for p ∈ Z. To see this is necessary, let v ∈ [0, 1] be a coordinate parameterizing the boundary of D g . Then if we define: we see that: Then in a neighborhood of ∂D g we may shift ψ → ψ + h(v), and (3.25) implies this is a single-valued coordinate transformation. This implies dψ → dψ + a, and so eliminates a in (3.22), and so after this coordinate transformation we may consistently identify the boundaries of ∂D g by the group action and obtain the compact space Σ g .

JHEP05(2018)116
For our solution, the fundamental domain for the compactification, D g , is chosen as in [4], so that vol(Σ g ) = 4π(g − 1) for g > 1 (3.26) and D dθdφ = 4π for g = 1. 14 Therefore, given that a = θdφ for the torus, and a = cosh θdφ in the higher genus case, (3.23) yields the following conditions Let's now turn to the Bolt solutions mentioned before. We can have Bolt solutions with topology O(−p) → Σ g imposing r b > s, λ(r = r b ) = 0, with r b a single zero of λ(r) and using the compactification above to get a genus g > 0 base manifold Σ g .
The absence of conical singularities at the location of the bolt requires This condition ensures that the metric near the Bolt takes the form (we take κ = −1 for simplicity) where again we have defined a new radial coordinater which parameterizes the distance to the bolt in this way:r ∝ √ r − r b . The metric (3.29) is regular atr = 0 if p 2|g−1| ∆ψ = 2π, hence we have retrieved exactly eq. (3.27).
Therefore, a M g,p boundary with g > 0 can be filled by the Bolt 15 solution, with topology O(−p) → Σ g and Euclidean time period (3.27).

Supersymmetry properties
The supersymmetry properties of the spherical NUTs and Bolts were analyzed in detail in [25], where it was shown that configurations satisfying preserve 1/2 of the supersymmetry, while those with are 1/4 BPS. In [25] the explicit form of the Killing spinor was provided as well. We will be interested in the 1/4 BPS ones. For our purposes, we wish to generalize the BPS conditions to solutions whose boundary is M g,p .
14 Notice here the slightly different conventions with respect to section 2.1, in order to match the conventions of [4]. Here vol(Σg) = 4π(g − 1) and the Chern number p appears in the periodicity of ψ. 15 We could as well consider toroidal NUT solutions. They however present conical singularities like the S 3 /Zp NUTs, and we do not consider them in our holographic checks.

JHEP05(2018)116
AdS 4 solutions with NUT charge and Σ g horizon in minimal N = 2 gauged supergravity in Lorentzian signature were presented in [44]. In the same paper, the integrability conditions for supersymmetry were analyzed, and necessary conditions for supersymmetry to be preserved were given. For the 1/4 BPS case they read Lorentzian solutions of this form were studied also in [45]: in the latter, it was shown that (3.32) is a sufficient condition for preserving 1/4 of supersymmetry, 16 following the procedure of [48]. We work here in Euclidean signature and, taking a pragmatic approach, we explicitly construct the Killing spinor for the Euclidean solutions under consideration. Therefore, we map the integrability condition of [44] via and we find that a necessary condition for the solutions to be 1/4 BPS is Notice that (3.34) reduces to (3.31) for κ = 1. In appendix A we construct the Killing spinor for the solutions satisfying the conditions (3.34), generalizing the procedure of [25], whose notation we follow closely, to the higher genus case. As one can see in appendix, the Killing spinor has only radial dependence, so that the compactification necessary to obtain a compact Riemann surface is allowed without breaking supersymmetry. We notice that the signs in this equation (3.34) can be reabsorbed with a change of coordinate {r → −r, ψ− → −ψ, φ → −φ} hence without loss of generality we focus on the set of parameters This condition denotes a four-dimensional subspace of BPS solutions, parameterized by p (Chern class of the bundle), g (genus of the base Riemann surface), Q (electric charge) and s (squashing parameter). As we will see in the next section, regularity imposes additional constraints on these quantities, and in particular, we will see that there are different bulk fillings for the same boundary data. These 1/4 BPS configurations with M g,p asymptotics will be the gravitational configurations of our interest. When solving the BPS equations, once the r component is satisfied, the remaining components of the Killing spinor equations reduce to the new minimal rigid supersymmetry equation for a 3d spinor υ. Here γ i are the Pauli matrices i = 1, 2, 3, ∇ (3) is the covariant derivative associated to the 3d boundary metric (3.22). The (global part of the) asymptotic gauge field is 16 See also [46] and [47], where the supersymmetry enhancement in the self-dual case was analyzed.

JHEP05(2018)116
which is mapped to the background R-symmetry gauge field (2.9) by taking into account that β = 2s/p, and the fact that, as already mentioned, the boundary metric in this section is defined with the following conventions with respect to the field theory one (2.6). Equations of the form (3.36) were shown to naturally arise on the boundary of supersymmetric AdS configurations in theories of gauged supergravity, by [30,31]. As noticed in [25], one needs to supplement the expression for the gauge field in (3.37) by a flat connection. The Killing spinor in appendix A is computed in a rotating frame: if one had to compute its expression in a static frame, a e iψ dependence would appear. The flat connection introduces an opposite phase which cancels with the latter, so the global form of the Killing spinor to be independent of ψ (see appendix A).

Moduli space of solutions
This subsection is devoted to the study of the moduli space of the solutions. The space of the configurations with g = 0 (κ = 1) was discussed already in [25], we extend it here to the cases with g > 0 as well. We focus on the cases This class of 1/4 BPS configurations include NUTs and Bolts, depending on the value of parameter Q. As already anticipated, for κ = 1 the former solutions (NUTs) have as conformal boundary a squashed S 3 , and for the latter (Bolts) the boundary is a squashed lens space S 3 /Z p . We start by analyzing the roots of λ(r). These are The largest root of λ(r), denoted with r 0 , can be either r 2 or r 4 , depending on the specific values of s and Q. We denote the two solutions with r + and r − respectively: We note that the requirement for a Taub-NUT solution is r 0 = r + = r − = s and gives Q = ± 1 2 (4s 2 − κ). For this value, the warp factor λ(r) has two coincident roots. This solution is well-defined for all value of the squashing s > 0 and its boundary is a p = 1 fibration over a 2d constant curvature Riemann surface Σ g : in particular, when g = 0, we retrieve the squashed three-sphere.
For the Bolt, imposing regularity near r = r 0 amounts to [37] where we have defined p = p |g − 1| for g = 1 , p = p for g = 1 . (3.43) Since r 0 > s for the Bolt solution, the numerator of (3.42) is positive. The derivative at the denominator is as well, since r 0 is the largest root of λ(r) and λ(r) goes as r 2 for r → ∞. Therefore in our considerations we will restrict to the case sp > 0, which encompasses both s > 0, p > 0 and s < 0, p < 0. We need to analyze the two cases (3.41) separately. If the largest root is r 4 , imposing r 0 = r + yields the following condition on Q: while if we impose r 0 = r − we obtain two possibilities Plugging these values of Q into r ± we obtain the following expressions for r + (using eq. (3.44)): and r − (using (3.45)) We can see then that we obtain up to four different branches denoted by the relations (3.44) and (3.45). 17 We will see shortly that each of the two solutions in (3.44) gives the same value I + for the on-shell action, which differs with respect to the one I − obtained for both the solutions in (3.45). Hence we collectively denote both signs in (3.44) with "Bolt + " and (3.45) with "Bolt − ." In order for these branches to exist, all the following conditions needs to be met: 1. the function under the square root f ± = (16 s 2 ± p) 2 − 128 κ s 2 should be positive, 2. respectively, r ± should be indeed the largest root of the warp factor λ(r) and 3. r ± > s in order to have a Bolt.
These requirements in particular set constraints on the value of the squashing parameter s. Indeed regular Bolt solutions might be present only in a limited interval of squashing parameter s: for instance, for g = 0 the Bolt + solutions for p = 1, 2 are regular only for 17 In the case g = 0 these branches were analyzed in [25], but we appear to find another additional Bolt+ solution with respect to their analysis, corresponding to the solution with Q − + . This branch gives a new solution only for p = 1, and the latter, despite being a genuinely different solution, has the same range of existence as the Bolt+ with Q + + . Hence the analysis of the on-shell action for g = 0 carries out unchanged with respect to [25]. More details in the appendix.

JHEP05(2018)116
a finite interval s ∈ [0, s 1 ], while they always exist for p ≥ 3 [25]. A detailed analysis of the moduli space for the different g cases is treated in appendix B. We remind the reader that spherical NUT solutions exist for all s > 0. Moreover, taking quotients of the latter generates the mildly singular NUT/Z p solutions mentioned in section 3.2, which have lens space boundary and are defined for all values of s.
We end this section by computing the total flux for the Bolt solutions. Indeed the field strength F of the Bolt has a non-trivial magnetic flux through the Bolt surface Σ g at r = r 0 , which lies inside the bulk (it is the point where the Bolt geometry caps off). The flux can be computed as where the ± refers to the two classes of Bolt + and Bolt − defined above. Notice that the form of the metric and field strength is fixed demanding regularity of the configurations, and the flux computed in (3.48) follows. Its value satisfies where q = p mod 2, which is consistent with A being a spin c gauge field. However, in order for the uplift to 11d to be consistent the flux at the Bolt is subject to a further condition.
In particular, depending on the geometry of the internal Sasaki-Einstein seven manifold Y 7 , the uplift is possible only for certain values of p: this is discussed in the next section.

Uplift to 11d
The 11d uplift ansatz for solutions of (Lorentzian) minimal gauged 4d supergravity was found in [49], where it was shown that any supersymmetric configuration of such 4d theory uplifts locally to a solution in eleven dimensions described by the following metric and 4-form G: where the ds 2 4 is the metric element in the four-dimensional spacetime, with volume form vol 4 , ξ is the contact form of Y 7 and Kähler Einstein base metric ds 2 6 . Here 4 is defined with respect to the 4d metric and 11 with respect to the eleven-dimensional one. The seven-dimensional metric on Y 7 satisfies R αβ = 6g αβ . The four-dimensional Newton constant G 4 is  (3.54) and the radius R * is . For a solution of 4d supergravity to uplift to a properly defined M-theory configuration, some parameters should be appropriately quantized. Because of this, there are some differences among the conditions we need to impose for NUTs and Bolts [25].
The Taub-NUT solution has topology R 4 and the gauge field A is globally a oneform on R 4 : the twisting ξ + 1 2 A is trivial. Therefore every NUT solution uplifts without restrictions to eleven dimensions. There is no quantization condition on G, which is itself a globally defined four-form on R 4 × Y 7 [25].
For the Bolt solution, characterized by topology O(−p) → Σ g , the situation is different. When Y 7 is a regular Sasaki-Einstein manifold, Y 7 is a U(1) bundle over the six-dimensional Kähler Einstein manifold B 6 . In this case ξ = dη + σ where dσ is the Ricci form on B 6 divided by four. Regularity imposes that the period of η is 2πI/4k, where k is a positive integer (which is the level of the Chern-Simons theory) and I = I(B 6 ) ∈ {1, 2, 3, 4} is the Fano index of B 6 . We will restrict here to the case k = 1 for simplicity. In order for the U(1) bundle to be well-defined at the Bolt itself, the flux should satisfy [25] 4 For the NUT and NUT/Z p , instead, the bundle is trivial. Eq. (3.55) for Bolts is a condition involving p and g and will be crucial point for the matching of the on-shell action with the partition function of the corresponding N = 2 superconformal field theory. For instance, the uplift on the round seven sphere S 7 (k = 1 and I(S 7 ) = 4), making use of (3.48), dictates the following constraint: which leads to p ± 2(g − 1) = 0 (mod 4) (3.57) Notice that there are cases for which Bolt solutions uplift for all p: for example the manifold Y 7 = M 3,2 , since I(M 3,2 ) = 1 [25]. It is worth mentioning another specific example we will focus later on: the manifold V 5,2 , whose Fano index is I(V 5,2 ) = 3. In this case the uplift condition (3.55) using (3.48) dictates the following constraint: which leads to p ± 2(g − 1) = 0 (mod 3) (3.59)

JHEP05(2018)116
We will see in section 4 the consequences of the condition (3.57) and (3.59) for the matching with the dual field theory computation.

On-shell action via holographic renormalization
We are now ready to compute the renormalized on-shell action for the branches of solutions that we have found in the previous sections, satisfying (3.35). In order to compute such quantity, we need to plug in the solution into the action (3.1), and perform the integral with extrema of integration r 0 (radial coordinate where the spacetime caps off) to r → ∞ (asymptotic boundary). However, this leads to a divergent quantity, therefore we need to adopt an appropriate renormalization scheme.
Techniques of holographic renormalization [50][51][52][53] have been extensively developed in the past years. Following these, one first regulates the action via the introduction of a radial cutoff r inf , which is sent to infinity after adding appropriate counterterms. The counterterms are function of the boundary 3-metric h ij and they read [52] where the term in Θ is the Gibbons -Hawking boundary term, defined as where n µ is a unit vector normal to the boundary and R 3 (h) is the Ricci scalar of the induced boundary metric. In what follows, we divide the action in two parts, in particular we single out the contribution of the vector fields, Plugging in the solution, we obtain the on-shell action which we denote with I. We compute where we have made use of The extremum of integration r 0 and the value of the gauge field contribution differ between NUTs and Bolts, whose on-shell actions are computed separately below.
Bolt branches. For the Bolt solution, instead, the expression for I F 2 is  We see from these results that the parameter s, which quantifies the squashing of the sphere, dropped out of the final result for both NUT and Bolts. Hence for this class of 1/4 BPS solutions the free energy does not depend on the squashing parameter, as was found in [25,54]. For g = 0 the on-shell action reduces to that found in [25]. For p = 0, instead, expressions (3.73)-(3.74) reduce to the entropy of supersymmetric 1/4 BPS black holes in gauged supergravity, found by [44,55].
Before proceeding further, we'd like to take into consideration another possible background with lens space S 3 /Z p boundary, which is the mildly singular AdS-Taub-NUT/Z p , also named "orbifold NUT." The free energy of this configuration is 1/p times that of the NUT solution [42]. However, if we want to have the same boundary data as the g = 0 Bolts, we need to consider orbifold NUTs where with the addition of ± p 2 − 1 units of magnetic flux, which is the same as that of the Bolt ± . Therefore the on-shell action will be supplemented by a term of the form F 2 [25] in the following way These orbifold configurations are acceptable solutions. However, as shown in the plot 3 of appendix B, their on-shell action is always higher or equal to the one of the Bolt +

JHEP05(2018)116
solutions, except in the case p = 2 where their flux vanishes and their free energy coincide. Hence when Bolts and NUT/Z p both exist, the Bolts are the favored configuration due to their lower free energy and we expect their on shell action to match the field theory computation.
We end this section by reminding the reader that the uplift of the solutions in 11d imposed the quantization condition (3.55) which constrains the values of allowed p, depending on the geometry of the internal Sasaki-Einstein manifold Y 7 . Reinstating the factors of N 3/2 in the above formulas, according to formula (3.53), we obtain , (3.76) . (3.77) Notice that the Bolt formula (3.77) reproduces for p = 0 the on-shell action of supersymmetric black holes of minimal gauged supergravity [55], which was computed in [27]. In the latter work (and more recently in [28,29] for the non-minimal case) it was moreover shown that the free energy of such configurations coincides with (minus) their entropy. The on-shell action computed here corresponds to an ensemble where the magnetic fluxes and the electric chemical potentials are kept fixed [56]. In other words, it refers to a canonical ensemble for the magnetic charges and a grand-canonical one for the electric ones. Since the partition function Z Mg,p computed in the field theory describes supersymmetric states with the same, fixed magnetic charges and it contains a sum over all electric charge sectors with fixed chemical potentials m i [9], it is natural to compare the expressions obtained in this sections (3.76)-(3.77) with those achieved in the field theory in section 2. The next section is devoted to this comparison.

Truncating to minimal supergravity
In general, the partition function we computed in section 2 is dual to a configuration in N = 2 gauged supergravity with non-trivial magnetic fluxes and profiles for the scalar fields. These profiles are governed by the asymptotic behavior (i.e., scalar modes) of these fields, which is determined by the background vector multiplets coupled to flavor symmetries in the dual field theory on M g,p , or equivalently, by the masses, m i , fluxes, s i , and R-charges, r i , that we choose. In order to compare to the supergravity computation of the previous section, we must suitably restrict these parameters so that the bulk vector multiplets are trivial, and we obtain a configuration in minimal supergravity.
As above, we focus on the ABJM model for concreteness; we will return to more general models in section 5. Then the bulk vector multiplets are trivial when we set all of the background vector multiplets on the boundary to be equal. Naively, this implies setting the parameters m i , s i , and r i all to be independent of i. However, since i m i = 0, this is only possible if m i = 0, which leads to a trivial solution. To get further, we recall JHEP05(2018)116 that, due to the shift symmetry (2.30), which allows us to shift the m i by integers, only the quantities [m i ] and −pm i + s i + (g − 1)r i are gauge invariant. Thus, we should impose that: [m i ] and n i ≡ −pn i + s i + (g − 1)r i are independent of i . (4.1) Following [7], we may interpret the fractional part of the masses, [m i ], as fixing the asymptotic behavior of the scalar fields in the bulk. 18 Then n i can be attributed to the net magnetic charge felt by the ith chiral field. Then we conjecture that (4.1) gives the truncation to minimal supergravity for the ABJM model.
Let us see what the condition (4.1) implies for large N behavior of the partition function. For the ABJM theory, recall that, for the solution with i [m i ] = 1, the expression for the large N partition function we derived in (2.69) can be written as: If we set all the n i equal, then this can only be satisfied for n i ∈ Z if we impose: Then we take n i = 1 4 (p + 2(g − 1)). Inserting these values into (4.2), we find: log Z ABJM Mg,p = πN 3/2 √ 2k 12 p + 4(g − 1) .
which satisfies: Then to set [m i ] andn i independent of i, we take: 1)) , (4.9)

Holographic comparison
At this point it is easy to compare the result obtained for the ABJM partition function on M g,p and the on-shell action of Bolt ± solutions, with uplift to M-theory on S 7 , corresponding to the case k = 1 in the field theory. The field theory result in (4.5) and (4.11) is summarized in subject respectively to the constraints (4.4) and (4.10) which read ± p + 2(g − 1) = 0 (mod 4). Moreover, the condition (4.13) maps exactly in the quantization condition (3.57) obtained in minimal supergravity. We recall that the latter condition is necessary for the uplift on S 7 to be globally defined for Bolt solutions. Notice that due to the condition (4.10) for ABJM, the Bolt + and Bolt − uplift for the same values of p. The dominant configuration, namely the one which has lowest free energy, is the Bolt + for p > 0, the Bolt − for p < 0.

NUTs and Bolts for S 3
Let us consider the special case of M g=0,p=1 = S 3 . One can see that the quantization condition (4.13) is not satisfied here, and indeed, the p = 1 Bolt ± solutions are not states in the ABJM theory, so our formalism cannot reproduce their free energy.
However, for the NUT-type solutions there are no further constraints for the uplift, therefore the AdS-Taub-NUT is a regular configuration which uplifts to 11d on S 7 . The NUT (g = 0) then provides a regular filling for the squashed S 3 and its on-shell action is (3.76) This coincides with the ABJM free energy on the S 3 at the conformal point, as computed in [15]. This lead us to a puzzle: why do we not recover the result of [15] for the S 3 partition function by our method? The precise mapping of parameters between our partition function and theirs is: 19 Z Then the superconformal R-charges, ∆ i = 1 2 , correspond to the case: To see what went wrong, recall the starting point of the computation of [15] was the integral formula, (2.16), in the special case g = 0, p = 1, and they solved this at large N by looking for a saddle point of this integrand. While their starting point, (2.16), and ours, (2.17), are equivalent at finite N , these two methods of taking the large N limit are evidently not equivalent. They can be summarized as follows: • Method 1 (the method we used above for general g and p, and used for p = 0 in [7]).
Extremize the functional, W[ρ, v α ], computing the twisted superpotential, obtaining eigenvalue distributions ρ W , v W α . Then the S 3 partition function is obtained by plugging this into the functional log Z S 3 [ρ, v α ], i.e.: • Method 2 (the method used in [15]). Directly extremize the functional log Z S 3 [ρ, v α ], computing the integrand of the S 3 partition function in (2.16), obtaining eigenvalue distributions ρ Z , v Z α . Then the S 3 partition function is obtained by: (4.20) In general, these two methods give different answers, and one must check which gives the dominant contribution to the partition function. We see that the second method gives a solution for ABJM at the superconformal point, while the first does not. However, as we will discuss in section 5 below, for certain other models, such as the V 5,2 model, both methods do lead to solutions. In general, we suggest the following interpretation of these methods: • Method 1, for general g and p, reproduces the on-shell action of the minimal SUGRA "Bolt" type solutions, whose boundary is M g,p , subject to the constraint (4.13).
• For M 0,1 = S 3 , Method 2 gives the on-shell action of the NUT-type solutions. 19 More precisely, as described in more detail in sections 5 and 6, this is true in the "physical gauge," in which case the R-charges ri can be tuned continuously, and appear in a combination mi → mi + ri ≡ ∆i.

JHEP05(2018)116
In other words, when one solution gives the dominant contribution to the free energy on the gravity side, the corresponding method gives the dominant contribution to the M g,p partition function. We note that the configurations with the lowest free energy (Bolt + for (4.13) and NUT for g = 0, p = 1) are defined for all values of the squashing. For more general g and p, one may ask if we can generalize Method 2 above by starting from the integral formula, (2.16), and finding a saddle, as done by [15] for S 3 . We comment on the case of S 3 /Z 2 in section 6 below. For g > 0, the integral formula (2.15) and (2.16) are more subtle, as one must carefully regulate the contribution the vector multiplet at points of enhanced Weyl symmetry, as discussed in [4,18]. Thus, we leave the investigation of the large N limit of this integral formula for future work. However, given the matching of the previous subsection, and the fact that a NUT-type solution does not exist for general g and p, we believe this gives strong evidence that Method 1 gives the correct large N behavior of the M g,p partition function for generic g and p.
Finally, let us mention that, for Method 1, we can naturally take p → 0: in this limit we retrieve the topologically twisted Witten index on Σ g ×S 1 studied in [3,4,7], which was used for the entropy matching of AdS 4 black holes. The Bolt solutions, indeed, present a non-trivial 2-cycle Σ g in the bulk at the location of the Bolt, which is threaded by magnetic flux. For p = 0 (which corresponds to the black hole case) this cycle extends all the way to the boundary. The NUT solutions, whose free energy is retrieved by Method 2, does not have this feature.

General quivers
In this section we discuss the more general quiver gauge theories of section 2.3.1, and some new issues that arise in taking the large N limit of their M g,p partition functions. We derive a general result relating the partition function of these gauge theories to the extremized values of the twisted superpotential, generalizing the "index theorem" of [13]. We point out a subtlety in the definition of the M g,p partition function in the presence of fractional R-charges, which is necessary for studying general quivers of the above type. We consider the explicitly the case of the case of M-theory compactified on the 7-manifold V 5,2 , and the truncation to minimal supergravity in this case. Finally, we discuss the truncation to minimal supergravity for general quivers of section 2.3.1, and describe a generalization of the "universal twist" of [27].

The M g,p partition function for general quivers
Using the ingredients provided in section 2.3, one can compute the M g,p partition function for more general quiver theories of the type discussed in section 2.3.1. It turns out that the final answer can be expressed in terms of the extremal twisted superpotential for the chosen theory a very simple form, given in (5.6) below. Let us describe this computation in more detail.
Let us first review the extremization of the twisted superpotential, which was described in [13] and computed in several examples in [57]. Let us suppose the theory has JHEP05(2018)116 a superpotential: where Φ i is the ith chiral multiplet, and we suppress color indices. Then this imposes the following constraints on the parameters appearing in the M g,p partition function: where m i , s i , and r i are the mass, flavor symmetry flux, and R-charge, respectively, for the ith chiral multiplet. In [13] it was found that a solution to the extremization of the twisted superpotential exists when we impose (translating to our notation): The extremization of W was worked out in several quivers of the above form in [57]. Suppose we can find a set of parameters satisfying (5.2) and (5.3). Then we may follow the procedure in section 2.3 to compute the large N behavior of the M g,p partition function. Namely, we plug the extremal eigenvalue distribution (5.4) into the Z Mg,p functional computed in section 2.3.4. We have worked this out in several of the examples considered in [57], and in all cases we have found the M g,p partition function has the following simple relation to the extremal twisted superpotential, W ext ([m i ]): where n i is as in (4.1), i.e.:

JHEP05(2018)116
Let us make several comments about this formula. First, in the case p = 0, it reduces to: (5.9) This agrees with the "index theorem" of [13].
Next, recall that, due to the constraint ( Fortunately, for the combination of derivatives appearing in the expression (5.6), one can check that this condition is indeed satisfied, as a consequence of (5.3) and (5.8).
We also note that for quivers of this type, the extremal twisted superpotential is typically a homogeneous function of the m i of degree 2 [13]. 20 Indeed, this was the case for the extremal twisted superpotential for the ABJM theory, in (2.58). In this case, the expression (5.6) simplifies to: In this form, the relation can be understood in terms the on-shell twisted superpotential, discussed in section 2.2. Namely, recall from (2.35) that we may write the various operators appearing in the sum over Bethe vacua in terms of the on-shell twisted superpotential and effective dilaton. In the large N limit, where the contribution from one vacuum is dominant, we may restrict our attention to the on-shell twisted superpotential and dilaton in that vacuum. Then we claim: Then inserting this into (2.35), one derives (5.12). We also note there is a symmetry of the partition function under m i → −m i , p → −p, W → −W, and k → −k, where k denotes the various Chern-Simons levels [18]. As a result there is an additional solution of the extremization of the twisted superpotential when:

JHEP05(2018)116
This implies: i q α in i = −p + 2(g − 1) , (5.15) wheren i is as in (4.7). When these conditions are satisfied, one finds: (5.16) Finally, we must mention an important caveat in the above computation. We assumed that it was possible to find a set of parameters satisfying (5.2) and (5.3), however, for certain theories one cannot find such a solution in integer R-charges r i . This forces us to consider the possibility of allowing fractional R-charges, which we turn to next.

Fractional R-charges and the R-symmetry background
So far we have always assumed that the R-symmetry used to place a theory supersymmetrically on M g,p assigned all chiral multiplets integer charges. This was necessary because our background includes an R-symmetry gauge field with non-zero flux, and so for the chiral multiplets to take values in well-defined line bundles, in general their R-charges must be integers. However, in some of the quiver gauge theories of section 2.3.1, some of the chiral multiplets necessarily have fractional R-charges, due to superpotential constraints, and so such theories naively seem beyond the scope of this partition function.
Fortunately, while integer R-charges are necessary to define the M g,p partition function for all g and p, for special choices of g and p we may allow certain fractional R-charges. However, there are certain subtleties in the definition of the M g,p partition function for fractional R-charges which will be important to understand for the large N computation. The following discussion is somewhat technical, and the main points are summarized at the end of this subsection.
To start, let us consider the case p = 0, corresponding to the partition function on Σ g × S 1 . Then the flux felt by a chiral multiplet of R-charge r includes a contribution r(g − 1). For this to give a well-defined vector bundle, we must then impose: For non-integer r, this restricts the possible values of g. However, when r is rational there are still infinitely many cases we may consider. Next take p = 0. Recall the R-symmetry gauge field, A R , is given by (2.9). This gauge field is part of a background gauge multiplet, V R with the vector component given by a complex gauge field [58]: 21 a For the background discussed in section 2.1, we can write this as: 21 We thank C. Closset for discussions about this background vector multiplet.

JHEP05(2018)116
where we recall η is the 1-form (2.3), and nγ is the flat U(1) gauge field on M g,p with first Chern class n ∈ Z p . As in (2.12), we may construct this gauge field as the pull-back of the gauge field with flux g − 1 on Σ g , i.e.: where a, defined in (2.8), satisfies 1 2π Σg da = 1. Now consider a chiral multiplet Φ with fractional R-charge, r. Then the gauge field on M g,p coupling to this chiral multiplet includes a contribution from the R-symmetry gauge field, which can be obtained by pulling back from Σ g , i.e.: Then in order for this to be well-defined, we must again impose (5.17). However, recall from (2.14) that we may equivalently write the gauge field (5.19) as: for n R ∈ Z and ν R ∈ R satisfying: The expression (5.20) corresponds to a particular choice, called the "A-twist gauge" in [18], where: For a general choice of (ν R , n R ) satisfying (5.23), the R-symmetry gauge field coupling to a chiral multiplet is 25) and the condition for this to be well-defined is: Thus, by changing the parameters (ν R , n R ), we may find backgrounds which allow more general choices of fractional R-charge. 22 For example, when g − 1 = 0 (mod p), we can take the "physical gauge" of [18]: and then we may take any r ∈ R. For a general choice of (ν R , n R ), one finds that the contribution of a chiral multiplet is given by:

JHEP05(2018)116
This reproduces (2.27) in the A-twist gauge, (5.24). Moreover, using (2.31), one can check that the partition function (5.28) does not depend on the parameters (ν R , n R ) satisfying (5.23) when the R-charges are integers. However, for fractional R-charges, the partition function will in general depend on these parameters. To be more explicit, let us write the superpotential as in (5.1), which imposes the constraints in (5.2), which we reproduce here: If we work in a general R-symmetry background, as in (5.28), it is natural to define the parameters:m and these satisfy: The conditions (5.31) carve out different subsets of the space of allowed masses and fluxes, depending on the choice of ν R and n R . However, many of these spaces can be related using the difference equation (2.31). For example, when it is possible to take all of the R-charges to be integer, that is, when it is possible to find a solution in integers γ i to: then, using the difference equation (2.31), we may shift the mass parameters and fluxes by: without changing the partition function. Then this modifies the condition (5.31) by shifting (ν R , n R ) → (ν R +1, n R +p). Since this is just a redefinition of parameters, these correspond to equivalent backgrounds on M g,p . Note that this argument relies only on the existence of some choice of integer R-charges, but the actual R-charge we pick need only satisfy (5.26). However, when (5.32) cannot be solved, the set of solutions to (5.31) for different ν R are not related by a redefinition of parameters. Then the M g,p partition function for different choices of ν R ∈ Z will give rise to inequivalent partition functions.
Changing the spin-structure. In the case of p even, we may consider yet another choice of R-symmetry background, which is: Then the R-symmetry gauge field is: As above, this does not affect the field strength of the R-symmetry gauge field. However, unlike the redefinitions above, this does change the R-symmetry gauge field itself by adding

JHEP05(2018)116
an additional flat connection, p 2 γ, with holonomy −1 around the S 1 fiber. In general, the presence of this additional holonomy means the Killing spinor equation, (2.10), no longer has a solution.
However, when p is even, H 1 (M g,p , Z 2 ) ⊃ Z 2 , and there is an additional choice for the spin structure on M g,p , differing in the sign the fermions incur as we wind around the S 1 fiber. Thus, if we accompany this shift of the R-symmetry gauge field by a change of spin structure, the Killing spinor equation is unchanged, and so this defines a valid background for p even. This ambiguity in the choice of spin structure and R-symmetry gauge field was discussed in the context of holography in [25].
Note that if we take the R-charges of all chiral multiplets to be even integers, then we find the contribution of a chiral multiplet for this R-symmetry background, (5.28), is equal to that in the A-twist gauge: where we used (2.31). Thus, for r ∈ 2Z, these backgrounds are equivalent. This follows because, for such a choice of R-symmetry, all the bosons have even R-charge, and the fermions have odd R-charge, and so the effect of the change of spin structure and Rsymmetry connection cancel for all fields, just as they did for the Killing spinor. More generally, if we can find a solution to: Then, if we shift:m this takes (ν R , n R ) → (ν R + n 2 , n R + np 2 ). In particular, if we can solve this for n = 1, all choices of ν R ∈ 1 2 Z are equivalent. However, when this is not possible, these different choices of ν R ∈ 1 2 Z in general lead to inequivalent backgrounds. We will see below that the background (5.34) will play a special role in the large N solution. 23 M g,p partition function for a general R-symmetry background. With these observations, let us return to the computation of the M g,p partition function of the quiver gauge theories of section 2.3.1. First, we note that, for finite N , the computation of the M g,p partition function of a gauge theory in a general background (ν R , n R ) is a straightforward modification of the procedure in the A-twist gauge, outlined in section 2.2. Namely, we simply make the following replacements in (2.32): (5.39) 23 Below we will consider the background (5.34) even in the case where p is not even. While this is somewhat formal, it can be justified in cases, as in (5.38), where we may use the shift symmetry, (2.31), to relate backgrounds with different values of νR.

JHEP05(2018)116
Then the M g,p partition function will depend on the parametersm i andñ i . In a general gauge, these parameters satisfy the constraints: Then the only difference between the M g,p partition function for different backgrounds, (ν R , n R ), is the set of constraints, (5.40), which we take the parameters to satisfy. Now let us return to the large N solution described in section 5.1 above. Recall that, for this solution to apply, we must find a set of parameters, m i , s i , r i , satisfying (5.2), i.e.: as well as: where: . which are clearly incompatible. However, given the discussion above, we see that the condition (5.41) is the condition appropriate for the A-twist gauge, while for a more general background, we should instead impose the constraint (5.40). Comparing to (5.42), we see that if we work in the background (5.34): where we must impose: then we may identify:m In this background, we may always solve the constraints and find a large N solution. When a solution also exists in the A-twist gauge background, this is related to the one above by a shift of mass parameters as in (5.38). This is the case, for example, for the ABJM theory. But even when no such solution exists in the A-twist gauge, we can still find a solution in the background (5.45). We will consider an example of such a theory in the next subsection. For notational simplicity, we will drop the tildes onm i andñ i below when the meaning is clear from context. To briefly summarize the above discussion, we have seen that in the presence of chiral multiplets with fractional R-charges, there are subtly different choices of supersymmetric JHEP05(2018)116 backgrounds on M g,p , which can be parameterized by the parameters ν R and n R , and which impose the constraints (5.40) on the parameters entering the partition function. The restrictions (5.3) and (5.8) which allowed us to find a large N solution for the M g,p partition in section 5.1 naturally force us to work with an R-symmetry background, (5.45), which is different than the A-twist gauge used in section 2.1. In some cases, such as for the ABJM theory, these two backgrounds are equivalent, but in general they are distinct.

Example: V 5,2 theory
To illustrate some of these features more concretely, let us now look at another example of a quiver gauge theory. We take the case where the bulk gravity theory comes from the compactification of M-theory on the manifold V 5,2 = SO(5)/SO(3), whose Fano index is I(V 5,2 ) = 3. The dual field theory, which we will call the V 5,2 theory, has a description [60] with gauge group U(N ) k=1 × U(N ) k=−1 with two pairs of bifundamental chiral multiplets, A i and B i , i = 1, 2, in the (N,N ) and (N , N ) representations, respectively, as well as an adjoint field, Φ α , for each of the two U(N ) factors, and superpotential: Note the presence of the Φ 3 superpotential term, which we saw in the previous section implies we may not work in the A-twist gauge background. Instead we work in the background (ν R = 1 2 , n R = p 2 +g −1), as discussed above. We denote the mass parameters by m A i , m B i , and m Φα , and similarly for the fluxes n. Then to find a solution at large N we impose: (p + 2(g − 1)), n Φα = 1 3 (p + 2(g − 1)) .
Note that to find a solution with n Φα ∈ Z, we must impose: This is an example of the general feature, discussed above, that in the presence of fractional R-charges, we may only define the M g,p partition function for certain choices of g and p.
The extremal twisted superpotential for this theory was computed in [57], with the result: Inserting the extremal eigenvalue distribution into the M g,p partition functional of section 2.3.4, we obtain the result: We have computed this explicitly using the prescription of section 2.3, but we note it is consistent with the general result (5.6), giving additional evidence for this conjecture.

JHEP05(2018)116
Now let us discuss the truncation to minimal supergravity. Similarly to the ABJM case, we look for solutions with: Then from (5.49), we see we must impose: In this case, the M g,p partition function simplifies to: As in the ABJM case, there is an additional solution with m → 1 − m, where: Let us now discuss the dual supergravity solutions. In the g = 0, p = 1 case, the NUT free energy, when this configuration is embedded in a 11d theory compactified on the V 5,2 manifold, is given by where we have used the fact that [61] Vol(V 5,2 ) = 27π 4 128 . (5.58) The result (5.57) agrees with the S 3 partition function of the dual quiver theory computed in [62]. Bolt solutions, instead, have free energy 59) and, in order to have a regular uplift to eleven dimensions on V 5,2 In this case, they are subject to the following constraint (see section 3.5): We see that the result (5.59) for the Bolt+ is reproduced by the M g,p partition function of the V 5,2 theory computed above in (5.55), Mg,p , and we also find the same constraint on p and g as in (5.50), with a similar relation for the Bolt− and (5.56).
However, something unexpected happens for p = 1. From the lower sign of (5.60) we see that p = 1 Bolt − are a solution of the theory, namely they lift up to M-theory compactified on V 5,2 . Their free energy is

JHEP05(2018)116
which is higher than that of the corresponding NUT solution (5.57) We are left then with the following puzzle: the NUT and p = 1 Bolt − have the same boundary data, in particular they are both regular 24 fillings for the squashed S 3 . However, our field theory result for these configurations exactly matches with (5.61). The field theory computation carried out previously in [62], which utilized a similar method to that of [15], matches (5.57). The two methods used here are the same as those discussed in section 4.3. As we observed there, it seems that the two methods reproduce different saddle points, which are realized in the dual supergravity picture in form of two solutions with different topology, namely the NUT and p = 1 Bolt configurations.
In this regard, it would be interesting to retrieve the M g,p partition function for the quiver dual to the M-theory reduction on the SE7 manifold M 32 [63,64]. In that case, there is no restriction for the uplift: Bolt solutions are present for every p. In particular, the Bolt + uplifts in the case with boundary S 3 , and has a lower free energy than the NUT solution. The S 3 free energy for this quiver was computed in [65,66] and once again it coincides with the on-shell action of the corresponding NUT solution. Given that this theory is chiral, the matrix model is not under good computational control with current methods and therefore computing the M g,p partition function is a difficult task. We leave this for future investigation.

Minimal supergravity for general quivers and the universal twist
Now consider a general quiver of the type studied in section 2.3.1. We would like to find a general analogue of the truncation to minimal supergravity discussed for the ABJM and V 5,2 models above. To do this, we make the following assumption. For general mass parameters, m i , we expect the partition function we have computed to be dual to a configuration in gauged supergravity. More precisely, in the case p = 0 it was argued in [7] that one must first extremize the partition function as a function of the mass parameters, m i , in order to compare to the supergravity solution. This reflects the attractor mechanism in the gravity theory, which fixes the behavior of scalar fields in the bulk as they approach the horizon [7,9,29] However, in the minimal case, there are no scalars, and so no such extremization is necessary.
We conjecture that a similar story holds for general p. Namely, we expect that for the case of generic m i , an extremization principle is necessary to compare the result to a solution in gauged supergravity. We leave the study of this general, non-minimal case for future work. However, we argue that the truncation to minimal supergravity corresponds to the case where we have chosen the m i such that no extremization is necessary. 24 We remind that the Bolt− solution exist only for the range of squashing s ∈ [0, ].

JHEP05(2018)116
Let us see how this works in more detail. We have argued in section 5.1 that the large N solution is given by: For a general theory, we may take an approach similar to that in [27]. As observed in [13], and as we will discuss in more detail in section 6, for a general quiver of the type above, one finds the relation: where log Z S 3 (∆ i ) is the partition function on the round sphere as a function of continuous trial R-charges, ∆ i ∈ (0, 2). Then F -maximization [21] implies that the twisted superpotential is extremized, as a function of the m i , when we set: where ∆ SC i are the superconformal R-charges. Now suppose that we also set: This implies that the second term in (5.63) vanishes. Moreover, the first term is already extremized as a function of the m i , by the argument above, and so applying the extremization procedure does not affect Z Mg,p (m i ). Thus, we conjecture that this gives a general solution to the truncation to minimal supergravity. In this case, we find the general relation: Using (5.64) and the fact that the S 3 free energy agrees with the free energy of the NUT solution, we may rewrite this as: Note the dependence on p and g is precisely the one we found in supergravity above. This can be thought of as a generalization of the "universal twist" of [27] to the case p = 0. Moreover, in order for the fluxes to be well-defined, we must impose: When ∆ SC i is rational, as in the cases above, this allows for families of solutions. Note the free energy, (5.68), and uplift condition, (5.69) agree with those found in the special cases JHEP05(2018)116 of the ABJM and V 5,2 models above. We conjecture these map to the M-theory uplift conditions for the corresponding supergravity backgrounds.
Finally, we note that there is a similar truncation corresponding to the Bolt − solution. Here we work in the background (ν R = − 1 2 , n R = − p 2 + g − 1), and find: (5.70) In cases where both truncations exist for the same g and p, such as for the ABJM theory discussed in section 4, these typically correspond to different choices of masses and fluxes in the field theory, leading to two different truncations to minimal supergravity. In this case, one of the solutions (namely, the Bolt + for p > 0 or the Bolt − for p < 0) is dominant over the other, and this will be picked out as the preferred solution when we perform the extremization of masses outlined at the beginning of this subsection.

Comparison with S 3 partition function
In this final section, we briefly discuss a curious relation between the extremal value of the twisted superpotential, W ext ([m i ]), and the large N limit of the S 3 partition function, as a function of trial R-charges ∆ i , which was first observed by [13], and which was used at the end of the previous section. Namely, they observed that: 25 where Z S 3 is the S 3 partition function, which was computed in the large N limit by [15], provided that we identify the trial R-charges ∆ i and the [m i ] via: In particular, recall one can find a solution to the extremization of W when one imposes the following conditions, for the superpotential constraints (5.1): which reproduces the constraint on the R-charges appearing in the S 3 partition function. Since the S 3 partition function appears as a special case of the M g,p partition function, namely S 3 = M g=0,p=1 , it is natural to ask if this relation can be understood in our framework. More precisely, it was argued in [18] that the M g=0,p=1 partition function in the "physical gauge", (5.27), reproduces the standard partition function on the round S 3 [20][21][22], namely:

JHEP05(2018)116
However, note that the identification between the mass parameters and trial R-charges here is not the one, (6.2), appearing in the relation (6.1). Moreover, we saw in section 4 that, in the case of the superconformal R-charges for the ABJM theory, ∆ i = 1 2 , we did not find a solution for the M g=0,p=1 partition function using our method, while the relation (6.1) above still holds in this case. Thus, we do not see a direct relation between our computation of the M g=0,p=1 partition function and the observation (6.1).
However, it turns out we can approach this relation from another point of view, through the lens space S 3 /Z 2 . The partition function on this space was first computed by Benini, Nishioka, and Yamazaki [67] (see also [42]), who considered more generally the case S 3 /Z p . In general, their background on S 3 /Z p is different from the one considered above and in [18], as the latter includes a flat R-symmetry gauge field, while the former does not. However, we claim that one actually recovers the partition function of BNY in the case p = 2 from the M g,p partition function if one takes the background of (5.34), i.e.: Then (5.26) implies we may take the r i ∈ R. To see this, we observe that the partition function of a chiral multiplet on the S 3 /Z 2 background of [42,67] can be written as: where s ∈ Z 2 is the torsion flux and r ∈ R is the R-charge, as one can straightforwardly show from the infinite product formula for the chiral multiplet given in [42,67]. This agrees with (5.28) up to a redefinition m → m/2. One may check this relation also holds for the Chern-Simons and gauge contributions, and so we find, for a general gauge theory: Recall the background (6.5) is precisely the one which allows us to take a large N limit of the M g,p partition function for a general quiver. Then, from (5.6) specialized to the case g = 0, p = 2, we find: Here we impose: Then a simple choice for the n i is to take them all to be zero, in which case we find the result: In fact, we can make a stronger statement. We claim that in the case n i = 0, and for arbitrary m i , the functionals computing −2W and 1 2πi log Z M g=0,p=2 for arbitrary eigenvalue JHEP05(2018)116 distributions are identical. For example, for the contribution of a bifundamental chiral multiplet, we have, from section 2.3: where recall g(u) = − 1 12 u(u − 1)(2u − 1), and in the present case: Then one finds the difference of the two functional gives: However, when we impose the constraint (2.39) and sum over all bifundamental chirals, we find that this vanishes. One can similarly check the other ingredients in the two functionals agree, and so we have: This immediately implies that their extremized values, (6.10), agree. Moreover, it implies that the two methods discussed in section 4.3 will give the same result on this space, since one is extremizing the same functional in both cases. This is indeed what happens for the supergravity solutions. For p = 2 the Bolt + has zero flux, and its on-shell action coincides with the NUT/Z 2 one. The two branches are actually continuously connected: they join at the point in phase space corresponding to s = 1/2 √ 2. For instance for ABJM we have: Given the proposed correspondence between the two methods and the free energy of the NUT and Bolt solutions, respectively, this suggests that these two supergravity solutions will agree for the case with conformal boundary S 3 /Z 2 , not just for minimal supergravity, but for general N = 2 gauged supergravity theories. It would be interesting to explore the special role played by S 3 /Z 2 in more detail. Finally, it was argued in [42] that the large N limit of the partition function on the BNY lens space, S 3 /Z p , is related to that of S 3 by: This relation arises because the bulk fillings for these boundary manifolds are simply Z p quotients of AdS 4 , and so their on-shell action is related by a factor of 1 p to the bulk dual of the S 3 partition function. Then if we combine (6.17) with (6.10), we arrive at the same relation (6.1), noted by [13].

JHEP05(2018)116 7 Discussion
In this paper we accomplished the task of computing the large N limit of the M g,p partition function for three-dimensional SCFTs with holographic M-theory duals and showing that it matches with the on-shell action of the minimal supergravity Bolt solutions. A few remarks are in order here.
First of all, in solving the supergravity equations of motion we have made some assumptions on the form of the solutions. For instance, for g = 0 we have imposed SU(2) × U(1) symmetry in the bulk, and we looked for configurations with a real metric. In principle we cannot exclude the existence of further M-theory solution with the same boundary data, which yield a different value for the on-shell action. It would be interesting to investigate further the phase space of solutions obtained when one releases one or more assumptions mentioned above.
One obvious extension of this work consists of incorporating vector multiplets in the gravity side, namely working with Bolt configurations which arise as solutions of N = 2 gauged supergravity with vector multiplets. Solutions of the U(1) Fayet-Iliopoulos STU model can be lifted up to 11 dimensions [68][69][70], and the analysis of Euclidean NUT or Bolt solutions in this framework is currently work in progress. Notice that Lorentzian solutions with non-trivial scalar profiles and NUT charge in theories of FI-gauged supergravity were found in [71,72].
With the addition of vectors and scalars, we expect a more intricate phase space to arise for the NUT and Bolt configurations. Phase transitions involving different solutions can arise 26 and this should be reflected in the field theory computation. It would also be interesting to study the quantum corrections to the on-shell action, and match them to the finite N result in the field theory, as was done recently in the context of AdS 4 black holes in [75].
As discussed in section 4, it would be interesting to study the large N limit of the integral formula (2.16) for the M g,p partition function, and compare this to the large N limit we found using the sum formula (2.17). In particular, we may hope to compare these different large N solutions to different (leading or sub-leading) saddles in the supergravity partition function, such as the NUT and Bolt solutions above.
Lastly, one can consider more general quiver gauge theories, including those with massive type IIA supergravity duals [76], as considered in [10,11,13]. It would also be desirable to obtain the partition function for the quiver theory dual to the M-theory reduction on M 32 . As mentioned at the end of section 5.3, for this theory all Bolt solutions are allowed, including in particular the p = 1 Bolt + , which is a regular filling for the squashed S 3 for a finite range of squashing parameters. It would be interesting to see if the field theory computation is able to reproduce its on-shell action I Bolt + p=1 = 3 4 I NUT , which is lower than the S 3 free energy computed in [65,66]. The question of finding the correct vacuum of the theory is of course entangled with the fact that, as already mentioned, we cannot rule out the presence of (perhaps less symmetric) branches of solutions with the same boundary data. A deeper understanding of these points would be desirable.
All these questions are left for future investigation and we hope to report back in the near future.

JHEP05(2018)116
• for R 2 (or, upon compactification g = 1) we choose τ 1 = sin ψdθ + cos ψdφ τ 2 = − cos ψdθ + sin ψdφ τ 3 = dψ − θdφ (A.8) • for H 2 (or, upon compactification g > 1) we have τ 1 = sin ψdθ + sinh θ cos ψdφ τ 2 = − cos ψdθ + sin ψ sinh θdφ τ 3 = dψ − cosh θdφ (A.9) We take the following basis of four-dimensional gamma matrices, which is the same used in [25]: where σ i , i = 1, 2, 3 are the Pauli matrices. The supersymmetry parameter is a Dirac spinor. It is convenient to write it in terms of its positive and negative chirality parts In what follows we moreover further distinguish the single parts of + and − as , so that the 4d Dirac spinor has this form Radial component of the Killing spinor equation. We write the warp factor λ(r) as where r α , α = 1, 2, 3, 4, are the four roots of the warp factor, and their explicit form is in (3.40). If we decompose the r component of the Killing spinor equation (A.1) into chiral parts, we obtain and In analogy to [25], we impose the following projections on the Killing spinor:

JHEP05(2018)116
Using the projection relation (A. 16) we find a solution to equations (A.14) and (A.15) of the form where υ is a two-component spinor which does not depend on r, and the two radial functions F (r) and G(r) read: 3d Killing spinor equation. Once we have inserted (A.17) the remaining components t, θ, φ of the Killing spinor equation (A.1) reduce to the following three-dimensional equation the spinor υ: Here ∇ (3) i and γ i are respectively the covariant derivative and gamma matrices for the 3d boundary metric ds 2 3 = τ 2 1 + τ 2 2 + 4s 2 τ 2 3 with spin connection where f (θ) is defined in (3.8). The 3d gauge field appearing in (A. 19) is the asymptotic value of the gauge field A, hence it has the form Choosing the vielbeinsē 1 = τ 1 ,ē 2 = τ 2 ,ē 3 = 2sτ 3 it is easy to see that the components of equation (A.19) are solved by where υ (0) is a constant. The full Killing spinor then reads The solution therefore preserves two supersymmetries out of eight, it is then 1/4 BPS. Notice that the Killing spinor only has radial dependence, hence no further supersymmetries are broken while taking the quotient to obtain a higher genus surface. Let us mention that if we had chosen the static vielbein, which for instance for κ = 1 reads

I(p)
Bolt+ , NUT+flux (g=0) Bolt+ NUT+flux+ Figure 3. This plot is showing the behavior of the on-shell action (in units of π/G 4 ) for the Bolt + solution and the NUT/Z p with the same units of flux for 1 < p < 12, as a function of p. One can see that the Bolt + has always lower free energy, so dominates the ensemble.
In the large N limit, t a becomes a continuous variable t, and we define: normalized by dtρ(t) = 1, and v α a → v α (t).

C.1 Contributions to twisted superpotential
We start with the contributions to the twisted superpotential. These were derived already in [9,13], however we review them here for completeness, and to establish the relation of our notations. We also employ a slightly different argument in some places than used in these works. Our starting point is the expression for the twisted superpotential at finite N in (2.20). Let us review the various ingredients in the U(N ) quiver gauge theories. For the CS terms, we have: If we take the ansatz (C.1), this becomes: Imposing the constraint α k α = 0, we may write this to leading order in N as: Next we compute the contribution of bifundamental and adjoint fields. Here we will take a slightly different approach from [7], and instead use the following argument, modified JHEP05(2018)116 from that of [15]. We will consider the large N limit of functions of the form: where u a andû b are the eigenvalues for two of the U(N ) gauge groups (or possibly the same one, in the case of adjoint fields), and f is some function, which may depend on other parameters which we suppress. Let us rewrite this as: We expect the second term to be subleading in N , but we will return to its contribution below. Focusing on the first term for now, and plugging in the ansatz (C.1), this can be written as: Now, expecting most of the contribution to come from the region near t = t , we make a change of variables t → ξ = N 1/2 (t − t ). Then this becomes, to leading order: where we defined δv(t) = v(t) −v(t), and we take the principle value to exclude the contribution from t = t . In the present case, a single bifundamental chiral multiplet contributes: W π (u a −û b + m) , (C.10) where we have defined the contribution of a chiral multiplet in the parity-preserving regularization (i.e., without the level − 1 2 CS term) as: 29 W π (u) = 1 (2πi) 2 Li 2 (e 2πiu ) + 1 4 [u]([u] − 1) + 1 24 , (C. 11) where, as in the main text, we define, for complex u: [u] = u − n, n ∈ Z, such that 0 < Re([u]) ≤ 1 . (C.12) Here we have used the freedom to change branch, taking W → W + nu + m, n, m ∈ Z, to make the twisted superpotential periodic under u → u + 1, which will be convenient below. 29 Here and below, we work on the principle branch of the polylogarithms, defined so that Lis(e 2πiu ) → 0 as u → i∞, which fixes the branch in the upper half u-plane, and extending continuously along vertical lines to the lower half plane. Then the resulting functions are periodic under u → u + 1, with branch cuts along Re(u) ∈ Z for Im(u) < 0.
To deal with this divergence, we make two observations. First, we only expect that certain choices of matter content and superpotential terms lead to theories with appropriate supergravity duals, with N 3/2 scaling of the degrees of freedom. Specifically, as in section 2.3.1, we impose: • The number of incoming and outgoing edges (bifundamental chirals) at any node (U(N ) gauge factor) in the quiver are equal.
• The sum of all masses of bifundamentals charged under a node is zero. 30 Then when we sum over the contributions of all chiral multiplets, we find the coefficient of ξ 2 max is: where I runs over the bifundamentals, with the Ith fundamental connecting the α I th and β I th gauge groups, and n I = v α I − v β I + m I − [v α I − v β I + m I ] ∈ Z, and we have used the 30 More precisely, it need only be integer for the argument below to work, and in the context of the twisted superpotential the integer part of the masses is irrelevant, so imposing the sum is precisely zero is unnecessary. However, when computing the partition function below, the integer part will be important, and we will need to impose this condition.

JHEP05(2018)116
conditions above to cancel the dependence on the v α and m I . Then the r.h.s. , while not necessarily zero, is 1 2 times an integer. Then there is one more fact we must use, which is that the twisted superpotential is only defined modulo changes of branch. For finite N , such changes of branch take the form: W → W + n a u a + n i m i + n, n a , n i , n ∈ Z .
(C. 20) For infinite N , we choose to focus on changes of branch which preserve Weyl symmetry of the gauge group, and consider only those terms which appear at leading order in N . Then there are two basic choices, which lead to the following terms when we plug in the ansatz (C. where we must take α n α = αn α so that this term is of order N 3/2 . Then we note that, after summing over the contributions from all chirals, the divergent terms from (C.17) are precisely of the form of the second change of branch above, so we may remove them. Then we find the contribution of a bifundamental chiral multiplet is: Note that while this is cubic in δv, the first consistency condition above implies that the cubic terms for each v α cancels when we sum over all chirals, and so the final functional is always a quadratic function of the v α . Finally, let us return to the subleading contribution in (C.7): N a=1 f (u a −û a ) → N dtρ(t)f (δv(t)) .

(C.24)
For the case of a bifundamental chiral multiplet, this gives: N dtρ(t)W π (δv+m) = N dtρ(t) 1 (2πi) 2 Li 2 (e 2πi(δv+m) )+ (C.25) 31 Note the condition that the twisted superpotential be finite fixes the choice of change of branch of the second type above, while the first type changes the twisted superpotential by a finite piece, and so is not specified uniquely.

JHEP05(2018)116
Naively this is subleading in N , however, note that as δv + m approaches an integer, the derivative of this expression with respect to δv diverges. This means the subleading term will start to compete with the leading term, and we will have to take it into account below. Specifically, let us write, for n ∈ Z: δv(t) + m = n + Ce −2πN 1/2 Y I (t) .
Using the limit: Next we consider the contributions to Z Mg,p . Here the starting point is the finite N expression for the partition function as a sum over Bethe vacua, (2.26). Specifically, we consider the contribution from a given Bethe vacuum, and take the eigenvalues u a in this vacuum to be given by the distribution ρ(t), v α (t), as above. Then we compute the functional:

JHEP05(2018)116
where the sum is taken over all bifundamental chirals charged under the αth gauge group, and this must hold for all gauge groups, and adjoints are counted twice. These two conditions can be summarized by: I∈α I = 0 , (C. 42) where in the sum we include the contribution of the vector multiplet, with I = g − 1.
Then, using also the condition I∈α m I = 0 imposed above, we find the coefficient of the term quadratic in ξ 2 max vanishes. Thus, we are left with the following finite contribution from a bifundamental: log Z Mg,p,bif = −2πN 3/2 dtρ(t) 2 G (δv + m) .

JHEP05(2018)116
Finally, we also have the contribution from the Hessian: (C. 51) Then, naively this has order N log N if all the components of the matrix stay finite, but near the tail regions there are divergences which cause it to contribute at leading order. As argued in [7], one finds: log H ≈ 2πN 3/2 tails t|δv(t)+m≈n dtρ(t)Y (t) . (C.52) Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.