Physical Attributes of Anisotropic Compact Stars in $f(R,G)$ Gravity

TL;DR

This paper investigates the physical properties and stability of anisotropic compact stars within $f(R,G)$ gravity, using specific models and graphical analysis to understand their behavior under modified gravity theories.

Contribution

It introduces a combined $f(R,G)$ gravity model and analyzes anisotropic star properties, energy conditions, and stability through 3D graphical methods.

Findings

Stars behave normally for positive $f(G)$ parameter $n$

Energy conditions are satisfied in the models

Stars remain stable under the considered conditions

Abstract

Modified gravity is one of the potential candidates to explain the accelerated expansion of the universe. Current study highlights the materialization of anisotropic compact stars in the context of $f(R,G)$ theory of gravity. In particular, to gain insight in the physical behavior of three stars namely, Her $X1$, SAX J $1808$-$3658$ and 4U $1820$-$30$, energy density, and radial and tangential pressures are calculated. The $f(R,G)$ gravity model is split into a Starobinsky like $f(R)$ model and a power law $f(G)$ model. The main feature of the work is a $3$-dimensional graphical analysis in which, anisotropic measurements, energy conditions and stability attributes of these stars are discussed. It is shown that all three stars behave as usual for positive values of $f(G)$ model parameter $n$.

Equations

$$G\equiv R^2-4R_{\mu \nu }{R}^{\mu \nu }+R_{\mu \nu \theta \varphi }{R}^{\mu \nu \theta \varphi },$$

$$S=\frac{1}{2K}\int d^{4}x\sqrt{-g}f(R,G)+S_M(g^{\mu \nu },\psi ).$$

$$R_{\mu \nu }-\frac{1}{2}{g}_{\mu \nu }R$$

$$V\equiv f_{R}R+f_{G}G-f(R,G),$$

$$d{s}^2=e^{a(r)}d{t}^2-e^{b(r)}d{r}^2-r^2(d{\theta }^2+{\sin }^2\theta d{\varphi }^2).$$

$$T_{\alpha \beta }^m=(\rho +p_t)u_{\alpha }{u}_{\beta }-p_t{g}_{\alpha \beta }+(p_r-p_t)v_{\alpha }{v}_{\beta },$$

$$\rho$$

$$p_r$$

$$p_t$$

$$\rho$$

$$8A^2{B}^2{r}^8+2n(1-6A^2{r}^4+2A(r^2+B{r}^4){)}^2+e^{2A{r}^2}(-1+3B{r}^2-2A^4{r}^8+A^3{r}^6\times$$

$$(3+4B{r}^2)-2A^2{r}^4(-2+B{r}^2+B^2{r}^4)-A{r}^2(8+6B{r}^2+B^2{r}^4)+2n(1-A^2{r}^4+$$

$$A(r^2+B{r}^4){)}^2+e^{A{r}^2}(2-6B{r}^2+18A^4{r}^8+3A{r}^2(3+B{r}^2{)}^2-3A^3{r}^6(5+8B{r}^2)+$$

$$6A^2{r}^4(-7+B^2{r}^4)-4n(1+6A^4{r}^8+3A(r^2+B{r}^4)+A^2{r}^4(-5+4B{r}^2+2B^2{r}^4)-$$

$$8A^3(r^6+B{r}^8))))+16e^{-2A{r}^2}(-3+3e^{A{r}^2}-B^2{r}^4+2A^3{r}^6(2+B{r}^2)-$$

$$A^2{r}^4(4+11B{r}^2+2B^2{r}^4)+A{r}^2(-3+4B{r}^2+5B^2{r}^4))\lambda -$$

$$16e^{-2A{r}^2}(2+A{r}^2)(-1+e^{A{r}^2}+B^2{r}^4(2+B{r}^2)-A(r^2+4B{r}^4+B^2{r}^6))\lambda +$$

$$2e^{-2A{r}^2}(-1+e^{A{r}^2}-3B{r}^2-B^2{r}^4+A{r}^2(2+B{r}^2))\times$$

$$(e^{A{r}^2}(r^2-2\lambda )+2(1+3B{r}^2+B^2{r}^4-A{r}^2(2+B{r}^2))\lambda )+2B{e}^{-2A{r}^2}{r}^2$$

$$\displaystyle(3-Ar^{2}+2Br^{2})\big{(}e^{Ar^{2}}(r^{2}-4\lambda)+4(1+3Br^{2}+B62r^{4}-Ar^{2}(2+Br^{2}))\lambda\big{)}\bigg{]},$$

$$p_r$$

$$p_t$$

$$p_r={\text{w}}_r\rho ,p_t={\text{w}}_t\rho .$$

$$0<{\omega }_r<1,0<{\omega }_t<1.$$

$$\frac{d\rho }{dr}=0,\frac{d{p}_r}{dr}=0,\frac{d^2\rho }{d{r}^2}<0,\frac{d^2{p}_r}{d{r}^2}<0.$$

$$d{s}^2=(1-\frac{2M}{r})d{t}^2-(1-\frac{2M}{r}{)}^{-1}d{r}^2-r^2(d{\theta }^2+{\sin }^{2}d{\phi }^2).$$

$$g_{tt}^{-}=g_{tt}^{+},g_{rr}^{-}=g_{rr}^{+},\frac{\partial g_{tt}^{-}}{\partial r}=\frac{\partial g_{tt}^{+}}{\partial r},$$

$$A=-\frac{1}{R^2}\ln (1-\frac{2M}{R}),$$

$$B=\frac{M}{R^3}(1-\frac{2M}{R}{)}^{-1},$$

$$C=\ln (1-\frac{2M}{R})-\frac{M}{R}(1-\frac{2M}{R}{)}^{-1}.$$

$$NEC:\rho +{\text{p}}_r\ge 0,\rho +{\text{p}}_t\ge 0,$$

$$WEC:\rho \ge 0,\rho +{\text{p}}_r\ge 0,\rho +{\text{p}}_t\ge 0,$$

$$SEC:\rho +{\text{p}}_r\ge 0,\rho +{\text{p}}_t\ge 0,\rho +{\text{p}}_r+2{\text{p}}_t\ge 0,$$

$$DEC:\rho >\mid {\text{p}}_r\mid ,\rho >\mid {\text{p}}_t\mid .$$

$$\frac{d{p}_r}{dr}+\frac{\nu \stackrel{ˊ}{}(\rho +p_r)}{2}+\frac{2(p_r-p_t)}{r}=0.$$

$$F_g+F_h+F_a=0,$$

$$F_g=-Br(\rho +p_r),F_h=-\frac{d{p}_r}{dr},F_a=\frac{2(p_t-p_r)}{r}.$$

Full text

Keywords: $f(R,G)$ Gravity, Spherically Symmetric Spacetime, Compact Stars.

PACS: 04.50.Kd, 98.80.-k.

1 Introduction

The spatial behavior of constituents of the universe is a complex phenomena. Visual observation of the remote universe enlightens us about its continuous expansion, and accentuates to find alternative rationale which may be helpful to explain the phenomena of expansion of the universe [1]. So, these alternative models are known currently as modified theories of general relativity (GR), some of which are $f(R)$, $f(G)$, $f(R,T)$, $f(G,T)$ and $f(R,G)$ theories of gravity, where $R$, $G$, and $T$ denotes Ricci scalar, the Gauss-Bonnet invariant and the trace of the energy momentum tensor respectively. The theory of GR explains the cosmological phenomena in weak field regimes while some modifications are required to address the strong fields in the scenario of continuous expansion of universe. The same thought was provoked and one of the modified theories, $f(R)$ gravity, was proposed in $1970$ by Buchdahl [2]. Hydrostatic equilibrium and stellar structure in $f(R)$ gravity have been investigated by considering the Lan e-Emden equation [3]. Some finite time future singularities in modified gravity were described with their solution using the addition of higher derivative gravitational invariants [4]. Harko [5] presented modified $f(R,T)$ theory in $2011$ by involving both matter and curvature terms. Furthermore, the next modification to the GR was Gauss-Bonnet gravity [6], also known as Einstein-Gauss-Bonnet gravity, which includes the Gauss-Bonnet term.

[TABLE]

where $R_{\mu\nu}$ and $R_{\mu\nu\theta\phi}$ are Ricci tensor and the Riemann tensor respectively. $f(G)$ gravity [7] is the further generalization of Gauss-Bonnet gravity. It has been shown that this kind of generalization may naturally approach an operative cosmological constant, quintessence and phantom cosmic acceleration [8]. Recently, another alternative theory with the title of $f(G,T)$ gravity has been proposed by Sharif and Ikram [9]. They studied the energy conditions for the well-known Friedmann -Robertson- Walker space time. So far some interesting work has been done in this modified theory [10]-[13]. Similarly, $f(\mathcal{T})$ gravity is also a modified theory, which has attracted the attention of theoretical physicists gain insight of acceleration and regular thermal expansion of late-time universe by assuming an alternative of cosmological constant, where $f(\mathcal{T})$ is a function of the torsion scalar [14]. An endeavor was made by combining $R$ and $G$ in a bivariate function $f(R,G)$ [15]-[22]. The function $f(R,G)$ provides a foundation for the double inflationary scenario [23]. It has been shown that $f(R,G)$ theory is consistent with the observational data [24, 25]. Besides its stability, this theory is perfectly suitable to explain the accelerating waves of the celestial bodies as well as the phantom divide line crossing and transition from acceleration to deceleration phases. An important feature of $f(R,G)$ gravity is that it reduces the risk of ghost contributions and gravitational action is regularized because of the term $G$ [26, 27]. Thus, modified theories of gravity seem interesting in explaining the universe in different cosmological contexts [28].

Study of compact stars has always been a strong topic of research [29]-[40]. Different properties like mass, radius and moment of inertia of neutron stars are studied and comparison has been developed with GR and alternative theories of gravity [41]. Some investigations of the structure of slowly rotating neutron stars in $R^{2}$ gravity have been done using two different hadronic parameters and a strange matter equation of state (EOS) parameter [42]. The mass-ratio for compact objects in the presence of cosmological constant has been derived [43]. In this article, we are interested in studying the compact stars and in examining their physical behavior using standard models of spherically symmetric space time by finding exact solution of field equations. We study the stability of anisotropic compact stars using generalized modified Gauss-Bonnet $f(R,G)$ gravity. The arrangement of this paper is as follows: In the second section, we find expressions for energy density, radial pressure and tangential pressure for anisotropic matter source using spherically symmetric metric in $f(R,G)$ gravity. Section 3 comprises the analysis of behavior of EOS parameters, anisotropy measure, the matching conditions, the energy bounds, Tolman Oppenheimer Volkoff (TOV) equation and stability analysis of three compact stars namely, Her $X1$, SAX J $1808$-$3658$ and 4U $1820$-$30$. The last section ends with the conclusion and summery of the paper.

2 Anisotropic Matter Configuration in $f(R,G)$ Gravity

We start with the action of modified Gauss-Bonnet gravity [44]

[TABLE]

Varying action (2) with respect to metric tensor yields the following modified field equations [25]

[TABLE]

where $f_{R}$ and $f_{G}$ are partial derivatives with respect to $R$ and $G$ respectively,

[TABLE]

and $T_{\mu\nu}^{(matt)}$ describes the ordinary matter. The most general spherically symmetric space time is [45]

[TABLE]

The energy-momentum tensor in case of anisotropic fluid is given by

[TABLE]

where $u_{\alpha}=e^{a/2}\delta_{\alpha}^{0}$, $v_{\alpha}=e^{b/2}\delta_{\alpha}^{1}$ are four velocities. Radial pressure and tangential pressures are $p_{r}$ and $p_{t}$ respectively, and the energy density is denoted by $\rho$. Using equations (5) and (6) in field equations (3), and after some manipulations we obtain,

[TABLE]

[TABLE]

[TABLE]

The set of three equations (7)-(9) involves five unknown functions $\rho,p_{r},p_{t},a,b$. Moreover, the equations are too much complicated and highly non-linear due to the involved bivariate function $f(R,G)$. So following Krori and Barua [45], we choose $a(r)=Br^{2}+C$, $b(r)=Ar^{2}$, $A$, $B$, and $C$ are constants. These constants will be determined by using some physical assumptions. For the present analysis, we propose $f(R,G)=f_{1}(R)+f_{2}(G)$. We further consider the Starobinsky like model $f_{1}(R)=R+\lambda R^{2}$, where $\lambda$ is an arbitrary constant and we take $f_{2}(G)=G^{n}$, $n\neq 0$. Also $f_{1R}=\frac{df_{1}}{dR}$, $f_{2G}=\frac{df_{2}}{dG}$ and a prime denotes the derivatives with respect to the radial coordinate. Using these assumptions, Eqs. (7)-(9) take the form

[TABLE]

[TABLE]

[TABLE]

[TABLE]

3 Physical Analysis and Graphical Representation

This section presents the physical properties of the solutions regarding EOS, anisotropic behavior, energy conditions, TOV and matching conditions along with stability analysis of three different compact stars, Her $X1$, SAX J $1808$-$3658$ and $4$U $1820$-$30$. Many EOS parameters have been considered in literature in different cosmological contexts. Staykov et al. [42] used two different hadronic parameters and a strange matter EOS parameter to study the structure of rotating neutron stars. Quadratic EOS parameters have been used to study the properties of compact stars [46, 47]. For the sake of simplicity, in this work we assume linear EOS [48]:

[TABLE]

It is important to note here that the EOS parameters are dependent on radius rather than a constant quantity as in the ordinary matter distribution. The non constant behavior is due to the usual matter and exotic matter contributions. Figs. $(1)$ and $(2)$ show the variation of EOS parameters $\omega_{r}$ and $\omega_{t}$ respectively. It is obvious from the figures that effective EOS in our model are the same as in normal matter distribution [49], i.e.

[TABLE]

3.1 Anisotropy Measure of the Compact Stars

The variation of energy density $\rho$, radial pressure $p_{r}$ and tangential pressure $p_{t}$ can be observed in Figs. $(3)$-$(5)$. The behavior of $\frac{d\rho}{dr}$ and $\frac{dp_{r}}{dr}$ in Figs. $(6)$ and $(7)$ reveals that $\frac{d\rho}{dr}<0$ and $\frac{dp_{r}}{dr}<0$. It shows that with the increase in radius of the compact star, both energy density and radial pressure decrease. We analyze the variation of $\frac{d\rho}{dr}$ and $\frac{dp_{r}}{dr}$ at the center of the compact star $r=0$ and found that

[TABLE]

Equation (15) gives the maximum value of $\rho$ and $p_{r}$ at center $r=0$. The anisotropy measurement $\triangle=\frac{2}{r}(p_{t}-p_{r})$ has been shown graphically in Fig. $(14)$. The anisotropy measurement is directed outward when $p_{t}>p_{r}$ which results $\triangle>0$ and directed inward when $p_{t}<p_{r}$ which results $\triangle<0$. It is depicted in Fig. $(14)$ that for the large values of $r$, $\triangle>0$ for all stars implying that the anisotropic force permits the construction of great massive configurations. It is worthy to mention that anisotropy measurement vanishes at the center of the star.

3.2 The Matching Conditions with Schwarzchild Exterior metric

The interior metric of the boundary surface remains same in the interior and exterior geometry of the compact star. It justifies the continuity of the metric components for the boundary surface of the star. Many choices for the matching conditions are possible and one can consider the vacuum outside a general spherically symmetric space time with suitable boundary conditions [50]. For the present analysis, we consider Schwarzchild solution for describing the exterior geometry. Many authors have considered Schwarzchild solution for this purpose giving some interesting results [51]-[53]. Therefore, the exterior metric given by Schwarzchild is

[TABLE]

The intrinsic metric (5) for the smooth match at the boundary surface $r=R$ with Schwarzchild exterior metric produces,

[TABLE]

where interior solutions and exterior solutions are represented by (-) and (+). By the matching of interior and exterior metrics, we obtain

[TABLE]

Using the approximate values of $M$ and $R$ for the compact stars under observation, the constants $A$ and $B$ are given in the following table [49]

[TABLE]

3.3 The Energy Bounds

The energy bounds have gained much importance in the discussion of some important issues in cosmology. In fact, one can investigate the validity of second law of black hole thermodynamics and Hawking-Penrose singularity theorems using energy conditions [54]. Many interesting results have been reported in cosmology using the energy bounds [55]-[60]. These energy conditions are defined as

[TABLE]

where the null energy conditions, weak energy conditions, strong energy conditions and dominant energy conditions are denoted by $NEC$, $WEC$, $SEC$ and $DEC$ respectively. In Fig. $(8)$, it is evident that all energy conditions are satisfied for Her $X1$. The energy conditions are also satisfied for other two stars but the graphs are not shown here.

3.4 Implementation of Tolman-Oppenheimer-Volkoff Equation

The TOV equation is expressed in the following generalized form

[TABLE]

It follows from equation (21),

[TABLE]

where

[TABLE]

Here $F_{g}$, $F_{h}$ and $F_{a}$ are gravitating force, hydrostatic force and anisotropic pressure force for compact stars. Using the values of $\rho$, $p_{r}$ and $p_{t}$ from equations (10)-(12), for compact star $4U$, Fig.$(9)$ shows the behavior of these forces. Fig. $(10)$ depicts that the TOV equation is satisfied for compact star $4$U. The TOV equation is satisfied for the other two stars as well and graphical analysis is not presented here.

3.5 Stability Analysis

Now we find the radial sound speed $\upsilon_{sr}$ and transverse sound speed $\upsilon_{st}$ to determine the stability of our model where

[TABLE]

For a stable model following conditions must hold.

[TABLE]

The conditions (25) are exhibited graphically in Figs. ($11$) and ($12$) and show that these conditions are true for considered compact stars. In $1992$, Herrera gave the important concept that for a potentially stable region $\upsilon_{sr}$ is greater than $\upsilon_{st}$ [61]. Therefore, behavior of $\upsilon_{sr}^{2}-\upsilon_{st}^{2}$ is shown in Fig.($13$). It can be observed that $\mid\upsilon_{sr}^{2}-\upsilon_{st}^{2}\mid<1$. So, proposed models of compact stars are stable (Fig. $13$).

4 Final Remarks

The current study deals with the physical attributes of the compact stars in the scenario of $f(R,G)$ gravity. In this regard, interior solutions are figured out for three compact stars namely, Her $X1$, SAX J $1808$-$3658$ and 4U $1820$-$30$, with the assumption that they have anisotropic internal structure. The analytic solutions of the interior metric in $f(R,G)$ are smoothly matched with Schwarzchild exterior metric. This matching has provided the values of unknown constants $A$, $B$, and $C$ in the form of observed values of radii and masses of the model compact stars. The nature of the stars has been discussed using values of these constants. The analysis of the physical attributes leads one to conclude the following results as regards the anisotropic compact stars in $f(R,G)$ gravity:

• •

The EOS parameters for compact stars are true as in the case of ordinary matter distribution in $f(R,G)$ gravity which shows that the compact stars are composed of ordinary matter. The matter density, and radial and tangential pressures get the maximum values at the center of the star and they are also decreasing functions. It confirms the fact that the matter components of compact stars are positive and remain finite in the interior of stars. Thus, all three compact stars in this present study are singularity free.

• •

It is also noticed that anisotropic force is directed outward for $p_{t}>p_{r}$ which means $\triangle>0$, while on the other hand the anisotropic force is directed inward for $p_{t}<p_{r}$ which shows that $\triangle<0$ . The graphical description of $\triangle>0$ for three different compact stars is revealed in Fig. ($14$).

• •

The energy conditions and the TOV are analyzed for all three compact stars under consideration. The behavior of energy conditions and TOV are shown in Fig.($8$) and Fig.($10$) for stars Her $X1$ and 4U respectively. Although these are verified for remaining two stars as well. Furthermore, the inequality $\mid\upsilon_{sr}^{2}-\upsilon_{st}^{2}\mid<1$ is satisfied as shown in Fig. ($13$) for the considered compact stars. So these models of stars are stable.

The important feature of the present study is the $3$-dimensional analysis and all the results show that the compact stars behave as usual for the $f(G)$ model parameter $0<n<200$. Thus it is conjectured that the stars behave in the same way for a polynomial form of $f(G)$, $n\geq 200$. Also the behavior of the stars can be checked for negative values of $n$ as well and considering some other forms of $f(R)$ gravity model. It is worth mentioning that the findings of the current study are in conformity with the results in [62] for $n=2$ and $f_{1}(R)=0$.