Black holes by gravitational decoupling

We investigate how a spherically symmetric fluid modifies the Schwarzschild vacuum solution when there is no exchange of energy-momentum between the fluid and the central source of the Schwarzschild metric. This system is described by means of the gravitational decoupling realised via the minimal geometric deformation approach, which allows us to prove that the fluid must be anisotropic. Several cases are then explicitly shown.


Introduction
The study of black holes represents one of the most active areas of gravitational physics, from both a purely theoretical and the observational point of view. The interest black holes generate is due not only to their exotic nature, but also because they constitute ideal laboratories to study gravity in the strong field regime, and test general relativity therein. However, confronting theoretical predictions with observations is an arduous and complicated task. A formidable step in this direction is the recent direct observation of black holes through the detection of gravitational waves, which opens a new and promising era for gravitational physics [1,2].
It is well known that general relativity predicts surprisingly simple solutions for black holes, characterised at most by three fundamental parameters, namely the mass M, angular momentum J and charge Q [3]. The original no-hair conjecture states that these solutions should not carry any a e-mail: jovalle@usb.ve b e-mail: casadio@bo.infn.it c e-mail: roldao.rocha@ufabc.edu.br d e-mail: adrian.sotomayor@uantof.cl e e-mail: zdenek.stuchlik@fpf.slu.cz other charges [4]. Therefore, as the observations of systems containing black holes improve, the degree of consistency of these observations with the predictions determined according to the general relativistic solutions (with parameters M, J and Q) will result in a direct test of the validity of general relativity in the strong field regime. There could in fact exist other charges associated with inner gauge symmetries (and fields), and it is now known that black holes could have (soft) quantum hair [5]. The existence of new fundamental fields, which leave an imprint on the structure of the black hole, thus leading to hairy black hole solutions, is precisely the scenario under study in this paper.
Possible conditions for circumventing the no-go theorem have been investigated for a long time in different scenarios (see Refs. [6][7][8][9][10][11][12][13][14][15] for some recent works and Refs. [16][17][18][19][20][21][22] for earlier works). In particular, a fundamental scalar field φ has been considered with great interest (see Ref. [23] and references therein). In this work, we will take a different and more general approach than most of the investigations carried out so far and, instead of considering specific fundamental fields to generate hair in black hole solutions, we shall just assume the presence of an additional completely generic source described by a conserved energy-momentum tensor θ μν . Of course, this θ μν could account for one or more fundamental fields, but the crucial property is that it gravitates but does not interact directly with the matter that sources the (hairless) black hole solutions we start from. This feature may seem fanciful, but can be fully justified, for instance, in the context of the dark matter. Achieving this level of generality in the classical scheme represented by general relativity is a non-trivial task, and the gravitational decoupling by Minimal Geometric Deformation (MGD-decoupling, henceforth) is precisely the method that was developed for this purpose in Ref. [24].
The MGD approach was originally proposed [25,26] in the context of the brane-world [27,28] and extended to investigate new black hole solutions in Refs. [29,30] (for some earlier works on the MGD, see for instance Refs. [31][32][33][34], and Refs. [35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50] for some recent applications). The MGDdecoupling has a number of ingredients that make it particularly attractive in the search for new spherically symmetric solutions of Einstein's field equations. The two main feature of this approach are the following [24]: • Extending simple solutions into more complex domains We can start from a simple spherically symmetric gravitational source with energy-momentum tensorT μν and add to it more and more complex gravitational sources, as long as the spherical symmetry is preserved. The starting sourceT μν could be as simple as we wish, including the vacuum indeed, to which we can add a first new source, saŷ where α (1) is a constant that traces the effects of the new source T (1) μν . We can then repeat the process with more sources, namelỹ and so on. In this way, we can extend straightforward solutions of the Einstein equations associated with the simplest gravitational sourceT μν into the domain of more intricate forms of gravitational sources T μν =T (n) μν , step by step and systematically. We stress that this method works as long as the sources do not exchange energymomentum among them, namely which further clarifies that the constituents can only couple via gravity. • Deconstructing a complex gravitational source. The converse of the above also works. In order to find a solution to Einstein's equations with a complex spherically symmetric energy-momentum tensor T μν , we can split it into simpler components, sayT μν and T (i) μν , provided they all satisfy Eq. (1.3), and solve Einstein's equations for each one of these parts. Hence, we will have as many solutions as are the contributions T (i) μν in the original energymomentum tensor. Finally, by a straightforward combination of all these solutions, we will obtain the solution to the Einstein equations associated with the original energy-momentum tensor T μν .
Since Einstein's field equations are non-linear, the MGDdecoupling represents a breakthrough in the search and analysis of solutions, especially when we deal with situations beyond trivial cases, such as the interior of self-gravitating systems dominated by gravitational sources more realistic than the ideal perfect fluid [51,52]. Of course, we remark that this decoupling occurs because of the spherical symmetry and time-independence of the systems under investigation.
In analogy with the well-known electro-vacuum and scalar-vacuum, in this paper we will consider a Schwarzschild black hole surrounded by a spherically symmetric "tensorvacuum", represented by the aforementioned θ μν . Following the MGD-decoupling, we will separate the Einstein field equations in (i) Einstein's equations for the spherically symmetric vacuum and (ii) a "quasi-Einstein" system for the spherically symmetric "tensor-vacuum". The MGD procedure will then allow us to merge the Schwarzschild solution for (i) with the solution for the "quasi-Einstein" system (ii) into the solution for the complete system "Schwarzschild + tensor-vacuum". Like the case of the electro-vacuum and (in some cases) scalar-vacuum, new black hole solutions with additional parameters q i besides the mass M can be obtained, each one associated with a particular equation of state for the "tensor-vacuum". Demanding the geometry is free of singularities and other pathologies, implies regularity conditions which show that not all of these parameters q i can be independent.
The paper is organised as follows: in Sect. 2, we first review the fundamentals of the MGD-decoupling applied to a spherically symmetric system containing a perfect fluid and an additional source θ μν ; in Sect. 3, new hairy black holes solutions are found by assuming the perfect fluid has sufficiently small support so that only θ μν exists outside the horizon; finally, we summarise our conclusions in Sect. 4.

MGD decoupling for a perfect fluid
Let us start from the standard Einstein field equations and assume the total energy-momentum tensor contains two contributions, namely where is the 4-dimensional energy-momentum tensor of a perfect fluid with 4-velocity field u μ , density ρ and isotropic pressure p. The term θ μν in Eq. (2.2) describes an additional source whose coupling to gravity is proportional to the constant α [53]. This source may contain new fields, like scalar, vector and tensor fields, and will in general produce anisotropies in self-gravitating systems. In any case, since the Einstein tensor satisfies the Bianchi identity, the total source in Eq. (2.2) must satisfy the conservation equation We next specialise to spherical symmetry and no timedependence. In Schwarzschild-like coordinates, a static spherically symmetric metric g μν reads where ν = ν(r ) and λ = λ(r ) are functions of the areal radius r only, ranging from r = 0 (the star center) to some r = R (the star surface), and the fluid 4-velocity is given by The metric (2.5) must satisfy the Einstein equations (2.1), which explicitly read where f ≡ ∂ r f and spherical symmetry implies that θ 3 3 = θ 2 2 . The conservation equation (2.4) is a linear combination of Eqs. (2.6)-(2.8), and yields We then note the perfect fluid case is formally recovered for α → 0. The Eqs. (2.6)-(2.8) contain seven unknown functions, namely: two physical variables, the density ρ(r ) and pressure p(r ); two geometric functions, the temporal metric function ν(r ) and the radial metric function λ(r ); and three independent components of θ μν . This system of equations is therefore indeterminate and we should emphasise that the space-time geometry does not allow one to resolve for the gravitational source {ρ, p, θ μν } uniquely.
In order to simplify the analysis, and by simple inspection, we can identify an effective densitỹ an effective radial pressurẽ and an effective tangential pressurẽ These definitions clearly illustrate that θ μν could in general induce an anisotropy, inside a stellar distribution sourced by T (m) μν alone. The system of Eqs. (2.6)-(2.8) may therefore be formally treated as an anisotropic fluid [54,55].
The MGD-decoupling 1 can now be applied to the case at hand by simply noting that the energy-momentum ten- The components of the diagonal metric g μν that solve the complete Einstein equations (2.1) and satisfy the MGD read [24] g μν =ĝ μν = g (1) μν (2.14) for μ = ν = 1, and g 11 =ĝ 11 + α g (1)11 , (2.15) so that only the radial component is affected by the additional source θ μν . This metric g μν is found by first solving the Einstein equations for the perfect fluid source T and then the remaining quasi-Einstein equations for the source θ μν , namelỹ where the divergence-free quasi-Einstein tensor with ν μ a tensor that depends exclusively on g μν to ensure the divergence-free condition. For the spherically symmetric metric (2.5), it reads We can then proceed by considering a solution to the Eqs. (2.16) for a perfect fluid [that is Eqs. (2.6)-(2.9) with α = 0], which we can write as is the standard General Relativity expression containing the Misner-Sharp mass function m(r ). The effects of the source θ μν on the perfect fluid solution {ξ, μ ρ, p} can then be encoded in the MGD undergone solely by the radial component of the perfect fluid geometry in Eq. (2.20). Namely, the general solution is given by Eq. (2.5) with ν(r ) = ξ(r ) and where f * = f * (r ) is the MGD function to be determined from the quasi-Einstein Eqs. (2.17), which explicitly read We also notice that the conservation equations for the additional energy-momentum tensor, ∇ μ θ μν = 0, yield which does not depend on the MGD function f * .
In the next section, we shall solve the above equations starting from the simplest vacuum solution given by the outer Schwarzschild metric, therefore in a region of space where the perfect fluid ρ and p vanish.

Black holes
When new paradigms beyond Einstein gravity are studied, the important question arises whether or not new black hole solutions exist. In order to address this point in general, we start from the results of the previous section and determine the MGD function f * for the vacuum Schwarzschild solution (2.27). Figure 1 schematically shows the kind of system we deal with. The MGD metric will therefore read where the MGD function f * can be determined by imposing restrictions on the energy-momentum θ μν to close the system of Eqs. (2.23)-(2.25).
In the following we shall explore specific equations of state for θ μν and impose basic constraints on the causal structure of the resulting space-time in order to have a welldefined horizon structure. In particular, we recall that for the Schwarzschild metric (2.27), the surface r H = 2 M is both a Killing horizon (determined by the condition e ν = 0) and an outer marginally trapped surface (the causal horizon, in brief, determined by the condition e −λ = 0). For the MGD Schwarzschild metric (3.1), the component g tt = e ν Fig. 1 Spherically symmetric source θ μν in the vacuum ρ = p = 0 is always equal to the Schwarzschild form in Eq. (2.27) and can only vanish at r = r H . This means that r H = 2 M is still a Killing horizon (which can also become a real singularity). However, the causal horizon is found at r = r h such that g rr (r h ) = e −λ = 0, or We should therefore require that r h ≥ 2 M, so that the surface r = r H is hidden behind (or coincides with) the causal horizon. Moreover, if r h > r H , the signature of the metric becomes (+ + −−) for r H < r < r h , which one might want to discard as well, since not only the expansion of outgoing geodesics vanishes for r → r + h , but also ingoing geodesics never cross r = r h : in this case the surface r = r h would act as a border of the outer space-time manifold. To summarise, the MGD metric (3.1) represents a proper black hole only if the causal horizon coincides with the Killing horizon, that is when r h = r H = 2 M, and this is therefore the condition we shall require in the following.

Isotropic sector
Let us start by considering the case of isotropic pressure, so that Equations (2.24) and (2.25) then yield a differential equation for the MGD function, namely whose general solution is given by where iso is a constant with dimensions of a length. Hence, the MGD radial component for an isotropic deformation of the Schwarzschild exterior becomes which is clearly not asymptotically flat for r M. 2 We therefore conclude that the additional source θ μν cannot contain an isotropic pressure if we wish to preserve asymptotic flatness. 2 In fact, it approaches the radial component of the de Sitter metric for r ∼ iso M.

Conformal sector
The energy-momentum tensor for a conformally symmetric source must be traceless. Since θ 2 2 = θ 3 3 , we therefore assume so that the system (2.23)-(2.25) becomes where f * is again MGD function and ξ the unperturbed Schwarzschild function. From Eq. (3.7), we find the radial deformation must satisfy the differential equation and it is important to highlight that the conservation equation (2.26) remains a linear combination of the system (3.8)-(3.7). The general solution for Eq. (3.10) is given by with c a constant with units of a length. Thus the conformally deformed Schwarzschild exterior becomes where = α c , and its behaviour for r M is given by The causal structure for this geometry is now more involved. We still have the Killing horizon of the Schwarzschild metric at r H = 2 M, but e −λ diverges for and there is a second zero of e −λ at We can thus rewrite the radial metric component as Note that r c < r H but, depending on the sign ands size of , the second zero could occur inside or outside the critical radius r c and the Killing horizon r H .
In order to clarify the nature of the above solution, we compute explicitly the effective densitỹ the effective radial pressurẽ 18) and the effective tangential pressurẽ The anisotropy is thus given by (3.20) The first thing we notice is that the density and pressures are regular on both r H and r 0 , but diverge at r = r c < r H , which is therefore a real singularity, albeit hidden inside the Killing horizon. We can then assume the black hole space-time is represented by the range r > r c , for which we must require that the region r c < r < r H has the proper signature, as discussed previously. This means that we must have r 0 ≤ r c , or > 0, (3.21) with = 0 of course representing the pure Schwarzschild geometry. We conclude that the conformal geometry in (3.12) represents a black hole solution with outer horizon r H = 2 M, and primary hairs represented by the parameter , which is constrained by the regularity condition (3.21). A solution similar to that in (3.12) was found in the context of the extradimensional brane-world [56].

Barotropic equation of state
If the additional source is a polytropic fluid, it should satisfy the equation of statẽ p r = Kρ , (3.22) with = 1 + 1/n, where n is the polytropic index and K > 0 denotes a parameter which contains the temperature implicitly and is governed by the thermal characteristics of a given polytrope. For instance, the ultrarelativistic degenerate Fermi gas has polytropic index n = 3, while the nonrelativistic degenerate Fermi gas is found for n = 3/2. (for more details, see for instance Refs. [57][58][59][60][61]). However, due to the unknown nature of the source θ μν , we will include the possibility that K < 0. From Eqs. (2.10) and (2.11) with ρ = p = 0, we then obtain By using Eqs. (2.23) and (2.24) in the expression (3.23) we obtain a first order non-linear differential equation for the MGD function, We immediately notice that the right hand side is well-defined for a generic only provided f * /(r − 2 M) > 0.
In order to proceed, we thus consider the simplest case = 1, so that Eq. (3.23) becomes a barotropic equation of state. This corresponds to an isothermal self-gravitating sphere of gas and is thus more appropriate for our purpose. It is worth mentioning that this self-gravitating sphere can also describe the collisionless system of stars in a globular cluster. The geometric deformation for = 1 and r > 2 M simplifies to where p > 0 is a length, and the MGD radial component of the metric reads (3.26) again for r > 2 M. We also note that asymptotic flatness at r → ∞ requires K ≤ −1, with K = −1 yielding the pure Schwarzschild metric. Next, we note that the effective density is given by (3.27) and diverges at r = 2 M unless −1 < K < 0. Of course, the effective pressurep r = Kρ also diverges at r = 2 M unless −1 < K < 0. The effective tangential pressure is given by (3.28) which also diverges at at r = 2 M unless −1/2 < K < 0. To summarise, the surface r = r H is a real singularity unless −1/2 < K < 0. However, this is not compatible with the asymptotically flat conditions, which requires K ≤ −1. We therefore conclude that the Killing horizon at r = r H = 2 M has become a real singularity, which is not hidden inside a larger horizon.

Linear equation of state
Now let us consider a generic equation of state in the form whose general solution for r > r H = 2 M is given by where is a positive constant with dimensions of a length, and with b = 2 and the condition A > 0 required by asymptotic flatness. Therefore the solution reads which again shows the horizon at r H = 2 M, beside a possible divergence at r = r c and a second zero at r = r 0 , like in the previous cases. The physical content of the system is again clarified by the explicit computation of the effective densitỹ the effective radial pressurẽ (3.36) and the effective tangential pressurẽ (3.37) Again, we see that the effective density and effective pressures diverge at which represents a true singularity at 0 < r c < r H for 0 < B < 2, that is for For B > 2 (that is, 0 < b < 2), this singularity occurs outside the Killing horizon, r c > r H , and this case cannot be considered any further. Secondly, the effective density and effective pressures satisfỹ showing thatp t < 0 when bothρ andp r are positive. We thus conclude that at least one of the thermodynamic variables will always be negative as long as the equation of state is linear. Moreover, the effective radial and tangential pressure are related bỹ Since A > 0, we conclude that the two pressures always have opposite signs and one of them will be negative. On the other hand, the effective density and effective radial pressure are related bỹ Since A > 0, the asymptotic behaviour in Eq. (3.43) demands 0 < A ≤ 1 in order to ensure that the density does not change its sign in the region 2 M < r < ∞ [the pressure (3.36) always has the same sign inside this region]. We conclude that the dominant energy conditionρ ≥|p r | cannot be satisfied with a linear equation of state of the form displayed in Eq. (3.29). Nonetheless, the effective density is positive everywhere if α < 0 3 . For 2 < b < 4 one has r c < 0 and there is no extra singularity beside the usual Schwarzschild one at r = 0. In this case, we must demand that no second zero r 0 > 0 of e −λ exists, otherwise the space-time signature would become unacceptable inside a portion of r > 0. This condition is immediately satisfied if α > 0, for any A > 0, that is for a > 1. We next notice that there is a second zero of e −λ at when α < 0. To have a proper black hole solution, this second zero r 0 must lie inside the relevant singularity. If 2 < b < 4, the relevant singularity occurs at r = 0 an we must have r 0 ≤ 0, (3.45) that is, if and |α| are small enough to satisfy Otherwise, if b < 0 or b > 4, the relevant singularity occurs at 0 < r c < r H , but r 0 > r c makes this case unsuitable. The final conclusion is thus that the linear equation of state (3.29) always produces black holes (with a Schwarzschild singularity at r = 0) if 2 < b < 4 and a > 1, provided α > 0 or α < 0 and Eq. (3.46) holds.

A particular solution with no extra singularity
The reader can see that Eq. (3.29) leads to a system very rich in possibilities, whose generic solutions 4 are given in Eqs. (3.34)-(3.37), and whose general analysis is detailed throughout Eqs. (3.38)- (3.46). The main feature of these solutions is that they do not satisfy the dominant energy condition. In this respect, let us recall that the energy conditions are a set of constraints which are usually imposed on the energy-momentum tensor in order to avoid exotic matter sources, hence they can be viewed as sensible guides to avoid unphysical situations. However, it is well-known that these energy conditions might fail for particular classical systems which are still reasonable [62]. In our case we are dealing with a gravitational source θ μν whose main characteristic is that it only interacts gravitationally with the matter that, by itself, would source the (hairless) black hole solution (2.27). Hence, one should not exclude a priori that such matter is a kind of exotic source (as indeed the conjectured dark matter is expected to be).
Of all the possible solutions, we shall here analyse the particular case b = 3 (with a > 1 for asymptotic flatness) as an example of space-time which does not contain any extra singularity beside the usual Schwarzschild one at r = 0. The radial metric component is obtained from Eq. (3.34) and reads −1) , (3.47) which makes it immediately clear that there is no second divergence. In fact, the effective density is given by 5 the effective radial pressure bỹ from which we obtain the asymptotic behaviours (3.52) Therefore, when α < 0, the pressurep r < 0 and the effective energy will always be positive for r > 2 M whenever a ≤ 3/2. Figs. 2 and 3 show the corresponding metric elements and density and pressures in (3.48)-(3.50) respectively for α = −0.7 and a = 1.4.

Conclusions
By making use of the MGD-decoupling approach, we have presented in detail how the Schwarzschild black hole is modified when the vacuum is filled by a generic spherically symmetric gravitational fluid, described by a "tensor-vacuum" θ μν , which does not exchange energy-momentum with the central source. For this purpose, we have separated the Einstein field equations into (i) the Einstein equations for the spherically symmetric vacuum ρ = p = 0 and (ii) the "quasi-Einstein" system in Eqs. (2.23)-(2.25) for the spherically symmetric "tensor-vacuum" θ μν . Following the MGD procedure, the superposition of the Schwarzschild solution found in (i) plus the solution for the "quasi-Einstein" system in (ii), has led to the solution for the complete system "Schwarzschild + tensor-vacuum." The quasi-Einstein system (2.23)-(2.25) was solved by providing some physically motivated equations of state for the source θ μν . In this respect, four different scenarios were considered, namely, (i) the isotropic θ 1 1 = θ 2 2 ; ii) the conformal θ μ μ = 0; iii) the polytropic α θ 1 1 = K (α θ 0 0 ) and iv) the generic linear equation of state in (3.29). In the isotropic case, we only found a metric which is not asymptotically flat for r → ∞, which means that the tensor-vacuum for a black hole cannot be isotropic as long as its interaction with regular matter is purely gravitational. On the other hand, the conformal case leads to the hairy black hole solution in Eq. (3.12), whose primary hairs is represented by the length , which is constrained by the regularity condition (3.21). Among all polytropic equations of state, we have only considered the barotropic = 1, which represents a tensor-vacuum made of an isothermal self-gravitating sphere of gas. This leads to the exterior solution in Eq. (3.26) endowed with the parameters {M, α, p , K }. Since the Killing horizon r = r H = 2 M becomes a real singularity, this solution may represent the exterior of a self-gravitating system of mass M and radius R > r H but not a black hole solution.
Finally, we have analysed the generic linear equation of state in Eq. (3.29), which includes both the conformal and barotropic fluids as particular cases. This leads to the solution in Eq. (3.34), showing that even a simple linear equation of state may yield hairy black hole solutions with a rich geometry described by the parameters {M, α, , a, b}, where {α, , a, b} represents a potential set of charges generating primary hairs. In this context, a particular black hole solution with primary hairs {α, a} was found in Eq. (3.47), whose main characteristic is the absence of other singularities in the region 0 < r < ∞.
All the black holes solutions mentioned above have the horizon at r H = 2 M and primary hairs represented by a number of free parameters. However, these parameters can be restricted by demanding (i) the correct asymptotic behaviour and (ii) regularity conditions for black hole solutions free of pathologies. In this respect, there are always a potential singularity r c and a possible second horizon r h in our solutions. In order to have a proper black hole, it is necessary that r c ≤ r H to avoid a naked singularity, and r h = r H to have a metric with a proper signature. We emphasize that r h > r H yields both g tt and g rr positive inside the region r H < r < r h . All these conditions yields restrictions on potential primary hairs. For instance, the linear equation of state (3.29) always produces black holes if 2 < b < 4 and a > 1, provided α > 0 or α < 0 and Eq. (3.46) holds.
We have shown that different characteristics of the gravitational source lead to different hairy black hole solutions. Therefore, the compatibility between some of these solutions and the observations could determine the main features of the tensor-vacuum, and eventually the fundamental field(s) that constitute it. Finally, we would like to emphasize that the non-existence of an isotropic tensor-vacuum that does not exchange energy-momentum with regular matter favours scenarios with Klein-Gordon type fields φ, which naturally induce anisotropy in the Einstein field equations. These scalar fields are found in a large number of alternative theories to general relativity.