Gravitational decoupled anisotropies in compact stars

Simple generic extensions of isotropic Durgapal--Fuloria stars to the anisotropic domain are presented. These anisotropic solutions are obtained by guided minimal deformations over a self gravitating isotropic system. When the isotropic and the anisotropic sector interacts in a purely gravitational manner, the conditions to decouple both sectors by means of the minimal geometric deformation approach are satisfied. Hence the anisotropic field equations are isolated resulting a more treatable set. The simplicity of the equations allows one to manipulate the anisotropies that can be implemented in a systematic way to obtain different realistic models for anisotropic configurations. Later on, observational effects of such anisotropies when measuring the redshift are discussed. To conclude, the application of the method over anisotropic solutions is generalized. In this manner, different anisotropic sectors can be isolated of each other and modeled in a simple and systematic way. Besides, a generic property of the minimal geometric deformation approach, its noncommutativity, is discussed. This property duplicates the solutions obtained through this approach; anisotropies applied in the reversed order give different physically acceptable configurations.


Introduction
The study of analytical solutions of Einstein field equations plays a crucial role in the discovery and understanding of relativistic phenomena. Theoretical grounds gives very few isotropic solutions under static and spherically symmetric assumptions. Worse yet, many less of these solutions have physical relevance passing elementary tests of astrophysical observations [1][2][3][4]. Furthermore, no astronomical object is constituted of perfect fluids only; hence isotropic approximation is likely to fail.
Anisotropic configurations have continuously attracted interest and are still an active field of research. Strong evidence suggests that a variety of very interesting physical phenomena gives rise to a quite large number of local anisotropies, either for low or high density regimes (see [5] and references therein). For instance, among high density regimes, there are highly compact astrophysical objects with core densities ever higher than nuclear density (∼ 3 × 10 17 kg/m 3 ) that may exhibit an anisotropic behaviour [6]. Certainly, the nuclear interactions of these objects must be treated relativistically. The anisotropic behaviour is produced when the standard pressure is split in two different contributions: i) the radial pressure p r and ii) the transverse pressure p t , which are not likely to coincide.
Anisotropies in fluid pressure usually arise due to the presence of a mixture of fluids of different types, rotation, viscosity, the existence of a solid core, the presence of a superfluid or a magnetic field [7]. Even are produced by some kind of phase transitions or pion condensation among others [8,9]. The sources of anisotropies have been widely studied in the literature, particularly for different highly compact astrophysical objects such as compact stars or black holes, either in 4 dimensions [10,11] as well as in the context of braneworld solution in higher dimensions [12][13][14].
The main purpose of the present article is to generalize anisotropic analogous solutions of a particular kind of isotropic compact objects by means of the so-called minimal geometric deformation approach (MGD hereinafter) [15,16]. This method was originally proposed in the context of the Randall-Sundrum braneworld [17,18] and was designed to deform the standard Schwarzschild solution [19,20]. It describes the 4D geometry of a brane stellar distribution, hence obtaining braneworld corrections to standard GR solutions. Therefore, it is a suitable method to obtain spherically symmetric and inhomogeneous stellar distributions that are physically admissible in the braneworld. The key point of this approach is that the isotropic and anisotropic sectors can be split. Thus, the decoupling of both gravitational sources can be done in a simple and systematic way establishing a new window to search for new families of anisotropic solutions of Einstein field equations.
In this work we applied a gravitational decoupling through the MGD approach to derive exact and physically acceptable anisotropic interior solutions analogous to the Durgapal and Fuloria superdense star [21]. There have been several proposals of anisotropic models analogous to Durgapal-Fuloria compact stars [22,23]; the MGD method seems to generalize them. The details of this method will be shown later, however the main lines goes as follows: Let us start with a well known spherically symmetric gravitational source T (0) µν . This source can be as simple as one would desire; one can start with any known perfect fluid or even with vacuum itself. Any classical solution works as a seed for this method. After this, one switch on a new source of anisotropy When gravitational sources are couple via gravity only, i.e. they do not exchange energy-momentum among each other, the set of equations can be split into two contributions. On the one hand, a well known sector is identified with the classical field equations of the chosen seed; the Durgapal-Fuloria solution for compact stars in our case. On the other hand, one is left with a simpler set of 'pseudo-Einstein' equations for the sources of the anisotropy, to be solved. Combining both sectors a full anisotropic and physically consistent solution of Einstein field equations is obtained. Of course one can switch on as many arbitrary sources of anisotropies T (i) µν as desired, as long as a strategy to solve the new sector can be found.
This method for decoupling non-linear differential equations can be applied in a systematic way and has a vast unexplored territory where it could give different novel perspectives. MGD does not only give consistent interior solutions for different isotropic perfect fluid in GR; it could also be conveniently exploited in relevant theories such as f (R)-gravity [24,25], intrinsically anisotropic theories as Hořava-aether gravity [26] or to study the stability of novel proposals of Black Holes, described by Bose Einstein gravitational condensate systems of gravitons [27][28][29]. This is a robust method to extend physical solutions into an anisotropic domain preserving the physical acceptability.
The paper is organized as follows: after this introduction, we present the Einstein field equations for an anisotropic fluid. In Section 3 we explain how the MGD approach is implemented to generate arbitrary anisotropic solutions. Section 4 is dedicated to apply this method to a particular seed, the Durgapal-Fuloria model for compact stars. In Section 5 we extend the method to seeds which are already anisotropic. The last two sections are dedicated to discuss the main results and summarize our conclusions.

Anisotropic effective field equations
The simplest approach to describe compact distributions modelling stellar structures, is to restrict the metric to be static and spherically symmetric. In the usual Schwarzschild-like coordinates the line element takes the standard form where the functions ν ≡ ν(r) and λ ≡ λ(r) depend on the radial coordinate only. The encoded metric is a generic solution of the Einstein field equations describing an anisotropic fluid sphere. The coupling constant between matter is given by κ = 8πG/c 4 . Along these lines we will work in relativistic geometrized units, G = c = 1. The observable features of the object will be determined by the exterior metric that will describe the geometry of the outer part. In the present article to maintain the treatment as simpler as possible, we will suppose a Schwarzschild vacuum outside. The corresponding anisotropic effective stress-energy tensor T µν is characterized by its diagonal components ρ, p r and p t , that are related to the geometric functions µ, ν through (3). Explicitly, The prime stand for derivatives w.r.t. r. There is another equation consequence of the Bianchi identities: the covariant conservation of the stress-energy tensor Since the discovery of the first interior stellar solution by Schwarzschild [30] and for several years, stars interior were supposed to be constituted by perfect fluids. It was not until 1933 when Lemaître [31] develop that spherically symmetry do not require the isotropic condition p r = p t , but only the equality of the two tangential pressures p θ = p φ = p t . The system of equations (4)-(7) governs the matter distribution within the star, which is assumed to be locally anisotropic (the radial and tangential pressure do not coincide). It is necessary to solve for five unknowns functions: two geometric functions, ν(r) and λ(r); and three effective scalar functions, ρ(r), p r (r) and p t (r).
However there are more unknowns than equations, hence the system is undetermined and constrains must be imposed. Some of them must be chosen by consistency of regularity, stability and (or) energy conditions of relativistic models; see for instance [5,7,32].
Throughout this article we will make use of the following representation for the effective energymomentum tensor The first term encodes a perfect fluid with isotropic pressure p = p r = p t , u µ is the normalized four-velocity field that accomplish u µ u ν g µν = 1. In our case, the perfect fluid will invariably be given by the Durgapal-Fuloria interior solution. Under this representation, the anisotropic sector is described by the θ-term. It describes additional gravitational sources responsible of the anisotropies. These source may contain new fields, whether scalar, vector or (and) tensor fields, coupled to gravity by means of a free dimensionless and constant parameter α. One of the simplest and most treated examples in the literature are the anisotropies that may arise due to extra interactions resulting from the presence of charge [33]; besides there are plenty of complex treatments of anisotropies generated by other sophisticated physical fields [34]. The effective stress-energy tensor (8) contributes at the level of Einstein equations with an effective energy density ρ, an effective radial pressure p r and an effective tangential pressure p t defined as Thus, each magnitude is written as a deviation from the GR solution due to the presence of the θ-term. The additive structure for the anisotropies allows the theory to have a straightforward limit to GR; setting α = 0 the standard Einstein equations for the perfect fluid are recovered.
Since the Einstein tensor is divergence free, under the representation taken in (8) the covariant conservation equation (7) yields This equation is a linear combination of (4) and (6), as commonly happens in perfect fluid solutions of Einstein equations. As this point let us remark the appearance of the anisotropy: there is not an a priori restriction for the components of θ µν ; however, if θ r r = θ ϕ ϕ when solving the equation system (4)-(6), we will be in the presence of the pressure anisotropy Therefore, an isotropic stellar distribution (perfect fluid) becomes anisotropic when the θ-term is turned on. In these lines we will follow a different approach to tackle the equation system (4)-(6); we will address this system by means of the MGD method. This theory decouples the Einstein field equations when deforming the metric of the corresponding GR solution [15,16,35,36].

Minimal geometric deformation approach
With the aim of approaching the system of equations (4)- (6) in an alternative manner, a briefly review on the MGD procedure will be presented. This method produces anisotropic corrections to standard GR solutions providing physically admissible non-uniform and spherically symmetric stellar distributions. The input (seed) is a known solution of Einstein equations: for instance the thermodynamic parameters satisfying (9), and the corresponding geometric functions λ(r) and ν(r). When a perfect fluid solution is taken as a seed, the isotropic condition p r = p t = p is automatically accomplished. The method will produce a drift in the effective pressures such that p r = p t . For doing so, one implements the most generic minimal geometric deformation over the metric without breaking the spherical symmetry of the initial solution; this is with e and f generic functions parametrizing the metric deformation. In Figure 1 a schematic picture exemplifies how this method extends GR solutions to anisotropic domains when releasing α. Even though the theory does not impose limits for the coupling strength, the physical acceptability of the new solution does so; if α is increased, the anisotropies become at some point unstable.
Although nothing prevent us from deforming the temporal component of the metric, it is general enough to start setting e * = 0; hence the effects of the θ µν source undergo in a deformation over the radial coordinate only. The peculiarity of the MGD method is that it entails in its formulation a decoupling of the equations of motion. As a consequence of taking the θ-sector as responsible of the minimal distortion of the metric, the system of equations (4)-(6) results quasi-decoupled: we obtain the Einstein equations for the chosen perfect fluid; and an effective 'pseudo-Einstein' system of equations governing the θ-sector. The only parameter that connects the two sectors is the temporal geometric function ν(r). At the same order as before, the temporal, radial and angular equations of motion relating the geometry of the spacetime to the thermodynamic characteristic of Figure 1: Any perfect fluid solution can be extended via the MGD approach [16]. When α = 0 we are in the space of isotropic solution of Einstein equations, delimited by the red dashed line. This solution has a smooth extension to the anisotropic domain releasing the α to be nonnull. the perfect fluid sector reduce to The definition of a perfect fluid entails in itself the covariant conservation of the stress-energy tensor, i.e.
The resulting equation is again a linear combination of the temporal and angular equations, (17) and (19), and yields It is worth noting that this system of equations is equivalent to Eqs. (4)- (6) if the coupling between the two sectors is set to zero; this is if the anisotropic sector vanishes. The temporal component of the metric must satisfy binding conditions in the anisotropic sector: these are the remaining 'pseudo-Einstein' field equations for the θ-sector Once again, one has the corresponding conservation equation that is a consequence of (7) and (20) satisfying separately. This equation is and it is explicitly written as This time the latter equation is not necessary linearly dependent of the 'pseudo-Einstein' equations, and there is no reason why it should be. At this point it makes explicit that the interaction between the two sectors is purely gravitational; this is, from (20) and (25) is clear that each sector is separately conserved and there is no exchange of energy-momentum between them.
To conclude this section, let us summarize. First we started with an indeterminate system of equations (4)- (6). Then, we performed a linear mapping of the radial geometric function of the metric (16) that results in a 'decoupling' of the Einstein field equations. We ended with two sets of equations: a perfect fluid sector {ρ ; p ; ν ; µ}, given by (17)-(21) where everything is known once a perfect solution of GR is chosen; and a simpler sector of three linearly independent equations that can be chosen from (22) to (26), for determining four unknown functions {f * ; θ t t ; θ r r ; θ ϕ ϕ }. Once the second sector is solved, we can identify directly the effective physical quantities introduced in (10), (11) and (12). At this point, is mandatory to recall that the underlying anisotropic effect which appears as a consequence of breaking the isotropic condition over the effective pressures, p t = p r , causes the appearance of the anisotropy Π(α; r) defined in Eq. (14).

Anisotropic Durgapal-Fuloria compact star
Let us proceed now to apply the MGD method with the aim of solving the Einstein field equations for the interior of anisotropic superdense stars. In the present work we will take as a seed the wellknown Durgapal-Fuloria solution {ν ; µ ; ρ ; p} modeling compact stars. As explained before, once the MGD method is applied the system of equations (4)- (6) is decoupled. Half of the decoupled equations (17)- (19) are already solved once the relativistic perfect fluid is chosen. For instance, the thermodynamic functions that characterize the Durgapal-Fuloria solution are with C an integration constant. The gravitational mass of a sphere of radius r is obtained integrating the density inside the corresponding volume; in spherical coordinates is This mass function has a well defined behaviour, vanishing at the center of the compact object, i.e. m(r = 0) = 0. It also determines the total mass evaluating the mass function at the surface, m(r = R) = M . A massive object deforms the surrounding spacetime; the Durgapal-Fuloria solution is defined by the following metric components e ν(r) = A 1 + Cr 2 4 , It is a custom in GR to write the radial component of the metric with the so-called compactness parameter, given by ξ = 2m/r. The spacetime results regular everywhere, even at the center where e λ(r=0) ≡ µ(r = 0) = 1; m vanishes faster than r as one can easily check from (29) inside (31). A is the second (and last) integration constant to be determined using boundary conditions over the surface r = R. In the present article the outer metric will be chosen to satisfy the Schwarzschild form-for simplicity, an uncharged compact star. Both constants A and C are positive; however they are expected to change as far as anisotropies begin to be considered. The remaining equations after the decoupling, (22) to (24), have to be solved if a generic anisotropic self-gravitating system is desired. The system of equation is as explain before underdetermined. A reasonable constrain is needed to close the system, but it is mandatory not to lose the physical acceptability of the solution. These issues will be discussed in what follows when three different anisotropic solutions (of many) are presented.

Pressure-like constraint for the anisotropy
In order to close the system of equations (22)-(24), additional information is needed. For instance, an equation of state for the source θ µν or some physically motivated constrain on f * (r). A first acceptable interior solution is deduced when forcing the associated radial pressure θ r r to mimic a physically acceptable pressure θ r r (r) = p(r) .
This means that one simple choice is to require that the stress-energy tensor for the perfect fluid coincides with the anisotropy in that direction. As a consequence of (32), the radial Einstein equations for the GR solution (18) and the radial 'pseudo-Einstein' equation (23) are equal. This gives immediately an expression for the radial component metric deformation The temporal component of the metric (30) remains non-deformed, so ν is computed directly. The resulting deformed component, the radial one in (16), then becomes It is explicit that when the α → 0 limit is taken, one gets the non-perturbed Durgapal-Fuloria solution; particularly for the radial component of the metric, e −λ(r) = µ(r).
With the above considerations, the metric can be written in terms of an effective mass function of the anisotropic sphere given by Expressed in this form, it is obtain one branch of MGD metrics that govern anisotropic interiors of GR solutions, whatever the GR solution is chosen. This branch corresponds to the pressure constrain imposed over the radial anisotropy. Therefore, the metric (2) is deformed to As we have closed the system of equations with the constrain (32), we can compute all the effective magnitudes that characterized the fluid; but first, the values of the integration constants A and C are needed to be fixed. This will be done by means of consistent matching conditions.

Matching conditions
A crucial aspect in the study of stellar distributions is the matching conditions at the star surface between the interior and the exterior spacetime geometries [37,38]. In our case, the interior stellar geometry is given by the MGD metric (36); while the outer part is assumed to be empty. Hence for r ≥ R the solution is given by the Schwarzschild vacuum solution. The continuity of the first fundamental form at the star surface Σ (defined by r = R) is given by ds 2 Σ = 0. This equation implies the continuity of the metric when crossing the surface and reads for the relevant components (tt and rr-components) as The superindices stand for the region from where we approach the surface, either from inside with a minus sign, or from outside using the plus sign. We must also take into account the Israel-Darmois matching condition at the stellar surface Σ that gives the continuity of the second fundamental form [G µν x ν ] Σ = 0; x ν is a unit vector. If we make use of the field equations (3), the continuity reads as [ T µν x ν ] Σ = 0. Using the full stressenergy tensor (8) and projecting in the radial direction x r = r, is written as [(T (P F ) rr + α θ rr ) r] Σ = 0. This leads to where the effective pressure comes from Eq. (11). On the r.h.s. we are in vacuum, hence the pressure must nullify. The equation system has been closed with the constrain (32), therefore p(R) = 0 is equivalent to request p(R) = 0 in (28). This equivalence makes the constant C not to vary from the perfect fluid solution once the anisotropies are considered. The value is With the constant fixed, we have fully determined the effective radial pressure of the anisotropic Durgapal-Fuloria solution A natural bound is obtained, α < 1. In Figure 2, it is shown the dependence of the pressure with a dimensionless radial coordinate r/R for different values of α. At first sight one can observe that the higher α is, the smaller the radial pressure becomes. The decreasing of the radial pressure is needed to produce the pressure anisotropy that is reflected in a change over the tangential pressure along the surface. The expression for the later pressure is written as p t (r; α) =p r + α 6 C 2 r 2 1 + 3 Cr 2 π 1 + Cr 2 1 + 9 Cr 2 2 .
The pressure (in both directions) must be a decreasing function along the radial coordinate. This condition restricts even more the values of α; higher values immediately triggers instabilities. In light of what was written in (38), the tangential pressure (41) determines another physical constrain for α: this pressure is meaningful as long as it remains positive everywhere p t (r) > 0; hence, so must be α > 0 to not contradict this statement in the surface where p r (R) = 0. From the latter equation, the anisotropy is directly computed; comparing with Eq. (14), we obtain Π(r; α) ≡ α 6 C 2 r 2 1 + 3 Cr 2 π 1 + Cr 2 1 + 9 Cr 2 2 .
Some comments are pertinent. The Durgapal-Fuloria solution is a fluid sphere with a solid crust.
In Figure 2 the density shows a discontinuity in the surface. The anisotropy smoothes this jump; the bigger the parameter α is, the lower the value of the density on the surface of the star. This behaviour immediately triggers the question on the profile of the effective mass function. This parameter has been defined in (35) and together with (33), is written as m(r; α) = 1 + α 2 (2 − 7 Cr 2 + C 2 r 4 ) (3 + Cr 2 )(1 + 9 Cr 2 ) m(r) .
An observer outside the star, would see a resulting mass M Schw surrounded by vacuum as it has been requested in (37). The continuity of the radial component of the metric (when crossing the star surface Σ) is direct: (36) identifies the Schwarzschild mass seen from outside with the effective mass of our solution ; i.e. M Schw ≡ m(R). Even more, a closer look at the mass function shows that the correction (αrf * )/2 in (35) vanishes at the surface r = R. This means that the effective total mass of the star is the same as the isotropic total standard mass; m(R) ≡ m(R) from (29). This issue is not surprising at all. The anisotropy mimics the radial pressure, hence the radial and tangential pressure start to drift apart in the region close to the solid surface. For this anisotropic behaviour to happen, both pressures must decrease in magnitude at the inner region. Of course this pressure discrepancy with respect to the isotropic solution makes the density to be disturbed. The equilibrium between gravitational collapse and pressure repulsion is modified; hence the mass function is redistributed to the center of the star. Despite this, the total mass of the anisotropic object remains unmodified. To conclude this section, let us determine the value of the remaining constant A from (37). The temporal component of the MGD metric (30) should match smoothly with the outer Schwarzschild region The constant A remains unchanged. This constant close one branch (α-dependent) of anisotropic solutions analogous to Durgapal-Fuloria; namely {ν ; λ ; ρ ; p r ; p t }. Of course, this solution is not unique. Different anisotropic solutions can be obtained starting from the Durgapal-Fuloria solution by means of requiring different constrains when closing the indeterminate system of equations. In next section we will consider a different constrain and we will see that a different anisotropic solution is obtained.

Density-like constrain for the anisotropy
Another useful constrain that gives an acceptable physical solution, is to impose that the anisotropy 'mimics' a density. This requirement is written as and closes the system of equations (22)- (24). The consequence of this ansatz is direct, the temporal Einstein equation for the perfect fluid (17) is identical to the temporal equation of motion for the θ µν -tensor (22). Equaling both equations, one notes immediately that the resulting equation has a total derivative structure. The integration is straightforward where K = 0 must be imposted for the invariants R, R µν R µν and R µνγσ R µνγσ to remain smooth and finite all over the inner region. Note that with this constraint the radial deformation is again totally determined by the solution of the perfect fluid. Eventually, one computes the relevant component of the metric; the minimally deformed component is written as in (16) (naming β to the coupling between sectors) as We can write the metric with the structure used in (36). The effective mass is written as a minimal deviation from the GR mass m(r) presented in (29) In the latter equality, if we make use of spherical coordinates and the corresponding relations, we have ρ r 2 dr = m/Ω 4 ; where the 4-dimensional solid angle is Ω 4 = dΩ and 2 Ω 4 /κ = 1. Of course, this is not surprising at all, the constrain for the anisotropy is to mimic the density, therefore the effective mass mimics the mass being proportional one to the other (unlike the previous case where the mass is exactly the same with respect to the standard GR solution).
Once the system is closed and the minimal deformation obtained, the remaining magnitudes are easily derived. As before this will be done by means of the smooth matching between the inner and outer region of the star.

Matching conditions
Here we will reproduce the same steps that we have done before in order to find the constants A and C; this time for the density ansatz (46). It is already known that the constant C is determined by the second fundamental form (38). Its value is In this case, the anisotropic sector has an influence on the integration constant. It is explicitly seen that in the limit of no coupling β → 0 the constant from Durgapal-Fuloria is recovered. Because of the ansatz where it has been required for the anisotropy to mimic the density, the effective value of the density is modified This is in complete accordance with the changes experienced by the mass. In what follows we will show that this solution only admits a minimal geometric deformation over the metric in only one 'direction'; the 'direction' to where the density and the mass is increase. The anisotropies restricted to the present constrain change the integration constants; for instance C is β dependent. An analysis over Eq. (50) shows that C increases when β becomes more negative. This behaviour makes Eqs. (49) and (51) to increase when β increase in modulus. Of course, theses both parameters can not increase without a limit; as in the previous case, anisotropies develop instabilities. In the first curve of Figure 3 it is seen how the density function increases in the inner region, while it slightly decreases its value over the surface's surroundings softening the crust. The mass function rises throughout the interior and the total effective mass is also increased. The remaining thermodynamic parameters are the effective radial pressure p r (r; β) = p − β C(3 + Cr 2 )(1 + 9 Cr 2 ) 7π(1 + Cr 2 ) 3 (52) and the effective tangential pressure p t (r; β) = p r − β C 2 r 2 π(1 + Cr 2 ) 2 . (53) As can be seen in the second graph of Figure 3, being the latter pressure different from the former, the anisotropy is developed. In the third and fourth set of curves, it is seen that both pressures are enhanced for the anisotropy to take place. When the mass increases, the enhancement of the density requires higher pressures for stability reasons. In order to get some insight in the underlying sign for β in the new solution, we will focus ourselves in the anisotropy, given by It is worth to note that this magnitude can not be negative (let us remind that C > 0) because if this were so, over the surface of the star where p r (R) = 0, we would have a negative tangential pressure p t (R) < 0 which is not physically acceptable. Therefore, positive pressures implies negative values for β. Now we have a physical domain for β. The constant A is found in an analogous manner than in the previous section, i.e. by means of Eq. (45). The value of this constant changes with β.
The usual constant of Durgapal-Fuloria is recovered in the limit of β → 0 as it should be.

Detectability and observational differences in anisotropic distributions
One of the many remarkable predictions of the theory of general relativity is the time dilation within a gravitational well. This results in footprints in the lines of the spectrum shifting towards the red. Although it is a useful quantity, particularly for astronomers, which allows to get some insight into compact stars physics, this effect is extremely difficult to deal with because of its complexity to be disentangle from the displacements and alterations due to various causes such as the Doppler, Zeeman and pressure effects among others. Theoretical derivations states that the redshift factor associated to a star comes when relating the proper time τ of the object with the observer clock t. This relation is given by the standard formula dτ 2 = g tt dt 2 that yields the following for the surface Therefore the relation between the emitted and observed frequency makes the redshift manifest [39].
Over the years, the study of anisotropies in compact objects has received considerable attention and this parameter is a simple way to contrast theory with observations. Formula (55) relates the measured redshift with the compactness parameter ξ = 2m/r of the star (that depends on the anisotropic coupling in this case α/β) introduced after Eq. (30). An observer outside would see the Schwarzschild metric, hence the redshift, that depends on ξ, is directly related with the total mass of the star (we are generating anisotropic contributions over fixed radius stars). Now, we want to investigate how this parameter evolve in our particular solutions.
Let us start with the fist solution derived in Section 4.1. The interest should focus in the mass function (44) evaluated over the surface. As we have explained after this equation the Schwarzschild mass remains unmodified with respect to the Durgapal-Fuloria mass, M Schw ≡ m(R) = m (DF ) (R). Hence, there is no observational evidence to differentiate an isotropic star to these anisotropic counterpart.
On the contrary, things change in the second case. The solution obtained in Section 4.1 is based on an increment of the total mass. Then a shift occurs when observing this anisotropic configuration The Schwarzschild mass coincides to the mass function over the surface; i.e. M Schw = m(R) in Eq.
(49). As it has been explained, the mass increases when β increases in modulus. Therefore, the compactness parameter is also increased; the star becomes more and more denser with β. Then, the parenthesis in (56) decrease and z(β) grows when |β| grows. This means that this anisotropic contributions increases the gravitational redshift as it is expected when the stars are more dense.

Anisotropizing an anisotropic Durgapal-Fuloria star
In section 3, we present a method to generate different anisotropic solutions of Einstein field equations using any well known perfect fluid as a seed. After this, we apply this prescription to the Durgapal-Fuloria perfect sphere. In section 4, with some reasonable constrains we found two novel physical anisotropic solutions analogous to the Durgapal-Fuloria compact star. The decomposition of Einstein equations (4)-(6) stem on the minimal geometric deformation (16); the anisotropic sector (22)-(24) is decoupled with respect to any perfect fluid sector (17)- (19). However, there is no need for the known sector to be a perfect fluid solution exclusively. Whatever solution of Einstein field equations, either a perfect or an anisotropic fluid, work as a seed for implementing the MGD decomposition. For instance, we can take any of the two previous founded solutions; e.g. the one obtained in section 4.1 given by {ν ; µ ; p r ; p t }. So as not to obscure how the method works, we will minimally deform the anisotropic solution along the radial component of the metric. While the temporal geometric function in Eq. (30) remains unchanged, the minimal distortion takes place only over the radial component of an anisotropic metric solution of Einstein equation; in this case (36). This deformation is caused by new generic sources of anisotropies (called ψ µν in order to avoid confusion with the deformed seed by θ µν ) that acts over the anisotropic energy momentum tensor (8) The Einstein field equations connecting the latter effective stress-energy tensor to the spacetime where the anisotropic radial component of the metric (34) has been used. The effective radial pressurep r = p r − βψ r r is computed from (60) The first integration constant C is obtain by means of the continuity of the second fundamental form, analogous condition to (38). Imposing the annulment of the latter effective radial pressure at the surface Σ, we get Both the integration constant C as well as the effective pressurep r recover the corresponding values: (39) and (40) in the limit of β → 0 or (50) and (52) when α → 0. Besides, this constant is required to plot the thermodynamic parameters.
In Figure 4 we present the corresponding evolution of the parameters of the theory: for a comparison, we include also the Durgapal-Fuloria isotropic solution (α = β = 0 in solid line). If for instance, one of the couplings move away from zero but the other remains null, the thermodynamic quantities behave as in Figures 2 (if α = 0) or 3 (if β = 0), as it is expected. After this, we plot the anisotropic seed to be minimally deformed by fixing the coupling α (red dashed-line). Finally, β drifts away the parameters again. We choose one smaller order or magnitude for the second deformation to make notorious the effect over the seed solution. An important statement is that as the seed is anisotropic and the corresponding tangential pressure is nonnull over the surface, then there is no restriction for β to be negative. β is allowed to be positive until either it decrease the tangential pressure until it becomes null, or the anisotropy becomes unstable.
Lest we forget, we include the expression of the two remaining parameters: first the effective densityρ = ρ + βψ t t , whose structure is analogous to (51), ρ(r; α, β) = (1 + β) ρ with ρ the seed density (43). And secondly, the anisotropic tangential pressurep t = p t − βψ ϕ ϕ that can be written asp t (r; α, β) =p r (r; α, β) + Π(r; α, β) . The anisotropy is now written as Π(r; α, β) = −β C 2 r 2 π(1 + Cr 2 ) 2 + (1 + β)α It is important to remark that in all expressions we recover both previous limits when either α or β are set to zero, and the standard Durgapal-Fuloria solution for α = β = 0. Let us conclude presenting the profile of the anisotropy over the surface Σ of the sphere. In Figure 5, we plot the function Π from Eq. (75) versus the couplings α and β. Nearby α closer to zero and β positive, the anisotropy becomes unphysical (Π < 0); thus deformations with this parameters are prohibited and excluded from the physical surface. In particular this procedure extends the range of physical values for β. Each time α grows, new possible values of β > 0 are released.
Besides should be noted that the anisotropy can be interpreted as a lineal combination of the first two studied cases Π(r; α, β) = Π(r; β) + (1 + β)Π(r; α) ; Each single parameter anisotropy have been computed in (42) and (54). In this case, the resulting Π(α, β) no longer coincides with β(ψ r r − ψ t t ) as mentioned in (14), because the seed is no more an isotropic fluid. Immediately, the linearity is translated to the stress-energy tensor; the components of the new ψ-sector can be written as a combination of the single minimal geometric deformations computed in Section 4 making the 'additive' character of the method manifest. If one starts with any perfect solution of GR, T (P F ) µν , a minimal deformation induced by an anisotropy subjected to a pressure structure, makes the stress energy tensor to become T µν = T (P F ) µν + α θ (pressure) µν . After this, a second minimal deformation acts over the already anisotropic solution, but now subjected to a density constrain. The new contribution is given by (77), therefore the effective energy-momentum tensor (58) is decomposed as This expression states the noncommutative structure of the MGD-decoupling method; the order in which the deformations take place matters. For instance, if we take as a seed the solution found in Section 4.2 where the deformation obeys a density-like constrain, and then we deform the anisotropic solution with a different constrain analogous to (32), the noncommutative character of respectively. Likewise, the effective energy-momentum tensor (78) becomes when the order of the two anistropizations is reversed. The reason of the noncommutativity is that the coefficients involved in the linear combinations of the components of the anisotropic tensor θ µν depends explicitly on the coupling constants. This is not surprising at all; we are dealing with nonlinear differential equations. The commutativity is likely to be lost in this kind of systems. From a perturbation theory point of view, one can think that the deformations over the metric (zero order) is due to the existence of the anisotropic term which acts at O(α); being α the coupling strength to the anisotropies. We must emphasize that, although the MGD approach seems like a perturbation technique, the method, in fact, is not, and this is easily visualized by noticing that the couplings do not necessarily have to be small, which is a crucial ingredient in perturbation theories. The deformation being a perturbation is just a well behaved limit of the theory, and means that we can softly deform the seed configuration. Being the theory noncommutative, successive and mixed perturbative deformations give different configurations depending on the order in which each of them are implemented. This provides infinite manners of deforming realistic configurations controlling rigorously the physical acceptability of the resulting anisotropic distribution.

Conclusions
In this paper we have presented different branches of solutions that models non-rotating and uncharged anisotropic superdense stars. Each anisotropic branch opens a possibility for new physically acceptable configuration obtained by guided deformations over the isotropic Durgapal-Fuloria stars, and exemplify some possible anisotropic distributions among the many that the MGD method generates. This prescription has been design to decouple the field equations of static and spherically symmetric self-gravitating systems. It associate the anisotropic sector with a deformation over the geometric potentials. In this work we have reported radial deformations only, but extensions to the temporal deformation may bring intriguing results. After the decoupling, one obtains a sector which solution is already known (seed sector) and the anisotropic sector which obeys a set of simpler 'pseudo-Einstein' equations associated to the metric deformation. It is worth to note that the minimal geometric deformations stem in an exclusively gravitational interaction between sectors; i.e. there is no exchange of energy-momentum among them.
When the equations are decoupled, no new information is introduced. Then we have an underdetermined system of equations; consistent constrains are needed. We have shown how intuitive constrain leads to new physical anisotropic solutions. Variations in the couplings between the seed and the anisotropic sector reveals consistent evolution of the thermodynamical parameter giving to the MGD method a new prove of validity. We also have discussed the observational features of the anisotropic sectors. When the anisotropy changes the compactness of the star, the observed redshift increases as it is expected. Not all anisotropic contributions have observational effects because some anisotropies only redistribute the thermodynamical parameters in the interior. However, when the anisotropy tweaks the compactness parameter, the star suffer a redshift. Therefore, observational data would bound the parameters of the model.
After presenting two branches of solutions that provides an infinite number of physical compact stars, we have proceed to generalize the method to deform anisotropic solutions. Any solution of Einstein equation admits a minimal deformation. Different anisotropic sources have additive effect, however these effects are noncommutative. The path to the final configuration matters, and normally deformations in reversed order produces different resulting configurations. Hence, the method provides a 'fine tunning structure' that generates an enormous amount of different physically acceptable anisotropic stars.