Charged anisotropic compact star core-envelope model with polytropic core and linear envelope

This manuscript is related to the construction of relativistic core-envelope model for spherically symmetric charged anisotropic compact objects. The polytropic equation of state is considered for core, while it is linear in the case of envelope. We present that core, envelope and the Reissner Nordstro¨\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ddot{o}$$\end{document}m exterior regions of stars match smoothly. It has been verified that all physical parameters are well behaved in the core and envelope region for the compact stars SAX J1808.4-3658 and 4U1608-52. Various physical parameters inside star are discussed herein, non-singularity and continuity at the junction has been catered as well. Impact of charged compact object together with core-envelope model on the mass, radius and compactification factor is described by graphical representation in both core and envelop regions. The stability of the model is worked out with the help of Tolman–Oppenheimer–Volkoff equations and radial sound speed.


Introduction
The phenomenon of gravitational collapse depicts gravitational interaction within gravitating sources termed as compact objects. This phenomenon takes place as a consequence of unbalance in internal gravitational pull and outward drawn pressure stresses. The collapse of massive stars results in the formation of gravitating objects such as neutron stars, white dwarfs and black holes. These stellar residues are extremely dense in nature depicting strong gravitational interaction. Schwarzschild [1] described gravitational field of spherically symmetric compact objects and presented the solutions of field equations for internal and external spacetime geometry of star. Chandrasekhar and Milne [2] explained the phenomenon of gravitational collapse responsible for creation of a e-mail: syedalimardanazmi@yahoo.com b e-mail: ifra.noureen@gmail.com (corresponding author) c e-mail: S2019109012@umt.edu.pk compact stars. Oppenheimer and Volkoff [3] studied gravitational analysis of compact stars and discussed role of pressure stresses within gravitating sources.
Einstein's general theory of relativity is based on the study of spacetime curvature. Symmetries in spacetime describe it's geometrical properties. In general relativity (GR), spherical symmetry is mostly considered as one can discuss the characteristics and properties of gravitating objects more conveniently. Spherical symmetry may also be considered due to weak magnetic field in compact stars [4]. Moreover, it is believed that consideration of spherical symmetry is essential for many solutions of Einstein's field equations. Hillebrandt and Steinmetz [5] discussed the stability of spherically symmetric compact objects. Das et al. [6] constructed model of spherically symmetric gravitational equations by solving partial differential equations with initial value problem. Herrera et al. [7] presented an algorithm for spherically symmetric solutions of Einstein's field equations. Spherically symmetric relativistic compact star models were constructed in to analyze the phenomenon of gravitational collapse of stars [8,9].
Takisa et al. [10] described that in spherically symmetric compact stars the physical parameters like density, pressure, temperature and redshift are positive throughout the region. In the study of compact object the discussion of inner fluid is of primary importance. Different types of fluid configurations can used to explore various aspects of compact stars such as perfect fluid, isotropic fluid, anisotropic fluid etc. In isotropic fluid, radial and tangential pressures of compact object are same. The property which makes anisotropic fluid different from perfect fluid is that in anisotropic fluid, radial and tangential pressures are not equal. In compact objects anisotropy in fluid pressure may also exist because of the presence of solid core and occurrence of various physical processes in star. Dynamics and structural features of gravitating systems in anisotropic environment are explored in detail in [11][12][13][14]. Maharaj and Chaisi [15] presented the con-struction of anisotropic models by using solutions calculated from isotropic models. The study of Ivanov [16] is based on the characteristics of anisotropy and how anisotropy is assists in study of compact star models. Lakshmi [17] presented model of Einstein's field equation of anisotropic compact star by using line element with spherically symmetric coordinates and discussed density distribution in anisotropic fluid sphere. Sah and Prakas [18] discussed anisotropic fluid distribution in spherically symmetric compact star.
The core and envelope are two different but adjacent layers of a star. They hold different matter distributions and have different macroscopic physical features comprising of mass of inertia, a measure of density, evaluation of mass, or radius. The core envelope models for enormous relativistic stars have been studied in the past several years. Durgarpal and Gehlot [19] provided exact interior solutions with different energy densities for the core and envelope regions. Das and Narlikar [20] determined gravitational redshifts from the core-envelope boundary to an object's surface under various conditions and their effect on the equation of state (EoS). Negi et al. [21] discussed the construction of core-envelope models for varying physical parameters, in particular distinct EoS for core and envelope region. Sharma and Mukherjee [22] demonstrated the Vaidya−Tikekar model for a star with a quark phase in its inner core, surrounded by less compact baryonic envelope. Ramesh and Thomas [23] described the core region of the star for anisotropic matter configuration. Paul and Tikekar [24] remarked in their work that interior of compact star based on core region and envelope region. Ward and Whitworth [25] described the role of core region of star in last stage during evolution of compact star, the occurrence of fusion process in the core region binds lighter elements into heavier. Also, it preserve equilibrium between gravity and pressure of the star. The occupation of heavier elements in core region gets enlarged and size of star also increases. When surrounding pressure of envelope is reduces, the temperature of core region cools down resulting in dense compact objects. Core region of the star is surrounded by envelope region. Hansraj et al. [26] found solutions of nonlinear Riccati equations by constructing core-envelope models. Takisa et al. [27] discussed the envelope matter and construct anisotropic compact star model. Gedela et al. established core-envelope models for superdense compact stars [28].
In compact star core region is of quark matter and envelope region is baryonic matter. Weber et al. [29] discussed the properties of quark matter in detail. Annala et al. [30] discussed that on the bases of the behavior of compact matter the quark matter is composed of three states (i) up, (ii) down and (iii) strange quark. The densities of all these three states of quark matter is equal. Pure quark matter possess more stability then other nuclear matter. Quark masses have strong internal coupling which vary with the environment. The baryonic matter is made up of baryons which itself is made up of subatomic particles such as electron, proton and neutron. Its supposed that at finite baryon number the baryonic matter is quark matter. Naoz et al. [31] showed the formation of early structure of baryon. It is believed that about 68% of universe is made up of dark energy, 27% is dark matter and remaining 5% is of other matters. Baryonic matter is a form of dark matter. Thats why baryonic matter is helpful while dealing with cold and dense stars.
Inclusion of charge with matter anisotropy unfolds many intersecting features in the study of gravitational interaction within compact objects. Bonnor [32] presented the model of anisotropic sphere, in his study electrical repulsion affects self-gravitating charged matter. The gravitational field in compact stars is not exerted by the electrical energy of the charge of that star. In fact, energy of the electric field of charge object is used to balance the gravitational field. Bhatia et al. [33] described the creation of electric field of charge as a consequence of movement of electrons from the core of the compact star. Rao et al. [34] have studied effects of charge on the pressure of fluid, mass and redshift of compact objects. Ray et al. [35] discussed the high value of radius of charged supermassive compact objects effected by repulsive force and limiting value of amount of charge that a star can bear.
In comparison to the heavy ion particles, the lighter weight electron particles arise and move to the top of surface of star. So, to stop the rise and escape of electrons from the surface of star, the electric field is generated by charge. After the lose of some electrons from the surface, the star become stable because of the electric field generated by charge and net charge of star become positive. In charged stars, mass increases as an effect of presence of electric repulsive forces. Their exist relation between mass and central density, so central density increase with the increase in mass of charged star. The reason of increase in the value of redshift of charged compact star, is due to the emitted light rays from the surface of star that are shifted towards the red [10].
Sun et al. [36] described that gravitational collapse occur when inside pressure is not enough to balance gravity in stars. The existence of charge counterbalance the pressure gradient and avoid collapse. Pant et al. [37] described that charge present in fluid can keep singularity away from system. The validity and reliability of charged anisotropic compact stars is verified by evaluating its physical conditions such as stability, mass, radius and adiabatic index. Kiess [38] explained the relation between mass density, charge density, radius and central density of compact object. The study of Takisa et al. [39] showed that charge present in the compact star can change its whole appearance and structure. Panahi et al. [40] presented a comparison of charge anisotropy with isotropy. Tamta and Fuloria [41] obtained solutions for Einstein Maxwell's field equations for charged anisotropic fluid. Ratanpal and Bhar [42] constructed models on anisotropic compact stars in the presence of charge. Malaver and Manuel [43] constructed charged anisotropic model without singularities in matter which fulfil the physical properties required by realistic star.
Maurya [44] presented family of solutions for relativistic charged compact stars with anisotropic sphere. He discussed that at the center of compact star the electric charge is zero. Hence charge increases while moving away from center and attain maximum value at the surface of star. As a result of this type of distribution of charge at the center and surface, the core become less stable as compare to the surface of a charged object. The solution for charged anisotropy by taking conformal symmetry in spherically symmetric spacetime was discussed in [10]. Prasad and Kumar [45] explained solutions of Einstein's-Maxwell field equations with uniform charge distribution in compact stars by Pseudo-spheroidal spacetime. The same authors gave three energy conditions for charged fluid sphere are null energy conditions, weak and strong energy conditions. An EoS is an equation which explains the state of matter under various physical conditions. The EoS provides information about fluid type, solid matter, interior structure and variables as density and pressure of compact stars. To understand the structure and properties of compact object establishment of EoS is very significant. It relates physical parameters of the system in a single relation. It explains the dependence of pressure on density in a star. Zeldovich [46] presented the EoS with physical parameters like density and pressure. Barreto and Rojas [47] explained anisotropic fluid distribution of compact star by EoS. Chan and Ryong [48] explained the state of matter of self gravitating stars by using EoS. Compact star models with different state equations have been studied in [49][50][51][52][53][54].
All the properties and conditions of a strong gravitational system as compact star can't be explained just by a single EoS so we use various EoS's. It can be of the form of linear, polytropic, quadratic or of other dependence. Sometimes when solutions are known EoS's are selected without physical motivational reasons just to clarify differential equations of models. Tooper [55] explained that how polytropic EoS deals with equilibrium of system and also relate pressure and density in the sphere. Nilsson and Uggla [56] constructed model with polytropic EoS for compact stars. In polytropic equation, the density and pressure of compact stars are linked together by a power law. Negi and Durgapal [57] examined many EoS's and their physical parameters. Sharma and Maharaj [58] applied a linear EoS for compact stars with quark matter. Seveso et al. [59] defined the results of relativistic EoS extended by polytropic equation for compact stars.
In astrophysics, static gravitational equilibrium is defined by Tolman-Oppenheimer-Volkoff (TOV) equation. The TOV equation determines the properties and structure of spherically symmetric objects in equilibrium state. In 1939, it was first time when the issue of obtaining exact solutions for spherically symmetric fluid body was discussed. That discussion was started with the research of Tolman [60] and Oppenheimer and Volkoff [3]. We have discussed the viability of our developed model with the help of TOV equations.
The manuscript is arranged as follows: in Sect. 2, Einstein Maxwell field equations are constructed followed by Line elements for ore, envelope and Boundary region. Viability conditions for construction of core-envelope model are explained in Sect. 3, constituting subsections related to details of constructed model. Physical analysis for coreenvelope model and graphical representation of physical parameters is given in Sect. 4. In last section the results are summarized, followed by the list of references.

Einstein Maxwell field equations
Spherically symmetric spacetime for charged anisotropic fluid in Schwarzschild coordinates defined as ( is employed in this work, written as where ν(r ) and λ(r ) are the metric potentials. The Einstein's field equations are of the form In above equation, we have used gravitational units i.e., 8π G = c = 1. Ricci tensor is defined by R i j and R represents the scalar curvature. The energy-momentum tensor for charged anisotropic matter distribution is of the form In Eq. (3), v i and χ i stands for 4-velocity and radial 4-vector respectively and are given as The electromagnetic energy momentum tensor in Eq. (2) is given by where, F i j denotes the electromagnetic field tensor, which satisfies following relations where φ i and J i are the 4-potential and 4-current density respectively, and can be defined by the following relations where, = (r ), and ς is the charge density. Now, let us assume the total charge contained inside a sphere of radius r is denoted by S(r ), and defined as Now, by simplifying Maxwell Eqs. (6) and (7), we obtain following differential equation where 'prime' denotes the derivative with respect to radial coordinate r . Integrating Eq. (10) w.r.t 'r ', we have Onset of Einstein-Maxwell field equations is From Eqs. (13) and (14), we can obtain the following equation.
where = p t − p r is anisotropy factor, it vanishes at the origin (r = 0) and in the scenarios where p t = p r . The force due to anisotropic pressure of fluid is attractive if p t < p r , and repulsive if p t > p r [44]. Implementation of the transformation x = r 2 , z(x) = e −λ(r ) and y(x) = e ν(r ) , to the system of Eqs. (12)-(15) implies Here, 'z x , y x ' and 'y x x ' denotes first and second derivative with respect to x, respectively. In order to construct a core-envelope model for spherically symmetric spacetime, the matching of interior spacetime to exterior spacetime is required. Further, the interior of compact objects is of two regions, the inner layer is known as core and outer layer is known as envelope. The spacetime for core region (R C ) and envelope region (R E ) are taken as where, we match interior spacetime from Eqs. (22) and (23) to the exterior Reissner-Nordström spacetime (R B ) (we consider r = R at boundary, where R is radius of star). The exterior metric is given by where M and Q denotes the total mass and total charge of anisotropic charged star, respectively.

Viability conditions for core-envelope model
Physical viability demands that the model should meet the following requirements in core, envelope and exterior region (developed by [27] are discussed below) • Geometrical constancy: All physical parameters should be well behaved throughout within the star, in its center, in core and envelope [27]. • Viable trends in density and pressure: Matter density (ρ) and pressure stresses (radial and transverse pressure) shall be continuous at the matching, shall have positive values at the core and envelope of the compact object. Moreover, the pressure to density ratio should be con-tinuous at the matching and this ratio shall be positive everywhere within the star [46]. • Compactification factor u(r ), gravitational red shift z(r ) and mass-radius relation m(r ) shall be continuous in core and envelope region. Further, these parameters shall either be increasing or decreasing along radial coordinate. • The metric functions of the core region should match with the gravitational potential of the envelope region [27].

The core-envelope model
For core region, we assume an anisotropic type of wellknown Tolman VII [61] for metric potential e −λ C with polytropic EoS given by where n denotes polytropic index, a, b and α are constants. Herein, we assume polytropic index n = 1 in Eq. (26). The fundamental assumption of theory of polytropes is that the pressure stresses of gravitating sources are density dependent. Chandrasekhar [62,63] worked on polytropic gaseous sphere to account density and mass of white dwarfs. The same author provided the detailed framework on thermodynamical properties of compact objects via polytropic EoS [64]. Tooper [65] discussed relativistic polytropes and constructed the general framework for the derivation of LEe. Kaplan and Lupanov [66] worked structure of polytropic spheres and discussed its stability. Kaufmann [67] constricted relationship in mass and radius for spherical polytropes by taking different values of polytropic index. Occhionero [68] and Kovetz [69] worked on the structure of rotating polytropes to ascertain refinements in the theory of polytropes presented by Chandrasekhar [64]. Eriguchi [70,71] explored hydrostatic equilibrium of polytropes by using numerical computations extended his method for larger values of poly-tropic index i.e., n = 4, 4.9. Sharma [72] approximated the solution to field equations for spherical polytropes by making use of Padé approximation. Moreover, theory of relativistic polytropes and impact of polytopic EoS in different physical situations is discussed in [73,74].
By substituting value of S 2 and using Eqs. (17), and (25) in Eq. (26), we get following equations as follows where Here P 1 = a 1 + 1 r and P 2 = −b 1 + 2 r are integrating constants.
Substitution of x = r 2 and use of Eqs. (25), (27) and (25) in the Eqs. (17)-(4) leads to the following equations: where For realistic modeling a single EoS may not be appropriate, since we are studying the structure of star in two layers. That is why two different EoS are employed for the discussion of the stellar model. We have chosen an EoS which linearly relates the energy density and radial pressure in the envelope. Linear EoS is consistent in the modeling of compact objects, one can achieve modeling of anisotropic spherical objects with quark matter distributions [74,75]. It is suitable in retrieving the mass-radius relationship and values of surface redshifts corresponding to realistic stars [76,77]. For envelope region, we assume the same anisotropic Tolman VII [61] with linear EoS where a, b, α and β are constants. Substituting values from Eqs. (33) and (17) in Eq. (34), we have Here making use of Eqs. (33), (34) and x = r 2 in system of Eqs. (17)-(4), we get following equations where

Matching at the core-envelope and boundary
The line elements given in Eqs. (22) and (23) of core and envelope region must match at boundary r = R C . Metric potential and radial pressure must also be continuous at boundary.

Matching conditions of interface metrics
The structure of numerous compact stars is made up of core and envelope layers with different pressures. To build such star model, following matching conditions are important.

Matching conditions of envelope and exterior
The line elements Eqs. (23) and (24) of envelope and exterior region must match at boundary r = R E (where R E is the radius of star).
In matching conditions (41)- (46) constants are M, R E , C, β and σ , defined as Here, M and R E is mass and radius of star, respectively. where The constants a, b, α, A and B are free parameters. The values of these constants are taken in such a way that all physical properties of assumed compact objects are well-behaved.

Physical analysis for core-envelope model
The viability conditions for well behaved core-envelope model are explained in following subsections. The realistic stars SAX J1808.4-3658 and 4U1608-52 are considered to show physical acceptability of developed model. Following table provides values of parameters that are calculated to arrive at masses and radii of core and envelope for 4U1608-52.

Geometrical nonsingularity
The metric potentials are stated as e ν = positiveconstant and e −λ = 1 at the center of star. According to this e ν and e −λ are invariable and non-singular at the center of star. Additionally, e ν and e −λ are continuous at the matching and increase or decrease about radial coordinate r .

Density and pressure trends
The fluid density (ρ) and pressures ( p r , p t ) for coreenvelope star model are continuous at the matching, having positive values throughout the region [78] (Figs. 1, 2, 3). Additionally, the Zeldovich's condition [46] that the ratios of pressure and density are positive and smaller then 1 in the star and continuous at matching (Fig. 4).

Mass-radius relation, red-shift and compactification factor
Earman [79] described the objectives of studying red-shift of compact stars in GR. The mass (m(r )) and red-shift (z(r )) are continuous at the matching and increases or decreases respectively with radial coordinate r for the compact stars S AX J 1808.4 − 3658 and 4U 1608 − 52 (Figs. 5, 6). Additionally, compactification factor (u(r )) is continuous at the matching and increase w.r.t r (Fig. 7) and be within Buchdahl limit [80].

Anisotropic constant
The plot of anisotropy factor is continuous as shown in (Fig. 8).

Causality conditions of radial sound speed
The radial sound speed fulfill the causal conditions of compact star and decrease monotonically at the exterior and continuous at matching. The v 2 r of core-envelope model of star is in Fig. 9.

Adiabatic index
In relativistic anisotropic sphere the adiabatic index r and the proportion of two specific heats take account of stability of stars, defined by [81], For a stable Newtonian sphere > 4 3 [62][63][64]. For coreenvelope model of stars adiabatic indexes are plotted in Fig. 10. The Fig. 10 presents the continuity of adiabatic index at matching.

Energy conditions
For a physically stable arrangement, the core and envelope of stars fulfill some inequalities called as energy conditions [82] are following: (i) Null energy condition (ii) Weak energy conditions for radial and tangential pressures (iii) Strong energy condition where, E = S r 2 . The change in energy conditions w.r.t r is continuous at matching and fulfill realistic conditions for core and envelope of star (Fig. 11).

TOV equation of core-envelope model
For static equilibrium in compact stars we assume the three force components of TOV equation. The three resulting forces of equilibrium state are; gravitational force (F g ), hydrostatic force (F h ) and anisotropic force (F a ) should be zero everywhere inside star and continuous at matching. The TOV equation [83] is defined as Here, gravitational mass is function or r is denoted by M g (r ) and defined as through the formula Tolman-Whittaker and Einstein's field equations.
Following is the equilibrium force equation which is equivalent to Eq. (57) Here, elements of Eq. (57) are F g , F h and F a . The TOV equation is fulfilled inside the star and forces (F g , F h and F a ) are continuous at matching and system is in equilibrium (Fig. 12).

Summary and conclusion
In this paper we have constructed a core-envelope model by taking a spherically symmetric compact object with anisotropic charged matter configuration. Polytropic EoS is chosen for inner layer termed as core region while a linear EoS is employed to discuss envelope region of compact stars. It is worth mentioning here that, the spacetime geometry of core, envelope and exterior region of star match smoothly. In order to develop a physically viable core-envelop model, we have accounted the values of masses, radii and compactification factors for two stars namely S AX J 1808.4 − 3658 and J 1808.4 − 3658 by incorporating electromagnetic effects. The viability conditions for core-envelop model are dis- Fig. 12 Variation of balancing forces with radial coordinate Table 1 Calculated values of parameters that give core radii (R C , mass and envelope radii R E ) for S AX J 1808.4-3658   Tables 1 and 2. Herein, we shall summarize the results determined from graphical representation of physical parameters and required conditions as follows: • The density profile and pressure components of both stars under consideration are displayed in Figs. 1, 2 and 3. It can be clear that these physical parameters are continuous at the junction and positively defined. • Graphical representation of pressure to density ratios given in Fig. 4 shows continuity at the matching and this ratio is positive everywhere within the stars. • The mass function is displayed in Fig. 5, that shows increasing trend in mass from center to the exterior boundary. • Variation of gravitational red-shift with radial coordinate is expressed in Fig. 6 for both stars that shows positive and monotonically increasing behavior from center to exterior boundary. • The compactification factor, i.e. the mass-to-radius ratio of a compact star is calculated by using the formula, u = M R . The compactification factor can be used to classify the compact objects in various categories, when u = 1 2 , it represents a black hole. From Fig. 7, it can bee seen that compactification factor increases with the radial coordinate r .
• The anisotropy factor has been plotted in Fig. 8, it is evident from the graph that anisotropy factor has minimum value near the center and it increases gradually as we move from the center towards boundary. • The radial sound speed v 2 r satisfies the casuality condition, i.e., continuous at the junction and decreases with r as shown in Fig. 9.
• The adiabatic index defined as r = (ρ + p r )/ p r ∂ p r /∂ρ, depicts variation in anisotropic pressure density with a given change in density. The adiabatic shall have value in range > 4 3 for stable stellar configuration, that is achieved in Fig. 10.
• For a physically stable model, the star's core and envelope must satisfy the energy conditions. It is evident that at the junction, the variation in energy conditions with r of the core and envelope is satisfied as given in Fig. 11. • In general relativity, the TOV equation constrains the structure of a spherically symmetric body in static gravitational equilibrium. From Fig. 12, the TOV condition is satisfied and the three balancing forces are continuous at the junction. This concludes that the system is in static equilibrium state.
From the Tables 3 and 4, we can see that by increasing the mass of star its central density also increase. In this paper, we verify all physical properties of stars by equilibrium forces of TOV equations.
Viability of above mentioned conditions represent a physically acceptable core envelope model.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no external data associated with the manuscript.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .