General Supersymmetric AdS 5 Black Holes with Squashed Boundary

We present a new family of asymptotically locally AdS5 squashed supersymmetric black hole solutions of Fayet-Iliopoulos gauged N = 2, D = 5 supergravity with two vector multiplets that have a natural uplift to type IIB supergravity. Our new family of black holes is characterized by three parameters, of which two control the horizon geometry while the latter regulates the squashing at the boundary. We evaluate the main physical properties of the family of solutions using holographic renormalization and find that the entropy is independent on the squashing and it is reproduced by using the angular momentum and the Page charges. In previously known solutions Page and holographic charges are equal, due to the vanishing of the Chern-Simons term that here, instead, is relevant. This result suggests that for asymptotically locally AdS5 solutions we should refer to the Page charges to describe the thermodynamics of the system. e-mails: alessandro.bombini@pd.infn.it , lorenzo.papini@pd.infn.it ar X iv :1 90 3. 00 02 1v 1 [ he pth ] 2 8 Fe b 20 19


Introduction
One of the main outstanding achievements of string theory is the description of asymptotically flat black holes in terms of microscopic constituents such as strings and branes [1][2][3][4][5][6], by which is possible to compute their entropy by microstate counting, reproducing correctly the Bekenstein-Hawking entropy as a statistical Boltzmann entropy [7,8].
When the AdS/CFT correspondence was conjectured [9][10][11] it seemed very natural to use it to provide a description of asymptotically Anti-de-Sitter (AdS) black holes in terms of microscopical states in a dual quantum field theory, thus extending the results obtained for asymptotically flat black holes. In particular, since the original correspondence is AdS 5 /CFT 4 , it appeared natural to attack the problem starting with supersymmetric asymptotically AdS 5 black holes. However, although such black holes were constructed fifteen years ago [12] and many generalization immediately followed [13][14][15][16], all the first attempts to provide an interpretation of their entropy in terms of a quantum field theory computation were basically unsuccessful [17] and this problem has remained a puzzle for many years.
The same problem in one lower dimension has been solved four years ago starting in [18]. There, a certain class of supersymmetric asymptotically locally AdS 4 black holes is analyzed and their entropy is obtained by a localization computation in the dual ABJM field theory. More precisely, the authors managed to reproduce the black hole entropy by Legendre transforming the large N contribution of the topologically twisted index of the ABJM theory, introduced in [19]. The field theoretic computation performed in [18] and related papers, as well as the corresponding interpretation on the gravity side discussed for example in [20,21], shade light on many non-trivial steps one must follow in order to reproduce the entropy; it was then reasonable to expect that these results could be inspiring to solve the problem also in five-dimensions.
A solution to the five-dimensional enigma has finally been proposed in [22][23][24] (see also [25] for previous progress) by reconsidering the field theory partition function dual to the black hole allowing for complex values of the fugacities. The approaches used in the above-mentioned papers are quite different and the relation between them remains as an interesting puzzle.
In each of these field theoretic computations it is crucial to understand which are the field theory states that contribute to the entropy. Further information on this may be collected by studying black holes which are not globally asymptotically AdS, but just locally. Such black holes may be obtained by deforming the conformal boundary and some solutions with these characteristics have been constructed in minimal five-dimensional gauged supergravity [26,27] and in five-dimensional Fayet-Iliopoulos gauged supergravity [28]. In all these papers, the authors considered a cohomogeneity-one ansatz with local SU(2) × U(1) × U(1) symmetry and managed to obtain supersymmetric black holes with a non conformally-flat boundary geometry containing a squashed three-sphere. These solutions are thus asymptotically locally AdS 5 (AlAdS 5 ) rather than asymptotically AdS 5 . In the squashed solution of minimal gauged supergravity [26,27] the geometry of the event horizon results completely frozen so that the entropy is uniquely fixed; instead the squashing at the boundary can assume any value. Therefore this AlAdS 5 solution behaves differently from the general asymptotically AdS 5 one of [12], since in the latter the entropy is controlled by the same parameter regulating the horizon geometry and can therefore vary. The minimal gauged supergravity solution has been generalized in [28], where the authors constructed a two-parameters family of squashsed black holes in Fayet-Iliopoulos five-dimensional gauged supergravity with an arbitrary number of vector multiplets. Of the two parameters on which the solution depends, one controls the horizon geometry, which is thus not completely frozen anymore, while the other one determines the squashing at the boundary. The entropy of the solution is regulated by only one parameter and behaves again differently with respect to the general solution of [13], where the entropy is controlled by three parameters. Although quite general and valid for any number of vector multiplets, the solution of [28] is obtained by imposing a particular ansatz on the scalar fields which constraints all their components orthogonal to the scalar vacuum expectation values in the supersymmetric AdS 5 vacuum to be the same. In principle each of such components is controlled by a different function; the authors impose their ansatz by requiring that all those functions are the same.
The main aim of the present paper is to construct general squashed solutions in five-dimensional Fayet-Iliopoulos gauged supergravity with two vector multiplets without imposing any ansatz. This means that we will abandon the restrictive conditions imposed in [28] and let the scalar fields unconstrained. We will look for solutions with two vector multiplets only because this is the case for which the solutions can be uplifted to be solutions of type IIB supergravity in ten dimensions, and are thus particularly relevant from a string-theoretical perspective. From an holographic point of view, the Fayet-Iliopoulos five-dimensional supergravity coupled to n V vector multiplets should be dual to a subsector of an N = 1 SCFT composed by an N = 1 energy-momentum tensor and U(1) n V flavour current multiplets. This description is made more precise in the particular n V = 2 case we are considering, where the supergravity theory can be regarded as a consistent truncation of type IIB supergravity and it is therefore holographically dual to a deformation of N = 4 super-Yang-Mills 1 .
In [28] the conditions, originally given in [13], to be imposed in order to have a supersymmetric solution to Fayet-Iliopoulos gauged supergravity are rearranged and partially solved. This process results into n V + 1 coupled ordinary differential equations. We explicitly examine the n V = 2 case we are interested in and we obtain three coupled ordinary differential equations which we simplify as much as we can using the particular form of the various functions in the U(1) 3 theory. However the equations remain very complicated and we could not find any new analytic solution. We therefore pass to a perturbative approach and we try to construct a near-horizon family of candidate black hole solutions and a near-boundary family of candidate AlAdS 5 solutions. Remarkably enough, we find numerically that the two families such obtained match in the bulk and we are able to construct the whole solution for a large part of the parameter space. We thus obtain a new three-parameter family of supersymmetric black holes, which posseses both non trivial and unconstrained gauge and scalar fields. Our solution generalizes the one-parameter one of [26,27] found in minimal gauged supergravity and also the twoparameters one of [28] since our scalar fields are unconstrained. From a technical point of view, the approach we adopt here to perform the perturbative and numerical analysis is similar in spirit to the one of [28]; however, the results and outcomes we find are more general since we abandon the simplifying ansatz imposed there on the scalar fields, so that we obtain an infinite set of new solutions, enlarging the family found there.
The horizon geometry and the horizon properties of our solution are controlled by two of the three total parameters, which are also responsible for the Page charges and the angular momentum. The last parameter regulates instead the squashing of the boundary geometry, but does not influence the horizon, whose geometry is completely independent on the squashing.
The paper is organized as follows. In sec. 2 we give an essential presentation of Fayet-Iliopoulos gauged supergravity and discuss its main features, including its uplift to string theory in the case n V = 2. We also write the U(1) 3 theory supersymmetry equations we need to solve and present the related first integrals. In sec. 3 we construct our family of solutions, first by analyzing the near-boundary and near-horizon regimes separately and then matching them numerically. In sec. 4 we present the main physical properties of our solution. To evaluate some of them we employ the technique of holographic renormalization. We explicitly compute the entropy and find that, as was also obtained in [28], it is remarkably reproduced by a simple formula containing the Page charges, instead of the holographic charges. We conclude in sec. 5, where we discuss our results. In app. A we report more details about the perturbative analysis of sec. 3, turning also the near-boundary solution in Fefferman-Graham coordinates. Finally in app. B we briefly review the main features of holographic renormalization and how this is used to compute the physical properties described in sec. 4.

N = 2, D = 5 Gauged Supergravity and its uplift
In this section we briefly introduce the N = 2 five-dimensional Fayet-Iliopoulos gauged supergravity and we discuss which conditions should be imposed to its various bosonic fields in order to obtain supersymmetric solutions. For the purposes of the present paper, we are interested only in cohomogeneity one supersymmetric solutions which possess a local SU(2)×U(1)×U(1) symmetry, which implies that the supersymmetry equations are just ordinary differential equations (ODE's) rather than partial differential equations. These supersymmetry equations were originally given in [13] and were recast into a sim-pler form in [28]. We recall that the general case is characterized by an arbitrary number of scalar fields n V while, in the case of interest for the present paper, we will consider n V = 2 since this is the case in which the uplift to type IIB supergravity is possible. The authors of [28] specialize to the particular case where all the components of the scalar fields orthogonal to the scalar vacuum expectation values in the supersymmetric AdS 5 vacuum are controlled by the same H function (up to some constants q I ); here, we will instead specialize to the U(1) 3 theory case with n V = 2 without imposing any simplifying ansatz. The same orthogonal components will thus be treated in full generality and let be free to assume any value; they will be controlled by three different functions H 1 , H 2 and H 3 linked by a unique constraint. Solutions to this particular theory are interesting because they have a natural uplift to type IIB supergravity, as we will discuss in sec. 2.4.

The theory
We now briefly review the main features of the five-dimensional N = 2 Fayet-Iliopoulos gauged supergravity we consider in this paper. Here we give some essential information about the theory and we refer to [13,30,31] for further details.
The bosonic sector of the theory is composed by the metric g µν , by n V + 1 Abelian gauge fields A I µ and by n V real scalar fields Φ I . For convenience it is customary to parametrize the latter using n V + 1 real fields X I which fulfil the constraint where C IJK is a constant, symmetric tensor satisfying The constraint (2.1) can be more easily written by introducing lower-index scalars X I defined as so that now (2.1) becomes We can also obtain the inverse relation by defining the tensor C IJK such that C IJK = δ II δ JJ δ KK C IJK . The Fayet-Iliopoulos gauging procedure introduces n V + 1 parameters V I in order to gauge a U(1) subgroup of the R-symmetry SU (2). In the bosonic sector the main consequence of this gauging is the introduction of the following scalar potential which is fully consistent with supersymmetry. The bosonic action of the theory in mostly-plus signature is: where g = det g µν is the determinant of the metric, F I = dA I are the n V + 1 field strengths and κ 2 = 8 π G N . The scalars appear in the action as contracted with the kinetic matrix Q IJ , which reads: and satisfies From the action (2.7) it is possible to derive the Einstein and Maxwell equations and also the equations of motion for the scalars, which we do not report because we will never need them in the present paper. They can be found in [13]. Assuming that C IJK V I V J V K > 0, there is a supersymmetric AdS 5 vacuum allowed by the theory which is characterized by a radius and by the constant values of the scalars X I . These are furthermore determined by the Fayet-Iliopoulos parameters V I as (2.11) We will always use this relation to trade the parameters V I with theX I in the rest of the paper. Using (2.5) we can therefore rewrite the scalar potential V as We are interested in AlAdS solutions to five-dimensional Fayet-Iliopoulos gauged supergravity. These have a holographic interpretation in terms of a dual four-dimensional N = 1 SCFT, including cases where the latter SCFT is not in a trivial state such as the conformal vacuum. However the holographic interpretation is actually under control in the cases where a consistent uplift of the theory to string theory or M-theory does exist. In section 2.4 we will briefly review how it is possible to embed the five-dimensional N = 2 Fayet-Iliopoulos supegravity we study in this paper as a consistent truncation of type IIB supergravity. We introduce the SU(2) left-invariant one formŝ σ 1 = cosψ dθ + sinψ sin θ dφ , σ 2 = − sinψ dθ + cosψ sin θ dφ , whereψ denotes a coordinate different than ψ to be introduced later. These one-forms satisfy dσ i = − 1 2 ijkσj ∧σ k . We choose the set of coordinates (y, ρ, θ, φ,ψ) to describe our solution and we assume the following ansatz for the five-dimensional metric where all the unknown functions f , w, a are dependent on the ρ coordinate only and for the rest of the paper the prime symbol will denote differentiation with respect to this coordinate. In the minimal theory, the functions f and w can be rewritten in terms of a only; in particular f assumes the form f min = 12 a 2 a 2 (a 2 a − a + 7aa a + 4(a ) 3 ) .
The gauge fields are given by where U I (ρ) are n V + 1 unknown functions to be further determined. The field strengths following from the above gauge fields are: so that their Hodge dual are 2 : Finally the scalar fields are only functions of the radial coordinate ρ, so that X I = X I (ρ).

The supersymmetry equations
Rearranging the supersymmetry equations originally given in [13], the authors of [28] showed that all the supersymmetric solutions of the N = 2, D = 5 Fayet-Iliopolous supegravity under consideration can be obtained by solving the following set of equations: where H I = H I (ρ) are n V + 1 functions which are defined such that and satisfy the constraintX Moreover, we have defined Our aim in this section is to explicitly rewrite the supersymmetry equations (2.19), (2.20) and all the various objects defined in [28] for the case n V = 2. We have the indices I, J, K running from 1 to 3 and We have also that the constraint on the scalars (2.1) becomes 25) and the kinetic matrix (2.8) for the scalars is given by The corresponding valueX I of the scalars in the AdS 5 vacuum is so the corresponding kinetic matrix in the same vacuum is just (2.28) Using these relations we can explicitly write the supersymmetry equations in the context of the U(1) 3 theory we are considering. The equations will depend only on a(ρ) and on three functions so that H 3 will never appear anymore throughout the paper. We define two particular combinations of H 1 and H 2 which will appear in the supersymmetry equations: We now proceed to rewrite the supersymmetry equations. In order to do this we let the index I run from 1 to 3, we use (2.30) to eliminate H 3 whenever it appears and we perform all the necessary contractions recalling (2.24). Doing so we obtain the following three equations (2.34) We have to solve these three equations (2.32), (2.33) and (2.34) in order to find new solutions in the U(1) 3 theory.
Once a, H 1 and H 2 are determined, all the other functions are fixed in terms of these. We now report the explicit expressions of them in the U(1) 3 case. These are straightforwardly obtained by the general relations reported in [28] by setting n V = 2 and performing the necessary contractions. We start with the function f that is determined as where h I are given in terms of H I as They obviously satisfy the constraint h 1 + h 2 + h 3 = 0 as the H I do. We proceed with the w function (2.38) The functions U 1 , U 2 and U 3 are given in terms of H I by and it is immediate to see that they satisfy so that only U 1 and U 2 are independent. Finally we consider the scalars; for the lowerindex scalars X I we have from equation (2.37) that Furthermore, only X 1 and X 2 are independent since it holds the following constraint For the upper-index scalars X I we obtain instead that and we have that only X 1 and X 2 are independent, since it results We now briefly discuss which conditions one should impose to reduce to previously known solutions of the N = 2, D = 5 Fayet-Iliopoulos gauged supergravity and to minimal gauged supergravity. To obtain the latter it is sufficient to take H 1 = H 2 = 0. Indeed, in this case the equations for H 1 and H 2 (2.32), (2.33) are trivially satisfied, while the equation for a (2.34) becomes the same given in [12]. All the physical relevant functions become the ones of minimal gauged supergravity. Indeed equation (2.35) gives f = f min , therefore the scalars (2.47) become just constants, the U 1,2 functions are equal and provide the same gauge field of [12] and equations (2.38), (2.43) reduce to the same form they take in minimal gauged supergravity. We can also easily obtain the U(1) 3 version of the general family of solutions given in [28]. Indeed the solutions of that paper are obtained by taking the simplifying ansatz H I (ρ) = q I H (ρ). Further conditions discussed in [28] fix the charges q I to assume in the U(1) 3 theory the values q 1 = q 2 = 1 6 , q 3 = − 1 3 (or cyclic permutations). Therefore to reduce to this class of solutions we have to take the limit H 1 = H 2 = 1 6 H, H 3 = − 1 3 H (or cyclic permutations). Doing so, equations (2.32) and (2.33) become equal and become the same as equation (2.58) of [28], while equation (2.34) reduces to (2.59) of the same paper. The Gutowski-Reall solution [13] is also recovered, since it is just a particular limit of the more general solutions of [28].

First integrals and conserved charges
The analysis performed in [28] is able to find three first integrals for the black hole solutions constructed there (cf. equations (3.45), (3.46), (3.47) in the above-mentioned paper). Two of them are straightforwardly obtained by considering the component of the Maxwell equation parallel toX I and by the orthogonal ones; while the last first integral is derived by manipulating the component t ψ of the Einstein equations using the Maxwell equation and the supersymmetric conditions. In particular in the special ansatz considered in [28] all the H I functions are equal (up to a constant) and therefore all the components of the Maxwell equation orthogonal toX I are also equal (up to a constant), so they globally provide only one non-trivial first integral. Thus a total number of three first integrals is obtained: two coming from the Maxwell equation and one from the Einstein equations.
In the general case of n V + 1 different H I functions, we find a total number of n V + 3 first integrals: one coming from the component of the Maxwell equation parallel toX I , n V + 1 coming from the orthogonal ones and finally one can be derived from the t ψ component of the Einstein equations. These first integrals are given by where the q I are defined such thatX I q I = 0 and they have to satisfy eq. (2.55) of [28]. Note that in the general n V case, the H I functions satisfy the constraintX I H I = 0, therefore one of the above first integrals is dependent from the others and thus there are n V + 2 independent first integrals. As it is for the three first integrals of [28], also the generalized n V + 3 ones we found have an interpretation in terms of conserved charges. We introduce the conserved Page electric charges as [32] where Σ ∞ denotes the three-dimensional ρ = ∞ hypersurface obtained by the general family Σ ρ which foliates a generic Cauchy surface (a hypersurface of constant time). By using the various definitions and relations given in the subsections above for the quantities appearing in (2.52), it is possible to show that the Page charges can be rewritten as with the overall factor introduced for convenience. The Page charges are therefore described by all the first integrals which derive from the Maxwell equation. The last first integral can be interpreted by considering a conserved angular momentum J referred to the symmetry generated by the killing vector ∂ ∂ ψ . This is given by the following generalized Komar integral [28]: Evaluating this angular momentum on the supergravity background we are considering, we find so that J is proportional to the last first integral. For the present paper, we are interested in the case n V = 2, therefore we should have 5 first integrals in total. By setting n V = 2 in (2.51) and performing the needed contractions, we find the following relations for our U(1) 3 case: Note that the q I are fixed to be q 1 = q 2 = 1 6 , q 3 = − 1 3 , as it is for the U(1) 3 version of the solution of [28]. Moreover, equations (2.53) and (2.55) are still valid with the I index which runs from 1 to 3. In (2.56) we already used the constraint (2.29) to eliminate H 3 in favour of H 1 and H 2 . As a consequence, we immediately see that K are not independent, but they satisfy the relation so that we have 4 independent first integrals in total. We will always use (2.57) to trade K 2 , so the set of independent first integrals we choose is 2 , K 3 ). As we shall see later in the paper, these first integrals will also help us to connect the parameters of the perturbative near-boundary solution we will construct with the parameters of the near-horizon one.
The Page charges P I and the quantity J defined by a Komar integral can be regarded as the electric charges and the angular momentum of the solution. However the procedure of holographic renormalization, which we will employ later in the paper, gives the possibility to define in a different manner analogous conserved quantities playing the same role; we will therefore compare them with the conserved Page charges and angular momentum we have defined in the present section. In particular, since the contribution provided by the Chern-Simons term to the holographic charges is different to the one for the Page charges, we should expect that those conserved quantities are not equal.

Closing remark: uplift to type IIB supergravity
Here we briefly review how we can embed N = 2, D = 5 supergravity with U(1) 3 gauge group in type IIB supergravity following [33]. Starting from type IIB on AdS 5 × S 5 , we can have a consistent truncation turning on the τ = C 0 + ie −Φ , C 4 and g M N fields, where x M = (x µ , y a ) with x µ being the AdS 5 coordinates and y a = (θ,ψ,φ 1 ,φ 2 ,φ 3 ) being the S 5 coordinates 3 .
We can write In this set of coordinates the round S 5 is here 5 ≡ and vol 5 are referred to the AdS 5 metric. Considering the only relevant part of the usual type IIB lagrangian and inserting the ansatz (2.59), with an appropriate field-redefinition, we land to eq. (2.7).

Constructing the solution
The three ODE's (2.32), (2.33), and (2.34) we obtained are difficult to solve analytically, therefore we resort to a numerical method to find new solutions. In order to do so, we series-expand the fields both in the near-boundary region ρ → ∞ and in the near-horizon one ρ → 0, and we fix the series coefficients by solving the ODE's order by order; after that we build our numerical solution by matching the two expansions in the bulk i.e. in a region where they overlap. Our perturbative and numerical approaches are similar in spirit to the ones adopted in [28], however here we will not impose any ansatz on the scalar fields and we will therefore look for more general solutions. In the near-boundary region we find expansions that are compatible with AlAdS solutions, in the near-horizon one we note that there are solutions which possess the characteristics of black holes and, using a numerical procedure, we establish that there are well-behaved solutions interpolating between these two regimes. For ease of notations, we will set = 1 in the whole section and we also change the label of the functions H 1 and H 2 to:

Near-boundary analysis
We now perturbatively solve the equations (2.32), (2.33) and (2.34) around ρ → ∞. This is the limit in which the solution approaches the conformal boundary. The unknown functions are a, H and K. We assume for them the following asymptotic expansions We furthermore assume a 0 = 0. We have included only odd powers of e ρ in the expansion for a: that is because any term weighted by an even power of e ρ would vanish because of the equations. For an analogous reason, the expansions for H and K involve only even powers of e ρ . We have obtained a perturbative solution for the three equations (2.32), (2.33) and (2.34) which is valid up to order O (e −10 ρ ) and is controlled by the following eleven parameters 4 : a 0 , a 2 = a 2,0 , c = a 2,1 , a 4 = a 4,0 , a 6 = a 6,0 , Using the near-boundary solution we found, we can perturbatively evaluate all the other relevant functions. However before computing them, we introduce the following which will be related to the squashing of the three-sphere at the boundary. We will trade the parameter c for v 2 when writing the main results of the present paper, since the latter has a clearer physical interpretation. We furthermore define the change of coordinates which trades y,ψ for t, ψ; in the latter set of coordinates it is easier to see that, at the conformal boundary, the metric is static. In the new set of coordinates, the metric and the gauge fields turn to: the σ i being defined in the same way as theσ i with ψ replacingψ. The functions f and w, which are independent on the change of coordinates (3.6), can be easily evaluated using (2.35) and (2.38). From their explicit expansions (A.3a, A.3b), we see that f goes to 1 in the near-boundary limit, while the w function has a e 2 ρ leading term. Both of these two near-boundary behaviours are fully consistent with an AlAdS solution.
We can proceed to compute the metric (3.7a) in order to verify that this is indeed static. We find that it can be rewritten as ds 2 = dρ 2 + e 2ρ ds 2 bdry + . . . , (3.8) with ds 2 bdry being the metric at the boundary, which reads so it is indeed static, as wanted. Looking at (3.9) we can also see that the last term is the metric of a three-sphere whose squashing is controlled by the parameter v. The change of coordinates affects also the supersymmetric Killing vector V , which becomes: We present the asymptotic form of the scalars X I and the near-boundary expansions of all the components of both the metric and the gauge fields in app. A.1.1, where we express the solution in Fefferman-Graham coordinates. The possibility to write our solution in Fefferman-Graham form will also confirm once more that it is indeed AlAdS 5 and the following analysis will show that four of the eleven free parameters, a 0 , c,H and K determine the bulk fields at the boundary and therefore play the role of a source in the dual quantum field theory. In particular, looking at (3.9) and recalling that v 2 is related to c as in (3.5), it is possible to predict that a 0 and c would control the metric at the conformal boundary, whileH andK should determine the scalar fields. The analysis in Fefferman-Graham coordinates will reveal that this is indeed the case. We refer to app. A.1.1 for further details on the Fefferman-Graham form of the metric.
We also recall that in sec. 2.3 we have introduced the four independent first integrals and K 3 . These are given by equations (2.56). Since the first integrals are obviously constants, we can evaluate them in the near-boundary region using the expansions (A.2a), (A.2b), (A.2c) we have found. In this way, we find that the three first integrals coming from the Maxwell equation, 2 , depend on the various near-boundary free parameters, among which the most subleading are a 4 , H 4 and K 4 . We can use these relations to express the latter free parameters with respect to the others and the first integrals. In the same fashion, we find that the first integral coming from the Einstein equations, K 3 , depends on a 6 ; therefore we also get an expression for a 6 with respect the other near-boundary free parameters and the first integrals. These relations for a 4 , H 4 , K 4 and a 6 are rather involved and are thus reported in app. A.3, in particular they are given by equations (A.29). The advantage we get from these equations is the following: we are able to evaluate all the first integrals both in the near-boundary region as well as in the near-horizon region, obtaining them as functions of, respectively, the near-boundary and the near-horizon parameters; combining the relations obtained in the near-horizon with the ones obtained in the near-boundary we will be able to express a 4 , H 4 , K 4 , a 6 as functions of only the remaining near-boundary parameters and the near-horizon ones. This will allow us to replace the most subleading parameters a 4 , H 4 , K 4 , a 6 with the others.

Near-horizon analysis
Having shown that solutions compatible with an AlAdS 5 behaviour exist in the nearboundary, we proceed to solve the ODEs (2.32), (2.33) and (2.34) in the region near to ρ → 0, which we identify with the interior region of our solution. It is reasonable to assume that the three unknown functions K, H and a can be written as a Taylor expansion near ρ = 0 We want to search for either new black hole solutions or new soliton solutions; in the first case we should have an event horizon at ρ = 0, in the second one the solution should close off smoothly in the same point. Looking at the metric (2.14), we see that both the types of solutions require α 0 = 0, which we therefore assume. Moreover, due to the symmetries of the ODE's we can take α 1 > 0 with no loss of generality 5 . We proceed to solve (2.32), (2.33) and (2.34) perturbatively up to O(ρ 13 ). We find that equations (2.32) and (2.33) fix uniquely the form of K and H with the coefficients ι 0 , η 0 , ι 1 , η 1 forced to vanish and all the others determined by the free parameters ι 2 , η 2 and by the coefficients of a. We find it convenient to define the new parameters so that the first terms in the expansions of K and H are: (3.14) Note that switching the parameters η ↔ ι we have that H ↔ K, as expected. The above expansions solve (2.32) and (2.33) without imposing any conditions on a. This function will then be constrained by the remaining equation (2.34) on which we now focus. The solution process of the latter equation brings us to distinguish between different cases. To solve the first non-trivial order of (2.34) we must satisfy the condition: which means that either α 2 = 0 or the parenthesis vanishes. An analogous condition was found in [28], while in the minimal gauged supergravity (which can be obtained by setting ι = η = 0) the corresponding condition forces to choose α 2 = 0 since the parenthesis cannot vanish in this theory [34]. Since we are interested in solutions which have a minimal gauged supergravity limit, we choose α 2 = 0 and proceed further. The next order yields: α 4 5819 α 6 − 5244 α 4 + 12 α 2 6912 η 2 + η ι + ι 2 + 65 + 32(36 η + 1)(36 ι + 1)(36 η + 36 ι − 1) = 0 , (3.16) this equation can be satisfied if α 4 = 0 or the parenthesis vanishes. In [28] setting the corresponding parenthesis to zero led to the black hole solution studied there, while setting α 4 to zero the only regular solution obtained was the one of [13]. In the minimal theory it is possible to obtain the black hole of [27] by setting α = 8 11 while the choice α 4 = 0 leads either to the regular soliton of [34] or to the black hole of [12]. The soliton found in [34] is a solution of minimal gauged supergravity of special kind: indeed, considering the expression of f min in (2.35), we can write f as in order to describe a soliton, the f function must start with a constant term in a small ρ-expansion, so that the solution closes smoothly. However from (3.17) we can argue that this is possible only if the numerator and the denominator have the same leading behaviour at small ρ. Plugging the expansions (3.11) in (3.17), we can easily check whether this is possible or not; in particular we note that the numerator goes as ρ 3 at small ρ, so we have to impose the same behaviour to the denominator. In the minimal case H = K = 0 this is easily achieved by taking α = ± 1 2 , since P = a a 3 + a a 7 a a + 4 (a ) indeed starts with a ρ 3 term if and only if α = 1 2 . This choice for α is the one taken in [34] and leads the author to find a soliton solution. In the general case we are considering in this paper, H and K are non vanishing, therefore recalling (3.11) it is evident that the denominator of (3.17) goes always as ρ, while the numerator begins with ρ 3 . We then conclude that there is no possibility to find a soliton solution in the U(1) 3 theory with non trivial scalars. We have furthermore verified that, even in our general framework where we do not have imposed any ansatz on the scalar fields, the choice α 4 = 0 leads only to the solution of [13] or to a singular solution 6 .
We therefore proceed to analyse the only new solution we find, which is obtained by setting the parenthesis in (3.16) to zero. This gives (3. 19) Here one has the possibility to choose either the plus or the minus sign; we leave this choice unspecified for now and proceed further. Setting ι as in (3.19), we continue to perturbatively solve the equation (2.34) finding that the solution is uniquely determined in terms of the free parameters α, η, α 3 and α 4 . We report in app. A.2 the first terms of the expansions of a, H and K. They are given by equation (A.25). In the same appendix, we also give more information about the near-horizon solution. Note that looking at (A. 25) it can appear that K can be obtained from H switching η ↔ ι, however one must keep in mind that ι is not a free parameter anymore, being fixed as in (3.19).
We now briefly show how to reduce our general solution to the U(1) 3 version of the one constructed in [28]. As discussed at the end of sec. 2.2, we need to impose H = K, which means η = ι. The condition (3.19) then becomes an equation for η which gives: which is consistent with the fact that it must be H limit = 1 6 H there as already stated above. As consequence, the expansions (A.25) fully reproduce the ones of [28].
The perturbative solution near to ρ = 0 is characterized by four free parameters: α, α 3 , α 4 and η. However due to the scale symmetry of the supersymmetry equations discussed at the end of sec. 2.2, we can rescale one of the parameters without changing the solution we found. In particular we note that α and η are left invariant under the action of these symmetries, while α 3 and α 4 can be rescaled. We can therefore argue that only three of the free parameters we found are physical and we choose to consider α 3 as an unphysical parameter. We will explicitly use the possibility to rescale α 3 to numerically match our small-ρ behaviour with the near-boundary one discussed in sec. 3.1, showing that an interpolating solution indeed exists for different values of the remaining physical free parameters.
We also define a new parameter ξ such that it is invariant under the scaling symmetry discussed at the end of sec. 2.2: We shall use this definition to trade α 4 with ξ wherever the former appears. From now on, our set of independent near-horizon free parameters will then be (α, η, ξ).
We now proceed to report the explicit near-horizon form of the metric, the gauge fields and the scalars, showing that they are compatible with a black hole which is a generalization of the U(1) 3 version of the solution presented in [28]. Here we present only the most relevant information about the near-horizon analysis and we refer to app. A.2 for further details. Although it is fixed as in (3.19), we will keep ι for compactness in the formulae below. We begin by presenting the form the five-dimensional metric takes at leading order in the small ρ expansion: where we have defined as ∆ and Θ the two following quantities (3.23b) We proceed to report the near-horizon expansions of the scalar fields X 1 and X 2 : while the last scalar field X 3 is easily determined using the constraint between them, given by eq (2.50). We do not report the expansions for the scalars with lower indices, X I , since they can be straightforwardly obtained from (2.3).
Finally we show the near-horizon behaviour of the gauge fields, which is

25b)
and, again, the third gauge field A 3 can be easily determined by the other two and will not be presented here.
Looking at the expansions of the various supergravity fields reported above, it is evident that the pertubative solution can be regarded as the near-horizon expansion of a black hole whose horizon is located at ρ = 0. Indeed, the metric has a divergent radial component which is O(ρ −2 ) while its spatial part stays finite as the limit ρ → 0 is approached. Furthermore, the supersymmetric Killing vector V , given by eq. (3.10), is everywhere timelike but on the horizon, where its norm −f 2 vanishes 7 . All the scalar fields stay regular as the limit ρ → 0 towards the horizon is approached and the same do the gauge fields, which are furthermore transverse to V in the gauge we have chosen 8 .
We now have to understand for which choice of parameters the solution has a well defined horizon at ρ = 0. In order to ensure regularity of the horizon we need that all the spatial diagonal metric components in eq. (3.22) retain their sign for every value of the radial coordinate ρ. This means we must assure that g ii > 0 for i = ρ, θ, φ, ψ 9 . From (3.22) it is easy to see that this translates into imposing the conditions ∆ > 0 and Θ > 0. Indeed, this ensures the positivity of g ii for every value of ρ. We are still left with the possible sign choice in eq. (3.19); both choices give a well defined black hole solution and we will analyze the parameter space for both of them, even if in the following we will report the numerical results only for the minus sign choice, which is the choice that leads to the largest space of regular solutions. In fig. 1 we report the parameter space in terms of α and , with the latter defined via so that the limit to H = K case of [28] is simply reproduced by the choice = 0. Note that we can trade η with using (3.26) only if α ≤ 8 11 , that is the maximum value considered for α in fig. 1; we have analyzed the parameter space for α > 8 11 and for generic values of η finding that no regular black hole horizon with real coefficients appears in this region 10 .
In fig. 1 we have reported the parameter space for both the sign choices. The region colored in red is the region where g ψψ > 0, while the region in blue is where g θθ > 0. We have colored in yellow the regions where we managed to find numerically a regular black hole solution with real coefficients, by interpolating the near-horizon solution of the present section with the near-boundary one of the previous section. In particular, the 7 This can be easily seen by eq. (A.26) in app. A.2, where the near-horizon expansion of f is presented. From this expansion it is evident that f ∼ ρ 2 at ρ ∼ 0 and therefore it vanishes at the horizon. 8 Indeed, it results V µ A I µ = 0 as it is easy to verify using eq. (3.25a). 9 Actually to ensure regularity of the horizon we should also require that g yy ≤ 0, where the equal sign holds only at the horizon ρ = 0. However this is already guaranteed by the fact that g yy = −f 2 with f being, as already stated, a real function which vanishes at the horizon. 10 Note that the regularity conditions ∆ > 0 and Θ > 0 must be combined with the existence condition of the square root in (3.19). We found that these three requests are never simultaneosly verified when α >  yellow dots are the point characterized by the most extreme values for the parameters α , for which we found a regular numerical solution. The points in the purple region that are not in the yellow one represent values of α , for which the horizon is well behaved but a full solutions seems not to exist. This is because we find divergences in the bulk when we try to numerically interpolate the near-horizon region with the near-boundary one. Notice that we have reproduced the results of [28] on the axis = 0. We have a nice explanation for the peculiar behaviour appearing at the point (α , ) = ( 2 3 , 0), reported in fig. 1a with a green dot, where in [28] it emerges a non-analytic behaviour of ∆: this is due to the peculiar structure of the ∆ function in the (α , ) plane.
Notice that, as it can be easily seen in fig. 1b, the region of existence of regular black hole solutions with the plus choice in eq. (3.19) is smaller than the one obtained with the minus choice; this is clearly visible from the form of the "yellow triangle" of solutions in the two cases. This is what is also found in the case = 0 of [28]. We also stress the fact that both g θθ and g ψψ quickly drop to be negative outside the region of the parameter space we have shown in the figure, so no regular horizon can be found there. There is a possible exception only for α ∈ (0, 1 2 ), where instead we have found a region of regular positive g θθ and g ψψ , but there H and K becomes complex.
We conclude this section by observing that we can use the near-horizon solution to obtain the dependence of the first integrals (2.56) on the near-horizon parameters α and η 11 . In order to do this, we have just to plug our near-horizon expansions for the supergravity functions into (2.56) and perform the computations. For conciseness we do not report the expressions such obtained here, but they are given by eq. (A.30) in app. A.3, together with the relations between the near-boundary parameters and the K I . As we mentioned at the end of sec. 3.1, confronting the near-boundary and the near-horizon expressions for the first integrals, we are able to write the most subleading near-boundary parameters a 4 , H 4 , K 4 and a 6 in terms of the remaining near-boundary ones and the near horizon ones. This allows us to eliminate these four parameters in all the expressions and we will proceed by doing it throughout the paper, as it simplifies many expressions.

The matching solution
In this section we proceed to match the near-horizon perturbative solution constructed in sec. 3.2 with the near-boundary one obtained in sec. 3.1, showing that a smooth interpolating solution indeed exist. This happens for all the points in the parameters space of fig. 1 that are inside the yellow region.
We begin by giving a brief explanation on how we construct the numerical solution. We use the near-horizon expansions of the previous section to set the initial conditions in the vicinity of the horizon, located at ρ 0, and then we numerically integrate the supersymmetry equations (2.32), (2.33) and (2.34) towards the near-boundary region, i.e. towards large values of ρ. We recall that we found in the near-horizon an unphysical parameter, α 3 , which may be rescaled at will; we use this possibility to set the appropriate rescaling such that the AlAdS behaviour of a holds in the near-boundary region. Obviously, in order to integrate the equations, we need to give numerical values to the near-horizon parameters α, and ξ. We tried many different values for the parameters α and in the whole possible region of regularity of the solution (which coincides with the region colored in purple in fig. 1) finding regularity in the interior only in the points (α, ) inside the yellow region. This means that for every point in the yellow region there is an interval of allowed values of ξ for which all the components of the metric, the scalars and the gauge fields are regular. The allowed interval of ξ depends on α and and is determined by regularity of the boundary geometry. All the points outside the yellow region lead to solutions which present fields that are not regular in the bulk; in particular for such solutions the function f turns out to have always a divergence at finite ρ. We shall therefore discard such solutions.
We found that the region of regularity of the solution corresponding to the minus sign quite subleading, it is instead absent.

choice in (3.19) is inside
while a similar, but smaller, range is found for the plus sign solution. From now on we will specialize on the minus sign choice, but all the characteristics of the solutions we will discuss are present also in the ones obtained choosing the plus sign. We constructed the full interpolating solution for many values of the near-horizon parameters inside the bounds reported in (3.27). As illustrative examples, we discuss in the following two different analyses performed on the solution: the first is made by fixing α and ξ and studying solutions with different , the second one is made by fixing α and and studying solutions for various values of ξ. The first analysis gives us the possibility to compare the characteristics of the new solutions we have found with the ones of [28], which are obtained by setting = 0. Since, as we know from [27,28,34], the parameter ξ is related to the squashing at the boundary, the second analysis allows us to study the new solutions (which present = 0) with different squashed boundary geometries.
We begin by presenting the solutions with different . We choose α = 0.84 and ξ = − 1 4 . In fig. 2 we show the numerical behaviour of the metric components and of the basic functions a, H, K. It is easy to see that in the near-horizon region their behaviour is in general different for the various choices of ; an exception are the f and a functions for which the differences are very small. We should notice that g θθ and g ψψ tend to a positive non-zero value for ρ = 0 and are always positive; this means that our solution has no Closed Timelike Curves (CTCs) in the whole region ρ ≥ 0. We have verified that the same happens for many different values of the parameters in the yellow region of regularity of fig. 1a. Also, since our solutions are rotating solutions, it is clear from the plot of g tt that an ergoregion emerges 12 . In fig. 2b we display both H and K, opportunely rescaled with a ρ −1 e −2ρ prefactor, to show that they are indeed different for all the choices of but the case = 0, where we have H = K.
In fig. 3 we show the numerical solutions for the scalar and gauge fields. From this picture is quite evident how the change in affects the global structure of the solution, since in the near-horizon region the fields get attracted to different asymptotic values, while, as for the metric components, they are all attracted to the same large ρ value.
We now proceed to study a particular solution for some different values of ξ. We choose for the other parameters the values α = 0.84 and = 0.008. 12 More precisely, the fact that g tt becomes positive implies that the vector ∂ ∂t becomes spacelike. If we regard this vector as the generator of time translations, then an ergoregion does exist in our solution. However we may also choose the vector (3.10) to be the generator and in this case there is no ergoregion. Similar features appear often in supersymmetric AdS black holes and were first discussed in [12].
The solution H, K.  and different values of , reported in the label. Each function is rescaled by its asymptotic behaviour at large ρ. We emphasize that both g θθ and g ψψ are positive in all the ρ ≥ 0 region, so our solution does not have any CTCs. Since instead g tt assumes positive values near the horizon, our solution does have an ergoregion.
We reported in fig. 4 the relevant functions a, H, K and the metric components. As opposed to the fixed ξ case, here the components go for ρ → ∞ to different values. Furthermore, their behaviour is very similar in the near-horizon region. This is because the effect of having a different ξ is almost negligible in the near-horizon, since the horizon geometry is controlled by α and , while in the near-boundary region the same effect is relevant, being the ξ parameter related to the squashing at the boundary. Again, we have an ergoregion, where g tt becomes positive, and no CTCs, since both g θθ and g ψψ are positive everywhere.   We then show fig. 5 where we reported the gauge and scalar fields for various ξ; again, we see that, in contrast with the fixed ξ case, their behaviour in the near-horizon region is similar for all the ξ, while they go to different values in the near-boundary region. The only exception is the value of A I t , that also differs in the near-horizon region. This is due to the fact that in the coordinates (t, ψ) we are using the time component of the gauge fields explicitly depends on the squashing v, as it is clearly visible by (3.25a). If we had used instead the coordinates (y,ψ), the time component A I t would vanish at the horizon (a) The solution a.     and different values of ξ, reported in the label. Each function is rescaled by its asymptotic behaviour at large ρ. We emphasize that both g θθ and g ψψ are positive in all the ρ ≥ 0 region, so our solution does not have any CTCs. Since instead g tt assumes positive values near the horizon, our solution does have an ergoregion.  and would not be influenced by ξ. We end this section by summarizing the main characteristics of the family of solutions we have constructed. Both the near-horizon analysis and the numerical one prove that our solutions are black hole solutions whose horizon geometry is controlled by two of the three near-horizon parameters, α and η. The last near-horizon parameter, ξ, is related to the squashing at the boundary, and is therefore related with the parameter v 2 controlling the squashing of the boundary three-sphere. Both the near-boundary analysis and the numerical one show that our solutions are AlAdS, a conformally flat boundary being obtained only when the S 3 is round (v 2 = 1). In the near-boundary region, the solution is controlled by eleven free parameters (a 0 , a 2 , a 4 , a 6 , v 2 , H 2 , H 4 ,H, K 2 , K 4 ,K), these should be connected with the near-horizon ones and should be dependent on them. We were able to trade four of them (a 4 , a 6 , H 4 , K 4 ) with the four linearly-independent first integrals (K 1 , K 2 , K 3 ) we have found in sec. 2.3; these first integrals are useful since they are immediately connected with the interior of the solution, being possible to express them with respect to the near-horizon parameters only. All the other near-boundary free parameters can be related to the near-horizon ones using a numerical procedure, as it was done in [27,28,34]. We have numerically shown that the near-boundary and nearhorizon behaviours we have found interpolate smoothly in the bulk, giving rise to regular solutions which are free from CTCs.

Physical properties of the solution
In this section we compute the relevant physical quantities that characterize the family of solutions we have built. These are the energy, the angular momentum, the holographic and Page charges, which can be computed using the near-boundary perturbative solution, the chemical potentials and the entropy, which instead can be derived by means of the near-horizon expansions. Once these quantities are known, we can perform some consistency checks, for example by verifying the quantum statistical relation.
In order to compute some of the above physical properties, we will use the technology of holographic renormalization [11,[35][36][37][38][39]. We perform such computations using the Fefferman-Graham radial coordiante r, introduced in app. A.1.1, instead of the usual one ρ. This is because the use of the Fefferman-Graham coordinate is standard in holography and may help to compare our results with other references. Moreover, we write the general Fefferman-Graham metric as where the five-dimensional coordinates split as x µ = (r, x i ) with x i = {t, θ, φ, ψ} and where h ij is the induced metric at the boundary of the spacetime. Similarly we define the boundary gauge fields A I i and the boundary field strengths F I ij . In this section we only report the results we got using holographic renormalization while we refer to app. B for a more detailed discussions about how these results are obtained. We remark that we performed holographic renormalization using a minimal subtraction scheme; all the physical quantities evaluated by means of this formalism are therefore refereed to this renormalization scheme. We also underline that the application of holographic renor-malization to five-dimensional Fayet-Iliopoulos gauged supergravity has been discussed in detail in [28].
We can compute via holographic renormalization the stress-energy tensor of our family of solutions 13 . Its expression is quite standard where R ij , R, B ij are the Ricci tensor, the Ricci scalar and the Bach tensor of the induced metric h ij . The other ingredients appearing in the expression (4.2) for the stress-energy tensor are the extrinsic curvature K ij of the induced metric h ij , its trace K and the superpotential W = 3 −1X I X I which derives from the scalar potential V given in (2.6). The conserved electric current also arises from holographic renormalization and it is given by where u i ∂ i = v 2 a 0 ∂ t is a unit timelike vector for the metric on the conformal boundary. The conserved charges (4.4) explicitly are (4.5) Note that the angular momentum precisely coincides with the generalized Komar integral (2.55). We can also compute the conserved electric charges as and it is fundamental for our discussion to remark that these differ on our solutions from the Page charges In fact, on our solutions, we have 2 ) , 2 ) , (4.9) We want to stress that this fact is a consequence of the squashing of the boundary: it is trivial to see from eqs. (4.6, 4.7) that the difference between holographic and Page charges is related to the Chern-Simons term which gives a different contribution to the two quantities. In usual non-squashed solutions like [12,13] the same term gives no contribution, since the field strength F I vanishes asymptotically, and the two different types of charges are therefore equal. This implies that we may have some relevant departure from the equation of [40] that relates the entropy and the charges, since we have no unique way to choose which charge is the correct one for reproducing the entropy; in fact, it will turn out later that the entropy is indeed reproduced in terms of the Page charges, instead of the holographic charges obtained in (4.8).
We remark that our results for the holographic and Pages charges, the energy E and the angular momentum J are valid for every AlAdS solution to five-dimensional Fayet-Iliopoulos gauged supergravity with n V = 2 which satisfies the supersymmetry equations (2.32, 2.33, 2.34). In particular all the conserved charges depend only on the first integrals K 1 , K 2 , K 3 , on the squashing at the boundary v and on the two scalar sourcesH,K. As we have already discussed in sec. 3.2 and as we will show explicitly in app. A.3, the first integrals are completely determined by the interior of the solution, since they can be expressed with respect the near-horizon parameters α and η (see (A.30) for their explicit expressions).
We now turn to examine the near-horizon properties of our solution. We can easily compute the entropy by looking at the horizon metric (3.22): This is the anticipated result: the entropy of the black hole solutions with squashed boundary can be reproduced by a simple combination of the conserved charges; however the formula does not involve the usual holographic electric charges but rather the Page charges, signaling a relevance for the entropy counting of the Chern-Simons term that was previously unnoticed by the non-squashed solutions of [12,13]. This means that the entropy formula found for example in [40] for asymptotically AdS 5 black holes retains its validity for these AlAdS 5 solutions if we identify the charges appearing therein with the Page ones instead of the holographic ones. The result we have found was already obtained in [28] and is somehow anticipated by the explicit form the different types of charges take: indeed the Q I depend on both the squashing v and the scalar sourcesH,K on which the horizon geometry is independent and therefore appear to be inadequate to describe an horizon quantity like the entropy. The Page charges P I , instead, depend only on the first integrals which are immediately related to the horizon geometry only and seem thus the right charges to describe the entropy. The angular velocity Ω of our family of solutions is easily read from the supersymmetric Killing vector (3.10) and it results: the electric potential is instead found to be 14 : Our next aim is to use all the quantities we have evaluated in this section to verify the extremal limit of the quantum statistical relation valid for general AlAdS spacetimes. The quantum statistical relation is given by where I is the Euclidean on-shell action, T is the temperature of the black hole and β = 1 T . We can take the extremal limit of this relation by recalling that for extremal black holes we have T = 0. Sending the temperature to zero, we obtain at the leading order the relation: which is valid for extremal black holes and should therefore be satisfied by our family of solutions.
To verify eq. (4.14) we are missing the Euclidean on-shell action. This can be computed using again holographic renormalization. In particular it is possible to define a renormalized Lorentzian on-shell action as where the regularized action is In the equation above, S bulk is the bulk action (2.7), while S GH and S ct are the Gibbons-Hawking term and the counterterms piece respectively. We show in app. B the explicit form of those pieces as well as how the whole computation of the renormalized action is performed; here we report instead only the final result, which is This action is evaluated in a minimal subtraction scheme, as all the other quantities computed via holographic renormalization reported in this section. We immediately notice that this depends only on the squashing at the boundary v and on the scalar sources H,K. Moreover, we should remark that the regularized action is gauge-dependent due to the Chern-Simons term in the bulk action (2.7); the result we reported above is valid when the gauge condition V µ A I µ = 0 is imposed at the horizon. This particular gauge is justified since it ensures regularity of the solution by avoiding divergences in the square norm of the gauge fields.
The Euclidean continuation of the Lorentzian regularized action can now be obtained by performing a Wick rotation on the time and by making the latter periodic of period β, so that we have Now we finally computed the last quantity involved in the quantum statistical relation. Plugging all the ingredients into eq. (4.14), we see that it is indeed verified. Note also that by recalling that Φ I = 0 for our family of solutions and by defining the holographic charge associated to the Killing vector (3.10) as: the quantum statistical relation assumes the form: which can also be seen as the BPS relation between the holographic charges, the anomalous contribution of [42,43] being already included.

Conclusions
In the present paper we have constructed a new family of supersymmetric AlAdS 5 black holes with a boundary geometry containing a squashed S 3 . These black holes generalize the solutions previously found in minimal gauged supergravity [26,27] and also the one of [28] for n V = 2 (since we have not imposed any ansatz on the scalar fields), and can be uplifted to be solutions of ten-dimensional type IIB supergravity. Our family of solutions depends on three-parameters of which two regulate the horizon geometry, the angular momentum and the Page charges while the remaining one determines the squashing at the boundary. The horizon properties are totally independent on the squashing; therefore if we set a particular horizon geometry by fixing the two former parameters, whatever the squashing the S 3 metric will flow to a fixed one at the horizon. This has somehow the flavour of the attractor mechanism for scalar fields in four dimensions [44]. Let us compare the number of parameters describing the horizon geometry we obtained with the number one should expect by a theoretical counting. For the n V = 2 model we are examining, we have five conserved charges, which are the energy, one angular momentum and three electric charges; however only four of these five total charges are independent since supersymmetry imposes one linear constraint among them. One can therefore expect to find black hole solutions with four parameters regulating the horizon geometry, but already in the solution of [13] one of these is constrained by the requirements to be imposed to avoid causal pathologies [45], so the independent parameters are three. We should therefore expect to be possible to find squashing solutions with three independent parameters regulating the horizon geometry, in addition to the one determining the squashing at the boundary. Our family of solutions presents an horizon geometry described by two parameters, generalizing the solutions constructed in [28] characterized by only one parameter regulating the horizon geometry due to the ansatz for the scalar fields adopted there. The black holes presented here are thus the most general squashed solutions found in the U(1) 3 theory we are studying; however we are still missing squashed solutions with an horizon geometry regulated by three parameters, which should be the most possible general ones according to the theoretical counting arguments reported above. It could be that the general three parameter solution breaks the SU(2) × U(1) 4 symmetry in the bulk and should thus be searched in a more general setup than Fayet-Iliopoulos gauged supergravity.
We have seen in sec. 4 that the entropy of our family of black holes is reproduced using the Page charges instead of the holographic charges. The two types of charges are different for AlAdS 5 solutions due to the presence of the Chern-Simons term, which does not vanish asymptotically like in the case of non-squashed solutions. The formula thus obtained for the entropy in terms of the conserved charges is in agreement with the typical one for asymptotically AdS 5 black holes reported for example in [40], provided the fact that for AlAdS 5 solutions the charges appearing there must be identified with the Page charges and not with the holographic ones. According to the extremization principle proposed in [25], the entropy of supersymmetric asymptotically AdS 5 black holes can be obtained by Legendre-transforming a certain function of chemical potentials conjugated to the black hole conserved charges. This has been further discussed in [22], where the authors identified the function to be Legendre-transformed as the on-shell action of the black hole and managed to reproduce the entropy of [15] using the extremization principle. The result for the entropy we have obtained in this paper suggests that the same extremization principle would work for our family of black holes if one takes into account the Page charges and their conjugate chemical potentials instead of the electric holographic charges. This distinction cannot be established looking at asymptotically globally AdS 5 solutions since the holographic charges and the Page charges coincide for them, due to the fact that the Chern-Simons term vanishes at the boundary.
AB wishes to thank the Institut de Physique Théorique, Université Paris Saclay, CEA, CNRS for the hospitality while the final part of this project was realized.
A More on the perturbative solution A.1 More on the near-boundary solution In this appendix we provide further details on the near-boundary analysis of sec. 3.1. We give some information and results about the main functions characterizing the solution in this region and we write it in Fefferman-Graham coordinates. These are the well-suited coordinates to be used to describe the near boundary behaviour of our solution, since they both allow to check that it is indeed AlAdS and to provide an holographic interpretation, helping us to understand the role of the different parameters controlling the expansions of the supegravity fields from a field theoretic point of view. All the functions presented in this appendix are evaluated in the coordinates (t, ψ) defined in (3.6) as We begin by reporting the first terms of the perturbative solutions for a, H and K: a(ρ) = a 0 e ρ + (a 2 + c ρ) e −ρ a 0 + a 4 + (2 − 16 a 2 − 5 c) c 12 + 54 (2 H 2 + 3H + K 2 )H+ From here we can see that, already at the displayed order, the coefficients of both H and K enter in the solution for a; this case is thus different from both the minimal gauged supergravity one [34] and the solution of [28], where K = H. Furthermore we have verified that the expansions reduce to the ones of [28] Looking at the expansions above we note that when ρ → ∞ we have the limits f → 1 and w → e 2ρ : both these behaviours are consistent with an AlAdS solution.
In the coordinates (A.1), the metric and the gauge fields become ds 2 = g ρρ dρ 2 + g θθ (σ 2 1 + σ 2 2 ) + g ψψ σ 2 3 + g tt dt 2 + 2g tψ σ 3 dt , (A.4a) The various components of the fields assume the following expressions Instead of giving the explicit near-boundary expansions for the various components of the metric, the gauge fields and the scalars in the radial coordinate ρ we used until now, we prefer to switch to the Fefferman-Graham radial coordinate r, since it provides a clearer and easier holographic interpretation of the solution.

A.1.1 Near-boundary solution in Fefferman-Graham coordinates
We begin by presenting the general expected forms the fields of our solution should take in Fefferman-Graham coordinates. The metric should assume in these coordinates the following form: having denoted the radial coordinate as r, while x i describe the various hypersurfaces at fixed r. Each of these hypersufaces has an induced metric h ij that presents the following large r expansion ij +h (4) ij log r 2 +h Obviously all the coefficients in the expansion above depend only on the x i coordinates and are thus independent on r. In a similar way, one can build in the radial gauge A I r = 0 the following large r expansion for the gauge fields For the large r ansatz we shall assume for the scalar fields, the situation is more involved, since this varies according to the mass of the scalar fields and to the conformal dimension of the corresponding dual SCFT operators. The scalar fields X I involved in the present paper have a mass that fulfils the relation m 2 2 = −4, so the corresponding SCFT operators have dimension ∆ = 2. From [38] we can therefore read that such scalars should present a Fefferman-Graham expansion of the form To move in Fefferman-Graham coordinates, we need to switch our radial coordinate ρ to the Fefferman-Graham radial coordinate r so that the g ρρ component of the metric becomes dr 2 r 2 . This means that we have to impose that which is equivalent to solve the ODE Finding an analytical solution of this equation turns out to be very hard, therefore we solve it perturbatively at large ρ obtaining a 2 0 r 2 = a 2 0 e 2ρ + This change of coordinates is such that with respect to (t, r, θ, φ, ψ) the metric becomes with all the components h ij having a large r expansion which is consistent with (A.7). The various coefficients of the metric component expansions assume the following expressions: .

(A.17)
Comparing with App. A of [34] we may notice that both the leading order terms and the next-to-leading order ones are the same as in minimal gauged supergravity. The effect of the supergravity vector multiplet fields is instead manifest at the following order where coefficients dependent onH, H 2 , H 4 ,K, K 2 and K 4 appear. In [28] the backreaction of the fields appeared at the same order, but it was dependent only on the three parameters of the function H introduced there, playing a role analogous toH, H 2 and H 4 . As it is well-known, the free parameters of the metric should appear in the coefficients h (0) and h (4) of the Fefferman-Graham expansion. In particular the parameters controlling h (0) should play the role of the source for the energy-momentum tensor T µν of the dual CFT, while h (4) should be its expectation value. The analysis of [34] shows that we should expect five free parameters appearing in the metric and these are found to be a 0 , c, a 2 , a 4 and a 6 , as it is possible to check looking at the various h (0) and h (4) coefficients above.
Let us now proceed to examine the scalar fields. They are indeed found to take the form (A.9) and their Fefferman-Graham coefficients are  We do not report the expansion for the third scalar X 3 since it is fixed by the constraint (2.50). We note furthermore that the expansion for the scalar X 2 coincides with the one for X 1 upon performing the switching which is equivalent to switch the coefficients of H with the ones of K and viceversa. Both X 1 and X 2 should have two free parameters in the coefficients (0)φ1,2 and (0) φ 1,2 which are interpreted respectively as the sources and the expectation values of the dual scalar operators. Looking at the expression of the coefficients above we clearly see that the free parameters areH,K, which are associated with the sources, and H 2 , K 2 , which are associated with the expectation values. We do not report the expansions of the lower-index scalars X I since they can be easily obtained using their definition (2.3).
Finally we look at the gauge fields. We find that they can be expressed in the Fefferman-Graham form (A.8) and their coefficients are: As for the scalars, we have not reported the third gauge field A 3 , since it is not independent from A 1 and A 2 and can be therefore computed by them. Furthermore we also note that the gauge field A 2 can be obtained from A 1 switching the parameters as in (A.20). Supersymmetry constraints the form of the gauge fields, relating them to the metric and the scalar fields. Indeed, looking at the expressions of the coefficients above, we see that the only free parameters appearing in the gauge fields are H 4 and K 4 ; all the other parameters were found before either in the metric or in the scalar fields. Furthermore the equations of motion in Fefferman-Graham coordinates do not determine the coefficients (0) A 1,2 and (2) A 1,2 ; these are rather fixed by supersymmetry.

A.2 More on the near-horizon solution
In this appendix we report additional information about the near-horizon solution we presented in sec. 3.2. In particular, here we show the small ρ expansions of all the functions not presented in the main text. We begin by reporting the first terms of the expansions of the functions a, H and K, which determine the whole near-horizon solution: Recall that ι is not an independent parameter, but it is fixed as in (3.19). Although at the displayed order the near-horizon expansion of the a function is equal to the one of minimal gauged supergravity [34], the presence of non-trivial scalars becomes evident at next orders where terms dependent on η appear. Using these expansions we can compute the functions f and w via eqs. (2.35) and (2.38) respectively. We obtain the following where ∆ is defined as in (3.23a). Looking at (A.26), we note that f begins as ρ 2 , so that it vanishes as ρ → 0. As explained in the main text, this is totally compatible with a black hole solution whose horizon is situated at ρ = 0. Note furthermore that terms controlled by η and ι appear already at the leading order in both the functions; therefore it is evident that these expansions differ from the ones of [28] and from the minimal gauged supergravity ones [27,34] providing a generalization of both of them. The last functions we explicitly report are the U I , which appear in the gauge fields. The functions U 1 and U 2 are easily computed by (2.39), (2.40) and we find while U 3 can be obtained by these using the constraint (2.42).

A.3 First integrals and equations for the subleading parameters
Here we report the explicit expressions for the near-boundary parameters a 4 , H 4 , K 4 , a 6 one can find by evaluating the first integrals K 1 , K 2 , K 3 in the near-boundary region and the form the latter assume in the near-horizon one.
By evaluating (2.56) in the near-boundary and solving the resulting equations for a 4 , H 4 , K 4 , a 6 we find the following relations 2 ) +H 6 4056K 2K + 8K 2 (153K 2 + 7) + 3687K 2 + 174(K  As a check of their correctness, we can evaluate these relations in the limit leading to the family of solutions of [28] and compare the results with the expressions for the same parameters reported in that paper. As explained in the main text, in order to take this limit we have to set H = K and to rescale the two functions by 1 6 . This is equivalent to set H 2 = K 2 = H limit 2 6 and to impose the same to H 4 , K 4 ,H,K. We have verified that, performing this limit, all the above expressions for the near-boundary free parameters reduce to the ones reported in [28]. Note that equations (A.29) are valid for every AlAdS solution with the characteristics discussed in sec. (3.1) and allow to trade the most subleading near-boundary parameters with the four independent first integrals 2 , K 3 . Although the number of arbitrary parameters remains the same, this procedure simplifies many expressions.
As we can see from equations (A.30), the first integrals are functions of the nearhorizon parameters α and η and are fully determined by these. As a consequence, we can determine the parameters a 4 , H 4 , K 4 , a 6 in terms of the remaining near-boundary ones and the near-horizon ones α, η. All the near-boundary parameters should be fully determined in terms of the near-horizon ones; if we had more first integrals at our disposal we could use them to get relations similar to (A.29) for the other near-boundary parameters and obtain the full dependence of the latter in terms of α, η, ξ . However we are not able to find any other first integral (if any exists), therefore we should resort to a numerical procedure to determine the dependence of the remaining near-boundary parameters on the near-horizon ones. This is done in [28] and also in [27,34].

B Holographic renormalization
In this appendix we report some details about how to compute some relevant physical quantities using holographic renormalization and we show their explicit expressions. We divide this appendix into two sections; in the first one we show the procedure to follow in order to obtain the renormalized Lorentzian on-shell action, while in the second one we explicitly compute the energy-momentum tensor, the holographic electric current and the scalar one-point functions. As already stated in the main text, we will always work with the Fefferman-Graham radial coordinate r, presented in app. A.1.1. We furthermore adopt a minimal subtraction scheme where the only counterterms added are those necessary to cancel the divergences. We will follow closely the notation of [28].

B.1 The renormalized on-shell action
As we explained in sec. 4, we can define a Lorentzian renormalized on-shell action as S ren = lim where the regularized action is S reg = S bulk + S GH + S ct . (B. 2) The first term coincides with the bulk action (2.7). Using the trace of the Einsten equations (2.10a) and combining the result with the Maxwell equation (2.10b), we can rewrite S bulk as: Recalling the form of the scalar potential (2.6), we can show that the following relation holds where p is the Ricci potential given in (2.23a). Looking at the expression of the latter and recalling the near-horizon expansion of a given in (A.25), we can immediately conclude that a 2 p → 0 at the horizon; therefore thanks to eq. (B.4) it is trivial to integrate the first term of (B.3) on M ρ and the only contribution will be given by the upper limit of integration. We set this upper limit to be r 0 : this is just a cutoff we use to regulate the large-distance divergences of the various pieces of the action. The second term of (B.3) can also be reduced to a boundary one via a De Rham theorem since Q IJ A I ∧ F J vanishes at the horizon in the gauge we have chosen and is regular everywhere outside the horizon. Therefore we rewrite the bulk action as S bulk = − 16π 2 3κ 2 a 2 p r 0 dt + which precisely coincides with eq. (4.2) shown in sec. 4. We find that the stress-energy tensor can be written as T ij dx i dx j = T tt dt 2 + T θθ σ 2 1 + σ 2 2 + T ψψ σ 2 3 + 2 T tψ dt σ 3 , (B. 15) where the components explicitly read 2 +K 18K + v 2 − 1 +K 6K (2) 2 + 3K (1) 2 +K v 2 − 1 ,
Similarly to the stress-energy tensor, we can define the holographic conserved currents as By varying the action with respect the boundary gauge fields, as prescribed by the formula above, we obtain: (B.22) Evaluating eq. (B.22) on our supergravity background, we find that the two non-vanishing components of the conserved currents are (B.25) We notice that j t 1 ↔ j t 2 ifH ↔K and K