Geometry of static \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$w=-1/5$$\end{document}w=-1/5 perfect fluid spheres in general relativity

We discuss the physical features of two recent classes of analytical solutions of the Einstein equations sourced by an exotic perfect fluid with equation of state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ P=-\rho /5$$\end{document}P=-ρ/5. These geometries depend on up to four parameters and are static and spherically symmetric. They describe compact spaces with naked central singularities.


Introduction
Recently, two new families of static and spherically symmetric solutions of the Einstein equations (without cosmological constant) were proposed by Semiz [1]. The matter source is a perfect fluid with constant barotropic equation of state P = −ρ/5, where ρ and P are the fluid energy density and pressure, respectively [1]. One would like to understand the physical nature of these solutions and assess whether they can be useful to model regions of stars, at least as toy models. The equation of state P = −ρ/5 is clearly unphysical, as one would be hard put to find realistic situations described by this fluid, but dark energy-like stars (and even phantom energy stars [2]) have been studied in the literature [3][4][5][6][7][8][9][10], as well as halos of exotic energy [11]. Although dark energy has pressure P < −ρ/3 and there are all indications that, if it is responsible for the present acceleration of our universe, it has equation of state P −ρ [12], our situation with P = −ρ/5 could still serve as a toy model for hypothetical objects formed by a negative pressure fluid. Moreover, from the mathematical point of view, simple solutions of the Einstein equations describing perfect fluids are relatively difficult to find. Although there are over one hundred analytical solutions of the Einstein equations sourced by perfect (and a e-mail: b.fazlpour@umz.ac.ir b e-mail: a.banijamali@nit.ac.ir c e-mail: vfaraoni@ubishops.ca (corresponding author) even imperfect) fluids that constitute potential candidates to model relativistic stars, or at least stellar regions [13,14], almost all of them turn out to be unphysical for one reason or another [14]. Here we examine the new solutions of [1] to understand their physical features (or lack thereof). These geometries are written in Buchdahl coordinates but it is more instructive from the physical point of view to rewrite them in terms of Schwarzschild-like coordinates, which we do here.
We follow the notation of Ref. [15]: the metric signature is −+++ and we use units in which the speed of light in vacuo c and Newton's constant G are unity, while κ ≡ 8π G to keep with Ref. [1].
Semiz's proposal consists of a four-parameter family of solutions of the Einstein equations with zero cosmological constant where R ab is the Ricci tensor of the metric g ab and R is the Ricci scalar. The matter source is a perfect fluid with stress-energy tensor where u a is the fluid 4-velocity and the equation of state is P = −ρ/5 [1]. These geometries are spherically symmetric and static in the appropriate coordinate range. There are two new classes of solutions in [1]: the most general family is parametrized by four constants (C 0 , C 1 , C 2 , C 3 ) and has line element with C 1 = 0, C 0 + C 1 r = 0, and where while dΩ 2 (2) ≡ dϑ 2 + sin 2 ϑ dϕ 2 is the line element on the unit 2-sphere. The energy density is [1] and is non-negative provided that which we assume in the following, while the limiting situation given by C 3 = 0 corresponds to vacuum. The solution for C 1 = 0 is not obtained continuously from Eqs. (3) and (5) in the limit C 1 → 0 but requires a separate discussion [1]. This second family is parametrized by the remaining three constants [1]: we begin by analyzing this second family (or "special solution" in the nomenclature of [1]) in the following section.
2 Special solution C 1 = 0 This 3-parameter (C 0 , C 2 , C 3 ) family of solutions is described by the line element [1] with C 0 = 0 and where In order to preserve the metric signature it must be g(r ) > 0 (if g(r ) becomes negative, the coordinates t and r switch their timelike and spacelike natures, as in the Schwarzschild geometry at the horizon r = 2m). We rewrite the line element (7) in terms of the areal radius R(r ) = g(r )/|C 0 |. This relation is inverted by first obtaining and solving for The argument of the square root in the right-hand side must be non-negative to keep r real, which gives the limitation on the range of the areal radius. The latter begins from zero at and then decreases until it vanishes again at The two coordinate charts r 1 ≤ r ≤ C 0 C 3 /2 and C 0 C 3 /2 ≤ r ≤ r 2 cover the same physical region 0 ≤ R ≤ R max . We restrict ourselves to r 1 ≤ r ≤ C 0 C 3 /2, in which d R/dr > 0, by choosing the negative sign in Eq. (2.5).
We write and, substituting the relation between differentials and using Eq. (14), the line element (7) becomes The equation ∇ c R∇ c R = g R R = 0 locating the apparent horizons (see, e.g., [16]) has R max as the only root, which is a single root and therefore there are no apparent horizons for R < R max (we discuss the physical meaning of the formal root R max below). The energy density [1] and the pressure P = −ρ/5 (which are always non-zero) diverge at the origin R = 0, which corresponds to r = r 1 , together with the Ricci scalar therefore there is a naked spacetime singularity at the origin R = 0.
The Misner-Sharp-Hernandez mass M MSH (R) defined in spherical symmetry by [17,18] for the geometry (16) and is negative in the region 0 < R < 4R max /5 around the naked singularity. This fact is not surprising: it has been argued that the Misner-Sharp-Hernandez mass (to which the Hawking-Hayward quasilocal mass [19,20] reduces in spherical symmetry [21]) is misbehaved near naked singularities, Cauchy horizons, or regions with the wrong asymptotic flatness [22,23]. This is the case, for example, for the inner region of the Reissner-Nordström black hole near the Cauchy horizon, for the entire Schwarschild spacetime with negative mass, and for the Fisher-Janis-Newman-Winicour-Buchdahl-Wyman scalar field solution of the Einstein equations [24-31] for the parameter values for which there is a naked singularity [32]. Let us come to the maximum value R max of the areal radius which, in spite of being a formal root of the equation ∇ c R∇ c R = 0, does not describe a horizon but is instead the antipode of the origin R = 0 in a compact space. To see this fact, it is instructive to study the behaviour of radial null geodesics in this geometry. Consider the outgoing (+) and ingoing (−) congruences of radial null geodesics with tangents l μ (±) = dx μ /dλ, where λ is an affine parameter along these curves. These tangents have components l μ (±) = l 0 , l 1 , 0, 0 and the normalization l since a null vector can be rescaled by a function, we can choose l 0 = 1 (which means choosing the coordinate time t as the affine parameter along these null geodesics), obtaining We then have the first order equations which integrate to where λ 0 is an integration constant. Unfortunately this relation cannot be inverted explicitly. Since (with the upper sign for outgoing and the lower one for ingoing radial geodesics), near the origin R = 0 it is d R/dt ∼ +∞ for outgoing and d R/dt ∼ −∞ for ingoing geodesics. Furthermore, d R/dt = 0 at R = R max . Outgoing radial null geodesics starting near the origin do so extremely fast but they slow down as they approach the maximum possible radius R max , which can only be reached with zero velocity (see Fig. 1). A null geodesic starting exactly at R max does so with zero velocity d R/dλ and remains there. Ingoing radial null geodesics starting near the maximum radius R max are slow and accelerate as they get closer to the central naked singularity, which they approach with infinite velocity We can also study radial timelike geodesics with tangents u μ = u 0 , u 1 , 0, 0 . The normalization u c u c = −1 gives with the upper sign for outgoing and the lower one for ingoing geodesics. The timelike Killing vector ξ a = (∂/∂t) a guaran- Fig. 1 The areal radius R versus the affine parameter λ along the radial null geodesics of the geometry (16), for the parameter values C 0 = C 2 = 2, C 3 = 3, and λ 0 = 1. Outgoing geodesics slow down as they approach R max , where they stop. Ingoing geodesics starting near R max do so extremely slowly but accelerate as they approach R = 0 tees the conservation of the energy per unit mass of the test particle E along these geodesic curves: where u 0 > 0 because these curves are future-oriented, hence E is strictly positive. Equation (28) then gives which tells us that: • For a given energy E determined by the initial condition R 0 ,Ṙ 0 , radial motion is only possible if (otherwise u 1 becomes imaginary). Ingoing radial motion stops at R min and a test particle cannot approach the origin, which is consistent with the fact that, according to Eq. (29), u 0 = E R/|C 0 | → 0 as R → 0. • Outgoing radial motion stops at R max , where u 1 vanishes for both outgoing and ingoing radial geodesics, and a particle starting initially at R max remains there irrespective of its initial energy. • Since R is limited by R max , the possible energies are bounded from below, particles with energy below, or equal to, the minimum threshold E min do not move.

Case
In this case we are left with only two parameters (C 0 , C 2 ). Now g(r ) = C 0 C 2 − r 2 , which requires The areal radius is with r = 0 corresponding to R = C 0 C 2 /|C 0 |, while r = C 2 sign(C 0 ) corresponds to the origin R = 0 of the physical radial coordinate. The areal radius R(r ) varies in the range (where, in the last equality, we used the fact that C 0 C 2 > 0) and is a decreasing function of r since d R/dr = −2r/|C 0 | is always negative in the allowed range. Inverting the relation between radial coordinates, one obtains which, in conjunction with yields the line element This geometry is the same as that of the previous case C 1 = 0, C 3 = 0 given by the line element (16), but now R max = |C 2 |. Again, the energy density is non-zero and the Ricci scalar diverges at the origin R = 0.

Case
For these parameter values, g(r ) = r (C 0 C 3 − r ) requires C 0 C 3 to be positive and, therefore, we have the range 0 ≤ r ≤ C 0 C 3 of the Buchdahl radius. Correspondingly, the areal radius varies in the interval beginning from zero at r = 0, increasing to the maximum and then decreasing until it vanishes again at r = C 0 C 3 . There are two coordinate charts 0 ≤ r ≤ C 0 C 3 /2 and C 0 C 3 /2 ≤ r ≤ C 0 C 3 covering the same physical region 0 ≤ R ≤ R max and we restrict ourselves to the former, in which d R/dr > 0. Equation (39) yields where we choose the lower sign for consistency with d R/dr > 0 and 0 ≤ r ≤ C 0 C 3 /2. Then g(r ) = |C 0 |R and give the line element which is the same as the line element (16), but with R max now given by Eq. (41). The Ricci scalar diverges at the origin R = 0 (which corresponds to r = 0 in the chart with d R/dr > 0), therefore there is a naked spacetime singularity there.

General solution C 1 = 0
The line element for the generic family of Semiz solutions is (3) [1]. The presence of four parameters with relatively wide ranges now makes it difficult to reach definite conclusions and we focus on special cases.

C 3 = 0 is Schwarzschild
When C 3 = 0, the energy density (5) and the pressure P = −ρ/5 vanish identically and this spacetime is empty.
Since the geometry is also spherically symmetric and asymptotically flat (as we are going to show) it must be the Schwarzschild one, according to the Jebsen-Birkhoff theorem [15]. In fact, we have f (r ) = 3 (C 1 C 2 + r ), the areal radius is and then dr = |C 1 |d R, yielding the line element where dt ≡ |C 1 |dt. This geometry is asymptotically flat: by introducing the constant (which is not necessarily positive) and rescaling the time coordinate according to t →t = |C 1 | t, the line element (49) is written as the Schwarzschild one describing a black hole if m > 0 and a naked central singularity if m < 0. 3 ; the areal radius is

Special case
We have which is positive for To proceed, let us consider the situation r ≥ −C 0 /C 1 , in which case R increases in the interval with R min ≤ R ≤ R max and We have again a compact space. Rewriting the line element (3) in terms of the areal radius produces a cumbersome expression that does not depend only on R but contains also r (R) because the relation R(r ) cannot be inverted explicitly.
3.3 The even more special case C 0 = C 2 = 0 In this case we have only two parameters (C 1 , C 3 ), f (r ) = r 3 + C 3 1 C 3 r 2 , and the areal radius is where the last inequality follows from C 1 C 3 ≤ 0. Since the areal radius is an increasing function of r in the interval 0, and then decreases, vanishing again at r = 3 |C 1 C 3 |C 2 1 . This compact space corresponds to the range of the areal radius, with R r/|C 1 | as r → 0 + . Equation (59) is inverted by first obtaining and then solving for where while the two remaining roots are imaginary. Substituting the relation between differentials and using yield the line element where Again, the many combinations of parameters and the cumbersome metric coefficients do not lend themselves to a straightforward and transparent analysis, but it is clear that also in this case we have a compact 3-space of finite extent. Using (3.20), the energy density (5) reduces to as R → 0 + . Hence, ρ and P are singular at the origin, together with the Ricci scalar R and As R → 0 + , we have the asymptotics and g 00 → ∞ while g 11 → 0 as R → 0.

Conclusions
We have studied the nature of the new classes of static and spherically symmetric solutions of the Einstein equations given recently in Ref. [1] when the matter source is a perfect fluid with equation of state P = −ρ/5. The analytical solutions of Ref.
[1] that we analyzed (except for the Schwarzschild solution obtained for C 3 = 0) describe compact spaces with naked central singularities. The "general" family of solutions (3) and (5) always reduces to Schwarschild when the parameter C 3 vanishes. In most other situations, the presence of three or four parameters and/or the cubic nature of the function R(r ) hamper a complete description of the geometry. However, in all cases analyzed, except for the empty spacetime associated with C 3 = 0, we find a compact space of finite volume (a feature mentioned in [1]). The fact that the geometry, together with the energy density and the pressure, is singular at R = 0 is not necessarily the death knell for these solutions. In fact, it is deemed acceptable for fluid solutions of the Einstein equations to only model limited regions of relativistic stars, a procedure that is reflected in the authoritative Ref. [13] and in the more specialized literature. Indeed, even Newtonian stars are rarely modelled with a single fluid, corresponding to the fact that different regions at different temperatures and densities are described by different equations of state unless the stellar material is well mixed, which only happens in certain types of stars. Therefore, there is in principle the (physically well motivated) possibility of excising the singularity and replacing it with a more realistic geometry sourced by matter with a different equation of state. However, if one wants to describe a stellar interior with this exotic fluid, one must match it with an asymptotically flat Schwarzschild exterior. The fact that the solutions of [1] describe compact spaces points to a possible analogy with the Oppenheimer-Snyder model of gravitational collapse to a black hole [33]. In this model, a compact, positively curved Friedmann-Lemaître-Robertson-Walker universe collapsing to a Big Crunch is matched to a Schwarzschild exterior on the surface of a 2-sphere of symmetry [33], satisfying the Darmois-Israel junction conditions [34,35]. However, in the Oppenheimer-Snyder model the matching is possible because the collapsing interior universe is filled by a dust with zero pressure everywhere. It is well known that the matching to a Schwarzschild exterior can only be done on a surface on which the pressure P(R) vanishes, otherwise the junction conditions are violated and there is a material layer on the matching surface, which is certainly not an ingredient of realistic stellar models. (This fact is highlighted in many studies of relativistic fluid balls [36][37][38][39][40][41][42][43] and fireballs [44].) However, for the fluid solutions of [1] under discussion, the pressure P(R) never vanishes. Therefore, the best that one could do is modelling a limited region of a stellar interior with the Semiz solutions for P = −ρ/5. To be physical, this region should correspond to a positive Misner-Sharp-Hernandez mass M MSH and, therefore, should be sufficiently far away from the singularity at R = 0. The excised region containing the origin should be modelled with a different, non-singular, solution of the Einstein equations. 1 Then, the w = −1/5 solution should be matched continuously with another "intermediate" solution with non-vanishing pressure on a surface of constant radius, and the pressure in this layer should then go to zero at larger radii to make it possible to match it to a Schwarzschild exterior, satisfying again the Darmois-Israel junction conditions. In the absence of a specific need for such an involved "star" model in astrophysics, we will not pursue this object further, limiting ourselves to pointing out the constraints for such a construction. Probably some of the phenomenology unveiled here for the geometries found in [1] also applies to other classes of perfect fluid solutions of the Einstein equations. Whether this is the case will be established in future work.