On the relativistic anisotropic configurations

In this paper we study anisotropic spherical polytropes within the framework of general relativity. Using the anisotropic Tolman–Oppenheimer–Volkov equations, we explore the relativistic anisotropic Lane–Emden equations. We find how the anisotropic pressure affects the boundary conditions of these equations. Also we argue that the behavior of physical quantities near the center of star changes in the presence of anisotropy. For constant density, a class of exact solution is derived with the aid of a new ansatz and its physical properties are discussed.


Introduction
In theoretical astrophysics there is a growing interest to discuss the stellar structures in which the matter content is an anisotropic fluid. The effect of anisotropy can be studied both in Newtonian gravity and general relativity. For the configurations with not extremely high densities, the presence of an anisotropy factor has been discussed in the Newtonian regime [1][2][3] and otherwise general relativity must be used.
In 1972 Ruderman [4] for the first time argued that the nuclear matter may have anisotropic features in the very high density regime (>10 15 gr cm 3 ) and thus the nuclear interaction needs to be treated from a relativistic point of view. After the pioneering work of Bowers and Liang [5] in 1974, the study of anisotropic distribution of matter got wide attention. Anisotropy arises due to the existence of solid stellar core or by the presence of type-IIIA superfluid [6,7], phase transitions, pion condensation in a star [8], electromagnetic field [9][10][11], rotation, etc. For an exhaustive review on the anisotropic pressure see [12], and a comprehensive study of this topic may be found in [13].
Most studies assume that the matter of the anisotropic star is described by an equation of state relating the radial pressure P r or transversal pressure P ⊥ to the energy density ρ. a e-mail: fshojai@ut.ac.ir The polytropic equation of state is one of the most important choices for a state equation which has a wide range of applications [14][15][16][17]. There are two possible cases of the polytropic equation of state, namely: the radial pressure is proportional to the baryonic density, or it is proportional to the total energy density. Here we consider the latter case, P r = Kρ γ , where γ = 1 + 1 n is the polytropic exponent, ρ is the total energy density and n to the energy density For n = 2, the exact solutions of Einstein's field equations for an anisotropic sphere have been obtained by Thirukkanesh and Ragel [18].
The anisotropic model has been investigated for uniform matter density in [19] and for variable density distribution in [20,21].
Anisotropic pressure can affect the stability of the stellar objects. Following the Bondi formalism, the stability of a locally anisotropic fluid was first studied in [22].A general study of this issue, following the Chandrasekhar formalism [23], was reported in [24][25][26]. Dev and Gleiser [2,27,28] have extended the variational formalism of Chandrasekhar for an isotropic star to study the stability of anisotropic general relativistic spheres against radial perturbations. Assuming some specific density profiles and the adiabatic exponent, they have found that there can exist stable relativistic anisotropic spheres leading to instability in isotropic spheres.
Here our main purpose is to obtain the exact solutions of general relativistic star equations in the presence of anisotropic pressure. Maharaj and Sharma [29] have obtained some exact solutions of compact anisotropic stars analytically by assuming a linear equation of state with an appropriate ansatz.
The investigation of anisotropic configurations was not restricted to static spherically symmetric cases. All possible static cylindrical symmetric solutions are analyzed in [30], satisfying the conformal flatness condition. Also in [31], the stationary case is studied in detail. A general study of this issue has been presented in [32].
Also the modeling of dense charged gravitating objects in strong gravitational fields has attracted much interest in recent years because of its relevance to relativistic astrophysics [33]. The investigation of these objects requires an exact solution of the Einstein-Maxwell system. An anisotropic fluid with a linear, quadratic, and polytropic equation of state in the presence of an electromagnetic field is studied in [34][35][36][37][38].
The importance of an equation of state in a stellar model has been emphasized by Varela et al. [39] who provided a mechanism of dealing with charged anisotropic matter in a general approach. Also the exact solutions for a charged anisotropic quark star have been reported by Sunzu et al. [40].
Anisotropy can exist not only in ordinary compact stars but also in hypothetical objects like boson stars [41] and gravastars [42]. Besides the isolated objects in astrophysics, anisotropy is present in globular clusters, galactic bulges, and dark halos [43,44].
In this paper we want to discuss the anisotropic polytropes in the context of general relativity. In the next section we derive the relativistic Lane-Emden equations in the presence of anisotropy factor. Then we obtain their boundary conditions and generalize the Chandrasekhar's theorem [45] to the anisotropic case. In Sect. 3, first we find a series expansion for the dimensionless Lane-Emden function, θ , and the anisotropy factor, , near the center of star. Then in order to obtain the exact solutions, we have assumed that the presence of anisotropy factor does not change the general form of the relativistic hydrostatic equilibrium equations. This ansatz leads us to find a new class of exact solutions for the anisotropic Lane-Emden equations for constant density. This procedure is different from the usual approach [2,18,46] which is to suppose an ansatz on the form of anisotropy factor and try to get the physical solution of the equations. At the end in Sect. 4, we present a summary of the main conclusions.

Relativistic anisotropic polytropes
Since we deal with the static and spherically symmetric configurations, we use the following metric: The components of the energy-momentum tensor for anisotropic matter are So the Einstein field equations become and the radial component of the Bianchi identities reduces to the equation of hydrostatic equilibrium: In fact this equation can be obtained from Eq. (3)-(5) and usually it is convenient to use it instead of (5). The set of equations (3), (4), and (6) contain five unknown functions, κ(r ), ν(r ), P r (r ), P ⊥ (r ), and ρ(r ), and we have to add two more equations to make the set of equations closed. To do this, we use the polytropic equation of state and introduce a new ansatz for the anisotropy factor, which is explained in the next section.
To obtain a unique solution of the Einstein equations, one must also specify the boundary conditions. Outside the star, the spacetime geometry is described by the Schwarzschild metric. Thus we have where R is the radius of star and M is its mass. By integrating (4) with initial condition κ(0) = 0 (in order to avoid a singularity, κ must tend to zero at least as r 2 for r → 0), we get By defining the auxiliary function Equation (4) becomes where M(r ) is the relativistic mass inside the radius r : By using (6) and (9), we can obtain the Tolman-Oppenheimer-Volkov (TOV) equation of hydrostatic equilibrium for the spherically static anisotropic star from (3): As is obvious, vanishing recovers the TOV equation for isotropic star. Moreover, in the non-relativistic limit (c → ∞), one can easily obtain the anisotropic Jeans equation for the radial pressure: The boundary conditions of (12) are These conditions imply that the radial pressure is finite at the origin and vanishes at the surface of star. The mass function M(r ) goes to zero at the origin and reduces to the total mass at the surface. For a specified and implicating the equation of state we obtain a closed system of equations, (11), (12), and (16). Then the metric coefficients are given by (6) and (10). As we have mentioned before, for many astrophysical applications, assuming the polytropic equation of state enables one to make an approximate model of a star. One can always parametrize this equation by introducing the Lane-Emden function defined by where ρ 0 and P 0 are the central values of relativistic density and pressure. With respect to this new function, Eq. (6) becomes where This equation yields where the integration constant ν 0 is the central value of ν. To find it, we use the fact that the value of the above expression is given by Eq. (7) at the surface of the star. This yields So we can write the metric coefficient ν(r ) in terms of θ and as .

Theorem The finite solutions of relativistic anisotropic Lane-Emden equations at the origin have to satisfy the following relations:
This is the generalization of Theorem 1 in [3] for the relativistic anisotropic case.

Analytical solutions of relativistic Lane-Emden equation
In this section we shall discuss how one can solve the Lane-Emden equations (25), (26) with the boundary conditions (28). These are a system of two equations for three unknowns, the function θ(ξ ), (ξ ), and η(ξ ). To solve these equations, we need another equation. One can assume either a special form for one of the above functions or a mathematical relation between these functions in order for the set of equations to be closed. Here we use the latter approach, which is explained in Sect. 3.2. First in Sect. 3.1 we find a series expansion of θ(ξ ), (ξ ), and η(ξ ) near the center of anisotropic star. For the isotropic case one can find the corresponding series expansion in [14] and the references in it. That is, where In Sect. 3.2 we shall solve the Lane-Emden equations assuming the density is constant throughout the star. A similar calculation for an isotropic star has been derived by Tooper [47]. In this case, the variables θ and η are separable and the exact analytical solutions of (25) and (26) for = 0, are [47]

Series expansion
Near the center of star, we assume that Inserting these expressions into the Lane-Emden equations (25) and (26), we get provided that the pressure and density gradient (and also (0) according to Eq. (28)) are finite at the center of star. Thus if r → 0, η ≈ ξ 3 , ≈ ξ , and θ ≈ 1. Comparing Eqs. (36)- (38) with (31) and (32), one can find how the anisotropy factor affects the behavior of quantities near the center of star.

An exact solution
Following [3], we assume that the presence of does not change the forms of the isotropic Lane-Emden equations and it only modifies the coefficients of these equations. To do this, noting Eq. (18), we introduce a constant m such that (hereafter we set c = 1) 2q 0 (n + 1)dθ − 4 rρ 0 θ n dr = 2q 0 (m + 1)dθ. (39) Substituting this into Eq. (18) and integrating it, the metric coefficient ν(r ) becomes which is identical to the isotropic one [14] in which the coefficient n has changed to m. Using Eqs. (9) and (40) and defining a dimensionless parameter ξ such that

. Equations (25) and (26) reduce to
Equations (42) and (43) are the extended forms of the relativistic Lane-Emden equations for anisotropic configurations, resulting by our ansatz described above.
Here we want to obtain the analytical solution of anisotropic Lane-Emden equations for constant density. For n = 0, a solution of Eqs. (42) and (43) satisfying the initial conditions (28) is in which k = 3 4(m+1) is a natural number and q 0 > 0. This solution vanishes at which is the dimensionless value of the radial coordinate at the boundary. Thus we find that With suitable choices of k and q 0 , the equations of (44) and (46) present a new class of solutions for an anisotropic constant-density star. This solution can be considered as a well-defined solution if the physical quantities like the anisotropy factor are finite inside the star. It can easily be shown that for all values of q 0 and k, the anisotropy factor (46) is finite. Also if q 0 ≥ 1, the anisotropy factor has a local maximum inside the star but this is not true for q 0 < 1. Here in Fig. 1 we have plotted | | P 0 as a function of ξ for some values of q 0 and k. For any choice, the plots are from the center to the corresponding radius of the star. The radial pressure is plotted in Fig. 2 for the same values of k and q 0 as in Fig. 1. Now it becomes clear that it is a positive regular function everywhere inside the star. It starts with its local maximum value at the center of the star and gradually decreases until it vanishes at the surface of star.
In [48], it is shown that all static anisotropic spherically symmetric solutions of Einstein's equations can be produced by two generating functions. One of them is proportional to the anisotropy factor, (ξ ) = −8π G (ξ ), and the other, z(ξ ), is defined by Using (40) and (46), the generating functions of our solution are Now let us to discuss the validity of the energy conditions. Since q 0 is positive, the weak, null, and strong energy conditions hold for every positive value of central pressure. A simple calculation shows that for cases with 0 < q 0 ≤ 1, the dominant energy condition is satisfied in the same way as the other conditions. But for cases with q 0 > 1, this condition is not satisfied within a radius of 2k 3 k } from the center of the star. In Fig. 3, we have plotted the dominant energy condition as a function of ξ for different values of q 0 and k. While the dashed lines in each cases represent − 1 q 0 and 1 q 0 , it is clear that for q 0 = 3 2 there is a violation of the dominant energy condition.
Other quantities of physical interest are the stellar radius and mass. To find these, one needs to find the radius of the star, R , measured by an external observer. This radius for an anisotropic polytropic sphere with constant density is For m = 0, Eq. (50) reduces to the corresponding radius in the case of isotropic star [14]. Now one can find the massradius relationship from (11): and using the value of ξ R (45), one can find the mass of an anisotropic star: To get the mass-radius relationship, one is required to evaluate the integrals of (50) and (51) with arbitrary radial coordinate which now appears as the upper limit of integrals. We thus get M(r ) r = 128(ρ 0 π G) 3 27 Here in order to plot the mass-radius relation we replace both M(r ) and r with their dimensionless parameters η(ξ ) and ξ . Figure 4 shows the dimensionless form of the massradius relation (54) as a function of ξ for specified values of k and q 0 .
Also the internal compactness of the star can be written as follows:

Concluding remarks
In this paper, we have discussed the relativistic Lane-Emden equations and their boundary conditions for an anisotropic star. We made a new proposal to find the exact solution of these equations for constant density. Noting Eq. (18), we assume that the effect of anisotropy is changing the coefficient of the first term from 2q 0 (n + 1) to 2q 0 (m + 1) or alternatively changing the coefficient of the third term from (q 0 θ + 1) to m(q 0 θ + 1) where m is different from the polytropic index in general. The first case is studied in Sect.