Cracking of Charged Polytropes with Generalized Polytropic Equation of State

We discuss the occurrence of cracking in charged anisotropic polytropes with generalized polytropic equation of state through two different assumptions; (i) by carrying out local density perturbations under conformally flat condition (ii) by perturbing anisotropy, polytropic index and charge parameters. For this purpose, we consider two different definitions of polytropes exist in literature. We conclude that under local density perturbations scheme cracking does not appears in both types of polytropes and stable configuration are observed, while with second kind of perturbation cracking appears in both types of polytropes under certain conditions.


Introduction
The theory of polytropes is vital in the evolution of mathematical models of compact objects and it attracts many researchers due to its simple form. In the study of polytropes, the main attraction is Lane-Emden equation, which led us towards the illustration of various astrophysical phenomena. Chandrasekhar [1] initially developed the theory of polytropes in Newtonian frame work with the help of laws of thermodynamics. Topper [2,3] used the hypothesis of quasi-static equilibrium form for the development of initial frame work of Newtonian polytropes. Kovetz [4] refined the work of Chandrasekhar [1] and reshaped the theory of polytropes. Abramowicz [5] was the first who presented higher dimension polytropes by developing Lane-Emden equation in higher dimension.
The study of electromagnetism and its effect on physical properties of astrophysical objects always fascinated the researches. Bekenstein [6] developed hydrostatic equilibrium equation (HEe) for the description of collapse of charged stars. Bonnor [7,8] presented the study on charged compact objects and described how electromagnetism affects the gravitational collapse. Bondi [9] used isotropic coordinates to analyzed the contraction of stars in the presence of charged. Ray et al. [10] examined the properties of stars with higher densities and concluded that approximately 10 20 coulomb charge can be hold by such stars. Herrera et al. [11] utilized structure scalars to illustrate compact objects having charged dissipative inner fluid distribution. Takisa and Maharaj [12] presented the mathematical model of charged compact objects with polytropic EoS.
The impact of anisotropy in the theory of general relativity is very important as we cannot study many physical phenomena without taking it into account. Cosenza et al. [13] presented a heuristic way for mathematical modeling of compact objects with anisotropic inner fluid distribution. Herrera and Santos [14] derived the anisotropic compact models in the frame work of general relativity. Herrera and Barreto [15] adopted a novel approach of effective variables for the description of physical variables involved in the anisotropic polytropic models. Herrera et al. [16] developed the governing equations in the presence of anisotropic stress for spherically symmetries. Herrera and Barreto [17,18] used the concept of Tolman-mass to check the viability of anisotropic polytropic models. Herrera et al. [19] adopted conformally flat approach to reduced physical parameters for simplification of Lane-Emden equations of polytropes.
The stability analysis of any developed model of stars plays a key role in general relativity. Any developed model cannot be used for the description of stars unless it is critically analyzed for stability. Bondi [20] developed HEe for stability analysis of neutral stars. Herrera et al. [21] proposed a novel way for the analysis of spherical symmetric models by means of cracking (overturning), which described the behavior of fluid distribution just after equilibrium state has been perturbed through density perturbation. Gonzalez et al. [22,23] provided an extension of Herrera et al. [21] by introducing local density perturbation (LDP). Azam et al. [24]- [28] used LDP for the analysis of various mathematical models of compact objects. Sharif and Sadiq [29] developed the model of charged polytropes. Azam et al. [30,31] developed the general frame work for charged polytropes with generalized polytropic equation of state (GPEoS) for spherical and cylindrical symmetries. They analyzed these models by means of Tolman-mass and Whittaker mass for spherical and cylindrical symmetries, respectively. Herrera et al. [32] have discussed the effect of small fluctuations of local anisotropy of pressure, and energy density on spherical polytropes. Sharif and Sadiq [33] have examined the effects of charge on spherical polytropes. Azam and Mardan [34] refined the work [33] for the analysis of charged polytropes.
The plan of work is as follows. In section 2, we provide some basic equations. Section 3 and 4 are devoted for the analysis of cracking through LDP and parametric perturbation respectively. In the last section we conclude over results.

Einstein-Maxwell Field Equations
We consider static spherically symmetric space time where ν(r) and λ(r) both depends only on radial coordinate r. The generalized form of energy-momentum tensor for charged anisotropic inner fluid distribution is given by where P t , P r , ρ, V i , S i and F mn represent the tangential pressure, radial pressure, energy density, four velocity, four vector and Maxwell field tensor for the inner fluid distribution. The Einstein-Maxwell field equations for line element Eq.(1) are given by Solving Eqs. (3)-(5) simultaneously lead to HEe where we have used ∆ = (P t − P r ). We take the Reissner-Nordsträm spacetime as the exterior geometry The junction conditions are very important in mathematical modeling of compact stars. They provide us the criterion for the collaboration of two metrics, which can results a physically viable solution [34,35]. For smooth matching of two space times, we must have and Misner-Sharp mass [36] leads to which has been used in the development of Lane-Emden equations [30].

Effect of Local Density Perturbation
In this section, we apply LDP [22,23] on charged conformally flat polytropes in equilibrium state. In LDP scheme it is assumed that all the physical parameters involved in the model and their derivative as function of density. Then the density is perturbed slightly and its effects has been observed on the HEe. Two different kinds of polytropes exist in literature will discussed here.

Case 1
We consider the GPEoS as so that the original polytropic part remain conserved, also the mass density ρ o is related to total energy density ρ as [7] Now taking following assumptions where P rc is the pressure at center of the star, ρ gc is the mass density evaluated at the center of CO, ξ, θ and v are dimensionless variables. We use conformally flat condition to find the expression of anisotropy factor ∆. The electric part of Weyl tensor is related to Weyl scalar given by [7,16] Using conformally flat condition, i.e., W = 0, along with Eqs.(3)-(5) in Eq.(13), we get Now differentiating Eq.(9) with respect to "r" and using the assumptions given in Eq.(12), we get The above equation along with Eq. (12) and (14) yields In order to observe the effects of LDP on conformally flat polytropes, we transform the HEe Eq.(6) by using Eqs.(10), (12) and (16) as Now we apply LDP to perturb all the physical variables in Eq. (17) and for this purpose we can write So the perturb form of the Eq. (17) can be written as where We will plot the force distribution δR 1 δρgc against the dimensionless radius ξ to observe possible occurrence of cracking. We say that cracking appears if force distribution changes it sign.

Case 2
Here, we consider the GPEoS as where mass density ρ o is replaced by total energy density ρ in Eq.(10) and they are related to each other as [7] We take following assumptions where c represents the quantity at center of the star, α 2 , α 3 and α 4 are same expressions as in Eq. (12) with α defined in Eq. (28). Carrying out the same process as in case 1, we get and the anisotropy factor ∆ comes out to be The HEe Eq.(6) will transform as proceeding in the same way, the perturb form of the Eq.(31) can be written as where We will plot the force distribution δR 2 δρgc against the dimensionless radius ξ to observe possible occurrence of cracking. We say that cracking appears if force distribution changes it sign.

Effect of Parametric Perturbation
In this section, we will study the stability of charged polytropes by perturbing the polytropic index and anisotropy factor. For this purpose we assume that our distribution satisfy the following relation where C is constant, producing the following form of HEe with h = 1 − 2C. Now we shall briefly review the main results for each case.

Conclusion and Discussion
In this work, we have applied two different perturbation schemes on two types of charged polytropes developed under the assumption of GPEoS [30]. In first perturbation scheme, we have used LDP scheme with conformally flat condition and sketched the force distribution function against radius of star. Here, we have assumed that all physical parameters involved in the model and their derivatives as function of central density. Figure 1-3 show the plot of force distribution function δR 1 δρgc against dimension less radius ξ. It is observed that for different values of central density and charge the system remain stable even after perturbation. These plots represent a picture of system just after perturbation. We observe rapid growth in the perturbed force near the center but it become consistent, stable and smooth as we move from center to outer surface of star. The zoom box inside the figure depicted perturbed forces near the center. It clearly indicated that the as the magnitude of charge increases gradually the peak of perturbed forces become less comparatively and it indicates that the presence charge stabilized the system after perturbation. Such behavior is observed in both types of polytropes under LDP scheme (see figs. 1-6).
In second perturbation scheme, we have perturbed the system through parameters involved in the model like charge, anisotropy and polytropic index. Such perturbation is carried out with the assumption that radial pressure remain same even after perturbation. Under this perturbation scheme, the perturbed force distribution δR 3 δn is plotted against radius of star (see 7-9). We get stable regions for small values of α as shown in 7 and 8. These plots also shows that the system is very sensitive towards the choice of parameters. Initially the perturbed force δR 3 δn strong behavior near the boundary of star (see Figure 7) but slightly change in parameters shifts the magnitude of force towards the center of star (see Figure 8). In figs. 7 and 8 stable configurations are observed. If the value of α is increased significantly then the system become unstable after perturbation and cracking (overturning) is observed. For smaller values of charge weak cracking near the center and strong overturning occurs in the outer region. For sufficiently high value of charge we observe weak overturning near the center and strong cracking in the outer regions as shown in fig. 9. Under parametric perturbation the second kind of polytropes did not show any stable region and remain unstable under different combinations of parametric values. Cracking and overturning is observed in polytropes of second type as shown in 10 and 11.
From the above discussion, it is concluded that charged polytropes developed under GPEoS remain stable if LDP scheme is applied under conformally flat condition. They are not sensitive towards the perturbation in central density and such behavior appeared in both types discussed in this work. When perturbation is carried out through parameters, then first kind of polytropes show stable behavior for small values of α in the presence of charge but cracking (overturning) also appeared in this case for larger values. Further the second kind of polytropes remain unstable under under parametric perturbation. Hence the proper choice of parameters is very crucial in the study of polytropic models.