Core-envelope model of super dense star with distinct equation of states

The aim of the present paper is to study an anisotropic spherically symmetric core-envelope model of a super dense star in which core is equipped with linear equation of state, consistent with the quark matter while the envelope is considered to be of quadratic equation of state by adopting the philosophy of Takisa et al. (Pramana J Phys 92:40, 2019). We demonstrate that all the physical parameters are realistic within the core as well as envelope of the stellar object and continuous at the junction. Our model is shown to be physically viable and substantiate with the strange stars SAX J1808.4-3658 and 4U1608-52. Further, We infer that if the mass of the star increases then central density results to higher values and core shrinks, which justifies the dominating effect of gravity for higher mass celestial objects.


Introduction
The composition of a super dense stellar material determines various macroscopic features including moment of inertia, degree of compactness, approximate values of mass and radius of a star. In order to construct a stellar model, an equation of state (EOS); a relation between pressure and energy-density of the stellar material is usually considered. However, it is believed that the interior configuration of a star has a deconfined quark core environed by an envelope consisting of baryonic matter [1]. Therefore, to use a single EOS may not be a suitable choice for realistic modelling of the whole star. The demand of more input parameters, viz, central density and density at the boundary suggests employing of two different EOS to construct a realistic stellar model. a e-mail: pant_rp@rediffmail.com b e-mail: satya235@gmail.com c e-mail: ravindra.bisht@yahoo.com d e-mail: neeraj.pant@yahoo.com A tentative insight of a quark phase in the interior of a compact star has been suggested by Witten [2] and Farhi and Jaffe [3] in independent works. Their works opened a new dimension of an entirely new class of compact stars known as strange stars. The matter compositions of these strange stars are in the form of a quark-gluon-plasma (QGP) state with a typical mass about 1 − 2M and radius about 10 km.
The precise form of stellar matter distribution is unknown if the matter density of a spherical stellar object is greater than the nuclear density [4][5][6]. The dubiety in the selection of a suitable type of EOS for stellar objects outside of the nuclear regime, core-envelope model or hybrid model may be a good choice in the modelling of these high dense stellar objects [7][8][9][10][11].
Empirically, in order to design the core envelope model, the Darmois-Isreali conditions should be satisfied at the junction of the core and the envelope. With this empirical aspect, many authors [1,[12][13][14][15][16][17][18][19][20][21][22] have developed stellar models, with an assumption of a core and an envelope, for highly dense relativistic objects. Eventually in most of the core-envelope solutions so far obtained have the continuity of only metric potentials and pressure at the junction of the core and the envelope. Negi et al. [23,24] developed models considering the choices of core and envelope as Tolman IV and Tolman V solutions respectively for perfect fluid. In their solution, apart from the continuity of metric potentials, pressure, one more physical parameter, i.e., energy density was continuous at the junction. In recent past, many attempts have been made by various authors to find exact analytic solutions of the Einstein field system using linear EOS with MIT bag model [25][26][27][28][29][30][31] as well as quadratic EOS [32][33][34][35][36][37][38][39]. On the other hand, various authors [40][41][42][43][44][45] have also explored new solutions of the Einstein field equations for anisotropic fluid under the Karmarkar condition [46].
More recently, Takisa et al. [47] considered anisotropic fluid for core-envelope stellar model of PSR J1614-2230 with core layer having a quark matter distribution with lin-ear EOS and envelope layer consisting of matter distribution with quadratic EOS. They successfully verified the continuity of some of the physical parameters, viz., energy density, pressures, anisotropy, radial sound speed and adiabatic index.
By motivation of the pioneer work of [47], we explore an anisotropic spherically symmetric core-envelope model of a super dense star in which core is equipped with linear EOS consistent with the quark matter while the envelope is considered to be of quadratic EOS with Tolman VII type metric potential (g rr ) such that e λ is 1 at origin thus geometrically non-singular. By virtue of this all physical parameters (density, pressures, red-shift, compactfication factor, anisotropic constants, causality condition, adiabatic index, energy conditions) are realistic within the core as well as the envelope of the stellar object and regular at the junction (interface). Our model is shown to be physically viable and substantiate with the strange stars SAX J1808.4-3658 and 4U1608-52 [48].

A system of the Einstein field equations
The inside of an anisotropic fluid sphere is described by the following spherically symmetric line element in Schwarzschild coordinates (x i ) = (t, r, θ, φ): where ν(r ) and λ(r ) are known as the metric potentials.
Assuming the matter inside of the fluid sphere anisotropic, the Einstein field equations (for the units G = c = 1 ) are given as where is the energy-momentum tensor, R i j is the Ricci tensor, R denotes the scalar curvature, ρ, p r and p t are the energy density, radial pressure appraised in the direction of the spacelike vector and transverse pressure in the orthogonal direction to p r respectively. In comoving coordinates v i = 1 g tt δ i t is the 4-velocity normalized in such a way that g i j v i v j = 1 and χ j = − 1 g rr δ i r is the unit spacelike vector in the radial direction, i.e., g i j χ i χ j = −1. For the geometry and matter accounted by the line element (1) and energy momentum tensor (3), the EFEs generate the following system of equations where . denotes the derivative with respect to the radial coordinate r . Using Eqs. (5) and (6) we get the measure of anisotropy (Δ) as The force due to the pressure anisotropy is repulsive if p t > p r , and attractive if p t < p r [49]. The existence of outward force p t > p r , allows the building of more compact distribution when using anisotropic fluid than when using isotropic perfect fluid, p t = p r [50].
In view of the following transformations x = r 2 , z(x) = e −λ(r ) and y(x) = e ν(r ) , the system of equations (4-7) becomes where ( ) and ( ) represent first and second derivatives with respect to x.
In relativistic stellar objects, the inside matter distribution may be compiled of two regions: an inner layer called as core and an outer layer named as envelope with distinct pressures. In order to make a core-envelope model for a given star, it is mandatory to classify space-time into a number of discrete regions. These regions comprise of the core (0 ≤ r ≤ R C , Region C), the envelope (R c ≤ r ≤ R E , Region E) and the exterior (R E > r , Region B). The corresponding line elements for the three regions can be taken as The exterior of a star for the Region B is the Schwarzschild exterior solution which is given in Eq. (14).

Conditions for a physically realistic core-envelope model
In order to make the model physically viable, one needs to verify the following conditions in core, envelope and exterior regions (conditions developed by [47] are further augmented): (i) Geometrical regularity: The metric potentials and matter variables should be defined at the center and should be well behaved throughout within the star [47]. (ii) Density and pressures trends: The matter density ρ, radial pressure p r and transverse pressure p t the core and envelope of the star should be continuous at the junction, positive and monotonically decreasing outward. Further, the pressure-density ratios should be positive and less than 1 throughout within the star (Zeldovich's condition [51]) and continuous at the junction. (iii) Mass-radius relation, Red-shift and Compactification factor: The mass function m(r ), compactification parameter u(r ) and gravitational red shift z(r ) for the core and the envelope of the star should be continuous at the interface and increasing and decreasing respectively with the radial coordinate r . (iv) Anisotropic constant: The radial pressure should coincides with the tangential pressure at the core of the star and should be continuous at the interface besides being asymptotic and increasing outward. within the star and all the three forces should be continuous at the interface resulting the system to be in static equilibrium. (ix) At the stellar boundary p r (R E ) = 0 [47]. (x) The metric potentials of the core region should match smoothly with the gravitational potentials of the envelope region [47]. (xi) The gravitational potentials of the envelope layer should connected smoothly over the boundary with the Schwarzschild exterior metric [47]. The core-envelope model For core region (0 ≤r ≤ R C ), we choose Tolman VII type of g 11 for the gravitational potential e λ satisfying a linear EOS where a, b, α, and β are constants. Substituting z value from Eq. (15) in Eq. (11) and using Eqs. 8, 9, 15 and 16 , we obtain the following differential equation: On integrating Eq. (17), we get where P 1 = 1(1 + 3α), P 2 = −b(1 + 5α) and P 3 is an integration constant. In view of Eqs. (15) and (18) and using x = r 2 , the system of Eqs. (8)-(11) becomes where d 1 = 4(a − 2br 2 )(ar 2 − br 4 − 1) + P 2 1 r 2 , d 2 = P 2 r 2 (−6ar 2 + 4br 4 + P 2 r 4 + 8) +P 1 (−2ar 2 + 2P 2 r 4 + 4), d 3 = 32π(ar 2 − br 4 − 1), d 4 = 4b(−ar 2 + br 4 + 1) + P 2 1 , d 5 = 2P 1 (a − 2br 2 + P 2 r 2 ) + P 2 (−2ar 2 + P 2 r 4 + 4). The metric potentials e −λ C , e ν E and matter variables are continuous and well behaved in the core region.

Matching conditions (interface at the core-envelope and boundary)
Gravitational potential and radial pressure must be continuous at the interface and at the boundary (Darmois-Isreali conditions). Therefore, we arrive with the following two set of conditions:

Continuity at the boundary
The envelope metric (13) must be connected smoothly over the boundary with the Schwarzschild exterior solution (14) at r = R E i.e., where R E is the radius of the star. The six matching conditions (31)-(36) along with the twelve constants, namely a, b, P 7 , P 3 , A, B, C, R C , R E , M, α, β form an undetermined system of equations. Solving the above system of equations, the mass M and radius R E of the star are obtained as and the remaining constants P 7 , P 3 , C, β are of the form where The remaining seven constants a, b, c, α, A, B and R I are free parameters. These constants are selected in such a way that all the physical properties of the considered stellar objects are well-behaved. and e −λ are continuous at the interface and monotonically increasing and decreasing respectively with the radial coordinate r as well (Fig. 1).

Density and pressures trends
The matter density ρ, radial pressure p r and transverse pressure p t for the core and envelope of stars are continuous at the interface, positive and monotonically decreasing outward [53] (Figs. 2, 3 and 4). Further, the stars satisfy Zeldovich's condition [51] i.e. the pressure-density ratios are positive and less than 1 throughout within the stars and continuous at the junction (Fig. 5).

Mass-radius relation, red-shift and compactification factor
The mass function m(r ) and gravitational red-shift z(r ) for the core and envelope of stars SAX J1808.4-3658 and . Also, compactification parameter u(r ) for the above stars is continuous at the interface and increasing in nature with r ( Fig. 8) and lies within the Buchdahl limit [52].

Anisotropic constant
In Fig. 9, the radial pressure coincides with the tangential pressure at the core of stars and continuous at the interface and increasing outward [53].

Causal nature of radial sound speed
The radial sound speed of a compact star model satisfies the causality condition at the center and monotonically decreasing outward with the continuity at the interface. The profile of v 2 r of both core and envelope of the stars are given in Fig. 10.

Adiabatic index
For a relativistic anisotropic sphere the stability counts on the adiabatic index Γ r , the ratio of two specific heats, defined by [54], Bondi [55] suggested that for a stable Newtonian sphere, the Γ value should be greater than 4 3 . The profiles of adiabatic indexes of the core and envelope of both stars are plotted in Fig. 11, which shows the adiabatic indexes are continuous at the interface and satisfy the Bondi condition [55].

Energy conditions
For a physically stable configuration, the core and envelope of the stars should satisfy the following inequalities simultaneously (which are known as energy conditions [29]): (i) null energy condition ρ + p r ≥ 0 (NEC) (ii) weak energy conditions ρ + p r ≥ 0, ρ ≥ 0 (WEC r ) and ρ + p t ≥ 0, ρ ≥ 0 (WEC t ) and (iii) strong energy condition ρ + p r + 2 p t ≥ 0 (SEC). From Fig. 12, it is clearly visible that the variation of energy conditions with r of the core and the envelope of both the stars are continuous at the interface and satisfying realistic conditions.

TOV equation of core-envelope model
Equilibrium state under three forces, i.e., the resulting forces; gravitational (F g ), hydrostatic (F h ) and anisotropic (F a ) should be zero throughout within the star and continuous at the interface. The TOV equation [56] is given as where M g (r ) is the gravitational mass within radius r and can be calculated as from the Tolman-Whittaker formula and EFEs. Equation (43) is equivalent to the following balanced force equation where F g , F h and F a respectively are components of Eq. (43) above. From Fig. 13, we can visualize that the TOV condition is satisfied within the stars and all the three forces are continuous at the interface, thereby, concluding that the system is in static equilibrium.

Conclusion
In this paper, we have explored an anisotropic spherically symmetric core-envelope model of stars SAX J1808.4-3658 and 4U1608-52 in which core and envelope are considered with the linear and quadratic equation of states respectively. Our core-envelope model is physically viable and substantiate with the substantially strange stars (i) SAX J1808.4-3658 with mass M = 0.7M and radius R E = 6.2 km Table 1 Values of constants that generate masses, core and envelope radii (R C , R E ) for two strange stars b (km) −2 A (km) − Table 2 The variation of physical parameters, i.e., adiabatic index, density, red-shift for (i) Model I with mass M = 0.7M and radius R E = 6. In [47], some of the physical properties like anisotropy is not increasing throughout the region of stars and the radial pressure and anisotropy both are not differentiable at the interface due to a sharp jerk at the junction region. However, in this paper we have verified that all the physical and geometrical parameters, including anisotropic, radial pressure, compactfication factor, energy conditions,and stability conditions using TOV equation balancing forces are with the viable trend and continuous with the smooth variations throughout inside the stars. Further, from Tables 1 and 2, we infer that if the mass of the star increases then central density results to higher values and core shrinks, which justifies the dominating effect of gravity for higher mass celestial objects.