Evolution of Self-Gravitating Fluid Spheres Involving Ghost Stars

Abstract

Exact solutions are presented which describe, either the evolution of fluid distributions corresponding to a ghost star (vanishing total mass), or describing the evolution of fluid distributions which attain the ghost star status at some point of their lives. The first two solutions correspond to the former case, they admit a conformal Killing vector (CKV) and describe the adiabatic evolution of a ghost star. Other two solutions corresponding to the latter case are found, which describe evolving fluid spheres absorbing energy from the outside, leading to a vanishing total mass at some point of their evolution. In this case the fluid is assumed to be expansion–free. In all four solutions the condition of vanishing complexity factor was imposed. The physical implications of the results, are discussed.

1. Introduction

In a recent paper [1] we have explored the possibility of finding fluid distributions of anisotropic fluids whose total mass is zero. Such configurations which we named ghost stars, are characterized by the presence of negative energy-density within the distribution, as consequence of which the total mass vanishes. The solutions presented in [1] are static, and may be regarded as either the final or the initial state of a dynamical process.

At this point the following remarks are in order

• As mentioned above the vanishing total mass of ghost stars is the result of the appearance of negative energy-density within some regions of the fluid distribution. Therefore such objects should be clearly differentiated from “evaporated” stars which have radiated all its energy away leading to the vanishing of its energy-density.
• The observational data collected from compact objects relies almost exclusively on two physical phenomena, namely: surface gravitational redshift and the influence of the gravitational field of the object on light rays passing close to it (its shadow). Both of which of course depend on the gravitational surface potential of the object (vanishing for any ghost star), which illustrates the potential observational consequences of such objects (see further comments on this issue in the last section).

In this work we endeavor to describe the evolution of fluid distributions describing the following possible scenarios

• The adiabatic evolution of a ghost star. The total mass remains zero all along the evolution.
• The non–adiabatic evolution of fluid distribution reaching at some point a zero total mass (with non-vanishing energy-density).

For the first scenario we obtain two solutions admitting a CKV. One of them corresponds to a CKV orthogonal to the four–velocity vector of the fluid, whereas the other admits a CKV parallel to the four–velocity.

The motivation behind the admittance of CKV is provided by the fact that it generalizes the well known concept of self–similarity, which play a very important role in classical hydrodynamics.

This explains the great interest aroused by this kind of symmetry (for very recent developments on this issue see, for example, refs. [2,3,4,5,6,7,8,9,10,11,12] and references therein).

This first couple of solutions represents fluid distributions contracting from arbitrarily large areal radius to compact fluid distribution with finite areal radius, or fluid distributions expanding from some finite value of the areal radius to infinity, always satisfying the vanishing total mass condition. The Darmois conditions at the boundary are satisfied, and as expected the energy-density is negative in some regions of the fluid distribution.

In the second scenario both solutions satisfy the vanishing expansion scalar condition. It is worth recalling that such a condition implies the existence of a cavity around the center, thereby helping in the modeling of voids observed at cosmological scales.

Newest results regarding expansion–free fluids may be found in [13,14,15,16,17,18,19,20,21,22,23,24] and references therein.

The two solutions satisfying the expansion–free condition represent fluid distribution contracting from arbitrarily large configurations to a singularity. They evolve absorbing energy from the outside, and at some point of their evolution their total mass vanishes. Thus the solutions pass through a ghost star state.

The presented solutions satisfy additional conditions which allow the full integration of field equations. These restrictions are the vanishing complexity factor condition as defined in [25,26] and the quasi–homologous evolution [27]. At any rate it is worth recalling that none of the above conditions is intrinsically related to the ghost star condition, and are assumed with the sole purpose to allow analytical integration of the Einstein equations.

A discussion on all the presented models is brought out in the last section.

2. The Relevant Equations and Variables

The system we are dealing with consists in a spherically symmetric distribution of collapsing fluid, bounded by a spherical surface $\mathrm{\Sigma }$. The fluid is assumed to be locally anisotropic (principal stresses unequal) and undergoing dissipation in the form of heat flow (diffusion approximation). We shall proceed now to summarize the definitions and main equations required for describing spherically symmetric dissipative fluids. A detailed description may be found in [27].

2.1. The Metric, the Energy–Momentum Tensor, the Kinematic Variables and the Mass Function

Choosing comoving coordinates, the general interior metric can be written as

$$d{s}^2=-A^{2}d{t}^2+B^{2}d{r}^2+R^2(d{\theta }^2+{sin}^2\theta d{\varphi }^2).$$

(1)

The energy momentum tensor in the canonical form, reads

$$T_{\alpha \beta }=\mu V_{\alpha }{V}_{\beta }+P{h}_{\alpha \beta }+{\mathrm{\Pi }}_{\alpha \beta }+q\left(V_{\alpha }{K}_{\beta }+K_{\alpha }{V}_{\beta }\right),$$

(2)

with

$$P={\displaystyle \frac{P_r+2P_{\perp }}{3}},h_{\alpha \beta }=g_{\alpha \beta }+V_{\alpha }{V}_{\beta },$$

$${\mathrm{\Pi }}_{\alpha \beta }=\mathrm{\Pi }\left(K_{\alpha }{K}_{\beta }-{\displaystyle \frac{1}{3}}{h}_{\alpha \beta }\right),\mathrm{\Pi }=P_r-P_{\perp },$$

where $\mu$ is the energy-density, $P_r$ the radial pressure, $P_{\perp }$ the tangential pressure, $q^{\alpha }=q{K}^{\alpha }$ the heat flux, $V^{\alpha }$ the four–velocity of the fluid, and $K^{\alpha }$ a unit four–vector along the radial direction. For comoving observers, we have

$$\begin{array}{ccc}\hfill V^{\alpha }& =& A^{-1}{\delta }_0^{\alpha },K^{\alpha }=B^{-1}{\delta }_1^{\alpha },\hfill \end{array}$$

(3)

satisfying

$$\begin{array}{c}\hfill V^{\alpha }{V}_{\alpha }=-1,V^{\alpha }{q}_{\alpha }=0,K^{\alpha }{K}_{\alpha }=1,K^{\alpha }{V}_{\alpha }=0.\end{array}$$

(4)

Both bulk and shear viscosity, as well as dissipation in the free streaming approximation can be trivially absorbed in $P_{\perp }$, $\mu ,P_r$ and q.

The Einstein equations for (1) and (2), are explicitly written in Appendix A.

The acceleration $a_{\alpha }$ and the expansion $\mathrm{\Theta }$ of the fluid are given by

$$a_{\alpha }=V_{\alpha ;\beta }{V}^{\beta },\mathrm{\Theta }={V^{\alpha }}_{;\alpha },$$

(5)

and its shear ${\sigma }_{\alpha \beta }$ by

$${\sigma }_{\alpha \beta }=V_{(\alpha ;\beta )}+a_{(\alpha }{V}_{\beta )}-{\displaystyle \frac{1}{3}}\mathrm{\Theta }{h}_{\alpha \beta }.$$

(6)

From the equations above we have for the four–acceleration and its scalar a,

$$a_{\alpha }=a{K}_{\alpha },a={\displaystyle \frac{A^{\prime }}{AB}},$$

(7)

and for the expansion

$$\mathrm{\Theta }={\displaystyle \frac{1}{A}}\left({\displaystyle \frac{\dot{B}}{B}}+2{\displaystyle \frac{\dot{R}}{R}}\right),$$

(8)

while for the shear we obtain

$${\sigma }^{\alpha \beta }{\sigma }_{\alpha \beta }={\displaystyle \frac{2}{3}}{\sigma }^2,$$

(9)

where

$$\sigma ={\displaystyle \frac{1}{A}}\left({\displaystyle \frac{\dot{B}}{B}}-{\displaystyle \frac{\dot{R}}{R}}\right),$$

(10)

in the above prime stands for r differentiation and the dot stands for differentiation with respect to t.

Next, the mass function $m(t,r)$ reads [28]

$$m={\displaystyle \frac{R^3}{2}}{R_{23}}^{23}={\displaystyle \frac{R}{2}}\left[{\left({\displaystyle \frac{\dot{R}}{A}}\right)}^2-{\left({\displaystyle \frac{R^{\prime }}{B}}\right)}^2+1\right].$$

(11)

Introducing the proper time derivative $D_T$ given by

$$D_T={\displaystyle \frac{1}{A}}{\displaystyle \frac{\partial }{\partial t}},$$

(12)

we can define the velocity U of the collapsing fluid as the variation of the areal radius with respect to proper time, i.e.,

$$U=D_{T}R,$$

(13)

where R defines the areal radius of a spherical surface inside the fluid distribution (as measured from its area).

Then (11) can be rewritten as

$$E\equiv {\displaystyle \frac{R^{\prime }}{B}}={\left(1+U^2-{\displaystyle \frac{2m}{R}}\right)}^{1/2}.$$

(14)

An alternative expression for m may be found using field equations, it reads

$$m=4\pi {\int }_0^r\left(\mu +q{\displaystyle \frac{U}{E}}\right)R^2{R}^{\prime }dr,$$

(15)

satisfying the regular condition $m(t,0)=0$.

From the above equation it follows that, since $R^{\prime }>0$ in order to avoid shell crossing singularities, the vanishing total mass condition ($m\left(r_{\mathrm{\Sigma }}\right)=0$) requires that the “effective energy-density” ($\mu +q{\displaystyle \frac{U}{E}}$) should be either zero (the trivial case), or changes its sign within the fluid distribution (ghost star).

2.2. The Complexity Factor

The complexity factor is scalar that measures the degree of complexity of a given fluid distribution [25,26]. It is identified with the scalar function $Y_{TF}$ which defines the trace–free part of the electric Riemann tensor (see [25,26] for details).

Using field equations we obtain

$$Y_{TF}=-8\pi \mathrm{\Pi }+{\displaystyle \frac{4\pi }{R^3}}{\int }_0^r{R}^3\left({\mu }^{\prime }-3q{\displaystyle \frac{UB}{R}}\right)dr.$$

(16)

In terms of the metric functions the scalar $Y_{TF}$ reads

$$\begin{array}{c}\hfill Y_{TF}={\displaystyle \frac{1}{A^2}}\left[{\displaystyle \frac{\ddot{R}}{R}}-{\displaystyle \frac{\ddot{B}}{B}}+{\displaystyle \frac{\dot{A}}{A}}\left({\displaystyle \frac{\dot{B}}{B}}-{\displaystyle \frac{\dot{R}}{R}}\right)\right]\\ \hfill +{\displaystyle \frac{1}{B^2}}\left[{\displaystyle \frac{A^{″}}{A}}-{\displaystyle \frac{A^{\prime }}{A}}\left({\displaystyle \frac{B^{\prime }}{B}}+{\displaystyle \frac{R^{\prime }}{R}}\right)\right].\end{array}$$

(17)

2.3. The Homologous and Quasi–Homologous Conditions

In the dynamic case, the discussion about the complexity of a fluid distribution involves not only the complexity factor which describes the complexity of the structure of the fluid, but also the complexity of the pattern of evolution.

Following previous works [26,27] we shall consider two specific modes of evolution as the most suitable candidates to describe the simplest pattern of evolution. These are, the homologous evolution (H) [26] characterized by

$$U=\tilde{a}\left(t\right)R,\tilde{a}\left(t\right)\equiv {\displaystyle \frac{U_{\mathrm{\Sigma }}}{R_{\mathrm{\Sigma }}}},$$

(18)

and

$${\displaystyle \frac{R_1}{R_2}}=\mathrm{constant},$$

(19)

where $R_1$ and $R_2$ denote the areal radii of two concentric shells ($1,2$) described by $r=r_1=\mathrm{constant}$, and $r=r_2=\mathrm{constant}$, respectively.

A somehow softer condition is represented by the quasi–homologous condition (QH) [27], which only requires the fulfillment of (18).

It can be shown, using the field equations (see [26,27] for details) that (18) implies

$${\displaystyle \frac{4\pi }{R^{\prime }}}Bq+{\displaystyle \frac{\sigma }{R}}=0.$$

(20)

Thus, H condition implies (19) and (20), whereas QH condition only implies (20).

2.4. The Exterior Spacetime and Junction Conditions

Since our fluid distribution is bounded we assume that outside $\mathrm{\Sigma }$ the space–time is described by Vaidya metric which reads.

$$d{s}^2=-\left[1-{\displaystyle \frac{2M\left(v\right)}{\mathbf{r}}}\right]d{v}^2-2d\mathbf{r}dv+{\mathbf{r}}^2(d{\theta }^2+{sin}^2\theta d{\varphi }^2),$$

(21)

where $M\left(v\right)$ denotes the total mass, and v is the retarded time.

The smooth matching of the full nonadiabatic sphere to the Vaidya spacetime, on the surface $r=r_{\mathrm{\Sigma }}=$ constant, requires the fulfillment of the Darmois conditions, i.e., the continuity of the first and second fundamental forms across $\mathrm{\Sigma }$ (see [29] and references therein for details), which implies

$$m(t,r)\stackrel{\mathrm{\Sigma }}{=}M\left(v\right),$$

(22)

and

$$q\stackrel{\mathrm{\Sigma }}{=}{P}_r,$$

(23)

where $\stackrel{\mathrm{\Sigma }}{=}$ means that both sides of the equation are evaluated on $\mathrm{\Sigma }$.

When Darmois conditions are not satisfied the boundary surface is a thin shell.

2.5. The Transport Equation

In the case of non–adiabatic evolution we have to resort to some transport equation to describe the evolution and spatial distribution of the temperature. Thus for example within the context of the Israel–Stewart theory [30,31,32] the transport equation for the heat flux reads

$$\begin{array}{ccc}\hfill \tau h^{\alpha \beta }{V}^{\gamma }{q}_{\beta ;\gamma }+q^{\alpha }& =& -\kappa h^{\alpha \beta }\left(T_{,\beta }+T{a}_{\beta }\right)\hfill \\ & -& {\displaystyle \frac{1}{2}}\kappa T^2{\left({\displaystyle \frac{\tau V^{\beta }}{\kappa T^2}}\right)}_{;\beta }{q}^{\alpha },\hfill \end{array}$$

(24)

where $\kappa$ denotes the thermal conductivity, and T and $\tau$ denote temperature and relaxation time respectively.

In the spherically symmetric case under consideration, the transport equation has only one independent component which may be obtained from (24) by contracting with the unit spacelike vector $K^{\alpha }$, we get

$$\tau V^{\alpha }{q}_{,\alpha }+q=-\kappa \left(K^{\alpha }{T}_{,\alpha }+Ta\right)-{\displaystyle \frac{1}{2}}\kappa T^2{\left({\displaystyle \frac{\tau V^{\alpha }}{\kappa T^2}}\right)}_{;\alpha }q.$$

(25)

Sometimes it is possible to simplify the equation above, in the so called truncated transport equation, when the last term in (24) may be neglected [33], producing

$$\tau V^{\alpha }{q}_{,\alpha }+q=-\kappa \left(K^{\alpha }{T}_{,\alpha }+Ta\right).$$

(26)

3. Exact Solutions

We shall now proceed to present exact analytical solutions describing two different scenarios,

• the evolution of a ghost star keeping its nature ($M=0$) all along its evolution
• the evolution of a compact object (not a ghost star) attaining the ghost star condition at some moment of its life.

We shall consider, both, the non–dissipative and the dissipative case.

In order to specify our models we need to impose some further restrictions. In this work such restriction will be

• The admittance of a conformal killing vector (CKV),

or

• The expansion–free condition.

In some cases the above conditions have to be complemented with additional restrictions such as

• The vanishing complexity factor condition.
• The homologous or the quasi–homologous approximation.

3.1. Solutions Admitting a CKV

In this subsection we shall consider spacetimes satisfying the equation

$${\mathcal{L}}_{\chi }{g}_{\alpha \beta }=2\psi g_{\alpha \beta },$$

(27)

where ${\mathcal{L}}_{\chi }$ denotes the Lie derivative with respect to the vector field $\mathit{\chi }$, which unless specified otherwise, has the general form

$$\mathit{\chi }=\xi (t,r,){\partial }_t+\lambda (t,r,){\partial }_r,$$

(28)

and $\psi$ in principle is a function of t and r. The case $\psi =constant$ corresponds to a homothetic Killing vector (HKV). The solutions described here are particular cases of solutions found in [8].

We shall consider two possible subclasses, both of which describe non–dissipative evolution

• ${\chi }^{\alpha }$ orthogonal to $V^{\alpha }$,
• ${\chi }^{\alpha }$ parallel to $V^{\alpha }$.

In the first case (${\chi }^{\alpha }$ orthogonal to $V^{\alpha }$), we shall obtain from the matching conditions, the QH condition and the vanishing complexity factor condition, with $M=0$, solution I.

In the second case (${\chi }^{\alpha }$ parallel to $V^{\alpha }$), we shall obtain from the matching conditions and the vanishing complexity factor condition, solution $II$.

Let us start by considering the case ${\chi }^{\alpha }$ orthogonal to $V^{\alpha }$, and $q=0$.

3.1.1. Solution I: ${\chi }_{\alpha }{V}^{\alpha }=q=M=0$

In this case we obtain from (27) (see [8] for details)

$$A=\alpha R={\displaystyle \frac{F\left(t\right)}{f\left(t\right)+g\left(r\right)}},B={\displaystyle \frac{1}{f\left(t\right)+g\left(r\right)}},$$

(29)

where f and g are two arbitrary functions of their arguments and $\alpha$ is a unit constant with dimensions of $\left[{\displaystyle \frac{1}{r}}\right]$.

Thus any model is determined up to three arbitrary functions $F\left(t\right),f\left(t\right),g\left(r\right)$, in terms of which the field equations read

$$\begin{array}{ccc}\hfill 8\pi \mu & =& {\displaystyle \frac{{(f+g)}^2}{F^2}}\left[{\displaystyle \frac{{\dot{F}}^2}{F^2}}-{\displaystyle \frac{4\dot{F}\dot{f}}{F(f+g)}}+{\displaystyle \frac{3{\dot{f}}^2}{{(f+g)}^2}}\right]\hfill \\ & +& 2g^{″}(f+g)-3g^{\prime 2}+{\displaystyle \frac{{\alpha }^2{(f+g)}^2}{F^2}},\hfill \end{array}$$

(30)

$$\begin{array}{c}\hfill 8\pi P_r={\displaystyle \frac{{(f+g)}^2}{F^2}}\left[{\displaystyle \frac{{\dot{F}}^2}{F^2}}+{\displaystyle \frac{2\dot{F}\dot{f}}{F(f+g)}}\right.\\ \hfill \left.-{\displaystyle \frac{3{\dot{f}}^2}{{(f+g)}^2}}+{\displaystyle \frac{2\ddot{f}}{f+g}}-{\displaystyle \frac{2\ddot{F}}{F}}\right]\\ \hfill +3g^{\prime 2}-{\displaystyle \frac{{\alpha }^2{(f+g)}^2}{F^2}},\end{array}$$

(31)

$$\begin{array}{ccc}\hfill 8\pi P_{\perp }& =& {\displaystyle \frac{{(f+g)}^2}{F^2}}\left[{\displaystyle \frac{{\dot{F}}^2}{F^2}}-{\displaystyle \frac{3{\dot{f}}^2}{{(f+g)}^2}}+{\displaystyle \frac{2\ddot{f}}{f+g}}-{\displaystyle \frac{\ddot{F}}{F}}\right]\hfill \\ & +& 3g^{\prime 2}-2g^{″}(f+g).\hfill \end{array}$$

(32)

Next, the matching conditions (22) and (23) on the surface $r=r_{\mathrm{\Sigma }}=constant$ read

$${\dot{R}}_{\mathrm{\Sigma }}^2+{\alpha }^2(R_{\mathrm{\Sigma }}^2-2M{R}_{\mathrm{\Sigma }}-\omega R_{\mathrm{\Sigma }}^4)=0,$$

(33)

and

$$2{\ddot{R}}_{\mathrm{\Sigma }}{R}_{\mathrm{\Sigma }}-{\dot{R}}_{\mathrm{\Sigma }}^2-{\alpha }^2(3\omega R_{\mathrm{\Sigma }}^4-R_{\mathrm{\Sigma }}^2)=0,$$

(34)

with $\omega \equiv g^{\prime }{\left(r_{\mathrm{\Sigma }}\right)}^2$.

Since (33) is just the first integral of (34), boundary conditions provide only one additional equation.

In order to specify a solution we still need to impose two additional conditions.

One of these conditions will be the quasi–homologous condition which implies because of (20) that the fluid is shear–free ($\sigma =0$), implying in its turn

$${\displaystyle \frac{\dot{B}}{B}}={\displaystyle \frac{\dot{R}}{R}}\Rightarrow F\left(t\right)=Const.\equiv F_0.$$

(35)

Thus the metric functions become

$$A={\displaystyle \frac{F_0}{f\left(t\right)+g\left(r\right)}},B={\displaystyle \frac{1}{f\left(t\right)+g\left(r\right)}},R={\displaystyle \frac{F_0}{\alpha \left[f\left(t\right)+g\left(r\right)\right]}}.$$

(36)

In order to determine $g\left(r\right)$, we shall further impose the vanishing complexity factor condition ($Y_{TF}=0$), producing

$$g\left(r\right)=c_{1}r+c_2,$$

(37)

with $c_1\equiv -\sqrt{\omega }$ (we have chosen the minus sign in the definition of $c_1$ to satisfy the condition $R^{\prime }>0$), and $c_2$ is another integration constant.

Then from the condition $M=0$, we obtain a solution to (33) which reads (see [8] for details)

$$R_{\mathrm{\Sigma }}^{\left(I\right)}={\displaystyle \frac{1}{\sqrt{\omega }cos\left[\alpha (t-t_0)\right]}}.$$

(38)

Finally using (38) in (36) we obtain the explicit form of $f\left(t\right)$ for this solution

$$f^{\left(I\right)}\left(t\right)={\displaystyle \frac{1}{\alpha }}{F}_0\sqrt{\omega }cos\left[\alpha (t-t_0)\right]+\sqrt{\omega }{r}_{\mathrm{\Sigma }}-c_2.$$

(39)

Thus, the corresponding physical variables read

$$\begin{array}{ccc}\hfill 8\pi \mu & =& -2\omega {cos}^2\left[\alpha (t-t_0)\right]\hfill \\ & +& \omega \left(1-{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right)\left[1-{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}+2cos\left[\alpha (t-t_0)\right]\right],\hfill \end{array}$$

(40)

$$8\pi P_r=-\omega \left(1-{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right)\left[4cos\left[\alpha (t-t_0)\right]+1-{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right],$$

(41)

$$8\pi P_{\perp }=\omega {cos}^2\left[\alpha (t-t_0)\right]-2\omega \left(1-{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right)cos\left[\alpha (t-t_0)\right].$$

(42)

The graphics of these physical variables for the solution I are given in Figure 1.

Solution I describes an expanding sphere, whose initial boundary areal radius grows from $1/\sqrt{\omega }$ at $\alpha (t-t_0)=0$, to infinity as $\alpha (t-t_0)\to \pi /2$, and a contracting sphere whose boundary areal radius decreases from infinity at $\alpha (t-t_0)=3\pi /2$ to $1/\sqrt{\omega }$ at $\alpha (t-t_0)=2\pi$. This picture repeating each time interval $\alpha (t-t_0)=2n\pi$, for any positive real integer n.

In order to determine the regions of the fluid distribution where the energy-density is negative (required to have a vanishing total mass) we shall write the condition $\mu =0$ from (40) in the form

$$z^2-2z(1+cosy)+2cosy(1-cosy)+1=0,$$

(43)

whose solution reads

$$z=1+cosy-\sqrt{3}cosy,$$

(44)

with $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ and $y\equiv \alpha (t-t_0)$.

The curve in Figure 2 is formed by all points in the plane $z,y$ where the energy-density vanishes. The curve divides the plane in two regions, corresponding to negative and positive values of $\mu$. As is apparent from the graphic of $\mu$ in Figure 1, these regions are denoted by A and B respectively, in Figure 2.

3.1.2. Solution II: $q=0;{\chi }^{\alpha }\left|\right|V^{\alpha }$

We shall next analyze the case when the CKV is parallel to the four–velocity vector in the absence of dissipation. In this case the Equation (27) produces

$$A=Bh\left(r\right),R=Br,{\chi }^0=1,\psi ={\displaystyle \frac{\dot{B}}{B}},$$

(45)

where $h\left(r\right)$ is an arbitrary function of its argument and ${\chi }^1=0$. It is worth noticing that in this case the fluid is necessarily shear–free, implying thereby that it evolves in QH regime.

Figure 2. $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ as function of $y\equiv \alpha (t-t_0)$ for the condition $\mu =0$. Regions A and B correspond to negative and positive values of $\mu$ respectively, for solution I.

Figure 2. $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ as function of $y\equiv \alpha (t-t_0)$ for the condition $\mu =0$. Regions A and B correspond to negative and positive values of $\mu$ respectively, for solution I.

Thus the line element may be written as

$$d{s}^2=B^2\left[-h^2\left(r\right)d{t}^2+d{r}^2+r^2(d{\theta }^2+{sin}^2\theta d{\varphi }^2)\right].$$

(46)

Next, using (45) and the field equations, the condition $q=0$ reads (see [8] for details)

$${\displaystyle \frac{{\dot{R}}^{\prime }}{R}}-2{\displaystyle \frac{\dot{R}}{R}}{\displaystyle \frac{R^{\prime }}{R}}-{\displaystyle \frac{\dot{R}}{R}}\left({\displaystyle \frac{h^{\prime }}{h}}-{\displaystyle \frac{1}{r}}\right)=0,$$

(47)

whose solution is

$$\begin{array}{c}\hfill R={\displaystyle \frac{r}{h\left(r\right)\left[f\left(t\right)+g\left(r\right)\right]}},\end{array}$$

(48)

implying

$$\begin{array}{c}\hfill B={\displaystyle \frac{1}{h\left(r\right)\left[f\left(t\right)+g\left(r\right)\right]}},A={\displaystyle \frac{1}{f\left(t\right)+g\left(r\right)}},\end{array}$$

(49)

where g and f are two arbitrary functions of their argument.

Thus the metric is defined up to three arbitrary functions ($g\left(r\right),f\left(t\right),h\left(r\right)$).

The function $f\left(t\right)$ will be obtained from the junction conditions (22) and (23).

Indeed, evaluating the mass function at the boundary surface $\mathrm{\Sigma }$ we obtain from (22) and (48)

$${\dot{R}}_{\mathrm{\Sigma }}^2={\alpha }^2{R}_{\mathrm{\Sigma }}^4\left[ϵ-V\left(R_{\mathrm{\Sigma }}\right)\right],$$

(50)

where

$${\alpha }^2\equiv {\displaystyle \frac{h_{\mathrm{\Sigma }}^2}{r_{\mathrm{\Sigma }}^2}},ϵ\equiv {\left(g^{\prime }\right)}_{\mathrm{\Sigma }}^2{h}_{\mathrm{\Sigma }}^2,$$

(51)

and

$$V\left(R_{\mathrm{\Sigma }}\right)={\displaystyle \frac{2\sqrt{ϵ}}{R_{\mathrm{\Sigma }}}}(1-a_1)+{\displaystyle \frac{a_1}{R_{\mathrm{\Sigma }}^2}}(2-a_1)-{\displaystyle \frac{2M}{R_{\mathrm{\Sigma }}^3}},$$

(52)

with $a_1\equiv {\displaystyle \frac{h_{\mathrm{\Sigma }}^{\prime }{r}_{\mathrm{\Sigma }}}{h_{\mathrm{\Sigma }}}}$.

On the other hand, from (23), using (48) we obtain

$$\begin{array}{c}\hfill 2{\ddot{R}}_{\mathrm{\Sigma }}{R}_{\mathrm{\Sigma }}-{\dot{R}}_{\mathrm{\Sigma }}^2-3ϵ{\alpha }^2{R}_{\mathrm{\Sigma }}^4-4{\alpha }^2\sqrt{ϵ}{R}_{\mathrm{\Sigma }}^3(a_1-1)\\ \hfill -{\alpha }^2{R}_{\mathrm{\Sigma }}^2{a}_1(a_1-2)=0.\end{array}$$

(53)

Thus assuming $M=0$, Equation (50) becomes

$${\dot{R}}_{\mathrm{\Sigma }}^2={\alpha }^2{R}_{\mathrm{\Sigma }}^4\left[ϵ-{\displaystyle \frac{2\sqrt{ϵ}(1-a_1)}{R_{\mathrm{\Sigma }}}}-{\displaystyle \frac{a_1(2-a_1)}{R_{\mathrm{\Sigma }}^2}}\right].$$

(54)

Solutions to the above equation in terms of elementary functions may be obtained by assuming $a_1=1$, in which case a possible solution to (54) is

$$R_{\mathrm{\Sigma }}^{\left(II\right)}={\displaystyle \frac{1}{\sqrt{ϵ}cos\left[\alpha (t-t_0)\right]}},$$

(55)

which exhibits the same time dependence as in solution I.

Thus, as in the previous model, solution $II$ describes an expanding sphere, whose initial boundary areal radius grows from $1/\sqrt{ϵ}$ at $\alpha (t-t_0)=0$, to infinity as $\alpha (t-t_0)\to \pi /2$, and a contracting sphere whose boundary areal radius decreases from infinity at $\alpha (t-t_0)=3\pi /2$ to $1/\sqrt{ϵ}$ at $\alpha (t-t_0)=2\pi$. This picture repeating each time interval $\alpha (t-t_0)=2n\pi$, for any positive real integer n.

Imposing further the vanishing complexity factor condition, then functions $h\left(r\right),g\left(r\right)$ are given by

$$h\left(r\right)=c_{6}r,g=c_{7}lnr+c_5.$$

(56)

The physical variables corresponding to this solution read

$$\begin{array}{ccc}\hfill 8\pi \mu & =& -2ϵ{cos}^2\left[\alpha (t-t_0)\right]+\hfill \\ & +& ϵln{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\left[2cos\left[\alpha (t-t_0)\right]+ln{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right],\hfill \end{array}$$

(57)

$$\begin{array}{c}\hfill 8\pi P_r=-ϵln{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\left[4cos\left[\alpha (t-t_0)\right]+ln{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}\right],\end{array}$$

(58)

$$8\pi P_{\perp }=ϵ{cos}^2\left[\alpha (t-t_0)\right]-2ϵln{\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}cos\left[\alpha (t-t_0)\right].$$

(59)

The behavior of these physical variables is depicted in Figure 3.

From the definition of the mass function (11), using (48), (49), (55) and (56), the condition $m\left(r_{\mathrm{\Sigma }}\right)=0$ implies

$$ϵ{cos}^2\left[\alpha (t-t_0)\right]\left(1-{\displaystyle \frac{{\alpha }^2}{c_6^2}}\right)+{\displaystyle \frac{ϵ{\alpha }^2}{c_6^2}}-c_7^2{c}_6^2=0.$$

(60)

The above equation is satisfied for any value of t if ${\alpha }^2=c_6^2$ and $ϵ=c_7^2{\alpha }^2$, which, as expected, are the same relationships which follow from (51) and (56).

In oder to determine the regions of the fluid distribution where the energy-density is negative (required to have a vanishing total mass) we shall write the condition $\mu =0$ from (57) in the form

$$lnz\left(2cosy+lnz\right)-2{cos}^{2}y=0,$$

(61)

whose solution reads

$$z=e^{-(1+\sqrt{3})cosy},$$

(62)

where $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$, $y\equiv \alpha (t-t_0)$.

The graphic of z as function of y is plotted in Figure 4. The curve contains all the points of the plane $z,y$ where the energy-density vanishes, and divides the plane in two regions (A and B) corresponding to negative and positive energy-density respectively.

3.2. Expansion–Free Models

We shall now present models for which the expansion scalar vanishes.

We recall that under such a condition the line element may be written as (see [24] for details)

$$\begin{array}{c}\hfill d{s}^2=-{\left\{{\displaystyle \frac{R^2\dot{R}}{\alpha }}exp\left[-4\pi \int qAB{\displaystyle \frac{R}{\dot{R}}}dr\right]\right\}}^{2}d{t}^2\\ \hfill +{\displaystyle \frac{{\alpha }^2}{R^4}}d{r}^2+R^2(d{\theta }^2+{sin}^2\theta d{\varphi }^2),\end{array}$$

(63)

where $\alpha$ is a unit constant with dimensions $\left[r^2\right]$.

Figure 3. $8\pi \mu /ϵ$, $8\pi P_r/ϵ$ and $8\pi P_{\perp }/ϵ$, as functions of $y\equiv {\displaystyle \frac{\alpha }{2}}(t-t_0)$ and $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ for solution II.

Figure 3. $8\pi \mu /ϵ$, $8\pi P_r/ϵ$ and $8\pi P_{\perp }/ϵ$, as functions of $y\equiv {\displaystyle \frac{\alpha }{2}}(t-t_0)$ and $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ for solution II.

Figure 4. $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ as function of $y\equiv \alpha (t-t_0)$ for the condition $\mu =0$. Regions A and B correspond to negative and positive values of $\mu$ respectively, for solution II.

Figure 4. $z\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$ as function of $y\equiv \alpha (t-t_0)$ for the condition $\mu =0$. Regions A and B correspond to negative and positive values of $\mu$ respectively, for solution II.

Besides, as is known, the expansion-free models are endowed with an internal vacuum cavity surrounding the center, accordingly the center of symmetry is not filled with fluid.

The solutions deployed below satisfy the vanishing complexity factor condition and the quasi-homologous evolution.

3.2.1. Solution III: $Y_{TF}=0$, $A=A\left(r\right)$, $R=R_1\left(t\right)R_2\left(r\right)$

For this model, we shall complement the expansion–free condition with the vanishing complexity factor condition $Y_{TF}=0$, and we shall assume that A only depends on the radial coordinate, and R is a separable function (i.e., $R=R_1\left(t\right)R_2\left(r\right)$).

From all these conditions the general form of the metric variables read (see [24] for details)

$$\begin{array}{ccc}\hfill R_1\left(t\right)& =& {\displaystyle \frac{{\nu }_0}{t+{\nu }_1}},\hfill \end{array}$$

(64)

$$\begin{array}{ccc}\hfill R_2\left(r\right)& =& {\nu }_2{A}^{{\nu }_3-1},\hfill \end{array}$$

(65)

$$\begin{array}{ccc}\hfill A& =& {\nu }_4{\left({\nu }_{3}r+{\nu }_5\right)}^{{\displaystyle \frac{1}{{\nu }_3}}},\hfill \end{array}$$

(66)

where ${\nu }_0$ and ${\nu }_3$ are dimensionless constants and ${\nu }_1,{\nu }_2,{\nu }_4$ and ${\nu }_5$ are constants with dimensions $\left[r\right],\left[r^2\right],[1/\left(r^{1/{\nu }_3}\right)]$ and $\left[r\right]$, respectively.

Let us choose for our model

$${\nu }_0=1,{\nu }_1={\nu }_5=0;{\nu }_3=2;{\nu }_2={\displaystyle \frac{1}{{\nu }_4^4}}.$$

(67)

Then the expression for the mass function evaluated at the boundary surface and the areal radius of the boundary surface become

$$m\left(r_{\mathrm{\Sigma }}\right)={\displaystyle \frac{{\alpha }^{1/3}{r}_{\mathrm{\Sigma }}^{1/3}}{\sqrt{2}{t}^{*}}}\left({\displaystyle \frac{1}{t^{*4}}}+1-{\displaystyle \frac{2}{t^{*6}}}\right),$$

(68)

and

$$R\left(r_{\mathrm{\Sigma }}\right)={\displaystyle \frac{{\alpha }^{1/3}{r}_{\mathrm{\Sigma }}^{1/3}\sqrt{2}}{t^{*}}},$$

(69)

where $t^{*}={\nu }_4^{2}t$ and we have put ${\displaystyle \frac{r_{\mathrm{\Sigma }}}{{\alpha }^2{\nu }_4^6}}=1$.

From (68) it follows that the total mass vanishes at $t^{*}=1$.

The physical variables for this model read

$$8\pi \mu ={\displaystyle \frac{1}{{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}}}\left({\displaystyle \frac{t^{*2}}{2x}}-{\displaystyle \frac{3}{t^{*4}}}-{\displaystyle \frac{3}{2x{t}^{*2}}}\right),$$

(70)

$$4\pi q=-{\displaystyle \frac{1}{{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}}}\left({\displaystyle \frac{\sqrt{2}}{t^{*3}\sqrt{x}}}\right),$$

(71)

$$8\pi P_r=-{\displaystyle \frac{1}{{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}}}\left({\displaystyle \frac{5}{2x{t}^{*2}}}+{\displaystyle \frac{t^{*2}}{2x}}-{\displaystyle \frac{3}{t^{*4}}}\right),$$

(72)

$$8\pi P_{\perp }={\displaystyle \frac{1}{{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}}}\left({\displaystyle \frac{3}{t^{*4}}}-{\displaystyle \frac{1}{x{t}^{*2}}}\right),$$

(73)

$$T={\displaystyle \frac{{\alpha }^{1/3}}{\sqrt{2x}{r}_{\mathrm{\Sigma }}^{2/3}}}\left[T_0\left(t^{*}\right)+{\displaystyle \frac{3\tau }{2\pi \kappa t^{*2}\alpha \sqrt{2x}}}+{\displaystyle \frac{r_{\mathrm{\Sigma }}^{1/3}ln\left(2x{r}_{\mathrm{\Sigma }}\right)}{4\pi \kappa {\alpha }^{2/3}{t}^{*}}}\right],$$

(74)

where the expression above for the temperature has been obtained using the truncated transport Equation (26), and $x\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$.

The “effective energy–density” appearing in the definition of the mass function (15) $\left(\mu +q{\displaystyle \frac{U}{E}}\right)$ evaluated at $t^{*}=1$ reads

$$8\pi {\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}\left(\mu +{\displaystyle \frac{qU}{E}}\right)={\displaystyle \frac{1}{x}}-3.$$

(75)

Figure 5 and Figure 6 depict the behavior of physical variables, the radial distribution of the “effective energy–density” at $t^{*}=1$, when the total mass vanishes, and the evolution of the total mass.

The model represents a contracting sphere with initial negative mass absorbing energy through the boundary surface. At $t^{*}=1$ the total mass vanishes becoming positive afterward.

It is worth mentioning that although the total mass tends to zero as $t\to \infty$, the fluid distribution does not characterize a ghost star in that limit, since in such a case the total mass tends to zero due to the fact that the integrand in (15) tends to zero as $t\to \infty$, and not because of change of sign of the effective energy-density as is the case for a ghost star.

3.2.2. $A=1$, $Y_{TF}=0$

For this model we shall assume that the fluid is geodesic, meaning

$$A(t,r)=1,$$

(76)

the above condition together with the expansion-free condition

$$B(t,r)={\displaystyle \frac{\alpha }{R^2}},$$

(77)

plus the condition $Y_{TF}=0$, produces

$$R={\displaystyle \frac{\alpha }{\left[t+\alpha b_2\left(r\right)\right]}},B={\displaystyle \frac{{\left[t+\alpha b_2\left(r\right)\right]}^2}{\alpha }},$$

(78)

where $b_2$ is an arbitrary function of its argument with dimensions $[1/r]$.

Figure 6. $8\pi P_{\perp }{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}$ as function of $t^{*}\equiv {\displaystyle \frac{t{r}_{\mathrm{\Sigma }}^{1/3}}{{\alpha }^{2/3}}}$ and $x\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$; $2^{1/2}m\left(r\mathrm{\Sigma }\right)/{\alpha }^{1/3}{r}_{\mathrm{\Sigma }}^{1/3}$ as function of $t^{*}$; $8\pi (\mu +q{\displaystyle \frac{U}{E}}){\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}$, evaluated at $t^{*}=1$, as function of x, for solution III.

Figure 6. $8\pi P_{\perp }{\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}$ as function of $t^{*}\equiv {\displaystyle \frac{t{r}_{\mathrm{\Sigma }}^{1/3}}{{\alpha }^{2/3}}}$ and $x\equiv {\displaystyle \frac{r}{r_{\mathrm{\Sigma }}}}$; $2^{1/2}m\left(r\mathrm{\Sigma }\right)/{\alpha }^{1/3}{r}_{\mathrm{\Sigma }}^{1/3}$ as function of $t^{*}$; $8\pi (\mu +q{\displaystyle \frac{U}{E}}){\alpha }^{2/3}{r}_{\mathrm{\Sigma }}^{2/3}$, evaluated at $t^{*}=1$, as function of x, for solution III.

For our model we shall choose

$$b_2=-{\displaystyle \frac{\sqrt{6}r}{\alpha }},\alpha =6r_{\mathrm{\Sigma }}^2,$$

(79)

producing

$$R={\displaystyle \frac{6^{(1/2)}{r}_{\mathrm{\Sigma }}}{t^{*}-x}},B={(t^{*}-x)}^2,$$

(80)

where

$$t^{*}={\displaystyle \frac{t}{\sqrt{6}{r}_{\mathrm{\Sigma }}}},r=x{r}_{\mathrm{\Sigma }},$$

(81)

and the mass function becomes

$$m={\displaystyle \frac{3r_{\mathrm{\Sigma }}}{\sqrt{6}{(t^{*}-x)}^9}}\left[{(t^{*}-x)}^8+{(t^{*}-x)}^4-6\right].$$

(82)

From (82) evaluated at $r=r_{\mathrm{\Sigma }}$ we see that $m\left(r_{\mathrm{\Sigma }}\right)=0$ at $t^{*}={\left(2\right)}^{1/4}+1\approx 2.19$, being negative before that time and becoming positive afterward.

For this metric, the physical variables and the shear read as follows:

$$\begin{array}{c}\hfill 8\pi \mu =-{\displaystyle \frac{1}{2r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^2}}-{\displaystyle \frac{9}{r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^6}}+{\displaystyle \frac{{\left(t^{*}-x\right)}^2}{6r_{\mathrm{\Sigma }}^2}},\end{array}$$

(83)

$$\begin{array}{c}\hfill 4\pi q=-{\displaystyle \frac{2\sqrt{6}}{3r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^4}},\end{array}$$

(84)

$$8\pi P_r=-{\displaystyle \frac{5}{6r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^2}}+{\displaystyle \frac{1}{r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^6}}-{\displaystyle \frac{{\left(t^{*}-x\right)}^2}{6r_{\mathrm{\Sigma }}^2}},$$

(85)

$$8\pi P_{\perp }=-{\displaystyle \frac{1}{3r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^2}}+{\displaystyle \frac{4}{r_{\mathrm{\Sigma }}^2{\left(t^{*}-x\right)}^6}},$$

(86)

$$\sigma ={\displaystyle \frac{3}{\sqrt{6}{r}_{\mathrm{\Sigma }}\left(t^{*}-x\right)}},$$

(87)

$$T=T_0\left(t\right)+{\displaystyle \frac{1}{\pi \kappa \sqrt{6}{r}_{\mathrm{\Sigma }}(t^{*}-x)}}\left[1-{\displaystyle \frac{2\tau }{\sqrt{6}{r}_{\mathrm{\Sigma }}(t^{*}-x)}}\right],$$

(88)

where, as in the previous case, the temperature has been calculated using the truncated transport Equation (26).

Figure 7 and Figure 8 depict the behavior of physical variables, the radial distribution of the “effective energy–density” at $t^{*}={\left(2\right)}^{1/4}+1$, when the total mass vanishes, and the evolution of the total mass.

The model represents a contracting sphere with initial negative mass absorbing energy through the boundary surface. At $t^{*}={\left(2\right)}^{1/4}+1$ the total mass vanishes becoming positive afterward. From the above it follows that

$$8\pi r_{\mathrm{\Sigma }}^2\left(\mu +q{\displaystyle \frac{U}{E}}\right)={\displaystyle \frac{5}{6{\left(t^{*}-x\right)}^2}}-{\displaystyle \frac{9}{{\left(t^{*}-x\right)}^6}}+{\displaystyle \frac{{\left(t^{*}-x\right)}^2}{6}}.$$

(89)

Evaluating the above expression at $t^{*}={\left(2\right)}^{1/4}+1$ we see that it changes of sign within the fluid distribution, thereby explaining the vanishing of the total mass at $t^{*}={\left(2\right)}^{1/4}+1$ (see Figure 8).

4. Discussion

We have presented a set of solutions of fluid spheres whose evolution involves ghost stars.

The first two solutions represent the adiabatic evolution of a ghost star, they admit a CKV which may be either orthogonal or parallel to the four-velocity vector. These solutions (I and $II$) describe either an expanding sphere, with an initial boundary areal radius growing from some finite value to infinity, or a contracting sphere whose boundary areal radius decreases from infinity to some finite final value. In both cases $m\left(r_{\mathrm{\Sigma }}\right)=0$ at all times and the energy–density is negative in some regions of the fluid distributions (see Figure 2 and Figure 4). These models are described by Equations (38)–(42) and (55)–(60) for models I and $II$ respectively. Their behavior is depicted in Figure 1, Figure 2, Figure 3 and Figure 4.

In both cases the vanishing complexity factor condition applies and the solutions match smoothly to the Minkowski space-time on the boundary surface of the fluid distribution.

The second couple of solutions satisfies the expansion-free condition and the vanishing complexity factor condition. In one case (solution $III$), these last conditions are complemented with the assumption that $A=A\left(r\right)$ and R is a separable function. In the last case (solution $IV$) we assume the fluid to be geodesic.

Model $III$ is described by Equations (64)–(75) and illustrated in Figure 5 and Figure 6. It describes a collapsing fluid whose total mass evolves from negative values to positive ones by absorbing radiation. At some point of the evolution ($t^{*}=1$) the total mass vanishes becoming positive afterward. At $t^{*}=1$ the effective energy–density is negative in some regions of the fluid distributions as shown in Figure 6.

This rather unusual scenario (compact object absorbing radiation), has been considered before [34] in order to explain the presence of gas in quasars, and modeled in [35].

Finally, the last model is geodesic (solution IV) and satisfies, besides the expansion–free condition, the vanishing complexity factor condition. Equations (80)–(88) provide full description of such model which depicts a collapsing fluid for which as $t\to \infty$, the energy-density and the radial pressure diverge and satisfy the equation of the state $\mu =-P_r$, while the heat flux vector and the tangential pressure vanish, and the temperature tends to $T_0$. As it happens in Solution $III$, at some point of its evolution $t^{*}={\left(2\right)}^{1/4}+1$ the total mass vanishes, and as expected the effective energy–density is negative in some regions of the fluid distributions as shown in Figure 8.

Unlike solutions $I,II$, solutions $III,IV$ do not satisfy Darmois conditions on the boundary surface, implying the presence of thin shells on these surfaces.

It is worth recalling that the very existence of ghost stars relies on the presence of negative energy–density (the effective energy–density in the non–adiabatic case). Negative energy-density (mass) is a condition which has been shown to be likely to appear in scenarios where quantum effects are relevant (see [36,37,38,39,40,41] and references therein).

An issue requiring much more research work concerns the possibility to observe a ghost star. We have in mind either a “permanent” ghost star as the case described by solutions I and $II$, or a compact object attaining momentarily the ghost star status, as the two models described by solutions $III$ and $IV$.

At present we contemplate three possible ways to establish (or dismiss) the very existence of a ghost star. On the one hand, by observing the shadow of such objects following the line of research open by the Event Horizon Telescope (EHT) Collaboration (see [42,43,44,45] and references therein). On the other hand the appearance of a ghost star, even if for a short time interval, involves radiating processes whose observation could help to identify a ghost star. Finally a ghost star exhibits null gravitational surface redshift, a fact that could be used to its identification. Models as $III$ and $IV$ seem to be particularly suitable for this purpose.

In relationship with this last point it would be very helpful to find an evolving model, with a positive energy flux at the boundary surface, leading asymptotically to a ghost star. Unfortunately, neither of the solutions presented here satisfy such a condition.