Anisotropic solutions by gravitational decoupling

We investigate the extension of isotropic interior solutions for static self-gravitating systems to include the effects of anisotropic spherically symmetric gravitational sources by means of the gravitational decoupling realised via the minimal geometric deformation approach. In particular, the matching conditions at the star surface with the outer Schwarzschild space-time are studied in great details, and we describe how to generate, from a single physically acceptable isotropic solution, new families of anisotropic solutions whose physical acceptability is also inherited from their isotropic parent.


Introduction
In a recent paper [1], the first simple, systematic and direct approach to decoupling gravitational sources in general relativity (GR) was developed from the so-called Minimal Geometric Deformation (MGD) approach. The MGD was originally proposed [2,3] in the context of the Randall-Sundrum brane-world [4,5] and extended to investigate new black hole solutions [6,7] (For some earlier works on MGD, see for instance Refs. [8,9,10,11], and for some recent applications Refs. [12,13,14,15,16,17,18]). The decoupling of gravitational sources by MGD (henceforth MGD-decoupling) is not only a novel idea, but also has a number of ingredients that make it particularly attractive in the search for new spherically symmetric solutions of Einstein's field equations, as discussed below.
The main and foremost feature of this approach is that we can start from a simple spherically symmetric gravitational sourceT µν 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ŷ and 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 energy-momentum among them, namely which further clarifies that the constituents shown in Eq. (1.1) can only couple via gravity. The converse 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, say T (i) µν , 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 energy-momentum 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 µν .
To summarise, the MGD-decoupling amounts to the following procedure: given two gravitational sources A and B, standard Einstein's equations are first solved for A, and then a simpler set of quasi-Einstein equations are solved for B. Finally, the two solutions can be combined in order to derive the complete solution for the total system A∪B. Since Einstein's field equations are non-linear, the above procedure represents a breakthrough in the search and analysis of solutions, especially when the involved situations are beyond trivial cases, such as the interior of self-gravitating systems dominated by gravitational sources more realistic than the ideal perfect fluid [19,20]. Of course, we remark that this decoupling occurs because of the spherical symmetry and time-independence of the systems under investigation.
Although decoupling gravitational sources in GR by a systematic way represents in itself a fact of significant theoretical importance, its practical relevance is no less. Indeed, this simple and systematic method could be conveniently exploited in a large number of relevant cases, such as the Einstein-Maxwell [21] and Einstein-Klein-Gordon system [22,23,24,25], for higher derivative gravity [26,27,28], f (R)-theories of gravity [29,30,31,32,33,34,35], Hořava-aether gravity [36,37], polytropic spheres [38,39,40], among many others. In this respect, the simplest practical application of the MGD-decoupling consists in extending known isotropic and physically acceptable interior solutions for spherically symmetric self-gravitating systems into the anisotropic domain, at the same time preserving physical acceptability, which represents a highly non-trivial problem [41] (for obtaining anisotropic solutions in a generic way, see for instance Ref. [42,43,44]).
The paper is organised as follows: in Section 2, we briefly review the effective Einstein field equations for a spherically symmetric and static distribution of matter with effective densityρ, effective radial pressurep r and effective tangential pressurep t ; Section 3 is devoted to the fundamentals of the MGD-decoupling; in Section 4, we provide detail on the matching conditions under the MGD-decoupling; in Section 5, we implement the MGD-decoupling to extend perfect fluid solutions in the anisotropic domain. In particular, three new exact and physically acceptable anisotropic solutions, generated from a single perfect fluid solution, are developed in order to emphasise the robustness of the approach; finally, we summarise our conclusions in Section 6.

Einstein equations for multiple sources
Let us start from the standard Einstein field equations 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 any additional source whose coupling to gravity is further proportional to the constant α [45]. This source may contain new fields, like scalar, vector and tensor fields, and will in general produce anisotropies in self-gravitating systems. We just recall that, since the Einstein tensor is divergence free, the total energy-momentum tensor (2.2) must satisfy the conservation equation In Schwarzschild-like coordinates, a static spherically symmetric metric reads ds 2 = e ν(r) dt 2 − e λ(r) dr 2 − r 2 dθ 2 + sin 2 θ dφ 2 , (2.5) 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 u µ = e −ν/2 δ µ 0 for 0 ≤ r ≤ R. The metric (2.5) must satisfy the Einstein equations (2.1), which explicitly read The conservation equation (2.4), which is a linear combination of Eqs. (2.6)-(2.8), yields where f ′ ≡ ∂ r f . We then note the perfect fluid case is formally recovered for α → 0.
The system (2.6)-(2.8) contains 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 θ µν . These equations therefore form an indefinite system. In the particular case where θ µν depends only on the density and the pressure, we need to prescribe only an additional equation to close the system Eqs. (2.6)-(2.8), just as we do for the perfect fluid in standard GR. However, at this stage we want to emphasise that it is not enough to know the space-time geometry to determine the gravitational source {ρ, p, θ µν } in general.
In order to simplify the analysis of the system (2.6)-(2.8), and by simple inspection, we can identify an effective densityρ = ρ + α θ 0 0 , (2.10) an effective isotropic pressurep and an effective tangential pressurep These definitions clearly illustrate that the source θ µν could in general induce an anisotropy, inside the stellar distribution. The system of Eqs. (2.6)-(2.8) could indeed be treated as an anisotropic fluid [46,47], which would require to consider five unknown functions, namely, the two metric functions ν(r) and λ(r), and the effective functions in Eqs. (2.10)-(2.12). However, we are going to implement a different approach, as explained below.

Gravitational decoupling by MGD
We shall implement the MGD in order to solve the system of Eqs. (2.6)-(2.9). In this approach, the system will be transformed in such a way that the field equations associated with the source θ µν will take the form of the "effective quasi-Einstein" Eqs. (3.12)-(3.14). Let us start by considering a solution to Eqs. (2.6)-(2.9) with α = 0, namely, a GR perfect fluid solution {ξ, µ, ρ, p}, where ξ and µ are the corresponding metric functions. The metric (2.5) now reads is the standard GR expression containing the mass function m. Now let us turn on the parameter α to consider the effects of the source θ µν on the perfect fluid solution {ξ, µ ρ, p}. These effects can be encoded in the geometric deformation undergone by the perfect fluid geometry {ξ, µ} in Eq. (3.1), namely where g and f are the deformations undergone by the temporal and radial metric component, respectively. Among all possible deformations (3.3) and (3.4), the so-called minimal geometric deformation is given by In this case, the metric in Eq. (3.1) is minimally deformed by θ µν and its radial metric component becomes whereas the temporal metric component e ν remains unchanged. (More precisely, ν(r) becomes ν(r, α), as can be seen, for instance, in Eq. (5.16) below.) Upon replacing Eq. (3.7) in the Einstein equations (2.6)-(2.8), the system splits into two sets of equations: i) the first one is given by the standard Einstein equations for a perfect fluid (the one with α = 0 we started from), whose metric is given by Eq. (3.1) with ξ(r) = ν(r): along with the conservation equation (2.4) with α = 0, namely ∇ ν T (m)µν = 0, yielding which is a linear combination of Eqs (3.8)-(3.10). And ii) the second set contains the source θ µν and reads (3.14) The conservation equation (2.4) then yields to ∇ ν θ µν = 0, which explicitly reads Eqs. (3.11) and (3.15) simply mean that there is no exchange of energy-momentum between the perfect fluid and the source θ µν , so that their interaction is purely gravitational.
As was previously remarked in Ref. [1], Eqs. (3.12)-(3.15) look very similar to the standard spherically symmetric Einstein field equations for an anisotropic fluid with energy-momentum ten- The corresponding metric would be given by However, the right-hand sides of Eqs. (3.12) and (3.13) are not the standard expressions for the Einstein tensor components G 0 0 and G 1 1 , since there is a missing −1/r 2 in both. Nonetheless, the system of Eqs. (3.12)-(3.15) may be formally identified as Einstein equations for an anisotropic fluid with energy-momentum tensor θ * µν , whose effective energy densityρ, effective isotropic pressurep r and effective tangential pressurep t are given respectively bỹ which can be written more concisely as With the transformation (3.20), the conservation equation (3.15) takes the standard form and therefore the interpretation of Eqs. (3.12)-(3.15) as standard spherically symmetric Einstein equations for the source θ * µν in Eqs. (3.17)-(3.19) is now complete. As was pointed out in Ref. [1], since Eqs. (3.12) and (3.13) do not contain the standard Einstein tensor components, we should expect that the conservation equation (3.15) for the source θ µν is not longer a linear combination of Eqs. (3.12)-(3.14). However, Eq. (3.15) still remains a linear combination of the system (3.12)-(3.14). The MGD approach therefore turns the indefinite system (2.6)-(2.8) into the set of equations for a perfect fluid {ρ, p, ν, µ} plus a simpler system of four unknown functions {f * , θ 0 0 , θ 1 1 , θ 2 2 } satisfying the three Eqs. (3.12)- (3.14) [at this stage we suppose that we have already found a perfect fluid solution, thus ν is determined], or the equivalent anisotropic system of Eqs. (3.17)- (3.19). Either way, the system (2.6)-(2.8) has been successfully decoupled.
The MGD-decoupling can be summarised in a more formal way as follows: consider a static spherically symmetric gravitational source T µν containing one isotropic perfect fluidT µν and n other gravitational sources T then the diagonal metric g µν , solution of Einstein equation G µν = −k 2 T µν , reads

(3.24)
This metric g µν is found by first solving the Einstein equations for the sourceT µν , and then by solving the remaining n quasi-Einstein's equations for the sources T where the divergence-free quasi-Einstein tensorG µν is related with the standard G µν bỹ with Γ ν µ (g) a tensor that depends exclusively on g µν to ensure the divergence-free condition. In the spherically symmetric representation it reads The explicit components ofG ν µ in terms of the metric in Eq. (3.16) are shown in the right-hand side of Eqs. (3.12)-(3.14).

Matching condition for stellar distributions
A crucial aspect in the study of stellar distributions is the matching conditions at the star surface (r = R) between the interior (r < R) and the exterior (r > R) space-time geometries [48]. In our case, the interior stellar geometry is given by the MGD metric where the interior mass function is given bỹ with m given by the standard GR expression in Eq. (3. 2) and f * the yet to be determined MGD in Eq. (3.7). The inner metric (4.1) should now be matched with an outer geometry where there is no isotropic fluid, that is p + = ρ + = 0. The exterior r > R may however not be a vacuum anymore since we can, in general, have new fields coming from the sector described by θ µν . The general outer metric can be written as where the explicit form of the functions ν + and λ + are obtained by solving the effective 4dimensional exterior Einstein equations The MGD will reduce the exterior Einstein equations (4.4) to the system (3.12)-(3.14), where the geometric function ν is given by the Schwarzschild solution instead of a perfect fluid solution.
Continuity of the first fundamental form at the star surface Σ defined by r = R reads and where M 0 = m(R) and f * R is the minimal geometric deformation at the star surface. Likewise, continuity of the second fundamental form reads where r µ is a unit radial vector. Using Eq. (4.8) and the general Einstein equations (2.1), we then find which leads to p − α θ 1 1 Σ = 0 . This matching condition takes the final form where p R ≡ p − (R). The condition in Eq. (4.11) is the general expression for the second fundamental form associated with the Einstein equations given by Eq. (2.1). By using Eq. (3.13) for the interior geometry in the condition (4.11), the second fundamental form can be written as where ν ′ R ≡ ∂ r ν − | r=R . Moreover, using now Eq. (3.13) for the outer geometry in Eq. (4.12) yields where g * R is the geometric deformation for the outer Schwarzschild solution due to the source θ µν , that is (4.14) Eqs. (4.6), (4.7) and (4.13) are the necessary and sufficient conditions for the matching of the interior MGD metric (4.1) to a spherically symmetric outer "vacuum" described by the deformed Schwarzschild metric in Eq. (4.14). Note that this exterior could be filled by fields contained in the source θ µν .
The matching condition (4.13) yields an important result: if the outer geometry is given by the exact Schwarzschild metric, one must have g * (r) = 0 in Eq. (4.14), which then leads to the conditionp It is important to remark that the star will therefore be in equilibrium in a true (Schwarzschild) vacuum only if the effective (in general anisotropic radial) pressure at the surface vanishes. In particular, if the inner geometric deformation f * (r < R) is positive and weakens the gravitational field [see Eq. (4.2)], an outer Schwarzschild vacuum can only be compatible with a non-vanishing inner θ µν if the isotropic stellar matter has p R < 0 at the star surface. This may be interpreted as regular matter with a solid crust [12]. Otherwise, the condition p R = 0 can be obtained by imposing that the right-hand side of Eq. (4.15) be proportional to p R , namely, α (θ 1 1 ) − R ∼ p R in Eq. (4.11), which leads to a vanishing inner deformation f * R = 0.

Interior solutions
Let us now solve the Einstein's field Eqs. (2.6)-(2.8) for the interior of a self-gravitating anisotropic system by the MGD decoupling. The physical variables {ρ,p r ,p t } defined by Eqs. (2.10)-(2.12) and the two metric functions {ν, λ} in Eq. (2.5) will be derived. The first step is to turn off α and find a solution for the perfect fluid Einstein Eqs. (3.8)-(3.10). In particular, we can simply choose a known solution with physical relevance, like the well-known Tolman IV solution {ν, µ, ρ, p} for perfect fluids [49], namely e ν(r) = B 2 1 + r 2 A 2 , (5.1) and The constants A, B and C in Eqs. (5.1)-(5.4) are determined from the matching conditions in Eqs. (4.5) and (4.8) between the above interior solution and the exterior metric which we choose to be the Schwarzschild space-time. This yields with the compactness M 0 /R < 4/9, and M 0 = m(R) the total mass in Eq. (3.2). The expressions in Eq. (5.5) ensure the geometric continuity at r = R and will change when we add the source θ µν . Let us then turn α on in the interior. The temporal and radial metric components are given by Eqs. (5.1) and (3.7) respectively, where the interior deformation f * (r) and the source θ µν are related through Eqs. (3.12)-(3.14). Hence, we need to prescribe additional information in order to close the system of Eqs. (3.12)-(3.14). We have many alternatives for this purpose, like imposing an equation of state for the source θ µν or some physically motivated restriction on f * (r). In any case, we must be careful in keeping the physical acceptability of our solution, which is not a trivial issue. In the following, this problem is addressed in such a way that three new, exact and physically acceptable interior solutions will be generated.

Solution I: mimic constraint for pressure
From the matching condition (4.15) we see that the Schwarzschild exterior solution will be compatible with regular interior matter as long as α (θ 1 1 ) − R ∼ p R . The simplest choice satisfying this critical requirement is θ 1 1 (r) = p(r) , (5.6) which, according to Eq. (3.9), can be written as From Eq. (3.13) we can see that the "mimic" constraint in Eq. (5.6) is equivalent to .
Hence the radial metric component reads where the expressions in Eqs. (3.7) and (5.1) have been used. The interior metric functions given by Eqs. (5.1) and (5.9) represent the Tolman IV solution minimally deformed by the generic anisotropic source θ µν . We can see that the limit α → 0 in Eq. (5.9) leads to the standard Tolman IV solution for perfect fluids. This is represented by a generic way in Fig. 1. Now let us match our interior metric in Eq. (2.5) with metric functions (5.1) and (5.9) with the exterior Schwarzschild solution (4.14) with g * (r) = 0. We can see that, for a given distribution of mass M 0 and radius R, we have four unknown parameters, namely {A, B, C} from the interior solution in Eqs. (5.1) and (5.9), and the mass M in Eq. (4.14). Since we have only the three matching conditions (4.6), (4.7) and (4.11) at the star surface, the problem is not closed. An obvious solution would be to set B = 1 in Eq. (5.1), corresponding to the time rescaling t →t = B t. However, we want to keep B near its expression in the Tolman IV solution of Eqs. (5.1)-(5.4) in order to see clearly the effect of the anisotropic source θ µν on the perfect fluid. We will therefore solve for {A, C, M} with respect to B as shown further below.
The continuity of the first fundamental form given by Eqs. (4.6) and (4.7) leads to and whereas continuity of the second fundamental form in Eq. (4.11) yields On the other hand, by using the condition in Eq. (5.11), we obtain for the Schwarzschild mass where the expression in Eq. (3.2) has been used. Finally, by using the expression in Eq. (5.15) in the matching condition (5.10), we obtain Eqs. (5.14)-(5.16) are the necessary and sufficient conditions to match the interior solution with the exterior Schwarzschild space-time. By using the mimic constraint in Eq. (5.6) along with Eq. (5.14), the effective isotropic pressurep r in Eq. (2.11) reads pressure are given respectively bỹ p t (r, α) =p r (r, α) + 3 α r 2 k 2 (A 2 + 3 r 2 ) 2 . (5.20) We can further solve Eq. (5.16) for A = A(B, α), with B as a free parameter around the value in Eq. (5.5), in order to obtain the effective radial pressurep r (r, α) in Eq. (5.17) shown in Figure 2 for two values of α and B 2 = 2/5. It appears that the anisotropy produced by θ µν decreases the effective radial pressure more and more for increasing α.

Solution II: mimic constraint for density
An alternative choice leading to physically acceptable solutions is the "mimic constraint" for the density, which yields the first order differential equation The solution is given by where the density ρ in Eq. (5.3) has been used. To avoid a singular behaviour at the center r = 0 we must impose C 0 = 0, so that Eq. (3.7) yields By using the expression in Eq. (5.24), the matching conditions in (4.6) and (4.7) lead to and whereas continuity of the second fundamental form in Eq. (5.12) yields . .
It is also interesting to investigate the effects of the anisotropy on the surface redshift [50], which in our case reads 33) and is displayed in Figure 4. In particular, we can see from the plot that the anisotropy increases the gravitational redshift at the stellar surface, hence a distant observer will detect a more compact distribution in comparison with the isotropic case. This result is in agreement with the fact that M > M 0 , as can be seen from Eq. (5.25).

Solution III: extending anisotropic solutions
So far we have seen how to generate exact and physically acceptable anisotropic solutions starting from a known isotropic solution. In order to emphasise the full potential of the MGD-decoupling, schematically represented by Eqs. Let us start by adding to the energy-momentum tensor in Eq. (2.2) one more anisotropic source with energy-momentum tensor ψ µν , that is where β is a new coupling constant. The Einstein equations now read where g * (r) represents the geometric deformation. After decoupling, the source ψ µν and the un- We can see that the anisotropic effects increases the compactness of the self-gravitating system and therefore the surface redshift z.
known deformation g * (r) satisfy the same Eqs. (3.12)-(3.14), which now read We then impose the mimic constraint for pressure as and following the same procedure as in Section 5.1, we obtain Continuity of the first fundamental form in Eq. (4.5) again leads to the expression in Eq. (5.10), where the mass M is now given by p t (r, α, β) =p r (r, α, β) + Π(r, α, β) , (5.47) with the anisotropy (see Fig. 5) One can easily check that the solution II of Section 5.2 is the particular case β = 0 of the solution III given by Eqs. (5.1), (5.43) and (5.45)-(5.47). Likewise, the solution I of Section 5.1 (with α → β) is the particular case α = 0 of the solution III. Eventually, we could set α = β in the solution III so that the anisotropic sector is determined by a single free parameter α.

Conclusions
By using the MGD-decoupling approach, we presented in detail how to extend interior isotropic solutions for self-gravitating systems in order to include anisotropic (but still spherically symmetric) gravitational sources. For this purpose, we showed that the Einstein field equations for a static and spherically symmetric self-gravitating system in Eq. (2.6)-(2.8) can be decoupled in two sectors, namely: the isotropic sector corresponding to a perfect fluidT µν shown in Eq. (3.8)-(3.10), and the sector described by quasi-Einstein field equations associated with an anisotropic source θ µν shown in Eqs. (3.12)-(3.14). These two sectors must interact only gravitationally, without direct exchange of energy-momentum. The matching conditions at the stellar surface were then studied in detail for an outer Schwarzschild space-time. In particular, the continuity of the second fundamental form in Eq. (4.15) was shown to yield the important result that the effective radial pressurep R = 0. The effective pressure (2.11) contains both the isotropic pressure of the undeformed matter sourceT µν and the inner geometric deformation f * (r) induced by the energy-momentum θ µν . We recall the latter could also represent a specific matter source, like a Klein-Gordon scalar field or any other form of matter-energy, but also the induced effects of extra-dimensions in the brane-world. If the geometric deformation f * (r) is positive and therefore weakens the gravitational field [see Eq. (4.2)], an outer Schwarzschild vacuum can only be supported if the isotropic p R < 0 at the star surface. This can in fact be interpreted as regular matter with a solid crust [12] as long as the region with negative pressure does not extend too deep into the star.
In order to show the robustness of our approach, three new exact and physically acceptable interior anisotropic solutions to the Einstein field equations were generated from a single perfect fluid solution. All these new solutions inherit their physical acceptability from the original isotropic solution. In particular, it was shown that the anisotropic source θ µν always reduces the isotropic radial pressurep(r) inside the self-gravitating system.
We would like to remark that the MGD-decoupling is not just a technique for developing physically acceptable anisotropic solutions in GR, but represents a powerful and efficient way to exploit the gravitational decoupling in relevant physical problems. The extension of GR solutions into the domain of more complex gravitational sources is a highly non-trivial theoretical problem. For instance, it is well known that a Klein-Gordon scalar field induces anisotropic effects when it is coupled with the gravitational field through the Einstein equations. Hence the MGD-decoupling represents a useful tool for extending GR isotropic solutions for self-gravitating systems into solutions of the Einstein-Klein-Gordon system. It could be implemented, for instance, to investigate the role played by a Klein-Gordon scalar field during the gravitational collapse. In this respect, it is worth mentioning that the MGD-decoupling can be generalised for time-dependent scenarios, as long as the spherical symmetry is preserved under slowly evolving situations, which means the stellar system is always in hydrostatic equilibrium [46].