Regularity condition on the anisotropy induced by gravitational decoupling in the framework of MGD

We use gravitational decoupling to establish a connection between the minimal geometric deformation approach and the standard method for obtaining anisotropic fluid solutions. Motivated by the relations that appear in the framework of minimal geometric deformation, we give an anisotropy factor that allows us to solve the quasi–Einstein equations associated to the decoupler sector. We illustrate this by building an anisotropic extension of the well known Tolman IV solution, providing in this way an exact and physically acceptable solution that represents the behavior of compact objects. We show that, in this way, it is not necessary to use the usual mimic constraint conditions. Our solution is free from physical and geometrical singularities, as expected. We have presented the main physical characteristics of our solution both analytically and graphically and verified the viability of the solution obtained by studying the usual criteria of physical acceptability.


Introduction
In 1916, Karl Schwarzschild obtained the first interior solution of the Einstein field equations [1]. This solutions describe a self-gravitating object sustained by a perfect and incompressible fluid which is embedded in a static and spherically symmetric vacuum space-time. Following the strategy of Schwarzschild, other interior solutions can be constructed providing suitable equations of state to close the system. However, in some cases the system obtained can not be analytically integrated and numerical models are required. Besides proposing an equation of state to relate thermodynamical quantities, we can use geometrical constraints on a e-mail: econtreras@usfq.edu.ec the functions. Indeed, following this program Tolman obtained a family of eight isotropic solutions [2].
For many years isotropic solutions were considered as well posed models to study stellar interiors. However, as was shown by Delgaty and Lake [3], very few of this solutions can be considered physically acceptable. More precisely, among 127 isotropic solutions they considered, only 16 passed all the test that any realistic solution must satisfy (for a list of physical conditions of interior solutions see, for example, [4]). Nevertheless, even when acceptable solutions can be found the perfect fluid model is inadequate when local anisotropy of pressure is assumed. Regardingly, anisotropic models have been considered as very reasonable for describing the matter distribution under a variety of circumstances [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. Now, as it is well known, assumption of local anisotropy in the fluid leads to the introduction of an extra unknown quantity in the system. In this sense, we need to imposse two conditions, either equations of state or geometric links between the metric variables in order to integrate the system. For example, in Ref. [6] Bowers and Liang besides assuming the Schwarzschild constraint on the density, they impossed that in order to avoid singularities in the Tolman-Oppenheimer-Volkoff (TOV) equation the anisitropy, p ⊥ − p r must satisfy, for example, the following constraint In the above expression ρ is the energy density of the fluid and p ⊥ and p r stand for the transverse and radial pressure of the fluid, respectively. The function f encodes the information of the anisotropy of the system which is no necessarily a linear function of the ra-arXiv:2001.08573v1 [gr-qc] 22 Jan 2020 dial pressure and C a parameter which measures the anisotropy strength. Finally, the exponent is constrained to n > 1 in order to aviod singularities in (1). In the same spirit, Cosenza et. al. [7] proposed a method which allowed them to find a family of non-isotropic models from any isotropic model which depends continuously on the constant C. The protocol consists in to take the energy density of any perfect fluid as the density of the anisotropic system and consider the following anisotropic function, f (p r , r) Following this procedure they were able to extend the Schwarzshild interior solution, Tolman IV, V and VI, and the Adler model to anisotropic domains. However, it is worth mentioning that in some cases numerical analysis were required to obtain the solution. The procedure just explained corresponds to the standard way in which problems have been solved in the presence of anisotropy.
Recently, the so-called Minimal Geometric Deformation (MGD) method [69][70][71][72] has emerged as an alternative to extend isotropic solutions in a straightforward and analytical way given the number of ingredients which convert it in a versatile and powerful tool to solve the Einstein's equations.
For example, the method has been used to obtain anisotropic like-Tolman IV solutions [35,39], anisotropic Tolman VII solutions [67] and a model for neutron stars [68]. In other contexts, MGD has been used to extend black holes in 3 + 1 and 2 + 1 dimensional space-times [42,53,56]. Moreover, in the context of modified theories of gravitation, the method has been used to obtain solutions in f (G) gravity [46], Lovelock [64], f (R, T ) [62], Rastall [71] and more recently interior solutions in the context of braneworld [70].
It is worth mentioning that in contrast to the standard strategy followed in Ref. [7] where the information of the isotropic solution entered via the energy density of a well known model, in the MGD method the isotropic solution is a sector of the total solution. More precisely, the isotropic solution is used as a seed to obtain anisotropic solutions of the Einstein equations as follows.
Let us consider the Einstein field equations and assume that the total energy-momentum tensor, T µν , can be decomposed as where T (m) µν is the matter energy momentum for a perfect fluid and θ µν an anisotropic source interacting with T (m) µν . Note that, since the Einstein tensor is divergence free, the total energy momentum tensor T (tot) µν satisfies ∇ µ T (tot)µν = 0. (6) It is important to point out that, as this equation is fulfilled and given that for a perfect fluid we also have ∇ µ T (m)µν = 0, then the following condition necessarily must be satisfied In this sense, there is no exchange of energy-momentum tensor between the perfect fluid and the anisotropic source and henceforth interaction is purely gravitational.
In what follows, we shall consider a static, spherically symmetric space-time with line element parameterized as where ν and λ are functions of the radial coordinate r only. Now, considering Eq. (8) as a solution of the Einstein equations, we obtain where the primes denote derivation with respect to the radial coordinate and we have defined Note that, at this point, the decomposition (5) seems as a simple separation of the constituents of the matter sector. Even more, given the non-linearity of Einstein's equations, such a decomposition does not lead to a decoupling of two set of equations, one for each source involved. However, contrary to the broadly belief, the decoupling is possible in the context of MGD. The method consists in to introduce a geometric deformation in the the metric functions given by where {g, f } are the so-called decoupling functions and α is a free parameter that "controls" the deformation.
It is worth mentioning that although a general treatment considering deformation in both components of the metric is possible (see Ref. [55]), in this work we shall concentrate in the particular case g = 0 and f = 0.
Doing so, we obtain two sets of differential equations: one describing an isotropic system sourced by the conserved energy-momentum tensor of a perfect fluid T (m) µν and the other set corresponding to quasi-Einstein field equations sourced by θ µν . More precisely, we obtain with for the perfect fluid and for the source θ µν that, whenever θ 1 1 = θ 2 2 , induce local anisotropy in the system as can be seen in Eqs. (13) and (14). It is worth noticing that the conservation equation which is a linear combination of Eqs. (21), (22) and (23). Note that unlike quasi-Einstein equations, which differ from the Einstein equations, this equation is completely analogous to an anisotropic TOV equation as can be seen in Ref. [7]. Now, given metric functions {ν, µ} sourced by a perfect fluid {ρ, p} that solve Eqs. (17), (18) and (19), the deformation function f can be found from Eqs. (21), (22) and (23) after choosing suitable conditions on the anisotropic source θ µν . It is worth mentioning that the case we are dealing with demands for an exterior Schwarzschild solution. In this case, the matching condition leads to the extra information required to completely solve the system.
Defining µ(r) = 1 − 2m(r) r in (16), the interior solution parameterized with (8) reads Now, outside of the distribution the space-time is that of Schwarzschild, given by In order to match smoothly the two metrics above on the boundary surface Σ, we must require the continuity of the first and the second fundamental form across that surface. Then it follows Note that, the condition on the radial pressure leads to Regardingly, if the original perfect fluid match smoothly with the Schwarzschild solution, i.e, p(r Σ ) = 0, Eq. (30) can be satisfied by demanding θ 1 1 ∼ p. Of course, the simpler way to satisfy the requirement on the radial pressure is assuming the so-called mimic constraint [35] for the pressure, namely in the interior of the star. Remarkably, this condition leads to an algebraic equation for f such that, in principle, any isotropic solution can be extended with this constraint. Another possibility is to use the mimic constraint for the density which leads to a differental equation for f which can be solved in some situations (see for example [39]). However, as far as we know, no physical requirements on the decoupler sector have been imposed up to now. More precisely, the anisotropic solution is obtained using information of the isotropic sector. In this respect, in this work we find an anisotropic solution assuming a regularity condition in the decoupler sector.
To be more precise, we impose the Bowers-Liang constraint, given by Eq.
In this sense, we propose the following condition on the decoupler sector reads, with f (θ 1 1 , r)r n−1 = ν /2 and C, as usual, is a constant that gauge the anisotropy strength. This ansatz is inspired by the relation between the components of the anisotropic energy-momentum tensor θ µν and the effective quantities given by Eqs. (12), (13) and (14).
Note that the function f in Eq. (32) is not the deformation function that appears in Eq. (16). It is clear that the replacement of (21), (22) and (23) in (32), leads to a differential equation for the deformation function, f , where the only required information is the metric function ν, which in the context of MGD is common for the three sectors involved. This work is organized as follows. In the next section we study the regularity condition on the decoupler sector induced by MGD. In section 3, we study the conditions for physical viability in interior solutions. Finally, the last section is devoted to final remarks.
2 Regularity condition on the decoupler sector As previously commented, regularity condition on the anisotropy of the system is an obligated requirement in the construction of well posed anisotropic solutions. In the context of MGD the anisotropy is induced by the decoupler sector sourced by θ µν which satisfy a conservation equation given by (24). Now, from Eq. (32) and after impossing the Consenza-Herrera-Esculpi-Witten anisotropy we obtain which can be formally solved to obtain where c 1 is a constant of integration. Of course, finding a solution of the above integral will depend on the particular form of the metric function ν.
As a particular case of application we shall consider the Tolman IV solution given by where A, B and d are constants. Next, replacing (35) in (34) we obtain which determines the decoupler sector completely and enusure the regularity of the anisitropy θ 2 2 −θ 1 1 . To complete the MGD program, the rest of the section is devoted to obtain the total like-Tolman IV anisitropic solution. From Eq. (16), the g rr component of the metric reads e −λ = A 2 + r 2 αc 1 r 2 A 2 + 2(C + 1)r 2 + Now, from (9), (10) and (11) we obtain the effective quantities In order to match the interior solution with the Schwarzschild exterior solution, we proceed to impose the continuity of the first and the second fundamental (see Eqs. (27), (28) and (29)) from where In this sense, the solution is parameterized by the mass M , the radius R, the parameter of anisotropy strength C, the MGD parameter α, and the constant of integration c 1 . We would like to point out that, unlike what usually occurs, in our solution we have not made use of any mimic constraint condition. To be more precise, we can see that imposing a physical condition on the decoupling sector (in this case to give the anisotropy factor) allows us to avoid the use of the condition.
It is worth mentioning that, although in our previous discussion, the physical condition was introduced on the decoupling sector at the level of quasi-Einstein equations, it could be also be introduced at the level of the conservation equation for the decoupling sector (24). This may open up new ways of approaching problems in MGD formalism.
In the next section we shall study the acceptability conditions for the like-Tolman IV anistropic solution obtained here.

Conditions for physical viability of interior solutions
The study of acceptability conditions of interior solutions is important because, as it is well known, in addition to solving Einstein's equations, we have to demonstrate that they are adequate to describe a real physical system. In this section we explore some of the conditions used for this purpose. To this end, we perform the physical analysis of the properties of the star solution by fixing the free parameters of the solution and providing plots. More precisely, in order to obtain a useful model for a compact anisotropic star we specify the mass M and the radius R of the star and impose some suitable conditions that the model should satisfy. The following conditions have been typically recognized as decisive for anisotropic fluid spheres.

Matter sector
A requirement on the matter sector to ensure acceptable interior solutions is that the density and pressures should be positive quantities. Besides, it also demanded that the density and pressures reach a maximum at the center and decrease monotonously toward the surface so thatp ⊥ ≥p r . In figures 1, 2 and 3 we show the behaviour ofρ,p r andp ⊥ respectively. It is worth noticing that the anisotropy factor C has an appreciable effect on the plots in the sense that the separation of each graphic parametrized with α increases as C decreases. To be more precise, in all the cases the the profiles in panel a) are almost indistinguishable in contrast to panel c) where the behaviour for different α is appreciable different.
To study the extra conditionp ⊥ >p r , in fig. 4 it is shown the anisotropy ∆ =p ⊥ −p r . Note that the anisotropy ∆ is a positive and increasing function, as expected.

Energy conditions
Another physical requirement we demand for interior solutions involve the energy conditions. As it is well known, an acceptable interior stellar should satisfy the the dominant energy condition (DEC), which implies that the speed of energy flow of matter is less than the speed of light for any observer. This condition reads In figure 5, it is shown that DEC is fulfilled by all the parameters considered here. A stronger condition is that the solution satisfies the strong energy condition (SEC) also, namelỹ As can be seen in fig. 6, the SEC is satisfied in the cases under consideration.

Causality
Causality is important to avoid superluminal motion. In other words, the causality condition demands that either the radial and tangential sound velocities, v r = dp r /dρ and v t = dp ⊥ /dρ respectively, are less than the speed of light. Given the behaviour of the radial and the traverse velocities illustrated in figures 7 and 8, we conclude that our model satisfy the causality condition requirement.

Adiabatic index
The adiabatic index, γ, serves as a criterion of stability of the interior solution. It can be shown that for anisotropic fluids the adiabatic index takes the form γ =ρ +p r p r dp r dρ , It is said that an interior configuration is stable whenever γ ≥ 4/3. In figure 9 we show the adiabatic index for different values of the free parameters involved. It is clear that the solution is stable regarding the adiabatic index criterion.

Stability against gravitational cracking
The appearance of non-vanishing total radial force with different signs in different regions of the fluid is a sign of instability. When this radial force, running from the center to the outside of the star, shifts from pointing to the center to pointing outward, the phenomenon has been called gravitational cracking [76]. In reference [77] it is stated that a simple requirement to avoid gravita- tional cracking is In figure 10 we show that in all the cases considered the solution is stable against gravitational cracking.

Final remarks
The minimal geometric deformation method has proven to be a simple and powerful tool for obtaining solutions of Einstein's field equations. The model studied in this article describing anisotropic fluid spheres meets all the requirements to be an acceptable solution both from a physical and mathematical perspective. In this work we established a connection between standard approaches to obtain anisotropic stellar solutions and the minimal geometric deformation method. The standard approach usually provides the anisotropy factor ∆, so we incorporate this information in the MGD method and in this way, the use of the mimic constraint condition becomes unnecessary.
Using the analytical Tolman IV perfect fluid solution in the MGD approach we get a new solution that represents the anisotropic extension of Tolman IV solution. This new analytical solution satisfies all the usual criteria of physical acceptability. We have evaluated the physical consistency of our anisotropic star structure by examining the structure of matter sector, energy conditions, causality, the adiabatic index and stability against gravitational cracking. Therefore, this could be used to model actual stellar compact structures, such as neutron stars.
As a continuation of the study presented here, we suggest to explore of anisotropic solutions using the conservation equation obtained from the decoupling sector. This and other aspects will be considered in future works.