Supercurrent anomalies in 4d SCFTs

We use holographic renormalization of minimal $\mathcal{N}=2$ gauged supergravity in order to derive the general form of the quantum Ward identities for 3d $\mathcal{N}=2$ and 4d $\mathcal{N}=1$ superconformal theories on general curved backgrounds, including an arbitrary fermionic source for the supercurrent. The Ward identities for 4d $\mathcal{N}=1$ theories contain both bosonic and fermionic global anomalies, which we determine explicitly up to quadratic order in the supercurrent source. The Ward identities we derive apply to any superconformal theory, independently of whether it admits a holographic dual, except for the specific values of the $a$ and $c$ anomaly coefficients, which are equal due to our starting point of a two-derivative bulk supergravity theory. We show that the fermionic anomalies lead to an anomalous transformation of the supercurrent under rigid supersymmetry on backgrounds admitting Killing spinors, even if all superconformal anomalies are numerically zero on such backgrounds. The anomalous transformation of the supercurrent under rigid supersymmetry leads to an obstruction to the $Q$-exactness of the stress tensor in supersymmetric vacua, and may have implications for the applicability of localization techniques. We use this obstruction to the $Q$-exactness of the stress tensor, together with the Ward identities, in order to determine the general form of the stress tensor and $R$-current one-point functions in supersymmetric vacua, which allows us to obtain general expressions for the supersymmetric Casimir charges and partition function.


Introduction
Supersymmetric quantum field theories on curved backgrounds admitting a notion of rigid supersymmetry have attracted considerable interest in recent years, mainly due to localization techniques. Such techniques utilize the Q-exactness of certain operators on compact manifolds admitting Killing spinors, and can be used to compute various observables exactly, for any value of the couplings [1].
All these approaches lead to consistent theories with rigid supersymmetry, and determine the classical Lagrangian, as well as the supersymmetry transformations of the classical fields. Moreover, one can apply the Noether procedure to this Lagrangian in order to derive the classical Ward identities that local gauge-invariant operators satisfy, reflecting the global symmetries of the theory. However, to obtain the quantum Ward identities one needs to compute the appropriate path integral on a general curved background in order to determine the global anomalies. Global anomalies for N = 1 supersymmetric theories in d = 4 have been classified using superspace cohomology arguments in [16], and have been computed explicitly in a number of examples, both in superspace [17,18], and in component language [19].
However, both the rigid supersymmetry transformations on a curved background and the general form of the Ward identities, including the quantum anomalies, can also be determined through holographic techniques, as was done for gauge transformations and the axial anomaly in [20], or for Weyl transformations and the trace anomaly in [21]. In particular, starting with a suitable gauged supergravity in asymptotically AdS d+1 space, one can derive the action of rigid supersymmetry on the sources of local gauge-invariant operators on the d-dimensional conformal boundary, as well as the quantum Ward identities these operators satisfy. The crucial point here is that the form of the Ward identities and of the global anomalies is determined by symmetries and is therefore universal. In particular, the same form of the Ward identities applies to a wider class of supersymmetric theories, which need not admit a holographic dual.
Of course, the holographic calculation does not determine the elementary Lagrangian on any curved background, except possibly for terms that are protected by non-renormalization theorems. Moreover, even when specifying the theory in terms of gauge-invariant operators and their sources, there are two possible caveats in deriving the rigid supersymmetry transformations and the quantum Ward identities using holographic techniques. The first is that the version of rigid supersymmetry induced holographically on the conformal boundary by bulk (minimal) gauged supergravity corresponds specifically to that obtained from conformal supergravity [22] on the boundary through the Festuccia-Seiberg argument [4,11] 2 . However, this reflects the fact that the holographic Ward identities are expressed in terms of the operators at the ultraviolet fixed point, and apply to both conformal and massive theories. The second is that the holographic calculation leads to specific values for the anomaly coefficients, which generically do not apply to theories without a holographic dual, or even to holographic theories at weak coupling. For instance, starting from two-derivative supergravity in AdS 5 necessarily leads to a = c in the resulting Ward identities on the conformal boundary [21]. However, starting from a higher derivative supergravity theory in the bulk may lift this degeneracy.
In this paper we focus on N = 2 superconformal theories in d = 3 and N = 1 superconformal theories in d = 4. Starting with minimal N = 2 gauged supergravity in four and five dimensions respectively, we use holography to derive the general form of the quantum superconformal Ward identities on an arbitrary curved background, including an arbitrary fermionic source for the supercurrent. A complementary analysis for N = 2 gauged supergravity coupled to hypermultiplets is carried out in [23]. The minimal N = 2 gauged supergravity describes holographically the current multiplet, consisting of the stress tensor, T i a , the R-symmetry current, J i , and the supercurrent S i . Their respective local sources are the vielbein e a i (0) , a U (1) gauge field A (0)i , and a chiral gravitino Ψ (0)+i , whose values specify the field theory background. Notice that in the presence of fermionic sources it is necessary to introduce the vielbein as a fundamental source, rather than the metric g (0)ij . In the absence of fermion sources, however, one may work exclusively with the metric g (0)ij and its conjugate symmetric stress tensor T ij . The superconformal Ward identities we obtain holographically are given in eq. (5.4). For N = 1 superconformal theories in d = 4 they are 3 where the local functions of the sources A M i , A R , A S , A W and A sW are given in eq. (5.3) and are related to the global anomalies of the theory. We have slightly simplified the notation here for the sake of brevity, but we refer to (5.4) and appendix A for the precise form of the Ward identities and a guide to our notation, respectively. For N = 2 superconformal theories in d = 3 the Ward identities are of the same form, except that the local terms A M i , A R , A S , A W and A sW are absent. For a generic N = 1 superconformal theory in four dimensions, the Weyl and R-current anomalies take the form where are respectively the Euler density and the supersymmetrized Weyl squared conformal invariant, while P [g (0) ] is the topological Pontryagin density. In order to specify the numerical factors, besides the a and c anomaly coefficients, one would need to appropriately normalize the gauge field A (0)i . The precise factors for the correctly normalized gauge field can be found e.g. in [25], which also corrected a sign in the original expressions for the anomalies in [19]. As anticipated, our computation reproduces these anomalies but with a = c. In particular, we obtain the full superconformal invariant W 2 , up to quadratic order in the supercurrent source. The terms A M i and A S that appear respectively in the diffeomorphism and supersymmetry Ward identities are related to the R-symmetry anomaly A R and take the form The appearance of these terms in the Ward identities simply reflects the fact that the R-current operator J i is the consistent current (see e.g. [26] for a recent overview). Writing the Ward identities in terms of the covariant (and gauge invariant) current J i cov , these terms get eliminated. In section 5 we confirm that the generating functional is invariant under local diffeomorphisms 4 and so indeed A M i does not represent an anomaly in the diffeomorphism invariance of the theory. In contrast, the term A S does lead to a non-invariance of the generating functional under local supersymmetry transformations, which will play a crucial role throughout our analysis. Finally, the terms A W and A sW represent anomalies under local Weyl and local superWeyl transformations, respectively.
Our main result concerns the consequences of the two fermionic anomalies, namely the supersymmetry anomaly A S and the superWeyl anomaly A sW , on the transformation of the supercurrent under rigid supersymmetry. In particular, we show that on backgrounds admitting a (conformal) Killing spinor ζ + , which satisfies the Killing spinor equation (6.9), the anticommutator of the corresponding supercharge with the supercurrent is generically anomalous, namely The supercurrent anomaly A i ζ , given explicitly in (6.20), is related to the supersymmetry and superWeyl anomalies as Moreover, we show through explicit examples that there are backgrounds admitting Killing spinors where all anomalies in the superconformal Ward identities are numerically zero, yet the supercurrent anomaly A i ζ is non-zero. In those cases, the supercurrent anomaly poses an obstruction to the Qexactness of the operator Although in two dimensions the anomalous term in the transformation of the supercurrent survives in flat space, leading to a central extension of the superVirasoro algebra, in four dimensions the anomalous term is manifest only in curved space. However, even in four dimensions, the consequences of this anomalous transformation of the supercurrent should be visible in flat space correlation functions. We show that the anomalous transformation (1.5) of the supercurrent under rigid supersymmetry provides a resolution to two paradoxes discussed recently in [27]. Firstly, considering the class 4 I am grateful to Davide Cassani for pointing this out to me. of supersymmetric backgrounds in eq. (6.1), which were originally obtained in [4,5], we revisit the argument of [28,29,30], according to which the supersymmetric partition function on this class of backgrounds is independent of the complex functions u(z,z) and w(z,z). We show that the u and w variations of the supersymmetric partition function are indeed a linear combination of certain components of the the operator (1.7). However, the fact that the supercurrent anomaly is generically non-zero on these backgrounds means that this operator is not Q-exact, and hence the supersymmetric partition function is not invariant under deformations of the functions u(z,z) and w(z,z). Although our derivation of the superconformal Ward identities is holographic, the argument we use to show the non-invariance of the supersymmetric partition function is a field theory argument, using only the supercurrent anomaly. However, the result matches precisely the answer obtained in [27], by explicitly evaluating the bulk on-shell action. We therefore conclude, that holographic renormalization is perfectly compatible with supersymmetry, and it is in fact the field theory assumption that the supersymmetric partition function is invariant that is not quite correct!
The second puzzle concerns the BPS relation among the conserved charges in supersymmetric states and the supersymmetric Casimir energy, which has been discussed extensively in the recent literature [31,32,33,34,35,27]. The BPS relation follows again from the Q-exactness of the operator (1.7), but the presence of a non-zero supercurrent anomaly implies that the BPS relation is also anomalous! We show that the supercurrent anomaly together with the Ward identities provide enough constraints to determine the one-point functions of the stress tensor and the Rcurrent in any supersymmetric state in terms two arbitrary scalar functions. This general solution for the supersymmetric one-point functions allows us to obtain general formulas for the conserved charges for any BPS state, including the supersymmetric Casimir charges of the global vacuum, as well as the supersymmetric partition function. In fact, we show that there is a one-parameter family of consistently defined charges, all of which satisfy the (anomalous) BPS relation. Moreover, the supersymmetric partition function is independent of the definition of the charges.
The rest of the paper is organized as follows. In section 2 we briefly review the minimal N = 2 gauged supergravity in four and five dimensions. Section 3 sets the stage for the subsequent holographic analysis and discusses the radial Hamiltonian formulation of the bulk supergravity dynamics. Section 4 contains our main technical results, where we carry out the procedure of holographic renormalization in order to derive all boundary counterterms, both bosonic and fermionic. Moreover, we show that there is a two-parameter family of supersymmetric renormalization schemes that also preserve parity, corresponding to adding the two possible superconformal invariants (4.19). In section 5 we use these boundary counterterms in order to define holographically the local current operators and derive the quantum superconformal Ward identities, including all quantum anomalies. In section 6 we discuss the anomalous transformation of the supercurrent under rigid supersymmetry, and we revisit the dependence of the supersymmetric partition function on the functions u(z,z) and w(z,z), parameterizing the supersymmetric backgrounds found in [4,5]. Finally, we use the anomalous transformation of the supercurrent, together with the Ward identities, in order to obtain general expressions for the one-point functions of the stress tensor and the R-current in any supersymmetric state, and to derive general expressions for the Casmir charges and the supersymmetric partition function. Our notation and a number of technical details are presented in three appendices.

Minimal N = gauged supergravity
Our focus in this paper is on the minimal gauged supergravities in four and five dimensions [36,37], which we briefly revisit in this section, keeping terms up to quadratic order in the gravitino. The field content of minimal gauged supergravity in both four and five dimensions comprises of the gravity multiplet only, consisting of the metric g µν , a complex Dirac gravitino Ψ µ (equivalently a pair of symplectic-Majorana gravitini) and the graviphoton A µ . 5 The minimal gauged supergravity action is obtained by gauging the U R (1) subgroup of the U Sp R (2) ∼ = SU R (2) R-symmetry group of N = 2 Poincaré supersymmetry and takes the form [36,37] where κ 2 = 8πG d+1 is the gravitational constant and are respectively the cosmological constant and U R (1) gauge coupling in D = d + 1 dimensions. As we verify in appendix C, the constants are fixed by supersymmetry. Notice that c 1 is zero unless d = 4, i.e. D = 5. It should be emphasized that although we have written the action (2.1) for generic dimension D = d + 1 in order to treat the cases D = 4 and D = 5 simultaneously, this action is not supersymmetric in dimensions other than four or five. The covariant derivative acts on the gravitino as where ω µαβ is the spin connection, Γ ρ µν is the Christoffel symbol of the metric g µν , and Γ αβ denotes the antisymmetrized product of two gamma matrices, as defined in appendix A, where we explain our index notation and gamma matrix conventions. 6 The bulk action (2.1) must be complemented by the Gibbons-Hawking boundary terms [38] which are required for the Dirichlet problem of the action (2.1) to be well posed on a finite cutoff surface ∂M ε . However, in order for the variational problem to be well posed on a non-compact asymptotically locally AdS space, additional boundary terms are required [39]. We will derive these 5 It can be checked that the number of on-shell bosonic and fermionic degrees of freedom match in both four and five dimensions, respectively as 2 + 2 = 2 + 2 and 5 + 3 = 4 + 4. 6 We hope that the distinction between the Christoffel symbols Γ ρ µν and the antisymmetrized products of gamma matrices will be clear from the context in the following. systematically in section 4. Moreover, the Gibbons-Hawking terms (2.5) are required for the action to preserve supersymmetry on the boundary.
The supersymmetry transformations that leave the action (2.1) off-shell invariant up to boundary terms are [36,37] where , c 4 = 1 and E α µ is the vielbein of the metric g µν , i.e. g µν = η αβ E α µ E β ν . These transformations are correct only to linear order in the gravitino, corresponding to the fact that the action (2.1) is specified to quadratic order in the gravitino. In appendix C we check explicitly that the supersymmetry transformations (2.6) leave the action (2.1) invariant up to boundary terms, which we compute.

Radial Hamiltonian formalism
The first step in order to construct the holographic dictionary for the supergravity theory (2.1) is to formulate the dynamics in Hamiltonian language, with the radial coordinate emanating from the conformal boundary playing the role of Hamiltonian time. This formulation allows one to systematically construct the covariant local boundary counterterms required in order to render the variational problem well posed on the conformal boundary, and to derive the renormalized observables in the dual field theory, as well as the Ward identities they satisfy, including all related global anomalies.
The radial Hamiltonian formulation of the dynamics singles out a radial coordinate r emanating from the conformal boundary ∂M in M and describes the evolution of the induced fields on the constant radius slices Σ r ∼ = ∂M under radial translations. This foliation of M need only hold in an open neighborhood of ∂M, since the variational problem and the holographic dictionary require knowledge of the space of asymptotic solutions only. Equivalently, the holographic dictionary and the boundary counterterms require only ultraviolet data.
The induced fields on the radial slices Σ r are obtained through a standard ADM decomposition of the bulk fields, given in eq. (A.7) in appendix A. This decomposition is somewhat more involved in the presence of fermion fields, since a Hamiltonian description of fermions requires an additional decomposition of the spinors according to radiality [40], defined in eq. (A.18). Further details on the ADM decomposition of spinor fields are provided in appendix A.
Radial Lagrangian Inserting the ADM decomposition (A.7) in the supergravity action (2.1) leads to a number of bulk terms that involve only up to first order radial derivatives on the induced fields, as well as a number of boundary terms that are exactly canceled by the Gibbons-Hawking terms (2.5). 7 We can therefore write S + S GH = Σr dr L, (3.1) where the radial Lagrangian is given by In deriving this Lagrangian we have decomposed the bulk vielbein as in eq. (A.9) and have partially fixed the SO(1, d) frame rotations by choosing the frame specified in eq. (A.13). As we shall see shortly, this does not lead to any loss of generality since one can obtain the first class constraint associated with frame rotations by invoking the fact that there is no torsion.
Canonical momenta As mentioned above, in order to define the symplectic variables associated with the gravitino we need to use the radiality projectors (A.18) to decompose its transverse components as [41] The variables Ψ +i and Ψ −i are symplectic conjugates and so one radiality can be treated as a generalized coordinate, while the other as the corresponding canonical momentum. We will adopt the convention that Ψ +i is the generalized coordinate and Ψ −i is (proportional to) the corresponding canonical momentum.
The canonical momenta conjugate to the vielbein e a i on Σ r , the gauge field A i , and the induced fermion fields Ψ +i and Ψ +i can be read off the radial Lagrangian (3.2) and are given by where we have defined for convenience the fermion bilinears Notice that the Lagrangian L does not contain any radial derivatives of the variables N , N i , a, Ψ r and Ψ r , and hence their conjugate momenta vanish identically. As we shall see momentarily, these non-dynamical Lagrange multipliers lead to a set of first class constraints, each associated with a local symmetry of the action (2.1). In particular, the lapse and shift functions N and N i reflect the diffeomorphism invariance of (2.1), a corresponds to local U (1) gauge transformations, while Ψ r and Ψ r reflect the local supersymmetry invariance. Notice, however, that there is no Lagrange multiplier associated with local frame rotations, although the action (2.1) is clearly invariant under such transformations. The reason for this is that the spin connection here is not an independent field and so there is no torsion. In a first order formalism, where the spin connection is treated as an independent field, the Lagrange multiplier associated with local frame rotations would be the spin connection component ω rab , while the phase space of the theory would include the generalized coordinate ω iab and its conjugate momentum. However, as can be seen in eq. (A.29), here ω rab is completely determined in terms of the vielbein e a i on Σ r , as well as the lapse and shift functions, and hence it is not an independent Lagrange multiplier. Moreover, ω iab and its conjugate momenta are not independent phase space variables.
Nevertheless, the fact that the torsion vanishes implies that the extrinsic curvature K ij , which can be identified with a certain component of the bulk Christoffel symbol (see eq. (A.26)), is symmetric, i.e. K [ij] = 0. It follows that canonical momenta (3.4) satisfy the algebraic constraint 6) or equivalently, and hence, on the reduced phase space where the fermion variables are set to zero the canonical variables e i a and π i a can be replaced by the induced metric γ ij and its conjugate momentum π ij . As we shall see later, this ensures that on a background where all the fermion sources are set to zero, the dual field theory possesses a symmetric and relativistic stress tensor.

Radial Hamiltonian
The radial Hamiltonian H corresponds to the Legendre transform of the Lagrangian (3.2) with respect to all dynamical variables, namely Introducing the gauge invariant momentum 9) and the gauge-covariant derivative (see also (A.32)) we find that the Hamiltonian takes the form , and we have defined F ± ≡ Γ ± F so that F = F + + F − . Since the momenta conjugate to the variables N , N i , a and Ψ r vanish identically, the form (3.11) of the Hamiltonian leads, via Hamilton's equations, to the first class constraints (3.13) We will see in the following that these constraints are directly related to the Ward identities of the dual field theory. An important observation is that the Hamiltonian constraint H = 0, as well as the supersymmetry constraints F ± = 0, are quadratic in the canonical momenta, while all other constraints are linear. This reflects the fact that the symmetries associated with the Hamiltonian and supersymmetry constraints are spontaneously broken by the radial slice Σ r and hence realized non linearly, while those associated with the remaining first class constraints are preserved. The non-linear constraints corresponding to the symmetries broken by the cutoff result in global anomalies in the dual quantum field theory. This is completely analogous to the way quantum anomalies appear in standard quantum field theory, through symmetry breaking by the ultraviolet regulator.
Hamilton-Jacobi formalism and flow equations The Hamiltonian determines the radial evolution of the induced fields e a i , Ψ +i and A i through Hamilton's equationṡ (3.14) In the gauge (A.8) only the Hamiltonian constraint contributes to these equations, which reaḋ Moreover, inverting the expression for π i Ψ in (3.4) we find that Ψ −i is given by These equations allow us to determine the radial evolution of the induced fields without invoking the second order field equations, or equivalently the other half of the Hamilton equations, by using instead the Hamilton-Jacobi expressions where S[e, A, Ψ + ] is Hamilton's principal function. In particular, inserting the expressions (3.17) for the canonical momenta in the first class constraints (3.13) leads to a set of functional differential equations for the functional S[e, A, Ψ + ]. Given a solution S[e, A, Ψ + ] of the these Hamilton-Jacobi equations, the flow equations (3.15) become first order equations for the induced fields. In the next section we will obtain the general asymptotic solution of the Hamilton-Jacobi equations.
In combination with the relations (3.17) and the flow equations (3.15) this solution allows us to determine the general asymptotic Fefferman-Graham expansions for the induced fields e a i , Ψ +i and A i , without solving the second order field equations. Since Hamilton's principal functional S[e, A, Ψ + ] coincides with the on-shell action evaluated with a radial cut-off r o , the asymptotic solution of the Hamilton-Jacobi equations also determines the divergences of the on-shell action, and hence the covariant boundary counterterms. Hence, the complete holographic dictionary can be obtained from a single object, namely the asymptotic solution of the Hamilton-Jacobi equations.

Recursive solution of the Hamilton-Jacobi equations
In this section we determine the asymptotic solution of the Hamilton-Jacobi equations, corresponding to the constraints (3.13), with the canonical momenta expressed as gradients of Hamilton's principal function S[e, A, Ψ + ]. As we will see, our algorithm requires solving only the Hamiltonian constraint H = 0, with all other constraints being automatically satisfied, up to the corresponding order in the asymptotic expansion. The resulting asymptotic solution for Hamilton's principal function will allow us to obtain the asymptotic Fefferman-Graham expansions for the bulk fields, and will be the basis for constructing the holographic dictionary in the next section.

Covariant expansion and recursive algorithm
Since the conformal dimensions of the operators dual to the fields in the Lagrangian (2.1) are known and fixed, we can obtain the general asymptotic solution of the Hamilton-Jacobi equation H = 0 by using the dilatation operator method [42]. From the leading asymptotic form of the fields in eq. (B.11) follows that the dilatation operator takes the form where S (n) = d d x L (n) are eigenfunctions of the dilatation operator, i.e. δ D S (n) = (d − n)S (n) . By construction, S (n) are covariant functionals of the induced fields e a i , A i and Ψ +i and higher order eigenfunctions, i.e. larger n, are asymptotically subleading relative to lower order ones. The formal expansion (4.2) is therefore a covariant asymptotic expansion. In order to utilize the expansion (4.2) we observe that the Hamilton-Jacobi relations (3.17) imply that for a generic variation of the induced fields [43] where ∂ i v i indicates a generic total derivative term. Applying this identity to local scaling transformations generated by the dilatation operator (4.1) gives where we have used the fact that L (n) are only defined up to a total derivative and hence can be absorbed in L (n) . The final ingredient necessary to set up the recursive algorithm for determining the terms S (n) in the expansion (4.2) is the leading asymptotic behavior of the fields in eq. (B.11). Together with the flow equations (3.15)-(3.16) and the expressions (3.17) for the canonical momenta, this leading asymptotic form of the fields determines the zero order solution S (0) to be The corresponding canonical momenta are 8 Inserting the expansion (4.2) in (3.12a) and using the identity (4.4) and the leading solution (4.5) of the Hamilton-Jacobi equation, the Hamiltonian constraint H = 0 reduces to a tower of linear equations for L (n) , n > 0, namely where the inhomogeneous terms up to order n = 4 are given by In writing these expressions we have again ignored terms higher than quadratic in the gravitino. From (4.7) and (4.8) follows that the first subleading term in the covariant expansion (4.2) is Varying this expression with respect to each of the induced fields gives the canonical momenta where we have used the identities (see e.g. (7.96) in [44]) and 12) in order to obtain the expression To proceed to the next order we insert the canonical momenta (4.10) in the expression for R (4) in (4.8). After some algebra we obtain and from (4.7) follows that If the boundary dimension d is not specified, the above procedure can be repeated indefinitely in order to determine the terms L (n) to arbitrary order. However, the recursion procedure stops at the term of zero dilatation weight. In particular, we only need to determine the terms L (n) for n ≤ d. Since we are interested in the cases d = 3 and d = 4 here, the above results suffice. However, as is clear from the above expressions, there is a significant difference between the two cases, due to the fact that for d = 4 one needs to go up to the term L (4) , which contains a pole when d = 4. This pole indicates the presence of a conformal anomaly and should be removed through the dimensional regularization prescription [42,43] 1 where r o is the radial cutoff in the Fefferman-Graham coordinates (A.7)-(A.8). More specifically, one replaces L (4) in the expansion (4.2) according to the rule The fact that the resulting asymptotic solution of the Hamilton-Jacobi equation depends explicitly on the radial cutoff is the holographic manifestation of the conformal anomaly [21] and implies that the implicitly covariant and local expansion in (4.2) is not a fully consistent asymptotic solution in this case. Alternatively, one could have started instead with a formal expansion that contains the logarithmically divergent term L (4) , in which case the recursion procedure would not produce any poles at d = 4. The two approaches are equivalent and an explicit comparison can be found in the recent review [45]. The punchline of the above analysis is that the general asymptotic solutions of the Hamilton-Jacobi equation for d = 3 and d = 4 take the form where S (d) is undetermined but finite (has dilatation weight zero), while S (0) , S (2) and S (4) = − 2 d 4 x R (4) correspond to the divergent part of the solution S and are given in equations (4.5), (4.9) and (4.14), respectively. The ellipses in (4.18) stand for terms with negative dilation weight that asymptotically go to zero, and hence are not relevant for the subsequent analysis. As we show next, this asymptotic solution of the Hamilton-Jacobi equation can be used to determine the Fefferman-Graham asymptotic expansions of bulk fields, and to derive the holographic dictionary.

Finite local counterterms and the supersymmetric renormalization scheme
The first implication of the asymptotic solution (4.18) of the Hamilton-Jacobi equation is that it determines the local and covariant boundary counterterms necessary to render the variational problem on the conformal boundary well posed, and consequently to renormalize the on-shell action. Hamilton-Jacobi theory identifies Hamilton's principal function S with the on-shell action, and so the divergences of the on-shell action are in one-to-one correspondence with the divergences of the asymptotic solution (4.18) of the Hamilton-Jacobi equation [46]. Hence, the counterterms can be defined as minus the divergent parts of the solution (4.18). The resulting renormalized action is invariant under boundary U(1) gauge transformations, diffeomorphisms, local Weyl transformations, as well as supersymmetry and superWeyl transformations, up to specific anomalies in the case d = 4. We will discuss these symmetries and derive the corresponding anomalies in section 5. In this subsection, however, we would like to specify the most general form of the boundary counterterms compatible with the above symmetries, including supersymmetry.
In the case d = 4, besides the divergent part of the solution (4.18), one may include in the boundary counterterms certain finite, local, and covariant terms that do not cancel any divergences, but may contribute to the value of the renormalized on-shell action, as well as to the value of certain one-point functions derived thereof. Such terms must preserve all local symmetries on the cutoff, since otherwise they would give rise to spurious anomalies, i.e. anomalies that are trivial cocycles and can be removed by the very same finite boundary terms. An example of a finite local boundary term that explicitly breaks local Weyl invariance (as well as supersymmetry) is the R 2 term, which leads to a total derivative contribution to the conformal anomaly. By construction, this contribution to the conformal anomaly is a trivial cocycle, that can be eliminated by removing the R 2 term from the action. Since we are interested in minimizing the possible contributions to the anomalies, in the following we will strictly insist that any candidate finite counterterm must be a local superconformal invariant, such that only non-trivial cocycles contribute to the anomalies.
Given the field content in the bulk, or equivalently the spectrum of local gauge invariant operators in the dual quantum field theory, the local superconformal invariants in any dimension can be classified. This amounts to a classification of the possible (non-trivial) contributions to the Weyl anomaly, since the Wess-Zumino consistency condition implies that the Weyl anomaly is itself a conformal invariant. For N = 1 supersymmetric theories in d = 4 this classification has been done using superspace techniques in [16]. Explicit calculations of the Weyl anomaly have been carried out in superspace in [17,18] and in components in [19]. The result that is relevant for us here is the general form of the bosonic part of the Weyl anomaly in eq. (4.3) of [19], which is the sum of the only two non-trivial superconformal invariants corresponding to our field content. 9 These superconformal invariants are E is the Euler density, which is a topological density and therefore does not receive fermionic corrections. The bosonic part of W 2 is a linear combination of the square of the Weyl tensor and F ij F ij . Notice that even though the square of the Weyl tensor and F ij F ij are separately Weyl invariant, only the specific linear combination in W 2 preserves supersymmetry. This means that there is a two-parameter family of supersymmetric renormalization schemes in this case. 10 The term S (4) in the solution (4.18) of the Hamilton-Jacobi equation is proportional to the linear combination which, as we will discuss in section 5, corresponds to the holographic Weyl anomaly for d = 4 and can be shown to be a superconformal invariant. Since the Euler density does not receive any contribution from the gravitino, it follows that the expression (4.14) determines the fermionic part of the superconformal invariant W 2 , up to quadratic order in the gravitino. To our knowledge, this is the first instance where this superconformal invariant, including the fermionic part, has been determined, at least in component language.
Combining the observations of this subsection with the general asymptotic solution of the Hamiton-Jacobi equation in (4.18), we conclude that the most general form of the local boundary counterterms compatible with supersymmetry (and parity), and which do not add trivial cocycles of the Wess-Zumino consistency conditions to the anomalies, is where s 1 and s 2 are arbitrary parameters, corresponding to a choice of renormalization scheme. In particular, the holographic renormalization scheme for d = 3 is unique, while for d = 4 there is a two-parameter family of supersymmetric schemes. Once a choice of scheme has been made, the renormalized on-shell action is defined through the limit The holographic dictionary identifies this quantity with the renormalized generating functional of the dual gauge-invariant operators, as we will review in section 5.

Fefferman-Graham asymptotic expansions from flow equations
Besides the boundary counterterms, the asymptotic solution (4.18) of the Hamilton-Jacobi equation also determines the general Fefferman-Graham expansions of the bulk fields. In particular, inserting the solution (4.18) in the expressions (3.17) for the canonical momenta, and these in turn into the Hamilton equations (3.15), results in a set of first order flow equations that can be integrated to obtain the asymptotic expansions of the fields. Up to the order necessary to obtain both the normalizable and non-normalizable modes of the asymptotic expansions, the first order equations take the forṁ where the ellipses stand for asymptotically subleading terms that only affect the asymptotic expansions beyond the normalizable mode. We have written these equations explicitly for the case d = 4, but kept the coefficients generic so that one can also apply them to the simpler case d = 3.
Notice that the unspecified quantities π j (4)b (x), π (4)i (x), π j Ψ(4) (x) correspond to the contribution to the flow equations from the undetermined finite part S (4) (or S (3) in the case d = 3) in the solution (4.18) of the Hamilton-Jacobi equation. As we will see in section 5, these quantities are directly related to the corresponding local operators in the dual field theory.
The first order equations (4.23) can be solved either recursively (see [45] for a step by step example), or by making an ansatz for the asymptotic form of the fields. The resulting asymptotic solutions for d = 4 and d = 3 are as follows.
Fefferman-Graham expansions for d = 4: The asymptotic solutions take the form where e a i (0) (x), A (0)i (x) and Ψ (0)+i (x) are arbitrary, while the remaining coefficients are given by In these expressions g (0)ij = e a i (0) e ja(0) is the boundary metric and the notations g i (2)j and g (2) stand respectively for g ik (0) g (2)kj and g ij (0) g (2)ij . Notice that the normalizable modes e a i (4) , A (4)i and Ψ (4)−i are related respectively to the undetermined canonical momenta π j (4)b (x), π (4)i (x), π j Ψ(4) (x) obtained from the finite part S (4) in (4.18). Since the latter are directly related to the dual operators, the expressions for e a i (4) , A (4)i and Ψ (4)−i in (4.25) provide the explicit relation between the dual operators and the coefficients of the Fefferman-Graham expansions. Finally, although the solution (4.18) allows one to determine the Fefferman-Graham expansions of the bosonic fields up to quadratic order in the gravitino, we only give explicitly these expansions to zero order in the gravitino in (4.25). These expansions generalize earlier results for the case of flat boundary metric and vanishing gauge field [38].
Fefferman-Graham expansions for d = 3: The asymptotic expansions for d = 3 do not contain any logarithmic terms and therefore take simpler form where e a i (0) (x), A (0)i (x) and Ψ (0)+i (x) are again arbitrary and the remaining coefficients are Notice that in this case the relation between the normalizable modes e a i (3) , A (3)i and Ψ (3)−i and the undetermined momenta π j is very simple. These asymptotic expansions are in agreement with those obtained in [47] in a similar setup.

Local operators, Ward identities and superconformal anomalies
In this section we use the boundary counterterms (4.21) in order to construct the holographic dictionary for the theory (2.1). In particular, we define the local operators dual to the bulk fields and derive the Ward identities these operators satisfy, including all bosonic and fermionic anomalies. These anomalies play a central role in our analysis of supersymmetric backgrounds in section 6.
The operators dual to the bulk fields correspond to the renormalized canonical momenta, namely These relations are the holographic version of the local renormalization group definition of local operators [48], which is why we omit the bracket notation · . Each of these operators is a function of arbitrary sources and so correlation functions can be computed by further differentiating these operators with respect to the sources. Notice that for d = 4, these operators depend on the value of the parameters s 1 and s 2 that parameterize the supersymmetric renormalization schemes. Specifically, the definition (5.1) of the local operators implies that they are identified respectively with the quantities π j (3)b , π (3)i and π j Ψ(3) appearing in the Fefferman-Graham expansions for d = 3, while for d = 4 they correspond to π j (4)b , π (4)i and π j Ψ(4) , up to terms proportional to the parameters s 1 and s 2 coming from the finite local counterterms parameterizing the renormalization scheme.
The Ward identities that the local operators (5.1) satisfy can be derived by inserting the covariant expansion (4.2) in the constraints (3.12) and (3.7) and isolating the terms of dilatation weight zero. This leads to the identities where the expressions on the r.h.s. vanish identically for d = 3, while for d = 4 are local in terms of the induced fields and take the form 11 Combining the definitions (5.1) of the dual operators with the identities (5.2) (taking into account the fact that the finite local counterterms in (4.21) parameterizing the supersymmetric renormalization scheme by construction obey the identities (5.2) without the anomaly terms) and removing the radial cutoff results in the Ward identities where the local terms on the r.h.s. of these identities are identical to the expressions on the radial cutoff in (5.3), except that the induced fields are appropriately replaced with the corresponding sources on the boundary of AdS. Analogous Ward identities, but with flat metric g (0)ij , were obtained for the N = 1 theory dual to the Klebanov-Strassler background in [49]. It cannot be overemphasized that the Ward identities (5.4) and the related anomalies (5.3) apply to generic N = 2 superconformal quantum field theories in three and generic N = 1 superconformal quantum field theories four dimensions, except that the form of the anomalies (5.3) in the four dimensional case holds for theories with a = c only, and the value of these anomaly coefficients may be different in a generic theory. Crucially, the form of the Ward identities (5.4) applies to any such theory, irrespectively of whether it admits a holographic dual. In particular, the alternative derivation of the Ward identities discussed in the following subsection is entirely field theoretic and allows one to extrapolate the results of section (6) beyond holographic theories.

Boundary supersymmetry and current anomalies
The Ward identities (5.4) reflect the local bulk symmetries surviving on the conformal boundary, including those broken by anomalies arising from the regularization and renormalization procedure. An alternative derivation of these identities, therefore, may be achieved by identifying the bulk symmetries acting on the boundary sources and considering how the renormalized generating functional transforms. From a bulk perspective, such symmetries correspond to all local bulk transformations preserving the Fefferman-Graham gauge (A.8). In the case of bulk diffeomorphisms these are known as Penrose-Brown-Henneaux (PBH) transformations [50,51,52]. In appendix B we derive the most general bulk symmetry transformations of the supergravity action (2.1) preserving the Fefferman-Graham gauge. We will refer to these as generalized Penrose-Brown-Henneaux transformations.
One of the main results of the analysis in appendix B is the identification of the local symmetries on the conformal boundary and their action on the field theory sources. Specifically, we find that the local symmetry transformations on the boundary are generated by the arbitrary bosonic parameters σ(x), ξ i o (x), λ ab o (x), as well as the fermionic variables o+ (x) and o− (x), all of which are functions of the boundary coordinates only. Removing the radial cutoff in the asymptotic relations (B.16) we determine that these transformations act on the field theory sources as where from (3.16) and (4.10) Notice that, although these transformations are derived holographically here, they should apply to any N = 2 superconformal theory in d = 3 and any N = 1 superconformal theory d = 4. Given the general variation of the renormalized on-shell action the Noether procedure for the transformations (5.5) leads to an alternative derivation of the Ward identities (5.4). In particular, inserting the transformations (5.5) for the sources in (5.7) and parameterizing the anomalies such that reproduces the Ward identities (5.4). However, the anomalies must be computed independently in this case, either holographically as above, or through field theory arguments. Note that the term A M (0)i in (5.4a) cancels out in the variation of the generating function (5.8), implying that the theory is invariant under local diffeomorphisms. Even though the field theory sources transform as tensors under the local symmetries, the corresponding operators do not, due to the anomalies. The simplest way to obtain the transformation of the currents under the local symmetries is by introducing a suitable Poisson bracket on the symplectic space of sources and local operators as discussed in appendix B. In particular, the general transformation of the current densities is given in (B.26), with the anomalous part of the transformation shown in terms of the anomaly densities in (B.27). For our purpose of studying supersymmetric backgrounds in section 6 it suffices to determine explicitly the transformations of the supercurrent under boundary supersymmetry and superWeyl transformations, generated respectively by o+ (x) and o− (x). The general transformation identities (B.26) for the currents and the explicit form of the anomalies in (5.3) determine Notice that, besides the tensorial part involving the currents, these transformations contain additional local terms due to the anomalies. This result will be central to our analysis of supersymmetric backgrounds in the next section. Another observation that will play an important role in the subsequent analysis is that the anomalous supersymmetry transformations (5.9) of the supercurrent can be obtained in a third way, namely from the bulk supersymmetry transformation of the gravitino. Finally, it should be stressed that the anomalous transformations (5.9) of the supercurrent are independent of the choice of supersymmetric renormalization scheme, i.e. of the choice of the parameters s 1 and s 2 in (4.21).

Fermionic anomalies and supersymmetry breaking
The anomalous transformation (5.9) of the supercurrent under local supersymmetry and superWeyl transformations is our main result. In this section we explore some of the consequences of this result for N = 1 superconformal theories in four dimensions. In particular, we show that even if a field theory background admits a Killing spinor and all anomalies are numerically vanishing, supersymmetry can still be broken by the anomalous transformation of the supercurrent.

Field theory backgrounds admitting conformal Killing spinors
In order to explore the consequences of the anomalous supercurrent transformation we consider a simple class of field theory backgrounds admitting Killing spinors that were originally obtained in [4,5] and were the subject of the recent analysis of [27]. Setting the AdS 5 radius to = 1, this class of rigid four-dimensional supersymmetric field theory backgrounds takes the form 12 ds 2 (0) = −dt 2 + (dψ + a) 2 + 4e w dzdz, where w(z,z) is a local function on the compact Riemann surface Σ 2 parameterized by the complex coordinates (z,z), a = a z (z,z)dz + az(z,z)dz is a local one-form on Σ 2 satisfying the condition da = iue w dz ∧ dz, (6.2) for a globally defined function u(z,z) on Σ 2 , and ∂ ψ is a nowhere vanishing Killing vector. The real constants γ and γ , as well as the function λ = λ(z,z) on Σ 2 are locally pure gauge, but their values are determined by the requirement that A (0) be a globally defined one-form. In particular, global considerations require that γ = 0 and γ be proportional to the first Chern class of Σ 2 . For an extensive discussion of the global properties the backgrounds (6.1) and explicit expressions for these parameters in terms of topological invariants we refer the reader to [27]. In the following it will be convenient to slightly reparameterize the background (6.1) in order to solve explicitly the constraint (6.2). In particular, writing the components of the one-form a in terms of a local function µ(z,z) as the condition (6.2) is automatically satisfied provided u(z,z) is expressed in terms of w(z,z) and µ(z,z) through the relation u = e −w ∂ z ∂zµ. (6.4) This allows us to parameterize the supersymmetric backgrounds (6.1) in terms of the two unconstrained functions w(z,z) and µ(z,z) as Conformal Killing spinors and vectors A Killing vector K i of a generic bosonic background specified by the sources e a (0)i and A (0)i corresponds to a diffeomorphism that leaves the sources invariant, up to local frame rotations and U (1) gauge transformations. The general local transformations of the sources (5.5) imply that in terms of the metric and the gauge invariant fieldstrength the conditions for K i to be a Killing vector are where L K denotes the Lie derivative with respect to K i . Similarly, a conformal Killing vector corresponds to a combined diffeomorphism and Weyl transformation that leaves the sources invariant up to the same transformations. In this case (5.5) give and so the Weyl factor σ is determined in terms of the vector K i as Killing spinors and conformal Killing spinors on a bosonic background are defined analogously. Namely, from the transformations (5.5) follows that a conformal Killing spinor ζ + satisfies 13 where Killing spinors are the subclass of such spinors having in addition ζ − = 0. It is a straightforward exercise to show that for any conformal Killing spinor ζ + , the spinor bilinear is a conformal Killing vector with Weyl factor In particular, if ζ + is a Killing spinor, then K i ζ is a Killing vector. In order to determine the Killing spinors that the backgrounds (6.5) admit, it is necessary to first make a choice of vielbein for the background metric, and to specify a basis for the gamma matrices. In the following we will choose the vielbein components as 13 Killing spinors correspond to special supersymmetry transformations and so the radiality of such spinors is unambiguously determined by the generalized PBH transformations (5.5). Our spinor conventions therefore differ from other common choices in the literature, and in particular from those in [27]. and the gamma matrix basis where σ a are the Pauli matrices. To avoid potential confusion with the complex conjugate coordinatez, in this section we denote frame indices with a dot rather than an overbar, as in the rest of the paper. With the above choice of vielbein and gamma matrices, one finds that the spinor ζ − of any conformal Killing spinor of the background (6.5) is given by Solving the conformal Killing spinor condition on a generic background of the form (6.5) with u = 0 we find that there is a unique solution given by Note that these spinors are globally well defined provided γ = 0 [27]. With this normalization the corresponding conformal Killing vector (6.11) takes the form In the special case u = 0, however, there are two independent solutions to the conformal Killing spinor conditions, namely In the following we will focus on the generic case admitting only the conformal Killing spinor (6.16).

Supercurrent anomalies and supersymmetric vacua
Restricting the local supersymmetry and superWeyl transformations of the supercurrent in (5.9) to the conformal Killing spinor (6.16) results in the transformation of the supercurrent under rigid supersymmetry, namely (6.20) where we have used the fact that on a bosonic background the Ward identity (5.4f) implies that the stress tensor T ij is symmetric, as well as the gamma matrix identity (6.21) and the conformal Killing spinor relation (6.9). Although all anomalies in the superconformal Ward identities (5.4) numerically vanish on the class of supersymmetric backgrounds (6.5), the anomalous term in the supercurrent transformation (6.20) under rigid supersymmetry is not vanishing for a generic background of the form (6.5).
The fact that the supercurrent transforms anomalously under rigid supersymmetry on the backgrounds (6.5) has important implications for the supersymmetric vacua of theories defined on such backgrounds. As we remarked earlier, the transformations (5.9) of the supercurrent under local supersymmetry and superWeyl transformations can alternatively be derived from the bulk supersymmetry transformation of the gravitino. It follows that for supersymmetric vacua, corresponding to supergravity solutions satisfying the bulk BPS equations, the expectation value of the supersymmetry variation of the supercurrent vanishes, namely This is commonly used to argue that the linear combination of the stress tensor and the supercurrent is Q-exact. However, the anomaly in the supercurrent transformation under rigid supersymmetry invalidates this argument and hence, the operator (6.23) is Q-exact only on the subclass of backgrounds of the form (6.5) where the supercurrent anomaly is numerically zero. A sufficient but not strictly necessary condition for the supercurrent anomaly to vanish is but clearly this is a too restrictive condition. By a careful analysis of the supercurrent anomaly in (6.20) one can determine the most general subclass of backgrounds of the form (6.5) for which the supercurrent anomaly vanishes, but we will not address this question here.

The w and µ deformations
The fact that the operator (6.23) is not in general Q-exact due to the supercurrent anomaly under rigid supersymmetry implies that the argument first presented in [28,29,30], showing that the supersymmetric partition function on backgrounds of the form (6.5) is independent of the functions w and µ, needs revisiting. In this subsection we reexamine this argument using the anomalous supercurrent transformation (6.20), and we provide a resolution to an apparent paradox in the recent analysis of [27]. We would like to emphasize that the subsequent analysis relies only on the superconformal Ward identities and the anomalous transformation of the supercurrent and, hence, it is applicable to any N = 1 superconformal field theory on backgrounds of the form (6.5), irrespectively of whether they admit a holographic dual.
The general variation (5.7) of the generating functional implies that under a bosonic variation Since the supersymmetric backgrounds (6.5) are parameterized by the two unconstrained functions w(z,z) and µ(z,z), we would like to examine the dependence of the supersymmetric partition function on these functions. As we now show, the variation of the partition function under infinitesimal deformations of the functions w(z,z) and µ(z,z) can be expressed as a linear combination of different components of the homogeneous part of the supercurrent transformation (6.20), i.e. the part proportional to the stress tensor and the R-symmetry current given in (6.23). This observation implies that the transformation of the supersymmetric partition function under w and µ deformations, for which the total variation of the supercurrent vanishes due to the identity (6.22), is determined entirely by the supercurrent anomaly!
The w deformation An infinitesimal deformation of the function w, keeping µ fixed, corresponds to the following transformations of the bosonic sources: with the variation of all other components vanishing.
The µ deformation Similarly, an infinitesimal µ deformation, keeping w fixed, corresponds to the source variations We next insert these transformations of the sources in the identity (6.25) and, by dropping various total derivative terms in the complex coordinates z andz, we express these variations in terms of certain components of the homogeneous transformation of the supercurrent under rigid supersymmetry. Specifically, for the w deformation we obtain The subscript 1 or 2 in the variations of the supercurrent denotes the relevant spinorial component. Moreover, from the second to the third equality we have used the R-symmetry Ward identity D (0)i J i = 0, since the R-symmetry anomaly vanishes on the backgrounds (6.5), as well as the identity (6.22) for the supersymmetric partition function. In the last step we have evaluated the indicated components of the supercurrent anomaly in (6.20) on the backgrounds (6.5). Remarkably, the final result agrees completely with that obtained in eq. (4.35) of [27], provided the "scheme" parameters ς and ς there are set to zero. 14 Since these parameters do not multiply superconformal invariants, they contribute trivial cocycles to the anomalies, and hence, should be set to zero. As we discussed in section 4, the supersymmetric schemes correspond instead to the parameters s 1 and s 2 in (4.21) and neither of these affects the supercurrent anomaly, which is scheme independent. Similarly, the µ variation of the partition function can also be related to a certain linear combination of the supercurrent variation under rigid superymmetry, namely Once again, this result agrees with the corresponding expression in eq. (4.36) of [27], after setting the parameters ς and ς there to zero. We would like to emphasize that the above argument for the dependence of the supersymmetric partition function on w and µ is fundamentally different from that in [27]. The derivation in [27] relies exclusively on the bosonic operators and their sources. In particular, the variation of the supersymmetric partition function is obtained by explicitly evaluating the bulk on-shell action on supersymmetric supergravity solutions and it is therefore inherently a holographic derivation. The argument we provide above, however, does not require the evaluation of the partition function at any stage. Instead, by including the supercurrent in the analysis, we express the w and µ variation of the supersymmetric partition function in terms of the supercurrent transformation under rigid supersymmetry. This relation is a direct consequence of the superconformal Ward identities and hence, our argument is in fact a field theory one! As we stressed a number of times, in our entire analysis we only rely on holography in order to derive the form of the superconformal Ward identities. These identities, however, apply more generally to SCFTs, with or without a holographic dual.
An immediate conclusion that follows from the above derivation is that, contrary to earlier claims, the supersymmetric field theory partition function is in fact not invariant under w and µ deformations around a generic background of the form (6.5). Moreover, holographic renormalization reproduces precisely this field theory property. In particular, the anomalous transformation of the supersymmetric partition function under w and µ deformations is scheme independent on both the field theory and holographic sides (provided the scheme preserves supersymmetry), and, since it is a direct consequence of a non-trivial supercurrent anomaly, it cannot possibly be removed by a local covariant boundary counterterm. This anomalous transformation is physical and so there is no need to add non-covariant counterterms as is done in [27], which explicitly break other symmetries of the theory! As we now explain, these observations also imply that the BPS relation between the conserved charges on supersymmetric vacua is also anomalous!

Conserved charges and the BPS relation
The conserved charges are a consequence of the Ward identities (5.4) [39]. The presence of anomalies in the Ward identities renders the argument leading to the conserved charges less trivial in general [53], but for the supersymmetric backgrounds (6.5) we are interested in here these anomalies numerically vanish.
Electric charge The R-symmetry Ward identity (5.4b) implies that the quantity where the integral is over a constant time slice of the geometry (6.5), is conserved for any background field A (0)i if ω = 1, or for any value of ω if the R-symmetry anomaly is numerically zero, as is the case for the backgrounds (6.5). Notice that we have introduced the notation · here in order to emphasize that the charges are associated with a given state of the theory and so involve the expectation value of the operators in a specific state. Moreover, the overall factor 1/ √ 3 in the definition of the electric charge is included in order to compensate the normalization of the background gauge field in (6.5) and get a canonical BPS relation between the charges later on. The Contracting this identity with a conformal Killing vector K i and utilizing the trace (5.4d) and R-symmetry (5.4b) Ward identities, one finds that the quantity is conserved provided both the Weyl anomaly A W (0) and A M (0)i vanish numerically, which is the case for the supersymmetric backgrounds (6.5). Note that these charges are not gauge invariant except when ω = −2, but they can be made gauge invariant for any ω by a minor modification [53].
Supercharges Finally, the supercharges associated with a conformal Killing spinor ζ + are These charges are conserved provided the supersymmetry and superWeyl anomalies vanish, and they are of course zero in a bosonic background.

BPS relation
The transformation of the supercurrent under rigid supersymmetry in (6.20) leads to the BPS relation between the conserved charges in supersymmetric vacua. In particular, contracting the identity (6.20) from the left with iζ + and taking the expectation value in a supersymmetric vacuum we obtain where and we have used the Killing spinor property (6.9). Moreover, from the conformal Killing spinor solution (6.16) we determine and which allow us to rewrite (6.35) in the form Integrating this identity over a constant time slice of the geometries (6.5) we obtain the relation where, Notice that, as the conformal Killing charges (6.33), the anomaly charge is in general gauge invariant only for ω = −2. However, the identity (6.40) is gauge invariant and holds for any ω. Identifying the mass and angular momentum respectively with the conserved charges (6.42) leads to the BPS condition In particular, the supercurrent anomaly gives rise to an anomaly in the BPS condition! This observation provides a resolution to another apparent paradox discussed in [27], as well as the earlier works [31,32,33,34,35], concerning the supersymmetric Casimir energy, which we will discuss momentarily. Interestingly, the anomalous BPS condition (6.43) seems to be compatible with the recent conjecture [57] for non-supersymmetric AdS vacua. 15 BPS states and Casimir charges In order to compute the supersymmetric Casimir charges it is necessary to determine the expectation values of the stress tensor and the R-current in general supersymmetric vacua. These expectation values must simultaneously satisfy the supersymmetry condition (6.22), as well as the trace, diffeomorphism and R-symmetry Ward identities. These conditions together account for 6 + 1 + 4 + 1 = 12 constraints, while the stress tensor and the R-current together have 10 + 4 = 14 components. It follows that the general form of the one-point functions for supersymmetric vacua must be parameterized in terms of two arbitrary functions. The supersymmetry condition (6.22) and the Ward identities can be solved systematically to obtain the following general solution for the supersymmetric one-point functions: 15 I thank Gavin Hartnett for suggesting this connection. and As anticipated, these supersymmetric expectation values are parameterized by the two arbitrary functions Υ(z, z) and Φ(z, z). This solution for the supersymmetric one-point functions can be generalized further by replacing where ν is an arbitrary constant and where W (0)ijkl is the Weyl tensor of the metric g (0)ij , which for backgrounds of the form (6.5) satisfies the identity The expressions (6.47) correspond respectively to the derivative of where W 2 (0) is the supersymmetric Weyl squared density defined in eq. (4.19), with respect to the metric g (0)ij and the gauge field A (0)i . Notice that although W 2 (0) vanishes on the backgrounds (6.5), its derivatives T i j loc and J i loc are not generically zero. However, it can be shown that they only contribute total derivative terms to the one-point functions of the currents, and so do not affect the value of the conserved charges. This fact also implies that the conserved charges are independent of the choice of supersymmetric renormalization scheme, parameterized by the two parameters s 1 and s 2 in (4.21).
The general solution (6.44) and (6.45) for the supersymmetric one-point functions allows us to compute explicitly the conserved charges in terms of the functions Υ(z, z) and Φ(z, z). We find As expected, these charges satisfy the BPS relation (6.43) for any value of the parameter ω and for arbitrary functions Υ(z, z) and Φ(z, z). Moreover, they apply to any BPS state, including states corresponding to supersymmetric black holes in the bulk.
We would now like to evaluate the conserved charges (6.50) explicitly in the global supersymmetric vacuum of the theory. This vacuum corresponds to a smooth horizonless geometry in the bulk, which allows one to determine the functions Υ vac (z, z) and Φ vac (z, z) explicitly. As was pointed out in [27], for smooth horizonless bulk solutions the Page R-charge vanishes identically, i.e. Q ω=1 e = 0, which implies that the time component of is a total derivative. This condition, together with the general supersymmetric solution (6.45), determines the function Φ(z, z) to be Φ vac (z, z) = 1 3 γR 2d − 1 24 (u 3 − 4uR 2d ). (6.52) In order to determine the function Υ(z, z) in the global supersymmetric vacuum we need an additional condition. In [27] this additional condition was that the Ricci potential be a globally defined one-form in the bulk. However, physically this condition is equivalent to requiring that the entropy vanishes, which is not a priori obvious how to express as a condition on the one-point functions of the stress tensor and the R-current. To circumvent this, below we determine the function Υ vac (z, z) by requiring the value of the free energy to agree with that computed in [27] (without the "new counterterms") through the global Ricci potential condition. We then manage to express the value of Υ vac (z, z) as a condition on the one-point functions of the stress tensor and the R-current. The resulting condition provides an expression for the supersymmetric entropy in terms of the local one-point functions. For the sake of presentation, however, we find it convenient to start with the zero entropy condition in terms of the one-point functions, and justify it a posteriori by the fact that it leads to the correct value of Υ vac (z, z).
In order to write down an expression for the entropy, or to be able to relate the free energy to the conserved charges, we need to identify the timelike Killing vector that would become null on the Killing horizon. The form of the graviphoton in (6.5) implies that for generic u(z, z) the null generator of the horizon is uniquely determined to be is constant. Given this timelike Killing vector, we conjecture that the entropy of any supersymmetric state can be expressed in terms of the expectation values of the stress tensor and the R-current as where is the unit vector orthogonal to χ. Notice that the integral is done over the surface orthogonal to χ ⊥ and the relative sign between the stress tensor and the R-current is opposite that appearing in the conserved charges (6.33). Obtaining a first principles derivation of this expression is highly desirable, but here we will be content with the fact that it gives the correct expression for Υ vac (z, z).
Having determined the functions Υ vac (z, z) and Φ vac (z, z) for the supersymmetric vacuum, we can evaluate explicitly the charges (6.50), which in this case are identified with the Casimir charges:

Supersymmetric partition function
We conclude with some remarks regarding the value of the supersymmetric partition function. On general grounds, the Euclidean partition function I is proportional to the Gibbs free energy, namely where β is the perimeter of the Euclidean time circle, Ω is the angular velocity, and Φ e is the electric potential. This relation applies to both supersymmetric and non-supersymmetric vacua and is renormalization scheme independent. In particular, the scheme dependence of the partition function compensates that of the mass M , while the remaining charges are scheme independent. This relation applies to any field theory, independently of whether it admits a holographic dual, but a general holographic proof was provided in [39] by combining the Noether current approach to entropy [58] and the holographic definition of the conserved charges. In order to evaluate the supersymmetric partition function we use the Casimir charges (6.59) and set S = 0. Moreover, from the Killing vector (6.53) follows that the angular momentum and electric potential are given respectively by Using these expressions we obtain This result is independent of the value of ω and of the choice of supersymmetric renormalization scheme. Moreover, it agrees with the free energy obtained in [27] (without the "new counterterms" and with opposite sign for the term proportional to (γ + 2γ ), due to our different orientation conventions, as discussed in footnote 14).

Appendices A Conventions and radial ADM decomposition
In this appendix we define our conventions and collect several identities that are extensively used in this paper.

A.2 Gamma matrix conventions and identities
The Gamma matrices Γ α with flat frame bundle indices satisfy the Clifford algebra The Gamma matrices with coordinate indices are defined using the inverse vielbein E µ α as and satisfy We do not need the explicit matrix representation of the Clifford algebra here, except for the calculation in section 6, where we specify the gamma matrix representation used in (6.14). As in [44], we will assume that the representation is Hermitian so that We also adopt the standard notation for totally antisymmetric products of Gamma matrices where the indices are antisymmetrized with weight one. Several identities these products satisfy in D-dimensions can be found e.g. in section 3 of [44]. Some of the identities that we utilize here are

A.3 Radial ADM decomposition
The radial ADM decomposition of the dynamical variables consists in choosing a suitable radial coordinate r emanating from the conformal boundary and describing the bulk space as a foliation by constant r-slices, which we denote by Σ r . Accordingly, all fields with tensor indices are decomposed along the radial and transverse directions as where γ ij , A i and Ψ i are dynamical induced fields on the slices Σ r , while N , N i , a and Ψ r , are non-dynamical Lagrange multipliers conjugate to first class constraints.
Fefferman-Graham gauge Since N , N i , a and Ψ r , are non-dynamical Lagrange multipliers, once the corresponding first class constraints have been derived, we can set them to any fixed value that is convenient. A choice that is particularly useful for reading off the holographic dictionary is the Fefferman-Graham (FG) gauge Whenever we gauge-fix the Lagrange multipliers in this paper we adopt this choice of gauge.
Vielbein and gamma matrices The vielbeins compatible with the ADM metric (A.7) can be expressed in terms of a (d + 1)-vector n α and a (d + 1) × d matrix e α j as The inverse vielbeins are then given by and satisfy the orthogonality relation In order to simplify the calculations we sometimes choose a convenient vielbein frame so that and e a i becomes the vielbein on the slice Σ r .
Gamma matrix decomposition and radiality The decomposition of the vielbein allows us to decompose the gamma matrices in radial and transverse components as follows: (A.14) These relations imply that Moreover, the Hermiticity property (A.4) translates into the identities while the Clifford algebra (A.3) becomes The fact that Γ anticommutes with all gamma matrices Γ i allows us to introduce the radiality projectors and the corresponding radially projected spinors [40] Note that the radiality projectors are independent of the spacetime dimension, but coincide with the usual chirality projectors for odd D = d + 1. Moreover, splitting the spinors according to their radiality is necessary not only for formulating a consistent Hamiltonian description of the fermion dynamics [41], but also in order to construct the asymptotic Fefferman-Graham expansions, since fermions of different radiality have different asymptotic behavior.
Levi-Civita tensor The Levi-Civita tensor in the (d + 1)-dimensional bulk is given by In the coordinate system (A.7) it takes the form where i 1 ···i d is the Levi-Civita tensor on Σ r and The Levi-Civita tensor is related with the antisymmetric products of the gamma matrices through the identity Γ µ 1 µ 2 ...µ D ∝ µ 1 µ 2 ...µ D .
(A. 22) In particular, for D = d + 1 with d even (see (3.31) in [44]) It follows that Radial decomposition of the Christoffel and spin connections We further need the ADM decomposition of the Christoffel symbol and of the spin connection. A straightforward calculation using the metric (A.7) determines the bulk Christoffel symbol: is the extrinsic curvature of the radial slice Σ r and K ≡ γ ij K ij . The torsion free spin connection is expressed in terms of the vielbein as Using the above decompositions of the vielbein and of the Christoffel symbol we determine the components of the spin connection (see also (88) and (89) in [41]): Choosing the convenient frame (A.13), these identities allow us to write It follows that the components of the covariant derivative acting on the gravitino take the form are the covariant derivatives of, respectively, a vector and a scalar spinor on the radial slice Σ r .

B Generalized PBH transformations
Penrose-Brown-Henneaux (PBH) transformations [50,51,52] are bulk diffeomorphisms that preserve the Fefferman-Graham gauge (A.8). In addition to diffeomorphisms, the supergravity action (2.1) is invariant under local supersymmetry and U (1) gauge transformations. Moreover, local SO(1, d) frame rotations not only leave the action (2.1) invariant, but also automatically preserve the FG gauge (A.8). However, as we discussed in appendix A, it is particularly useful to make a specific choice of frame defined by the conditions (A.13). In combination with the FG gauge (A.8), these conditions amount to requiring that We will refer to this gauge fixing condition as the strong Fefferman-Graham gauge. It is equivalent to the FG gauge (A.8), but also partially gauge fixes the frame rotations. In the remaining of this appendix we will determine the most general local bulk transformations that preserve the gauge conditions (B.1). Such transformations generalize the PBH diffeomorphisms and have a central role in constructing the holographic dictionary, and in particular, in determining how supersymmetry acts on the boundary. We start by considering the transformation of all bulk fields under infinitesimal diffeomorphisms ξ µ , frame rotations λ αβ = −λ βα , U (1) gauge transformations θ, as well as local supersymmetry transformations . The transformation of the fields under local supersymmetry transformations is given in (2.6). Under bulk diffeomorphisms while under local frame rotations and under U (1) gauge transformations Requiring that the total transformation of all the fields preserves the conditions (B.1) leads to the following set of differential equations for the local parameters: The last equation can be simplified by decomposing the gravitino and the spinor parameter using the radiality projectors (A.18), namely This leads to the two spinor equationṡ Moreover, combining the second and third equations in (B.5) using λr b = −λ ar η ab , gives The above equations can be solved to determine the general form of the local transformations preserving the strong FG gauge conditions (B.1). The general solution contains a number of integration functions and takes the form (see also section C of [47]) In particular, although the gauge fixing conditions (B.1) allow λ ar (r, x) to be arbitrary, preserving the asymptotic form of the vielbein e a i requires that λ ar (r, x) ∼ λ a or (x), up to subleading terms that remain arbitrary, but do not affect any physical observable. Moreover, the expressions for ± look formal since ± enter on the r.h.s. of the last two equations in (B.9) as well. However, this dependence on the r.h.s. is asymptotically subleading, and so are the last two terms in the expression for Ξ ± in (B.10). The expressions for ± , together with that for Ξ ± in (B.10), therefore, allow one to recursively determine the asymptotic solution of the spinors ± . To leading order asymptotically the general solution (B.9) takes the form Notice that each local transformation is sourced by an arbitrary function of the transverse coordinates. As is discussed in the main part of the paper, these parameters correspond to the local symmetries of the dual field theory in a background of arbitrary sources. To determine how the local sources transform under these local symmetries we next consider how the generalized PBH transformations act on the induced fields on Σ r , in the limit that r → ∞.
Under local bulk transformations that preserve the strong FG gauge (B.1) the induced fields on the radial slices Σ r transform as δ ξ,λ,θ, e a i = σė a i + ξ j ∂ j e a i + (∂ i ξ j )e a j − λ a b e b i + Using the leading asymptotic form of the induced fields in eq. (B.11), and of the local symmetry parameters in eq. (B.12), we find that to leading order asymptotically δ ξ,λ,θ, e a i ∼ σ e a i + ξ j o ∂ j e a i + (∂ i ξ j o )e a j − λ a o b e b i + 1 2 ( + Γ a Ψ +i − Ψ +i Γ a + ), (B.16a) Taking the limit r → ∞, these expressions result in the transformations of the field theory sources given in eq. (5.5). Finally, projecting the gravitino transformation in (B.14) with the P − radiality projector gives From eq. (3.16) follows that this corresponds to the leading asymptotic transformation of the canonical momentum π i Ψ .

B.1 Ward identities as generators of generalized PBH transformations
In section 5 we argued that the generalized PBH transformations are directly related with the Ward identities (5.4), which were obtained from the first class constraints (3.13). Specifically, varying the renormalized on-shell action (4.22) with respect to arbitrary PBH transformations of the form (5.5) provides an alternative derivation of the Ward identities. As we now demonstrate, the connection between the Ward identities and generalized PBH transformations admits an elegant formulation on the symplectic space of sources and local operators, which can be employed in order to obtain the transformations of the local operators under PBH transformations. Inserting the Fefferman-Graham expansions in the symplectic form associated with the bulk action (2.1) [59,60], it can be shown that the field theory sources and the corresponding operators defined through (5.1) parameterize a symplectic manifold equipped with the symplectic 2-form The symplectic form (B.19) allows us to introduce the Poisson bracket [61] {A, B} PB = d d x δA δe a for any functions A and B on the space of sources and local operators. It should be emphasized that, although this result is obtained here holographically from the bulk theory, it applies to any local quantum field theory where the local operators are defined as in the local renormalization group [48]. Indeed, the subsequent analysis in this appendix applies to any local quantum field theory [45], independently of whether it admits a holographic dual.
The Ward identities (5.4) imply that the local functions where e (0) ≡ −g (0) and These anomalous contributions to the transformation of the local operators are central to the analysis of the symmetries preserved by the quantum field theory vacua. Finally, the Poisson bracket of the generating functions (B.23) can be used to determine the infinite dimensional superalgebra of generalized PBH transformations, as is done e.g. in eq. (3.2.35) of [61] for local Weyl transformations and diffeomorphisms. However, in order for the algebra to close in the present context, it is necessary to know the bulk supergravity action and the corresponding supersymmetry transformations to all orders in the gravitino.
Combining these transformations we find that the total supersymmetry variation is a pure boundary term provided the constants c 1 , c 2 , c 3 , c 4 take the values , c 4 = ±1 while the gauge coupling is given by The signs in these expressions are correlated so that there are only two possible choices: either the upper sign or the lower sign must be chosen for all constants. Choosing the upper sign leads to the action (2.1) and the supersymmetry variations (2.6). Note also that the constant c 1 is non-zero only for d = 4, i.e. for D = 5.