General relativistic model for mixed fluid sphere with equation of state

We generate a general frame work to solve the Einstein system with an equation of state that describe static spherically symmetric anisotropic matter distribution in terms of a generating function. It is examined for a Van der Waals type equation of state with a physically reasonable form of generating function. The model satisfies all the required major physical properties of a realistic star. It is shown to be stable in the low-density regime that may represent a liquid–gas mixed fluid sphere.


Introduction
Many astrophysical and nuclear physics phenomena are of interest in self-gravitating stellar objects when its structure of matter approaches densities ρ ∼ 0.16 fm −3 at finite temperature (< 20 MeV). Exact solutions of the Einsteins field equations for static spherically symmetric manifolds are important in the description of such anisotropic relativistic spheres. The interior spacetime matches smoothly with the exterior Schwarzschild model [1]. Many exact solutions to the field equations have been generated by different approaches with generalized forms for one of the gravitational potentials that does have an equation of state (EoS) (linear [2][3][4][5][6][7], quadratic [8][9][10], polytropic [11][12][13][14][15][16], Van der Waals [17], etc.) and without [18][19][20][21][22][23][24][25] a particular barotropic EoS relating the pressure to the energy density. However, among large number of such work reported over the years, relatively few of these solutions correspond to non-singular metric functions with physically acceptable energy momentum tensor.
In this work we generate a general frame work describing static spherically symmetric anisotropic matter distribution in terms of an equation of state p r = p r (ρ) and gravitational potential as generating function. We treat Van der Waals type a e-mail: chalmusragel@esn.ac.lk b e-mail: thirukkanesh@esn.ac.lk of EoS and a physically reasonable form of gravitational potential as a particular example to model mixed fluid astrophysical objects that is of physical significance. There are many phenomena in stellar objects where first-order phase transitions occur (e.g., pion condensation, kaon condensation, deconfinement and chiral phase transitions) [26][27][28], and the Van der Waals fluid is a good model to illustrate principal features of such scenario [29,30]. At low-temperature high-density limits, as the densities grow the separation distance between particles decreases and quantum degeneracy pressure due to the exclusion principle supports against further collapse and the average particle energy increases. This leads to the deconfinement of the nucleons and mixed phases are expected with various first-order phase transitions, such as liquid-gas transition in the low densities [31][32][33], meson condensations and hadron-quark deconfinement transitions at high-density region [34][35][36][37].
Meson condensations (pions and kaons) have been extensively investigated as new states of high-density hadronic matter, which may be realized in hybrid stars. Meson condensed states are strongly interacting systems of mesons and baryons whose dynamics is controlled by the underlying chiral symmetry, and quark confinement/deconfinement play an important role in connecting hadron dynamics and quarkgluon dynamics, leading to a unified description. The breaking of chiral symmetry triggers the onset of kaon condensation in both hadronic matter phase and the color-flavor locked (CFL) phase [38]. The authors have shown that a Van der Waals type EoS with appropriate choice of parameters can be used to model a neutral strange quark star with ultra-dense matter in spherically symmetric mixed phases in color singlet domain [17].
On the other hand, in the low-density regime (ρ s ≈ 0.16 fm −3 ) at finite temperature (T < 20 MeV), there appears a "liquid-gas" mixed phase equilibrium [27,31], where the coexisting nuclear phase (liquid) and nuclear/ electron phase (gas) are non-congruent and are separated in space forming pasta structures. Such frustrated system of "pasta phase" appears due to competition between the Coulomb interaction and the strong force, and is constituted by different geometrical configurations [39][40][41]. Nuclear pasta forms in core-collapse supernova events [32,33,42] and in the inner crust of neutron stars [43], and observation of the signals of the mixed phase has been suggested in the spectra of the gravitational waves [44]. Pressure anisotropy [45][46][47] is expected to play an important role in such mixed fluid and we show Van der Waals EOS could represent a low-density scenario of a star.

Anisotropic model with EoS
The interior of static spherically symmetric star can be described by the metric For anisotropic imperfect neutral matter distribution the energy momentum tensor can be taken as In the above equation the energy densityρ,, radial pressure p r and tangential pressure p t are measured relative to the comoving fluid 4-velocity u i = e −ν δ i 0 . For the metric (1) and energy momentum tensor (2) the Einstein field equations becomes where primes denote differentiation with respect to r . To obtain the Einstein system of Eqs. (3)-(5), we assumed 8π G = 1 and the speed of light c = 1. The system of Eqs.
(3)-(5) describes the behaviour of the gravitational field inside an anisotropic star. The mass contained within the sphere of radius r is given by Using the transformation suggested by Durgapal and Bannerji [48] x = r 2 , the system (3)-(5) becomes 4x Zÿ y where A is an arbitrary constant and dots denote differentiation with respect to the variable x. Under this transformation (7) the mass function (6) becomes A physically realistic fluid matter distribution is expected to satisfy the barotropic equation of state With the inclusion of (12), the general solution to the system (8)-(10) can be written as in terms of the gravitational potential Z , where the quantity Δ = p t − p r is the measure of anisotropy, and d is a constant of integration. Therefore the line element (1) takes the form Hence, any solution describing static spherically symmetric anisotropic matter distribution with any form of equation of state can be easily determined by the generating function Z (r ) if p r = p r (ρ) is specified. In particular, for a linear EoS where α and β are real constants the metric (20) reduces to The metric (22) with Tolmann IV type potential for Z has been shown to well model a strange star candidate [6].

Particular model with Van der Waal's EoS
In principle we can express any model explicitly from the system (13)-(18) that describe the behaviour of an anisotropic star if the form of the equation of state p r = p r (ρ) and the gravitational potential Z are specified. Consequently, in this paper, we consider the modified Van der Waals equation of state where α, β and γ are real constants together with the form for the gravitational potential where a is a positive real constant. The gravitational potential Z is regular at the origin and well behaved in the stellar interior for a wide range of values for the parameters a. This form of gravitational potential were previously used to study an isotropic perfect fluid sphere [48] and shown to satisfy all the physical requirements for realistic neutron star [49]. Later this model were generalized by Thirukkanesh and Maharaj [19] and Thirukkanesh and Ragel [50]. Therefore the form chosen in (24) is physically viable. Substituting (23) and (24) into (19) we obtain where Hence an exact model for the system (13)-(18) is as follows where y is given by (25). The solution (26)

Physical analysis
In the analysis of physical viability, the following conditions have been generally recognized to be crucial for anisotropic fluid spheres [51,52]: (i) the density ρ and the radial pressure p r should be positive inside the star; (ii) the density ρ, the radial pressure p r and the tangential pressure p t should be monotonically decreasing from center to surface; (iii) the radial pressure p r must vanish but the tangential pressure p t need not necessarily vanish at the boundary of the sphere. However, the radial pressure should be equal to the tangential pressure at center of the fluid sphere; (iv) inside the static configuration the speed of sound should be less than the speed of light, i.e., 0 ≤ V 2 Sr = dp r dρ ≤ 1 and 0 ≤ V 2 St = dp t dρ ≤ 1; (v) the energy-momentum tensor has to obey the conditions ρ − p r − 2 p t ≥ 0 and ρ + p r + 2 p t ≥ 0; (vi) At the boundary of the star r = R, the interior metric matches smoothly with the exterior Schwarzschild where M is the total mass of the sphere.
Now we show that the model generated satisfy the above physical conditions and hence, physically represent a realistic star. From condition (iii), at the boundary of the star r = R, we have p r (R) = 0 which yield For suitable choice of parameters α, β and γ we can get the positive radius R. The matching condition (vi) imply where M = m(R) = 3a R 3 4(1+a R 2 ) . The condition (33) does not impose any restrictions on the parameters. However the condition (34) impose restriction on the parameters A as Due to complexity of the solution, we show graphically that the matter variable are well behaved throughout the interior of the star and the plotted radial dependence of physical quantities of the model illustrates that the model satisfy the conditions (i)-(v). The graphs have been plotted for a particular choice of parameter values a = 0.001, α = 45, β = − 0.1655 and γ = 1 with stellar boundary R = 11.5 km. The choice of parameters have been such that the model is physically viable satisfying the conditions (i)-(v). Figures 1 and 2 show that the energy density and the radial pressure decreases monotonically with increasing radius, and the density is non-zero at the surface and the radial pressure vanishes, which is the requirement for realistic star. However, as seen in Fig. 3, the tangential pressure is non-zero at the surface, which is physically possible [51], and such a scenario of anisotropic fluid has been explained by assumption of particles moving on circular orbits [53,54] and the tangential pressure of a surface layer is related to surface tension [55]. As demonstrated in Fig. 4, the anisotropy Δ is repulsive  and increases monotonically towards the surface boundary remaining finite and continuous in the interior. Moreover, Figs. 5 and 6 illustrate that inside the static configuration the square of speed of sound (both radial and tangential) remain between 0 and 1, satisfying another physical requirement. Figures 7 and 8 show that the energy-momentum tensor of the model is physically admissible, satisfying the condition (v). Therefore the model, for the chosen parameters satisfy all major physical conditions (i)-(vi) listed above of a realistic star. Moreover, Figure 9 illustrates that the Page 5 of 6 306   St − V 2 Sr ≤ 0. For the model generated by these parameter values, we calculate the mass of the stellar object to be M = 0.6818M , central density ρ(0) = 2.4123 × 10 14 g cm −3 and surface density ρ(R) = 1.9646 × 10 14 g cm −3 , which shows that the stellar object has a density less than the nuclear saturation density 2.8×10 14 g cm −3 , with model parameters that could   represent the "liquid-gas" mixed phase equilibrium in the low-density regime.

Conclusion
In this paper, a general frame work is generated to describe static spherically symmetric anisotropic matter distribution with an equation of state p r = p r (ρ) and gravitational poten-