BPS black holes from massive IIA on S$^6$

We present BPS black hole solutions in a four-dimensional $\mathcal{N}=2$ supergravity with an abelian dyonic gauging of the universal hypermultiplet moduli space. This supergravity arises as the SU(3)-invariant subsector in the reduction of massive IIA supergravity on a six-sphere. The solutions are supported by non-constant scalar, vector and tensor fields and interpolate between a unique $\textrm{AdS}_{2} \,\times\, \textrm{H}^2$ geometry in the near-horizon region and the domain-wall DW$_{4}$ (four-dimensional) description of the D2-brane at the boundary. Some special solutions with charged AdS$_{4}$ or non-relativistic scaling behaviours in the ultraviolet are also presented.


Motivation and outlook
The search for BPS black hole solutions in four-dimensional N = 2 gauged supergravities with an embedding in string/M-theory has recently captured new attention in light of the gravity/gauge correspondence.
An interesting program started with the classification of asymptotically AdS 4 black holes in N = 2 supergravity coupled to vector multiplets in the presence of U(1) Fayet-Iliopoulos (FI) gaugings and non-constant scalars [1,2]. The case with three vector multiplets (STU model), a square root prepotential and all the FI parameters identified, corresponds to the U(1) 4 -invariant subsector [3,4] of the maximal SO(8)-gauged supergravity [5]. This supergravity arises from the reduction of eleven-dimensional supergravity on a seven-sphere [6], and has a maximally supersymmetric AdS 4 solution dual to the three-dimensional ABJM superconformal field theory [7] at low Chern-Simons (CS) levels k and −k . When uplifted to eleven dimensions, this solution corresponds to the Freund-Rubin AdS 4 × S 7 vacuum [8] describing the near-horizon geometry of the M2-brane. A charged version of this AdS 4 vacuum corresponds to the ultraviolet behaviour of the BPS black holes constructed in [1,2] (see [9,10,11] for M-theory models also containing hypermultiplets). In contrast, the infrared behaviour approaches an AdS 2 × S 2 geometry with the scalars determined by the attractor mechanism [2,12,13]. The holographic interpretation is an RG flow across dimensions, more specifically, between a CFT 3 and a CFT 1 . Using supersymmetric localisation techniques, a counting of microstates of BPS black holes in AdS 4 was performed in the dual field theory [14,15] -identified as a deformation of the ABJM theory by a topological twist [16] -and it was shown to match the Bekenstein-Hawking entropy [17,18].
This correspondence also has a realisation on the D3-brane of the type IIB theory, once the latter is reduced on a five-sphere to a five-dimensional maximal SO(6)-gauged supergravity [19]. In this case, solutions interpolating between AdS 5 and AdS 3 × Σ 2 geometries, with Σ 2 being a Riemann surface, have a holographic interpretation in terms of RG flows between a CFT 4 and a CFT 2 [20,21]. The field theory dual is a topologically twisted N = 4 super Yang-Mills theory (SYM). spin gravity multiplet vector multiplet universal hypermultiplet 2 g µν 1 A 0 µ A 1 µ 0 χ , ϕ φ , a , ζ ,ζ Table 1: Bosonic fields in the N = 2 and SU(3)-invariant sector of the maximal supergravity multiplet in four dimensions.
The present paper continues this program and classifies BPS black hole solutions in an N = 2 subsector of the four-dimensional maximal ISO(7)-gauged supergravity constructed in [22]. This supergravity arises from a reduction of the massive IIA theory on a six-sphere [23,24]. We focus on the SU(3)-invariant subsector which is described by an N = 2 supergravity coupled to a vector multiplet and the universal hypermultiplet (see Table 1). Because of the massive IIA origin, this setup differs from the M-theory and type IIB cases discussed before. For instance, the massive IIA theory has a DW 4 domain-wall solution (instead of an AdS 4 vacuum) as the four-dimensional description of the near-horizon limit of the D2-brane [25]. Such a DW 4 solution is the non-conformal analog of the AdS 4 (AdS 5 ) vacuum in the M-theory (type IIB) models, and thus controls the ultraviolet behaviour of generic BPS flows.
In this paper we present a two-parameter family of BPS black hole solutions that feature a unique AdS 2 ×H 2 geometry in the infrared and flow to a charged version of the DW 4 solution describing the D2-brane in the ultraviolet. The scalar fields in the vector multiplet and hypermultiplet are non-constant along the flow and enter the black hole horizon as dictated by the attractor equations. For specific values of the parameters, the solutions flow to either an N = 2 charged AdS 4 vacuum or to a non-relativistic metric in the ultraviolet, instead of the generic charged DW 4 solution. It would be very interesting to understand these flows from a dual field theory perspective using the massive IIA on S 6 /SYM-CS duality [23,26].
2 N = 2 supergravity with abelian gaugings from massive IIA Massive IIA ten-dimensional supergravity admits a consistent truncation on the six-sphere to maximal D = 4 supergravity with a dyonic ISO(7) gauging [23,22,24]. Within this truncation, there is a subsector that is invariant under the action of an SU(3) subgroup of the ISO(7) gauge group, and is given by an N = 2 supergravity coupled to a vector multiplet and the universal hypermultiplet [22]. The dynamical (bosonic) degrees of freedom of this N = 2 subsector are summarised in Table 1.
We follow closely the N = 2 supergravity conventions of [27] except for a change of gauge in the ansatz for the vector and tensor fields, to be discussed below. The two real scalars in the vector multiplet (see Table 1) can be grouped into a complex one describing the special Kähler manifold M SK = SU(1,1)/U(1) in terms of holomorphic sections X M (z) = (X Λ (z), F Λ (z)) . Here M is a symplectic Sp(4) vector index, whereas Λ = 0, 1 runs over the first (electric) or second (magnetic) half of components. It proves convenient to define a symplectic product of vectors where Ω M N is the antisymmetric invariant matrix of Sp (4). In terms of it, the Kähler potential associated to M SK can be expressed as K = − log(i X,X ) . In the model of [23] the sections take the form and satisfy the relation F Λ = ∂F/∂X Λ for a prepotential F of the form whereas the Kähler potential yields a Kähler metric of the form The generalised theta angles and coupling constants for the vector fields entering the Lagrangian are encoded in a complex matrix that depends only on the scalar z Extracting R ΛΣ ≡ Re(N ΛΣ ) and I ΛΣ ≡ Im(N ΛΣ ) from (2.6), we introduce a scalar matrix M M N (z) that restores symplectic covariance and will be relevant later on when presenting the BPS equations. It takes the form is a redefined (non-holomorphic) set of symplectic sections with Kähler covariant derivatives given by Consider now the universal hypermultiplet M QK = SU(2,1)/(SU(2)×U(1)). The four real scalars spanning this quaternionic Kähler geometry are collectively denoted q u = (φ, a, ζ,ζ), with metric The specific N = 2 models that we focus on involve an abelian R × U(1) gauging of two isometries of this quaternionic manifold. The relevant Killing vectors k α (where α = R or U(1)) are and can be derived from an SU(2) triplet of moment maps P x α of the form The gaugings under consideration in this work are of the dyonic type first introduced in [28] and further explored in [29]. These gaugings involve both electric A µ Λ and magnetic A µ Λ vector fields as gauge connections in the covariant derivatives. The vector fields can be arranged into an Sp(4) symplectic vector A µ M = (A µ Λ ,Ã µ Λ ) in terms of which the covariant derivatives for the scalars in the hypermultiplet read Following [27], we have introduced Killing vectors of the form K M u ≡ Θ M α k α u in (2.11) in order to restore symplectic covariance.
The embedding tensor Θ M α in (2.11) is constant and specifies the linear combinations of electric and magnetic vectors that enter the gauge connection. Consistency requires a quadratic constraint on the embedding tensor of the form Θ α , Θ β = 0 [30]. This constraint can be viewed as an orthogonality condition between the charges Θ M α in (2.11), and guarantees that a dyonic gauging involving electric and magnetic vectors can always be rotated back to a purely electric one by a change of symplectic frame. This change of symplectic frame is usually assumed in the literature in order to have a description involving electric vectors solely. However, a formulation in terms of a prepotential F might be no longer available after changing the symplectic frame. In this work, we stay with the prepotential in (2.4) and do not perform any symplectic rotation to an electric frame. As a result, we deal with dyonic gaugings involving non-zero magnetic charges Θ Λα .
Consistency of the gauge algebra in the presence of magnetic charges requires one to introduce auxiliary two-form tensor fields B µν α that modify the field strengths of the dynamical vectors. For abelian gaugings, the latter are given by [30] H µν Lastly, the tensor fields come along with their own set of tensor gauge transformations, which are intertwined with the ordinary vector gauge transformations. We will discuss the gauge fixing of this symmetry in the next section. Using differential form notation, the bosonic Lagrangian that describes the dynamics of the dyonic gaugings of N = 2 supergravity reads [27] where the last line is a topological term that is non-zero whenever magnetic charges Θ Λα are present. Together with the Einstein-Hilbert term, and due to the abelian gauging in the hypermultiplet sector, the Lagrangian also contains a scalar potential V g given by where, as for the Killing vectors entering (2.11), we have now introduced a symplectic vector of momentum maps P x M ≡ Θ M α P x α in order to restore symplectic covariance [27]. Therefore, the Lagrangian (2.13) becomes completely specified in terms of the geometric data for M SK and M QK presented previously (Killing vectors, etc.), as well as a constant embedding tensor Θ M α encoding the gauging of the theory.
The model of [23] The N = 2 dyonically gauged supergravity we explore in this work appears from the reduction of massive IIA supergravity on the six-sphere [23,22,24]. These gaugings are determined by an embedding tensor Θ M α of the form where g and m are constant parameters identified with the inverse radius of the six-sphere and with the Romans mass parameter, respectively, and are assumed to be positive. The parameter g sources the electric part of the embedding tensor whereas the parameter m activates the magnetic one. By setting m = 0 , the gauging is of electric type and the resulting N = 2 supergravity model has an uplift to the massless IIA theory (and thus also to M-theory).
From the explicit form of the embedding tensor in (2.15) it follows that the R factor in the gauge group R × U(1) is gauged dyonically by the vectors A 0 andÃ 0 , whereas the U(1) factor is gauged only electrically by the vector A 1 . This can be seen from the covariant derivatives (2.11) of the scalars in the universal hypermultiplet which, for our specific model, take the form As a result, the shift symmetry associated with the Killing vector k R = ∂ a in (2.9) is gauged with the linear combination α − ≡ g A 0 − mÃ 0 of the graviphoton and its magnetic dual, whereas that of the k U(1) Killing vector is gauged using the vector A 1 in the vector multiplet, and the scalars ζ andζ are charged under it. The model also contains a tensor field that modifies the electric field strengths according to (2.12), resulting in where we have relabelled the tensor field as B 0 ≡ B R . Therefore, the scalar a in (2.16) is a Stückelberg field, and the tensor field B 0 becomes massive. Since the U(1) factor of the gauge group is gauged electrically only, the tensor field B U(1) decouples from the system and can be consistently set to zero. When particularised to the embedding tensor in (2.15), the generic N = 2 supergravity Lagrangian in (2.13) becomes It is important to note that the dyonic nature of the gauging implies the introduction of the magnetic vectorÃ 0 and the tensor field B 0 which, however, does not affect the counting of degrees of freedom. These fields are not dynamical, as can be seen from the variations of the Lagrangian (2.18) with respect to them, which produce two first-order differential relations The former is a duality relation between the tensor field and the scalars in the universal hypermultiplet, whereas the later is the duality relation between the graviphoton and its magnetic dual. As anticipated below (2.12), the introduction of the tensor field comes along with an additional tensor gauge symmetry given by a one-form gauge parameter Ξ 0 . Up to a total derivative, the Lagrangian (2.18) is invariant under the tensor gauge transformation Finally, plugging the embedding tensor (2.15) into the expression of the scalar potential in (2.14), and making again use of the scalar geometry data, one obtains (2.21) The full set of equations of motion that follows from the N = 2 supergravity Lagrangian (2.18) is presented in appendix A.

BPS equations in dyonically gauged N = 2 supergravity
The generic Lagrangian (2.13) of dyonically gauged N = 2 supergravity has recently been considered in [27] to study static BPS flow equations with spherical S 2 (κ = 1) or hyperbolic H 2 (κ = −1) symmetry. In this section we make extensive use of the results derived therein, and simply fetch the main results and equations needed to find BPS solutions in our model.

Field ansatz and gauge fixing
The most general metric compatible with sphericity/hyperbolicity and staticity is given by where we have partially-fixed diffeomorphisms by imposing that the radial component of the metric is the inverse of the temporal one. The functions U (r) and ψ(r) are assumed to depend solely on the radial coordinate r , and the same holds for the scalar fields z(r) and q u (r) . As we show below (see eq. (4.6)), the existence of a regular horizon in the infrared (IR) imposes that the scalars ζ andζ must vanish there. Furthermore, we will impose boundary conditions in the ultraviolet (UV) such that ζ andζ vanish at r → ∞ . Then, by looking at the equations of motion in (A.4) and at the form of V g in (2.21), it is consistent to take From now on we restrict our study to configurations where this relation is imposed, which allows us to simplify the forthcoming discussion. This restriction also implies an enhancement of the residual symmetry of the SU(3)-invariant sector of maximal supergravity to an SU(3) × U(1) symmetry as a consequence of turning off the scalar fields charged under the U(1) factor of the gauge group. Let us consider now the ansatz for the vector and tensor fields. For the vectors, staticity and spherical/hyperbolic symmetry of the associated field strengths imply that with p Λ being the constant magnetic charges of the electric gauge fields. We work in the gauge in which the radial components A r Λ (r) dr are set to zero. The ansatz for the magnetic vector and the tensor field are given bỹ where e 0 can be identified with a constant electric charge of A 0 upon the use of the duality relation between electric and magnetic vectors in (2.19 and we can use the last one to set a = 0 . Furthermore, the U(1) current sourcing the righthand-side of the Maxwell equation (A.2) for the A 1 vector vanishes whenever ζ =ζ = 0 . This allows to introduce the dual magnetic vector toÃ 1 satisfying such that the charge e 1 is a constant of motion. Combining (3.7) with the second equation in (2.19) we can then write duality relations between electric and magnetic vectors of the form In ref. [27], the ansatz for the tensor field was of the form B 0 dφ . By performing the gauge transformation (2.20) the vector charges in the two gauge choices are related as p 0 (r) [27] = p 0 (3.3) + 1 2 m b0(r) and e0(r) [27] = e0 (3.4) + 1 2 g b0(r) . We prefer to work with the spherically/hyperbolic symmetric form for B 0 in (3.4), which is consistent with constant charges for the vector fields.
Note that we do not need to solve forÃ t 1 as it does not enter any equation of motion. On the other hand, the integration constant e 1 makes and appearance in the first order equations (3.8). These read (3.9) The second expression in (3.9) allows to integrate out A t 1 since it appears only via radial derivatives. On the other hand, the temporal components of the electric and magnetic fields A t 0 andÃ t 0 enter the equations of motion of the remaining fields via the combination α − t = g A t 0 − mÃ t 0 . Summarising, the spherical/hyperbolic and static ansatz we have imposed reduces the equations of motion to a system of two first-order differential equations (for b 0 and α − t ) and five second-order differential equations (for φ , ϕ , χ , U and ψ ), together with a first-order constraint coming from the radial component of the Einstein equations. The equations of motion of ϕ and χ are displayed in (A.6) and (A.7). The equations of motion of U , ψ and φ simplify to (3.10)

First-order BPS equations
The equations of motion obtained from the Lagrangian (2.13) with the spherical/hyperbolic and static ansatz plugged in can be obtained from the effective one-dimensional action where the primes denote derivatives with respect to the radial coordinate r . As pointed out in the previous section (see footnote 1), the ansatz for the tensor fields in [27] differs from the one in (3.4) by a tensor gauge transformation (2.20). Consequently, our symplectic vector Q M containing the vector charges is given by where the fourth row and column are zero due to our restriction (3.2). The matrix H is non-invertible. This seems at odds with the appearance of H −1 in the effective action (3.11) but, as discussed in detail in [27], the matrix H −1 is defined to satisfy the condition H H −1 H = H , which is weaker than H −1 H = I . Finally, the one-dimensional potential V 1d is given by with V BH = − 1 2 Q T M Q being the black hole potential in N = 2 ungauged supergravity, that depends on the charges and on the scalar matrix M(z) in (2.7).
The authors of [27] also identified a real function 2|W | that solves the Hamilton-Jacobi equation for the effective action (3.11) provided a charge quantisation condition where Q x ≡ P x , Q . The complex function W is given by in terms of the central charge Z = Q, V and a superpotential L = Q x P x , V . Using |W | , and up to a total derivative, the effective action (3.11) can be written as a sum of squares yielding a set of BPS first-order equations. To integrate the BPS equations it is convenient to keep the phase β in (3.16) as a dynamical variable, although by its very definition is not independent of the other functions in (3.11). The set of BPS equations following from the effective action (3.11) then reads [27]: where A r = Im(z ′ ∂ z K) = − 3 2 e ϕ(r) χ ′ (r) is the U(1) Kähler connection in M SK . The system (3.17) must be supplemented with the charge quantisation condition in (3.15), the expression of the phase β as a function of the other scalars dictated by (3.16), and with a set of additional constraints

Black holes and BPS flows
In this section we present the attractor equations for the the near-horizon region of BPS black holes in the N = 2 supergravity model we are investigating. Then we find BPS black hole solutions for which the scalar fields both in the vector multiplet and the universal hypermultiplet vary along the radial coordinate. The generic solutions interpolate between a unique AdS 2 × H 2 geometry in the near-horizon region and the domain-wall DW 4 (fourdimensional) description of the D2-brane at r → ∞ . However, special behaviours at r → ∞ also occur when the boundary conditions at the horizon are fine tuned. All the plots presented in this section have been generated with g = m = 1 , which can always be achieved by a rescaling of the fields.

Near-horizon region and attractor equations
The near-horizon geometry of an extremal four-dimensional black hole is given by The functions e U (r) and e ψ(r) in the metric (3.1) take the form where L AdS 2 and L Σ 2 are the curvature radii of the AdS 2 and Σ 2 factors of the AdS 2 × Σ 2 near-horizon geometry. In the parameterisation (4.1) we have shifted the radial coordinate r to place the horizon at r h = 0 . Using the equations for U ′ and ψ ′ in (3.17), and plugging in the functions (4.1), one obtains e −U (Z + i κ L 2 Σ 2 L) = 0 . Since this equality has to hold for any value of the radius in the AdS 2 × Σ 2 fixed point, it follows that Assuming that the scalars enter the horizon as constants, i.e. z ′ = q u′ = 0 , it follows from (3.17) that β ′ = 0 and Q ′ = 0 . Moreover, it can be shown from (4.2) and the first relation in (3.18) that K u , V = 0 . All these consequences of the AdS 2 × Σ 2 form of the metric imply that the BPS equations (3.17) can be rewritten as the set attractor equations derived in [27] where it is understood that all scalars and b 0 are evaluated at the horizon. As for the general BPS equations, the charge quantisation condition (3.15) and the additional constraints (3.18) must be imposed. The latter constraint also imposes H Ω A t = 0 , implying that in the AdS 2 × Σ 2 region g A t 0 = mÃ t 0 and A t 1 = 0 .
Let us characterise the near-horizon geometries in the model arising from the reduction of the massive IIA theory on the six-sphere. First of all, since Q ′ (r h ) = 0 , it follows from (3.12 The (quadratic) charge quantisation condition (3.15) reduces in this case to where we have made use of the first constraint in (3.5). Here we are reinstating temporarily the scalars ζ andζ to show explicitly how the attractor equations set them to zero. This is seen from the last expression in (4.3), which in particular does not involve the charges Q . In our specific model this equation imposes and fixes all the values of the scalars at the horizon but φ h in terms of the gauging parameters. Substituting (4.6) into the charge quantisation condition (4.5) gives Plugging these results into the first and second equations in (4.3) produces a set of algebraic relations. The system has a solution only if κ = −1 (hyperbolic horizon) and the scalars, charges and radii take the values (4.8) These are related to each other by an overall change in the sign of the charges Q h → −Q h . Moreover, using the definition of the phase β given in (3.16), one finds that β h = π 3 ∓ π 2 . From now on we select the first of these solutions, namely, the one with β h = −π/6 .

Asymptotically AdS 4 solutions with charges
The same configuration of the scalar fields that we have found in the analysis of the attractor equations can be seen to extremise the scalar potential V g in (2.21). In absence of charges, this configuration supports an AdS 4 × S 6 solution of massive IIA supergravity preserving N = 2 supersymmetry and SU(3) × U(1) symmetry [23]. As a consequence of the spherical/hyperbolic symmetry, the metric functions depend explicitly on κ and take the form with L 2 AdS 4 = 3 and V * g being the value of the potential (2.21) at the extremum. Here we are denoting the radial coordinate asr since, as we show below, it is shifted by a constant with respect to the one used in the previous section.
Since the set of BPS equations (3.17) requires the quantisation condition (3.15) to be satisfied, it is clear that this solution is not captured in the present setup. However it can be shown that, in the presence of charges, there is a Reissner-Nordström-AdS like solution with the same value for the scalars and with (4.10) Substituting (4.10) into the BPS equations (3.17), one finds a one-parameter family of solutions with charges which yields a function f (Q) in (4.10) of the form (4.12) This one-parameter family of solutions corresponds to an asymptotically AdS 4 geometry with non-trivial charges turned on. Near the origin,r = 0, the solution gives rise to a naked singularity. The family admits a non-extremal generalisation by adding to the metric function e 2U in (4.10) a mass term of the form −2M/r. With this the metric is a solution of the secondorder equations of motion in appendix A (but not of the BPS equations), and the the geometry in the IR is regularised by a horizon. This indicates that the naked singularities of (4.10) are of the good type in the classification of [31]. There is a particular case of the BPS solution (4.11) with which connects with the attractor solution in (4.8). It corresponds to an extremal Reissner-Nordström black hole solution with AdS 2 × H 2 geometry in the IR. This choice of e 1 charge yields a function f (Q) in (4.10) of the form with the horizon located atr 2

BPS flows from the DW 4 to AdS 2 × H 2
We have shown that the attractor equations (4.3) select a unique configuration of charges and scalar fields, given in (4.8), such that a horizon with hyperbolic symmetry exists. This AdS 2 × H 2 geometry in the IR can be reached from a charged AdS 4 geometry in the UV yielding the extremal BH solution in (4.13)-(4.14) with constant scalars. In this section we construct numerically more BPS solutions, and show that the analytic BH-AdS geometry corresponds to a very special point within a two-dimensional parameter space of configurations. These solutions generically interpolate between an AdS 2 × H 2 geometry in the IR and a DW 4 domain-wall geometry governed by the D2-brane in the UV (see Figure 1). To understand how the UV geometry is dictated by the D2-brane, let us recall the form of such a solution in massless IIA supergravity. This is given by a metric (in Einstein frame) and a dilaton eΦ of the form In addition, there is a four-form fluxF (4) = 5 g e φ e 2(ψ−U ) sinh θ dt ∧ dr ∧ dθ ∧ dφ that is electrically sourced by the D2-brane. The dependence with the radius of the different functions is given by e 2U ∼ r but this correction is suppressed as one approaches the boundary at r → ∞ [25]. This can be seen from the potential of the corresponding four-dimensional gauged supergravity or from the fermion mass terms entering the supersymmetry transformations obtained upon reduction on S 6 . In both cases the Romans' mass parameter appears dressed up with a function of the scalars that suppresses its contribution near the boundary. In the presence of non trivial Q charges, as it is the case in this work, a similar effect occurs: the charges are dressed up with functions of the scalars that make their induced corrections subleading near the boundary. Furthermore, perturbing the BPS equations around the DW 4 geometry shows that only relevant deformations are turned on [25]. For this reason, the D2-brane solution of the massless IIA theory generically governs the UV asymptotics also in the massive setup with finite charges. In addition, having a solution whose UV is governed by the DW 4 configuration necessarily implies a running of the dilaton e φ belonging to the universal hypermultiplet. This implies that all the solutions that we describe in this section contain running hyperscalars.
In order to solve the BPS equations, we shoot numerically from the extremal horizon. To impose appropriate boundary conditions, we first identify the irrelevant perturbations around the unique AdS 2 × H 2 solution given by the metric and fields in (4.1) and (4.8). Expanding the BPS equations (3.17) near the horizon at r = 0 , one finds the following regular corrections to the metric and field functions: We have performed a numerical scan of 10 6 points in the (c 1 , c 2 )-plane within the range −100 ≤ c 1,2 ≤ 100 . The result is depicted in Figure 1, which we explain now in some detail. The shaded region corresponds to regular BPS configurations that interpolate between the AdS 2 × H 2 solution (4.1) and (4.8) in the IR, and flow to the DW 4 solution (4.16) in the UV. All these configurations have the same behaviour at large r given by (4.16) together with where we are omitting corrections that fall off at r → ∞ and a coefficient that depends on the specific choice of (c 1 , c 2 ) . Importantly for the D2-brane interpretation, the two dilatons e ϕ and e φ become identified asymptotically and the axion χ goes to zero faster than the dilatons as r increases. The BPS solution with (c 1 , c 2 ) = (−1, −1) is represented in Figure 1 by a (red) triangle, and the profiles for the corresponding fields are shown in Figure 2. Note that, despite this solution having c 1 = c 2 , the function b 0 still flows non trivially as it receives a correction at a larger order than the one given in (4.17). Even though b 0 ∼ r 1/2 asymptotically, the solution approaches the DW 4 , thus indicating that this mode does not carry infinite energy at the boundary. Nonetheless, it is possible to tune the values of (c 1 , c 2 ) to find solutions such that b 0 approaches a constant when r → ∞ . We have denoted the locus of such parameters with the (grey) dashed line in Figure 1. The shaded region of regular solutions in Figure 1 is delimited. The upper (red line) and lower (brown line) boundaries yield configurations that do not approach the DW 4 solution (4.16) but acquire non-relativistic behaviours in the UV. For instance, the (blue) circle approaches a Lifshitz spacetime with z = 2 whereas the (green) square approaches a conformally Lifshitz spacetime with (z, θ) = (1.86, −0.705) . Lastly, the (black) rhombus at the origin of the parameter space (c 1 , c 2 ) = (0, 0) is special and produces the asymptotically AdS 4 solution with constant scalars in (4.13)-(4.14). This is the only point in Figure 1 satisfying c 1 + 3 c 2 = 0 , or equivalently, setting to zero the irrelevant deformations in (4.17) for the dilaton e φ in the universal hypermultiplet. Moving slightly away from this point into the shaded region modifies the UV behaviour of the solution making it flow to the DW 4 . We show this behaviour in Figure 3 where we have produced the plot by setting (c 1 , c 2 ) = (0, −10 −8 ) . One sees that the logarithmic derivatives of the metric functions coincide quite accurately with the ones dictated by the asymptotically AdS 4 solution in (4.13)-(4.14) (red, dashed line) up to a value of the radial coordinate beyond which the functions in our ansatz transition to that of the DW 4 asymptotics (4. 16).

Non-relativistic UV asymptotics
As previously mentioned, the solutions associated with the points at the boundary of the shaded region in Figure 1 have a non-relativistic scaling in the UV. An example of this behaviour is given by the (blue) circle in that figure, for which the BPS solution asymptotes a scaling solution with broken Lorentz symmetry e 2U ∼ r 2 , e 2(ψ−U ) ∼ r , β ∼ 0 , b 0 ∼ r , (4.20) and constant scalars at large values of the radial coordinate. This corresponds to a nonrelativistic metric of the Lifshitz type with dynamical exponent z = 2. Along the boundary line that joins the (blue) circle and the (black) rhombus from above (red line), the scaling solution (4.20) receives some logarithmic corrections that we have not investigated in detail. A different non-relativistic scaling in the UV occurs for solutions associated with the points in the boundary line connecting the (blue) circle and the (black) rhombus in