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 (TOV) 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 behaviour 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 anisotropy factor has been discussed in 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 very high density regime (> 10 15 gr cm 3 ) and thus the nuclear interaction needs to be treated under 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.
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 ρ. The polytropic equation of state is one of the most important choices for a state equation which has a wide range of applications [12,13,14,15]. This equation of state is written as P r = Kρ γ , where γ = 1+ 1 n is the polytropic exponent and n is known as the polytropic index. For n = 2, the exact solutions of Einstein's field equations for an anisotropic sphere have been obtained by Thirukkanesh and Ragel [16].
The anisotropic model has been investigated for uniform matter density in [17] and for variable density distribution in [18,19].
Anisotropic pressure can affect on the stability of the stellar objects. Dev and Gleiser [2,20,21] have extended the variational formalism of Chandrasekhar [22] for isotropic star to study the stability of anisotropic general relativistic spheres against the 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 to obtain the exact solutions of general relativistic star equations in the presence of anisotropic pressure. Maharaj and Sharma [23] 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. Cylinderically symmetric models have been discussed in [24].
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. The investigation of these objects requires an exact solution of the Einstein-Maxwell system. Anisotropic fluid with a linear, quadratic and polytopic equation of state in the presence of electromagnetic field is studied in [25,26,27,28,29].
The importance of an equation of state in a stellar model has been emphasized by Varela et al [30] who provided a mechanism of dealing with charged anisotropic matter in a general approach. And also the exact solutions for a charged anisotropic quark star have been reported by Sunzu et al. [31].
Anisotropy can exist not only in ordinary compact stars but also in the hypothetical objects like Boson stars [32] and Gravastars [33]. Besides the isolated objects in astrophysics, anisotropy is present in globular clusters, galactic bulges and dark halos [34,35].
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 [36] to the anisotropic case. In section 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,16,37] 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 section 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 energy-momentum tensor for anisotropic matter are: T j k = 0 if j = k So the Einstein's field equations become: and the radial component of Bianchi identities reduces to the equation of hydrostatic equilibrium: In fact this equation can be obtained from equations (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 anisotropy factor which is explained in the next section.
To obtain a unique solution of the Einstein's 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 TOV equation of hydrostatic equilibrium for the spherically static aniotropic star from (3) : As it 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 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 implying an 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, equation (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 relation (7) at the surface of the star . This yields to: So we can write the metric coefficient ν(r) in terms of θ and ∆ as: Substituting e −κ from (9) and dν dr from (18) into equation (3), and using (10), the above equation reduces to: Now defining the dimensionless parameters ξ and η : , equations (10) and (23) take the following form: which are the general relativistic Lane-Emden equations for anisotropic polytropes. In the non-relativistic limit, according to (19), q 0 ∼ 0 and one can recover the Newtonian anisotropic Lane-Emden equation [3]. Also setting ∆ = 0, the isotropic Lane-Emden equation is obtained.
In equation (25) if ξ tends to zero, we will have: So we arrive at the following theorem: 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. 1 The plus sign holds if −1 < n < ∞, and the minus sign if −∞ < n < −1

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 mathemtaical relation between these functions in order to the set of equations be closed. Here we use the latter approach which is explained in Section 3.2. First in Section 3.1 we find a series expansion of θ(ξ), ∆(ξ) and η(ξ) near the center of anisotropic star. For isotropic case one can find the corresponding series expansion in [12] and the references in it. That is: In Section 3.2 we shall solve the Lane-Emden equations assuming the density is constant throughout the star. A similar calculation for isotropic star has been derived by Tooper [38]. In this case, the variables θ and η are separable and the exact analytical solutions of (25) and (26) for ∆ = 0, are [38]:

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 equation (28) ) are finite at the center of star. Thus if r → 0, η ≈ ξ 3 , ∆ ≈ ξ and θ ≈ 1.
Comparing equations (36)- (38) with (31) and (32), one can find that how the anisotropy factor affects the behaviour of quantities near the center of star.

An exact solution
Following [3], we assume that the presence of ∆ does not change the forms of isotropic Lane-Emden equations and it only modifies the coefficients of these equations. To do this, noting equation (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 equation (18) and integrating it, the metric coefficient ν(r) becomes: which is identical to the isotropic one [12] in which the coefficient n has changed to m. Using equations (9) and (40) and defining a dimensionless parameter ξ such that: (25) and (26) reduce to: Equations (42) and (43) are the extended form of relativistic Lane-Emden equations for anisotropic configurations, resulted by our ansatz describing above.
Here we want to obtain the analytical solution of anisotropic Lane-Emden equations for constant density. For n = 0, a solution of equations (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 solution for an anisotropic constant-density star. This solution can be considered as a well-defined solution if the physical quantities like the anisotropy factor be finite inside the star. It can be easily 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 Figure 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 Figure 2 for the same values of k and q 0 as in Figure 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. Now let us to discuss the validity of 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 same 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 star. In Figure 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 = 1 2 there is a violation of the dominant energy condition. Other quantities of physical interest are the stellar radius and mass. To find these, it is needed 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, relation (47) reduces to the corresponding radius in the case of isotropic star [12]. Now one can find the mass-radius relationship from (11): then according to relation (41) : and using the value of ξ R (45), one can find the mass of an anisotropic star: To get the mass-radius relationship, it is required to evaluate the integrals of (47) and (48) 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 mass-radius relation (51) as a function of ξ for specified values of k and q 0 .
Also the internal compactness of the star can be written as follows: and correspondingly the surface redshift z(r) is given by :

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 equation (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 polytropic index in general. The first case is studied in section (3.2) with details. For the second case using equation (18), one can introduce m such that: So the resulted Lane-Emden equations will be : where the plus sign holds if −1 < n+1 m < ∞ and the minus sign if −∞ < n+1 m < −1. In this case the dimensionless parameter ξ is introduced as r=( ± n+1 m q 0 4πGρ 0 ) 1/2 ξ. These equations for n = 0 yield the same solution (44) and (46) in which k = 3 4 m.