Stability and improved physical characteristics of relativistic compact objects arising from the quadratic term in pr=αρ2+βρ-γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_r = \alpha \rho ^2 + \beta \rho - \gamma $$\end{document}

We investigate the stability and enhancement of the physical characteristics of compact, relativistic objects which follow a quadratic equation of state. To achieve this, we make use of the Vaidya–Tikekar metric potential. This gravitational potential has been shown to be suitable for describing superdense stellar objects. Pressure anisotropy is also a key feature of our model and is shown to play an important role in maintaining stability. Our results show that the combination of the Vaidya–Tikekar gravitational potential used together with the quadratic equation of state provide models which are favourable. In comparison with other equations of state, we have shown that the quadratic equation of state mimics the colour-flavour-locked equation of state more closely than the linear equation of state.


Introduction
Compact objects such as neutron stars and the more specialised strange quark stars are ideal candidates for applying Einstein's theory of general relativity and for studying matter under extreme conditions. Such compact systems also draw attention to the characteristics and interactions of quarks and how this manifests itself within the physical properties of these stars. The pioneering works of Witten [1] and more recently that of Weber [2] include in-depth studies into the stability and properties of quark matter and the possibility for the existence of strange-matter stars. From the point of view of the microphysics of compact matter, the basic features of the theory of quantum chromodynamics (QCD) have a e-mail: Thirukkanesh@esn.ac.lk b e-mail: bogadi.robert@gmail.com (corresponding author) c e-mail: megandhreng@dut.ac.za d e-mail: moyos@dut.ac.za been used to assist in formulating equations of state (EoS) [3]. In the context of massive compact objects such as quarkstars, the large masses, high energy densities and pressures involved require relativistic treatments using Einstein's theory of general relativity (GR) or the more recent higher-order gravity theories such as Einstein-Gauss-Bonnet gravity [4].
Equations of state remain a key aspect in the study of compact objects. In work involving GR, equations of state have been used in setting gravitational potentials and also ensuring that physical viability is maintained. In particular, the linear equation of state, p = α(ρ − ρ s ) has been popular and successful in numerous studies [5,6]. The colour-flavour-locked equation of state has also recently been actively utilised [7], allowing for a more in-depth description of the microphysics of highly compact material which has surpassed the nuclear saturation density. If strange stars have a non-homogeneous, shell-type internal structure, as proposed for neutron stars, then different regions within the quark matter might invite the use of say region-specific equations of state as opposed to using a single EoS for the entire description of the star. This has already been considered in so-called hybrid stars [8] in which the high pressures and energy densities of the quark core might be underestimated by the simplistic linear equation of state. Thus for neutron stars of the order of 2M or greater, a quadratic equation of state offers the possibility for augmenting the pressures and densities within the core of these more massive stars. The quadratic EoS has been investigated and found to offer exact solutions to the Einstein-Maxwell field equations which are physically acceptable [9][10][11]. Physical characteristics within a star could include additional processes such as the production of hyperons and generation of condensates, and an environment for the production of hyperons and generation of condensates would be favoured by higher core densities and pressures. A recent method in which a quadratic EoS is employed for taking into account the distinguishing features of the core of massive relativistic stars is the core-envelope method [11,12]. Although we do not employ this method in our study, it nevertheless promotes the use of the quadratic equation of state for superdense, massive objects.
In considering highly compact matter from the viewpoint of GR, it is necessary to employ gravitational potentials suited to superdense matter. The Vaidya-Tikekar potential which is known to generate superdense stellar models, is favourable in this case [13]. Within the ansatz is the spheroidal parameter K which was originally set at K = −2 by Vaidya and Tikekar. This however can be adjusted to improve the computed physical characteristics for the star, should certain quantities such as pressure anisotropy and sound speed show unreasonable behaviour or if the system shows marked instability. Other metric potentials such as that of Finch and Skea [14] have also been employed in compact systems governed by a quadratic equation of state [15] although the extension back to the linear regime might not be possible. Such a restriction does not arise in our study.
Pressure anisotropy is also a key feature of our model. Since the pioneering work of Bowers and Liang [16], the benefits of considering systems which incorporate anisotropy have been well-noted [17,18]. In this paper, the method used and solutions obtained are similar to those of Maharaj and Takisa [20] with the emphasis being on application to stars of known masses and well-predicted radii. This is achieved via the Vaidya-Tikekar potential. Realistic estimates of energy density, anisotropic pressures and stability parameters are sought through variation of the equation of state parameters in the quadratic EoS of the form, p = αρ 2 + βρ − γ .

The field equations
We assume the spacetime manifold to be static and spherically symmetric. This assumption is consistent with the development of models used to study the physical behaviour of relativistic astrophysical objects such as neutron stars and other similarly compact objects. The interior geometry of a spherically symmetric static star is described by the line element (1) in Schwarzschild coordinates (x a ) = (t, r, θ, φ). We take the energy momentum tensor for an anisotropic neutral imperfect fluid sphere to be of the form where ρ is the energy density, p r is the radial pressure and p t is the tangential pressure. These quantities are measured relative to the comoving fluid 4-velocity u i = e −ν δ i 0 . For the line element (1) and matter distribution (2) the Einstein field equations can be expressed as where primes denote differentiation with respect to r and our choice of units are such that 8π G/c 4 = 1. The system of equations (3)-(5) governs the behaviour of the gravitational field for an anisotropic imperfect fluid. The mass contained within a radius r of the sphere is then given by The field equations can be cast in a different but equivalent form by introducing the transformation where C is a positive real constant. This transformation was first suggested by Durgapal and Bannerji [21]. Under this transformation, the system (3)-(5) is described by 4x Zÿ y where dots denote differentiation with respect to the variable x. The mass function (6) as computed under the above transformation is given by For a physically realistic relativistic star we expect that the matter distribution should satisfy a barotropic equation of state p r = p r (ρ) and in this study we proceed with a quadratic equation of state of the form where α, β and γ are real constants. Then it is possible to write the system (8)- (10) in the simpler form where the quantity Δ = p t − p r is a measure of the pressure anisotropy. A similar approach has been followed based on a linear equation of state [22]. The above system of equations governs the behaviour of the gravitational field for an imperfect fluid source.

Exact models
In the system (13)-(17), there are six independent variables (ρ, p r , p t , Δ, y, Z ) and only five independent equations. This suggests that it is possible to specify one of the quantities involved in the integration process. Equation (17) is the master equation in the integration process. In this treatment we specify the gravitational potential Z so that it is possible to integrate (17). The explicit solution of the field equations (13)-(17) then follows. We choose where K is a real constant. The gravitational potential Z was originally used by Vaidya and Tikekar [13] and more recently by Bhar [27] to study superdense stars, and the form of (18) has been found to be physically reasonable. Substituting (18) into (17) we obtaiṅ On integrating (19) we then obtain The solution (21) -(25) may now be applied to modelling an anisotropic, relativistic star according to a quadratic equation of state. Applying the solution to (11), the mass function takes the form The surface redshift is given by The proper radius is given by

Physical constraints
In assessing the physical viability of solutions to Einstein's field equations, the following conditions are tested, and implemented in the case of boundary conditions, for compact, anisotropic fluid spheres [17,18]: (i) non-negative values for the energy density ρ and the radial pressure p r inside the star; (ii) monotonically decreasing profiles of the energy density ρ, the radial pressure p r and the tangential pressure p t from centre to surface; (iii) vanishing radial pressure p r at the surface boundary; (iv) subluminal sound speeds within the static configuration, i.e., 0 ≤ V 2 rs = dp r dρ ≤ 1 and 0 ≤ V 2 ts = dp t dρ ≤ 1; (v) stability condition, i.e., −1 ≤ V 2 ts − V 2 rs ≤ 0; (vi) constraints for the energy-momentum tensor: ρ − p r − 2 p t ≥ 0 and ρ + p r + 2 p t ≥ 0; (vii) smooth matching of the interior metric with the exterior Schwarzschild metric at the boundary of the star r = R, where M is the total mass of the sphere.

Physical application
Three well-studied systems, Cen X-3, Vela X-1 and PSR J1614-2230 [19], were selected with evenly spaced masses above 1.4M as shown in Table 1. Stars of lower mass were avoided as the lower central densities and pressures expected for such systems would be less suited to our gravitational model. A value of −2 was initially investigated in our study and significantly lower densities and pressures were obtained. In addition, the anisotropy profiles for K = −2 displayed some unphysical characteristics, especially for lower mass stars. The anisotropy profiles for our systems utilising K = −5 are more typical in both magnitude and curve progression as shown in Fig. 3.

Discussion
We now provide further discussion of the trends and physical viability of our models. In Fig. 1 we observe that the density profiles are smooth, monotonically decreasing functions of the radial coordinate. As the quadratic parameter α increases, the magnitude of the density increases at each interior point of the compact object. Figure 2 (left panel) displays the radial pressure as a function of the radial coordinate. As expected, the radial pressure decreases monotonically from some finite value at the centre and vanishes at the boundary of the star. It is interesting to note that an increase in α supports configurations with larger radii. This confirms the trend in the density profiles -larger values of α lead to higher densities. We observe that the tangential pressure (Fig. 2, right panel) is well-behaved. Higher contributions from the quadratic term lead to higher tangential pressures. The anisotropy parameter is plotted in Fig. 3. For each interior point of the stellar configuration, Δ > 0. This means that tangential pressure dominates the radial pressure everywhere inside the star. A positive anisotropy parameter leads to a repulsive contribution from the anisotropic force which may lead to more massive and stable configurations. Figure 4 (left and right panels) indicate that our model obeys the causality requirements. An interesting feature is the minima which appear in the tangential sound speed. This could model an additional anisotropy which distinguishes material closer to the surface from that of the core. The cracking method due to Herrera ascertains that if the tangential speed of pressure wave v t exceeds the radial speed v r , then the model is potentially stable. This requirement can be articulated mathematically as v 2 t − v 2 r < 0. We observe from Fig. 5 that the interior of our stellar model is stable. We further observe that this stability is enhanced with an increase in the quadratic contributions from the EoS. Our model obeys the energy conditions as revealed in Fig. 6. In addition, we have plotted the stability index, Γ in Fig. 7. We note that Γ > 4/3 everywhere inside the stellar fluid.
We have compared the quadratic EoS to the colour-flavour locked (CFL) EoS and the linear EoS in Fig. 8. Similar plots and comparisons have been found in the literature for linearly approximated EoS's [23] and the more recently pursued CFL EoS [24][25][26]. It is interesting to note that the quadratic EoS mimics the CFL EoS more closely than the linear EoS. This is a new and interesting observation. Thirukannesh et al. [7] recently showed that the linear EoS closely resembles the CFL EoS. Our study shows that a better approximation to the CFL EoS is the quadratic EoS. The surface redshift is plotted in Fig. 9. We observe that the redshift increases with an increase in radius and that our results compare well with other studies [27][28][29]. Furthermore, we observe higher surface redshifts for larger values of α. This ties in with our earlier observation that larger values of α lead to higher densities.

Conclusion
We have shown that the quadratic equation of state is wellsuited to the description of highly dense and massive neutron stars, strange-matter stars and possibly the more exotic hybrid stars for which higher core densities and pressures are sought. Our results have shown the effectiveness of the quadratic term in achieving the higher densities and pressures while maintaining stability of the core. Pressure anisotropy is a key feature of our models and the anisotropy parameter displayed radial profiles which were favourable and comparable to other studies. The Vaidya-Tikekar potential used, with the spheroidal parameter (K = −5) used, resulted in the possibility of an additional inhomogeneity as suggested by the minima in the tangential sound speed profiles. According to the sound speed stability criterion, our models appear to be stable except possibly near the surface which could result in a crust which exhibits cracking. The phenomenon of cracking and its relationship to inhomogeneity and anisotropy has been studied by Herrera et al. [18] and the necessity for pressure anisotropy appears to be of central importance in the stability of compact, relativistic objects. The possible instability near the surface as shown in our results may however be reduced by increasing the effect of the quadratic term in the EoS used in this study. Hence the benefit of the quadratic term is highlighted. Comparison of the quadratic equation of state with the linear and the CFL equations of state show that the quadratic term is to some extent a perturbation, featuring more prominently in the core, and that the values of α used serve as adjustments where linearity is likely to fail. The quadratic term also shows enhancements of the surface redshift. Lastly, the adiabatic index shows that the configurations are stable and that increasing α results in a reduction in the rate of decline in stability, albeit the index remaining above Chandrasekhar's limit (Γ ≥ 4/3). The apparent reverse trend in stability as compared to the potential instability analysis as given by the sound speed anisotropy results in Fig. 5, are no doubt a result of the different masses and radii of the candidate stars used in our investigation. In light of this, it is proposed that for a future investigation of the quadratic equation of state with the Vaidya-Tikekar potential, a single candidate star be chosen and studied with respect to variation in the spheroidal parameter.
Acknowledgements RB and MG acknowledge support from the office of the Deputy Vice-Chancellor for Research and Innovation at the Durban University of Technology.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data used to plot the graphs may be reproduced from Table 1, and the relevant equations given in the article.] 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 .