Anisotropic Tolman VII solution by gravitational decoupling

TL;DR

This paper constructs an anisotropic extension of the Tolman VII solution using gravitational decoupling, providing an exact, physically acceptable model for compact objects with differing radial and tangential pressures.

Contribution

It introduces a novel anisotropic Tolman VII solution via gravitational decoupling, expanding the modeling of compact objects with anisotropic pressures.

Findings

Effective density and perfect fluid density are explicitly compared.

Radial and tangential pressures differ, indicating anisotropy.

Anisotropy depends on the coupling constant.

Abstract

Using the gravitational decoupling by the minimal geometric deformation approach, we build an anisotropic version of the well-known Tolman VII solution, determining an exact and physically acceptable interior two-fluid solution that can represent behavior of compact objects. Comparison of the effective density and density of the perfect fluid is demonstrated explicitly. We show that the radial and tangential pressure are different in magnitude giving thus the anisotropy of the modified Tolman VII solution. The dependence of the anisotropy on the coupling constant is also shown.

Equations

$${\hat{G}}_{\mu \nu }=-k^2{\hat{T}}_{\mu \nu };\text{to find the metric}{\hat{g}}_{\mu \nu },$$

$${\tilde{G}}_{\mu \nu }=-k^2{\tilde{T}}_{\mu \nu };\text{to find the metric}{\tilde{g}}_{\mu \nu },$$

$$G_{\mu \nu }=-k^2{T}_{\mu \nu };\text{to find the metric}{g}_{\mu \nu },$$

$$R_{\mu \nu }-\frac{1}{2}R{g}_{\mu \nu }=-k^2{T}_{\mu \nu }^{(\mathrm{tot})},$$

$$T_{\mu \nu }^{(\mathrm{tot})}=T_{\mu \nu }^{(\mathrm{perfectfluid})}+T_{\mu \nu }^{(\mathrm{othersource})}.$$

$$T_{\mu \nu }^{(\mathrm{perfectfluid})}=\left(\rho +p\right)u_{\mu }{u}_{\nu }-p{g}_{\mu \nu },$$

$$T_{\mu \nu }^{(\mathrm{othersource})}=\alpha {\theta }_{\mu \nu }.$$

$$T_{\mu \nu }^{(\mathrm{tot})}=\left(\rho +p\right)u_{\mu }{u}_{\nu }-p{g}_{\mu \nu }+\alpha {\theta }_{\mu \nu }.$$

$${\nabla }_{\nu }{T}^{(\mathrm{tot})\mu \nu }=0.$$

$$d{s}^2=e^{\nu (r)}d{t}^2-e^{\lambda (r)}d{r}^2-r^2(d{\theta }^2+si{n}^2\theta d{\varphi }^2),$$

$$-k^2\left(\rho +\alpha {\theta }_0^0\right)=-\frac{1}{r^2}+e^{-\lambda }\left(\frac{1}{r^2}-\frac{{\lambda }^{\prime }}{r}\right),$$

$$-k^2\left(-p+\alpha {\theta }_1^1\right)=-\frac{1}{r^2}+e^{-\lambda }\left(\frac{1}{r^2}-\frac{{\nu }^{\prime }}{r}\right),$$

$$-k^2\left(-p+\alpha {\theta }_2^2\right)=\frac{1}{4}{e}^{-\lambda }\left(2{\nu }^{\prime \prime }+{\nu }^{\prime 2}-{\lambda }^{\prime }{\nu }^{\prime }+2\frac{{\nu }^{\prime }-{\lambda }^{\prime }}{r}\right).$$

$$-p^{\prime }-\frac{{\nu }^{\prime }}{2}\left(\rho +p\right)+\alpha {\left({\theta }_1^1\right)}^{\prime }-\frac{{\nu }^{\prime }}{2}\alpha \left({\theta }_0^0-{\theta }_1^1\right)-\frac{2}{r}\alpha \left({\theta }_2^2-{\theta }_1^1\right)=0.$$

$${\rho }_{\mathrm{eff}}(\mathrm{radial})=\rho +\alpha {\theta }_0^0,$$

$$p_r=p-\alpha {\theta }_1^1,$$

$$p_t=p-\alpha {\theta }_2^2.$$

$$\pi =p_t-p_r=\alpha ({\theta }_2^2-{\theta }_1^1).$$

$$d{s}^2=e^{ϵ(r)}d{t}^2-e^{\gamma (r)}d{r}^2-r^2(d{\theta }^2+si{n}^2\theta d{\varphi }^2),$$

$$e^{\gamma (r)}={\left(1-\frac{2m(r)}{r}\right)}^{-1},$$

$$e^{ϵ}\mapsto e^{\nu }=e^{ϵ+\alpha g},$$

$$e^{-\gamma }\mapsto e^{-\lambda }=e^{-\gamma }+\alpha f,$$

$$g\mapsto 0,$$

$$f\mapsto f^{\ast }.$$

$$k^2\rho =\frac{1}{r^2}-e^{-\gamma }\left(\frac{1}{r^2}+\frac{{\gamma }^{\prime }}{r}\right),$$

$$k^{2}p=-\frac{1}{r^2}+e^{-\gamma }\left(\frac{1}{r^2}+\frac{{\nu }^{\prime }}{r}\right),$$

$$k^{2}p=\frac{1}{4}{e}^{-\gamma }\left(2{\nu }^{\prime \prime }+{\nu }^{\prime 2}+\frac{2{\nu }^{\prime }}{r}\right)+\frac{1}{4}{\gamma }^{\prime }{e}^{-\gamma }\left({\nu }^{\prime }+\frac{2}{r}\right).$$

$$p^{\prime }=-\frac{{\nu }^{\prime }}{2}(\rho +p).$$

$$k^2{\theta }_0^0=-\frac{f^{\ast }}{r^2}-\frac{f^{{\ast }^{\prime }}}{r},$$

$$k^2{\theta }_1^1=-f^{\ast }\left(\frac{1}{r^2}+\frac{{\nu }^{\prime }}{r}\right),$$

Full text

I Introduction

Einstein’s gravitational field equations are partial nonlinear differential equations – their solution is very difficult with exception of some simplified situations. Immediately after Einstein introduced the general relativity(GR), K. Schwarzschild solved the vacuum Einstein equations Schwarzschild (1916) describing exterior of a spherically symmetric and static sphere. The simple internal spherically symmetric solution with special uniform distribution of matter has been found by Schwarzschild Schwarzschild (1916), and generalized for spacetimes with non zero cosmological constant in Stuchlík (2000); Böhmer (2004). Some important internal solutions of the Einstein equations were found by R. Tolman for perfect fluid spheres with fluid described by polytropic equations of state Tolman (1939). The polytropes in spacetimes with non-zero cosmological constant were extensively discussed in Stuchlík (2000); Stuchlík et al. (2016). Interesting properties of the spherically symmetric polytropes were discussed in Stuchlík et al. (2016); Novotnỳ et al. (2017); Hod (2018a, b); Stuchlík et al. (2017).

Anisotropic pressure in stellar distribution implies unequal radial and tangential pressure. The possible reasons for anisotropies in fluid pressure are presence of mixture of different fluids, different kinds of phase transitions Sokolov et al. (1980), viscosity, rotation, magnetic field, superfluid Kippenhahn et al. (1990) or existence of a solid core. From our general notion, it is understandable that anisotropic solutions could represent a realistic description of astrophysical interest. Anisotropy of pressure in a perfect fluid sphere was first described by G. Lemaitre Lemaître (1933). In 1974, Bower and Liang Bowers and Liang (1974) described importance of locally anisotropic equations of state for relativistic spheres. Ruderman has shown in his significant work that nuclear matter may be anisotropic in very high density regions ($\rho>10^{17}\mathrm{kg/m^{3}}$) Ruderman (1972). A few articles devoted to generating anisotropic solutions are available in literature Thirukkanesh et al. (2018); Herrera et al. (2008); Chaisi and Maharaj (2006).

There are eight perfect fluid solutions of the Einstein field equations presented by Tolman Tolman (1939). Among these solutions, the Tolman IV and the Tolman VII solutions are physically interesting and could depict some neutron star configurations. For detailed analysis of the Tolman VII solution see Sotani and Kokkotas (2018); Kiess (2017); Bhar et al. (2017); Singh et al. (2016); Raghoonundun (2016); Raghoonundun and Hobill (2016a, b); Bhar et al. (2015); Raghoonundun and Hobill (2015); Papazoglou and Moustakidis (2016); Kiess (2012); Neary and Lake (2001); Neary et al. (2001).

Ovalle introduced an anisotropic version of the isotropic Tolman IV solution in Ovalle et al. (2018a), by using the so called minimal geometric deformation (MGD) method developed by himself Ovalle (2017); Ovalle et al. (2018a). Our objective in the present paper is to generate a new anisotropic solution from the isotropic Tolman VII solution (Tolman, 1939), by using the Ovalle MGD method.

The MGD method is the first simple, systematic and direct method of decoupling gravitational sources in GR. Initially, MGD method was proposed Ovalle (2008, 2010a) in the context of the Randall-Sundrum braneworld model Randall and Sundrum (1999a, b) –references for earlier works on the MGD method,  Ovalle (2009, 2010b); Casadio and Ovalle (2012); Ovalle and Linares (2013); Ovalle et al. (2013); Casadio et al. (2014), and for some recent applications see  Ovalle et al. (2015); Casadio et al. (2015); Cavalcanti et al. (2016); Casadio and da Rocha (2016); Ovalle et al. (2017); da Rocha (2017a, b); Fernandes-Silva and da Rocha (2018); Casadio et al. (2017); Ovalle et al. (2018a); Heras and Leon (2018); Fernandes-Silva et al. (2018); Estrada and Tello-Ortiz (2018); Contreras and Bargueño (2018); Morales and Tello-Ortiz (2018); Gabbanelli et al. (2018); Sharif and Sadiq (2018a, b); Graterol (2018); Ovalle et al. (2018b); Contreras (2019); Contreras et al. (2019); Maurya and Tello-Ortiz (2019); Fernandes-Silva et al. (2019); Contreras (2019); Ovalle (2019); Ovalle et al. (2018c); Sharif and Saba (2018)). The notable feature of the MGD method is that it preserves the spherical symmetry, as well as the physical acceptability. MGD method thus opens up a new window to search for physically acceptable anisotropic solutions.

We review first the formulation of decoupling the Einstein field equations for different gravitational sources. We have to solve

[TABLE]

and

[TABLE]

instead of looking for solution for the total energy momentum tensor $T_{\mu\nu}=\hat{T}_{\mu\nu}+\tilde{T}_{\mu\nu}$,

[TABLE]

where the constant $k^{2}=8\pi G/c^{4}$, and in the geometric units ($c=G=1$) there is $k^{2}=8\pi$. $T_{\mu\nu}$ is the total energy momentum tensor, $\hat{T}_{\mu\nu}$ is the energy momentum tensor for one gravitational source and $\tilde{T}_{\mu\nu}$ is the energy momentum tensor for another gravitational source. After solving Eqs. (1) and (2), we can find metric $g_{\mu\nu}$ by combining $\hat{g}_{\mu\nu}$ and $\tilde{g}_{\mu\nu}$.

Generally, we can extend this formalism for any number of gravitational sources. In that case we have to solve Einstein’s field equations for each source term, and combine the separately found metrics in order to get the metric related to the total energy momentum tensor. The number of physically acceptable solutions is not large. Delgaty and Lake examined physically acceptability of 127 known isotropic solutions  Delgaty and Lake (1998). They found that only 16 of them has physical relevance.

The paper is organized as follows. In Section II we discuss decoupling of Einstein’s field equations. In Section III we introduce the MGD method. In Section IV we study the condition for matching the interior solution to the exterior one. Section V is devoted to description of the Tolman VII perfect fluid solution. In Section VI we derive the new anisotropic solution and demonstrate its interesting properties. In Section VII conclusions are presented.

II Decoupling of Einstein’s field equations

We shortly review the method of the gravitational decoupling developed by Ovalle Ovalle (2017). In the framework of the decoupling method, Einstein’s field equations are expressed in the form

[TABLE]

where,

[TABLE]

The energy momentum tensor for perfect fluid reads

[TABLE]

where, $\rho$,  $p$ and $u_{\mu}$ are density, pressure and four velocity of perfect fluid, respectively. We consider the contribution from other gravitational source($\theta_{\mu\nu}$) modified by an intensity parameter $\alpha$

[TABLE]

The total energy momentum tensor for perfect fluid coupled with another gravitational source causing anisotropy in the self gravitating system then reads

[TABLE]

The term $\theta_{\mu\nu}$ in Eq. (7) stands for any source like scalar, vector, or tensor field, causing anisotropies in the fluid. As a consequence of Bianchi identity the conservation law holds

[TABLE]

The line element of spherically symmetric spacetime in Schwarzschild coordinates $\left(t,r,\theta,\phi\right)$ takes the form

[TABLE]

where, $\nu(r)$ and $\lambda(r)$ are functions of the radial coordinate ($r$) which ranges from the compact object centre ($r=0$) to its surface ($r=R$). Four-velocity of the static fluid is given by $u^{\mu}=e^{-\nu/2}\delta_{0}^{\mu}$, at radii $0\leq r\leq R$. The general metric given by Eq. (10) obeys the Einstein field equations having the following form

[TABLE]

[TABLE]

[TABLE]

The conservation Eq. (9), which is linear combination of Eqs. (11), (12) and (13), reads

[TABLE]

Here ′ denotes differentiation of the function with respect to $r$. We can easily identify from Eqs. (11), (12) and (13) the effective density

[TABLE]

the effective isotropic pressure

[TABLE]

and the effective tangential pressure

[TABLE]

It is evident that the $\theta_{\mu\nu}$ source introduces an anisotropy into perfect fluid which is given by

[TABLE]

Now, we have five unknowns, namely the metric coefficients $\nu(r)$, $\lambda(r)$, the effective density ($\rho_{\mathrm{eff}}$), the effective radial pressure ($p_{r}$), and the effective tangential pressure ($p_{t}$). To solve the Einstein equations to obtain these functions, we have to proceed using the MGD method.

III Minimal geometric deformation

MGD method is a very strong technique to decouple Einstein’s field equations for different source terms. This method has been developed in a simple and elegant way recently Ovalle (2017), we shortly review this method. We consider a perfect fluid solution of energy density $\rho$ and pressure $p$ that is described by the metric coefficients $\left(e^{\epsilon},e^{\gamma}\right)$. The line element of the corresponding spherically symmetric solution reads

[TABLE]

where the radial metric coefficient takes the form

[TABLE]

with $m(r)$ being so called GR mass function. We introduce the MGD transformations in the form Ovalle (2017)

[TABLE]

[TABLE]

where $g$ and $f$ are deformations of temporal and radial metric coefficients respectively, arising as an effect of introduction of the anisotropy. The minimal geometric deformation is given by the conditions

[TABLE]

and

[TABLE]

Then temporal component of the metric $e^{\nu}$ remains unchanged, while the additional gravitational source ($\theta_{\mu\nu}$) causes a deformation of the radial component according to Eq. (22). If we incorporate the deformed metric into Eqs. (11), (12) and (13), we see that each equation can be decomposed into two equations – one holds for the perfect fluid, the other one is involving $\theta_{\mu\nu}$. The set of equations for perfect fluid ($\alpha=0$) is fully determined by the metric coefficient $\epsilon(r)=\nu(r)$ and takes the form

[TABLE]

[TABLE]

[TABLE]

The conservation equation takes the form

[TABLE]

The equations involving the $\theta_{\mu\nu}$ term read

[TABLE]

[TABLE]

[TABLE]

The conservation equation for the $\theta_{\mu\nu}$ term takes the form

[TABLE]

and it is a linear combination of Eqs. (29), (30) and (31).

We see that Einstein’s equations (11), (12) and (13) are decoupled by deforming the radial metric component according to Eq. (21). The conservation Eqs. (28) and (32) for the MGD solution have to be satisfied simultaneously with the general conservation law given by Eq. (9). We conclude that both systems, perfect fluid and other gravitational source, conserve independently, i.e., these two systems cannot exchange energy-momentum, and their interaction is solely gravitational Ovalle (2017).

IV Matching condition

The matching of the interior ($r<R$) and exterior ($r>R$) solutions of any mass distribution at the boundary of the interior ($r=R$), has to be considered very carefully in the case of mixed sources treated in the framework of the MGD method Ovalle (2017). The interior metric of our consideration is given by Eq. (19) along with Eqs (20), (21), (22) and (23), and its line element can be written as

[TABLE]

where the interior mass function reads

[TABLE]

$m(r)$ is the standard GR mass function in Eq.(20), and the function $f^{*}$ has to be calculated later. We assume there is no matter outside, i.e., $\rho^{+}=p^{+}=0$. In general, additional gravitational source $\theta_{\mu\nu}$ could influence the exterior geometry, and in such a case the general exterior metric can be written as

[TABLE]

where the functional form of $\nu^{+}(r)$ and $\lambda^{+}(r)$ can be determined by solving Einstein’s equations for the exterior geometry

[TABLE]

Continuity of the first fundamental form implies

[TABLE]

From above we obtain equations

[TABLE]

and

[TABLE]

Here $M_{0}=m(R)$ and $f^{*}_{R}$ is the minimal geometric deformation at the boundary of the fluid distribution. Considering the Israel-Darmois matching condition at the surface ($r=R$), we get second fundamental form that reads

[TABLE]

where $r^{\mu}$ is the unit radial vector. Proportionality of $G_{\mu\nu}$ and $T_{\mu\nu}$ implies

[TABLE]

Eq. (41) can be rewritten with the help of Eq. (12) in the form

[TABLE]

As there is no fluid (matter) assumed outside the configuration, we can write

[TABLE]

where $p_{R}=p^{-}(R)$. With the use of Eq. (30) for the inner and the outer geometry, we can write

[TABLE]

here $g^{*}_{R}$ is the deformation of the outer geometry due to the matter term $\theta_{\mu\nu}$, and $M_{s}$ is the Schwarzschild mass function. So, for matching the interior geometry with the exterior one, Eqs (38), (39) and (44) are necessary and sufficient condition. If the outside geometry is the Schwarzschild vacuum one, Eq. (44) reduces to

[TABLE]

which implies that the effective radial pressure at the surface should vanish.

V Interior perfect fluid Tolman VII solution

Let us summarize properties of the well known interior Tolman VII perfect fluid solution Tolman (1939) for which we apply the MGD method. Its metric coefficients are given by

[TABLE]

[TABLE]

where A,B,C are constants of the solution, yet to be determined. The line element for this solution is described by Eq. (19) with $e^{\epsilon(r)}=e^{\gamma(r)}$. Using Eqs (25) and (26), we calculate radial profiles of the pressure and density of the perfect fluid

[TABLE]

and

[TABLE]

where $z=\left[\log\sqrt{\frac{\frac{2r^{2}}{A^{2}}-\frac{A^{2}}{4R^{2}}+b_{1}}{C}}\right]$, $b_{1}=\sqrt{1+\frac{4r^{4}}{A^{4}}-\frac{r^{2}}{R^{2}}}$. The constants A, B and C are calculated according to the matching conditions given by Eqs (37) and (40), under assumption of the outside Schwarzschild vacuum solution, and read

[TABLE]

[TABLE]

[TABLE]

where $\frac{M_{0}}{R}\leq 4/9$ , $b_{\mathrm{2}}=\sqrt{1-\frac{2M_{0}}{R}}$ , $b_{\mathrm{3}}=\left(\frac{R^{5}}{-2M_{0}+R}\right)^{\frac{3}{2}}$, and $M_{\mathrm{0}}=m(R)$ is the total mass of the fluid configuration given by Eq. (20). For the realistic fluid configurations, the perfect fluid density should be positive everywhere within $0<r<R$. For such configurations, the condition $\rho^{2}<\frac{3R^{2}}{5(R-2M_{0})}$ has to be satisfied for all values of $r$. The configurations that do not satisfy this condition are not physically acceptable. Examples of the density and pressure radial profiles of the Tolman VII solution are shown in Fig. 1 . We can see they resemble profiles attained for (neutron) stars and could be thus applied in astrophysical context. For detailed study and application of the Tolman VII solution see Tolman (1939).

VI Anisotropic Tolman VII solution by gravitational decoupling

Now let us turn-on the parameter $\alpha$ to get an anisotropic solution. According to the matching condition (Eq. (45)), the effective pressure should vanish at the boundary which implies from Eq. (12) that $p_{R}\sim\alpha(\theta_{1}^{1})_{R}$ . If we choose the deformation in the radial component of metric

[TABLE]

it essentially means that

[TABLE]

Eq.(54) is denoted as ‘mimic’ constraint by Ovalle Ovalle (2017), and it requires to effective pressure to vanish at the boundary. Using Eqs. (22) and (48), we get the radial metric component

[TABLE]

So, Eqs (46) and (55) represent the Tolman VII solution being minimally deformed by the gravitational source $\theta_{\mu\nu}$; the original Tolman VII perfect fluid solution can be recovered while $\alpha\rightarrow 0$. Now, we have a new anisotropic solution ((46) and (55)). Our task is to match the new solution with the exterior Schwarzschild vacuum metric. Due to the matching conditions (Eqs. (38) and (39)) we arrive to the relations

[TABLE]

and

[TABLE]

By using Eq. (20), the above equation gives the Schwarzschild mass($M_{s}$) due to the relation

[TABLE]

Considering the Schwarzschild vacuum outside, the second fundamental form (Eq. (43)) reads

[TABLE]

and as a consequence of the ‘mimic’ constraint (Eq.  (54)), it reduces to the condition

[TABLE]

Using Eq. (48) this implies for the constant $C$ the relation

[TABLE]

Using the expression for the Schwarzschild mass given in Eq. (58), we obtain from Eq. (56)

[TABLE]

from which we can determine $B$ while $C$ is given by Eq. (61). Eqs. (61) and (62) are necessary and sufficient conditions for matching the anisotropic interior solution with exterior Schwarzschild vacuum.

By using the ‘mimic’ constraint (Eq. (54)) in Eq. (16), from Eq. (48) we get the radial profile of the radial pressure in the form

[TABLE]

The radial profiles of the effective density and the tangential pressure are then given by the relations

[TABLE]

and

[TABLE]

here $\delta\rho$ is change in density and $\pi(r,\alpha)$ is measure of anisotropy, being defined as

[TABLE]

[TABLE]

where the parameters are determined by the relations

$Y_{1}=4[160r^{8}R^{4}+A^{8}(6r^{4}-8r^{2}R^{2}+3R^{4})+A^{4}(-64r^{6}R^{2}+44r^{4}R^{4})]$ ,

$Y_{2}=16A^{2}r^{2}R^{2}[20r^{4}R^{2}+A^{2}(-3r^{2}+R^{2})]$ ,

$Y_{3}=[-80r^{6}A^{4}-3A^{8}(r^{2}-R^{2})+4A^{4}(8r^{4}R^{2}-5r^{2}R^{4})]$ ,

$Y_{4}=8A^{4}r^{2}R^{4}$ ,

$Y_{5}=4r^{4}R^{2}+A^{4}(-2r^{2}+3R^{2})$ .

The radial profiles of the radial pressure, the tangential pressure and the effective density are shown in Fig. 2, and the effect of the anisotropy in dependence on radius and the coupling constant ($\alpha$) is shown in the Fig. 3.

VII Conclusions

We started from the Tolman VII perfect fluid solution and using the framework of MGD, we get a new anisotropic solution. We demonstrate that the effect of anisotropy is increasing with increasing coupling constant($\alpha$). We also demonstrate that with increasing radius the radial profile of the anisotropy effect increases reaching its maximum value and then it is decreasing to vanish at the edge of the configuration. We calculate the Schwarzschild mass, $M_{s}$ (from Eq. 58) which is same as the mass of the perfect fluid stellar distribution, $M_{0}$. So, in the case considered here, the anisotropy effect is not affecting the mass of the stellar distribution. The reason behind this is evident from the Fig. 2 where we see that the effective density crosses the perfect fluid density profile. The effect of trapping of null geodesics in the anisotropic solution has the same character as in the perfect fluid solution Neary et al. (2001), because the corresponding effective potential of the null geodesics depends on the metric components $g_{tt}$ and $g_{\phi\phi}$ on, which remain unchanged in comparison with the isotropic solution – only $g_{rr}$ changes in the case considered in our paper.

Acknowledgements

S.H. and Z.S. would like to acknowledge the institutional support of the Faculty of Philosophy and Science of the Silesian University in Opava, the internal student grant of the Silesian University Grant No. SGS/12/2019 and the Albert Einstein Centre for Gravitation and Astrophysics under the Czech Science Foundation Grant No. 14-37086.