Class I polytropes for anisotropic matter

In this work we study class I interior solutions supported by anisotropic polytropes. The generalized Lane--Emden equation compatible with the embedding condition is obtained and solved for a different set of parameters in both the isothermal and non--isothermal regimes. For completeness, the Tolman mass is computed and analysed to some extend. As a complementary study we consider the impact of the Karmarkar condition on the mass and the Tolman mass functions respectively. Comparison with other results in literature are discussed.


I. INTRODUCTION
Polytropic equations of state have played a remarkable role in astrophysics (see [1,2] and references therein), and have been extensively used to study the stellar structure. For example, white dwarfs have been modeled considering Newtonian polytropes [3,4] while in more compact configurations (e.g. neutron stars, super Chandrasekhar white dwarfs) [5,6] general relativistic polytropes have been extensively studied (see, for example [7][8][9][10][11][12][13][14][15][16][17], and the references therein ). For an isotropic fluid, the theory of polytropes is based on the polytropic equation of state which reads 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, * Electronic address: anthony.ramos@yachaytech.edu.ec † Electronic address: cynthia.arias@yachaytech.edu.ec ‡ Electronic address: ernesto.fuenmayor@ciens.ucv.ve § Electronic address: econtreras@usfq.edu.ec 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 isotropy is a very common assumption in the study of compact objects, there is strong evidence that suggests that for certain ranges of density, a large number of physical phenomena can cause local anisotropy. A possible source is related to intense magnetic fields observed in compact objects such as white dwarfs, neutron stars or magnetized strange quark stars [18][19][20][21][22]. Another source is 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 [23][24][25], and in the superposition of two isotropic fluids, to name a few. It is important to note that, although the degree of anisotropy may be small, the effects produced on compact stellar objects may be appreciable [26][27][28][29][30][31][32]. 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, this has been recently proven in [33]. Thus, after assuming that the fluid pressure is anisotropic, the two principal stresses (say P r and P ⊥ ) are unequal and the polytrope arXiv:2103.05039v1 [gr-qc] 8 Mar 2021 equation of state reads Note that, the introduction of P ⊥ as a new variable yields an additional degree of freedom and therefore Eq. (2) is not enough to integrate the Lane-Emden equation. Thus, in order to decrease the number of degrees of freedom the introduction of an additional equation of state is mandatory. Alternatively, we can impose certain conditions on the metric variables as the vanishing of the Weyl tensor, implemented in [6] to obtain a 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.
Other approaches as the Randall-Sundrum model [34] or 5-dimensional 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). In this regard, it is clear that embedding of four-dimensional spacetimes into higher dimensions is a very useful tool in order to generate astrophysical interior models. Even more, models embedded in five dimensional spacetimes satisfy the so-called Karmarkar or class I condition [35] (for recent developments see, , for example) which allows to choose one of the metric functions as the one which generates the total solution. In this work we consider stellar interiors supported by anisotropic fluids fulfilling the polytropic equation of state for the radial pressure and additionally, we implement the class I condition to close the system and construct the corresponding generalized Lane-Emden equation. The manuscript is organized as follows. In the next section we review the main aspects of anisotropic polytropes and obtain the generalize Lane-Emden equation. Next, we introduce the class I condition and solve the Lane-Emden equation numerically. We conclude this section with a brief discussion about the incidence of the class I condition on the Tolman mass. 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 [63] The matter content of the system is described by the energy-momentum tensor where ρ is the total energy density, is the four velocity of the fluid and s µ is defined as with the properties s µ u µ = 0, s µ s µ = −1.
Replacing (3), (5), (6) and (7) in (4), we have where primes denote derivative with respect to r. Furthermore, we shall consider that the fluid distribution is surrounded by the Schwarzschild vacuum, 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. 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 as 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 As already mentioned in the introduction, in order to integrate equation (19), we shall need additional information. In this work we use the polytropic equation of state for the radial pressure and the class I condition in order to close the system. In the next section we will implement the polytropic equation of state to construct the general Lane-Emden equation.

B. Relativistic polytrope for anisotropic fluids
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. Let us first consider the case γ = 1. Now, defining 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), we may rewrite (22) as with P rc = Kρ γr c . Now, replacing (22) and (24) in (19) we obtain where α = P rc /ρ c . Let us now introduce the following dimensionless variables from where (26) reads which corresponds to the generalized Lane-Emden equation for this case.
Let us now consider the isothermal case: n = ±∞, γ = 1. In this case, the definition of the variable ψ becomes, and we can write equation (22) as Replacing (30) and (31) in TOV equation (19), we obtain Introducing dimensionless variables, Eq. (32) becomes which corresponds to the Lane-Emden equation for the isothermal case.

III. THE CLASS I CONDITION
Embedding of four-dimensional spacetimes into higher dimensions is an invaluable tool in generating both cosmological and astrophysical models [43]. As it is well known, the Karmarkar condition is necessary and sufficient to ensure class one solutions [35] which for spherically symmetric space-times reads (R θφθφ = 0) leading to with e λ = 1. Now, using (24), (25) and (37) in (20) we obtain from where, after using the set of dimensionless variables (27) and (28), we arrive at Finally, replacing (39) in (29), the class I generalized Lane-Emden equation for γ = 1 reads We can repeat the same procedure for the γ = 1 (isothermal) case. Again, after using (30), (31) and (37) in (20) we obtain, from where, by means of the definitions (33) and (34) we have Finally, introducing (42) in (35) we arrive to the class I generalized Lane-Emden equation for γ = 1 (isothermal case), To complement the discussion, we proceed to calculate the Tolman mass (for the γ = 1 case) [58], defined by or, alternatively [26] m T = e (ν+λ)/2 m + 4πr 3 P r .
For the numerical calculations it is convenient to change to the following dimensionless variables in terms of which equation (51) reads, where G(x) becomes, This procedure, where the integration of the TOV equation is carried out, allows us to find the metric function ν (using the appropriate boundary conditions) and finally Tolman's mass, (56) and (57), is totally equivalent to the one developed in [59] for the different cases exposed. Now, the above results will be integrated numerically, with the corresponding boundary conditions, in order to obtain the solution for each case considered. Figure 1 shows the behaviour of e ν as a function of x for different values of the parameters involved. Note that, the function is monotonously increasing as required. In figures 2 and 3 it is shown the integration of Eqs. (43) (isothermal case) and (40), respectively, for the values of the parameters indicated in the figure legend. Note that for the isothermal case (γ = 1) the configurations are unbounded and as a consequence, it is meaningless to define a surface potential or the total mass. However, for the γ = 1 case ψ is monotonously decreasing as expected, and the radial pressure P r vanishes at the surface as required by the continuity of the second fundamental form. The parameter y ("the surface potential") which measures the degree of compactness is plotted in Fig. 4 as function of the polytropic index n for different values of α. Note that the plots have discontinuities in some point n d , where we observed that for n < n d , the solution is well-behaved, and for n ≥ n d , the numerical solution starts to have either imaginary components or strange oscillatory behavior. Figures 5 and 6 display the Tolman mass (normalized by the total mass), for the case γ = 1, as function of x for the selection of values of the parameters indicated in the legend. The behavior of the curves is qualitatively the same for a wide range of values of the parameters. We observe that as we move from the less compact configuration (curve a) to the more compact one (curve d), the Tolman mass tends to concentrate on the outer regions of the sphere. Therefore, it may be inferred from the figures 5 and 6 that more stable configurations correspond to smaller values of y since they are associated to a sharper reduction of the Tolman mass in the inner regions, as usually happens if the anisotropy is fixed by the imposed condition. This is the usual strategy adopted by the fluid distribution to keep the equilibrium; it tries to concentrate the Tolman mass in the outer regions. This behavior was already observed for different families of anisotropic polytropes discussed in [5]. The same result is obtained from the Fig

IV. KARMARKAR CONDITION AND TOLMAN MASS
In this section we would like to establish some relationships that may be useful to analyze the anisotropy inherited by the Karmarkar condition. In order to relate the influence of the Karmarkar condition with the mass function m and the Tolman mass, m T , we will take an instructive path to complete our analysis of the fluid distribution of the compact body. From (18)  usual definition of the mass function, Note that, we can write condition (36), for our static and spherically symmetric case, in order to define a new variable, K, in the following form, where the Karmarkar condition is satisfied when K = 0. In this way we can keep track of how this condition influences our anisotropic distribution of fluid of the compact object. Using the components of the Riemann tensor appropriate to the symmetry we may obtain from (59) the following expression, which, after using Einstein equations (8), (9) and (10), can be written as Now, we proceed in two different ways in order to interpret the influence of the Karmarkar condition over our interior solution. First, from (8) and (9), we obtain the derivatives of the metric functions ν and λ in the form Now, replacing these expressions in the last term of (61) we obtain which relates the mass function m with the K condition, defined in (60). In contrast, if we use (44) and Einstein equation (8), we can write the last term of (61) in terms of the Tolman mass function, m T , as Finally, using (45) in the previous equation we have which constitutes a relation between the mass function, the Tolman mass and the K condition. So, we have been successful in relating the components of the energymomentum tensor, the mass function and the Tolman mass with the defined K condition for the symmetry considered in our development. This could allow us to discover the influence over these mass functions by imposing the Karmarkar condition and then be able to find out how this metric (embedding) condition influences the distribution of our compact fluid. In order to illustrate the effect of the Tolman mass on K we shall consider a simple example, namely, an anisotropic solution with vanishing radial pressure [62], with which ensure the fulfilling of the matching conditions (12), (13) and (14). Next, replacing (67), (68), (69), (70), (71), (72) and (73) in (66) and using (45) to express the mass function in terms of m T we arrive at In

V. CONCLUSIONS
In this work, we have developed a general method to construct locally anisotropic polytropes, and we have applied it to the specific case of a class I solution for sherically symetric interior space-times. As it occurs in the specific case of conformally flat polytrope spheres [6], our models are necessarily anisotropic (principal stresses unequal) and therefore are not continuously linked to isotropic polytropes (just like it ocurres in other reported anisotropic distributions [60]). Since the inclusion of pressure anisotropy implies an additional degree of freedom, the integration of the ensuing Lane-Emden equation, in the general case, requires additional information.
Usually this information is provided through equations of state or other type of assumptions about the nature of the anisotropy. Here, we have supplied such additional information by assuming that the polytrope has the specific anisotropy derived from de Karmarkar condition. This represents an heuristic assumption whose validity will be confirmed (or not) from its application to specific problems. Then, the geometrical constraint, chosen to work with, produce a specific type of anisotropy leading to a Lane-Emden equation that required a numerical analysis, given the high nonlinearity of the equations. We could differentiate two possible cases depending on whether γ = 1 or γ = 1. We have integrated the Lane-Emden equations for these cases, for a very large set of values of the parameters, even though, only a very specific set for each case is exhibited, this because the qualitative behavior of the system does not change much for a wide range of the parameters. The main reason to present such models was not to describe any specific astrophysical scenario, but to illustrate the applications of our approach. However, the protocol here developed could be useful for dealing with several phenomena in which anisotropic polytropes appear, with special emphasis in the physics of compact objects and so on, the obtained models exhibit some interesting features which deserve to be commented. We do not know if distributions that are usually represented as polytropes, and that also constitute very compact objects, like neutron stars (that are described by a Fermi gas), or Super-Chandrasekhar white dwarfs could be adequately modeled in this way. For such configurations it is evident that general relativistic effects as well as the inclusion of pressure anisotropy are unavoidable, and therefore its study could be carried on using the approach presented here. Here too, care must be exercised with the fact that some of the physical phenomena present in such configurations (e.g. very strong magnetic fields) could break the spherical symmetry.
It is interesting to highlight that polytropes with γ = 1 are used to model systems with non-degenerate isothermal cores which play an important role in the analysis of the Schönberg--Chandrasekhar limit. Any distribution of pressure anisotropy would affect the Schönberg-Chandrasekhar limit not only by affect-ing the structure of the polytrope, but also by the modifications of the virial theorem, introduced by the anisotropy. In the last section we have followed a similar strategy previously used to express the Tolman mass in terms of traces of the Riemann tensor (matter sector) and the only non-vanishing scalar associated to the electric part of the Weyl tensor in static and spherically symmetric spacetimes, E (see [61], for example). In this sense, it is clear how the different parts affect the Tolman mass. For example, for conformally flat solutions the Weyl condition reads E = 0 which leads to a geometric constraint which allows to close the system. Similarly, here we have obtained expressions relating the scalar K with the traces of the Riemann tensor, the mass function m (64) and the Tolman mass (66) and class one solutions are obtained when the Karmarkar condition is fulfilled, namely K = 0. In this sense, expressions (64) and (66) may be useful to investigate the role played by Karmarkar embeding condition in our fluid distribution corresponding to the compact object. These aspects must be approached with the necessary care and detail in future works. We do not know if there are any real physical phenomena that can produce a distribution of local anisotropy similar to the one inherited through the metric Karmarkar condition described in our models. So, it could be interesting to be able to provide a physical meaningful picture for a source of the anisotropy that we have obtained (by means of magnetic field or charge distribution) of the type one could expect to find in some astrophysical settings. Although it is certainly an interesting question, the answer lies beyond the scope of this work. However, we believe that our approach provides a useful tool to carry out this and other investigations.
Eur. Phys. J. C manuscript No. (will be inserted by the editor) Class I polytropes for anisotropic matter.
A. Ramos a,1 , C. Abstract In this work we study class I interior solutions supported by anisotropic polytropes. The generalized Lane-Emden equation compatible with the embedding condition is obtained and solved for a different set of parameters in both the isothermal and nonisothermal regimes. For completeness, the Tolman mass is computed and analysed to some extend. As a complementary study we consider the impact of the Karmarkar condition on the mass and the Tolman mass functions respectively. Comparison with other results in literature are discussed.
For an isotropic fluid, the theory of polytropes is based on the polytropic equation of state which reads 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 equaa e-mail: anthony.ramos@yachaytech.edu.ec b e-mail: cynthia.arias@yachaytech.edu.ec c e-mail: ernesto.fuenmayor@ciens.ucv.ve d e-mail: econtreras@usfq.edu.ec tion 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 isotropy is a very common assumption in the study of compact objects, there is strong evidence that suggests that for certain ranges of density, a large number of physical phenomena can cause local anisotropy. A possible source is related to intense magnetic fields observed in compact objects such as white dwarfs, neutron stars or magnetized strange quark stars [16][17][18][19][20]. Another source is 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 [21][22][23], and in the superposition of two isotropic fluids, to name a few. 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][28][29][30]. 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 in-arXiv:2103.05039v1 [gr-qc] 8 Mar 2021 2 homogeneity and shear, this has been recently proven in [31]. 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 Note that, the introduction of P ⊥ as a new variable yields an additional degree of freedom and therefore Eq.
(2) is not enough to integrate the Lane-Emden equation. Thus, in order to decrease the number of degrees of freedom the introduction of an additional equation of state is mandatory. Alternatively, we can impose certain conditions on the metric variables as the vanishing of the Weyl tensor, implemented in [6] to obtain a 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. Other approaches as the Randall-Sundrum model [32] or 5-dimensional 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 selfgravitating spheres embedded in a 5-dimensional flat space-time (embedding class one). In this regard, it is clear that embedding of four-dimensional spacetimes into higher dimensions is a very useful tool in order to generate astrophysical interior models. Even more, models embedded in five dimensional spacetimes satisfy the so-called Karmarkar or class I condition [33] (for recent developments see, , for example) which allows to choose one of the metric functions as the one which generates the total solution. In this work we consider stellar interiors supported by anisotropic fluids fulfilling the polytropic equation of state for the radial pressure and additionally, we implement the class I condition to close the system and construct the corresponding generalized Lane-Emden equation. The manuscript is organized as follows. In the next section we review the main aspects of anisotropic polytropes and obtain the generalize Lane-Emden equation. Next, we introduce the class I condition and solve the Lane-Emden equation numerically. We conclude this section with a brief discussion about the incidence of the class I condition on the Tolman mass. The last section is devoted to final remarks and conclusions.
2 The polytrope for anisotropic fluid

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 is described by the energy-momentum tensor where ρ is the total energy density, is the four velocity of the fluid and s µ is defined as with the properties s µ u µ = 0, s µ s µ = −1.
Replacing (3), (5), (6) and (7) in (4), we have where primes denote derivative with respect to r. Furthermore, we shall consider that the fluid distribution is surrounded by the Schwarzschild vacuum, 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. From the radial component of the conservation law, 1 We are assuming natural units G = c = 1 3 one obtains the generalized Tolman-Oppenheimer-Volkoff equation for anisotropic matter which reads, Alternatively, using where the mass function m is as usually defined as 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 As already mentioned in the introduction, in order to integrate equation (19), we shall need additional information. In this work we use the polytropic equation of state for the radial pressure and the class I condition in order to close the system. In the next section we will implement the polytropic equation of state to construct the general Lane-Emden equation.

Relativistic polytrope for anisotropic fluids
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.
Let us first consider the case γ = 1. Now, defining 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), we may rewrite (22) as with P rc = Kρ γr c . Now, replacing (22) and (24) in (19) we obtain where α = P rc /ρ c . Let us now introduce the following dimensionless variables from where (26) reads which corresponds to the generalized Lane-Emden equation for this case.
Let us now consider the isothermal case: n = ±∞, γ = 1. In this case, the definition of the variable ψ becomes, and we can write equation (22) as Replacing (30) and (31) in TOV equation (19), we obtain Introducing dimensionless variables, Eq. (32) becomes which corresponds to the Lane-Emden equation for the isothermal case.

The class I condition
Embedding of four-dimensional spacetimes into higher dimensions is an invaluable tool in generating both cosmological and astrophysical models [41]. As it is well known, the Karmarkar condition is necessary and sufficient to ensure class one solutions [33] which for spherically symmetric space-times reads (R θφθφ = 0) [33] leading to with e λ = 1. Now, using (24), (25) and (37) in (20) we obtain from where, after using the set of dimensionless variables (27) and (28), we arrive at Finally, replacing (39) in (29), the class I generalized Lane-Emden equation for γ = 1 reads We can repeat the same procedure for the γ = 1 (isothermal) case. Again, after using (30), (31) and (37) in (20) we obtain, from where, by means of the definitions (33) and (34) we have Finally, introducing (42) in (35) we arrive to the class I generalized Lane-Emden equation for γ = 1 (isothermal case), To complement the discussion, we proceed to calculate the Tolman mass (for the γ = 1 case) [56], defined by or, alternatively [24] m T = e (ν+λ)/2 m + 4πr 3 P r .

5
Using (18) and (50) in (46) we arrive at where For the numerical calculations it is convenient to change to the following dimensionless variables in terms of which equation (51) reads, where G(x) becomes, This procedure, where the integration of the TOV equation is carried out, allows us to find the metric function ν (using the appropriate boundary conditions) and finally Tolman's mass, (56) and (57), is totally equivalent to the one developed in [57] for the different cases exposed. Now, the above results will be integrated numerically, with the corresponding boundary conditions, in order to obtain the solution for each case considered. Figure 1 shows the behaviour of e ν as a function of x for different values of the parameters involved. Note that, the function is monotonously increasing as required. In figures 2 and 3 it is shown the integration of Eqs. (43) (isothermal case) and (40), respectively, for the values of the parameters indicated in the figure legend. Note that for the isothermal case (γ = 1) the configurations are unbounded and as a consequence, it is meaningless to define a surface potential or the total mass. However, for the γ = 1 case ψ is monotonously decreasing as expected, and the radial pressure P r vanishes at the surface as required by the continuity of the second fundamental form. The parameter y ("the surface potential") which measures the degree of compactness is plotted in Fig. 4 as function of the polytropic index n for different values of α. Note that the plots have discontinuities in some point n d , where we observed that for n < n d , the solution is well-behaved, and for n ≥ n d , the numerical solution starts to have either imaginary components or strange oscillatory behavior. Figures 5 and 6 display the Tolman mass (normalized by the total mass), for the case γ = 1, as function of x for the selection of values of the parameters indicated in the legend. The behavior of the curves is qualitatively the same for a wide range of values of the parameters. We observe that as we move from the less compact configuration (curve a) to the more compact one (curve d), the Tolman mass tends to concentrate on the outer regions of the sphere. Therefore, it may be inferred from the figures 5 and 6 that more stable configurations correspond to smaller values of y since they are associated to a sharper reduction of the Tolman mass in the inner regions, as usually happens if the anisotropy is fixed by the imposed condition. This is the usual strategy adopted by the fluid distribution to keep the equilibrium; it tries to concentrate the Tolman mass in the outer regions. This behavior was already observed for different families of anisotropic polytropes discussed in [5]. The same result is obtained from the Fig. 7 where we plot the Tolman mass as a function of x for the selection of values of the parameters indicated in the legend.
For completeness, we plot the Tolman mass, as function of x, for several polytropic indexes and compactness parameters represented by y. The result i shown in Fig. 8 for the selection of values of the parameters indicated in the legend.
It is worth mentioning that all the numerical results show well-behaved solutions for small alpha values and a specific range of values for n.

Karmarkar condition and Tolman mass
In this section we would like to establish some relationships that may be useful to analyze the anisotropy inherited by the Karmarkar condition. In order to relate the influence of the Karmarkar condition with the mass function m and the Tolman mass, m T , we will take an instructive path to complete our analysis of the fluid distribution of the compact body. From (18) we obtain   the usual definition of the mass function, Note that, we can write condition (36), for our static and spherically symmetric case, in order to define a new variable, K, in the following form, where the Karmarkar condition is satisfied when K = 0. In this way we can keep track of how this condition influences our anisotropic distribution of fluid of the com- pact object. Using the components of the Riemann tensor appropriate to the symmetry we may obtain from (59) the following expression, which, after using Einstein equations (8), (9) and (10), can be written as Now, we proceed in two different ways in order to interpret the influence of the Karmarkar condition over our interior solution. First, from (8) and (9), we obtain the derivatives of the metric functions ν and λ in the Now, replacing these expressions in the last term of (61) we obtain 3 − 16π 2 r 4 T 0 0 T 1 1 = e −(ν+λ) K, which relates the mass function m with the K condition, defined in (60). In contrast, if we use (44) and Einstein equation (8), we can write the last term of (61) in terms of the Tolman mass function, m T , as Finally, using (45) in the previous equation we have which constitutes a relation between the mass function, the Tolman mass and the K condition. So, we have been successful in relating the components of the energymomentum tensor, the mass function and the Tolman mass with the defined K condition for the symmetry considered in our development. This could allow us to discover the influence over these mass functions by imposing the Karmarkar condition and then be able to find out how this metric (embedding) condition influences the distribution of our compact fluid. In order to illustrate the effect of the Tolman mass on K we shall consider a simple example, namely, an anisotropic solution with vanishing radial pressure [60], with 8 which ensure the fulfilling of the matching conditions (12), (13) and (14). Next, replacing (67), (68), (69), (70), (71), (72) and (73) in (66) and using (45) to express the mass function in terms of m T we arrive at In Fig. 9, we show the parametric plot of K(x) as a function of the normalized Tolman mass, η(x) in the interval x ∈ (0, 1) for different values of y. Interestingly, the solution deviates from the K = 0 condition as y increases.

Conclusions
In this work, we have developed a general method to construct locally anisotropic polytropes, and we have applied it to the specific case of a class I solution for sherically symetric interior space-times. As it occurs in the specific case of conformally flat polytrope spheres [6], our models are necessarily anisotropic (principal stresses unequal) and therefore are not continuously linked to isotropic polytropes (just like it ocurres in other reported anisotropic distributions [58]). Since the inclusion of pressure anisotropy implies an additional degree of freedom, the integration of the ensuing Lane-Emden equation, in the general case, requires additional information. Usually this information is provided through equations of state or other type of assumptions about the nature of the anisotropy. Here, we have supplied such additional information by assuming that the polytrope has the specific anisotropy derived from de Karmarkar condition. This represents an heuristic assumption whose validity will be confirmed (or not) from its application to specific problems. Then, the geometrical constraint, chosen to work with, produce a specific type of anisotropy leading to a Lane-Emden equation that required a numerical analysis, given the high nonlinearity of the equations.
We could differentiate two possible cases depending on whether γ = 1 or γ = 1. We have integrated the Lane-Emden equations for these cases, for a very large set of values of the parameters, even though, only a very specific set for each case is exhibited, this because the qualitative behavior of the system does not change much for a wide range of the parameters. The main reason to present such models was not to describe any specific astrophysical scenario, but to illustrate the applications of our approach. However, the protocol here developed could be useful for dealing with several phenomena in which anisotropic polytropes appear, with special emphasis in the physics of compact objects and so on, the obtained models exhibit some interesting features which deserve to be commented. We do not know if distributions that are usually represented as polytropes, and that also constitute very compact objects, like neutron stars (that are described by a Fermi gas), or Super-Chandrasekhar white dwarfs could be adequately modeled in this way. For such configurations it is evident that general relativistic effects as well as the inclusion of pressure anisotropy are unavoidable, and therefore its study could be carried on using the approach presented here. Here too, care must be exercised with the fact that some of the physical phenomena present in such configurations (e.g. very strong magnetic fields) could break the spherical symmetry.
It is interesting to highlight that polytropes with γ = 1 are used to model systems with non-degenerate isothermal cores which play an important role in the analysis of the Schönberg--Chandrasekhar limit. Any distribution of pressure anisotropy would affect the Schönberg-Chandrasekhar limit not only by affecting the structure of the polytrope, but also by the modifications of the virial theorem, introduced by the anisotropy.
In the last section we have followed a similar strategy previously used to express the Tolman mass in terms of traces of the Riemann tensor (matter sector) and the only non-vanishing scalar associated to the electric part of the Weyl tensor in static and spherically symmetric space-times, E (see [59], for example). In this sense, it is clear how the different parts affect the Tolman mass. For example, for conformally flat solutions the Weyl condition reads E = 0 which leads to a geometric constraint which allows to close the system. Similarly, here we have obtained expressions relating the scalar K with the traces of the Riemann tensor, the mass function m (64) and the Tolman mass (66) and class one solutions are obtained when the Karmarkar condition is fulfilled, namely K = 0. In this sense, expressions (64) and (66) may be useful to investigate the role played by Karmarkar embeding condition in our fluid distribution corresponding to the compact object. These aspects must be approached with the necessary care and detail in future works.
We do not know if there are any real physical phenomena that can produce a distribution of local anisotropy similar to the one inherited through the metric Karmarkar condition described in our models. So, it could be interesting to be able to provide a physical meaningful picture for a source of the anisotropy that we have obtained (by means of magnetic field or charge distribution) of the type one could expect to find in some astrophysical settings. Although it is certainly an interesting question, the answer lies beyond the scope of this work. However, we believe that our approach provides a useful tool to carry out this and other investigations.