A stellar model with diffusion
in general relativity
Abstract
We consider a spherically symmetric stellar model in general relativity whose interior consists of a pressureless fluid undergoing microscopic velocity diffusion in a cosmological scalar field. We show that the diffusion dynamics compel the interior to be spatially homogeneous, by which one can infer immediately that within our model, and in contrast to the diffusionfree case, no naked singularities can form in the gravitational collapse. We then study the problem of matching an exterior Bondi type metric to the surface of the star and find that the exterior can be chosen to be a modified Vaidya metric with variable cosmological constant. Finally, we study in detail the causal structure of an explicit, selfsimilar solution.
1 Introduction
The structure of singularities formed in the gravitational collapse of bounded matter distributions is a widely investigated problem in general relativity. The question of whether such singularities are naked, i.e., visible to faraway observers, or whether on the contrary they are safely hidden inside a black hole, has been the subject of innumerable works in the physical and mathematical literature. Nevertheless the problem remains poorly understood, even in the idealized case in which the collapsing body is spherically symmetric. A complete solution is only available for the gravitational collapse of a spherically symmetric dust cloud and can be found in the pioneering works by OppenheimerSnyder [1] and Christodoulou [2]. (In [3, 4] Christodoulou analyses the question of existence and stability of naked singularities in the gravitational collapse of a massless scalar field. Notwithstanding the importance of these works, our focus is on matter models that describe material bodies, such as perfect fluids and kinetic particles [5].)
The OppenheimerSnyder model consists of a collapsing spatially homogeneous and isotropic dust interior, described by the contracting FriedmannLemaître solution, matched at a comoving boundary with a Schwarzschild exterior. Dropping the homogeneity assumption of the interior leads to the class of LemaîtreTolmanBondi solutions. These inhomogeneous stellar models were studied numerically by Eardley and Smarr [6], and analytically by Christodoulou [2]. It was shown that, in contrast to the case studied by Oppenheimer and Snyder, the spatially inhomogeneous collapse leads to the formation of naked singularities. See the prologue of [7] for an historical review on the gravitational collapse problem in general relativity.
In this paper we initiate the study of the gravitational collapse of matter subject to diffusion. We believe that the inclusion of diffusion dynamics in the gravitational collapse problem is meaningful both from a mathematical and physical point of view. From one hand it is well known that diffusion terms introduce a regularizing effect in the equations, which might prevent the formation of naked singularities in general relativity. On the other hand the physical relevance of diffusion phenomena is unquestionable and the applications in general relativity have been discussed in [8, 9, 10].
We begin our study with the simplest possible model, namely a spherically symmetric dust cloud undergoing diffusion in a cosmological scalar field. In this case the regularizing effect due to diffusion is overwhelming: The interior of the dust cloud is forced to be spatially homogeneous. By this fact one can easily infer that, in contrast to the diffusionfree scenario described above, naked singularities cannot form in the gravitational collapse of a spherically symmetric dust cloud in the presence of diffusion.
Another interesting property of our model is that, in contrast to the diffusionfree case [1, 11, 2], the exterior of the star cannot be static. The simplest spherically symmetric solution of the Einstein equations that can provide a suitable exterior region for our stellar model is given by a Vaidya type metric which includes a variable cosmological constant (a generalization of the radiating version of the Schwarzschild(Anti)deSitter family of solutions).
A detailed analysis of our model is given in the following sections. We conclude this introduction by outlining the diffusion theory of matter in general relativity. There exist two versions of this theory: a kinetic one [8], which is based on a FokkerPlanck equation for the particle density in phasespace, and a fluid one [9], which is the (formal) macroscopic limit of the kinetic theory. In the present paper we apply the fluid theory. We recall that the energymomentum tensor and current density of a perfect fluid are
(1) 
where is the restframe energy density, the pressure, the fourvelocity and the particle density of the fluid. The diffusion behavior is imposed by postulating the equations
(2a)  
(2b) 
where is the diffusion constant, which measures the energy gained by the particles per unit of time due to the action of the diffusion forces. The second equation entails the conservation of the total number of fluid particles.
By projecting (2a) along the direction of and onto the hypersurface orthogonal to , we obtain the following equations on the matter fields:
(3a)  
(3b) 
The system (3) on the matter variables must be completed by assigning an equation of state between the pressure, the energy density and the particles number density. In this paper we assume that the fluid is pressureless (dust fluid). As the energymomentum tensor of the fluid is not divergencefree, see (2a), we have to postulate the existence of an additional matter field in spacetime to reestablish the (local) conservation of energy. The role of this additional matter field is that of the solvent matter in which the diffusion of the fluid particles takes place. The simplest choice is to assume the existence of a vacuum energy scalar field with energymomentum tensor , which leads to the following Einstein equations for the spacetime metric (in units ):
(4) 
The evolution equation on the scalar field determined by (4), the Bianchi identities, and the diffusion equation (2a) is
(5) 
2 The stellar model
Throughout the rest of the paper we assume that spacetime is spherically symmetric. The particle number density and the fourvelocity at each point are also spherically symmetric and satisfy (2b). Under these assumptions one can cover an open neighborhood of the center of symmetry by a coordinate system such that the metric and the fourvelocity take the form
Here are functions of and is the standard metric on . The center of symmetry is defined by the timelike curve . This coordinate system is called comoving and it is defined up to a transformation of the time and radial coordinates.
2.1 The interior
We assume that the interior of the star defines a region of spacetime covered by comoving coordinates. In particular, we assume that there exist and such that and within . The timelike hypersurface given by
will be identified with the surface of the star. To identify the region as the interior, we assume that the matter fields are nowhere vanishing on . The freedom in the choice of the comoving coordinates will be used to impose
(6) 
so that is the proper time of observers at rest with respect to the boundary of the star, and
(7) 
so that the comoving radius coincides initially, i.e., at time , with the radius function of the group orbits. It will now be proved that, when the fluid is pressureless, the diffusion dynamics compel the interior of our model to be spatially homogeneous. We denote , , for any function , and similarly denotes the derivative of any function of one variable.
Theorem 2.1.
Proof.
Equation (3b) for reads , which together with the boundary condition (6) implies . Hence , as claimed. Equation (2b) now reads
(9) 
Moreover (5) implies that and are functions of only, related by
(10) 
The Einstein equations (4) are:
(11)  
(12)  
(13)  
(14) 
where the only nonzero component of the energy momentum tensor is . The last Einstein equation gives
(15) 
where is an arbitrary positive function of the radial variable. Substituting into (9) we obtain
Integrating with the boundary condition (7) we obtain
(16) 
where we denoted . Substituting (15) and (16) into (12) we obtain the equation
The latter implies that there exists a constant such that
Finally (11) entails that is a function of only. ∎
From now on we assume that the interior fluid is pressureless. By rescaling the radial coordinate we may assume that , hence the spacetime metric in the interior of the star takes the RobertsonWalker form
(17) 
where for , and the equations (3)(5) reduce to the following. The conservation of the particle number density (2b) gives , while satisfy
(18a)  
(18b)  
(18c)  
(18d) 
where
(19) 
The initial data for the system (18) are given on and consist of a quadruple , with , such that the constraint equation (18a) is satisfied at time , i.e.,
In particular, the set of admissible initial data comprises a threedimensional manifold, which we denote by . By a regular solution of (18) in the interval with initial data we mean a triple of functions , satisfying (18b)(18d) and such that , for , and . Let be the maximal time of existence of a regular solution with a given set of initial data. We say that the regular solution is global if .
For , i.e., in the absence of diffusion, the interior reduces to a dust cloud in a spacetime with cosmological constant . In particular, for and , the exterior must be Schwarzschild, and we recover the wellknown OppenheimerSnyder model [1]. The analysis of the latter model is greatly simplified by the fact that a closed formula solution is known for the Einstein equations, see eg. [12]. In contrast to this, the general solution to the system (18) is not known, and so the analysis of the stellar interior in the presence of diffusion is more complicated. This analysis has been carried out in [13], where a complete characterization of the qualitative behavior of solutions to the system (18) has been obtained using dynamical systems methods. The relevant properties for the present study are summarized in the following theorem:
Theorem 2.2.
The manifold of initial data can be written as the disjoint union of two threedimensional submanifolds , such that

For initial data in the submanifold , the corresponding regular solution of (18) is global and holds for all . We call these interior models “expanding”.

For initial data in the submanifold , the corresponding regular solution of (18) exists only up to a finite time and as . Moreover there exists a constant such that
We call these interior models “collapsing”.
The asymptotic behavior of the scalar factor claimed in (ii) corresponds to the property proved in [13] that collapsing (dust) solutions behave like the FriedmannLemaître diffusionfree solution in the limit toward the singularity. It implies that the spacelike singularity at is a curvature singularity, as the Kretschmann scalar satisfies
In particular, spacetime is inextendible beyond the spacelike hypersurface and no outgoing light ray can emanate from the singularity. We obtain the following important corollary:
Corollary 2.1.
There exist no naked singularities in the spherical collapse of dust clouds undergoing diffusion in a cosmological scalar field.
We remark that the above result is independent of the spacetime exterior and applies to all collapsing models. In the absence of diffusion the formation of (locally and globally) naked singularities in the gravitational collapse of a dust cloud occurs for spatially inhomogeneous models, see [2]. In the present case the action of the diffusion forces prevents the formation of naked singularities by an explicit and compelling regularizing effect: it forces the dust interior to be spatially homogeneous.
Let us recall the definition of some important geometrical/physical quantities associated to the metric (17). Let
(20) 
The region in the interior where is the region of trapped surfaces, while defines the regular interior region. The hypersurface defines the apparent horizon. We define a local mass function through the identity
(21) 
which reduces to the standard definition of MisnerSharp mass when [11]. By (18a) we obtain the important identity
Moreover by (18c)(18d) we have
hence the mass function is given by
(22) 
Thus the local mass is conserved only in the absence of diffusion and it is otherwise linearly increasing in time. The latter behavior is intimately connected with the irreversibility of the diffusion process. In fact, using that the entropy of a dust fluid equals the energy per particle, i.e., (see [9]), and due to the conservation of the total number of particles, i.e., , we have
hence by (22) and recalling the definition (19) of the parameter we obtain that the entropy is linearly increasing in time:
2.2 The exterior
We write the metric on the exterior using Bondi coordinates:
(23) 
where are functions of and . The time coordinate is the ingoing (advanced) null coordinate for and the outgoing (retarded) null coordinate for . We assume that the comoving boundary between the interior and the exterior region is expressed as in Bondi coordinates, i.e.,
so that in (23). As the boundary is assumed to be timelike, we require that
(24) 
i.e., the first fundamental form of has the signature . Moreover, letting be the transformation of the time coordinate on the boundary, we require , which means that the time orientation of spacetime does not change across the boundary. We recall that two metrics and may be matched on if and only if they satisfy the junction conditions that they induce the same first and second fundamental form on , see [11] and references therein.
Theorem 2.3.
Proof.
The first and second fundamental forms induced by the interior metric (17) on are given respectively by
while the same quantities induced by the exterior metric (23) are
Hence the junction condition gives immediately (25) as well as
(29) 
The junction condition is equivalent to and
Using (25) in the latter equation we obtain
In the term we use
and so doing we obtain (27). Replacing (25) and (27) into (29) gives (26). ∎
In principle, any exterior metric that satisfies the conditions in the theorem can be matched to the interior. However a particularly simple and natural choice can be made as follows. We assume that in the whole exterior, so that (27) gives
(30) 
Using (20) and (21) in (30) we obtain
(31) 
Moreover the junction condition (d) in Theorem 2.3 becomes
Using (18b) and (31) we may write the latter as
(32) 
If we now require (32) to be valid for all , and not only on , and solve the resulting differential equation on subject to the boundary condition (31) at we obtain
where
(33) 
Hence we obtain the following corollary.
Corollary 2.2.
The exterior metric can be chosen to be the Vaidya type metric with variable cosmological constant given by
(34) 
where is given by (33) and
(35) 
is the exterior mass function.
Conversely, if we now start with a metric of the form (34) and impose the matching conditions (31)(32), we obtain the identities and , that is to say, the mass function and the cosmological scalar field must be continuous through the boundary in order for the metric (34) to be an admissible exterior.
The generalized Vaidya metric (34) solves the Einstein equation (4) with cosmological scalar field and the energymomentum tensor
where, by (18c), (22) and (26),
(36)  
(37) 
When , i.e., in the absence of diffusion, the functions and become constant and so . In this limit our model reduces to the OppenheimerSnyder model with cosmological constant [14]. The energymomentum tensor describes a cloud of dust particles moving along the null directions orthogonal to . Observe that the dust particles in the exterior do not undergo diffusion in the cosmological scalar field. This outcome of the model is consistent with the wellknown fact that particles moving along null directions are not subject to diffusion, see [15]. Requiring the weak energy condition to hold in the whole exterior region forces us to restrict to the outgoing Vaidya solution.
Proposition 2.1.
. Moreover the weak energy condition , for all , holds only for the outgoing Vaidya metric (34).
We remark that the property implies that the surface of the star is not radiating, i.e., the star is thermally isolated from the exterior. In particular the null dust particles in the exterior are not emanated by the star, but rather they are created spontaneously by the decaying vacuum energy field in the exterior (see [9] for a similar behavior occurring in cosmology). Note also that for we have for . Hence if we consider a spatially homogeneous exterior, the interior can be chosen to be the ingoing Vaidya metric, leading to a void model with diffusion which generalises the wellknown EinsteinStrauss model [16].
It is worth to point out the following important differences between our model and the analogous diffusionfree model studied in [17, 18]. Firstly, in the diffusionfree model the Vaidya solution can arise as the exterior of a fluid ball only if the pressure of the fluid is positive (dust is not allowed). Secondly the matching boundary in the diffusionfree model cannot be comoving. The latter property can be naturally regarded as being the cause of energy dissipation from the surface of the star as seen by comoving observers, see the discussion in [17, 18].
In view of Proposition 2.1 we restrict from now on to the generalized outgoing Vaidya metric (34), which we rewrite in standard notation as
(39) 
The apparent horizons of the metric (39) are the hypersurfaces in the exterior where
As , see (27), and owing to (20), the interior and exterior apparent horizons intersect on the boundary.
3 Selfsimilar solutions
In this section we study in detail the causal structure of a particular solution to our model. Recall that the interior is assumed to be pressureless. Assuming in addition that it is selfsimilar, we obtain the following solution of (4):
(40a)  
(40b)  
(40c)  
where is the real (positive) solution of the polynomial equation  
(40d) 
The Penrose diagrams of this solution are given in Figure 1. It has been shown in [13] that, for an open set of initial data, solutions of both types and have an intermediate behavior that is close to the special solution (40).
Note that we shifted the origin of time so that the singularity appears at . This is a curvature singularity as the Ricci scalar curvature satisfies
Toward the future the solution is forever expanding () without acceleration (). There is only one apparent horizon for this solution, located at
(41) 
Let the boundary of the star. We distinguish three cases:

; the interior has an apparent horizon in this case.

; the boundary of the star coincides with the apparent horizon.

; the interior has no apparent horizon.
Now, let the matching surface as seen by exterior observers. We assume that the exterior metric is given by (39). From (25), (26) we obtain
and (39) becomes
(42) 
where
The mass and the cosmological scalar field in the exterior are given by
Note also that the metric (42) has a curvature singularity at , for its Kretschmann scalar is given by
The apparent horizons of the metric (42) are the hypersurfaces where , where
(43) 
Proposition 3.1.
The following holds:

In case (i) there is no apparent horizon in the exterior region and for all .

In case (ii) the apparent horizon of the metric (42) coincides with the apparent horizon of the interior, as well as with the matching surface:
Moreover for all .

In case (iii) there exists such that the metric (34) has an apparent horizon at . Moreover for and for .
Proof.
For the proof it suffices to notice that the function in (43) attains its maximum at and . ∎
The model under discussion is selfsimilar, with the lines , , on the plane being tangent to the homothetic vector field in the exterior. We call such curves homothetic curves.
Theorem 3.1.
There exists such that the homothetic curve is spacelike for , null for and timelike for . In case (iii) there holds . Moreover the homothetic curve is the first ingoing radial null geodesics that escapes to nullinfinity.
Proof.
The metric induced on the hypersurface in the exterior is given by
where is the function defined in (43). The hypersurface is spacelike if , timelike if and null if . A straightforward calculation shows that
As when , there exists one (and only one) such that for , and for . The first part of the result follows. Moreover and so when the exterior apparent horizon exists we have . As to the last statement, the equation for the ingoing radial null geodesics is
Consider a solution such that , for some . Let . As
then is decreasing for . It follows that , for all . Hence , for all and therefore
We conclude that for large enough the ingoing light ray must hit the boundary in cases (i)(ii) and the apparent horizon in case (iii). In the latter case the ingoing light ray will enter the regular region and therefore hit the surface of the star in finite time. In conclusion, all ingoing light rays escaping to infinity other than the null homothetic curve must satisfy , for all , which concludes the proof of the theorem. ∎
The result of the previous theorem leads to the Penrose diagram shown in Figure 2.
Acknowledgements
A. A. is supported by the projects EXCL/MATGEO/0222/2012, PTDC/MATANA/1275/2014, by CAMGSD, Instituto Superior Técnico, by FCT/Portugal through UID/MAT/04459/2013, and by the FCT Grant No. SFRH/BPD/85194/2012.
The authors would like to thank J. Natário for numerous discussions and the anonymous referee for valuable comments that helped to improve the paper. Furthermore, A. A. thanks the Department of Mathematics at Chalmers University, Sweden, for the very kind hospitality.
References
 [1] J. Oppenheimer, Snyder : On Continued Gravitational Contraction. Physical Review 56, 455–459 (1939)
 [2] D. Christodoulou: Violation of Cosmic Censorship in the Gravitational Collapse of a Dust Cloud . Commun. Math. Phys. 93, 171–195 (1984)
 [3] D. Christodoulou: Examples of naked singularity formation in the gravitational collapse of a scalar field, Ann. of Math. 140, 607–653, (1994)
 [4] D. Christodoulou: The instability of naked singularities in the gravitational collapse of a scalar field. Ann. of Math. 149, 183217 (1999)
 [5] H. Andréasson: The EinsteinVlasov System/Kinetic Theory. Living Rev. Relativ. (2011) 14: 4
 [6] D. M. Eardley, L. Smarr. Time functions in numerical relativity: marginally bound dust collapse, Phys. Rev. D 19, 2239 (1979)
 [7] D. Christodoulou: The Formation of Black Holes in General Relativity. European Mathematical Society Monographs in Mathematics, (2009)
 [8] S. Calogero: A kinetic theory of diffusion in general relativity with cosmological scalar field. JCAP 11/2011, 016 (2011)
 [9] S. Calogero: Cosmological models with fluid matter undergoing velocity diffusion. J. Geom. Phys. 62, 2208–2213 (2012)
 [10] S. Calogero, H. Vetten: Cosmology with matter diffusion. JCAP 11/2013, 025 (2013)
 [11] F. Fayos, J. M. M. Senovilla, and R. Torres: General matching of two spherically symmetric spacetimes. Phys. Rev. D 54, 4862 (1996)
 [12] J. Griffiths and J. Podolsky. Exact spacetimes in Einstein’s General Relativity. Cambridge Monographs on Mathematical Physics, (2009).
 [13] A. Alho, S. Calogero, M. P. Ramos, A. J. Soares: Dynamics of RobertsonWalker spacetimes with diffusion. Ann. Physics 354, 475–488 (2015)
 [14] K. Nakao: The OppenheimerSnyder spacetime with a cosmological constant. Gen. Rel. Grav. 24, 10691081 (1992)
 [15] R. M. Dudley: Lorentzinvariant Markov processes in relativistic phase space. Ark. Mat. 6, 241268 (1966)
 [16] A. Einstein, E. G. Strauss, The influence of the expansion of space on the gravitation fields of surrounding the individual stars. Rev. Mod. Phys. 17, 120124 (1945)
 [17] F. Fayos, X. Jaen, E. Llanta, J. M. M. Senovilla: Matching the Vaidya and RobertsonWalker metric. Class. Quant. Gra. 8, 20572068 (1991)
 [18] W. B. Bonnor, A. K. G. de Oliveira, N. O. Santos: Radiating Spherical Collapse. Physics Reports 181 5, 269–326 (1989)