Integration of the Lane-Emden equation for relativistic anisotropic polytropes through Gravitational Decoupling: a novel approach

In this work we propose a novel approach to integrate the Lane-Emden equations for relativistic anisotropic polytropes. We take advantage of the fact that Gravitational Decoupling allows to decrease the number of degrees of freedom once a well known solution of the Einstein field equations is provided as a seed so after demanding the polytropic equation for the radial pressure the system is automatically closed. The approach not only allows to extend both isotropic or anisotropic known solutions but simplifies the computation of the Tolman mass whenever the Minimal Geometric Deformation is considered given that the $g_{tt}$ component of the metric remains unchanged. We illustrate how the the method works by analyzing the solutions obtained from Tolman IV, Durgapal IV and Wymann IIa isotropic systems as a seed for the integration.

The theory of polytropes is based on the polytropic equation of state, which can be written, in the case of a Newtonian-isotropic fluid like, where P denotes the isotropic pressure, ρ 0 stands for the mass (baryonic) density and K, γ, and n are usually called the polytropic constant, polytropic exponent, and polytropic index, respectively. Once the equation of state (1) is assumed, the whole system is described by the Lane-Emden equation that may be numerically solved for any set of the parameters of the theory. When the constant K is calculated from natural constants, the polytropic equation of state may be used to model a completely degenerate Fermi gas in the non-relativistic (n = 5/3) and relativistic limit (n = 4/3). In this case, Eq.
(1) provides a way of modeling compact objects such as white dwarfs and allows to obtain in a rather direct way the Chandrasekhar mass limit. Otherwise, if K is a free parameter, the models can be used to describe either an isothermal ideal gas or a completely convective star. These models related to isothermal ideal gas are relevant in the so-called Schönberg-Chandrasekhar limit [2].
Although local pressure isotropy is a very common assumption in the study of stellar objects, it is well known that many physical processes can produce deviations of the isotropy and/or fluctuations of the local anisotropy in pressure. These facts may be caused by a large variety of physical phenomena of the kind we expect to find especially in compact objects, so we find strong evidence that suggests that for certain ranges of density, a large number of physical phenomena can cause local anisotropy and therefore we must take it into account to describe realistic models. Among all the possibilities we would like to mention a few which might be particularly related to our primary interest. A possible source of anisotropy is related to intense magnetic fields observed in compact objects such as white dwarfs, neutron stars or magnetized strange quark stars [19][20][21]. Another source constitutes the high indexes of viscosity expected to be present in neutron stars, in highly dense matter produced by opacity of matter to neutrinos in the collapse of compact objects [22,23], and the superposition of two isotropic fluids. It is important to note that, although the degree of anisotropy may be small, the effects produced on compact stellar objects may be appreciable [24][25][26][27]. So, the assumption of an isotropic pressure is a very stringent condition, specially in a situation in which the compact object is modeled as a structure with high density (as neutrons stars, for example). Besides, the isotropic pressure condition becomes unstable by the presence of physical factors such as dissipation, energy density inhomogeneity and shear as it has been recently proven in [28]. These facts explain the renewed interest in the study of fluids not satisfying the isotropic pressure condition and justify the interest of extending the theory of polytropes to anisotropic compact fluids.
Thus, after assuming that the fluid pressure is anisotropic, the two principal stresses (say P r and P ⊥ ) are unequal and the polytrope equation of state reads As it happens in the Newtonian case, the introduction of the tangential pressure P ⊥ leads also in the relativistic case to an underdetermination of the problem requiring to impose an additional condition. Thus, in order to decrease the number of degrees of freedom the introduction of an additional condition is mandatory. In this regard, we can impose certain conditions on the metric variables as the vanishing of the Weyl tensor [29], implemented in [6] to obtain the conformally flat polytrope for anisotropic matter. This condition has its own interest, since it has been seen that highly compact configurations may be obtained with the specific distribution of anisotropy created by such a condition [29]. Other approaches as the Randall-Sundrum model [30] or 5dimensional warped geometries have served as an inspiration for other type of conditions, relating radial derivatives of the metric functions in spherically symmetric spacetimes, that produce self-gravitating spheres embedded in a 5-dimensional flat space-time (embedding class one). Even more, models embedded in five dimensional spacetimes satisfy the so-called Karmarkar or class I condition [31] (for recent developments see, [15,[32][33][34][35][36][37][38][39][40][41][42][43][44][45], for example) which allows to choose one of the metric functions as the one which generates the total solution. We can also consider an extra equation of the state as polytropic equation of state for the tangential pressure as reported in [17,46]. Alternatively, the use of a new concept of complexity based on the scalar Y T F appearing in the orthogonal splitting of the Riemann tensor has increased in recent years [47]. It is important to mention that, very recent studies have considered stellar anisotropic models with a certain complexity (or fulfilling the vanishing com-plexity condition Y T F = 0) and it has been established a relationship between families of solutions that have different complexities with the possible occurrence of cracking [48,49]. In this work we take an alternative route to integrate the Lane-Emden equation which consists of considering stellar interiors supported by anisotropic fluids fulfilling the polytropic equation of state for the total obtained radial pressure in combination with the Gravitational Decoupling (GD) [50] through the Minimal Geometric Deformation approach (MGD) (see  for recent developments). It is worth mentioning that the GD has been used in order to extend know solutions of Einstein field equations by coupling different sources, to decoupling a complex energy momentum tensor in simpler components, to find new solutions in theories beyond general relativity, and to find both static and stationary black hole solutions. In this work, we have demanded that the interaction between both the perfect ("seed") and the decoupling fluids (the θ extra sector) is such that the total effective source fulfills a polytropic equation of state. As we shall explain later, the use of GD by MGD allows to extend well-known solutions (the so-called "seed" sector) by deforming the seed metric and adding an unknown source to the matter sector in a suitable manner.
The advantage of such a procedure is that instead of providing two conditions to close the system of three differential equations in the static and spherically symmetric case, only one extra condition is required (because from the start it must be provided the seed solution). This procedures is clearly convenient in the integration of the relativistic Lane-Emden equation in the sense that it is sufficient with supplying the polytropic equation of state for the radial pressure to close the Einstein's system of field equations. More precisely, it is not necessary to propose any ansatz for the anisotropy or give some geometric condition but only a well-known stellar interior configuration. This work is organized as follows. In the next section we study the basic equations of general relativity as well as a summary of the theory of relativistic polytropes. In Section III, we review the main aspects of GD through the Minimal Geometric Deformation (MGD) formalism. We dedicate Section IV to obtaining the Lane-Emden equations for polytropes by gravitational decoupling and study some specific models. Finally, the last section is devoted to final remarks and conclusions.

A. The field equations and conventions
Let us consider a static and spherically symmetric distribution of anisotropic matter which metric, in Schwarzschild-like coordinates, is parametrized as where ν and λ are functions of r. The metric (3) has to satisfy Einstein field equations 1 The matter content of the system (describing an anisotropic fluid) is represented by the energymomentum tensor where ρ is the total energy density, is the four velocity of the fluid which satisfies u µ u µ = 1, and s µ is defined as with the properties s µ u µ = 0, s µ s µ = −1.
Furthermore, we shall consider that outside the fluid distribution the space-time is given by the Schwarzschild solution, namely where M represents the total energy of the system. In order to match the two metrics smoothly on the boundary surface r = r Σ = constant, we require continuity of the first and second fundamental forms across that surface. As a result of this matching we obtain the well known result, where the subscript Σ indicates that the quantity is evaluated at the boundary surface. 1 We are assuming natural units G = c = 1 From the radial component of the conservation law, one obtains the generalized Tolman-Oppenheimer-Volkoff equation for anisotropic matter which reads, Alternatively, using where the mass function m is as usually defined by we may rewrite Eq. (16) in the form where measures the anisotropy of the system. For the physical variables appearing in (19) the following boundary conditions apply which corresponds to a stellar configuration surrounded by the Schwarzschild vacuum [85] As already mentioned in the introduction, in order to integrate equation (19), we shall need additional information. In this work we propose a novel scenario; namely, we satisfy Einstein's system of field equations using the polytropic equation of state for the radial pressure together with the Gravitational Decoupling [50], by means of the Minimal Geometric Deformation approach (once the seed solution is given). In the next section we will implement the polytropic equation of state to construct the generalized anisotropic Lane-Emden equation which will arrive after introducing the polytropic equation of state in the TOV equation given by (19) followed by some redefinitions of the parameters involved.

B. Relativistic anisotropic polytropes
When considering the polytropic equation of state within the context of general relativity two different possibilities arise leading to the same equation in the Newtonian limit [5]. The first one, preserves the original polytropic equation of state (1) and the second case allows another (natural) possibility that consists in assuming that the relativistic polytrope is defined by, In this case the baryonic density ρ 0 is replaced by the total energy density ρ in the polytropic equation of state. The general treatment is very similar for both cases and therefore, for simplicity, we shall restrict here to the case described by (22). It can be shown that the relationship between the two densities is given by [5], As it is well known from the general theory of polytropes, there is a bifurcation at the value γ = 1. Thus, the cases γ = 1 and γ = 1 have to be considered separately. In the context of our work we want to focus in bounded (fluid) compact star models so the γ = 1 will not be considered.
Let us define the variable ψ by where ρ c denotes the energy density at the center (from now on the subscript c indicates that the variable is evaluated at the center). Now, we may rewrite (22) as with P rc = Kρ γ c . Replacing (22) and (24) in (19), the TOV equation can be written as where α = P rc /ρ c . Let us now introduce the following dimensionless variables from where (26) reads (see [5,7] for details) with Equation (31) in combination with (32) corresponds to the generalized Lane-Emden equation for an anisotropic fluid characterized by a polytropic equation of state.

III. GRAVITATIONAL DECOUPLING
In this section we introduce the GD by MGD (for more details, see [50]). Let us start by considering the Einstein field equations (4) sourced by certain T (tot) µν which can be written as where T (s) µν represents the matter content of a known solution, namely the seed sector, and θ µν describes an extra source coupled through the parameter β. It is essential to point out that the additional term βΘ µν is not considered a perturbation, i.e., the coupling parameter β could indeed be larger than unity (such coupling is introduced in order to control the effect of the unknown anisotropic source). Note that, since the Einstein tensor fulfills the Bianchi's identities, the total energy-momentum tensor satisfies the conservation equation so that, whenever ∇ µ T µν(s) = 0, the condition is automatic and as a consequence, there is no exchange of energy-momentum between the seed solution and the extra source θ µν (the interaction is entirely gravitational).
We are restricting ourselves to the spherically symmetric static anisotropic fluid case with internal metric (in Schwarzschild coordinates) given in (3). In this case the source can be expressed like, Note that, given the symmetry of the system, the components of the extra source, θ µν , depend on the radial coordinate only. To be more precise, although we are assuming anisotropic systems, it means that both the radial and the tangential pressures are different so that θ 11 = θ 22 = θ 33 . Even more, although an explicit dependence on the angular coordinates could represent an anisotropic systems too, this break our assumption that the model we are considering here is a spherically symmetric one. Then, equations (4), (33), (36) and (37) lead to the fact that the total energy-momentum tensor (T (tot) µν ) satisfies the system of field equations (8), (9) and (10), where now, in the left hand side we have Because, in general, θ 1 1 = θ 2 2 , we find that the system represents an anisotropic fluid. Is clear that the non-linearity of Einstein's equations avoids that the decomposition (33) leads to two sets of equations; one for each source involved. Nevertheless, the decoupling is possible in the context of MGD as we shall demonstrate in what follows.
Let us introduce a geometric deformation in the metric functions given by where f is the so-called decoupling function and β is the same free parameter that "controls" the influence of θ µν on T (s) µν . It is worth mentioning that although a general treatment considering deformation in both components of the metric is possible, in this work we shall concentrate in the particular case where the deformation is only implemented on the g rr component. Now, replacing (41) in the system (8)-(10), we are able to split the complete set of differential equations into two subsets: one describing a seed sector sourced by the conserved energy-momentum tensor T and the other set corresponding to quasi-Einstein field equations sourced by θ µν , where we have defined ρ θ = βθ 0 0 , P θ r = −βθ 1 1 and P θ ⊥ = −βθ 2 2 . Note that, as the seed sector is sourced by a conserved energy-momentum tensor, namely ∇ µ T µν(s) = 0, the the components of θ µν satisfy the conservation equation ∇ µ θ µ ν = 0, given by Remarkably, although the quasi-Einstein equations differ from Einstein equations, the expression given in (48) is completely analogous to the anisotropic Tolman-Opphenheimer-Volkoff (TOV) equation. Next, to complete the process, we require to satisfy the matching conditions (12), (13) and (14) on the boundary surface Σ.
To conclude this section, we would like to emphasize the importance of GD by MGD as a useful tool to find solutions of EFE. As it is well known, in static and spherically symmetric spacetimes sourced by anisotropic fluids, EFE reduce to three equations and five unknowns, namely {ν, λ, ρ, P r , P ⊥ }. In this sense, two auxiliary conditions must be provided: metric conditions, equations of state, complexity of the system, etc. However, in the context of MGD a seed solution is given, namely, a metric {ν, µ} that solve Eqs. (42)- (44). Now, note that as we are only deforming the radial metric by Eq. (41) , the decoupling sector system given by Eqs. (45)-(47) has a metric given by the pair {ν, f }. In this regard, it worth emphasising that through the MGD the number of degrees of freedom automatically reduces from five to four, namely {f, θ 0 0 , θ 1 1 , θ 2 2 }, given that both the seed and the decoupling sector share the same temporal metric so that only one extra condition is required. In general, this condition is implemented in the decoupling sector given by Eqs. (45), (46) and (47) by some equation of state (or condition) which leads to a differential equation for the decoupling function f . For example: i) The mimic constrain for the density, namely θ 0 0 = ρ (s) , which ensures that both secctors (seed and decoupling) have the same energy density profile. ii) The mimic constrain for the pressure, namely −θ 1 1 = P (s) r , which ensures that both sector have the same profile for the radial pressure. Another possibility is to supply a condition for the total solution (which of course entails a relation between the two sectors). For example, if we consider the simple barotropic condition P r = ρ, this encodes the relation so in the context of GD the system is closed enough. In this work we propose the polytropic equation of state which encodes the relation between both sectors. As we shall see later, we demonstrate that this strategy allows to write the Lane-Emden equation and the anisotropy present in the total solution in terms of the pair (η, ξ), defined in (27) and (30), and the metric variables of the seed sector.

IV. GRAVITATIONAL DECOUPLING AND ANISOTROPIC POLYTROPES
In this section we will outline the procedure to obtain the structural Lane-Emden equations by means of the GD. In order to do so, we note that integration of Eqs. (31) and (32) provides (ψ, η) as a function of ξ whenever the anisotropy ∆ = P ⊥ − P r is supplied. In this work we take advantages of the GD to write ∆ in terms of (η, ξ) and the seed metric functions (ν, µ) as we shall explain in what follows.
First, note that the total anisotropy can be written as with ∆ (s) = P (s) ⊥ − P (s) r and ∆ θ = P θ ⊥ − P θ r . Now, by replacing (43), (44), (46) and (47) in (52) we obtain where, from now on, primes indicate derivation respect to the variable ξ. Next, from (18) and (41) and using (27), (28) and (30), the decoupling function f and its derivative can be written as with η given by (30). Finally, replacing (54) and (55) in (53) we arrive at ∆ = πρ c α(n + 1)8π × µ e −µ + 2α(n + 1) (η − η ξ) ξ 2 ν + 2 ξ where we have used the dimensionless variables defined by equations (27)- (30). Note that, the final expression of the total anisotropy only depends on (η, ξ) (dimensionless mass-function and radius) and the known metric functions of the seed sector (ν, µ), as expected. At this point, it is worth emphasizing the system is closed in the following sense. As we stated before, for any anisotropic and spherically symmetric system, the problem reduces to solving three Einstein field equations for five unknowns, namely the metric functions {ν, λ} and the variable of the matter sector {ρ, P r , P ⊥ }. In this regard, two additional conditions are required in order to close the system. In this work, one of this conditions is supplied by the polytropic equation of state relating the radial pressure with the energy density given by Eq. 22. The other condition corresponds to that the temporal metric ν is a known function once a seed solution is specified.
To complement the discussion, we proceed to calculate the Tolman-Whittaker mass [86], m T , which in the context of GD is a simple task and this fact in itself reinforces the use of this tool novelty linked to the study of polytropes. Indeed, using the definition of the active gravitational mass, it is noted that the only relevant information that should be provided is the mass function that is encoded in e −λ . This is so, given that the g tt metric component is known since it corresponds to the metric function of the seed solution. Thus, replacing (18) in (57) and using (30), we arrive at from where with For the numerical calculations it is convenient to change to the following dimensionless variables in terms of which the Tolman mass can be written as We observe from (62) that each anisotropic polytropic model is characterized by a unique y.
In what follows we shall proceed to integrate the Lane-Emden equations by considering different seed sectors and be able to obtain, through the MGD formalism, relativistic anisotropic polytropes as the total solution. To achieve this, we will introduce the final anisotropy function given in (56) into the generalized TOV or Lane-Emden equation for an anisotropic fluid characterized by a polytropic equation, given by (31).

A. Tolman IV
In this section we shall use the well known Tolman IV solution [87] as a seed so we can close the system and find an anisotropic total solution that also complies with the fact of satisfying a polytrope equation of state, this, carried out within the framework of GD by means of MGD, as was previously explained. Then, we have where a, b and c are constants, that have the same dimensions of r. In order to describe our model, we have fixed the values of a = 1.41421, b = 2.23607 (in dimensions of a length) after a numerical probing in the parameters space.
In Fig. 1 we perform the integration of system of equations (31) and (32) for the values of the parameters indicated in the figure legend where we represent ψ (energy density) as a function of the dimensionless variable ξ for different values of the polytropic index n for the polytrope model obtained with the Tolman IV seed-solution. We observe that ψ is monotonously decreasing, as expected for well behaved general relativistic polytropes and the radial pressure P r vanishes at the surface as required by the continuity of the second fundamental form.
In Fig. 2 it is shown the mass-function η as a function of the variable ξ for different values of the index n. The results also depends on the "rigidity parameter" (α = P rc /ρ c ) related with the relativistic limit of the Lane-Emden equation (α → 0 corresponds to the Newtonian limit). Note that for each polytrope model (given by n), η is an increasing function that grows until the numerical evaluation stops at the surface of the compact object (ξ Σ ), where naturally coincides with the total Schwarzschild mass M . In Fig. 3 its shown the "surface potential" y which measures the degree of compactness as a function of the polytropic index n for different values of α. The parameter y is a relevant variable since it measures the compactness of the configuration, it is related with the redshift and it will also be useful in the analysis of the behavior of the active gravitational mass (Tolman-mass). Note that to produce the numerical behaviour for η/η Σ T (represented in Fig. 4) we must feed the data of y obtained from previous Fig. 3. We see that by increasing the index of the polytrope, the Tolman mass generally increases and is concentrated towards the outer surface of the object. In terms of stability, and keeping in mind the physical meaning of the Tolman mass, it may be said that more stable configurations correspond to smaller values of the Tolman mass concentrated towards the center of the sphere, since these are associated to a sharper reduction of the active gravitational mass in the inner regions, thereby providing a clear physical picture of the described scenario. Fig. 5 shows the dependence of the local pressure anisotropy with the dimensionless radial variable ξ and the polytrope index. Note that ∆ starts at zero, where P r = P ⊥ , is a monotonically increasing function as expected, and the radial pressure P r vanishes at the surface as required by the continuity of the second fundamental form.It is observed that the anisotropy function decreases with the increase of the polytrope index n, and this fact can be related to the Tolman mass distribution, existing a relationship between the stability of compact objects and their local pressure anisotropy, an issue that has been widely reported [5,7]

B. Durgapal IV
Now, we apply the procedure described in section 4 to the Durgapal IV [88] interior solution as a seed in the GD framework. The starting point is e ν = a(cr 2 + 1) 4 e −µ = 7 − 10cr 2 − c 2 r 4 7(cr 2 + 1) 2 + bcr 2 (cr 2 + 1) 2 (1 + 5cr 2 ) 2/5 , where a, b (with units of the inverse of a length squared) and c (dimensionless) are constants. In this case, the fixed constant values used in the Durgapal IV seed are a = 0.6254, b = 2.1479. and c = 1.372281. In figure 6 it is shown the integration of Eqs. (31) and (32), with previous use of the anisotropy function (56), for the Durgapal IV seed-solution, plotted for the values of the parameters indicated in the figure. Note that ψ is monotonously decreasing as expected. Again, the numerical analysis depends on, both, the polytropic index n and the "rigidity" parameter α (related to the ratio of pressure and density at the center of the compact object). Each n defines a specific stellar object that, for this model (Durgapal IV seed-solution) presents different boundaries, ξ Σ , consistent with the fact that the radial pressure P r vanishes at the surface (Σ) as required by the continuity of the second fundamental form. However as we have mentioned, the ensuing qualitative behaviour, namely, larger values of ψ (for bigger n) for smaller values of ∆ (everywhere throughout the sphere), is maintained for a wide range of values. It is straightforward to obtain this behavior by comparing with Fig. 10. Similar behaviors have been reported in [5].
The mass-function η as function of ξ for different values of n is presented in Fig. 7. The numerical evaluation has been stopped when η reaches ξ Σ , where the mass function, the Tolman mass, and the total mass become equal. Of course, each polytrope model, with different boundary (fact observed in Fig. 6) in turn also present a different total contained energy. The mass function η, related with the radial metric function λ, constitutes a parameter that allows us to analyze the energy contained within the sphere, although the Tolman mass has been shown to be more useful to describe the active gravitational mass [5,24].
The parameter y is plotted in Fig. 8 as function of n for different α and the normalized Tolmann mass η/η Σ T is shown in Fig. 9. Again, the active gravitational mass increases from the center to the surface, showing some peculiar behaviors for some values of the polytropic index (orange line: n = 3). This fact has been reported before (see references [7] for a more detailed discussion).
Finally, Fig. 10 is dedicated to exposing the behavior of the local anisotropy of the pressure as a function of the dimensionless parameter ξ. For this anisotropy function we obtain the usual behavior. We note that the anisotropy is an increasing function and when the index of the polytrope increases the anisotropy given in some layer of the stellar object decreases.

C. Wyman IIa
In this section we use the Wyman IIa solution [89] as a seed, namely with a, b and c are constants. For this model, the fixed constant values are a = 1.5297, b = 1.0. and c = 2.148. In Fig. 11 we show the behavior of the matter sector through integration of the Lane-Emden equation and we get ψ as function of ξ for the polytropic model represented by Wyman IIa seed-solution. Note that the density is positive inside the compact star, reaches its maximum at the center and decreases monotonously outwards, as expected. Unlike the two previous models, it is observed that all the configurations obtained for different values of the polytropic index, have the same boundary surface where P rΣ = 0.
In Fig. 12 we plot the η-mass as function of ξ for different values of n which grows appropriately from zero to its total limit value and larger values of the massfunction are obtained for smaller values of n.
The surface potential y and the normalized Tolman mass, for a selection of values of the parameters are plotted in figures 13 and 14 respectively. Also, the conclusions extracted from these figures are basically the same as the ones reached from the previous models. At any rate, maximal values of y correspond to bigger values of α, which represents a consistent and expected result for all models.
In Fig. 15 we show the behavior of the anisotropy function ∆ as a function of the redefined variable (ξ) for different values of the duplet of parameters (n, α) indicated in the caption of the corresponding figure. An increasing behavior is manifested and the ∆-function decreases with increasing n. In its turn, more compact configurations correspond to smaller values of the anisotropy and smaller values of the Tolman mass in the inner regions. In other words, for this case, smaller values of n (corresponding to more compact configurations) reach stability by reducing the active gravitational mass in the inner regions. Therefore, it may be inferred from this figure that more stable configurations correspond to smaller values of n since they are associated to a sharper reduction of the Tolman mass in the inner regions. In the same way, smaller polytrope indices imply an increase in the local pressure anisotropy, so that this function generates more stable situations for more compact spheres.

V. DISCUSSION
In order to obtain a general relativistic polytrope model we must integrate the hole system of EFE. If we are assuming that the fluid pressure is anisotropic, we need to provide further information about the anisotropy inherent to the problem under consideration. The introduction of a new variable yields an additional degree of freedom and therefore Eq. (2) is not enough to integrate the Lane-Emden equation. For doing that, there are different ways to deal with this problem. It can be assumed an ansatz for the anisotropy allowing a specific modeling [5]. Certain physical conditions can be imposed on the metric variables such as the vanishing of the Weyl tensor (conformally flat polytrope models) [6] and the Karmarkar Class I polytrope [15]. Also, the vanishing complexity polytrope model has been developed [47].
Since polytropes represent a variety of fluid systems with a wide range of applications in astrophysics (e.g. Fermi fluids), we have described hereby a general framework for modeling general relativistic polytropes (with local pressure anisotropy), by means of the Gravitational Decoupling process, through the Minimal Geometric Deformation approach. Thus, we have built another method to obtain the generalized Lane-Emden equations, that allows us to find a specific model. Assuming the polytropic equation of state for the radial pressure allows us to obtain an expression for the total anisotropy of the system in terms of the dimensionless defined mass-function (η), given in (30), and the metric variables of the seed sector. This represents an enormous advantage, since now we are able to obtain a great variety of general relativistic polytropes models choosing well-behaved known solutions as seeds. In this work we have studied three polytrope models considering the Tolman IV, Durgapal IV and Wyman IIa seed-solutions. For each case we have found and integrated numerically the full set of equations; the generalized Lane-Emden equations for anisotropic matter. Such extensions of the polytrope solutions to the general relativistic case are mandatory if one has to deal with ultra compact objects such as neutron stars, where general relativistic effects cannot be neglected. However it should be stressed that the above mentioned models are presented here with the sole purpose to illustrate the method. The natural way to obtain models consists in providing the specific information about the kind of anisotropy present in each problem, but having said this, it is important to emphasize that the obtained models exhibit some interesting features which deserve to be commented.
Thus, we observe in Figs. 1, 6 and 11 that bounded configurations exist for a range of values of the parameters involved and only in the Wyman IIa seed-solution model all the radii of the border coincide. However due to the existence of a larger number of parameters than in the isotropic case, the conditions for the existence of finite radius distributions are more involved than in the latter case. Figs. 3, 8 and 13 show the behavior of the "surface potential" and for all the models considered, y grows as n increases and has bigger values for bigger rigidity parameter (α). It should be emphasized that the relationship between the maximal values of y (maximal surface redshift) and the local anisotropy of pressure, has been discussed in detail in the past [90][91][92]. The explanation for such interest is easily understood, if we recall that the surface redshift is an observable variable, which thereby might provide information about the structure of the source, that is related with the polytrope model (the index n) that describes the fluid and also its degree of anisotropy. This is so, since each anisotropic polytropic model is characterized by a unique y as can be seen in Eq. (62). The polytropic models favoring higher redshifts are clearly exhibited in the Fig. 13, for the Wyman IIa solution, where a peculiar behavior is observed for the highest value of α. This "anomalous" behavior could be related to the extreme (maximal) value of this parameter (α = 1), that corresponds to the stiff equation of state P r = ρ, which is believed to describe ultradense matter [93].
The correspondence mentioned above, between the y-parameter and the polytropic index n, suggests that a relationship can be established between compactness and the local anisotropy distribution of the fluid, which in turn is linked to interesting facts that appear in the search for the stability of the system, except for particular cases where such a correspondence is broken. In order to delve deeper into this feature, in Figs. 4, 9 and 14 we have investigated the behaviour of the Tolman mass (that could be related with the stability) within the sphere for each model. These facts have been extensively reported in [5-7, 24, 94] where is discussed the efficiency to diminish the Tolman mass in the inner regions and to concentrate it in the outer ones depending on the anisotropic factor which brings out the role played by the anisotropy in the stability of the fluid configuration.
In principle, a somewhat speculative argument, states that it could be possible to investigate the characteristics of the fluid of the star by determining n by means of the surface potential. This would be valid for such configura-tions where the general relativistic effects as well as the inclusion of pressure anisotropy, are unavoidable, and are modeled resorting to a polytropic equation of state. So, a potential application of the approach presented here could apply to the study of super-Chandrasekhar white dwarfs that have masses of the order of 2.8M . Nevertheless, care must be exercised with the fact that some of the physical phenomena present in such configurations (e.g. intense magnetic fields) could break the spherical symmetry.