Study of polytropes with generalized polytropic equation of state

The aim of this paper is to discuss the theory of Newtonian and relativistic polytropes with a generalized polytropic equation of state. For this purpose, we formulated the general framework to discuss the physical properties of polytropes with an anisotropic inner fluid distribution under conformally flat condition in the presence of charge. We investigate the stability of these polytropes in the vicinity of a generalized polytropic equation through the Tolman mass. It is concluded that one of the derived models is physically acceptable.


Introduction
The theory of polytropes has a significant role to play in studying the inner structure of astrophysical compact objects (CO) very precisely. In Newtonian gravity, various useful physical phenomena have been addressed with polytropic equations of state (EoSs). Chandrasekhar [1] presented the basic theory of Newtonian polytropes emerging through the laws of thermodynamics for polytropic spheres. Tooper [2,3] provided the basic formalism of polytropes for a compressible fluid under the assumption of quasi-static equilibrium. He extended his work for an adiabatic fluid sphere and provided the fundamental framework to derive the Lane-Emden equation (LEe) for relativistic polytropes. Kovetz [4] redefined some anomalies in the theory of slowly rotating polytropes presented by Chandrasekhar [1]. Abramowicz [5] extended the general form of the LEe for spherical, planar, and cylindrical polytropes for higher dimensional spaces.
In general relativity (GR), polytropes have been discussed by many researchers by means of the LEe which can be derived from the hydrostatic equilibrium configuration of a e-mail: azammath@gmail.com; azam.math@ue.edu.pk b e-mail: syedalimardanazmi@yahoo.com; ali.azmi@umt.edu.pk c e-mail: ifra.noureen@gmail.com; ifra.noureen@umt.edu.pk d e-mail: aziz3037@yahoo.com; aziz.rehman@umt.edu.pk the CO. Cosenza et al. [6] presented a heuristic procedure to construct anisotropic models in GR. Herrera and Santos [7] provided a comprehensive study to discuss the possible causes for the existence of local anisotropy in self-gravitating systems both in the Newtonian and the relativistic regime. Herrera and Barreto [8] investigated relativistic polytropes under the assumption of a post-quasi-static regime and presented a new way to describe various physical variables like pressure, mass, and energy density by means of effective variables. They considered two possible configurations of polytropes in the framework of GR and found that only one was physically viable. Anisotropy plays a very vital role in the theory of GR applied to a discussion of spherical CO. Herrera et al. [9] developed a full set of governing equations for spherically symmetric dissipative fluids with anisotropic stresses to study self-gravitating systems in the framework of GR. Herrera and Barreto [10,11] used the Tolman mass (measure of active gravitational mass) to find the stability of both Newtonian and relativistic polytropes with an anisotropic inner matter configuration. Herrera et al. [12] analyzed in detail anisotropic polytropes with a conformally flat condition, which was useful in reducing the parameters involved in a relativistic modified LEe. Recently, Herrera et al. [13] discussed the stability of anisotropic polytropes by means of cracking.
The presence of charge on the stars is an important physical ingredient in studying the dynamics of astrophysical objects. To discuss the effect of charge on CO is always of great value in GR. Bekenstein [14] introduced the idea of hydrostatic equilibrium to observe gravitational collapse in charged CO. Bonnor [15,16] observed that electric repulsion can delay the gravitational collapse in spherically symmetric CO. Bondi has [17] provided a detailed explanation of the contraction in isotropic radiating CO with the help of Minkowski coordinates. Koppar et al. [18] provided a new technique to derive a charge generalization of a known static fluid solution for spherical symmetry. Ray et al. [19] examined the high density CO, which can hold an amount of charge of approximately 10 20 coulomb. Herrera et al. [20] analyzed charged dissipative spherical fluids. The physical meaning of structure scalars is analyzed for charged dissipative spherical fluids and for neutral dust in the presence of a cosmological constant. Takisa and Maharaj [21] studied the models of charged polytropic spherical symmetries. Sharif and Sadiq [22] presented a general formalism to study the effects of charge on anisotropic polytropes and developed a modified LEe. Azam et al. [23][24][25][26] analyzed different charged CO models to formulate their stability in the scenario of linear and quadratic regimes. It is found that the stability of these objects depends on the electromagnetic field as well as the choice of EoS.
In the development of the CO models the selection of the EoS is a very crucial issue. A polytropic EoS, P r = Kρ 1+ 1 n 0 , has been used by many researchers [5][6][7][8][9][10][11][12][13][14] for development and discussion of spherical CO. Chavanis [27] proposed the generalized polytropic EoS P r = α 1 ρ o + Kρ to discuss various cosmological aspects of the universe. He developed a model of the early universe to elaborate the transformation from a pre-radiation era to the radiation era for positive indices n > 0. In this continuity, he also produced models which described the late universe by considering negative indices in the case of a generalized polytropic EoS [28]. Freitas and Goncalves [29] used a generalized polytropic EoS to study primordial quantum fluctuations and build a universe with constant density at the origin. In this work, we will use a generalized polytropic EoS to discuss the theory of Newtonian polytropes and relativistic charged polytropes.
The plan of this paper is as follows. In Sect. 2, we shall provide the basics of the Einstein-Maxwell field equation and the hydrostatic equilibrium equation. In Sect. 3, Newtonian polytropes will be discussed and Sect. 4 is devoted to the observation of relativistic polytropes. Energy conditions and the conformally flat condition are discussed in Sect. 5. Section 6 is devoted to a stability analysis of polytropes. In the last section, we shall conclude and present our results.

The Einstein-Maxwell field equations
We consider a static spherically symmetric space-time where ν, λ are functions of r only. The energy-momentum tensor for an anisotropic matter distribution is given by where P t , P r , and ρ, respectively, represent the tangential pressure, the radial pressure, and the energy density of the inner fluid distribution. The four-velocity V i and four-vector S i satisfy the following equations: The electromagnetic energy-momentum tensor is defined by where is the Maxwell field tensor and it satisfies the following relations: where ψ i is the four-potential, μ 0 is the magnetic permeability, and J i is the four-current. Also, the four-potential and four-velocity satisfy the following relations in co-moving coordinates: where ψ is the scalar potential and σ is the charge density. The Maxwell equation (6) yields where the prime denotes differentiation with respect to r . From the above, we have where q(r ) = 4π r 0 μe λ 2 r 2 dr represents the total charge inside the sphere.
The Einstein-Maxwell field equations for the line element of Eq. (1) are given by and solving Eqs. (10)-(11) simultaneously leads to the hydrostatic equilibrium equation [22] d P r dr − 2 r + qq 8πr 3 + (4πr 4 P r − q 2 + mr) where we have used = (P t − P r ). The junction conditions are commonly used to relate inner and outer regions of the CO over a boundary surface. The choice of the outer geometry totally relies on the inner fluid distribution of the CO. For spherical symmetry, if the matter in the interior of a star is charged and anisotropic, we take the Reissner-Nordsträm space-time as the exterior geometry, Junction conditions play a very significant role in the study of the CO. These conditions determine whether the combination of two space-time metrics provides a physically viable solution or not when a hyper-surface divides the space-time into interior and exterior regions. The smooth matching of interior and exterior regions through the first and second fundamental forms [30][31][32][33] yields the following relations on the boundary surface: and the Misner-Sharp mass leads to [34,35] In spherical symmetry the Misner-Sharp mass can be described as a measure of the total energy within a sphere of a radius r at a time t. The theory of polytropes is based on the assumption of hydrostatic equilibrium and polytropic EoS. In this work, we discuss Newtonian and relativistic polytropes by using a generalized polytropic EoS, which is a combination of linear and polytropic EoSs.

Newtonian Polytrope
In the framework of Newtonian gravity, the polytropic EoSs are very useful in describing a huge variety of situations like inner pressure, types of inner fluid distribution etc. In this section, we formulate Newtonian polytropes through a hydrostatic equilibrium equation in the scenario of a generalized polytropic equation of state. Polytropes are supposed to be in hydrostatic equilibrium; we have to keep in mind that the appearance of cracking occurs when the system is taken out of equilibrium. We consider the following basic equations [8]: and in spherical coordinates the Poisson equation is given by here ρ o is the mass density and φ is taken to be the Newtonian gravitational potential. Equations (17)- (18), along with the generalized polytropic EoS [27,28] lead to a modified LEe (for γ = 1). Here K is the polytropic constant and n is the polytropic index. We take where ρ gc is the density evaluated at center c, and using Eqs. (20) in (17), we get From Eqs. (18), (20), and (21), we obtain where A 2 2 and A 2 3 are given by Inserting ξ = r A 2 in the above equation, we obtain the modified form of the LEe, with the boundary conditions The boundary of the sphere is defined by ξ = ξ n , such that θ(ξ n ) = 0.

The relativistic polytropes
This section deals with the relativistic configuration of polytropes with generalized EoS. The generalized polytropic EoS is the linear combination of the linear EoS 'P r = α 1 ρ o ' and the polytropic EoS 'P r = Kρ 3 ) matter. The polytropic part describes the cosmology of the early universe for positive values of the polytropic index, whereas it elaborates the late time universe with negative values of the polytropic index [27,28]. In order to discuss the cosmic behavior, ρ o was taken to be the Planck density, but for relativistic discussion we will consider it as the mass density and total energy density in case 1 and case 2, respectively. Here, we shall present the general formalism for relativistic polytropes with generalized polytropic EoS and investigate the presence of charge in the following two cases.

Case 1
Here, we consider the generalized polytropic EoS as so that the original polytropic part remains conserved; also the mass density ρ o is related to the total energy density ρ by [12] ρ = ρ o + n P r .
Now making the assumptions α = P rc ρ gc , α 2 = 1 + (n + 1)(α 1 + αθ), α 3 = (n + 1)α, where P rc is the pressure at center of the star, ρ gc is the mass density evaluated at the center of the CO, ξ , θ and v are dimensionless variables. Using the above assumptions along with the EoS (26), the hydrostatic equilibrium equation Eq.
Now differentiating Eq. (16) with respect to r and using the assumptions given in Eq. (28), we get Here, we consider the generalized polytropic EoS where the mass density ρ o is replaced by the total energy density ρ in Eq. (26), and they are related to each other by [12] We make the following assumptions: where c represents the quantity at the center of the star, α 2 , α 3 , and α 4 obey the same relations as in Eq. (28) with α defined in Eq. (33). Using the above assumptions along with the EoS (31), the hydrostatic equilibrium equation (13) becomes Now differentiating Eq. (16) with respect to r and using the assumptions given in Eq. (33), we get Equation (34)

Energy conditions and conformally flat condition
The energy conditions in GR are designed to obtain maximum possible information without enforcing a particular EoS. The energy conditions are provided in the sense that energy density cannot be negative because if it allows for random +ve and −ve energy regions, the empty space would become unstable. The energy conditions satisfied by all the spherically symmetric models are [22,36] ρ + q 2 8πr 4 > 0, For case 1 the conditions given in Eq. (36) turn out to be and for case 2 the conditions given in Eq. (36) emerge as We observe that the coupled equations (29)-(30) and Eqs. (34)-(35) form a system of differential equations. These systems involve three variables and we want some additional information to study a polytropic CO. We use the conformally flat condition to reduce one variable in the above said system of equations. The electric part of the Weyl tensor is related to the Weyl scalar given by [12,22] Now using the conformally flat condition, i.e., W = 0, along with the field equations, Eqs. (10)-(11), in Eq. (39), we get The above equation along with Eqs. (28) and (30) for case 1 yields Similarly, the anisotropy parameter for case 2 turns out to be Using Eqs. (41) in (29), we obtain the first equation of a coupled differential system corresponding to case 1, The above equation coupled with Eq. (30) provides a modified LEe [see the appendix, Eq. (55)] which describes the conformally flat polytropes for case 1.

Stability analysis
The stability of the model can be discussed by means of the Tolman mass, which measures the active gravitational mass of the CO. The important feature of the Tolman mass is that it can be evaluated by integrating over the region occupied by matter or electromagnetic energy [37][38][39]. The modified form of the Tolman mass for an anisotropic spherically symmetric metric (1) is given by [22] m T = M r 3 r 3 +r 3 r r e (λ+ν)/2 2 where represents the values calculated at the boundary of the CO. In order to get an expression for ν, we will solve the Einstein-Maxwell field equations (10)-(11) simultaneously; we have The integration of the above equation yields We define the dimensionless variables by For case 1, using Eqs. (28), (46), and (48) in (45), we get where y = (n+1)αν ξ and is given by For case 2, inserting Eqs. (33), (46), and (48) in (45), we obtain and is given by In order to study the physical viability of these models, we have calculated the Tolman mass whose behavior describes some physical features of these models, particularly the stability of the model. Figures 1, 2 and 3 have been plotted for case 1 (relativistic polytropes) corresponding to different parametric values n, α, y, Q, and α 1 = 0.5 [12,22]. Figure  1 shows that the Tolman mass is gradually decreasing and does not show any abnormal behavior as the value of α and charge Q increases. The curves become steeper with higher values of α and charge Q, but they remain positive, which is consistent with the results of Herrera et al. [12] for neutral polytropes. It is noted that for smaller values of y, the sphere gets more compact, which corresponds to an equilibrium configuration. This behavior is shown by the migration of the Tolman mass toward the boundary surface. In terms of stability, it happens due to a sharper reduction of the active gravitational mass in the inner regions of the sphere, which may correspond to more stable configurations for small values of y. The smooth behavior of Fig. 2 describes the solution of Eq. (30) which shows that the stability of the model can be enhanced through decreasing the values of y. Figure 3 represents the anisotropy of the model, which is larger near These results are consistent with [12,22] for charged and neutral polytropes.

Conclusion and discussion
In this work, we have developed the general framework to discuss the Newtonian and relativistic polytropes with a generalized polytropic EoS in the presence of an electromagnetic field under the conformally flat condition. The generalized polytropic EoS P r = α 1 ρ o +Kρ is the union of linear and polytropic EoSs. It is widely used in cosmology to explain different eras of the universe with the help of mathematical models. In the cosmic scenario ρ o was taken to be the Planck density but for a relativistic discussion we have taken it as a mass density and the total energy density in case 1 and case 2, respectively. In order to discuss the physical features of Newtonian and relativistic ploytropes, we have formulated the general formalism to obtain a modified LEe. The solutions of the LEe are called polytropes, which depend on the density profile function of the dimensionless radius ξ and the order of the solution is constrained by the value of the polytropic index n. The LEes are helpful to explain a relativistic CO as they produced a simple solution to describe the internal structure of the CO. But the cost of this simplicity is a power law relationship between the pressure and the density, which should be valid throughout the inner matter of the CO. The LEes have been initially developed for relativistic charged ploytropes with an anisotropic factor involved in them [see Eqs. (53) and (54)], which have three unknown variables. The conformally flat conditions are used to simplify these equations by eliminating anisotropy factor, which corresponds to a modified LEe [see Eqs. (55) and (56)], for charged polytropic spheres in the context of a generalized polytropic EoS. The energy conditions are helpful in order to check the physical viability of the CO in GR without enforcing an EoS. The energy conditions for both cases of relativistic polytropes have been developed in the presence of charge.
The stability of the relativistic polytropes is analyzed with the help of the Tolman mass. For case 1, the model behaved well as shown in Figs. 1, 2 and 3. The models developed here are anisotropic (in pressure) and helped us to discuss highly charged CO solutions. Figure 1 represents the physical behavior of the Tolman mass for different values of charge and model's parameters [12,22]. The Tolman mass, which is a measure of the active gravitational mass, gradually shifts toward the boundary surface but remains stable and does not present any fluctuations. Figure 1 represents the slow migration of the Tolman mass toward the boundary as the charge increases, corresponding to higher stability of the model. It is worth mentioning here that our results for the Tolman mass are considerably consistent with work of Herrera et al. [12] for polytropic EoSs. Figure 2 describes the solution of Eq. (30), which shows that the stability of the model can be enhanced through decreasing the values of y. This figure has been plotted for the indicated values of the triplet n, α, Q and shows a stable behavior for small values of α. Also, the configurations of the solution remain stable for an increasing value of charge, although the solution curves start bending toward the boundary of the CO, as shown in Fig. 2. It is also evident that the behavior of the parameter shows dispersal in the beginning, i.e., near the center of the CO, and it becomes more consistent as we move toward the boundary with the increase of charge. So, we can say that charge is the stabilizing factor here for relativistic polytropes with generalized polytropic EoS. Thus, the discussion of charge is essential for the study and exploration of relativistic polytropic COs as regards their physical viability.
The inclusion of anisotropy plays a vital role in the study of the astrophysical CO. The inner fluid distribution is assumed to be anisotropic for a spherically symmetric CO. Figure 3 represents the graphs of , which is the ratio of the anisotropy and the central mass density. The curves remain stable even after an increase of the charge but rapidly move toward the boundary of the sphere, and this phenomenon shows that the anisotropy is affected by the presence of an electromagnetic field. Hence, the presence of anisotropy with charge may lead toward gravitational collapse of the polytropic CO and can also be a source of the cracking phenomenon. The presence of cracking can be tested by means of methods provided by Herrera [40] and Gonzalez [41]. In this work, methods for the development of a generalized LEe have been presented.
We observe that in case 2 polytropes are not physically viable due to a violation of the energy conditions.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .