Compact Anisotropic Models in General Relativity by Gravitational Decoupling

Durgapal's fifth isotropic solution describing spherically symmetric and static matter distribution is extended to an anisotropic scenario. To do so we employ the gravitational decoupling through the minimal geometric deformation scheme. This approach allows to split Einstein's field equations in two simply set of equations, one corresponding to the isotropic sector and other to the anisotropic sector described by an extra gravitational source. The isotropic sector is solved by the Dugarpal's model and the anisotropic sector is solved once a suitable election on the minimal geometric deformation is imposes. The obtained model is representing some strange stars candidates and fulfill all the requirements in order to be a well behaved physical solution to the Einstein's field equations.


I. INTRODUCTION
In 1915 Albert Einstein [1][2][3][4] stunned the scientific community by presenting one of the greatest achievements of theoretical physics and at the same time one of the most beautiful theories known to date, which also has a great experimental support [5]. We are talking about General Relativity (GR from now on). Shortly after its publication, Schwarzschild [6] presented the first solution to the Einstein field equations. Solution that describes the neighborhood of a compact object that is spherically symmetric and static and free of material content, i.e. vanishing pressure and density. Later decades witnessed the great development that this wonderful theory experienced. Proof of this is reflected in the work of Professor Tolman [7], who found several analytical solutions which describe stellar interiors subject to a spherically symmetric and static geometry whose matter distribution corresponds to a perfect fluid (equal radial and tangential pressures p r = p t ). However, subsequent investigations showed that configurations with spherical symmetry do not necessarily meet the condition p r = p t . Of course, the introduction of anisotropies (unequal radial and tangential pressures, p r = p t ) allowed a better understanding about highly dense objects. The pioneering works by Ruderman [8] and Bowers and Liang [9] about anisotropic matter distribution opened the doors to one of the most researched branches nowadays. At present it remains a great challenge to find solutions to the Einstein field equations that also meet the requirements of physical admissibility [10]. Despite the great difficulty that it presents, an extensive variety of works available in the literature  (and reference contained therein) have successfully addressed the study and understanding of the role played by the anisotropy in stellar interiors. As Mak and Harko have argued [12], anisotropy can arise in different contexts such as: the existence of a solid core or by the presence of type 3A superfluid [54], pion condensation [55] or different kinds of phase transitions [56]. The presence of anisotropy introduces several features in the matter distribution, e.g. if we have a positive anisotropy factor ∆ = p t − p r > 0, the stellar system experiences an a repulsive force (attractive in the case of negative anisotropy factor) that counteracts the gravitational gradient. Hence it allows the construction of more compact objects when using anisotropic fluid than when using isotropic fluid [11][12][13]. Furthermore, a positive anisotropy factor enhances the stability and equilibrium of the system. On the other hand, matter tensor containing anisotropy, must be consistent with physical requirements for astrophysical applications [36].
Under the above background our motivation in the present paper is to extent isotropic solutions to the Einstein field equations to an anisotropic scenario. Specifically we have extended Durgapal's fifth model [57] to the anisotropic domain. To do so, we have applied a novel and systematic approach which opens up new possibilities for studies of compact stellar configurations that include anisotropic matter. This method was originally proposed in the context of the Randall-Sundrum braneworld [58,59] and was designed to deform the standard Schwarzschild solution [60,61]. Basically this scheme works decoupling gravitational sources through the so called minimal geometric deformation (MGD hereinafter) [62,63].
So, this work is organized as follows: Section II is devoted to the construction of the isotropic extension of the Durgapal's fifth model to an anisotropic scenario, in section III we analyze all the necessary requirements that an anisotropic solution of the Einstein field equations must meet to be physically admissible. Finally in section IV some conclusions are reported.
As was pointed out earlier, one of the essential features of anisotropic models is the inequality between radial and tangential pressures ı.e. p r = p t . To reach this condition in the present model, the starting point is to introduce an extra gravitational source which in principle can be e.g a scalar, vectorial or tensorial field. This extra source is coupled to the energymomentum tensor associated to the seed solution through a dimensionless coupling constant α. Explicitly it reads whereT µν corresponds to a perfect fluid given by (2)-(3). So, putting together the expressions (1) and (4) the Einstein field equations are The primes denote differentiation with respect to the radial coordinate r. From now on relativistic geometrized units are employed, that is c = G = 1. Bianchi's identities invokes the following conservation law being the above expression a linear combination of the equations (5)- (7). Moreover, ρ, p r and p t represent the effective density, the effective radial pressure and the effective tangential pressure respectively, that are given by The presence of the θ-term clearly introduces an anisotropy if θ r r = θ ϕ ϕ . Thus the effective anisotropy is defined as Solve the system of equations (5)- (7) is not an easy task. In order to tackle it we will employ the gravitational decoupling via the MGD approach [62]. This method consists in deforming the metric potentials e ν(r) and e λ(r) through a linear mapping given by where h(r) and f (r) are the corresponding deformations. It's worth mentioning that the foregoing deformations are purely radial functions, this feature ensures the spherical symmetry of the solution. The so called MGD corresponds to set h(r) = 0 or f (r) = 0, in this case the deformation will be done only on the radial component, remaining the temporal one unchanged (it corresponds to set h(r) = 0). Then the anisotropic sector θ µν is totally contained in the radial deformation (15). After replacing (15) into the system of equations (5)- (7), it is decoupled in two systems of equations. The first one corresponds to α = 0 it means perfect fluid matter distribution along with the conservation equatioñ it is a linear combination of the equations (16)- (18). The other set of equations corresponds to the θ sector The corresponding conservation equation ∇ ν θ µν = 0 then yields to this expression is a linear combination of the quasi-Einste-in equations. At this point it is remarkable to note that both the isotropic and the anisotropic sectors are individually conserved, it means that both systems interact only gravitationally.
From the equation (1) we have these expressions together with (2) and (3) The previous choice of the function f (r) is in accordance with the behavior of the metric potentials, that is: positive, regular and increasing monotone functions with increasing radius within the compact object. Figure 1 shows the increasing behaviour of the deformation function f (r) with increasing radius. Then from equations (20), (21) and (22) we obtain the following components for the source θ µν Hence, the anisotropic Durgapal's fifth model is described by the following metric potentials e ν(r) = A 1 + Cr 2 5 , (1 + 6Cr 2 ) 1/3 + αCr 2 1 + Cr 2 , (31) and thermodinamically characterized by whereρ(r) andp(r) are given by the expressions (2) and (3). From the latter equations, the anisotropy factor is directly computed, yielding to ∆ (α; r) = α 4π

A. Junction conditions
In order to find the arbitrary constants A, B and C we must match our interior solution where dΩ 2 ≡ sin 2 θdφ 2 + dθ 2 . For this purpose we will use the Israel-Darmois junction conditions [77,78]. These conditions require the continuity of the metric potentials e ν(r) and e λ(r) across the surface Σ of the compact object defined by r = R (It is known as the first fundamental form). Then we have The other condition is the second fundamental form [T µν r ν ] Σ = 0, where r ν is a unit vector projected in the radial direction. So the second fundamental form leads to regarding that the Schwarzschild exterior solution describes a vacuum space-time. The equation (39) provides us the following expression for the constant B The equations (37), (38) and (40)

III. PHYSICAL ANALYSIS
In this section we will analyze the basic requirements that all anisotropic solution to Einstein field equations must fulfill in order to be an admissible physical model describing stellar interiors [80].

A. Regularity
The departure point is to analyze if the model is free from physical and geometric singularities and non zero positive values of e ν(r) and e λ(r) i.e (e λ(r) )| r=0 = 1 and (e ν(r) )| r=0 > 0.
It is clear from expressions (30) and (31)   Another important aspect is related to the behaviour of the effective thermodynamic observables ρ(r), p r (r) and p t (r) within the stellar configuration. They must be positive and monotonically decreasing functions as they approach the to boundary of the compact object and their maximum values must be attain at the center r = 0. Evaluating (33) and (34) at then Moreover, to satisfy Zeldovich condition at the interior p r /ρ at center must be ≤ 1. Therefore In order to ensure the positiveness of p r and p t inside the distribution, from (41) the constant C must be positive. On the other hand (42) and (43) For the coupling constant α we have throughout the study we have fixed α = 0.3. The bound given for (45) leads to p t > p r ⇒ ∆ > 0 which represents a force due to the anisotropy directed outwards. Therefore we should have more massive and compact configurations [11,12]. We can observe from figures 5, 6 and 7 that p r , p t and ρ are monotonically decreasing functions with increasing radius and attain their maximum values at center of the star. Besides panels a, b and c in figure   8 shown how the effective radial pressure p r and the effective tangential p t pressure drift apart. The Zeldovich's condition is shown in figures 9 and 10. Clearly it is satisfied for the present model. Table II

B. Causality condition
An admissible anisotropic solution to Einstein field equations must satisfies causality condition ı.e both the radial v r and tangential v t sound speeds inside the object are less than the speed of light c (in relativistic geometrized units the speed of light becomes c = 1 ). Explicitly it reads v r (r) = dp r (r) dρ(r) ≤ 1 (46) v t (r) = dp t (r) dρ(r) ≤ 1. (47) As shown in figure 11 both speeds fulfill the above requirement. In addition, they have their maximum value at the center of the object (high density region) and decrease monotonically towards the surface (lower density region).

C. Energy conditions
Within the anisotropic matter distribution the energy should be positive. In order to ensure it, the energy-momentum tensor has to obey the null energy condition (NEC), the respectively.

D. Effective mass-radius ratio and redshift
In the study of spherically symmetric and static distributions associated with a perfect fluid, the maximum limit of the mass-radius ratio must satisfy the following upper bound u =M /R < 4/9 (in the units c = G = 1) [83]. However, in the presence of an anisotropic matter distribution this limit is more general [84]. However, it can be obtained from the effective mass defined as [12] M ef f = 4π Explicitly the effective mass M ef f reads Then the compactness factor u becomes So, the surface redshift (see Fig. 15) can be calculated using the compactness factor u given by (50), as follows The presence of a positive anisotropy factor ∆ > 0 does not impose an upper limit on the surface redshift Z s . In distinction with what happens in the case of isotropic distributions, where the maximum value that the surface redshift Z s can reaches is Z s = 4.77 [9]. Therefore, the surface redshift for anisotropic matter configurations is greater than its isotropic counterpart.

E. Equilibrium condition
The equilibrium of the system lies on the Tolman-Oppenheimer-Volkoff (TOV) equation [7,85]. In this case the equilibrium of the anisotropic fluid sphere is under three different forces. This forces are the hidrostatic force F h , the gravitational force F g and the anisotropic repulsive force F a introduced by the presence of a positive anisotropy factor ∆. In fact, the hidrostatic force F h and the anisotropic repulsive force F a counterbalance the gravitational force F g . Therefore, the collapse of the compact object to a point singularity may be avoided during the gravitational collapse. In conclusion, the presence of anisotropies within the stellar configuration enhance the equilibrium of the system.
Then the TOV equation describing the equilibrium condition for an anisotropic fluid distribution is given by We can observe from figures 16, 17 and 18 that the model is in equilibrium under the mentioned forces.

F. Stability conditions
Another relevant aspect in the study of anisotropic fluid spheres which is enhanced by the presence of anisotropies in the matter distribution is the stability. In order to verify if the present model is stable we analyze both the relativistic adiabatic index Γ [86,87] and the square of sound speeds behaviour inside the compact object [88].
It is well known from the studies about Newtonian isotropic fluid spheres that the collapsing condition correspond to Γ < 4/3. In distinction with the relativistic counterpart this condition becomes [89,90] where ρ 0 , p r0 and p t0 are the initial density, radial and tangential pressure when the fluid is in static equilibrium. The second term in the right hand side represents the relativistic corrections to the Newtonian perfect fluid and the third term is the contribution due to anisotropy. It is clear from (53) that in the case of a non-relativistic matter distribution and taking p r be equal to p t ı.e ∆ = 0, the bracket vanishes and we recast the collapsing Newtonian limit Γ < 4/3. On the other hand, Heintzmann and Hillebrandt [87] showed that in the presence of a positive an increasing anisotropy factor ∆ = p t − p r > 0, the stability condition for a relativistic compact object is given by Γ > 4/3, this is so because positive anisotropy factor may slow down the growth of instability. Due to the fact that gravitational collapse occurs in the radial direction, it is enough to analyze the adiabatic index in such direction. So, the adiabatic index is given by [91] Γ r = ρ + p r p r dp r dρ .
We can see from figure 19 that the model is in complete agreement with the condition Γ r > 4/3. Then the model is stable.
Based on the Herrera's cracking concept [92] Abreu et. al. [88] established another alternative to study the stability of a self-gravitating anisotropic fluid sphere. This approach states that the region is potentially stable where the radial speed v r of sound is greater than the transverse v t speed of sound. This implies that there is no change in sign v 2 t − v 2 r . The later assumption is equivalent to 0 ≤ |v 2 t − v 2 r | < 1. We note from figure 11 that radial speed of sound is always greater than transverse speed of sound and also from figure 20 0 ≤ |v 2 t − v 2 r | < 1 everywhere inside the star. On the other hand, figure 21 shows that there in no change in sign. These features represent that the proposed physical model is stable.