Gravitational Collapse and Spacetime Singularities
Gravitational Collapse and Spacetime Singularities
†
Issued as a paperback
Gravitational Collapse and
Spacetime Singularities
PA N K A J S . J O S H I
Tata Institute of Fundamental Research,
Mumbai, India
CAMBRIDGE UNIVERSITY PRESS
Cambridge, New York, Melbourne, Madrid, Cape Town, Singapore, São Paulo
© P. S. Joshi 2007
Cambridge University Press has no responsibility for the persistence or accuracy of urls
for external or third-party internet websites referred to in this publication, and does not
guarantee that any content on such websites is, or will remain, accurate or appropriate.
To my parents,
Arunadevi Shantilal Joshi
and
Shantilal Ramshankar Joshi
Contents
Preface page ix
1 Introduction 1
3 Spherical collapse 60
3.1 Basic framework 62
3.2 Regularity conditions 69
3.3 Collapsing matter clouds 71
3.4 Nature of singularities 79
3.5 Exterior geometry 87
3.6 Dust collapse 90
3.7 Equation of state 129
vii
viii Contents
References 255
Index 269
Preface
ix
x Preface
for numerous discussions and work as cited, on the themes described here.
Special thanks are due to R. Goswami and I. H. Dwivedi for their ideas and
help and for our studies together. A. Mahajan and S. Khedekar helped with
the manuscript.
1
Introduction
1
2 Introduction
when the quantum effects are incorporated together with the gravitational
force, the classical spacetime singularity may be resolved, and may no longer
exist in the full theory.
Gravitational collapse is thus a key phenomena for many astrophysical
processes for stars or other larger systems in the universe. In particular, the
very advanced stages of collapse of a massive star are occurrences in nature
where the effects of both gravity and the quantum would be combined. Even
if the final spacetime singularity, as predicted by classical general relativity,
may be resolved, possibly through quantum gravity effects, such a collapse
will necessarily give rise to spacetime regions of ultra-high mass densities and
curvatures, where the physical effects will be extreme.
The important physical issue would then be whether such extreme gravity
regions formed in the gravitational collapse of a massive star are visible to
external observers in the universe. An affirmative answer here would mean
that the physical phenomena of the gravitational collapse of a massive star
could provide a very good laboratory to study quantum gravity effects in the
cosmos, and this may help towards generating clues for an, as yet, unknown
theory of quantum gravity. A laboratory similar to that provided by the
early universe is then created in the later stages of the continual collapse of a
massive star. An additional feature would be that, whereas the early universe
was a unique event that happened only once, the collapse phenomena would
continue to occur whenever a sufficiently massive star in the universe died
on exhausting its nuclear fuel. If such ultra-strong gravity regions become
visible to external observers in the spacetime, an opportunity to observe the
quantum gravity effects in the universe is provided.
The answer to this is determined by the causal structure of spacetime
in the vicinity of a spacetime singularity. This is actually decided by the
dynamics of the gravitational collapse of the matter cloud, as it evolves from
a regular initial data, defined on an initial surface, from which the collapse
develops. This dynamical evolution is governed by the Einstein equations.
In other words, it is only the study of the collapse dynamics of the matter
clouds that would decide the visibility or otherwise of the ultra-strong gravity
regions. If, as the collapse evolves, the event horizons of gravity develop much
before the spacetime singularity forms, then these extreme gravity regions are
hidden away from the external universe, and a blackhole forms as the collapse
outcome. On the other hand, if such horizons are delayed or fail to develop
during collapse, as governed by the internal dynamics of the collapsing cloud,
then the scenario where the extreme gravity regions are visible to external
observers occurs, and a visible naked singularity forms.
The importance of gravitational collapse processes in relativistic astro-
physics was realized when Datt (1938) and Oppenheimer and Snyder (1939)
used general relativity to study the dynamical collapse of a homogeneous
spherical dust cloud under its own gravity. This model gave rise to the
4 Introduction
concept of a blackhole. The term blackhole itself was popularized only later
in the 1960s. The above work established that, under idealized conditions, a
collapsing cloud of matter with zero pressure will necessarily give rise to a
blackhole. Such a blackhole is a region of spacetime from which no light or
matter can escape away to faraway external observers, and which necessarily
covers the spacetime singularity or the regions of extreme physical condi-
tions from the external universe. Specifically, in order to create a blackhole
as the final state of gravitational collapse of the star, an event horizon must
develop in the spacetime earlier than the time when the final spacetime sin-
gularity forms. Such an event horizon is a one-way membrane such that light
or matter can fall into the region covered by it, but cannot escape away. If
the event horizon developed prior to the formation of the singularity, neither
the singularity nor the collapsing matter that has fallen within it would be
observable to an external observer, and a blackhole is said to have formed
as the final endstate of the collapsing star. All the matter of the star is then
supposed to be crushed into the infinite density singularity at the center of
the blackhole.
How early and when the horizon will actually develop in a realistic col-
lapse is determined by the dynamics of the collapsing matter, the physical
conditions within the star, and the dynamical evolution of the cloud as gov-
erned by the Einstein equations of gravity. Investigations in high energy
astrophysics have already used the concept of a blackhole quite extensively.
However, the actual understanding of the phenomena of gravitational col-
lapse, and the conditions under which it can lead to the blackhole formation,
or otherwise, within the framework of general relativity has progressed only
relatively recently.
Further to the early studies mentioned above, it was generally assumed
that the final endstate of collapse of a massive star will be a blackhole only.
However, several important questions remained unanswered. For example,
what would be the effects of non-zero pressures, which would be certainly
important in the later stages of collapse, towards determining the collapse
endstate, or, how will an inhomogeneous cloud collapse, say with a physically
realistic density profile that is higher at the center and decreases slowly as
one moves away from the center of the star? Early work on gravitational
collapse focused only on simple models with idealized conditions, assuming
a totally homogeneous density within the star, zero pressures, and so on,
which would not be physically realistic. For example, a realistic star must
have non-zero internal pressures, and its density would be typically higher
at the center, as compared with its outer layers.
These physical issues and important questions have been crucial to the
foundations of blackhole physics. But, not much attention could be paid to
them, mainly due to the complexity of the equations of general relativity.
This is because, in general, the Einstein equations are non-linear, second
Introduction 5
order partial differential equations that are quite difficult to solve. Therefore,
the only model available until the late 1960s for the dynamical gravitational
collapse of a massive matter cloud was that of a homogeneous, pressureless
spherical cloud. In addition, not much attention was paid to these issues
by the general relativists of the 1940s and 50s, who, by and large, did not
consider such ultra-high energy phenomena to be physically realistic or of
much astrophysical significance.
As indicated above, it was only the discovery in the 1960s of very
high energy astrophysical phenomena that generated a keen theoretical inter-
est in the continual gravitational collapse processes. However, mathematical
difficulties and the complexity of gravity theory did not allow much progress.
Then, the cosmic censorship hypothesis was introduced by Penrose (1969),
which conjectured that the outcome of any generic gravitational collapse of
a massive star must lead necessarily only to a blackhole formation as the
collapse final state. This hypothesis thus suggested that the extreme and
ultra-strong gravity regions, or the spacetime singularity, must always nec-
essarily be covered within an event horizon of gravity, and that the external
observers should never be able to see the singularity. This assumption means
that whatever the physical conditions and forces within the massive stars
may be (for example, they may be inhomogeneous in their density distribu-
tion, the pressures may be non-zero, or they may not be totally spherical
and so on), the outcome of their continual collapse must give rise to a black-
hole only. In other words, this amounts to an assumption of the nature of
the allowed dynamical evolutions of the collapsing clouds, namely that the
Einstein equations must permit only those evolutions that create the event
horizon necessarily much prior to the formation of the final singularity or
the ultra-strong gravity regions. Then, the singularity would be necessarily
hidden within the horizon, which is a one-way surface, not allowing it to
be seen by any external observers.
The cosmic censorship conjecture thus implies that no ultra-strong grav-
ity regions forming in continual collapse will be visible to outside observers.
That is, no naked singularity will develop in the collapse, and the event hori-
zon developing in the dynamical collapse will always manage to cover these.
Hence, the outcome of any gravitational collapse is necessarily a blackhole,
and external observers can never see any ultra-strong gravity regions forming
in the collapse, as indicated in Fig. 1.1.
As yet, a specific mathematical formulation for cosmic censorship that has
been properly defined does not exist. Then, a proof of the same would have to
be obtained within the framework of Einstein’s gravity theory. The cosmic
censorship assumption nevertheless provided a major impetus to develop-
ments in blackhole physics, and two parallel streams of developments took
place. On one hand, the theoretical properties of blackholes were devel-
oped extensively, using cosmic censorship as the basic assumption, thus
6 Introduction
Spacetime singularity
Event horizon
Initial surface
creating the laws of blackhole thermodynamics and related aspects (see for
example, Hawking and Ellis, 1973). On the other hand, efforts to estab-
lish the censorship hypothesis continued, as it was clear all along that this
assumption was absolutely fundamental to the theory and applications in
blackhole physics, and so it needed a rigorous formulation and proof within
the framework of general relativity. It is widely recognized that a proof
of the censorship conjecture would place blackhole physics and its applica-
tions on a sound footing, whereas its failure would actually throw blackhole
dynamics and related applications into serious doubt. Hence, the validity,
or otherwise, of the cosmic censorship conjecture has remained an issue of
crucial importance for all these years. The efforts to prove it have not suc-
ceeded for the past three decades, and there are even serious difficulties
in formulating any rigorous mathematical version or a statement for this
conjecture.
The theme that the only way out of this impasse is to study rigorously
the dynamical gravitational collapse phenomena within the framework of
Einstein’s theory of gravity is proposed and developed here. This has been
investigated extensively in the last couple of decades, and some of the
issues that have been addressed include: what is the outcome of a continual
gravitational collapse under physically realistic conditions, as governed by
the Einstein equations? Will it be necessarily a blackhole as hypothesized
by the censorship conjecture, or would it give rise to a naked singularity,
where ultra-strong gravity regions forming in collapse are visible to external
observers? In the latter case, would it be possible to observe the quantum
Introduction 7
gravity effects taking place in these visible ultra-strong gravity regions? Some
of these issues are discussed here.
A detailed study of the collapse phenomena may be the only way towards
any possible physically realistic formulation of censorship, if one exists. Such
a study and investigation of collapse could also lead to novel physical insights
and possibilities emerging out of the intricacies of the gravitational force. It
would appear that beyond the studies so far, mainly of static and stationary
solutions modeling blackholes, investigating dynamical evolutions as permit-
ted by the Einstein equations would offer new insights into the nature of
gravity. This is an arena that has been explored less, and which needs to be
investigated carefully in detail.
To this end, gravitational collapse scenarios with non-zero pressures and
more realistic equations of state for classes of general matter fields are con-
sidered here. A general formalism is developed to treat the spherical collapse
from regular initial data. These considerations also point to why it has not
been possible so far to make any definite progress on the censorship conjec-
ture. It is seen that it is first necessary to acquire a deeper and more extensive
understanding of the dynamical evolutions and gravitational collapse pro-
cesses in general relativity. Recent work on studying and understanding
the final fate of dynamical gravitational collapse in gravitation theory is
discussed. General matter fields are considered so as to include important
physical features in the collapse, such as inhomogeneities in matter distri-
bution, non-vanishing pressures, different forms for the equations of state
of the collapsing matter, and other such aspects. It is seen that in spherical
gravitational collapse, given the matter initial data on an initial surface from
which the collapse develops, there are the rest of the free initial data such
as the velocities of the collapsing shells, and the classes of the dynamical
evolutions as permitted by the Einstein equations, which lead to the final
state that is either a blackhole, or a naked singularity that is a visible ultra-
strong density and curvature region forming in the collapse not covered by
an event horizon. The nature of the outcome depends on the regular initial
data from which the collapse evolves, and the allowed dynamical evolutions
in the spacetime, as permitted by the Einstein equations.
After the basics of the structure and properties of spacetimes and the
essentials of relativity theory are summarized in Chapter 2, the above issues
are discussed in Chapter 3. Collapsing dust clouds, which generalize and
include as a special case the Oppenheimer–Snyder dust collapse models, and
which give an idea of the possible outcomes of gravitational collapse in terms
of a blackhole or a naked singularity, are also discussed in Chapter 3. The
Oppenheimer–Snyder dust collapse scenario is included here as a special case
when the cloud is homogeneous. It is seen, however, that a more realistic den-
sity profile with a density higher at the center and decreasing as one moves
away from center, gives rise to a naked singularity as the collapse endstate
8 Introduction
Spacetime singularity
Event horizon
Outgoing null geodesics
Initial surface
Fig. 1.2 If the collapsing cloud is inhomogeneous, with a density higher at the
center, the trapped surface formation and event horizon in the collapse are
delayed to give rise to a naked singularity, where the ultra-strong gravity regions
are visible to outside observers.
(see Fig. 1.2). In general, it is seen that the collapse outcome depends on
the nature of the initial matter profiles and the evolutions allowed by the
Einstein equations. The structure of this spacetime in the homogeneous den-
sity case gives rise to the basic notion and concept of a blackhole. The dust
collapse picture provides a concrete background to the possible final states
of a continual gravitational collapse.
Chapter 4 then studies several useful aspects of spacetime structure,
singularities and collapse, as related to the cosmic censorship hypothesis pos-
sibilities and the structure of naked singularities developing in gravitational
collapse. It is pointed out that while the cosmic censorship does not hold in
general relativity in the obvious sense of ruling out naked singularities from
all physically realistic gravitational collapse models, any definite formulation
of this hypothesis will depend on a detailed analysis of stability and gener-
icity aspects related to collapse scenarios, and the naked singularities and
blackhole phases developing as final outcomes of the gravitational collapse.
Several possibilities towards any plausible formulation are discussed.
In light of the results available so far and the emerging scenario, the key
physical issue is the possible final state of a massive star. The basic problem
to be addressed is: what will the final outcome of the gravitational collapse
of a massive star be when it collapses freely at the end of its life cycle on
exhausting its nuclear fuel under the force of its own gravity? Under realistic
astrophysical conditions, will it turn into a blackhole, or does it terminate as
a naked singularity? Are there any observable consequences in the latter case?
Introduction 9
10
2.1 The manifold model 11
(1) The sets uα form a cover for M , that is, any given p in M must be in a
uα for some value of α, and
M= uα . (2.1)
α
Rn
f(u)
M
u f f( p)
p
u f
q f(u)
f(q)
Fig. 2.1 All events p and q in the manifold have neighborhoods which are
homeomorphic to subsets in Rn . The points p, q ∈ M have coordinates
of φα (p) and φβ (q). Whenever the neighborhoods in M intersect, there should
be a smooth change of coordinates.
change, the change is necessarily smooth. That is, if {xi } and {y i } are
local coordinates of p ∈ M in uα and uβ respectively, then the functions
xi = xi ( y 1 , . . . , y n ) are C ∞ -functions from Rn to Rn . A maximal or complete
atlas is chosen for the spacetime manifold M, that is, if {uα , φα } is an atlas
for M , one selects for M the atlas that consists of all other atlases that are
compatible with {uα , φα }. This implies that their union with {uα , φα } is also
a C ∞ -atlas.
This implies that one has included all possible, mutually compatible coor-
dinate systems for the manifold M. A C r -manifold is defined in a similar way,
where it is required that the transition functions φα ◦ φ−1 β are r-times
continuously differentiable, where a continuous function is denoted by C 0 .
The Euclidian plane R2 , or Euclidian space Rn , is, in itself, a manifold as
it is covered by a single chart Rn , where φ would be the identity map with
the coordinate range −∞ < xi < ∞ for i = 1, . . . , n. Another example of
such a manifold is the two-sphere S 2 defined by
The six hemispherical open sets Oi± for i = 1, 2, 3 are given by Oi± =
{(x1 , x2 , x3 ) ∈ S 2 | ±xi > 0}, which cover S 2 . Each Oi± is mapped
onto the open disk {(x, y ) ∈ R2 | x2 + y 2 < 1} by the projection maps such
as f1+ (x1 , x2 , x3 ) = (x2 , x3 ). The overlap functions fi± ◦ ( fj± )−1 are C ∞ -
functions in their domain of definition. Thus, S 2 is a two-dimensional,
C ∞ -manifold that cannot be covered by a single coordinate system. Similarly,
the sphere S n in n-dimensions is also a differentiable manifold.
2.1 The manifold model 13
Fig. 2.2 The tangent space Tp at a point p ∈ M , which gives the set of all
directions at that point.
with at least one V i being non-zero, and applying this to the coordinate
functions x1 , . . . , xn gives V i = 0 for all i, a contradiction). Therefore, the
vectors (∂/∂xj ) span the vector space Tp , the space of all tangent vectors at
p (see Fig. 2.2). The vector space structure here is defined by
V = V i ei , (2.8)
where the numbers V i are called the components of V with respect to the
basis ei . In a coordinate basis, V i = dxi /dt. Again, {∂/∂xi } forms a basis
of Tp which means the dimension of Tp is n.
For the tangent space Tp at p ∈ M , the vector space of all dual vectors at
p, also called covariant vectors or one-forms at p, can be naturally defined. A
one-form ω at p is a real-valued linear functional on Tp , denoted by ω (X ) ≡
ω, X , and the linearity condition implies
Given a tangent space basis {ea }, a unique set of one-forms {ea } is given
by the condition that the one-form eb maps a vector V into V b , that is, the
bth component of V in the basis ea . Therefore,
eb , V = V b , (2.10)
where the right-hand side is the Kronecker delta function. The linear
combinations of one-forms ω and η are defined by
ωa = ω, ea . (2.13)
Therefore, the set of all one-forms at the event p forms a vector space at p,
the dual of Tp , and is denoted by Tp∗ . The basis ea is a dual basis to ea . If
ω ∈ Tp∗ and V ∈ Tp , then
The commutator for any two coordinate basis vector fields vanishes. If f
and g are any two smooth functions, it can be seen that [V , W ]( f + g ) =
[V , W ]( f )+[V , W ](g ) and that [V , W ](αf ) = α[V , W ]( f ) for any α ∈ R.
It can be shown that
[V , W ]( f g ) = f [V , W ](g ) + g [V , W ]( f ), (2.16)
[V 1 + V 2 , W ] = [V 1 , W ] + [V 2 , W ]. (2.17)
16 The spacetime manifold
∂f
df, V = V a . (2.19)
∂xa
The local coordinate functions (x1 , . . . , xn ) are used to define a set of one-
forms (dx1 , . . . , dxn ), which is a basis dual to the coordinate basis because
a ∂ ∂xa
dx , b = = δba . (2.20)
∂x ∂xb
Also,
∂ ∂f
df = df, dxa = dxa , (2.21)
∂xa ∂xa
Let B be the set of all ordered basis {ei } for Tp . If {ei } and {ej } are in
B , then
ej = aij ei . (2.24)
where {xi } and {xj } are local coordinates in Ui and Uj respectively. The
Möbius strip is an example of a non-orientable manifold.
2.1.4 Tensors
Tensors are geometric objects on a spacetime, and are invariant under the
change of coordinates. The stress–energy tensor represents various matter
fields existing on a spacetime, such as the electromagnetic field, dust, and so
on. On the other hand, the global geometry and curvature of the manifold
are described by tensor fields such as the metric tensor and the curvature
tensor. In the general theory of relativity, the form of physical laws remains
unchanged under the general transformation of coordinates (principle of gen-
eral covariance). So, physical fields are represented by various tensor fields on
a spacetime, and the laws governing these are tensor equations, valid under
arbitrary coordinate transformations. In an inertial coordinate frame, these
reduce to special relativity laws.
A tensor T of type (r, s) at p ∈ M is a multilinear real-valued functional
on the Cartesian product, as given by
where
T i1 ...ir j1 ...js ≡ T (ei1 , . . . , eir , ej1 , . . . , ejs ), (2.29)
and {ei } and {ei } are basis vectors at p for Tp and Tp∗ respectively.
The space of all tensors of type (r, s) at p is called the tensor product
Tsr ( p), denoted by,
with r-factors of Tp and s-factors of Tp∗ . The dimension of Tsr is nr+s , with n
being the dimension of the manifold. This is a vector space of all (r, s) tensors
over real numbers with the addition of tensors and scalar multiplication
defined in a natural manner. In particular, a vector is a tensor of type (1,0)
and a one-form is a tensor of type (0,1). Using the basis vectors {ei } and
{ei } for the tangent space and cotangent space at p, the set
forms a basis for the tensor product Tsr ( p) with all indices running from 1
to n. Then, any tensor T ∈ Tsr can be expressed as
Using the equations above for a transformation under the change of basis
vectors, it is seen that the contraction C 1 1 is independent of the basis used
and so invariant under a change of coordinates. Similarly, T can be contracted
over any pair of contravariant and covariant indices.
In the space of all tensors of type (r, s) at p, the addition of two tensors
T and T , and multiplication by a real number α, are defined as
(T + T )(ω 1 , . . . , ω r , X 1 , . . . , X s ) = T (ω 1 , . . . , ω r , X 1 , . . . , X s )
+ T (ω 1 , . . . , ω r , X 1 , . . . , X s ),
(αT )(ω 1 , . . . , ω r , X 1 , . . . , X s ) = αT (ω 1 , . . . , ω r , X 1 , . . . , X s ). (2.37)
The outer product of the two tensors T and S of type (r, s) and (r, s ) can
now be defined in terms of their components to give a new tensor T ⊗ S ,
(T ⊗ S )i1 ...ir+r j1 ...js+s = T i1 ...ir j1 ...js S ir+1 ...ir+r js+1 ...js+s . (2.38)
This allows new tensors to be constructed out of vectors and dual vectors.
A tensor field of type (r, s) on M is an assignment of a tensor of the same
type at all p in M . It is C k -differentiable if all the components of T have the
same differentiability as functions of the coordinates.
If T is a (0, 2) type tensor, it acts on the pairs of vectors V , W to produce
a real number T (V , W ) = Tij V i W j . Then T is called symmetric if
T (V , W ) = T (W , V ). (2.39)
2.2 The metric tensor 21
If {ei } is a basis for the tangent space, this implies T (ei , ej ) = T (ej , ei ),
that is
Tij = Tji . (2.40)
1
T(ij) = (Tij + Tji ), (2.42)
2!
1
T[ij] = (Tij − Tji ). (2.43)
2!
1 i
T[ijkl] = i
[T + Tklj i
+ Tljk − Tkjl
i
− Tlkj
i
− Tjlk
i
]. (2.44)
3! jkl
the metric tensor acts on pairs of vectors to give a number, and is symmetric
in its indices. Choosing a coordinate basis,
where gij = g (∂/∂xi , ∂/∂xj ). If V and W are any vectors, this gives
g (V , W ) = gij V i W j . This is written conventionally in the form of a distance
between two infinitesimally separated points in the spacetime as
Xi = gij X j , X i = g ij Xj . (2.48)
In particular,
g ik gkm = g i m = δ i m , (2.50)
i i
and the Kronecker delta δm transforms as components of a tensor, so δm and
i
gm are identical tensors.
The tensors T i j , T j i , or T ij are treated as representations of the same
geometric object because these are uniquely associated. Such an isomor-
phism between the covariant and contravariant arguments is equivalent to
the procedure of ‘raising’ and ‘lowering’ of indices as pointed out above. The
multilinear map
g : Tp × Tp → R (2.51)
If the metric has the form gij = (+1, . . . , +1) then it is called positive definite.
Then g (X, X ) = 0 implies X = 0. It is called a Lorentzian metric if the
form is
gij = diag(+1, . . . , +1, −1), (2.53)
with (n − 1) terms being positive. The metric is indefinite in the sense that
the magnitude of a non-zero vector could be positive, negative or zero. The
vector X ∈ Tp is called timelike, null, or spacelike, as defined by
The equations (2.56) and (2.55) are equivalent to the expression ds2 =
gij dxi dxj , which represents the infinitesimal arc length along γ .
2.3 Connection
The notion of parallel transport of a given vector X in Euclidian spaces can
be defined by requiring that, in going from a point p to q , both the magnitude
and direction of X must not change. If both these for the tangent vector
remain unchanged along a curve, it is called a straight line, along which
the tangent is parallel transported. In Euclidian space, if a vector is parallel
transported from points p to q along two different curves, the result is the
same, independent of the path taken. However, this is not necessarily the
case for a general affine manifold. For a differentiable manifold, the notion
of the parallel transport of vectors is defined by introducing a connection on
M.
Let X be a vector field on M , the derivative operator ∇X on M then gives
the rate of change of vectors or tensor fields along X at all p ∈ M . If Y is
another vector field at p, then the operator ∇X maps Y into a new vector
field Y → ∇X Y such that the following are satisfied:
∇( f Y ) = df ⊗ Y + f ∇Y . (2.57)
∇Y = Y i ; j e j ⊗ ei . (2.58)
∇∂/∂xi Y ≡ ∇i Y , (2.62)
using the rules defining the connection given above, and the relation
∂ ∂f
X( f ) = Xi i
( f ) = Xi i , (2.63)
∂x ∂x
one obtains
∂Y k ∂
∇X Y = X i
+ Γk ij Y j . (2.64)
∂xi ∂xk
Comparing this with (2.61)
∂
∇X Y = Y k
;i X
i
, (2.65)
∂xk
26 The spacetime manifold
where
∂Y k
Y k ;i ≡ + Γk ij Y j . (2.66)
∂xi
It can be seen that the components of the vector ∇X Y are given as Y k ;i X i ,
and
∂Y i
Y i, j ≡ (2.67)
∂xj
can be defined.
Then, taking the transformation of the coordinates {xi } → {xi } when
the basis vectors transform as ei → ei , it can be seen that Y i , j does not
transform like the components of a tensor. Similarly, consider the connection
coefficients in the new coordinate system,
Γk i j = ek , ∇ei ej . (2.68)
Transforming the dashed vectors to the original coordinate system and using
the conditions (2) and (3) above gives, in a coordinate basis,
∂xk ∂xi ∂xj k ∂ 2 xk
Γk i j = Γ ij + . (2.69)
∂xk ∂xi ∂xj ∂xi ∂xj
It follows that, because of the presence of the second derivative terms above,
the coefficients Γi jk also do not transform like the components of a tensor.
Consider, however,
∂ i j ∂
∇X Y = (Y i
;jX
j
) = (Y ; j X ) , (2.70)
∂xi ∂xi
which implies
i ∂xj ∂xi j
j i
Y ; j X = Y ; j j X . (2.71)
∂x ∂xi
Since the above is true for an arbitrary vector X j , Y i ; j are components of a
tensor.
Further, if Γi jk and Γi¯jk are components of two different connections on M ,
then using the coordinate transformations, it can be seen that the quantities
C i jk = Γi¯jk − Γi jk (2.72)
T i jk = Γi jk − Γi kj . (2.75)
Γi jk = Γi kj , (2.76)
and
∇X (S ⊗ T ) = ∇X S ⊗ T + S ⊗ ∇X T (2.78)
for any vector field X and tensor fields S and T . Furthermore, ∇ must agree
with the usual notion of a directional derivative, that is,
∂f
∇X f = df, X = Xf = X i . (2.79)
∂xi
Finally, ∇ must commute with contractions, that is,
As shown earlier,
with
∇X T i1 ...ir j1 ...js = T i1 ...ir j1 ...js ;a X a . (2.82)
Now, by considering the expansion for ∇i (ej ⊗ ek ) it is seen that
∇a ei = −Γi ac ec , (2.83)
with
∂ωk
ωk; j ≡
− Γi jk ωi . (2.85)
∂xj
In general, the covariant derivative of a tensor T can be written as
∂T i1 ...ir j1 ...js im
T i1 ...ir j1 ...js ;a = + Γ he T i1 ...e...ir j1 ...js
∂xa m
− e
Γ hjn T i1 ...ir
j1 ...e...js . (2.86)
n
This can also be seen by specializing to the frame of free fall. In such a
case, all the connection coefficients vanish and the metric is locally that of
special relativity. Then, gij = ηij and the partial derivatives of gij vanish.
Therefore, from the above equation for gij;k , one obtains gij;k = 0. As this
is a tensor equation, it must hold in all frames in general, and one can again
proceed as earlier.
2.4 Non-spacelike geodesics 29
DT
= T i...l k...m;h X h , (2.93)
∂t
where X is the tangent to γ . Then, γ is called a geodesic if its tangent vector
is parallel transported along it. That is, if X is the tangent vector field along
γ , then ∇X X is proportional to X . Therefore, there exists a function f such
that
∇X X = f X. (2.94)
Writing the components, this implies (X i ; j X j )ei = f X i ei always holds,
and so X i ; j X j = f X i along the geodesic curve. But, it is always possible
to reduce f to zero by a suitable choice of the curve parameter t along γ , so
the equation for the geodesic is written as
X i ; j X j = 0, (2.95)
where the X i values are components of the tangent vector to the geodesic.
The parameter t is called an affine parameter along γ , which is an
affinely parametrized geodesic. If {xi } denotes a local coordinate system,
the components X i are written as X i = dxi /dt and the equation for
geodesics is
d2 xi j
i dx dx
k
+ Γ jk = 0. (2.96)
dt2 dt dt
The affine parameter along the geodesic is determined up to an additive and
multiplicative constant. Thus, if t is an affine parameter, then so is t = at + b
and X i ; j X j = 0. Here b = 0 gives a new choice of the initial point γ (0) and
a= 0 implies a renormalization of the vector X .
A geodesic in (M, g ) is called timelike, spacelike, or null if its tangent
vector is timelike, spacelike, or null respectively. Here, the timelike or null
geodesics that represent the paths of particles or photons in the spacetime
are mainly considered. Since the tangent to a geodesic is parallel transported,
a timelike or null geodesic remains the same always and it cannot become
spacelike. In a Riemannian manifold with a positive definite metric, such
geodesics give the curves of shortest distance between its points. However, in
a spacetime with a Lorentzian metric, the non-spacelike geodesics maximize
the distance between the points, as defined by (2.56). If there is a timelike
geodesic between the points p and q , there is no shortest distance geodesic
30 The spacetime manifold
between them because, by introducing null geodesic pieces, one could always
join these points by curves of arbitrary small lengths. On the other hand,
any maximal length curve between p and q must necessarily be a timelike
geodesic.
The geodesic equations above are n equations in n variables xi with i =
1, . . . , n. Thus, the existence theorems for differential equations ensure that,
given xi and dxi /dt, that is, given any initial point p and the value of the
tangent vector X i , a unique geodesic through p with this value of tangent
exists. This can be used to define the exponential map Ep : Tp → M from
the tangent space at p into the spacetime. Under this map, any given tangent
vector X i in Tp is mapped to a point in M , a unit affine parameter distance
away along the unique geodesic determined by p and X i . It is clear that the
exponential map may not be defined on all of Tp because all the geodesics in
M passing through p may not extend to all the values of the affine parameter.
In such a case, M is called geodesically incomplete. On the other hand,
if the exponential map is defined on all of Tp for all points p, then M is
called geodesically complete. Then, all geodesics in M extend for all values
of their affine parameter. Also, the map Ep may not be one–one, because
the geodesics might cross each other. However, it can be shown (Bishop and
Critendon, 1964) that for a sufficiently small neighborhood Np of p there
is a neighborhood of the origin in Tp which is diffeomorphically mapped
onto Np by the exponential map which is one–one and well-defined on this
neighborhood. In such a case, the exponential map can be used to define
the normal coordinates on the neighborhood Np of p. Since Tp is an n-
dimensional vector space equivalent to Rn , the coordinates of any r ∈ Np can
be chosen to be the n-coordinates of the vector X p that is mapped onto it.
This coordinate system has the property that the geodesics are mapped into
straight lines and the connection coefficients vanish at p. Therefore, this
coordinate system turns out to be quite convenient for calculations at the
point p. The neighborhood Np can have a further property that any two
points in it can be joined by a unique geodesic contained totally within Np .
Such a neighborhood of p is called a convex normal neighborhood.
The geodesic equations are derived from the parallel transport required for
the tangent vector. If one requires that the curve must extremize the length,
namely that δl = 0, and works out
b
dxa dxb
1/2
δl = δ
gab dt dt
dt = 0, (2.97)
a
using the variational methods, it turns out that the resulting equations
are precisely the geodesic equations (2.96). So, the geodesics extremise the
lengths of curves between any two spacetime points. If the events p and q
are timelike related, and if there is a maximum length timelike curve from p
2.4 Non-spacelike geodesics 31
and to write the Lagrange equations, which are the equations of the space-
time geodesics. Then, by comparison with the geodesic equations, it is also
possible to evaluate the quantities Γi jk for the spacetime.
In the Minkowski spacetime, the surface t = 0 is a three-dimensional sur-
face with the time direction always normal to it. Any other surface of
t = const. is also a spacelike surface in the same way. In general, let S
be an (n − 1)-dimensional manifold. If there is a C ∞ -map φ : S → M that
is locally one–one (there is a neighborhood N for every p ∈ S such that φ
restricted to N is one–one) and φ−1 is also C ∞ as defined on φ(N ), then
φ(S ) is called an immersed submanifold of M . If φ is globally one–one, then
φ(S ) is called an embedded submanifold of M . It may also be required that
φ be a homeomorphism with the induced topology on φ(S ) from M . Lower
dimensional embedded submanifolds in M represent well-behaved surfaces
in the spacetime.
A hypersurface S of any n-dimensional manifold M is defined as an (n− 1)-
dimensional embedded submanifold of M , and Vp is denoted by the (n − 1)-
dimensional subspace of Tp of the vectors tangent to S at any p ∈ S . It
follows that a vector n ∈ Tp exists that is unique up to the scale, and that
is orthogonal to all the vectors in Vp . This is called the normal to S at p. If
the magnitude of n is either positive or negative at all points of S without
changing the sign, then n could be normalized so that gab na nb = ±1. If
gab na nb = −1, then the normal vector is timelike everywhere and S is called
a spacelike hypersurface. If the normal is spacelike everywhere with a positive
magnitude, S is a timelike hypersurface, and S is a null hypersurface if the
normal na is null at S .
The timelike geodesics could be used to define the synchronous coordinate
system in the neighborhood of a spacelike hypersurface in the spacetime as
below. Let ∇a be the metric connection that satisfies ∇a gbc = gbc;a = 0. Let
S be a spacelike hypersurface, then, for every p ∈ S , let γ be the unique
timelike geodesic with a tangent na , that is, the congruence of these curves
at points of S is orthogonal to S . Then, in the neighborhood of that portion
of S , the coordinates q → (x1 , . . . , xn−1 , t) are assigned for any point q in the
future of p along γ , where t is the parameter along γ and x1 , . . . , xn−1 are
the spatial coordinates of p. In particular, if the geodesics in the congruence
are parametrized by the proper time t with the magnitude of the tangent
given by −1 along γ , then the spacelike surfaces are given as {t = const.}
surfaces. Note that S can be labeled as the {t = 0} spacelike surface.
32 The spacetime manifold
The synchronous coordinates have the important property that when the
congruence {γ} is orthogonal to S0 , it will also be orthogonal to subsequent
surfaces St given by t = const. Clearly, {γ} is orthogonal to S0 by construc-
tion. To see that this holds for any St for t within the domain of construction,
let X a be any basis vector for the tangent space at a point of St . Then,
na ∇ a X b = X a ∇ a nb . (2.100)
This implies
∂
R(X, Y )Z = Ri jkl X kY lZ j . (2.105)
∂xi
Now, in order to evaluate (2.104), note that
Similarly,
[∇Y (∇X Z )]i = Z i ; jk X j Y k + Z i ; j X j ;k Y k . (2.107)
Finally,
∇k ∇l Z i − ∇l ∇k Z i = Ri jkl Z j . (2.111)
The above equation could also be taken as the defining equation for
the components of the curvature tensor. As shown by the left-hand side
of (2.111), the Riemann curvature tensor provides the measure of non-
commutation of a tangent vector when parallel transported along different
curves to arrive at the same spacetime point.
In place of the vectors X ,Y , and Z the basis vectors ei can now be chosen.
Then,
∂ Γi lj ∂ Γi kj
Ri jkl = − + Γi ka Γa lj − Γi la Γa kj . (2.114)
∂xk ∂xl
As pointed out earlier, given the metric tensor g on M , there is a unique,
torsion-free connection on M defined by the condition ∇X g = 0, which is
equivalent to the vanishing covariant derivative of the metric tensor, gij;k = 0.
34 The spacetime manifold
∇X (g (Y , Z )) = Xg (Y , Z )
= ∇X (gij Y i Z j )
= g (∇X Y , Z ) + g (∇X Z, Y ). (2.115)
Evaluating Y (g (Z, X )) and Z (g (X, Y )) and adding the first and subtract-
ing the second from (2.115) gives
The covariant derivatives of the Riemann tensor satisfy the Bianchi identities
given by Ri j[kl;a] = 0, which is the same as
A straightforward proof would involve writing down each term above explic-
itly, substituting from (2.114), and then taking a summation. There are
certain additional symmetries that are valid when the connection is the one
induced by the metric. In this case,
The Riemann tensor Rijkl defined by the metric has the symmetry
which means R(ij)kl = 0. Also, in this case the Riemann tensor is symmetric
in the pairs of the first two and last two indices,
The symmetry properties of the Weyl tensor follow from those of the Riemann
tensor discussed above; it possesses the same symmetries as the Riemann
tensor. Also, it can be verified that the following identically vanishes,
g ik Cijkl = 0. (2.129)
The Weyl tensor is that part of the curvature tensor for which all contractions
vanish for any pair of contracted indices,
C i jil = 0. (2.130)
If the Weyl tensor vanishes throughout the spacetime with Cijkl = 0 at all
points, then it can be shown that the metric gij must be conformally flat.
This means that a conformal function Ω(xi ), 0 < Ω < ∞, exists such that
where ηij is the flat Minkowskian metric. The Weyl tensor is conformally
invariant in the sense that under a conformal transformation gij → g ij =
Ω2 gij ,
i i
Cjkl = Cjkl . (2.132)
It is possible to show that a necessary and sufficient condition for the space-
time metric to be conformally flat is that the Weyl tensor must vanish
everywhere.
2.5 Spacetime curvature 37
The geodesic deviation equation, which is also called the Jacobi equation,
is now derived. This characterizes the coming together, or moving away,
of spacetime geodesics from each other as a result of the spacetime cur-
vature. Consider a smooth one-parameter family of affinely parametrized
non-spacelike geodesics, characterized by the parameters (t, v ), where t is
the affine parameter along a geodesic and v = const. characterizes different
geodesics in the family with t, v ∈ R. Such non-spacelike geodesics span a
two-dimensional submanifold on which t and v could be chosen as coordi-
nates. The vectors T = ∂/∂t and V = ∂/∂v are then coordinate vectors for
which [T , V ] = 0. Then, since the torsion tensor is vanishing,
∇T V = ∇V T , (2.133)
DV j = V i ∇i T j . (2.134)
D2 V j = DV i ∇i T j + V i D(∇i T j )
= (T k ∇k V i )(∇i T j ) + V i T l ∇l ∇i T j . (2.135)
∇l ∇i T j − ∇i ∇l T j = Rj kli T k . (2.136)
D2 V j = (V k ∇k T i )(∇i T j ) + ∇i ∇l T j V i T l + Rj kli T k V i T l
= V k ((∇k T i )(∇i T j ) + (∇k ∇l T j )T l ) + Rj kli T k V i T l
= V k (∇k (T i ∇i T j )) + Rj kli T k V i T l
= Rj kli T k V i T l . (2.137)
The equation
D2 V j = −Rj kil T k V i T l (2.138)
where τ is the proper time along the timelike trajectory describing the par-
ticle worldline. The simplest second rank tensor constructed from these two
quantities is
T ij = ρ0 ui u j . (2.139)
The component T 00 of this energy–momentum tensor is
dx0 dx0
T 00 = ρ0 . (2.140)
dτ dτ
In a special relativistic frame, this can be interpreted as the relativistic energy
density of matter. It can also be shown that requiring this tensor to have
zero divergence in such a frame gives the conservation of the energy and
momentum.
A perfect fluid is characterized by an additional scalar quantity, the pres-
sure p = p(xi ). In the limit as the pressure vanishes, this must reduce to the
dust form of matter. Furthermore, one also demands the conservation laws
in a special relativistic frame, and that these should reduce to the classical
equations of continuity and the Navier–Stokes equations in the appropriate
limits. Then this energy–momentum tensor is given in a general frame as
T ij = (ρ + p)ui uj + pg ij , (2.141)
T ij ; j = 0. (2.142)
The above equation for the stress–energy tensor contains much information
on the matter fields in a spacetime. For example, if the spacetime contains
a Killing vector K i then the above equation could be integrated to give a
conservation law. The conserved vector in such a case is defined as P i =
T ij Kj , and P i ;i = 0 as a consequence of (2.142) and the Killing equation
Then, the integration of P i ;i over a compact region implies that the total
flux over a closed surface of the energy–momentum is zero in the direction
of the Killing vector (Hawking and Ellis, 1973). Even when the spacetime
40 The spacetime manifold
does not admit a Killing vector, given any point p a Riemannian normal
coordinate system at p could be set up so that the metric components have
the Minkowskian values and the connection coefficients Γi jk vanish at p. A
small enough neighborhood of p could then be chosen so that the values of
gij and Γi jk differ by an arbitrarily small amount from the values at p. Using
this fact, it could be shown that isolated test particles should move along
timelike geodesics (Fock, 1939; Dixon, 1970).
Furthermore, all matter fields are assumed to obey the postulate of local
causality, that is, the equations governing the matter fields are such that for
any p ∈ M , there is an open neighborhood U of p in which a signal can be
sent between any two points of U if and only if there is a non-spacelike curve
joining these points. This principle is valid in special relativity and is also
accepted in general relativity.
The above principles effectively imply that it is the spacetime metric, and
the quantities derived from it, that must appear in the equations for phys-
ical quantities, and that these equations must reduce to the flat spacetime
equations when the metric is Minkowskian. This is the basic content of the
general theory of relativity, where the spacetime manifold is now allowed
to have topologies other than R4 , and the metric gij can be non-flat glob-
ally. In general relativity, the matter fields expressed by the stress–energy
tensor are related to the non-flat nature of the spacetime by means of the
Einstein equations, which are the basic equations satisfied by the spacetime
metric. In Einstein’s theory, one does not discuss the physical interaction
of matter fields in a fixed background metric prescribed in advance. Actu-
ally, the gij values are treated as dynamical variables that depend on the
matter content of the spacetime and that are to be solved from the Einstein
equations.
The Newtonian theory gives an important indicator towards obtaining this
relationship between the matter content and the spacetime geometry, where
the gravitational field is described by a potential φ. The tidal acceleration
between nearby particles is given in terms of the separation between them
and the second derivatives of φ. In a curved spacetime manifold, such tidal
accelerations are described by the Jacobi equation (2.138), in terms of the
Riemann curvature tensor. Furthermore, the Poisson equation
∇2 φ = 4πρ (2.144)
must be recovered in the Newtonian limit. Both in special and general theo-
ries of relativity, the matter content is described by the stress–energy tensor
Tij , and the mass–energy density ρ corresponds to the quantity Tij V i V j .
Therefore, each side of Poisson’s equation corresponds to the Riemann ten-
sor as expressed in the Jacobi equation and Tij V i V j respectively. Another
important indicator for the comparison is provided by the Bianchi identities.
2.6 The Einstein equations 41
∇a Ra l + ∇i Ri l − ∇l R = 0. (2.146)
∇i Gij = 0, (2.147)
Then, the contracted Bianchi identities actually imply the local conservation
of energy and momentum through the Einstein equations. Taking the trace
of the above,
R = −8πT. (2.150)
Substituting this back in the above gives an alternative form of the Einstein
equations,
Rij = 8π (Tij − 12 T gij ). (2.151)
The definition of the Ricci tensor suggests that the Einstein equations
depend on the derivatives of the metric up to the second order. These equa-
tions are highly non-linear in gij , however, they are linear in the second
derivatives of gij . In fact, the quantities Rij and Rgij are the only second-
rank symmetric tensors that are linear in the second derivatives of the metric,
and involve only up to the second derivatives of gij . Actually, the Einstein
equations are a coupled system of non-linear second order partial differ-
ential equations for gij . This makes the task of solving these extremely
difficult. Several symmetry assumptions on the spacetime generally need
to be imposed in order to work out the metric components as a solution
to the Einstein equations. Some solutions that are useful in the context of
gravitational collapse are discussed later.
42 The spacetime manifold
Given the T ij values, the field equations may be viewed as a set of differ-
ential equations to determine the gravitational potential values, g ij , that fix
the resulting geometry. A particularly important case here is that of vacuum
solutions when T ij = 0. On the other hand, one could arbitrarily specify
the ten metric potentials and then compute the Einstein tensor Gij . Then,
the field equations determine the energy–momentum tensor Tij . However,
in that case, the resulting T ij values turn out to be unphysical most of the
time in that the energy conditions ensuring the positivity of mass–energy
density may be violated. Such a violation of the energy conditions is rejected
on physical grounds as all observed classical fields obey positivity of energy
density, which is closely connected with the physical features of gravitation
theory.
In general, the field equations are ten equations, connecting a total of
twenty quantities, which are the ten components of gij and the other ten
components of Tij . Therefore, the field equations are the conditions placing
constraints on the simultaneous choice of these twenty quantities. If part of
the gravitational potentials and the matter contents are determined from
physical conditions, then such conditions can fully determine the matter and
geometry. In particular, if the vacuum equations are considered,
then there are ten equations to determine the ten quantities gij . However,
the Bianchi identities
∇j Gij = 0 (2.153)
place four differential constraints on these equations, which are not all
independent. Therefore, there is an indeterminacy in that there are fewer
equations than unknowns to be determined. Furthermore, there is an intrin-
sic gauge freedom available in the general theory of relativity that does not
allow a complete determination of the metric potentials. This is given by the
coordinate freedom that allows a transformation from one set of coordinates
xi to any other set of coordinates xi . However, this coordinate freedom
could be used to impose conditions on the metric components. For example,
choosing the normal coordinates gives g00 = 1 and g0α = 0, α = 1, 2, 3 in this
coordinate system. This leaves six other components to be determined from
the field equations. The issue is closely linked to the Cauchy problem in gen-
eral relativity where the basic problem is that, given initial data on a regular
spacelike hypersurface, one would like to determine its unique evolution in
the future or past.
The Einstein equations can admit a cosmological constant. Note that the
most general second rank tensor which can be constructed out of Gij and gij ,
which is divergence free and involves the derivatives of the metric tensor up
to second order only, is the linear combination Gij + Λgij (Lovelock, 1972),
2.7 Exact solutions 43
It is seen that for an empty spacetime with Tij = 0, the Einstein equations
then become
Rij = −Λgij . (2.155)
If Λ = 0, then one does not obtain the Newtonian theory in the limit of
slow motions and weak fields. However, if the magnitude of Λ is very small,
then such departures will be quite negligible and approximate agreement
with the Newtonian theory is obtained.
The range of the coordinates r, θ, φ are 0 < r < ∞, 0 < θ < π , and
0 < φ < 2π . Two such coordinate neighborhoods are needed to cover all
of the Minkowski spacetime.
In coordinates (t, x, y, z ) it is obvious that the geometry is flat, because
all the metric components are constants, so all the connection coefficients
are vanishing. In other coordinate systems, such as spherical coordinates
(t, r, θ, φ), the connection coefficients Γi jk do not all vanish (for example,
Γ1 22 = r); however, all the Riemann curvature tensor components do still
vanish.
The Lorentz transformations on the Minkowski spacetime are defined as
the set of metric preserving isometries that are linear and homogeneous.
Physically, these represent the change of reference frame from one to another
inertial observer, given by the coordinate change
xi → xi = Li j xj . (2.158)
contained inside the chronological past or the future of the origin and there
are timelike geodesics outside this past or future cone that do not meet any
of these surfaces.
This leaves the causal structure of M invariant, because the null geodesic
paths of the metric ηij and the unphysical metric gij are the same under a
conformal mapping, up to a reparametrization. Therefore, the past of any
non-spacelike curve γ is unchanged and there is a natural correspondence
between ideal points in both the spacetimes. Since light cones are unaltered
by a conformal transformation, the boundary attachment obtained in this
manner is coordinate independent.
In the Minkowski spacetime, the advanced and retarded null coordinates
can be introduced, where
v = t + r, u = t − r, (2.161)
which gives a reference frame, based on null cones, which is most suitable to
analyze the radiation fields.
Under this transformation of coordinates, the metric becomes
with −∞ < v < ∞ and −∞ < u < ∞. Now, the information at future
null infinity corresponds to taking the limit as v → ∞, which amounts to
moving in the future along u = const. light cones. Similarly, past null infinity
corresponds to u → ∞. This procedure could be made precise in a coordinate
independent way. The Minkowski spacetime M can be compactified by means
48 The spacetime manifold
and then by adding the closure to add the null infinities. New coordinates p
and q can be introduced by
and the metric ḡij on the unphysical spacetime M , after the conformal
transformation, is given by
It is possible to see now that the metric above, with the coordinate ranges as
given above, is a manifold embedded as a part of the Einstein static universe.
To see this, let
T = p + q, R = p − q, (2.167)
and then (2.166) becomes, in (T, R, θ, φ) coordinates,
empty exterior solution where the Ricci tensor vanishes and is matched at the
boundary to an interior solution inside the body. In (t, r, θ, φ) coordinates,
the metric can be given as
−1
2m 2m
ds = − 1 −
2
dt + 1 −
2
dr2 + r2 dΩ2 , (2.170)
r r
where dΩ2 = dθ2 + sin2 θ dφ2 . The coordinate t is timelike and r, θ, φ are
spacelike coordinates. The radial coordinate r has the property that the
two-sphere given by t = const., r = const. has the two-metric given by
2Φ 2GM
g00 1 + 2
=1− 2 , (2.172)
c c r
m2
Rijkl Rijkl = (2.174)
r6
diverges, so the point r = 0 is a real spacetime singularity. It is not possible
to extend the spacetime across this singularity in a continuous manner. A
maximal extension of the Schwarzschild manifold, which covers both the
patches above, was obtained by Kruskal (1960) and Szekeres (1960).
The vacuum Schwarzschild geometry exterior to the collapsing matter
arises for a dust cloud collapse. In the case of homogeneous collapse, a
Schwarzschild blackhole is necessarily produced in a complete gravitational
collapse, and the event horizon fully covers the resulting spacetime singular-
ity of infinite curvature and density. In other words, as the collapse evolves,
the event horizon develops earlier than the formation of the singularity. This
situation is significant for the cosmic censorship hypothesis and blackhole for-
mation scenarios. The interior metric in this case is exactly that of a closed
Friedmann model in a time reversed sense.
In the extended Schwarzschild manifold, r = 2m is a null hypersurface
and each point on this surface is a two-sphere of area 16πm2 . Note that the
metric component g00 = (1 − 2m/r) is positive for r > 2m; however, g00 < 0
for r < 2m. Therefore, it is no longer possible to use t as a time coordinate in
that range, as the coordinates t and r reverse their roles there, and spacetime
is no longer static. Therefore, the r = 2m surface is also called the static
limit. The vector ∂/∂t with components ξ i = δ i 0 = (1, 0, 0, 0) gives the time
translation, leaving the gij unchanged as it does not involve the time coor-
dinate. Therefore, ξ is a Killing vector that leaves the spacetime geometry
unchanged. One obtains ξ 2 = gij ξ i ξ j = g00 , and for the Schwarzschild metric
ξ 2 vanishes on r = 2m. Hence, at the static limit the timelike Killing vector
becomes null.
The Schwarzschild geometry is an illustration of the basic principle that
Einstein used to formulate gravitation theory, namely that matter tells the
spacetime how to curve. To see this, consider the Schwarzschild solution in a
spacelike surface t = const. and in the equatorial plane θ = π/2. The metric
of this two-dimensional curved surface is given by
dr2
ds2 = + r2 dφ2 . (2.175)
(1 − 2m/r)
Consider a static observer along a Killing direction, for whom the four-
velocities are ui = ξ i /|ξ|. Suppose a static source with four-velocity ui1 emits
a photon with four-momentum pi (so pi; j pj = 0 with a suitable parametriza-
tion) and is observed by a static observer with four-velocity ui2 . Now, take
the directional derivative of ξi pi along the geodesic tangent pi ,
The first term vanishes because ξ i is a Killing vector and the second term
vanishes because of the geodesic equation. Therefore, the ratio of energies
measured at these two points by static observers is given by
i 1/2
E1 u pi
= i 1 1/2 . (2.177)
E2 (u pi )2
Since ξ 2 = g00 , this is the gravitational red-shift formula for a static source
and observer in terms of the metric components. It is now seen that if the
observer remains at a finite radius but the source approaches r = 2m, the
red-shift tends to infinity. Therefore, as a particle falls into the blackhole,
approaching r = 2m, the light rays emitted by it are infinitely red-shifted as
observed by a distant static observer in the outside spacetime.
collapsing dust ball. This is the same line element as that of the closed
homogeneous and isotropic Friedmann models given by
dr2
ds = −dt + R (t)
2 2 2 2 2
+ r dΩ , (2.179)
1 − r2
dt r
= . (2.180)
dr r − 2m
t = r + 2m ln | r − 2m | + const. (2.181)
In the region r > 2m we have dr/dt > 0 and hence r increases with increasing
t. Thus the above describes the congruence of outgoing radial null geodesics.
The in going null trajectories are given by the negative sign solutions, which
are given by
t = −r − 2m ln | r − 2m | + const. (2.182)
54 The spacetime manifold
In the region below r = 2m, the coordinates r and t change their space-
like and timelike nature and hence the light cones tip over. As a result, no
observer in the region r < 2m can remain at a constant value of r, but they
must move inwards to fall within the intrinsic curvature singularity at r = 0.
Each point in the (t, r) plane represents a two-sphere of area 4πr2 .
This gives an idea of the phenomena happening in the region r < 2m,
namely that any material particle or photon here must fall in the spacetime
singularity, and that it cannot escape to larger values of r to communicate
with external observers in the spacetime. Hence, this is called the black-
hole region in the spacetime. However, the Schwarzschild picture gives the
impression that from the outside r > 2m, no photons or particles could
fall in this blackhole and they will take infinite time to reach the surface
r = 2m. It turns out that this is actually a coordinate defect arising due
to the coordinate singularity at r = 2m, as can be seen by going to the
Eddington–Finkelstein coordinates, where the idea is to choose a new time
coordinate such that the in going null geodesics become straight lines in the
spacetime. It is clear from the Schwarzschild consideration above that the
appropriate change may be given by
for the r > 2m region. The solution with such a coordinate change is now
regular at r = 2m and the coordinate range is now 0 < r < ∞. This is
called an analytic extension of the Schwarzschild solution. A time-reversed
solution is obtained if a different time coordinate is introduced,
in which case the outgoing null geodesics are straight lines. A simpler way to
write the metric in the new coordinate system is to introduce the advanced
null coordinate v defined by
Spacetime singularity
Event horizon
r=0
r = 2m
Apparent horizon
Fig. 2.3 Homogeneous dust cloud collapse. The trapped surfaces form when the
star enters a r = 2m radius. The event horizon forms prior to the singularity,
creating a blackhole as the collapse endstate.
where ki is a null vector radially directed outwards. The metric is best given
in the null coordinates (u, r, θ, φ) as
2m(u)
ds = − 1 −
2
du2 − 2du dr + r2 dΩ2 , (2.188)
r
the star, and is matched by the Schwarzschild metric in the exterior of the
radiation zone.
The form for the energy–momentum tensor σ is defined to be the energy
density of the radiation as measured locally by an observer with a four-
velocity vector v i . Thus, σ is the energy flux as well as the energy density
measured in this frame,
σ ≡ Tij v i v j , (2.189)
with v i vi = −1. Working out the connection coefficients, the Ricci tensor in
null coordinates is given by
2 dm(u) 0 0
Rij = − δ iδ j . (2.190)
r2 du
This implies that the Ricci scalar Ri i = R = 0, and hence the Einstein
equations give
1 dm(u) 0 0
Tij = − δ iδ j , (2.191)
4πr2 du
which is the energy–momentum tensor of a radiating field in the geometric
optics form. From (2.189) and (2.191),
1 dm(u)
σ=− , (2.192)
4πr2 du
which is the expression for the energy density of radiation.
In the case when m(u) = const., the relationship of the null coordinates
used here with the Schwarzschild coordinates (t, r, θ, φ) is not difficult to see.
In such a case, the transformation given by Finkelstein (1958) can be used
to diagonalize (2.188),
given by
dm
L∞ (u) = lim 4πr2 q = − , (2.195)
r→∞,U =0 du
that is, it is the negative rate of change of mass of the radiating body.
The surface r = 2m(u) has many interesting properties, as pointed out
by Lindquist, Schwartz, and Misner (1965). Unlike the Schwarzschild case,
where r = 2m is a null hypersurface, for the Vaidya radiating star metric,
this is a spacelike hypersurface. The induced metric on this hypersurface is
given by
dm
ds |r=2m(u) = −2
2
du2 + r2 dΩ2 . (2.196)
du
This induced metric has the signature (+, +, +) since dm/du < 0. As a
result, the position of light cones on this surface is such that for all timelike
vectors in the forward light cone at all points on this surface, dr/du > 0.
Therefore, no timelike trajectory from the outside region r > 2m(u) can
come and cross this surface to enter inside the r < 2m(u) region. The solution
is of type D in the Petrov classification of spacetimes and possesses a normal
shear-free congruence with a non-zero expansion.
The Vaidya solution has been used extensively in the context of the cosmic
censorship hypothesis, and to study gravitational collapse endstates (see for
example, Papapetrou, 1985; Joshi, 1993, for details). In that case, imploding
radiation shells are considered, rather than the outgoing case considered here.
Then, the function m is taken to be non-decreasing and the advanced null
coordinate t + r is used. For the case of a linear mass function m(u) = λu,
an interesting causality and horizon structure arises, depending on the rate
of collapse ṁ(u), that is, as decided by the magnitude of λ. In this case,
the model is self-similar in that it admits a homothetic Killing vector. When
the collapse is fast enough, the horizon forms well before the singularity to
fully cover it, and a blackhole forms as the collapse endstate. However, for a
slow enough collapse when λ is smaller than a critical value, the horizon and
trapped surface formation is delayed and a naked singularity develops as the
collapse final state (see Fig. 2.4).
The above scenario may be of physical interest if, in the very late stages
of collapse, the star converts itself into some kind of an imploding radiation
ball. On the other hand, it can be argued that real collapsing stars may not
have the form of the energy–momentum tensor of a radiation fluid in the very
late stages of their collapse. In any case, not much is known today on the
actual equation of state for the collapsing star in the very late stages of its
gravitational collapse. The above consideration really shows what is actually
possible within the framework of general relativistic gravitational collapse,
subject to various physical reasonability and regularity conditions, such as
the validity of energy conditions and evolution from regular initial data.
2.7 Exact solutions 59
Event horizon
Singularity
at
r=0 Outgoing
r = 2m
null geodesic
Incoming matter
60
Spherical collapse 61
Here, the quantities A, B, C , and D are the functions of t and r that are
to be determined. Introducing a new radial coordinate r = D1/2 and a new
time coordinate t by requiring that dt = F [A dt − B dr ], where F (t, r ) is
a suitable integrating factor, the line element reduces to
where the primes have been dropped. Here, ν = ν (r, t) and ψ = ψ (r, t), and
the quantities eν = 1/F 2 A and eψ = C − B 2 /A appearing in the metric
are always positive.
In general, the spherical symmetry of a spacetime can be defined using
Killing vectors. There must be three linearly independent spacelike Killing
3.1 Basic framework 63
[X 1 , X 2 ] = X 3 , [X 2 , X 3 ] = X 1 , [X 3 , X 1 ] = X 2 , (3.4)
and their orbits must be closed. Using these properties, the line element
above could be derived rigorously for a spherically symmetric spacetime.
A general formalism for a spherically symmetric gravitational collapse
including pressure, was developed by Misner and Sharp (1964). The
spherically symmetric spacetime is written in comoving coordinates as
where dΩ2 = dθ2 + sin2 θ dφ2 is the usual metric on the two-sphere, and φ
and λ are functions of t and r. The stress–energy tensor is that of a perfect
fluid given by
Tij = (ρ + p)ui uj + pgij . (3.6)
The spatial velocities ui are vanishing here and the spatial coordinates of a
given particle remain constant throughout the collapse. A function m(r, t)
is introduced by the definition
−1
2m
λ
e = 1 + Ṙ −
2
R2 , (3.7)
r
where a dash denotes a derivative with respect to r, and for any function f,
˙ −φ ∂f
f =e . (3.8)
∂t
The coordinate t here gives the proper time along the particle world lines.
Integrating the conservation equation T ij ;j = 0 and solving the Einstein
equations, the Misner–Sharp equations for the spherically symmetric collapse
can be written as
(see for example, May and White, 1966) can be opted for in order to get an
idea of the evolution of the collapsing system.
The general spherically symmetric metric, in a spacetime of dimension N ,
can be written in the form
N −2
N
Gab = [gab {2RR + (N − 3)Q} − 2RR,ab ] , (3.14)
2R 2
1
N
G22 = − R2 + (N − 3){2RR + (N − 4)Q} , (3.15)
2
i−2
N
Gii (i > 2) = sin2 (θk ) N G22 , (3.16)
k=1
where
and
= g ab ab . (3.18)
Here, ab is the Ricci tensor evaluated from the two-metric gab , and is the
scalar curvature also evaluated from the two-metric gab .
The spacetime geometry within the spherically symmetric collapsing cloud
is described by comoving coordinates (t, r, θ i ), which are specified below.
The matter field is chosen to be of general type I, which is a broad class
including most of the physically reasonable matter forms, including dust,
perfect fluids, massless scalar fields, and such others. This class of matter
is specified by the requirement that the energy momentum tensor for the
matter admits one timelike and three spacelike eigenvectors (Hawking and
Ellis, 1973; Stephani et al., 2003).
The coordinates (t, r, θ i ) are then chosen to be those along these eigenvec-
tors, which makes the coordinate system comoving, that is, the coordinate
system moves with the matter. The freedom of coordinate transformations
3.1 Basic framework 65
of the form t = f (t, r) and r = g (t, r) can be used to make the gtr term in
(3.12) and the radial velocity of the matter vanish. In this case, the general
metric in the comoving coordinates (t, r, θ i ) must have three general arbi-
trary functions of t and r, and this can be written in the form (Landau and
Lifshitz, 1975)
In this comoving frame, the energy–momentum tensor for any matter field
that is type I is given in a diagonal form
i
T t t = −ρ(t, r), T r r = pr (t, r), T θ θi = pθ (t, r). (3.20)
The quantities ρ, pr , and pθ are the energy density, and radial and tangential
pressures ascribed to the matter field respectively.
The matter cloud can be chosen to have a compact support at an ini-
tial spacelike surface of t = ti , with 0 < r < rb , where rb would denote
the boundary of the cloud. Outside this boundary, the interior solution has
to be matched through suitable junction conditions, with another suitable
spacetime metric to complete the full spacetime geometry.
The matter fields are taken to satisfy the weak energy condition, that is,
the energy density measured by any local timelike observer is non-negative.
This ensures the physical reasonability for the collapsing matter fields consid-
ered here. Another energy condition frequently used is the dominant energy
condition, which demands that, for any timelike observer, the local energy
flow is non-spacelike. These two conditions are frequently regarded as the
main and important energy conditions that are physically reasonable (Hawk-
ing and Ellis, 1973; Wald, 1984). All classical observed matter fields satisfy
these conditions. For the energy conditions to be satisfied, for any timelike
vector V i,
Tik V i V k ≥ 0 (3.21)
ρ ≥ 0, ρ + pr ≥ 0, ρ + pθ ≥ 0, (3.22)
|pr | ≤ ρ, |pθ | ≤ ρ. (3.23)
Now, with the above metric, the following quantities can be evaluated
R,a R,a = Ṙ2 e−2ν − R 2 e−2ψ , (3.24)
R = e−2ν [R̈ + Ṙ(−ν̇ + ψ̇ )] − e−2ψ [R + R (ν − ψ )], (3.25)
66 Spherical collapse
and
e2ν ν R e2ψ ψ̇ Ṙ
R,ab = R̈ − ν̇ Ṙ − δa0 δb0 + R − ψ R − δa1 δb1 + Qab ,
e2ψ e2ν
(3.26)
where
Qab = Ṙ − ν Ṙ − ψ̇R δa0 δb1 + δa1 δb0 . (3.27)
Using the above equations, the Einstein equations Gik = Tik now take the
form (in the units 8πG = c = 1)
(N − 2)F (N − 2)Ḟ
ρ= −2
, pr = − , (3.28)
2R N R 2RN −2 Ṙ
(N − 2)(pθ − pr ) R pr
ν = − , (3.29)
ρ + pr R ρ + pr
Ġ H
−2Ṙ + R + Ṙ = 0, (3.30)
G H
F
G − H = 1 − N −3 , (3.31)
R
where
G(t, r) = e−2ψ (R )2 , H (t, r) = e−2ν (Ṙ)2 . (3.32)
The arbitrary function F = F (t, r) here has the interpretation of the mass
function for the cloud, which gives the total mass in a shell of comoving radius
r, at an epoch t. The energy condition ρ ≥ 0 implies F ≥ 0 and F ≥ 0.
Since the area radius vanishes at the center of the cloud, from (3.28), it is
evident that in order to preserve the regularity of density and pressures at
any non-singular epoch t, F (t, 0) = 0, that is the mass function vanishes at
the center of the cloud.
As seen from (3.28), there is a density singularity in the spacetime at R = 0
and at R = 0. However, the latter ones are due to shell-crossings (Yodzis,
Seifert, and Muller zum Hagen, 1973, 1974), which basically indicate the
breakdown of the coordinate system used. These are not generally regarded
as genuine spacetime singularities, because they can be possibly removed
from the spacetime to extend the manifold (Clarke, 1986, 1993). Hence, only
the shell-focusing singularity at R = 0, which is a genuine physical singularity
where all matter shells collapse to a zero physical radius, will be considered
here. This will be discussed in more detail later.
Note that, in general, for a general matter field with non-vanishing pres-
sures as considered here, there are a variety of dynamical time evolutions
possible from the given matter density and pressure profiles as prescribed
on an initial surface (that are called matter initial data here), from which
3.1 Basic framework 67
Note that not all the initial data functions above are mutually independent,
as, from (3.29),
r
(N − 2)(pθ0 − pr0 ) pr0
ν0 (r) = − dr. (3.37)
0 r(ρ0 + pr0 ) ρ0 + pr0
Therefore, apart from the matter initial data, which describe the initial den-
sity and pressure profiles at the initial epoch t = ti , the rest of the initial
data that are free is ψ0 (r), which essentially describes the velocities of the
collapsing matter shells as discussed later.
Note that in the above, when the pressures are taken to be vanishing, the
scenario reduces to a dust collapse. In this case, the only two free initial
data functions are the initial density and the velocity profile for the cloud.
As discussed later, in this case, for any given initial density profile ρ(r)
for the matter cloud, there are classes of initial velocity profiles ψ (r) that
take the final collapse outcome to either a blackhole, or a naked singularity,
depending on the choice of initial velocities. This is subject to the validity
of the energy conditions and regularity conditions for gravitational collapse.
The Oppenheimer–Snyder homogeneous dust collapse corresponds here to a
special choice of the initial density profile that is a constant function, and
another special choice of the initial velocity function. In fact, if the initial
density profile is chosen to be constant, but the initial velocities are allowed to
take different functional forms, then an initially constant and homogeneous
density can also later inhomogenize, and result in a final evolution that is a
3.2 Regularity conditions 69
and the mass function. This gives the freedom of choice of two free functions.
Selection of these free functions, subject to the weak energy condition and
the given regular initial data for collapse at the initial surface, determines the
matter distribution and metric of the spacetime throughout, and thus leads
to a particular time evolution of the initial matter and velocity distributions
with which the collapse began. As will be shown, it turns out that, given the
matter initial profiles in terms of ρ0 , pr0 , and pθ0 , there are the rest of the
initial data at t = ti , and the classes of solutions to the Einstein equations
that are explicitly constructed here, which give either a blackhole or a naked
singularity as the endstate of collapse. This outcome depends on the nature
of the rest of the initial functions, and the classes of the dynamical evolutions,
as allowed by the Einstein equations.
An important point to be noted here is that, in the description above no
mention has been made so far of the equations of state that the matter must
obey. Typically, these are of the form, pr = pr (ρ) and pθ = pθ (ρ). If these
are specified, then there is no freedom left, and there are seven equations for
seven variables. If this were to be incorporated right away, the only way to
proceed to find the collapse endstate would be to assume a specific equation
of state that the matter must satisfy. Then, one has to examine the collapse
problem accordingly, and examine the nature of the final singularity resulting
from the dynamical evolution, as governed by the Einstein equations. There
have been many collapse studies in the past using this approach, such as for
example, for the dust equation of state for perfect fluids, and others. The
limitation of such an approach, however, is that there is very little existing
knowledge on what a realistic equation of state should be that the matter
has to satisfy at the extreme high densities that a continual collapse realizes
in its advanced stages. For example, even for neutron star densities that are
relatively low compared with those of the later stages of a continual collapse,
there is a great deal of uncertainty on the equation of state for such neutron
matter. As a result, the neutron star mass limits are uncertain. Therefore,
specific or special assumptions used on the equation of state may turn out
to be physically unrealistic or restrictive, and untenable in the final stages
of the collapse. In fact, diametrically opposite views exist on the possible
equation of state in the very late stages of the collapse. For example, while
there are many arguments suggesting that pressures must play an important
role in the later stages of the collapse, the other view is that in such late
stages the matter must necessarily be dust-like (see for example, Hagedorn,
1968; Penrose, 1974a, 1974b).
Under this situation, the approach taken here is that no specific or partic-
ular equation of state is assumed to begin with. All further considerations are
carried out in a general manner, in terms of the allowed initial matter profiles
and the allowed dynamical evolutions of the Einstein equations, towards
determining the blackhole and naked singularity endstates for collapse.
3.3 Collapsing matter clouds 71
Subsequently, in Section 3.7, the role that the equation of state plays towards
further fine tuning the blackhole and naked singularity outcomes as collapse
endstates is discussed. The advantage of such an approach here is that, first
all the collapse equations are written in general, and only then are different
subcases distinguished, depending on the corresponding equation of state
under consideration. As shall be shown, various important subcases and
equations of state such as dust, perfect fluids, and others are included as
special cases of the treatment given here.
Exterior
Time
Interior
Fig. 3.1 The collapse evolves from regular initial data from an initial space-
like hypersurface. The cloud has a compact support on the spacelike slice and
the interior collapse spacetime is matched to a suitable exterior metric at the
boundary.
exist. No such families exist, originating from the singularity, when the end
product of collapse is a blackhole. This is a definition that is independent of
the coordinates used and the slicing chosen. In the case of a blackhole devel-
oping, the resultant spacetime singularity will be hidden inside an event
horizon of gravity, remaining unseen by external observers. On the other
hand, if the collapse ends in a naked visible singularity, there is a causal
connection between the region of singularity and faraway observers, thus
enabling, in principle, a communication from the super-dense regions close
to the singularity to faraway observers.
Given the matter initial profiles in terms of the initial functions as given
by ρ0 (r), pr0 (r), pθ0 (r) at the initial epoch t = ti from which the collapse
commences, the purpose now is to construct and examine possible evolutions
(classes of solutions to the Einstein equations) of such a matter cloud from
the given initial data, to investigate its final states (see Fig. 3.1).
While constructing the classes of solutions that give the collapse evolutions,
given the matter initial data at t = ti , as much generality as possible is
preserved. The mass function F (t, r) for the collapsing cloud must have the
following general form:
F (t, r) = r3 M(r, v ), (3.39)
where M > 0 is at least a C 1 -function of r for r = 0, and at least a C 2 -
function for r > 0. It is to be noted that F must have this general form,
which follows from the regularity and finiteness of the density profile at the
initial epoch t = ti , and at all other later regular epochs before the cloud
collapses to the final singularity at R = 0. This requires, from the Einstein
equation for energy density, that F must behave as r3 close to the regular
center. Hence, note that since M is a general (at least C 2 ) function, (3.39)
is not really any ansatz or a special choice, but quite a generic class of the
mass profiles for the collapsing cloud, consistent with and allowed by the
regularity conditions. Therefore, no special choice of F is made, but it is
3.3 Collapsing matter clouds 73
3M + r [M,r + M,v v ]
ρ(r, v ) = , (3.40)
v 2 (v + rv )
and
M,v
pr (r, v ) = − . (3.41)
v2
The regular density distribution at the initial epoch is given by
That is, A(r, v ),v ≡ ν /R , and since at t = ti one has R = r, this gives
[A(r, v ),v ]v=1 = ν0 (r). The main interest here is in studying the shell-focusing
singularity at R = 0, which is the physical singularity where all the matter
shells collapse to zero physical radius. Therefore, assume that there are no
shell-crossing singularities in the spacetime where R = 0, and that the func-
tion A(r, v ) is well-defined. From (3.38), the form of ν (t, r) can be generalized
and chosen as the class of functions given by
Here, b(r) is another arbitrary function of the shell radius r. By using the
regularity condition on the function v̇ at the center of the cloud,
where b0 (r) is the energy distribution function for the shells. Using (3.43) in
(3.29),
Rpr
pθ = RA,v (ρ + pr ) + 2pr + . (3.48)
R
3.3 Collapsing matter clouds 75
It can be seen that in general, both the density and radial pressure blow up
at the singularity, so the above equation implies that the tangential pressure
would also typically blow up at the singularity.
Now, using (3.39), (3.43), and (3.46) in (3.31),
R1/2 Ṙ = −eν(r,v) (1 + r2 b0 )Re2r A(r,v) − R + r3 M. (3.49)
The negative sign on the right-hand side of the above equation corresponds
to a collapse scenario where Ṙ < 0. Defining a function h(r, v ) as
e2r A(r,v) − 1
h(r, v ) = = 2A(r, v ) + O(r2 ), (3.50)
r2
(3.31) becomes
v 1/2 v = − e(rA+v) vb0 + e2v (vh + M). (3.51)
Note that the variable r is treated as a constant in the above equation. The
above equation gives the time taken for a shell labeled r to reach a particular
epoch v from the initial value v = 1. Expanding the function t(v, r) around
the center of the cloud, provided the functions within the integral above are
sufficiently regular, gives
where
Hence, it can be seen that the time taken for a shell labeled r to reach the
spacetime singularity at R = 0 (which is the singularity curve) is given as
1
v 1/2 dv
ts (r) = . (3.56)
0 e(rA+ν) vb0 + e2ν (vh + M)
76 Spherical collapse
v (N −3)/2
Q(r, v ) = . (3.57)
e(rA+ν) v (N −3) b0 + e2ν (v (N −3) h + M)
Let A be the rectangular area in the (r, v ) plane defined by the lines
r = 0, r = , v = 0, v = 1. (3.59)
d d2
φ1 (r) = [ts (r)], φ2 (r) = [ts (r)], (3.60)
dr dr2
and this implies that the singularity curve ts (r) would be a well-defined
C 2 -function near the center.
Several well-studied collapse models such as dust collapse models and
others, satisfy these or stronger conditions, and so the singularity curve is
well-defined and expandable. Some examples of such singularity curves that
form special cases for the consideration above follow.
In the case of a dust collapse, the expression for the singularity curve is
given by (Goswami and Joshi, 2004a, 2004b)
1
v 1/2 dv
ts (r)dust = , (3.61)
0 M(r) + vb0 (r)
where M(r) and b0 (r) are well-defined C 2 -functions of the comoving coor-
dinate r and are well-defined at r = 0. It can then be easily seen that the
3.3 Collapsing matter clouds 77
where
Again, since M(r) and b0 (r) and L(r) are well-defined C 2 -functions of the
coordinate r, it can be seen that the above singularity curve is well-defined
and differentiable at r = 0.
To generalize this, and to give another explicit example, it can be shown
that (Goswami and Joshi, 2006) given any initial data of the form (3.36),
classes of dynamical evolutions that give rise to a well-defined and differ-
entiable singularity curve always exist. By the freedom of choice of free
functions, the evolution functions M(r, v ) and A(r, v ) can be chosen in the
following way,
The Einstein equations here imply that, for the above class of evolutions,
the radial pressure remains static. However, the tangential pressure blows up
along with the density at the singularity v = 0 and is given by
It is now easy to check that the above class of evolutions admits a well-
defined and differentiable singularity curve, as both the functions M(r, v )
and A(r, v ) are well-defined and are C 2 at r = 0 and v = 0.
The point is that, as long as in the construction of classes of the solutions
the functions M and A are taken to be sufficiently regular, then the sin-
gularity curve also becomes regular and expandable, at least up to the first
order. If the singularity curve is not regular even to that extent, it means
that it has no well-defined tangent at the origin, and so the collapse evolu-
tion is not regular enough and may not be of physical interest. In any case,
the purpose here is not to analyze all possible collapse evolutions from a
given initial distribution of matter and velocities. The aim really is to show
that, given any such regular distribution from which the collapse develops,
there exist classes of dynamical evolutions that take the other collapse to
either the naked singularity or the other collapse endstates, depending on
the choice of allowed functions. It is also shown that well-known classes of
collapse models, such as dust, perfect fluids, and others, form special classes
of the treatment given here.
Once a singularity curve that is at least C 2 exists, the function can be
Taylor expanded near the center as
In other words, a continual collapse in finite time ensures that the above
condition holds. Also, from (3.51) and (3.53), for small values of r, along
constant v surfaces,
v 1/2 v = (vb00 + vh0 + M0 (v ))X (v ) + O(r). (3.71)
It is now clear that the value of X (0) depends on the functions b0 , M, and
h, which in turn depend on the initial data at t = ti , the dynamical variable
v , and the evolution function A(r, v ). Therefore, a given set of initial matter
distributions and the dynamical profiles, including the energy distribution of
the shells, completely determine the tangent at the center to the singularity
curve.
3.4 Nature of singularities 79
F 2 FF F2
K = 12 − 32 + 48 . (3.72)
R4 R2 R5 R R6
g ik R ,i R ,k = 0. (3.73)
F
=1 (3.75)
RN −3
in a general N -dimensional spacetime, which gives the boundary of the
trapped surface region of the spacetime. For the usual spacetime with N = 4,
this becomes F = R. If the neighborhood of the center gets trapped prior
to the epoch of singularity, then it is covered and a blackhole results, other-
wise it could be naked when non-spacelike future directed trajectories escape
from it.
Therefore, the important point is to determine if there are any future
directed non-spacelike paths emerging from the singularity. To investigate
this, and to examine the nature of the central singularity at R = 0, r = 0,
consider the equation for the outgoing radial null geodesics that is given by
dt
= eψ−ν . (3.76)
dr
It would be desirable to examine if there are any families of future directed
null geodesics emerging from the singularity, thus causing a naked singularity
phase as the collapse endstate. The singularity occurs at v (ts (r), r) = 0,
that is, at R(ts (r), r) = 0. Therefore, if there are any future directed null
geodesics terminating in the past at the singularity, R → 0 as t → ts along
these curves. Now, writing (3.76) in terms of the variables (u = rα , R),
dR 1 −(α−1) Ṙ ψ−ν
= r R 1 + e , (3.77)
du α R
or equal to four, that is, when F (ts (r), r) > 0 for r = 0, all singularities for
r > 0 are necessarily covered since F/R → ∞, and hence dR/du → −∞.
This is when both the pressures pr and pθ are positive with the energy con-
ditions are satisfied. Therefore, in such a case only the central singularity at
R = 0, r = 0 could be naked.
Hence, the central singularity at r = 0, R = 0 needs to be examined
to determine if it is visible or not, and to investigate if any solutions to
the outgoing null geodesics equation exist that terminate in the past at the
singularity, and that go to faraway observers in the future. The conditions
under which this can happen are to be determined. Also, note that since
the singularity curve and the evolution functions are regular, the limit of the
functions H , G, and F/R at r → 0, t → ts0 can be calculated. From (3.46),
as A(r, v ) is a well-defined function, G(ts0 , 0) = 1. Also, from (3.51) at this
point, H ≈ r2 /v . Calculating this limit on the t = ts0 plane from (3.71), at
the point (ts0 , 0), H = 0. Hence, from (3.31), F/R = 0 in this limit.
Let x0 now be the tangent to the outgoing null geodesics in the (R, u)
plane at the central singularity
R dR
To find out whether the null geodesic equation admits any solution of x0
that is positive and finite at the central singularity, the values of H , G, and
F/R at (ts0 , 0) can be used in (3.78). Also, (3.71) can be used to get the
value of v v 1/2 on the v = 0 surface at r = 0 (that is, on the point (ts0 , 0)).
Therefore, solving (3.78),
3
M0 (0)X (0),
3/2
x0 = (3.80)
2
and the equation of the radial null geodesic emerging from the singularity is
given by R = x0 u in the (R, u) plane, or in (t, r) coordinates it is given by
It now follows that if X (0) > 0, then x0 > 0, and a radially outgoing
null geodesic emerges from the singularity, giving rise to a naked central
singularity. However, if X (0) < 0, then a blackhole solution exists, as there
will be no such trajectories (see Fig. 3.2). If X (0) = 0, then the next higher
order non-zero term in the singularity curve equation will have to be taken
into account, and a similar analysis has to be carried out by choosing a
different value of α.
It can be seen that the above is both a necessary and sufficient condition
for an outgoing radial null geodesic emerging from the singularity to exist.
82 Spherical collapse
Initial
surface
0(r), p0(r), pr0(r)
Fig. 3.2 Given the regular matter initial data, the collapse can evolve either to
a blackhole or a naked singularity final state, depending on the choice of the rest
of the free functions, such as the velocities E(r) of the collapsing shells and the
allowed dynamical evolutions as given by the Einstein equations.
Assume that such a geodesic does exist and that in the (R, u) plane it is
given by the equation R = x0 u, with x0 > 0. Then, at the central singularity
(R = 0, u = 0), the tangent to such a geodesic must be x0 . Also, this tangent
must be the root of the equation
⎛ ⎞
1/2
dR 3 R vv − F/R
⎝ √ 1√
du
−
5 u
+ 1/2 √ ⎠ = 0 (3.82)
(R/u) G G+ H
Along this curve t → ts0 as r → 0. And, as X (0) > 0, this curve is outgoing
in the sense that t increases as r increases along the curve. The quantity
(−gtt dt2 + grr dr2 ) along this curve can be calculated in the vicinity of the
central singularity. Using (3.71), (3.51), and (3.46), for this curve,
3
5 3
2/3
= M0 (0)X (0) r −dr2 + dr2 = 0. (3.85)
3 2
That is, in the vicinity of the central singularity the curve considered above
is null. Therefore, for any given positive value of the tangent to the singu-
larity curve at the central singularity, it is always possible to find a null and
outgoing curve terminating in the past at the central singularity, making the
singularity naked.
Note that basically it is the geometry of the trapped surfaces and the
apparent horizon that decides the visibility, or otherwise, of the spacetime
singularity. Different kinds of collapse evolutions lead to different trapped
surface configurations, thus leading to visibility, or otherwise, of the final
singularity. For example, while in a homogeneous dust collapse, the trapped
surfaces and apparent horizon form early enough to cover the singularity;
when inhomogeneities are included, the trapping is naturally delayed so as
to allow the singularity to be visible.
Hence, some remarks on the nature of the apparent horizon and its relation
to the visibility, or otherwise, of the singularity are made below. To find the
equation of the apparent horizon near the central singularity, let the time
corresponding to a shell labeled by r entering the apparent horizon, in terms
of the variable v , be vah (r). Then, from (3.75), it can easily be seen that
vah (r) is the root of the equation
r2 M(r, v ) − v = 0. (3.86)
Now, using (3.52), the equation for the apparent horizon in the (t, r) plane
can be written as
vah (r)
v 1/2 dv
tah (r) = ts (r) − . (3.87)
0 e(rA+ν) vb0 + e2ν (vh + M)
It is obvious that the necessary condition for the existence of a locally naked
singularity is that the apparent horizon curve must be an increasing function
84 Spherical collapse
Event horizon
Outgoing null
geodesics
Initial surface
Fig. 3.3 Local versus global visibility. If the singularity is only locally visible,
the light rays come out, but they then fall back again at the center without
coming out of the cloud. On the other hand, for a globally visible singularity,
the outgoing rays reach the boundary of the cloud and can reach faraway external
observers.
of the star, thus not allowing the ultra-dense regions to be visible to faraway
external observers (see Fig. 3.3).
This issue has still not been studied for the general class of models such
as the ones considered here that may be of interest. However, for the case
of dust collapse models, this has been studied in some detail (Joshi and
Dwivedi, 1993a), and it is shown that, whenever the singularity is locally
naked, one can always choose the classes of the mass and energy functions
suitably, as one moves away from the center for the larger values of the radial
coordinate r, in such a manner that the singularity becomes globally visible.
The point is, while the local visibility of the central singularity is basically
decided by the conditions near the center, the global visibility depends on the
overall behavior of these functions within the matter cloud, away from the
center. This behavior can still be freely chosen at larger values of r. In other
words, for the dust collapse models, once the singularity is locally visible,
there are always classes of mass and energy functions which can be chosen
in order to make it globally visible. Another important related point is, as
such, there is no scale in the problem, and the size of the collapsing cloud
could be quite large. In such a case, even if the singularity is only locally
visible, it can still be seen for a long enough time by the observers. Therefore,
in principle, a locally naked singularity is also as serious a violation of the
cosmic censorship as a globally visible singularity, and there may not be a
qualitative difference in the two cases in many situations of physical interest
(Penrose, 1979).
86 Spherical collapse
fundamental forms are matched. That is, the metric coefficients and the
extrinsic curvature are matched at the boundary of the cloud.
Whereas the procedures used below are standard, the particular case
treated here will be described in some detail so as to give the exact pic-
ture of the overall collapse scenario emerging. Note that since the matching
is for the second fundamental form Kij , there is no surface stress energy, or
surface tension at the boundary (see for example, Mazur and Mottola, 2004).
The metric just inside Σ is
which describes the geometry of the collapsing cloud. The metric in the
exterior of Σ is given by
2M (rv , V )
ds2+ =− 1− dV 2 − 2dV drv + rv2 dΩ2 , (3.89)
rv
where V is the retarded outgoing null coordinate and rv is the Vaidya radius.
Matching the area radius at the boundary results in
R(rb , t) = rv (V ). (3.90)
Then, on the hypersurface Σ, the interior and exterior metrics are given by
and
2M (rv , V ) drv
ds2Σ+ =− 1− +2 dV 2 + rv2 dΩ2 . (3.92)
rv dV
1
nV+ = − , (3.95)
2M (rv , V ) drv
1− +2
rv dV
2M (rv , V ) drv
1− +
rv dV
nr+v = . (3.96)
2M (rv , V ) drv
1− +2
rv dV
1
Kab = Ln gab . (3.97)
2
That is, the second fundamental form is the Lie derivative of the metric
with respect to the normal vector n. The above equation is equivalent to
1
Kab = gab,c nc + gcb nc,a + gac nc,b . (3.98)
2
−
Now, setting Kθθ − Kθθ
+
Σ
= 0 on the hypersurface Σ,
2M (rv , V ) drv
1− +
rv dV
RR e−ψ = rv . (3.99)
2M (rv , V ) drv
1− +2
rv dV
Simplifying the above equation using (3.93) and the Einstein equations,
F Re−ν √
M (rv , V ),rv = + √ H ,t + Re2ν ν e−ψ. (3.102)
2R G
90 Spherical collapse
r=0
Generalized
Vaidya
metric
Outgoing
Interior radiation
geometry
Fig. 3.4 For a general matter field, the interior collapsing metric is matched to
an exterior that is a generalized Vaidya geometry.
Any generalized Vaidya mass function M (v, rv ) that satisfies (3.102) will
then give a unique exterior spacetime with required equations of motion
given by other matching conditions, (3.100), (3.101), and (3.90).
To see that the set of all such functions M (v, rv ) is non-empty, the
examples of a charged Vaidya spacetime M = M (V ) + Q(V )/rv , and the
anisotropic de Sitter spacetime M = M (rv ) are two different solutions of
(3.102) (see for example, Joshi and Dwivedi, 1999; Giambo, 2005). This
gives two unique exterior spacetimes, both of which are subclasses of the
generalized Vaidya metric (see Fig. 3.4).
It is of course also possible to treat this problem from the perspective of
only pure general relativity, and the theoretical cosmic censorship aspect. In
this case, there is no need to cut off the cloud and match it at the boundary.
Suitable fall off conditions can just be imposed for the collapsing matter,
so that faraway the metric becomes Minkowskian at the spatial infinity (see
for example, Choptuik, 1993, for the case of numerical scalar field collapse
models).
Space Space
(a) (b)
Fig. 3.5 The nature of the initial density profile affects the final outcome of
collapse. For a constant density (a) a blackhole develops, but for a profile slightly
higher at the center of the cloud (b) the final state is a naked singularity.
R2
ds2 = −dt2 + dr2 + R2 (dθ2 + sin2 θ dφ2 ). (3.103)
1+f
F
= (t, r) = , (3.105)
R2 R
where is the energy density, R is the area radius of the cloud, which is a
function of both t and r, and
F
Ṙ2 = + f. (3.106)
R
Here, the dot and the prime denote partial derivatives with respect to the
parameters t and r respectively.
As the continual gravitational collapse situation is of concern here, it is
required that Ṙ(t, r) < 0. The quantities F and f are arbitrary functions
of r, as obtained from the integration of the Einstein equations, and, as
seen from the above equations, these have the interpretation of the mass and
velocity functions respectively. The quantity 4πR2 (t, r) gives the proper area
of the mass shells and the area of such a shell at r = const. goes to zero when
R(t, r) = 0. Integration of (3.106) gives
The function t0 (r) is a constant of integration. So, there are in all three
arbitrary functions of r, as given by f (r), F (r), and t0 (r). However, the
94 Spherical collapse
R(0, r) = r. (3.111)
The time t = t0 (r) corresponds here to the value R = 0 where the area of
the shell of matter at a constant value of the coordinate r vanishes. It follows
that the singularity curve t = t0 (r) corresponds to the time when the matter
shells meet the real physical singularity, which corresponds to the entire
cloud shrinking to a zero radius. Therefore, the range of the coordinates is
given by
0 ≤ r < ∞, −∞ < t < t0 (r). (3.113)
The case considered by Oppenheimer and Snyder (1939), which discussed the
homogeneous dust cloud collapse, is a special subcase of the above general
dust solution, corresponding to = (t) only throughout, and with a special
value of the velocity function as given by f (r) = 0. The above comoving
coordinates, tailored to the collapsing matter, which have a somewhat spe-
cial physical significance, have been used here, and these provide a physical
insight into the collapse processes. In the case of a homogeneous density, the
interior of the cloud has a special, more simple form, which is the same as
the time reversed case of the Friedman solution in cosmology.
It is seen, in general, that when the matter density is allowed to be inho-
mogeneous and can change as a function of the coordinate radius r, then,
unlike the collapsing Friedmann case (where the physical singularity R = 0
occurs at a constant epoch of time, say, at t = t0 ), the singular epoch now is
a function of r, as a result of the inhomogeneity in the matter distribution.
The Friedmann case could be recovered from the above equations by setting,
for example, t0 (r) = t0 (r) = 0. This corresponds to the case of a simulta-
neous singularity. It follows that, for an inhomogeneous matter distribution,
different shells of matter arrive at the singularity at different times as given
by ts (r), and a singularity curve rather than a constant singular epoch is
obtained, as was the case for a constant density distribution.
The function f (r) above can be used to classify the spacetime as bound,
marginally bound, or unbound, depending on the range of its values, which are
where (0, r) = ρ(r) is the initial density distribution, from which the grav-
itational collapse of the cloud develops. For physical reasonableness, the
weak energy condition would be assumed throughout the spacetime, that is,
Tij V i V j ≥ 0 for all non-spacelike vectors V i . This implies that the energy
density is positive everywhere ( ≥ 0), including the region near the center
of the cloud r = 0.
The partial derivatives of R, as given by R and Ṙ, play an important role
in the analysis, and, from the above equations,
1/2
3 1
R = rα−1 (η − β )X + Θ − η − β X 3/2 G(−P X ) P+
2 X
≡ rα−1 H (X, r), (3.116)
1/2
Λ1/2 1 3
Ṙ = 2
−βX 2
+ P + Θ − η − β X 3/2 G(−P X )
2rX X 2
−N (X, r)
≡ , (3.117)
r
where, using the notation given by Joshi and Dwivedi (1993a),
F f
X = (R/rα ), η = η (r) = r , β = β (r) = r , p = p(r) = rf /F,
F f
(3.118)
F
P = prα−1 , Λ= , (3.119)
rα
and
√
t0 Λ 1+β−η (η − 32 β )G(−p)
Θ ≡ α−1 = + . (3.120)
r (1 + p)1/2 r3(α−1)/2 r3(α−1)/2
The function β (r) is defined to be zero when f is constant and zero. The
factor rα has been introduced here for the sake of convenience in examin-
ing the structure of the naked singularity. The exact value of the positive
constant α ≥ 1 is to be determined and chosen later, and it will depend
on the different models of the spacetime that allow naked singularities. The
functions H (X, r) and N (X, r) are defined by (3.117) and (3.118).
96 Spherical collapse
all the matter shells reduce to a zero size, and therefore this is a genuine
physical singularity. Hence, the occurrence of such shell-focusing singulari-
ties at the center of the collapsing dust cloud are investigated here, and their
nature and structure for the TBL spacetimes are examined. It is clear that
a shell-focusing singularity occurring at r > 0, R = 0 is totally spacelike and
therefore the discussion would be confined to the singularity at r = 0.
The equation for the shell-focusing singularity R(t0 , r0 ) = 0 is as given
earlier, and occurs at r = r0 at the coordinate time t = t0 . The singularity
is called a central singularity if it occurs at r = 0. This central shell-focusing
singularity can be naked, though gravitationally weak (Newman, 1986), for
the class of TBL models considered by Christodoulou (1984), for which the
energy density, which is assumed to be positive everywhere and taken to be
non-zero at r = 0, and the metric functions are even and smooth functions of
t and r. Translated in terms of the parameters defined above, this corresponds
to the case for which η (0) = 3, β (0) = 2, and p(r) is an even smooth function
of r. In terms of the functions F (r) and f (r), this amounts to the conditions
It was, however, pointed out by Waugh and Lake (1988, 1989) and Ori and
Piran (1987, 1990) that this class of gravitationally weak naked singularities
excludes the self-similar TBL models, where the singularity was shown to
be gravitationally strong along the Cauchy horizon, which is a null geodesic
coming out of the singularity. It was further pointed out that the strong
curvature naked singularity is not necessarily confined to the self-similar
models only, but that there are many classes of TBL models that are non-self-
similar, and the naked singularity occurring there is gravitationally strong
(Dwivedi and Joshi, 1992). In these classes, the collapse terminates in a naked
singularity that is gravitationally strong for a wide range of TBL models that
are non-self-similar in general, and which include all the self-similar models as
a special subclass. In the notation used here, these models are characterized
by the conditions η (0) = 1, with F (r) and f (r) being analytic at r = 0.
Rather general differentiability conditions can be required only on the
functions F (r) and f (r) in that they are assumed to be at least C 1 at the
center r = 0, ∞ > η (0) > 0, and β (0) is finite. In order to ensure the
metric to be C 2 , which is the basic requirement in general relativity, these
differentiability conditions imposed here on the functions f and F are suffi-
cient. Frequently, much stronger conditions are imposed on these functions
and the metric, requiring them to be smooth and analytic. However, these
are considered here to be too strong in general. Anyway, the models with
stronger conditions become special cases of the above. Furthermore, from the
physical reasonableness, one would require F (0) = 0 (otherwise there will
be a massive singularity present already at r = 0), which implies η (0) > 0.
98 Spherical collapse
where dΩ2 = dθ2 + sin2 θ dφ2 . The value of the Schwarzschild radial coordi-
nate is rs = R(t, rc ) at the boundary r = rc . Also, m(rc ) = M , where M is the
total Schwarzschild mass enclosed within the dust ball of coordinate radius
r = rc .
The apparent horizon, which indicates the boundary of all trapped surfaces
that define the region that allows no light to escape, forming in the spacetime
during the gravitational collapse in the interior dust ball, lies at R = F (r).
The corresponding time t = tah (r) for the apparent horizon is given by
r3/2 G(−p)
t = tah (r) = √ − F G(−f ). (3.124)
F
The emissions from the shell-focusing singularity R(t0 , r0 ) = 0 for all r0 > 0
would lie in the region above t = tah , that is, t0 > tah for all r0 > 0, with t0
being the time when the singularity at r = r0 occurs. Hence, all radiations
would be future trapped from the shell-focusing singularities at r > 0. At
r = 0 however, t(0) = tah (0) and the singularity could be at least locally
visible. Any light ray terminating at this singularity in the past goes to the
future infinity if it reaches the surface of the cloud r = rc earlier than the
apparent horizon at r = rc . In such a case the singularity would be globally
naked. The nature of this central spacetime singularity has to be examined
in terms of its visibility or otherwise.
3.6 Dust collapse 99
As seen from the Einstein equations above in the dust case, once the
mass function is specified at the initial epoch, the energy function f (r) fully
specifies the velocity distribution of the in-falling shells. Also, the energy
condition then implies that F ≥ 0 and the collapsing matter cloud condition
implies that Ṙ < 0. For dust clouds it follows from the equations of motion
that once Ṙ < 0 at the initial epoch from where the collapse commences,
then, at all epochs, the same condition holds, and thus there is a continual
collapse without any reversal until the shell-focusing singularity at R = 0
is reached. In other words, there is no bounce possible in the dust collapse
once the collapse has initiated with Ṙ < 0.
The scaling independence of the comoving coordinate r can again be used
to give
where
This means the coordinate r has been scaled in such a way that, at the
initial epoch, R = r, and at the singularity, R = 0. Therefore, R = 0 both
at the regular center r = 0 of the cloud, and at the spacetime singularity
where all matter shells collapse to a zero physical radius. The regular center
is then distinguished from the singularity by suitable behavior of the mass
function F (r), so that the density remains finite and regular there at all
times until the singular epoch. The introduction of the parameter v then
allows the spacetime singularity to be distinguished from the regular center,
with v = 1 all through the initial epoch, including the center r = 0, and
which then monotonically decreases with time as the collapse progresses to
the value v = 0 at the singularity R = 0.
From the equations of motion, it is evident that to have a regular solution
over all space at the initial epoch, the two free functions F (r) and f (r) must
have the following forms
where M(r) and b(r) are at least C 1 -functions of r for r = 0, and at least
C 2 -functions for r > 0. This is dictated by the condition that the density
and energy distributions must be regular at the initial epoch and should not
blow up. This is because if the mass function F did not go as a power of at
least r3 closer to the origin, then, as implied from the equations of motion,
the density will be singular at the origin r = 0 as it will diverge there. This
cannot be accepted as regular initial data for the collapse. Similarly, (3.107)
implies that f (r) is determined once the velocity profile is specified, and vice
100 Spherical collapse
versa, for a given initial density distribution. Since the center of the cloud
is taken to be at rest in any spherically symmetric distribution, the leading
term in the energy profile must be at least r or higher. Then again, (3.107)
implies the behavior for f (r) as above.
Here, the negative sign implies that v̇ < 0, that is, the matter cloud is
collapsing. Integrating the above equation with respect to v ,
1
v 1/2 dv
t(v, r) = . (3.129)
v M(r) + vb(r)
3.6 Dust collapse 101
X2 (v ) X3 (v )
t(v, r) = t(v, 0) + rX (v ) + r2 + r3 + ···, (3.130)
2 6
where the function X (v ) is given by
1 1
v 1/2 (M1 + vb1 )dv
X (v ) = − , (3.131)
2 v (M0 + vb0 )3/2
where
Therefore, the time taken for the central shell to reach the singularity is
given by
1
v 1/2 dv
ts0 = √ . (3.133)
0 M0 + vb0
From the above equation it is clear that for ts0 to be defined,
In other words, the continual collapse condition implies the positivity of the
above term. Hence, the time taken for other shells to reach the singularity
can be given by the expansion
Also, from the equation for v̇ above and (3.131), for small values of r along
constant v surfaces,
v 1/2 v = (M0 + vb0 ) (X (v ) + rX2 (v ) + · · · ) . (3.136)
Now, it can easily be seen that the value of X (0) depends on, and is com-
pletely characterized by, the functions M (r) and b(r), which in turn specify
fully the initial mass and energy distributions for the collapsing matter. Spec-
ifying these functions is equivalent to specifying the regular initial data for
the collapse on the initial surface t = 0. In other words, a given set of density
and energy distributions completely determines the slope to the singularity
curve at the origin, which is the central singularity. Also, it is evident that,
given any one of these two profiles, the other one can always be chosen in
such a manner that the quantity X (0) will be either positive or negative.
102 Spherical collapse
dt
= eψ. (3.138)
dr
Choosing α = 5/3, and using (3.129) together with the collapse condition
Ṙ < 0,
dR 3 R v v 1/2 1 − FR
= +
1/2
. (3.140)
du 5 u R e−ψ R e−ψ R +
Ṙ
u
3.6 Dust collapse 103
If the null geodesics terminate at the singularity in the past with a definite
tangent, then at the singularity dR/du > 0 in the (u, R) plane, and this
must have a finite value.
In the case under consideration, all singularities for r > 0 are covered
since F/R → ∞ in the limit of the approach to the singularity, and hence
dR/du → −∞. Therefore only the singularity at the central shell could be
naked.
In order to see the possible emergence of null geodesics from the central
singularity, (3.141) needs to be analyzed. The limits of the concerned func-
tions in (3.141) are calculated at the central singularity. In the TBL case, in
the limit of t → ts , r → 0, e−ψ R → 1. Also, from the equation for v̇ above
and (3.137), in this limit Ṙ → 0. It then follows in general, from the Einstein
equations discussed above, that the term F/R goes to zero in this limit.
It would be interesting to find out when there will be future directed null
geodesics emerging from the central singularity with a well-defined and defi-
nite positive tangent in the (t, r) or (R, u) plane, thus making the singularity
visible. The tangent to the null geodesic at the singularity can be defined as
R dR
Using (3.141) and (3.137), together with the required limits as above,
3
M0 X (0).
3/2
x0 = (3.142)
2
The necessary and sufficient conditions for a naked singularity to exist, that
is, for the null geodesics with a well-defined tangent to come out from the
central singularity can be deduced. Suppose X (0) > 0, then, from (3.142),
x0 > 0 always; then in the (R, u) plane the equation for the null geodesic
that comes out from the singularity is given by
R = x0 u. (3.143)
In other words, (3.143) is a solution of the null geodesic equation in the limit
of the central singularity. Therefore, given X (0) > 0, a solution of radially
outgoing null geodesics emerging from the singularity can be constructed.
This makes the central singularity visible. In the (t, r) plane, the above null
geodesics near the singularity will be given as
It follows that X (0) > 0 implies x0 > 0 and radially outgoing null geodesics
emerge from the singularity, giving rise to the central naked singularity.
104 Spherical collapse
On the other hand, if X (0) < 0, then the singularity curve is a decreas-
ing function of r. Hence, the region around the center gets singular before
the central shell, and after that it is no longer in the spacetime. Now, if
there would have been any outgoing null geodesic from the central singular-
ity, it must then go to a singular region or outside the spacetime, which is
impossible. Hence, when X (0) < 0, there is always a blackhole solution.
If X (0) = 0 then the next higher order non-zero term in the singularity
curve equation will have to be taken into account, and a similar analysis by
choosing a different value of α in (3.140) will have to be carried out.
It has thus been shown above that X (0) > 0 is the necessary and sufficient
condition for null geodesics to emerge from the central singularity with a
definite positive tangent. It should be noted, however, that in general the
dependence of R on r along the outgoing null geodesics from the singularity
does not necessarily have to be of a power law form. However, in order
to satisfy the regularity and physical relevance, examining the trajectories
that come out with a regular and well-defined tangent is physically more
appealing, which is the case examined here.
Shown below for completeness is that if null geodesics of any form come
out at all, then those with a definite tangent also must emerge from the
central singularity. Towards this end, consider the equation of the apparent
horizon. From the equations F = R and (3.130), this is given by
X2 (0)
tah (r) = ts0 + rX (0) + r2 + · · · − O(r3 ). (3.145)
2
that non-dust forms of matter, and the effects of pressures must be suitably
taken into account. In such a case, as pointed out above, it is seen that even
if only homogeneous initial density profiles (with non-zero initial pressures)
are considered, then the pressure by itself can also cause sufficient distortions
in the formation of the apparent horizon so as to cause a naked singularity
as the endstate of the collapse, rather than a blackhole. However, if the
equation of state is homogeneous, together with the initial data having a
homogeneous density profile, then no naked singularity will appear.
Note that the formalism given here brings out the role of the initial
data in causing the blackhole and naked singularity endstates for the four-
dimensional dust collapse in a clear and transparent manner. To be specific,
it is seen, using (3.131), that given any initial density distribution for the
cloud, a suitable energy profile can be chosen so that the evolution could end
in either a blackhole or a naked singularity, depending on the choice made.
In other words, there is a non-zero measure of energy distributions that will
take the given density profile to a blackhole. The same holds for a naked
singularity to evolve from the same initial density. The converse is also true,
namely that given any initial energy distribution, the density profiles that
give rise to either of these endstates can be chosen.
condition for the existence of positive real roots into the actual constraints
on the initial density distribution in the cloud.
This issue is investigated here, and it would be expected that the degree
of inhomogeneity in the matter distribution plays a role in determining the
final fate of the collapse. This is seen by working out the explicit conditions
for collapse to an endstate that is either a blackhole or a naked singularity,
depending on the initial conditions chosen. It turns out that for the class
under consideration here, these outcomes are characterized in terms of the
existence of real positive roots of a quartic equation. This enables the black-
hole or naked singularity configuration as the endstate of the gravitational
collapse to be related in terms of the initial density distribution ρ(r) and
radius rb of the massive body.
As used earlier, the quantity R(t1 , r1 ) denotes the physical radius of a shell
of collapsing matter at a coordinate radius r1 and on the time slice t = t1 .
The quantities F and f are arbitrary functions of r. In further discussion,
the class of solutions is restricted to f (r) ≡ 0, which comprise the marginally
bound TBL models. Similar considerations can be developed for the models
with f > 0 or f < 0. As the collapsing cloud is of interest here, Ṙ(t, r) < 0.
The epoch R = 0 denotes a physical singularity where the spherical shells of
matter collapse to a zero radius, and where the density blows up to infinity.
The time t = t0 (r) corresponds to the value R = 0 where the area of the
shell of matter at a constant value of coordinate r vanishes. This singularity
curve t = t0 (r) corresponds to the time when the matter shells meet the
physical singularity. In the case of a finite cloud of dust, there will be a cut
off at r = rb , where the metric is matched smoothly with a Schwarzschild
exterior, as discussed earlier.
In TBL models, the freedom to relabel the dust shells arbitrarily exists. It
is given by r → g (r) for any shell with r = const. on any t = const. epoch.
Therefore, at any constant time surface, say at t = t0 , R(r, t0 ) can be chosen
to be an arbitrary function of r. This arbitrariness reflects essentially the
freedom in the choice of units. The choice of scaling at t = 0 as given by
R(r, 0) = r is made for the sake of convenience in the calculation. With this
scaling, the Ṙ equation (with f = 0), can be integrated to obtain
3
R3/2 (r, t) = r3/2 − F (r)t, (3.146)
2
and the energy equation is given by
4/3
(r, t) = , (3.147)
2 G(r) 2 G (r)
t− t−
3 H (r) 3 H (r)
where G(r) = r3/2 , G (r) = (3/2)r1/2 , and H (r) = F (r).
108 Spherical collapse
Now, write F (r) ≡ rλ(r), and assume λ(0) ≡ λ0 = 0 and is finite, which
is the class of models considered by Dwivedi and Joshi (1992). This means
that, near the origin, F (r) goes as r in the present scaling, and the density
at the center behaves with time as (0, t) = 4/3t2 . This is a general class
of models that includes all self-similar solutions, as well as a wide range of
non-self-similar models, which are found to be quite adequate for the purpose
of the present investigation. The central density becomes singular at t = 0,
and the singularity is interpreted as having arisen from the evolution of the
dust collapse, which had a finite density distribution in the past on an earlier
non-singular initial epoch.
To check if the singularity could be naked, the possibility that future
directed null geodesics would come out of the singularity at t = 0, r = 0
needs to be examined. The equations of the outgoing radial null geodesics in
the spacetime, with k as the affine parameter, can be written as
dK t
+ Ṙ K r K t = 0, (3.148)
dk
dt Kt
= r = R , (3.149)
dr K
where K t = dt/dk and K r = dr/dk are tangents to the outgoing null geo-
desics. The partial derivatives R and Ṙ that occur in the above equations
can be worked out from the solution given above, and these are most suitably
written as
1−η t
R = ηP − √ + η Ṙ, (3.150)
λ r
λ 1−η t
Ṙ = √ +η , (3.151)
2rP 2 λ r
where
rF
R(r, t) = rP (r, t), η = η (r) = . (3.152)
F
The functions η (r) and P (r) are introduced because they have a well-defined
limit in the approach to the singularity.
If the outgoing null geodesics terminate in the past with a definite tangent
at the singularity (in which case the singularity would be naked), then using
(3.150) and l’Hospital’s rule,
t dt
X0 = lim = lim = lim R , (3.153)
t→0, r→0 r t→0, r→0 dr t=0, r=0
3.6 Dust collapse 109
2 2P 3/2
X− √ =− √ . (3.154)
3 λ 3 λ
2
X0 < √ . (3.155)
3 λ0
V (X0 ) = 0, (3.156)
where !
λ0
V (X ) ≡ Q + X − X. (3.157)
Q
In order to be the past end point of the outgoing null geodesics, at least
one real positive value of X0 must satisfy (3.156), subject to the constraint√
implied by (3.155), as stated above, and which is given by X0 < 2/3 λ0 .
In general, it was shown by Dwivedi and Joshi (1992) that if the equation
V (X0 ) = 0 has a real positive root, the singularity would be naked. Whenever
this is not realized, the evolution will lead to a blackhole. Such a singularity
could be either locally or globally naked depending on the global features of
the function λ(r).
The sense in which the terms naked singularity and blackhole are used
needs to be clarified. When there are no positive real roots to (3.156), the
central singularity is not naked, because there are no outgoing future directed
null geodesics from the singularity in this case. Furthermore, as discussed
earlier, the shell-focusing singularity R = 0 for r > 0 is always covered.
Hence, in the absence of positive real roots, the collapse will always lead to a
blackhole. On the other hand, if there are positive real roots, it follows that
the singularity is at least locally naked. The global visibility aspect of such
a locally naked singularity will be discussed in the next subsection. Such a
locally naked singularity would also be globally naked when the outgoing
trajectories could reach arbitrarily large values of r, that is, the signals reach
faraway observers. Otherwise, the collapse outcome would still be a blackhole
110 Spherical collapse
when these trajectories that emerge from the singularity fall back to the
singularity again without emerging from the horizon. This is a violation of
weak censorship only. The occurrence of either of these situations will really
depend on the nature of the function λ(r). The conditions under which this
locally naked singularity could be globally naked have been discussed, for
instance, by Joshi and Dwivedi (1993a). The occurrence of positive real roots
implies the violation of strong cosmic censorship, though not necessarily
of weak cosmic censorship. In other words, a blackhole and locally naked
singularity are not mutually exclusive alternatives. It can also be shown that
whenever there is a positive real root as above, a family of outgoing null
geodesics always terminates at the singularity in the past.
The condition for the occurrence of a naked singularity is now examined
in some detail. The condition V (X0 ) = 0 can be written as
2
Y3 Y − − α(Y − 2)3 = 0, (3.158)
3
√ 3/2
where Y = λ0 X0 , and α = λ0 /12 have been set. Note that F (r), and
hence λ0 , are always necessarily positive. Using standard results, it can be
shown that this quartic equation has positive real roots if and only if α > α1
or α < α2 , where
26 √
α1 = + 5 3 ≈ 17.3269
3
and
26 √
α2 = − 5 3 ≈ 6.4126 × 10−3 . (3.159)
3
To derive this condition, first note that if it has a real root, it must be posi-
tive, as negative values of Y do not solve this equation. Writing the general
quartic as ax4 + 4bx3 + 6cx2 + 4dx + e = 0 one defines H = ac − b2 , I =
ae − 4bd + 3c2 , J = ace + 2bcd − ad2 − eb2 − c3 , and ∆ = I 3 − 27J 2 . If
∆ < 0, the quartic has two real and two imaginary roots. If ∆ > 0, all roots
are imaginary unless H < 0 and (a2 I − 12H 2 ) < 0, in which case they are
all real. The application of such a procedure to the quartic in (3.158) leads
to the conditions on α given above.
It follows, however,
√ from earlier discussions, that along any such outgoing
null geodesics λ0 X0 = Y < 2/3. Then, (3.155) implies that the larger
range of α for the existence of roots, that is α > 17.3269, is ruled out in
the sense that no outgoing trajectories can meet the singularity with this
larger value of the tangent X0 . This is seen by writing α as a function of
Y , which shows that if α > α1 , then Y > 2. It thus follows that a naked
singularity arises if and only if α < α1 , or equivalently, if and only if λ0 <
0.1809. Whenever the limiting value λ0 does not satisfy this constraint, the
gravitational collapse of the dust cloud must end in a blackhole. The physical
3.6 Dust collapse 111
interpretation for the quantity λ0 can be obtained from the equation for the
time evolution of the energy density. If the collapse starts at a time −t0 < 0,
and ρc is the initial energy density at the center, then ρc = 4/3t20 . If ρc is
the initial density gradient at the center, then λ0 = 16ρ3c /3ρ2 c . Putting in
the units gives
8πG 16ρ3c
λ0 = . (3.160)
c4 3ρ2
c
and the naked singularity results for the values given by β < 0.0113. The
occurrence of one or the other outcome is governed by the conditions of a
combination of the initial central density and the initial density gradient at
the center. The cosmic censorship hypothesis of Penrose could, in the present
context, be translated to the conjecture that values of β smaller than 0.0113
do not occur in a realistic collapse.
Now, the value of the parameters β or λ0 , given the initial density profile
for the collapsing massive star, need to be calculated. Given the initial central
density, the initial density gradient dρ/dr = ρc can be evaluated at the center
as follows. First, note that the expression for dρ/dr can be written, and it is
√ of r → 0, this always goes to a finite quantity that is
seen that in the limit
proportional to 1/ λ0 . Now, given the initial data in the form of the density
distribution ρ(R) for the body at an initial non-singular epoch of time, in
terms of the physical radius R, integration can give the form of the mass
function F (R). Then, the above provides a functional relationship r(R) that
can be inverted, in principle, to express R in terms of r. The mass profile
can then be written explicitly as F (r), and λ0 is evaluated as the limit of
F (r)/r as r → 0.
This model can be generalized in many ways, for instance by considering
the most general class of functions F (r) and f (r). Also, it has been shown
earlier that the pattern of a transition from the blackhole configuration to
the naked singularity configuration persists in models with more general
equations of state. It is desirable to cast results for these models in terms of
constraints on the initial density distribution. For a discussion on non-self-
similar collapse, see Lake (1991).
1/2
a
3
i 2
L(λ) = (X ) ds, (3.162)
0 i=0
which is a finite quantity. The X i values are the components of the tan-
gent vector to the curve in a parallel propagated tetrad frame along the
curve. Each such incomplete curve defines a boundary point of the space-
time, which is a singularity. In order to know if this is a genuine physical
singularity, one would typically like to associate it with unboundedly growing
spacetime curvatures. If all curvature components and the scalar polynomi-
als formed out of the metric and the Riemann curvature tensor remain finite
and well-behaved in the limit of approach to the singularity along an incom-
plete non-spacelike curve, it may be possible to remove the singularity by
extending the spacetime when the differentiability requirements are lowered
(Clarke, 1986).
There are several ways to formalize such a requirement. For example, a
parallel propagated curvature singularity is the end point of at least one non-
spacelike curve on which the components of the Riemann curvature tensor are
unbounded in a parallel propagated frame. On the other hand, a scalar poly-
nomial singularity has a scalar polynomial in the metric and Riemann tensor
taking an unbounded large value along at least one non-spacelike curve with
this singular end point. This includes the cases such as the Schwarzschild
singularity where the Kretschmann scalar Rijkl Rijkl blows up in the limit
as r → 0. Will genuine curvature singularities occur in general relativity?
The answer, for the case of parallely propagated curvature singularities, is
provided by a theorem of Clarke (see for example, Clarke, 1993). This theo-
rem shows that for a globally hyperbolic spacetime that is C 0− -inextendible,
when the Riemann tensor is not very specialized in the sense of not being
type D and electrovac at the singular end point, then the singularity must be
a parallely propagated curvature singularity. Similarly, a class of physically
114 Spherical collapse
for some fixed constant A along the trajectory in the limit of approach to the
singularity as t → 0. This provides a sufficient condition for all the two-forms
µ(k ), defined along the singular null geodesic, to vanish as the singularity
is approached, and implies a very powerful curvature growth establishing
a strong curvature singularity. For timelike geodesics this will imply that
all the volume forms defined by the Jacobi fields along these trajectories
must vanish in the limit of approach to the singularity, or they must vanish
infinitely many times in this limit.
3.6 Dust collapse 115
F (K t )2 F (K t )2
Ψ ≡ Rab K a K b = = . (3.167)
R2 R R2 R
For radial null geodesics, using l’Hospital’s rule and the TBL equations as
given earlier, and the fact that in the limit of approach to the singularity
r → 0, X → X0 ,
4η0 Λ0
lim k 2 Ψ =
√ 2 , (3.168)
k→0 H0 X 02
2 (3α − η0 ) 1 + f0 − N0
where η and Λ are as defined in (3.118) and (3.119). Hence, it is seen from
the definition of Λ in (3.120) that
So, as discussed above, it can be concluded that the condition for strong
curvature is satisfied along non-spacelike geodesics as well if α = η0 , and if
such families meet the naked singularity in the past.
The Kretschmann scalar Rijkl Rijkl along the geodesics in the TBL
spacetimes goes as
K ∝ r2(η0 −3α). (3.172)
Hence, the singularity is a scalar polynomial singularity as long as α > η0 /3.
3.6 Dust collapse 117
will end in a naked singularity when V (X ) = 0 has two real positive and
two complex roots. Let x1 , x2 (x1 > x2 ) be two such positive roots of this
equation. The equation of geodesics, in the form r = r (x), X = x2 is given
by (Joshi and Dwivedi, 1993a)
(x − x2 )n2
r = r(X ) ≡ r(x) = D f1 (x), (3.173)
(x − x1 )n1
where
Ax + B
f1 (x) = exp − dx . (3.174)
x2 + D1 x + D2
Here, n1 , n2 , A, B, D1 , and D2 are constants given by
√
Λ0 3 √
x + 4
x − x + Λ0 = (x − x1 )(x − x2 )(x2 + D1 x + D2 ), (3.175)
2
3x3 n1 n2 Ax + B
√ √ = − + 2 , (3.176)
Λ0 3
x + 2 x − x + Λ0
4 x − x1 x − x2 x + D1 x + D2
and D is the constant that labels the different geodesics. The constants n1 ,
and n2 are positive. For the case Λ0 = 7/17,
The event horizon is represented by the geodesic for which Xc = Λ(rc ). Since
it is outgoing, dR/d(rα ) is positive at r = 0 and ejected into the region
R > F where dR/dr is positive. Therefore, all the geodesics that reach
r = rc , where the TBL metric is matched with the Schwarzschild exterior,
with Xc > Λ(rc ) would escape to infinity, while others would become ingoing.
It follows that the geodesics that reach future infinity with their past end
point at the singularity are given by (3.180) with Xc > Λc = Λrc . Hence,
in the case when a family of geodesics terminates at the singularity with a
tangent X = X0 and η (r)Λ(r) < αX0 , for rc ≥ r > 0, the singularity would
be globally visible, as there would always be some geodesics that escape to
infinity.
Consider the case now when the equation V (X ) = 0 has two positive roots
X1 and X2 (X1 > X2 ). In such a case, families of curves would emerge from
the singularity with the tangent either X1 or X2 . Let η (r)Λ(r) < αX2 for
rc ≥ r > 0, then it ensures that some geodesics would cross the boundary
of the cloud with Xc > Λ(rc ) making the singularity globally naked. The
same holds even when more than two positive roots exist. Therefore, if the
family of geodesics do terminate at the singularity with tangent X0 , then
the condition η (r)Λ(r) < αX0 for rc ≥ r > 0 implies the global visibility of
the singularity.
cosmological term Λ. This will give rise to a vacuum energy density, corre-
sponding to a positive sign of Λ. This would represent a spatially uniform
time independent energy density distribution, and its positive value acts as
a globally repulsive force field. Also, the cold dark matter models with a
substantial component supplied by Λ can fit the current observational data
(Ostriker and Steinhardt, 1995).
Therefore, it is useful and interesting to consider the Einstein equations
with a non-zero Λ, to investigate the gravitational collapse, and also the
implications in cosmology. With such a perspective, the structure of a sin-
gularity in the spherically symmetric dust models with a non-vanishing
cosmological constant is discussed here. Dust models with a cosmological
term are discussed in the literature (see for example, Krasinski, 1997). These
models can be matched with the Schwarzschild–de Sitter spacetime at the
boundary of the cloud (Omer, 1965; Lake and Roeder, 1979; Lake, 2000.)
Here, the TBL collapse models are analyzed with a cosmological constant,
and provide the general solution to the Einstein equations with dust as the
source term.
The weak energy condition, Tij V i V j ≥ 0, for all non-spacelike vectors
i
V , is assumed in the spacetime for the matter to be physically reasonable,
especially in the case of the gravitational collapse. The Einstein equations
are analyzed to check whether there are globally regular solutions, as have
been shown to exist in the case of homogeneous density dust solutions with
a positive cosmological constant (Markovic and Shapiro, 2000). The general
solution to the Einstein equations in the case corresponding to the marginally
bound models in the TBL spacetimes, which is when the energy function f =
0, is also derived, and the condition for avoiding shell-crossings is discussed.
The structure of the singularity can be investigated by studying the outgoing
radial null geodesics near it. Then, it can be examined how the given initial
conditions, in the presence of the cosmological term, affect the final state
results in terms of the formation of a blackhole or a naked singularity. It
can also be seen why the naked singularity turns out to be stable to the
introduction of the cosmological term.
Here, the gravitational collapse of a dust cloud in the presence of a cosmo-
logical term is discussed in some detail, and how it may affect the dynamical
evolution of the collapse is examined. In particular, it is found that there
is a non-trivial detailed dynamical structure, especially for the f < 0 case,
and this provides a strong motivation for the present study, namely how the
introduction of a cosmological term non-trivially changes the structure of
the collapse. In most cases, there is no rebounce, and then the nature of the
singularity forming as the endstate of the collapse is studied. At a more tech-
nical level, the change in the dynamical structure of the collapse, when there
is a non-vanishing Λ term present, is reflected in the change in the nature
and structure of the singularity, which is far from obvious. There are some
3.6 Dust collapse 121
important changes in the collapse dynamics with the introduction of the cos-
mological constant. In particular, the analysis of the roots equation referred
to earlier shows that part of the naked singularity spectrum is covered when
a positive cosmological term is introduced, but it is never removed altogether
irrespective of its value. In this sense, the naked singularity is stable to the
introduction of the cosmological term.
In fact, there has been hope through the last three decades or so that it
may be possible to avoid singularities by the introduction of a cosmological
term with a positive sign (see for example, Hawking and Ellis, 1973, pp. 139,
362), especially in the cosmological case. Hence, it is important to investigate
the collapse dynamics, especially on a large or medium scale in the universe,
in the presence of a cosmological term. In fact, as shown by Deshingkar
et al. (2001), for a certain range of initial data, there is a rebounce possible
due to the presence of the positive cosmological term. Such a possibility will
be especially important in the case of the collapse of either the clusters or
superclusters of galaxies, where the densities may be sufficiently low at the
initial epoch, and such a rebounce can become reality.
Towards the effects of a non-vanishing cosmological term on the final fate
of an inhomogeneous collapsing dust cloud, it is shown that, depending on
the nature of the initial data from which the collapse evolves, and for a
positive value of the cosmological constant, it is possible to have a globally
regular evolution where a bounce develops within the cloud. The initial data
causing such a bounce is characterized in terms of the initial density and
velocity profiles for the collapsing cloud. In the other cases, the result of
the collapse is either the formation of a blackhole, or a naked singularity,
as the endstate of the collapse. It is also found that a positive cosmological
term can cover a part of the singularity spectrum that is visible in the corre-
sponding dust collapse models for the same initial data. The basic set of the
Einstein equations and the regularity conditions for the collapse is discussed,
and then the possibility when there is a non-vanishing cosmological term
present, when an initially collapsing cloud rebounces at a later epoch so that
a singularity does not form, is investigated. Therefore, there is an occurrence
of three phases, the collapse, the reversal, and the subsequent dispersal. A
general solution with a non-zero Λ term for the marginally bound case is
also given.
The model for a self-gravitating, spherically symmetric, inhomogeneous
dust cloud with cosmological constant is given by the metric
R 2
ds2 = −dt2 + dr2 + R2 (dθ2 + sin2 θ dφ2 ), (3.181)
1 + f (r)
with the stress–energy tensor given by Tji = (r, t)δti δjt , with (r, t) being
the energy density of matter. The weak energy condition is assumed for the
matter, which implies that (r, t) ≥ 0. This is equivalent to the strong energy
condition in this case, as the principal pressures are zero.
The Einstein equations in the presence of a cosmological constant can be
written as
F (r) Λ
Ṙ2 = + f (r) + R2 , (3.182)
R 3
F
(t, r) = 2 . (3.183)
R R
The collapse situation with Ṙ < 0 is mainly considered. The two free func-
tions that characterize the dust cloud are those representing the total mass
F (r), and the energy f (r), inside the shell labeled by the comoving coordi-
nate r. The cosmological term Λ can, in principle, be of either sign, however
the recent observations indicated above seem to favor the positive sign. The
above equation for Ṙ can, in principle, be integrated, and after integration, a
constant of integration is obtained. The constant of integration can be fixed
by using the scaling freedom. Here, this is fixed by setting R = r on the
initial hypersurface (t = 0). The two free functions F (r) and f (r) can be
fixed by prescribing the initial density and velocity profiles through (3.183)
and (3.182) respectively. Assume these free functions have the form
where the choice of the first non-vanishing term is made in order to have
regular initial data on the initial surface t = 0 from which the collapse
evolves. From (3.183), it is clear that it is possible to have both shell-crossing
(R (tsc (r), r) = 0), and shell-focusing singularities (R(tsf (r), r) = 0) in these
spacetimes, depending on the dynamics of the shells given by (3.182). Here,
only cases where there are no shell-crossings are considered, as these are gen-
erally not considered to be genuine singularities, and the main interest here
is in studying the nature of the physical singularity corresponding to R = 0
where the matter shells shrink to zero radius. This puts some restrictions on
the initial data.
The metric exterior to the collapsing cloud has the Schwarzschild–de Sitter
form,
ds2 = −gdt2 + g −1 dr2 + r2 (dθ2 + sin2 θ dφ2 ), (3.185)
where g = 1 − (2M/r) − (Λr2 /3). Matching the solutions (3.181) and (3.185)
at the boundary rb of the collapsing dust cloud, 2M = F (rb ) is obtained.
The boundary can be made to bounce from an initially collapsing phase
3.6 Dust collapse 123
V (R, r)
Ṙ2 = . (3.186)
3R
V = 3F + 3f ζr + Λζ 3 r3 , (3.189)
which is a cubic equation with three roots in general. From the theory of
cubic equations, if all the three roots of the equation above are real, then
at least one of them has to be positive and at least one negative. Note that
V (R = 0, r > 0) = 3F > 0. Hence, any regular region between R = 0 and the
first zero of (3.189), that is, ζ1 > 0, always becomes singular during collapse
and so two real positive roots for (3.189) are required for the possibility of
a rebounce. The region between the two real positive roots is a forbidden
region. So, starting on the right side, the collapsing shells bounce back and
then there is continuous expansion. Since one of the real roots has to be
negative and the region between two real positive roots is forbidden, it is not
possible to have oscillating solutions in these spacetimes.
In the case of Λ = 0, it is well-known that a rebounce is not possible, and
the collapse necessarily results in a singularity. The cubic then reduces to a
linear equation and the solution is given as R = Rmax (r) = −F/f . So, only
in the case where f < 0 can Ṙ = 0 for positive R. This corresponds to the
maximum possible physical radius for a given shell, that is, even if an initial
124 Spherical collapse
expansion is the starting point, a given shell will reach the maximum radius
Rmax (r), and then it will recollapse.
In the case when Λ < 0, one, and only one, root is always positive.
The other two roots are negative if 9F 2 < −(4f 3 /Λ), or else they are
complex conjugates. Therefore, any initial configuration becomes singular
in this case. The real positive root in this case gives an upper bound on
the radius R = Rmax (r) of a shell labeled r. This upper bound occurs as
the negative (attractive) contribution from the Λ term keeps on increasing
with an increasing R, while the contribution from gravitational attraction
keeps on decreasing, and so at some point for any value of f , the attraction
due to Λ starts to dominate. Hence, even if there is an initially expanding
configuration, finally there must always be a collapse in this case.
For Λ > 0 and f (r) ≥ 0, it can easily be seen that it is never possible to
have ζ̇ 2 = Ṙ2 /r2 = 0, so a singularity always forms if an initial collapse is
always started from.
For the case when Λ > 0 and f (r) < 0, one root of the cubic above is
always negative. If
4f 3
F2 > − , (3.190)
9Λ
then the other two roots are complex conjugates. So, the singularity always
forms in such a case in an initially collapsing configuration. On the other
hand, if the initial data is such that
4f 3
F2 < − , (3.191)
9Λ
then the other two roots are real and positive. Denote these roots by ζ1
and ζ2 with ζ1 < ζ2 , and the region between these two roots is forbidden.
The entire space of allowed dynamics is given by the two disjoint regions
[0, ζ1 ] and [ζ2 , ∞]. If the initial scale factor ζ0 lies in the first section, i.e. if
ζ0 < ζ1 , then the singularity is always the end point of the collapse. Here, rζ1
represents the upper bound for the physical radius of a shell in this region.
If ζ0 lies in the second section, that is, ζ0 > ζ2 , then there will be a rebounce
from the initial collapsing configuration. After the rebounce, the physical
radius of the shell keeps on increasing forever. There is no upper limit for
the maximum value of ζ in this region and rζ2 gives the lower bound for the
physical radius of a shell, that is, the shell rebounces at R = rζ2 .
From the above discussion, it can be seen that a rebounce is possible only
in the case when Λ > 0 and f < 0, and when the following two conditions
are satisfied:
4f 3
F2 < − (3.192)
9Λ
3.6 Dust collapse 125
and ζ0 > ζ2 . With the scaling used here, the condition (3.192) can be
written as
!
f 1 9F 2 Λ
1>2 − 2 cos arccos − . (3.193)
r Λ 3 4f 3
and
ζ̇ 2 + ξ ζ̈
g ij rij = 6 (3.195)
ζ2
stay finite, as ζ (t, r) > 0 for the regular models. For further details see
Deshingkar et al. (2001) and Arun Madhav et al. (2005).
There are several features here that are worth noting, and which have
an interesting physical significance as far as the dynamics of the collapse is
concerned. These also illustrate the effects that a non-vanishing cosmological
term may have towards determining the final fate of a collapsing cloud of
matter. First, with a negative value of the cosmological term, all the solutions
become closed and a singularity always forms in the future even if one starts
with initial expansion. This is to be expected because such a value will
only contribute in a positive manner to the overall gravitational attraction
of matter, and it just acts as a constant positive energy field helping the
collapse. Next, there is a range of initial data where the collapse necessarily
ends in a singularity however large the positive value of the cosmological
term is. This is contrary to the belief sometimes expressed that a positive
cosmological constant can always cause a rebounce, provided it is sufficiently
large. Finally, it is clear from the above that there can be a rebounce only if
the initial density is sufficiently low for a given positive value of Λ. This is so
because the cosmological term becomes dominant with increasing distances,
and gravity dominates at higher densities. Therefore, as the cloud is more
disperse, with a lesser density but a larger size, the effect of the cosmological
term is greater.
It can therefore be seen that while a bounce and a regular solution occur
for a specific range of the initial data, for the majority of the regular initial
data space, the collapse results in a spacetime singularity where the densities
and curvatures blow up. While for Λ = 0 the structure of this singularity
and when it will be naked or covered is known in detail, the above discussion
gives some understanding of the effects of a non-zero Λ towards the structure
of the singularity forming in the collapse. To see this more clearly for a dust
126 Spherical collapse
collapse, the case f (r) = 0 is now analyzed explicitly in some detail. While
f = 0 has been chosen for the simplicity and clarity of the considerations,
similar behavior would also be expected in other cases. Equation (3.182) can
now be written as
−1/2
F (r) Λ 2
t − tc (r) = ± + R dR, (3.196)
R 3
Using the scaling freedom in this solution, R(t = 0, r) = r can be set. This
determines tc (r) as
m
2 r3/2 ∞
(−1)m (2m − 1)!! Λr 3
tc (r) = tsf (r) = 1+ .
3 F (r) m=1
2m (2m + 1)m! 3F (r)
(3.198)
From the above expressions,
1/2 1/2
F rF 3F + ΛR3/2 r
R = R+ 1− . (3.199)
3F 3F 3F + Λr3/2 R
As seen below, the behavior of collapsing shells near the center would
depend only on the first non-vanishing derivatives of the density and veloc-
ity profiles. Therefore, the local visibility conditions are unaffected by the
boundary conditions such as the initial choices of the mass function and the
actual value of the Λ term. On the other hand, the global behavior of the
trajectories emerging from the singularity can change due to the addition of
a non-zero cosmological term.
To analyze the nature of the central singularity, the method developed by
Joshi and Dwivedi (1993a), which gives a necessary and sufficient condition
for the local visibility of the singularity is followed. The main idea here is to
see if outgoing radial null geodesics in the spacetime, meeting the singularity
in their past with a well-defined real positive tangent vector in a suitable
plane are possible.
The equation for these null geodesics in the spacetime (3.181), for f (r) = 0
is given by
dt
= R . (3.200)
dr
dR R F Λ 2
= 1− + R , (3.201)
du αrα−1 R 3
R dR R
X0 = lim = lim = lim
u→0, R→0 u u→0, R→0 du u→0, R→0 αr α−1
F Λ
1− + R2 ≡ U (X0 , 0), (3.202)
R 3
where λ = F/u and λ0 is the limit as r → 0. If the above has a real positive
root X0 then there will be at least one null geodesic emerging from the
singularity at R = 0, u = 0 with the root X = X0 as a tangent in the (u, R)
plane.
When n < 3, α < 3 and therefore λ0 = 0, and the above equation
reduces to
F
X0 = −
3/2
n , (3.204)
2F0 1 + ΛF0 /3
which always has a real positive root. Apart from an additional Λ term, this
equation is analogous to the corresponding equation obtained for the TBL
models, and reduces to the same case for Λ = 0. Therefore, when either
the first or the second derivative of density is non-zero, the singularity is
always at least locally visible. The addition of the cosmological term changes
only the value of the tangent (X0 ) to the outgoing radial null geodesics
from the singularity, but not the visibility property itself. Therefore, the
corresponding dust naked singularity spectrum is stable to the addition of a
positive cosmological constant.
In this case, as in the earlier studies on TBL models, the smaller root
will be along the apparent horizon direction and a family of geodesics will
come out along this direction. In this case, set α = 3, and the first term in
(3.201) blows up, and the second term goes to zero such that the product is
F0 . For the value n = 3, which corresponds to the critical case in the TBL
models as discussed
above, λ0 = F0 is obtained. Introducing X = F0 x2 and
5/2
ξΛ = F3 /F0 1 + ΛF0 /3, the root equation above becomes
2x4 + x3 − ξΛ x + ξΛ = 0. (3.205)
This equation is similar to the corresponding case in the TBL models, with
a modification in the definition of ξ . From the theory of quartic equations,
√
this equation admits a real positive root for ξΛ < ξcrit = −(26 + 15/ 3)/2.
Therefore, for a given central density F0 , and the inhomogeneity parameter
as given by F3 , the naked dust singularity can be partly covered by a positive
Λ, because there is an additional positive term in the denominator in ξΛ . But
it is interesting to note that, however large, a finite Λ term cannot completely
cover the corresponding visible part of dust models.
In a similar manner, it is easy to see that a negative Λ will open up some
covered part in the dust collapse final state spectrum. As such, the discussion
above does not use the positivity of Λ, so the discussion goes through even
3.7 Equation of state 129
for negative Λ. For all values of n > 3, λ0 = ∞, and there cannot be a real
positive root in these cases. The final singularity is then hidden behind the
event horizon.
It can therefore be seen that the gravitational collapse of a dust cloud
with a non-zero cosmological term can develop into a blackhole, a naked
singularity, and even a globally regular solution, as the final outcome of
a collapse. Each of these outcomes is determined by the choice of initial
parameters, given in terms of the density and velocity profiles of the cloud.
Although for simplicity, this was restricted to the f (r) = 0 case while ana-
lyzing the structure of the singularity, the results can be extended to the
general case. The main aim here has been to examine how the presence of a
non-zero cosmological term affects the dynamical evolution of the collapse.
It has been found that there is a non-trivial detailed dynamical structure,
especially for the f < 0 case, and this provides a strong motivation for
the present study, namely how the introduction of a cosmological term non-
trivially changes the collapse outcomes. In most cases, there is no rebounce,
and then the nature of the singularity forming as the endstate of the col-
lapse is studied. At a more technical level, this is reflected in the change in
the structure of the roots that characterize the outgoing null geodesics, as
discussed above.
In particular, this analysis shows that it is possible to cover a part of the
naked singularity spectrum in the corresponding critical branch of solutions
in the TBL models, when a positive cosmological term is introduced. How-
ever, the naked singularity of dust collapse is never removed, irrespective
of the value of Λ. In this sense, the naked singularity remains stable to the
introduction of a cosmological term.
It can be seen that there are certain important changes in the collapse
dynamics with the introduction of a cosmological term, and the study here
brings this out, allowing the implications of a non-zero Λ towards the final
outcome of gravitational collapse to be understood. These results are of
interest in view of the recent observational claims about a non-vanishing
cosmological constant.
and metric variables. Hence, these functions are fully determined only as a
result of the time development of the system from the initial data, provided
the relation between the density and pressures, that is, a given ‘equation of
state’ is known.
In principle, it is possible to choose these functions freely. For example,
one could specify the matter and velocity profiles at the initial epoch and
also the dynamical evolutions, such as F (v, r) and ν (v, r), only subject to
an energy condition and regularity. Note that v here plays the role of a time
coordinate and this then fully determines the collapse evolution. One can
then calculate the energy density, and the radial and tangential pressures for
the matter. However, in such a case, the resultant ‘equation of state’ could be
quite strange in general. If any equation of state of the form pr = f (ρ) and
pθ = g (ρ) is given, then it is clear from (3.28) and (3.29) that there would be
a constraint on the otherwise arbitrary functions M and A, which specify
the required class, if the solutions of the constraint equations exist. As for the
physics of very high density matter, presently, there is very little idea on what
kind of an equation of state the matter would follow, especially at very high
densities. Closer to the collapse endstates ultra-high energy densities and
pressures are certainly important. Hence, if the possibility that the property
of the matter fields or the equation of state could freely be chosen as above
is allowed for, then the analysis is certainly valid and gives several useful
conclusions on possible collapse endstates. In such a case, it is also possible
that the chosen equation of state will be in general such that the pressures
may explicitly depend not only on the energy density, but also on the time
coordinate.
From such a perspective, it is now shown and pointed out below that the
analysis, as given above, does in fact include several well-known equations
of state and useful classes of collapse models as special cases. It is also seen
that the energy conditions are satisfied throughout the collapse.
Discussed above was the idealized class of dust collapse models where the
pressures were taken to be vanishing. This has been studied extensively and
has yielded many important insights into collapse evolutions. In this special
case, the Einstein equations can be solved completely and some remarks
are made below on the N -dimensional generalization of the usual TBL dust
collapse models. The metric in this case is given as
R2
ds = dt −
2 2
dr2 − R2 (t, r)dΩ2N −2 . (3.206)
1 + r2 b0 (r)
(N − 2)F
=ρ (3.207)
2R(N −2) R
3.7 Equation of state 131
and
F (r)
Ṙ2 = + f (r). (3.208)
R(N −3)
In the case of dust, the mass function must be F = F (r), and hence the
regularity condition implies that
The energy condition here gives, in the range 0 < r < rb , M(r) ≥ 0 and
3M + rMr ≥ 0. In this case now, the function X (v ) that was discussed
earlier is given as
1 1 v (N −3)/2 (M1 + v (N −3) b1 )dv
X (v ) = − , (3.210)
2 v (M0 + v (N −3) b00 )3/2
and the time taken for the central shell to reach the singularity is
1
v (N −3)/2 dv
ts0 = . (3.211)
0 M0 + v (N −3) b00
It is now seen clearly that any given sets of density and velocity profiles at
the initial epoch completely determine the tangent to the singularity curve
at the central singularity. Also, (3.80) becomes
(N −1)/2 N − 1
x0 = M0 X (0). (3.212)
2
It therefore follows that, given any specific density profile of the collapsing
dust cloud, a velocity profile can be chosen so that the endstate of the col-
lapse would be either a naked singularity or a blackhole, depending on the
choice made, such that the energy conditions are also satisfied throughout
the collapse. The converse also holds, namely a given velocity profile can be
chosen for the cloud at the initial epoch, and then there are density pro-
files that will lead the collapse to either the naked singularity or blackhole
final states. As seen, these conclusions hold irrespective of the number of
dimensions of the spacetime. In this sense, this treatment unifies and gen-
eralizes the earlier results of the dust collapse. The interesting point that
comes out here is, given an initial density profile for the collapsing cloud, the
space of velocity profile functions is divided into the regions that lead the
collapse either to a blackhole or a naked singularity evolution, depending on
the choice made, and the converse holds similarly.
While the dust equation of state discussed above is fairly standard and
extensively used, it is widely believed that pressures would play an impor-
tant role in gravitational collapse considerations (Jhingan and Magli, 2000;
132 Spherical collapse
where k ∈ [−1, 1] is a constant. The case k = 0 gives the dust case, and
k = 1 is the stiff fluid case. Presently only the case of positive pressures is
considered. In this case k > 0, and the energy conditions give
From the above equation of state and the Einstein equations it can
immediately be seen that the function M is now the solution of the
equation
(N − 1)kM + krM,r + Q(r, v )M,v = 0, (3.220)
where
Q(r, v ) = (k + 1)rv + v. (3.221)
Now, (3.220) has a general solution of the form
F (X, Y ) = 0, (3.222)
where X (r, v, M) and Y (r, v, M) are the solutions of the system of equations
dM dr dv
− = = . (3.223)
3kM kr Q
Therefore, it can easily be seen that (3.220) admits classes of solutions when
v > 0. Also, solving the equation for the central shell r = 0, with boundary
conditions ρ → ∞ as v → 0,
m0
M(0, v ) = , (3.224)
v3k
where m0 is a constant. By choosing m0 > 0, the central shell can be made
to satisfy the energy condition ρ(t, 0) > 0 for all epochs. Then, by the
continuity of the density function, it can be said that an -ball exists around
the central shell, for
√ which v (t, r) > 0 and also ρ(t, r) > 0. But at the
central singularity, vv ≈ X (0). This implies that classes of solutions exist
that satisfy the energy conditions and also admit a naked singularity as the
collapse endstate. For a further discussion on perfect fluid collapse, and the
details of blackhole and naked singularity formation, see Goswami and Joshi
(2006) and references therein.
It is seen from the above discussion that several well-known classes of
collapse models and equation of states form subcases of the considerations
134 Spherical collapse
given here. Along with these models, the above analysis would work for any
other models with other equations of state, if solutions to the constraint
equations on M and A are permitted. Hence, it follows that the above
provides an interesting framework for the study of dynamical collapse, which
is one of the most important problems in gravitation physics today.
4
Cosmic censorship
135
136 Cosmic censorship
The spacetime model worked with must satisfy several physical reasonabil-
ity conditions in order to represent the physical universe as observed. One
of these is that the spacetime should have a regular causal structure. The
physical condition that no material particle signals can travel faster than the
velocity of light fixes the causal structure for the Minkowski spacetime. Also,
in general relativity, which uses the framework of a general spacetime mani-
fold, it was noted that locally the causality relations were the same as in the
Minkowski spacetime, which is the background manifold for special relativ-
ity. However, globally, there could be important differences in the causality
structure and other physical properties, due to a different spacetime topol-
ogy, as compared with that of the Euclidean space R4 of the Minkowski
spacetime (Geroch, 1967).
The strong gravitational fields that arise either in cosmology or gravita-
tional collapse, as signaled by the occurrence of spacetime singularities in
general relativity, also may have significant causality and topology implica-
tions for the structure of the spacetime in the vicinity of these ultra-strong
gravity regions. The causal structure of the spacetime has been studied in
detail, especially from the perspective of the occurrence of spacetime singu-
larities in collapse and cosmology (see for example, Geroch, 1970b; Penrose,
1972; Hawking and Ellis, 1973). Here, in this section, this causal structure
and its relationship with the spacetime topology are discussed. The known
results and definitions required for later chapters are also reviewed and some
results that are of intrinsic interest for the global aspects of spacetimes, and
which emphasize the close interplay between the causal structure and topo-
logy, are given. Basic causal structure ideas and definitions are given, and
topological properties of several spacetime sets and various causality con-
ditions are stated in order to arrive at a unified causality statement for a
reasonable model of spacetime. The results in the present section are valid
for manifolds of arbitrary dimensions greater than or equal to two.
An event p chronologically precedes event q , denoted by p q , if there
is a smooth, future directed timelike curve from p to q . If such a curve is
non-spacelike, namely timelike or null, p causally precedes q , or p < q . The
chronological future I + ( p) of p is the set of all points q such that p q . The
chronological past of p is defined similarly. Therefore,
I + ( p) = {q ∈ M | p q}, (4.1)
I − ( p) = {q ∈ M | q p}. (4.2)
+
I (S)
p
q
Fig. 4.1 The future of a set S is the union of the futures of all the events in
this set.
The relations and < are transitive, and for events p, q , and r, p q and
q < r or p < q and q r implies p < r (Penrose, 1972). It is seen from this
that
I + ( p) = J + ( p), (4.5)
and also
I˙+ ( p) = J˙+ ( p), (4.6)
and the chronological (causal) pasts of the subsets of the spacetimes are
defined similarly. Such dual definitions or results will often be taken as
granted.
Suppose there is a future directed timelike curve from p to q . There is a
local region containing q that the timelike curve must enter in which special
relativity is valid. Therefore, there is a neighborhood N of q such that any
point of N can be reached by a future directed timelike curve from p. It is
thus seen that for any event p ∈ M , the sets I + ( p) and I − ( p) are open in M .
The above also implies that the sets I ± (S ) are open, as they are the union
of open sets in M . However, the sets J ± ( p) are neither open nor closed in
general (see for example, Hawking and Ellis, 1973, or Joshi, 1993). Therefore,
4.1 Causal structure 139
The boundary of these points define the null cone in the tangent space
Tp . Suppose now that W and Z are any two non-null vectors in Tp ,
then,
g (W , Z ) = 1
2
[ g (W + Z, W + Z ) − g (W , W ) − g (Z, Z )] . (4.10)
4.1 Causal structure 141
has two distinct roots λ1 and λ2 , as g (X, X ) < 0 and g (Y , Y ) > 0. The
knowledge of the null cone then implies that λ1 and λ2 can be determined
in principle. But λ1 λ2 = g (X, X )/g (Y , Y ). Hence, the null cone gives the
ratio of the magnitudes of a timelike and a spacelike vector. Therefore, each
term in the equation above is determined up to a factor, and so g (W , Z ) is
determined up to a factor.
The local causality principle implies that, over small regions of space and
time, the causal structure is the same as in the special theory of relativity.
However, on a larger scale global pathological features, such as the violation
of time orientation, possible non-Hausdorff nature or non-paracompactness,
disconnected components of spacetime, and such others, may show up.
Such pathologies are to be ruled out by means of reasonable topological
assumptions, and one would like to ensure that the spacetime is causally
well-behaved. This is carried out by means of introducing various causal-
ity conditions such as the non-occurrence of closed timelike or non-spacelike
curves (causality), and the stability of this condition under small perturba-
tions in the metric (stable causality). In fact, Carter (1971) pointed out that
there is an infinite hierarchy of such causality conditions for a spacetime.
It would appear reasonable to demand that physically realistic spacetimes
do not allow either closed timelike or closed non-spacelike curves, as this
would give rise to the phenomenon of entering one’s own past. However,
general relativity and Einstein’s equations as such do not rule out such a
possibility on their own. For example, the Gödel universe (Gödel, 1949) has
closed timelike curves through each point of the spacetime. Again, the global
topology of M can cause closed timelike curves. For example, for the cylinder
M = S 1 × R, obtained from the Minkowski spacetime by identifying t = 0
and t = 1 hypersurfaces with the metric given by ds2 = −dt2 + dx2 , the
circles x = const., are closed timelike curves. In fact, for all p ∈ M here,
I + ( p) = I − ( p) = M . Such examples could be discarded as mathematical
pathologies in the spacetime topology, and the Gödel universe may be termed
unrealistic because it is a rotating model that does not correspond to the
observed universe.
More difficult to rule out are the Kerr solutions of a spinning gravitational
source (Kerr, 1963), which contain closed timelike curves if the rotation is
sufficiently fast in comparison with the value of the mass parameter. These
could possibly represent the final fate of a massive collapsing star that is
rotating. If a star failed to get rid of enough spin during the process of col-
lapse, it would give rise to a time machine in the spacetime. Wormholes in
142 Cosmic censorship
q
+
t=1
I ( p)
∩
–
I (q)
p t=0
Σ t = –1
Fig. 4.2 A globally hyperbolic spacetime has a fixed topology M = Σ×R, and is
foliated entirely by spacelike surfaces Σ. For any events p and q, the intersection
I + (p) ∩ I − (q) is compact. Such a spacetime obeys strong cosmic censorship and
allows no naked singularities.
t2 − x2 − y 2 − z 2 = const. (4.13)
are not, as the past or future null cones of the origin are boundaries of the
domain of dependence for these spacelike surfaces. Again, if a point from the
Minkowski spacetime is deleted, the resulting M admits no Cauchy surface
and M is not globally hyperbolic.
It was shown by Geroch (1970a) that a globally hyperbolic M has a unique
topological structure, and it admits no topology change in that M is then
homeomorphic to R × S , where S is a three-dimensional submanifold and
for each t ∈ R, {t} × S is a Cauchy surface for M .
The basic idea of the proof for the above involves introducing a finite
measure µ on M so that µ(M ) = 1. Then, a function h : M → R is
introduced by
µ(J + ( p))
h+ ( p) = . (4.14)
µ(J − ( p))
The function h− is defined similarly. The sets h± = const. are seen to
be Cauchy surfaces for M . Such functions are called causal functions on
a spacetime.
It can be asked if the converse is true in some sense, that is, whether the
direct product spacetimes (M, g ) with M = S × T , where each S × {t} is
spacelike and each {x} × T is timelike, are always globally hyperbolic. The
answer is no; refer to Clarke and Joshi (1988), who discuss a spacetime that
is a direct product as above, but not globally hyperbolic.
Next, let S be a partial Cauchy surface. Then N = D+ (S ) ∪ D− (S ) = M
and N must be a proper subset of M . The boundary of N in M can be
146 Cosmic censorship
divided into two regions H + (S ) and H − (S ) that are respectively called the
future and past Cauchy horizons of S ,
The hope was then that, when more general solutions were considered
with a lesser degree of symmetry requirements, such singularities would be
avoided.
This issue was sorted out when a detailed study of the structure of a
general spacetime and the associated problem of a spacetime singularity was
taken up by Hawking, Penrose, and Geroch (see for example, Penrose, 1968;
Geroch, 1971; Hawking and Ellis, 1973). It was shown by this work that a
spacetime will admit singularities within a very general framework, provided
that it satisfies certain reasonable assumptions, such as the positivity of
energy, a suitable causality assumption, and a condition such as the existence
of trapped surfaces. It thus follows that the spacetime singularities form a
general feature of the relativity theory. In fact, these considerations also
ensure the existence of singularities in other theories of gravity that are
based on a spacetime manifold framework and satisfy the general conditions
stated above.
with the range of coordinates given by −∞ < t < ∞, 0 < r < ∞, 0 < θ < π ,
but with 0 < φ < a, with φ = 0 and φ = a identified and a = 2π . There is a
conical singularity at r = 0 through which the spacetime cannot be extended
and the singular boundary is related to the timelike two-plane r = 0 of the
Minkowski spacetime.
The important question that arises is whether such singularities develop
even when a spacetime of generality is considered and if so, under what con-
ditions. In order to consider this question, it is first necessary to characterize
more precisely what is meant by a spacetime singularity.
While trying to characterize a spacetime singularity, the first point to note
is that by very definition, the metric tensor has to be well-defined at all the
regular points of the spacetime. Since this is no longer true at a spacetime
singularity such as those discussed above, a singularity cannot be regarded
as a regular point of the spacetime, but must be treated as a boundary
point attached to it. This situation causes difficulty when one attempts to
characterize a singularity by the criterion that the curvatures must blow up
near the singularity. The trouble is, since the singularity is not a part of the
spacetime, it is not possible to define its neighborhood in the usual sense in
order to discuss the behavior of curvature quantities in that region.
Characterizing the singularity could be attempted in terms of the diver-
gence of the components of the Riemann curvature tensor along non-spacelike
trajectories of the spacetime. The trouble with this is that the behavior of
such components will, in general, change with the change of frames used, and
this approach is not really of much help. The curvature scalars or the scalar
polynomials in the metric and the Riemann tensor could be used, and they
could be required to achieve unboundedly large values. This is the behavior
encountered in the Schwarzschild and the Friedmann models. However, it is
possible that such a divergence of curvature scalars occurs only at infinity for
a given non-spacelike curve. In general, it looks reasonable to demand that
some sort of curvature divergence must take place along the non-spacelike
curves that encounter a spacetime singularity. However, a general charac-
terization of the singularity in terms of the curvature divergence runs into
various difficulties. For example, for the plane wave vacuum solutions, the
polynomials in the curvature scalars vanish, but the curvature tensor is still
4.2 Spacetime singularities 153
of such surfaces cannot contain any conjugate points, and these must leave
the hypersurface before encountering a conjugate point. This puts strong
constraints on the nature of such surfaces, and the singularity theorems result
from an analysis of these limits.
Consider a congruence of timelike geodesics in the spacetime. This is a
family of curves such that precisely one timelike geodesic trajectory passes
through each point p. Choosing the curves to be smooth, this defines a
smooth timelike vector field on the spacetime. On the other hand, a given
smooth vector field on the spacetime specifies a congruence of curves in the
manifold.
Let V i denote the timelike tangent vector to the congruence. Choosing the
parameter to be the proper time along such timelike trajectories, this can be
normalized to be a unit tangent vector,
V i Vi = −1. (4.17)
and also
hij hij = hi i = δ i i + V i Vi = 3. (4.21)
For the given congruence of timelike geodesics, the expansion, shear, and
rotation tensors are respectively defined as
Furthermore, note that σij and ωij are purely spatial quantities in the sense
that
σij V i = ωij V i = 0. (4.26)
Also, note that
σ i i = hij σij = θ − 13 hij hij = 0. (4.27)
The covariant derivative of V is then expressed as
V k ∇k ∇j Vi = V k ∇j ∇k Vi + Rilkj V l V k . (4.29)
dθ
= V k ∇k V i ;i = −(V k ;i V i ;k ) − Rlk V l V k , (4.31)
dτ
where τ is the affine parameter along the geodesic.
Using (4.28) in the above and the anti-symmetry properties of the tensor
ωij , after some simplification,
dθ
= −Rlk V l V k − 13 θ2 − σij σ ij + ωij ω ij (4.32)
dτ
can be obtained, which can be written as
dθ
= −Rlk V l V k − 13 θ2 − 2σ 2 + 2ω 2 . (4.33)
dτ
Equation (4.33) above is called the Raychaudhuri equation (Raychaudhuri,
1955), which describes the rate of change of the volume expansion as the
timelike geodesic curves in the congruence are moved along.
The second and third terms on the right-hand side involving θ and σ are
always positive. Consider now the term Rij V i V j , which by the Einstein
equations can be written as
The term Tij V i V j above represents the energy density measured by a time-
like observer with the unit tangent V i , which is the four-velocity of the
observer. For all reasonable classical physical fields this energy density is
generally taken as non-negative, and it is assumed that for all timelike
vectors V i ,
Tij V i V j ≥ 0, (4.35)
is satisfied. Such an assumption is called the weak energy condition. On the
other hand, it is also considered reasonable to believe that the matter stresses
will not be so large as to make the right-hand side of (4.34) negative. This
will be satisfied when
Tij V i V j ≥ − 12 T (4.36)
is satisfied. Such an assumption is called the strong energy condition, and it
implies that for all timelike vectors V i ,
Rij V i V j ≥ 0. (4.37)
By continuity, it can be argued that the same will then also hold for all null
vectors.
Both the strong and weak energy conditions will be valid for well-known
forms of matter, such as the perfect fluid, provided that the energy density ρ
is non-negative and that there are no large negative pressures that are bigger
or comparable to ρ, when converted into physical units.
An additional energy condition often required by the singularity theo-
rems is the dominant energy condition, which states that in addition to the
weak energy condition, the pressure of the medium must not exceed the
energy density. This can be equivalently stated as, for all timelike vectors
V i , Tij V i V j ≥ 0 and the vector T ij Vi is a non-spacelike vector. Such a
condition would be satisfied provided that the local speed of sound does not
exceed the local speed of light.
With the strong energy condition being satisfied, the Raychaudhuri equa-
tion implies that the effect of matter on the spacetime curvature causes a
focusing effect in the congruence of the timelike geodesics due to gravita-
tional attraction. This, in general, causes the neighboring geodesics in the
congruence to cross each other, which gives rise to caustics or conjugate
points. This separation between nearby timelike geodesics is governed by the
geodesic deviation equation,
D2 Z j = −Rkil
j
V k Z iV l, (4.38)
Suppose now γ is a timelike geodesic. Then, two points p and q along γ are
called conjugate points if a Jacobi field along γ exists that is not identically
zero, but vanishes at p and q . From the derivation of the Raychaudhuri
equation given above, it is clear that the occurrence of conjugate points
along a timelike geodesic is closely related to the behavior of the expansion
parameter θ of the congruence. In fact, it can be shown that the necessary
and sufficient condition for a point q to be conjugate to p is that for the
congruence of timelike geodesics emerging from p, θ → −∞ at q (see for
example, Hawking and Ellis, 1973). The conjugate points along the null
geodesics are also similarly defined. Consider, for example, a congruence
of null geodesics emanating from a point p. If infinitesimally nearby null
geodesics of the congruence meet again at some other point q in the future,
then p and q are said to be conjugate to each other.
Similarly, let S be a smooth spacelike hypersurface, that is, it is an
embedded three-dimensional submanifold. Consider a congruence of time-
like geodesics orthogonal to S . Then, a point p along a timelike geodesic γ
of the congruence is said to be conjugate to S along γ if a Jacobi vector field
along γ exists that is non-zero at S , but vanishes at p. This means that there
are two infinitesimally nearby geodesics orthogonal to S that intersect at p.
Again, the equivalent condition for this to happen, in terms of the parameter
θ, is that the expansion θ for the congruence orthogonal to S tends to −∞
at p. If V i denotes the unit timelike tangent vector field of the congruence
of timelike geodesics, where V i denotes the normal to S , then the extrinsic
curvature χij of S is defined as
χij = ∇i Vj , (4.39)
dθ θ2
≤− , (4.41)
dτ 3
discussed. Using the causal structure analysis, it is shown that there must be
maximal length timelike curves between certain pairs of events in the space-
time. As pointed out above, a causal geodesic that is both future and past
complete must contain pairs of conjugate points if the spacetime satisfies
the generic condition and an energy condition. This is then used to draw
the necessary contradiction in order to show that M must be non-spacelike
geodesically incomplete. This result can be stated as in the following.
A spacetime (M, g ) cannot be timelike and null geodesically complete if
the following are satisfied:
(1) Rij K i K j ≥ 0 for all non-spacelike vectors K i ;
(2) the generic condition is satisfied, that is, every non-spacelike geodesic
contains a point at which K[i Rj]el[m Kn] K e K l = 0, where K is the
tangent to the non-spacelike geodesic;
(3) the chronology condition holds;
(4) in M , either a compact achronal set without edge or a closed trapped
surface exists, or a point p exists such that for all past directed null
geodesics from p, θ must eventually be negative.
The main idea of the proof is that it can be shown that the following three
cannot hold simultaneously: (a) every inextendible non-spacelike geodesic
contains pairs of conjugate points, (b) the chronology condition holds, (c) an
achronal set S exists in the spacetime such that E + (S ) or E − (S ) is compact,
where E + (S ) and E − (S ) denote the future and past of the edge of set S .
If this is shown, then the theorem is proved, as (3) is the same as (b),
(4) implies (c), and (1) and (2) imply (a). First, note that (a) and (b)
imply the strong causality of the spacetime (see Prop. 6.4.6 of Hawking and
Ellis, 1973). Next, it can be shown that if S is a future trapped set and if
strong causality holds on I + (S ) then a future endless trip γ exists such that
γ ⊂ Int D+ (E + (γ )). Now, T = J − (γ ) ∩ E + (S ) is defined, T turns out to be
past trapped, and hence λ, a past endless causal geodesic in Int(D− (E − (T )),
exists. Then, a sequence {ai } receding into the past on λ, and a sequence
{ci } on γ to the future is chosen. The sets J − (ci ) ∩ J + (ai ) are compact and
globally hyperbolic, so a maximal geodesic µi from ai to ci exists for each i.
The intersections of µi with the compact set T have a limit point p and a
limiting causal direction. The causal geodesic µ with this direction at p must
have a pair of conjugate points. This is then shown to be contradictory to
the maximality property of the geodesics stated above.
There is a more general way in which the singular points in a spacetime
can be defined using the terminal indecomposable pasts (TIPs) and terminal
indecomposable futures (TIFs), as discussed earlier (Penrose, 1974a, 1979).
The spacetime is assumed to be strongly causal. Here, a curve means a map
γ from an interval [0, a) of the real line into M , where a could possibly be
infinity. Therefore, the curve starts at an initial point γ (0) with a definite
4.3 Blackholes 161
p q
I –(p) I –(g)
I –(q)
Fig. 4.3 The event p is a regular point in M and the set I − (p) defines a PIP. For
a non-spacelike curve γ going to infinity, I − (γ) gives an ∞-TIP. If λ is a finite
length curve going into a spacetime singularity, then I − (λ) is a singular TIP.
tangent but has no end point, as the interval is open at a. Such a curve will
be called extendible if it is possible to extend the map γ to an end point γ (a)
in M , and otherwise it is called inextendible. Of particular interest here are
the inextendible non-spacelike curves. The TIPs and TIFs are generated by
future directed and past directed timelike curves respectively, and they give
all the boundary points of spacetime that include both the singularities and
points at infinity. Such a boundary point is called ∞-TIP, which is a point
at infinity if it is generated by some timelike curve of infinite proper time
length in the future. A singular TIP is one that is not generated by any such
timelike curve of infinite length (see Fig. 4.3). Similarly, ∞-TIFs and singular
TIFs can be defined. The existence of a singular TIP defines a singularity of
spacetime, giving a class of future directed inextendible timelike curves that
have a finite proper time length, but no future end point.
As pointed out by Clarke (1986), the basic requirement for the ideal end
point of a timelike curve to be called a singularity, rather then a regular
boundary point, is that there should be no extension of the spacetime possible
in which the curve in question could be continued. If such an extension
existed, then the singularity would be similar in some sense to the coordinate
singularity in the Schwarzschild geometry at r = 2m. Therefore, the question
of singularity depends on what type of extension is allowed for the spacetime.
Therefore, a boundary point is called a C k -singularity of the spacetime if
there is no C k -extension of M that removes it. Clarke then defines the index
k as a measure of the strength of the singularity in the sense that the smaller
the k , the stronger the singularity.
4.3 Blackholes
Assuming cosmic censorship in the form of an asymptotic predictability in
the spacetime, much of the theoretical basic blackhole physics is carried out
in an asymptotically flat and predictable model.
162 Cosmic censorship
large near I . The null geodesics of these two conformally related spacetimes
are completely identical as point sets, as mentioned earlier. However, the
affine parameters along the geodesics in M and M are related as dv̄ =
Ω2 dv . Therefore, the affine parameter along the null geodesics in M must
blow up near I + in the future, and similar behavior holds near the past
infinity.
The definition above is quite stringent in that it assumes that every null
geodesic has two end points: in the future and in the past at I . Although this
is satisfied in the Minkowski spacetime, this does not hold in spacetimes such
as the Schwarzschild and Reissner–Nordström cases, which contain event
horizons and a blackhole region. Here, the future directed null geodesics
that enter the blackhole must end in the singularity at r = 0 and cannot
have an end point at I + . One would like to include these spacetimes in
the general asymptotically flat class, and so a spacetime M is defined to
be weakly asymptotically simple and empty (or a WASE spacetime) if an
asymptotically flat spacetime M in the above sense exists such that there
is a neighborhood of I in M that is isometric to an open set in M . This
definition covers the Schwarzschild and the Reissner–Nordström cases and
also the Kerr solutions.
Some important general properties of blackholes are now discussed. For
a detailed treatment on blackhole physics and results such as the blackhole
uniqueness theorems and details, see Hawking and Ellis (1973) and Wald
(1984). As stated earlier, the fundamental motivation for the concept of
a blackhole comes from a spherically symmetric homogeneous dust collapse
that has two important features. First, for a star undergoing a complete grav-
itational collapse, a region of trapped surfaces forms below r = 2m, from
which no light rays escape to an observer at infinity. Therefore, a blackhole
forms in the spacetime. Second, the ultimate fate of the star undergoing
the collapse is an infinite curvature singularity at r = 0, which is com-
pletely hidden within the trapped surface region and the blackhole. Hence,
no emissions or light rays from the singularity could go out to any observer
at infinity, and the singularity is causally disconnected from the outside
spacetime.
The question now is whether these conclusions can be generalized for a non-
spherically symmetric collapse, and whether they are valid at least for small
perturbations from exact spherical symmetry. It is known from Hawking
and Ellis (1973), using the stability of a Cauchy development in general
relativity, that the formation of trapped surfaces, and hence of a blackhole,
is a stable property when departures from spherical symmetry are taken into
account. Considering a spherically symmetric collapse evolution from given
initial data on a partial Cauchy surface S , the formation of trapped surfaces
T in the form of all the spheres with r < 2m in the exterior Schwarzschild
geometry is found. The stability of the Cauchy development then implies
164 Cosmic censorship
that, for all initial data sufficiently near to the original data in the compact
region J + (S ) ∩ J − (T ), the trapped surfaces must still occur. Then, the
curvature singularity of a spherical collapse also turns out to be a stable
feature as implied by the singularity theorems, which show that the closed
trapped surfaces always imply the existence of a spacetime singularity under
reasonable general conditions.
There is no proof available, however, that the singularity will continue
to be hidden within the blackhole and remain causally disconnected from
outside observers, even when the collapse departs from the homogeneous
dust cases, or is not exactly spherical. If the singularity became visible to
external observers, the predictability in the spacetime will be undermined
because new information could come from the singularity where the densities
and curvatures could be arbitrarily large.
Hence, in order to generalize the notion of blackholes to gravitational
collapse situations other than exactly spherically symmetric homogeneous
dust cases, it becomes necessary to rule out such naked singularities by means
of an explicit cosmic censorship assumption. This could be stated as follows:
if S is a partial Cauchy surface, then there are no naked singularities to the
future of S that can be seen from the future null infinity I + . This is true for
the spherical homogeneous dust collapse, where the breakdown of physical
theory at the spacetime singularity does not disturb the prediction in the
future for the outside asymptotically flat region.
This assumption is made precise by considering the spacetimes (M, g )
that admit a weakly asymptotically simple and empty conformal completion
(M , ḡ ). Then, (M, g ) is said to be future asymptotically predictable from a
partial Cauchy surface S if
B = M − J − (I + ). (4.43)
4.3 Blackholes 165
I+ I+
p
r
q
Fig. 4.4 No non-spacelike curves from the blackhole region B can reach the
future infinity and the events faraway in the spacetime, so it is cut off from the
outside universe. The events p and q are in a blackhole, whereas the event r that
is outside the blackhole can send signals to infinity.
is called the event horizon. As pointed out earlier, this horizon must be an
achronal surface generated by null geodesics that could have past end points
in M , but that have no future end points. For the Minkowski spacetime,
J − (I + ) = M and there is no blackhole. However, for the Schwarzschild case,
J − (I + ) is the region for the spacetime exterior to r = 2m, and the event
horizon is given by the null hypersurface r = 2m, which is the boundary of
the blackhole region 0 < r < 2m.
By definition, no spacetime singularities are visible at null infinity in an
asymptotically predictable spacetime. In fact, this is true also for trapped
surfaces in this case, provided either the weak or strong energy condition
is satisfied, which implies that all trapped surfaces must be fully contained
within the blackhole region and not visible from I + (Wald, 1984). Specifi-
cally, let (M, gij ) be future asymptotically predictable from a partial Cauchy
surface S , and Rij K i K j ≥ 0 for all null vectors K i . Then, if T is a closed
trapped surface in D+ (S ), T ∩ J − (I + ) = ∅.
Let (M, gij ) be an asymptotically flat spacetime with the associated
unphysical conformal spacetime (M , ḡij ). Suppose p ∈ I + and q ∈ M ∩
J − ( p). Let γ be the future directed null geodesics generator of I + through
p, and let r ∈ γ be any point. Then, q ∈ I − (r) in M . Therefore,
J − (I + ) = I − (I + ) in M , and hence J − (I + ) is open in M . Therefore, the
blackhole region B = M − J − (I + ) is closed in M . This implies that the
event horizon is contained in B . As such, the blackhole region B need not
166 Cosmic censorship
2GM
v02 = . (4.47)
R
When the radial velocity of the particle is less than v0 , it must fall back to
the body, otherwise it escapes to infinity. Therefore, if the mass distribution
and radius of the body were such that c2 = 2GM/R, where c is the velocity
of light, then for any larger mass or smaller radius of the body, even light
will not escape. This was realized by Laplace in 1798, who pointed out that,
for a star with a density the same as the Sun but radius 250 times larger, no
light could escape from its surface.
It is clear that a blackhole could not be observed directly, but gravitational
effects exhibited by such an object must be looked for. Although there is no
conclusive evidence available for the existence of blackholes at the moment,
presently the best candidates seem to be the binary stars in which one of the
partners is visible and the other is supposed to be a blackhole. Such a black-
hole would suck matter from its visible component, in the process forming
an accretion disk around the blackhole. Before the in-falling matter spirals
down the blackhole, the inner, hot regions are believed to produce intense
bursts of X-rays formed by synchrotron radiation. Therefore, the discovery of
the X-ray source Cygnus XI in 1971, which shows rapid variations, indicates
possible evidence for blackholes. Further to this, several other X-ray binaries
have been proposed as possible candidates for blackholes.
4.4 Higher spacetime dimensions 169
where i−1
N −2
2
dΩ2N −2 = sin (θ ) (dθi )2
j
(4.49)
i=1 j=1
Take the matter field to satisfy the weak energy condition, that is, the
energy density measured by any local observer is non-negative, and so for
any timelike vector V i ,
In the case of a finite collapsing cloud, there is a finite boundary 0 < r < rb ,
outside which the cloud is matched to a Schwarzschild exterior. The range of
the coordinates for the metric is then 0 < r < rb , and −∞ < t < ts (r), where
ts (r) corresponds to the singular epoch R = 0. Solving the N -dimensional
Einstein equations, the generalized Tolman–Bondi–Lemaitre (TBL) metric
can be obtained as
R2
ds = −dt +
2 2
dr2 + R2 (t, r)dΩ2N −2 , (4.52)
1 + f (r)
where f (r) is an arbitrary function of the comoving radius r, and f (r) > −1.
The equations of motion are then given by
(N − 2)F F (r)
= ρ, Ṙ2 = + f (r). (4.53)
2R(N −2) R R(N −3)
Here, F (r) is an arbitrary function of the comoving coordinate r and has
the interpretation of mass function of the dust cloud, as discussed earlier,
and f (r) is the energy function, which specifies the velocity profiles for the
collapsing shells for any given mass function. The energy condition then
implies F ≥ 0. Consider the shell-focusing naked singularity at R = 0. The
continual collapse condition is given as Ṙ < 0.
Using the scaling independence of the comoving coordinate r,
where
v (ti , r) = 1, v (ts (r), r) = 0, v̇ < 0. (4.55)
This means the coordinate r has been scaled in such a way that at the initial
epoch, R = r, and at the singularity R = 0. As noted previously, R = 0 both
4.4 Higher spacetime dimensions 171
where F (r) and f (r) are the mass function and the energy function
respectively.
Following Goswami and Joshi (2006), now assume that the initial density
and energy functions ρ(r) and f (r) are smooth and even, ensuring their
analytic nature. Note that the Einstein equations do not impose any such
restriction, which has to be physically motivated, and this implies a certain
mathematical simplicity in the arguments to deal with a dynamical collapse
situation. It follows that both M(r) and b(r) are now smooth C ∞ -functions,
which means the Taylor expansions of these functions around the center must
be of the form
M(r) = M0 + r2 M2 + r4 M4 + · · · , (4.57)
and
b(r) = b0 + r2 b2 + r4 b4 + · · · , (4.58)
that is, all odd terms in r vanish in these expansions, and the presence of
even terms only would ensure smoothness.
To predict the final state of the collapse for a given initial mass and velocity
distribution, the singularity curve resulting from the collapse of successive
matter shells, and the apparent horizon developing in the spacetime are
studied below (see also Goswami and Joshi, 2004a). A decreasing apparent
horizon in the (t, r) plane then is a sufficient condition for a blackhole, as it
shows the entrapment of the neighborhood of the center before the central
singularity.
Consider the continual collapse of the dust cloud to a final shell-focusing
singularity at R = 0, where all matter shells collapse to a zero physical
radius. With the regular initial conditions as above, the Einstein equation
172 Cosmic censorship
Here, the negative sign implies that v̇ < 0, that is, the matter cloud is
collapsing. Integrating the above equation with respect to v gives
1
v (N −3)/2 dv
t(v, r) = . (4.60)
v M(r) + v (N −3) b(r)
X2 (v ) X3 (v )
t(v, r) = t(v, 0) + rX (v ) + r2 + r3 + ··· . (4.61)
2 6
Now, from the equation for the form of the mass function above and (4.58),
X (v ) = 0. The function X2 (v ) is then given as
1
v (N −3)/2 (M2 + v (N −3) b2 )dv
X2 (v ) = − 3/2
, (4.62)
v (M0 + v (N −3) b0 )
where
Note that the value of X2 is determined fully by the initial values of the mass
function F (r) and the energy function b(r). Therefore, the time taken for
the central shell to reach the singularity is given by
1
v (N −3)/2 dv
ts0 = . (4.64)
0 M0 + v (N −3) b0
From the above equation it is clear that for ts0 to be defined,
that is, the continual collapse condition implies the positivity of the above
term. Hence, the time taken for other shells close to the center to reach the
singularity can now be given by
X2 (0)
ts (r) = ts0 + r2 + ··· . (4.66)
2
In order to determine the visibility, or otherwise, of the singularity at R =
0, the causal structure of the trapped surfaces, and the nature and behavior
4.4 Higher spacetime dimensions 173
F
= 1. (4.67)
RN −3
Broadly, it can be stated that if the neighborhood of the center gets
trapped earlier than the singularity, then it is covered, otherwise it is naked
with families of non-spacelike future directed trajectories escaping away from
it. For example, it follows from the above equation that along the singularity
curve t = ts (r) (which corresponds to R = 0), for any r > 0, F (r) goes
to a constant positive value, whereas R → 0. Hence, it follows that trap-
ping already occurs before the singularity develops at any r > 0 along the
singularity curve ts (r).
Now, there is a need to determine when there will be families of non-
spacelike paths coming out of the central singularity at r = 0, t = ts (0),
reaching outside observers, and when there will be none. The visibility,
or otherwise, of the singularity is decided accordingly. By determining the
nature of the singularity curve and its relation to the initial data, it is pos-
sible to deduce whether the trapped surface formation in the collapse takes
place before or after the central singularity. It is this causal structure that
determines the possible emergence, or otherwise, of non-spacelike paths from
the singularity, and settles the final outcome in terms of either a blackhole
or a naked singularity.
From (4.67),
vah (r) = [r2 M(r)]1/(N −3). (4.68)
Using the above equation in (4.60),
vah (r)
v (N −3)/2 dv
tah (r) = ts (r) − . (4.69)
0 M(r) + v (N −3) b(r)
Event horizons
Singularities
Apparent horizons
(a) (b)
Fig. 4.5 Typical apparent horizon behavior in a blackhole (a) and a naked sin-
gularity (b) collapse. Whereas in the blackhole case the apparent horizon is
decreasing away from the center, in the naked singularity situation it typically
increases in time away from the center.
4.5 Formulating the censorship 175
1
Tij ≡ Gij . (4.71)
8π
176 Cosmic censorship
weak (Tipler, Clarke, and Ellis, 1980; Israel, 1986a, 1986b; Newman, 1986;
Newman and Joshi, 1988; Nolan, 1999). Some of these suggestions are dis-
cussed here. In order to avoid the difficulties associated with the question as
to which forms of matter and equations of state should be considered rea-
sonable, it is also suggested that one examines first a purely vacuum version
of the censorship. That will show whether pure gravity allows naked singu-
larities. In fact, Geroch and Horowitz (1979) have detailed several possible
approaches to the censorship formulations, and pointed out difficulties in
each case.
To summarize the situation, while the cosmic censorship hypothesis is a
crucial assumption underlying all blackhole physics, gravitational collapse
theory, and many important related areas in gravity physics, no proof and
formulation is available today. The first major task here is actually to formu-
late rigorously a satisfactory version of the hypothesis. The proof of cosmic
censorship would confirm the already widely accepted and applied theory
of blackholes, while its overturn would throw the blackhole dynamics into
serious doubt. Therefore, cosmic censorship turns out to be one of the most
important issues for general relativity and gravitation theory today. Even if
true, a proof for this conjecture does not seem possible unless some major
theoretical advances of the mathematical techniques and the understanding
of the global structure of the Einstein equations are made.
This situation leads to the conclusion that the first and foremost task here
is to carry out a detailed and careful examination of gravitational collapse
scenarios that possibly give rise to a naked singularity formation. Until this is
done properly, trying out different formulations for censorship may not help,
because without really knowing what is involved in the collapse processes,
one may be in for a complete surprise as far as the final fate of a gravitational
collapse is concerned. It is only such investigations of general collapse situ-
ations that could indicate which theoretical advances to expect for a proof,
and which features to avoid while formulating the cosmic censorship.
Under this situation, while no mathematical formulation or proof of cen-
sorship is currently available, another alternative is to ask if it would be
possible to impose suitable physical constraints on the gravitational collapse
so as to ensure the validity of the censorship. In other words, physically real-
istic collapses should satisfy the censorship. Then again, the precise physical
conditions under which the censorship is supposed to be holding are to be
specified. The advantage then is that even if a certain set of physical condi-
tions did not work towards proving the censorship, there is still the option
of trying out another set of physical constraints to continue further efforts.
Eventually, this may lead to an appropriate mathematical formulation of
the censorship conjecture to be established. Many natural looking physical
conditions can be proposed and tried out, with these being indicated as a
remedy to rule out naked singularities.
180 Cosmic censorship
This is one of the basic conditions assumed in the classical gravity descrip-
tion, and it should be satisfied by the matter fields constituting the star, at
least until the collapse has proceeded to such an advanced stage so as to enter
a phase governed by quantum gravity. This is the stage when the classical
description starts breaking down in one way or another.
As noted earlier, if completely arbitrary matter fields are allowed for, it
is quite easy to produce naked singularities. For example, one possibility is
to start with a geometry allowing families of future directed non-spacelike
geodesics, which are future endless, but which terminate in the past at the
singularity. Then, define the matter fields to be given by Tij ≡ (1/8π )Gij .
It is thus obvious that in gravitational collapse consideration must be given
to scenarios where matter fields do satisfy reasonable physical conditions. It
would be hoped that a suitable energy condition would be one of these, as all
observed classical fields do obey such a condition. A further motivation would
be the energy conditions that have been used extensively in the singularity
theorems in general relativity, and which predict the existence of singularities
in gravitational collapse and cosmology.
It would be nice to see if the censorship is obeyed once the matter fields
have been assumed to satisfy suitable energy conditions. It turns out, how-
ever, that there are several classes of collapse models where collapsing matter
does satisfy a proper energy condition, but the collapse leads to an endstate
that is a naked singularity.
Actually, there are classes of collapse models where satisfying the energy
condition appears to be aiding the naked singularity formation as the final
state of collapse, in turn making the naked singularity physically more
4.5 Formulating the censorship 181
where the metric functions are taken to depend on X = t/r due to self-
similarity. The outgoing null geodesics (Joshi and Dwivedi, 1992, 1993a)
can be worked out from the naked singularity, which turn out to be related
to the density and pressure distributions in the spacetime via the Einstein
equations. These are then given by
This is one of the most important physical constraints necessary for any
possible version of a cosmic censorship statement. In general, the regularity
conditions on the initial data for the collapse can come in many forms. If
realistic collapse scenarios of matter clouds such as gravitationally collapsing
massive stars are to be modeled, then the densities, pressures, and other
physical quantities must be finite and regular at the initial spacelike surface
from which the collapse develops. That is, the initial surface should not admit
any density or curvature singularities in the initial data so as to represent a
collapse from regular matter distribution.
Generally, this is ensured by imposing the usual differentiability conditions
on the functions involved, together with requirements of finiteness and regu-
larity. It is known from the gravitational collapse analysis so far, that regular
182 Cosmic censorship
This has been one of the most useful physical requirements, which was
explored rather thoroughly in order to develop a formulation for the cos-
mic censorship conjecture. The idea has been that any singularity that will
develop from a realistic collapse has got to be physically serious in various
aspects, including powerful divergences in all important physical quantities
such as densities, pressures, curvatures and others, at least at the classical
level. A typical condition for the singularity to be gravitationally strong is,
in addition to the divergences such as above, the gravitational tidal forces
must diverge and all physical volumes are crushed to zero size in the limit of
approach to the naked singularity. A sufficient condition for this to happen is
1
Rij V i V j ∝ , (4.75)
k2
4.5 Formulating the censorship 183
physically important equations of state are also included. The collapse can
be studied for a general form of matter, the so-called type I matter fields (all
the known physical forms of matter, such as dust, perfect fluids, and massless
scalar fields are included in this class), subject to the weak energy condition.
The result is, given an arbitrary but regular distribution of matter on the ini-
tial surface, that there are always evolutions available from this initial data
that would either result in a blackhole or naked singularity, depending on the
allowed choice of free functions available from the Einstein equations. More
specifically, in a spherically symmetric collapse with a type I general matter
field, given the distribution of density and the radial and tangential pressure
profiles on the initial surface from which the collapse develops, it is then
possible to choose the free function describing the velocities of the in-falling
shells in such a manner so as to have a blackhole or a naked singularity as
the final end product, depending on this choice.
An interesting feature is that, while there are several widely used and
familiar equations of state available that result in the formation of naked
singularities as the final fate of the collapse, there is still not a single equation
of state available so far that ensures necessarily that the end product will be
a blackhole only. Under this situation, it is quite possible that the physics
that causes collapse endstates may not be directly related to the equation of
state or the form of matter collapsing.
A more promising alternative could be to work with a rather different
representation of the matter, such as the one given by the Einstein–Vlasov
statistical description (see for example, Rendall, 2005).
the very late stages of the collapse, once the densities and curvatures are
sufficiently high so that quantum gravity effects become important in the
process of a continual collapse. Such quantum effects would come into play
much before the actual formation of the classical naked singularity, which
itself may possibly be smeared out by the quantum effects.
The key question then is that of the possible visibility, or otherwise, of
these extreme strong gravity regions, which do develop in the vicinity of
the classical naked singularity. It is then the causal structure, that is, the
communicability, or otherwise, of these extreme strong gravity regions where
quantum gravity should prevail, which would make the essential difference as
far as the physical consequences of a naked singularity are concerned, rather
than various purely classical aspects such as red-shift.
As seen above, the physical conditions are not able effectively to rule out
naked singularities, which in turn may lead to some possible formulation of
the cosmic censorship conjecture, either a physical or a mathematical one.
With each of the above conditions, there are counter-examples that obey
such a physical constraint, but which produce a naked singularity as the
endstate of a dynamical collapse.
There are three further possibilities that are under active investigation
today towards a possible formulation of the censorship conjecture, and which
may offer a better hope for a cosmic censorship hypothesis. These are now
briefly discussed below.
Outgoing
light rays
Quantum gravity
fuzzy region
Collapsing Collapsing
star star
Fig. 4.6 Quantum effects may resolve the spacetime singularity, either covered
or naked. If the collapse outcome is a naked singularity, there may be the oppor-
tunity to observe the quantum gravity effects taking place in the ultra-strong
gravity regions in the universe.
gravity whenever such an event occurs in the universe (see for example, Vaz
and Witten, 1994, 1995, and Wald, 1997). Then, the gravitational collapse
phenomena could provide a possibility of actually testing the laws of quantum
gravity, and every time a massive star collapses in the universe, there is
potentially a laboratory to test the laws of quantum gravity (see Fig. 4.6).
For an earlier discussion on quantum effects near singularities, see Kodama
(1979) and Hiscock, Williams, and Eardley (1982).
In the case of a blackhole developing in the collapse of a finite sized object
such as a massive star, such strong gravity regions have got to be necessarily
hidden behind an event horizon of gravity, which would be well before the
physical conditions became extreme. Then the quantum effects, even if they
caused qualitative changes closer to the singularity, will be of no physical
consequence. This is because no causal communications are then allowed
from within such horizons. On the other hand, if the causal structure were
that of a naked singularity, then communications from such a quantum grav-
ity dominated extreme curvature ball would be visible in principle, either
directly or via secondary effects such as shocks produced in the surrounding
medium. Some of these issues will be discussed further in the next chapter.
There has been a suggestion that all naked singularities, whenever they are
produced by matter fields such as dust, perfect fluids, and such others, should
188 Cosmic censorship
It would appear that this is the key issue on which any possible future
formulation and proof of the censorship has to crucially depend. Even if
naked singularities do develop in collapse models, if they were not generic
and stable in some suitably well-defined sense, that would make a good case
for censorship. For example, most of the current classes of naked singu-
larities discussed here are within the framework of a spherically symmetric
collapse. While there are some indications that naked singularities do develop
in non-spherical collapses as well, as discussed in the next chapter, such a
non-spherical collapse remains largely uncharted territory, and it would be
essential to examine it rather thoroughly.
The key question to be resolved here is, whilst it is known that physically
reasonable initial data do give rise to naked singularities, will the initial data
subspace, which gives rise to the naked singularity as the collapse endstate,
4.5 Formulating the censorship 189
taken seriously. The key difficulty to addressing this question is that there
are no well-defined criteria or formalisms available in the general theory of
relativity to test the stability and genericity. As opposed to Newtonian the-
ory, there is no well-defined notion for stability available in Einstein’s theory.
Under this situation, there are a variety of ways in which the above question
can be asked and treated, and there is no unique answer available in a general
manner. Many times, it is the physics of the situation that guides the path
adopted, as pointed out above. Nevertheless, this is an important and basic
issue to try to answer.
The singularity theorems establish the existence of spacetime singularities
in the form of incomplete non-spacelike geodesics, both for gravitational
collapse and cosmology. These theorems, however, give no information on the
nature of the singularities, such as whether they occur in the past or future,
the possible growth of curvature in the limit of approach to the singularity,
and whether they will be covered by an event horizon hidden from all outside
observers, or whether they may be causally connected to external observers
in the universe. Hence, while dealing with gravitational collapse scenarios it
becomes inevitable to assume the cosmic censorship in the form of a future
asymptotic predictability of the spacetime, in order to make any progress in
the theory and applications of blackholes. It has to be ensured, through the
censorship conjecture, that the singularities of the collapse are necessarily
hidden inside an event horizon. This assumption essentially relates to the
causal structure in the later stages of a gravitational collapse.
Here, the occurrence of naked singularities arising in the gravitational
collapse of massive matter clouds such as stars have been discussed. The
examination of such dynamical collapse scenarios, concerned with the evo-
lution of regular initial data from a well-behaved initial value surface, imply
the following basic conclusions. First, such a naked singularity forms in the
dynamical evolution of several forms of matter, such as the collapse of in-
flowing radiation, dust, or perfect fluids. Second, a non-zero measure set
of non-spacelike trajectories, in the form of families of non-spacelike curves,
is emitted from the naked singularity, as opposed to a single null geodesic
escaping, which corresponds to a single wave front emerging. Finally, such a
singularity is physically significant in the sense that it is a powerfully strong
curvature singularity as the curvatures diverge rapidly along all the trajec-
tories meeting the naked singularity in the past. It can also be noted that
in dust collapse (see Fig. 4.7), and in similar models, these are seen to arise
from a non-zero measure set of initial data sets (Saraykar and Gate, 1999;
Mena, Takavol, and Joshi, 2000).
How seriously such a naked singularity is to be taken, is it generic and
stable, and what are the implications towards the formulation and proof of
the cosmic censorship hypothesis? It may be noted that the gravitational
collapse situations investigated so far have been spherically symmetric. Is it
192 Cosmic censorship
0.96
0.94
M3
0.92
0.9
33 34 35 36 37
E5 Naked singularity
Blackhole
Fig. 4.7 The blackhole and naked singularity phases in a dust collapse, in terms
of the initial density and velocity profiles (from Mena, Tavakol, and Joshi, 2000).
possible that the naked singularities occurring are artifacts of this assumed
symmetry? As discussed later, there are results that show that naked sin-
gularities occur in some non-spherical collapse models as well, and so their
existence need not be due to the assumed spherical symmetry only. Further-
more, as was shown by the singularity theorems, the singularities developing
in spherical situations still persist, even when small perturbations are taken
into account. It is possible that a similar situation may arise here also, and
so the detailed investigation of the spherically symmetric scenario becomes
quite important.
In any case, if the cosmic censorship is generically correct, it has to hold in
spherical symmetry as well. Hence the results here on the spherical collapse
show that only a substantially fine-tuned version of the censorship holds if
one exists.
then for this property to be stable, it would be expected that all the nearby
spacetimes would have the same property. For example, one says that the
spacetime is causal, and not admitting any closed non-spacelike curves is
a stable property, if for all nearby metrics ḡ on the same manifold, the
spacetime (M, ḡ ) also remains causal. Similarly, for a spacetime having naked
singularities, one would like to check if these naked singularities are stable. If
(M, g ) has naked singularities, but all nearby (M, ḡ ) have none, then these
are not stable.
The issue then is to make this notion of ‘nearby’ metrics more precise, and
for that one has to have a topology on the space of all metrics ḡ for a given
spacetime M . One way to do this is to define a C 0 -topology, by defining the
nearby metrics to be those that are nearby in their values at the spacetime
coordinates, and the open balls can be defined by the requirement that for
a given > 0, the ball consists of all metrics ḡ , such that |gij − ḡij | < .
The issue then would be, should it be required only for the values to be
nearby, or must the derivatives also be nearby? In general relativity, the
metric tensor is required to be at least C 2 . If the first derivatives are also
required to be close, then a C 1 -topology should be used, and if the second
derivatives are required to be close, then a C 2 -topology has to be used. A
concerned property, such as the spacetime having, or not having, a naked
singularity should then be examined in such a context with the appropriate
topology defining the nearby metrics. There is no uniqueness in defining such
a topology on the space of all metrics, and a C ∞ -topology could be used on
the space of all metrics. The main difficulty in a stability and genericity
analysis in general relativity is then that of identifying a suitable physically
relevant topology for the space of all metrics, and then it has to be examined
if there are open sets with non-zero measures with a given property.
As there is no uniqueness in the above procedures, the typical approach
has been to examine perturbations of a given situation, as directed by the
physical requirements of the model. For example, in a gravitational collapse,
for the homogeneous dust collapse, one may try and see if the final state of
the collapse remains unchanged if perturbations that are inhomogeneities in
a given homogeneous profile are introduced, or one may want to perturb the
original solution to a new solution that allows for small non-zero pressures
for a given initial data from which the collapse evolves. Very many such
analyzes have been carried out for gravitational collapses, as indicated in the
previous chapter and in the discussion here.
The method used allows the genericity and stability aspects related to the
occurrence of naked singularities in a gravitational collapse to be considered.
The main motivation for studying a higher-dimensional collapse is that, while
the censorship may fail in the four-dimensional manifold of general relativity,
it can possibly be restored due to the extra physical effects arising from the
transition itself to a higher-dimensional spacetime continuum.
There have been several investigations in recent years on the spherically
symmetric collapse of dust in higher-dimensions. The recent revival of inter-
est in this problem is motivated, to an extent, by various higher-dimensional
theories, including the string paradigm to unify the forces of nature, and the
brane-world scenarios. Different specialized subcases of the general problem
of dust collapse in higher spacetime dimensions have been considered. For
example, the marginally bound case in a general spacetime dimension was
studied by Ghosh and Beesham (2001). The same case was also studied with
an added and physically motivated assumption on initial density profiles,
that the first derivatives of the initial density distribution for the collapsing
cloud must be vanishing (Sil and Chaterjee, 1994; Patil, Ghate, and Saraykar,
2001; Banerjee, Debnath, and Chakraborty, 2003). Also, the non-marginally
bound case with the geometric assumption that spacetime is self-similar, was
examined by Ghosh and Banerjee (2003) for a five-dimensional model.
These studies do provide an idea of what is possible in a gravitational
collapse as far as its endstate spectrum is concerned. It is obvious from the
discussion so far that any possible proof of the censorship must be inspired by
additional physical inputs into the current framework of thinking, with one
of these being a possible transition to higher spacetime dimensions. Any such
alternatives would be worth exploring due to the fundamental significance
of the censorship in blackhole physics. If naked singularities did develop in
realistic gravitational collapses of massive objects, they may have properties
that would be rather different from those of blackholes, both theoretically as
well as observationally, and a comparison of these two cases may prove quite
interesting.
The effect of dimensions on the final fate of the evolution of the matter
cloud that collapses from a given regular initial data is now examined. A
spherically symmetric dust collapse is considered with N ≥ 4 dimensions.
There have been suggestions that if the in-falling velocity of the matter shells
is so high that the effects of pressures are negligible, then dust may be a
good approximation in the final stages of a collapse. Dust collapse is worth
investigating in any case, as it has continued to serve as a basic paradigm in
blackhole physics.
To focus the discussion, consider a model initial density profile given by
ρ2 ρ3
ρ(ti , r) = ρ0 + rρ1 + r2 + r3 + · · · , (4.76)
2! 3!
4.6 Genericity and stability 197
∞
2ρn
M(r) = Mn rn , Mn = , (4.77)
n=0
(N − 2)(N + n − 1)n!
b2
b(r) = b0 + rb1 + r2 + ··· . (4.78)
2!
First, consider the marginally bound class of collapse models for a trans-
parent understanding of the problem. This is the case when the energy
function b(r) above vanishes identically for the collapsing shells. In this case,
the first non-vanishing coefficient Xn (0), where n > 0, could be worked out
as discussed in Chapter 3. This is given by
n! Mn
Xn (0) = − . (4.79)
N −1 M0
3/2
2/(N −1)
N − 1
lim lim v = M0 X2 (0) r4/(N −1) . (4.80)
t→ts r→0 4
198 Cosmic censorship
Also, in the same limit, the function F/RN −3 has the form
F r2 M0
lim lim = . (4.81)
t→ts r→0 RN −3 v (N −3)
Therefore, it is clear from (4.80) and (4.81) that if N > 5, then for the
limt→ts , limr→0 , F/R → ∞ and thus the endstate of collapse will always
be a blackhole, as discussed previously. It follows that for a marginally
bound dust collapse, with ρ1 = 0, that is, when the initial density pro-
file is sufficiently differentiable and smooth, the cosmic censorship holds in a
higher-dimensional spacetime with N = 6, or higher.
In the above, the spacetime dimension was taken to be six or higher. Now
consider the case when the spacetime dimension is five, but still with an
analytic initial density profile. In the case of a five-dimensional marginally
bound collapse with ρ1 = 0, the tangent to the outgoing radial null geodesic
at the singularity in the (R, u) plane can be written as
1 − F/R2
x20 = M0 X2 (0) . (4.82)
1 + F/R2
The sufficient condition for the existence of an outgoing null geodesic from
the singularity is that x0 > 0, which in the above case amounts to
M2
ξ≡ < −2. (4.83)
M20
But again, the outgoing null geodesic should be within the spacetime, that
is, the slope of the geodesic must be less than that of the singularity curve,
dt dt
lim lim ≤ . (4.84)
t→ts r→0 dr null dr sing
M2
ξ= ≤ −8. (4.85)
M20
Therefore, from (4.83) and (4.85), it can be seen that, for an outgoing null
geodesic from the singularity to exist, ξ ≤ ξc = −8, in which case the result
is a naked singularity, otherwise a blackhole results as the collapse endstate.
Note that this situation has an interesting parallel to the four-dimensional
collapse scenario, where a similar critical value exists. However, it is for
the coefficient ρ3 , when both ρ1 and ρ2 are vanishing. Therefore, with the
increase of the spacetime dimension by one, the criticality separating the
4.6 Genericity and stability 199
blackhole and naked singularity phases shifts at the level of the second density
derivative from the earlier third density derivative.
An interesting observation that could be made here is that for ξ < −2
there is an increasing apparent horizon at the singularity. The apparent
horizon is given by R = F , so it initiates at the central singularity r =
R = 0 and, in the above case, it is increasing in time (as opposed to the
Oppenheimer–Snyder case of a homogeneous dust collapse). Therefore, for
the range −2 > ξ > −8, no trapped surface is formed before the singularity
epoch, but there is still a blackhole as the collapse endstate. This confirms
that the absence of a trapped surface until before the singularity is necessary,
but not a sufficient condition for the formation of a naked singularity. This
is relevant, especially for numerical collapse simulations, where the criterion
for a naked singularity formation is often taken to be just the absence of
trapped surfaces on an evolving sequence of spacelike surfaces, in a particular
slicing of the spacetime. The above example shows that the mere absence of
trapped surfaces cannot be taken as proof that the collapse terminates in a
naked singularity.
It is useful to note here that the well-known Oppenheimer–Snyder class of
collapse solutions is a special case of a marginally bound dust collapse in four
dimensions, in which case the initial density profile is homogeneous, that is,
Mn (n > 0) = 0 for all n. The point is, if the initial density is homogeneous,
but if the collapse is not marginally bound, then the non-zero energy function
f could inhomogenize the collapse at later epochs. In the present case, as
f = 0, at all later epochs the density also remains a function of time only,
that is, it is homogeneous at all later times as well, and it can clearly be
seen that the final outcome of this class of collapse is always a blackhole.
Furthermore, from (3.68) and (3.53) it can be seen that
2
ts0 = . (4.87)
(N − 1)M0
1/2
out, this is a function space consisting of all possible mass functions F and
energy functions f . It is worth knowing how, for example, a naked singularity
endstate would be affected when one moves from a given density and energy
profile, which gave rise to this state, to a nearby density or energy profile in
this space of all initial data. The issue of how given density and energy distri-
butions determine the final collapse state has been discussed quite extensively
in the usual four-dimensional dust collapse models, although a somewhat
different methodology was used. These results were completed to give a full
and general treatment of the four-dimensional case, and the typical result
is that given any density profile, one could choose the energy profile (and
vice versa), so that the collapse endstate would be either a blackhole or a
naked singularity, depending on this choice.
As can be seen from the considerations here, these results are generalized to
the case of a higher-dimensional collapse situation, and the method allows a
more definite statement on the genericity of naked singularity formation to be
made. As seen from the discussion above and in Chapter 3, the quantity X (0)
is fully determined from the initial data functions and their first derivatives.
Once it is positive, the collapse ends in a naked singularity and a negative
value gives the blackhole final state. It follows by continuity that, given a
density profile, if the energy profile chosen is such that the collapse ends in a
naked singularity, that is if X (0) > 0, then there is a whole family of nearby
velocities such that this will continue to be the case, and then the naked
singularity forms an open subspace in the initial data space. The same of
course holds for blackhole formations, and both these are neatly separated
open regions in the initial data space. But if, on physical grounds, it is taken
that both ρ1 and b1 must vanish, then the dust treatment gives X (0) = 0,
and the cosmic censor may be restored.
It may be argued that if all the assumptions such as those discussed above
can be suitably motivated physically, then it may be possible to restore
cosmic censorship in a higher-dimensional spacetime for the gravitational
collapse of dust. These conditions will now be discussed in some detail. That
the equation of state must be dust-like in the final phases of the collapse is
a strong assumption, but it is not a possibility that can quite be ruled out,
as discussed earlier. After all, very little is known on the equations of state,
especially what it would be like in the advanced stages of the collapse. Also,
it is not ruled out that, in the late stages of the collapse, the configuration
is like a marginally bound one, especially in the vicinity of the singularity.
The introduction of pressures may, or may not, change such a scenario.
In this case then, it may be possible to recover the censorship if one moves
to a higher-dimensional spacetime arena. This is subject to the validity of
several extra physical inputs, as described above. On the other hand, once
more general situations of either a non-marginally bound case, or with a
more general form of matter, or without restrictive extra assumptions on
202 Cosmic censorship
the nature of the initial density profiles are moved to, then generically both
the blackhole and naked singularity phases could result as endstates of the
collapse in a higher-dimensional spacetime scenario as well. In this way, a
dynamical collapse in general relativity offers a rich spectrum of possibilities
to investigate.
exists in the sense that the solution exists for arbitrary large times for an
asymptotic observer. For a discussion of such an approach in the context
of self-gravitating scalar fields, see Malec (1995). The problem of the global
existence of solutions is discussed by Malec, and an explicit example of an
initial configuration that results in a naked singularity is found at the center
of symmetry.
Scalar field collapse has also been numerically studied, as mentioned above.
A family of scalar field solutions was considered where a parameter p char-
acterized the strength of the field. The numerical calculations showed that,
for blackhole formation, there is a critical limit p → p∗ and the mass of the
resulting blackholes satisfy a power law Mbh ∝ ( p − p∗ )γ , where the critical
exponent γ has a value of about 0.37. It was then conjectured that such
a critical behavior may be a general property of a gravitational collapse,
because similar behavior was found in some other cases, including imploding
axisymmetric gravitational waves. Also, the case of the collapse of radiation
with an equation of state p = ρ/3 was considered, assuming self-similarity
for the solutions. It is still not clear if the critical parameter γ will have the
same value for all forms of matter chosen, and further investigations may be
required to determine this issue. As the parameter p moves from the weak
to the strong range, very small mass blackholes can form. This has relevance
to the censorship because, in such a case, one can probe and receive mes-
sages from arbitrarily near to the singularity, and this is naked singularity
like behavior. Attempts have also been made to construct models analyti-
cally that may reproduce such a critical behavior assuming self-similarity,
and solutions were constructed that have dispersal, together with solutions
with blackholes or naked singularities.
The causal structure of spacetime near a singularity and the nature of tra-
jectories emerging from it are analyzed below in some detail. If null and other
timelike paths also emerge from the singularity, then in principle, particle
and energy emission from such ultra-dense regions is allowed. Such emissions
are basically governed by the nature of non-spacelike paths near the singu-
larity. These trajectories are examined and it is shown that if a null geodesic
emerges, then families of future directed non-spacelike curves that also neces-
sarily escape from the naked singularity exist. The existence of such families
is crucial to the physical visibility of the ultra-dense regions.
Here, no underlying symmetries are assumed for the spacetime, and some
earlier considerations on the nature of causal trajectories emerging from a
naked singularity are generalized and clarified. Singularities are the regions
where the physical conditions such as densities and curvatures are at their
extreme. While the big-bang singularity of cosmology is visible in principle,
and gave rise to the universe as a whole, it cannot actually be seen. On
the other hand, when a massive star dies and collapses continually under
gravity, the eventual spacetime singularity can be either hidden within an
event horizon of a blackhole, or it could be visible to outside observers,
depending on how the collapse of the cloud evolves. While a naked singularity
forming in the collapse could provide an opportunity for the physical effects
taking place in these extreme regions to be observable to outside observers
in the universe, the actual visibility of such extreme gravity regions will
depend on the nature and structure of non-spacelike paths emerging from
the singularity.
If a continual collapse leads to a naked singularity formation then, even if
quantum gravity resolves it eventually, the point is that the causal structure
of spacetime in the vicinity of the ultra-dense regions allows them to be seen
by an external observer. Therefore, the quantum effects taking place in the
regions with arbitrarily high matter densities and curvatures can be seen by
the external observers. Any physical effects emerging will be again governed
by the existence of families of non-spacelike paths from the vicinity of the
singularity. It is therefore important to understand the structure of such
families within a gravitational collapse framework.
The important physical issue then is whether such a naked singularity
forming in the gravitational collapse could radiate away energy and par-
ticles. This depends crucially on the existence and structure of families of
non-spacelike trajectories emerging from its vicinity. Also, the actual physical
appearance and size of the singularity will be determined by the non-radial
null trajectories, and the energy emission, if any, will be governed by the
timelike curves and other non-spacelike trajectories escaping from the sin-
gularity. For this reason, several authors have considered the possibility of
non-radial null geodesics emerging from a naked singularity in the context
of spherically symmetric dust collapse models (Deshingkar and Joshi, 2001;
4.6 Genericity and stability 205
Mena and Nolan, 2001, 2002; Deshingkar, Joshi, and Dwivedi, 2002). Also,
families of non-spacelike and timelike geodesics have been worked out in
self-similar perfect fluid collapses (Joshi and Dwivedi, 1992, 1993b), and for
Vaidya radiation collapse models (Dwivedi and Joshi, 1989, 1991). Most of
these considerations have been in the framework of spherically symmetric
spacetimes, at times together with other symmetry conditions, such as self-
similarity of the models imposed, and within the framework of a specific
matter model.
A general consideration of the nature of non-spacelike trajectories near a
naked singularity will be of much interest from such a perspective. Here, the
non-spacelike trajectories from a naked singularity in general are examined,
and it is shown that if a radial null geodesic emerges, then large families of
non-spacelike curves also necessarily emerge from the singularity. It is thus
seen that the existence of a radial null geodesic is sufficient to ensure the
existence of families of timelike and non-spacelike trajectories escaping, and
in this sense a single photon escaping in a radial direction from the singu-
larity is never an isolated phenomenon. This generalizes and clarifies earlier
considerations in this direction, without assuming any symmetry conditions
on the underlying spacetime or assuming a specific matter model for the
collapse.
When the collapse ends in a naked singularity, the causal structure near the
singularity is such that a null geodesic trajectory γ emerges from it, as shown
earlier. Specifically, γ is future directed, which in the past terminates at the
singularity, and is therefore a past incomplete null geodesic. To examine in
general the possible existence and nature of non-spacelike curves emerging
from this naked singularity, the causal boundary construction developed by
Geroch, Kronheimer, and Penrose (1972), where the spacetime M is taken
to satisfy a suitable causality condition such as strong causality, which rules
out the existence of closed timelike curves, is used here. In this procedure,
a boundary is attached to the regular spacetime manifold, which includes
spacetime singularities as well as the points at infinity.
Note that a boundary attachment to the spacetime manifold is essential to
treat the regular spacetime events, together with its singularities and points
at infinity in a unified manner. There are different ways to attach a boundary
to the spacetime, and they do not necessarily all give the same result. Here
the rather basic approach, given by Geroch, Kronheimer, and Penrose (1972),
is used as it depends essentially only on the causal structure of the spacetime,
which is much more fundamental than, for example, the differential structure
of the spacetime manifold. Also, from a physical point of view, each ideal
point here is directly associated with the region of spacetime that it can
influence, or that it would be influenced by.
An open set W in the spacetime is called a future set if it contains its
own future, that is, I + (W ) ⊂ W . Furthermore, a future set W is called an
206 Cosmic censorship
depending on the nature of the initial data from which it evolves and the
dynamical evolutions of the collapsing cloud, as allowed by the Einstein
equations.
However, as remarked above, if the visibility and other related physical
characteristics of a naked singularity that formed in the gravitational collapse
are to be explored, then it is important to examine and understand the
structure of the families of non-spacelike curves from the singularity. Again,
if non-spacelike curves emerge from the singularity, but do not go out of the
boundary of the collapsing cloud, then the singularity will be only locally
visible, and outside observers would not be able to see it. It is necessary
therefore to understand the structure of non-spacelike curves from a naked
singularity in general.
It is now possible to do this. Since the set S = I + (γ ) is a TIF, it fol-
lows from Geroch, Kronheimer, and Penrose (1972) that in this case, a past
inextendible timelike curve λ must exist, such that S = I + (λ). In the case
of the collapse ending in a naked singularity and a radial null geodesic γ
escaping from it, it can thus be seen that the set I + (γ ) is a TIF, and so, by
the above result, there is a timelike curve λ generating this TIF, in the sense
that S = I + (λ). Since both the non-spacelike trajectories γ and λ represent
the same ideal or boundary point of the spacetime that is the naked singular-
ity, and since I + (γ ) = I + (λ) by definition, it follows that the future directed
timelike curve λ must terminate in the past at the naked singularity. In other
words, it has been shown that a timelike curve λ, which escapes away to the
future, and which terminates in the past at the naked singularity exists.
It follows that if p ∈ λ and q is any other event such that q ∈ I + ( p),
then there are timelike curves from the naked singularity to q . This proves
the existence of families of infinitely many future directed non-spacelike tra-
jectories escaping away from the naked singularity. In general, if λ is any
other future directed non-spacelike curve such that I + (λ ) = I + (λ), then it
follows that they all represent the same TIF, which is the naked singularity,
and that λ terminates in the past at this singularity.
Therefore, it can be seen that there is an infinity of future going non-
spacelike curves that emerge from the singularity, if a single null geodesic has
emerged. These include timelike curves as well as non-radial non-spacelike
geodesics. It is seen that the usual method employed to show the existence
of a naked singularity in the collapse, which establishes the existence of a
radial null geodesic escaping away, is sufficient to lead to the existence of
infinite families of future going non-spacelike curves from the naked singu-
larity, as shown here. In the present consideration, it is no longer required to
have any special symmetry assumptions on the spacetime, such as spherical
symmetry, self-similarity, or others, or any specific form of matter model
such as the dust equation of state, which are usually assumed in such
discussions.
208 Cosmic censorship
In particular, this also clarifies and generalizes the earlier results on dust
collapses and other models mentioned above, which have focused on non-
spacelike null geodesics. The null geodesics of the spacetime have, of course,
a special role to play as far as the visibility of the singularity is concerned.
From such a perspective, the existence of radial versus the non-radial families
of null geodesics from a naked singularity will be briefly discussed. Suppose
a radial null geodesic emerges from the naked singularity S developing in a
continual collapse. In this case, as seen above, a timelike curve λ generating
the TIF set I + (λ) that represents the boundary point S exists. All other
future directed non-spacelike curves γ that satisfy I + (λ) = I + (γ ) generate
the same TIF representing the boundary point S , and they give the families
of particle or photon trajectories escaping away from the naked singularity.
The boundary of this future set, which is a TIF, is a three-dimensional null
hypersurface that is ruled by the radial, as well as non-radial, null geodesics
generators γ , which are all incomplete when extended in the past, and which
all have the property that I + (λ) = I + (γ ). This shows that the existence
of a radial null geodesic is sufficient also to give families of non-radial null
geodesics emerging from the singularity. This generalizes the earlier results
on the existence of non-radial null geodesics from the singularity for the
spherically symmetric dust collapse, when a radial null geodesic emerges
from the naked singularity.
It is seen that once a singularity is naked, it gives rise to infinitely many null
as well as timelike curves to escape away from. In this sense, the emission of
paths representing particle or photon trajectories from the naked singularity
is a generic phenomena. This is essential and is a necessary condition for
the naked singularity to give rise to any physical effects that may possibly
be observed by external observers. In the present consideration, the global
visibility of the singularity, that is, the situation in which once the families
of non-spacelike curves emerge from the naked singularity when they will
actually cross the boundary of the cloud to escape to an outside observer,
has not been discussed. It is known, however, in several cases including
spherical dust collapses, that whenever a singularity is locally naked, then
the rest of the free functions in the model can be chosen so as to make it
globally visible.
It is known, for example, in the case of dust collapses, that once the
singularity is locally naked, the choice of a suitable behavior of the mass
function (which is a free function, subject only to some physical conditions
such as an energy condition and regularity of the initial data) away from the
center, allows the null rays to emerge from the boundary of the cloud. It
may also be noted that in various classes of self-similar collapses, once the
singularity is locally naked it becomes necessarily globally visible (Joshi and
Dwivedi, 1992, 1993b). In any case, as there is no scale in the problem, once
the singularity is locally visible, an observer within a large enough blackhole
4.6 Genericity and stability 209
will still be able to see it for a long enough time. In such a scenario, the
escape of rays outside the boundary of the cloud would not be crucial. A
discussion on global visibility in a more general context of perfect fluids has
recently been given by Giambo (2006).
5
Final fate of a massive star
The considerations on gravitational collapse so far have been with the moti-
vation to address the physical questions such as the role of collapse in
astrophysics and cosmology. Many of the cosmic processes, such as the
birth of stars, the formation of galaxies, and others, are not well under-
stood today, but it is clear that gravitational collapse will play a major role
there. Hence, understanding the dynamics of the collapse is important, as
has been attempted here in various cases.
The important question of the final fate of massive stars at the end of
their life cycle, when they have used all their nuclear fuel, and when grav-
ity becomes the sole and key governing force, has drawn much attention
for many decades. The importance of this issue was highlighted by Chan-
drasekhar (1934), who pointed out that the life history of a star of small
mass must be essentially different from that of a star of large mass, and that
while a small mass star can pass into a white dwarf stage, a star of large mass
cannot go to this state, and one is left speculating on other possibilities. The
question as to what happens when a massive star, heavier than a few solar
masses, collapses under its own gravity has been a fundamental key problem
in astronomy and astrophysics. If the star is sufficiently massive, beyond
the white dwarf or neutron star mass limits, then a continued gravitational
collapse must ensue without achieving any equilibrium state, when the star
has exhausted its nuclear fuel. To understand the possible endstates of such
a continual gravitational collapse, the dynamical collapse scenarios must be
studied within the framework of a gravitation theory such as Einstein’s the-
ory. The theory of singularities discussed earlier then implies that, under
rather general physical conditions, a spacetime singularity must develop.
Considerations here show that, according to Einstein’s theory of gravity,
such a star in a continual collapse can end up either in a blackhole final state,
from which no communications in the form of light or particles would come
out, or it can go to a naked singularity, where the collapse proceeds but the
210
Final fate of a massive star 211
trapped surfaces do not form early enough to cover the singularity. In such
a case, the extreme gravity regions can communicate and send out physical
effects to the external universe, and it is also suggested that a huge amount of
energy could possibly be released, in principle, during the final stages of the
collapse from the regions close to the classical singularity (Joshi, Dadhich,
and Maartens, 2000). According to these suggestions, an enormous amount of
energy may probably be generated, either by some kind of a quantum gravity
mechanism, or by means of an astrophysical process where the region could
simply turn into a fireball-like situation, creating shocks into the surround-
ing medium. One way to study the structure of these extreme regions is to
examine the complete spectrum of non-spacelike geodesics through which
this energy could escape. Even if a fraction of the energy so generated is able
to escape to a distant observer, it cannot be ruled out that an observational
signature may be generated. It therefore becomes important to look into
these possibilities in some detail, and to consider likely observable signatures
of each of the blackhole and naked singularity final states of the collapse,
from the perspective of a faraway observer.
The theoretical properties and possible observational signatures of a
blackhole and naked singularity would be significantly different from each
other, and this could be of potential interest from the perspective of astro-
physical observations. An immediate distinction is, for example, in the case
where the collapse ends in a blackhole, an event horizon develops well before
the occurrence of the singularity, and thus the regions of extreme physical
conditions are always hidden from the outside world. But if the collapse devel-
oped into a globally naked singularity, then the energy of the region neigh-
bouring the singularity can escape via the available non-spacelike geodesics
paths or via other non-geodetic, non-spacelike trajectories to a distant
observer.
In this chapter, various aspects on gravitational collapse that have emerged
from the analysis so far are considered. An important issue is the final fate
of a non-spherical collapse. In earlier treatment, it was shown that strong
curvature naked singularities arise in a variety of situations involving dust,
perfect fluids, and other forms of matter. There are many interesting ques-
tions that are under active investigation at the moment. For example, could
naked singularities generate bursts of gravity waves? What kind of quantum
effects will take place near a naked singularity? What will be the generic
outcome for the case of a non-spherical collapse? Many of these issues could
have interesting physical implications. The possibility that the ultra-high
energy astrophysical phenomena, such as the gamma ray bursts, may have a
strong connection to the physics and dynamics of the gravitational collapse
of massive stars cannot be ruled out. In fact, many of the current gamma
ray burst models involve a collapsar, emphasizing the role of the massive
star collapse. Another intriguing possibility is that a naked singularity may
212 Final fate of a massive star
possibly provide some kind of observational signatures for the quantum grav-
ity effects taking place in ultra-strong gravity regions. This would then be
an exciting prospect in view of the current lack of knowledge on quantum
gravity.
In Section 5.1, the life cycle of massive stars is discussed. The large mass
stars, several times heavier than the Sun, follow a characteristically different
life cycle from stars of about one solar mass. Such massive stars typically
live a much shorter life than small mass stars, and undergo a catastrophic
gravitational collapse at the end of their life cycle. In Section 5.2, how a
physically realistic collapse must evolve and the different dynamical forces
at work to govern the final fate of such a collapse are discussed. To give an
explicit example, it is shown how the spacetime shearing forces, caused by the
inhomogeneities in the matter distribution of the star, delay the formation of
trapped surfaces during the evolution of the collapse. This allows the ultra-
strong regions of gravity and the spacetime singularity to become visible
to faraway observers, and it is seen how a naked singularity forms rather
naturally as a collapse endstate.
For any realistic consideration of gravitational collapse, departures from
spherical symmetry have to be taken into account. While this problem is
much more complex and not much work is still carried out on this, the Szek-
eres collapse models are discussed in Section 5.3. These are not spherically
symmetric, and have no Killing vectors. It is shown that these again admit
both blackhole and naked singularity final states. While much analytic and
numerical work remains to be carried out on non-spherical collapse, this
model explicitly illustrates and shows that the naked singularity final state
for the collapse is not necessarily limited to spherical symmetry, and that
these can result in a non-spherical collapse as well.
If both blackhole and naked singularity final states do occur in physi-
cally realistic gravitational collapse scenarios, the main issue would then be
the outcome that nature may prefer for the final state of a collapsing star.
While, mathematically, many evolutions may be possible, in physically real-
istic situations only certain evolutions may be actually realized. From such a
perspective, in Section 5.4 the paradoxes that are associated with the black-
hole formation as the collapse final state, when a massive star undergoes
gravitational collapse at the end of its life cycle, are discussed. Finally, in
Section 5.5 it is discussed how the quantum effects, which will be prominent
in the late stages of the collapse, can play a crucial role in creating a massive
emission of mass and energy in the very late stages of the collapse. It could
then be asked if a ‘quantum star’ comes into being in the very final stages of
the collapse, which allows for a powerful emission of mass and energy. Such
a burst-like emission can have observational features on the one hand and,
on the other, it would help resolve the naked singularity formation, thus
allowing the classical singularity to be removed.
5.1 Life cycle of massive stars 213
The equation of state for an ideal electron Fermi gas was approximated
by Chandrasekhar (1931, 1934), who showed that there is a maximum mass
limit for the mass of a spherical, non-rotating star to achieve a white dwarf
stable state. This is given by
2
2
Mc ∼ 1.4 M , (5.1)
µe
In any case, it does not seem likely that all such massive stars would be
able to throw away almost all, or a very major part, of their mass in such
a process. That is because, for stars having tens of solar masses, this would
amount to throwing away almost ninety percent of the mass of the star. In
a supernova, typically only the outer layers of the star are blown off, and
no suitable mechanism that could achieve such a high degree of efficiency
for the mass ejection from the star is envisaged as yet. In the case of a
massive star, in the course of its normal life of nuclear burning, not a very
large portion of its mass can be removed through radiation. Then, once the
gravitational collapse initiates, which is a catastrophic process, a star that
has lived millions of years collapses gravitationally in a matter of seconds.
Now, during such a catastrophic collapse, if the shock that is produced could
not blow off almost all the outer layers, these layers would fall on the newborn
neutron star that has momentarily formed at the core in the process of the
collapse, and the further collapse continues again as the pull of gravity would
exceed the balancing neutron degeneracy pressures.
It is thus seen that the evolution of a massive star, when it exhausts its
nuclear fuel, causes an inevitable continual collapse at the end point of a
stellar evolution. The general relativity theory then implies that a spacetime
singularity must form of necessity in such a scenario. The basic ingredients
of such spacetime singularities were discussed earlier. The cosmic censorship
conjecture then asserts that any such singularities forming in gravitational
collapse must be covered necessarily within an event horizon of gravity, invis-
ible to any external observer, and that in general the final state of such an
evolution must then be a blackhole. This requires an appropriate formulation
of the concepts of a blackhole as well as the cosmic censorship, within the
framework of an appropriate spacetime geometry. To decide on this issue,
the dynamical gravitational collapse has to be analyzed to see how the struc-
ture of the horizons evolve. It is within such a perspective that the analysis
here on gravitational collapse has taken place, and it can be seen that the
dynamical gravitational collapse processes are entirely fundamental in this
way to the basics and applications of blackhole physics.
permitted by general relativity, nature might select only a few of these for
the actual evolution of a massive star. Therefore, further examination of the
physical forces operating within a star as the collapse proceeds is required.
This issue is discussed here in the context of a spherical collapse by means
of an explicit example. This also throws some light on why naked singularities
form at all in a gravitational collapse. The key physical features that possibly
cause the development of a naked singularity, rather than a blackhole, as the
endstate of a gravitational collapse are investigated, and it is seen that suffi-
ciently strong shearing effects near the singularity can delay the formation of
the apparent horizon. This exposes the singularity to an external observer,
in contrast to a blackhole, where it is hidden behind an event horizon due to
the early formation of the trapped surfaces. The final outcome of a gravita-
tional collapse in general relativity is an issue of much importance from the
perspective of blackhole physics and its astrophysical implications, and one
needs to understand the key physical characteristics and dynamical features
in the collapse that give rise to a naked singularity, rather than a blackhole.
Here, the treatment of Joshi, Dadhich, and Maartens (2002) is followed.
It is seen that it is the inhomogeneity and related shearing effects within
the cloud that, if sufficiently strong near the central worldline of the col-
lapsing cloud, would delay the formation of the apparent horizon so that
the singularity becomes visible, and communication from the extreme strong
gravity regions to outside observers becomes possible. When the inhomo-
geneity and related shear forces are weak, or for the extreme case of no shear
in a fully homogeneous collapse, the collapse necessarily ends in a blackhole,
because an early formation of the apparent horizon leads to the singularity
being hidden inside an event horizon.
For the spherical gravitational collapse of a massive matter cloud, the
interior metric in comoving coordinates is
where 4πF (ti , rb ) = M , which is the total mass of the collapsing cloud, and
where r > rb is a Schwarzschild spacetime. As earlier, the rescaling freedom
in r is used to set
R(ti , r) = r , (5.5)
so the physical area radius R increases monotonically in r, and with Ri = 1,
there are no shell-crossings on the initial surface from which the collapse
develops. Of interest here is the central shell-focusing singularity at R =
0, r = 0, which is a gravitationally strong singularity, as opposed to the
shell-crossing singularities that are weak, and through which the spacetime
may sometimes be extended.
The evolution of the density and radial pressure are given by
F Ḟ
ρ= , pr = . (5.6)
R2 R R2 Ṙ
The central singularity at r = 0, where density and curvature are infinite,
is naked if there are outgoing non-spacelike geodesics that reach outside
observers in the future and terminate at the singularity in the past. Outgoing
radial null geodesics of the metric (5.2) are given by
dt
= eψ−ν . (5.7)
dr
First consider the case of a homogeneous density collapse, ρ = ρ(t). Writ-
ing f = e−2ψ R2 − 1, the Einstein equations give f − e−2ν Ṙ2 = −F/R.
Then, as discussed in Chapter 3, the geodesic equations can be written as
!
dR f + F/R R
= 1− , (5.8)
du 1+f αrα−1
where u = rα (α > 1). If there are outgoing radial null geodesics terminating
in the past at the singularity with a definite tangent, then at the singularity,
dR/du > 0. Now in the case of homogeneous density, the entire mass of the
cloud collapses to the singularity simultaneously at the event (t = ts , r = 0),
so that F/R → ∞. Then in that case, from (5.8), dR/du → −∞, so that
no radial null geodesics can emerge from the central singularity. It can be
shown that all the later epochs t > ts are similarly covered.
It has thus been shown that for a spherical gravitational collapse with
homogeneous density, the final outcome is necessarily a blackhole. The
pressures, however, can be arbitrary on which no conditions are imposed.
This conclusion does not require homogeneity of the pressures pr and
pθ , and is independent of their behavior. This generalizes the well-known
Oppenheimer–Snyder–Datt result of the special case of dust, where the
218 Final fate of a massive star
As the density grows without bound, trapped surfaces develop within the col-
lapsing cloud. These can be traced explicitly via the outgoing null geodesics,
and the equation of the apparent horizon, t = tah (r), which marks the bound-
ary of the trapped region and is given by (see also the discussion on apparent
horizons in Section 3.6)
If the apparent horizon starts developing earlier than the epoch of the sin-
gularity formation, then the event horizon can fully cover the strong gravity
regions including the final singularity, which will thus be hidden within a
5.2 Evolution of a physically realistic collapse 219
blackhole. On the other hand, if trapped surfaces form sufficiently later dur-
ing the evolution of the collapse, then it is possible for the singularity to
communicate with outside observers.
For the sake of clarity, consider a marginally bound collapse, f = 0,
although the conclusions can be generalized to hold for the general case.
Then, (5.9) can be integrated to give
3
R3/2 (t, r) = r3/2 − (t − ti )F 1/2 (r), (5.14)
2
and (5.12) and (5.13) lead to
3 1/2
2 r
ts (r) = ti + , (5.15)
3 F (r)
2
tah (r) = ts (r) − F (r). (5.16)
3
The central singularity at r = 0 appears at the time
2
t0 = ts (0) = ti + √ , (5.17)
3ρc
1 r
(3F − rF ) .
2
σ 2 ≡ σab σ ab = 2 (5.18)
2 4
6R R F
Here, a generic inhomogeneous mass profile can be chosen to have the form
F (r) = F0 r3 + F1 r4 + F2 r5 + · · · (5.19)
The time dependent factor in square brackets on the right of (5.21) decreases
monotonically from 1 at t = ti to 0 at t = t0 . Therefore, the qualitative role
of the shear in singularity formation can be seen by looking at the initial
shear. The initial shear σi = σ (ti , r) on the surface t = ti grows as rn ,
n ≥ 1, near r = 0. A dimensionless and covariant measure of the shear is
the relative shear, |σ/Θ|, where
Ṙ Ṙ
Θ=2 + (5.23)
R R
is the volume expansion. It follows that
σ
−nFn n
= √ r [1 + O(r)]. (5.24)
Θ i 3 6F 0
It can now be seen how such an initial shear distribution determines the
growth and evolution of the trapped surfaces, as prescribed by the apparent
horizon curve tah (r), given by (5.22). If the initial density profile is assumed
to be smooth at the center, then ρi (r) = ρc + ρ2 r2 + · · · , with ρ2 ≤ 0,
which corresponds to F (r) = F0 r3 + F3 r5 + · · · , with F2 ≤ 0. Now suppose
that ρ2 (and hence F2 ) is non-zero. Then, (5.22) implies that the apparent
horizon curve initiates at r = 0 at the epoch t0 , and increases near r = 0
with increasing r, moving to the future (see Fig. 5.1). Note that as soon
as F2 is non-zero, even with a very small magnitude, the behavior of the
apparent horizon changes qualitatively. Rather than going back into the past
from the center, as would happen in the homogeneous collapse case with
F2 = 0, it is now future pointed. This leads to a locally naked singularity as
the collapse endstate. The singularity may be globally naked and visible to
faraway observers, depending on the nature of the density function at large
values of r.
5.2 Evolution of a physically realistic collapse 221
1.2e-06
1e-06
8e-07
F3 = –30
6e-07
tah
4e-07
F3 = –15
2e-07
F3 = –2
0
F3 = –0.1
–2e-07
0 0.001 0.002 0.003 0.004 0.005
r
In the general case, where there is not a necessarily smooth initial density,
this condition is violated for n = 1, 2, as follows from (5.22). The apparent
horizon curve initiates at the singularity r = 0 at the epoch t0 , and increases
with increasing r, moving to the future, i.e. tah > t0 for r > 0 near the cen-
ter. The behavior of the outgoing families of null geodesics has been analyzed
in detail in these cases, and it is known that the geodesics terminate at the
singularity in the past, which results in a naked singularity, as discussed in
Section 3.6. In such cases, the extreme strong gravity regions can commu-
nicate with outside observers. For the case where n = 3, (5.25) shows that
it is possible to have a blackhole if F3 ≥ −2F0 , or a naked singularity if
5/2
When the dust density is homogeneous, the apparent horizon starts devel-
oping earlier than the epoch of the singularity formation, and the singularity
then is fully hidden within a blackhole. There is no density gradient at the
center, and no shear present within the cloud. On the other hand, if a den-
sity gradient is present at the center, then the trapped surface development
is delayed due to the presence of the shear, as seen above, and, depending
on the ‘strength’ of the density gradient and shear at the center, this may
allow the singularity to be visible. It is the rate of decrease of the shear, as
222 Final fate of a massive star
Db Eab = 13 Da ρ , (5.26)
where Da is the covariant spatial derivative. Note that the magnetic Weyl
tensor vanishes for spherical symmetry. In turn, the gravito-electric field is
a source for shear, or equivalently, the shear is a gravito-electric potential,
2 2
uc ∇c σab + Θσab + σac σ c b − σ 2 hab = −Eab . (5.27)
3 3
Da ρ = −4σDa σ − 2ΘDb σab − 3Db (uc ∇c σab ) − 3σa b Dc σbc − 3Db (σac σ c b ),
(5.28)
5.2 Evolution of a physically realistic collapse 223
where the shear constraint Db σab = 2/3Da Θ has been used. Equation (5.28)
explicitly links the behavior of the density gradients and the shear near the
center, which was discussed above. The free gravitational field that mediates
this link can also provide a covariant characterization of the singularity for-
mation. By (5.24) and (5.27), the relative gravito-electric field E/Θ2 (where
E 2 = 1/2E ab Eab ) near r = 0 is given at t = ti by
E −7nFn n
= √ r [1 + O(r)] . (5.29)
Θ2 i 18 6F0
then (5.6) shows that ρi and (pr )i are constants. The density and pressure
may however develop inhomogeneities as the collapse proceeds, depending
on the choice of the remaining functions, including in particular the initial
velocities of the collapsing shells. The collapse may then end up in either a
blackhole or a naked singularity (for a discussion on this for the case of a
dust collapse, see Joshi and Dwivedi, 1993a). In fact, it can be shown that
zero shear implies a blackhole for these models. By (5.3), (5.6), and (5.30),
224 Final fate of a massive star
the shear-free condition leads to R /R = 1/r, and (5.6) then shows that
ρ = ρ(t), i.e. the density evolution is necessarily homogeneous. As shown
above, the collapse thus necessarily ends in a blackhole. For the class of
models given by (5.30), whenever the collapse ends in a naked singularity,
the shear must necessarily be non-vanishing. Although this class of models is
somewhat special, the result indicates that the behavior of the shear remains
a crucial factor even when pressures are non-vanishing.
It would appear that the only way a singularity becomes visible is when
suitable modifications occur in the geometry of trapped surfaces and the
apparent horizon, and a delay in the trapped surface formation is necessary.
As has been shown here, the shear provides a rather natural explanation for
the occurrence of locally naked singularities in certain collapse models. Suffi-
ciently strong shearing forces in spherical collapsing dust do affect the delay
in the formation of the apparent horizon, thereby exposing the strong gravity
regions to the outside world, and leading to a naked singularity formation.
When shear decays rapidly near the singularity, the situation is effectively
like the shear-free case, with a blackhole endstate. The important point is
that naked singularities can develop in quite a natural manner, very much
within the standard framework of general relativity, as governed by shearing
effects.
In the case of a spherical dust collapse, shear and density inhomogeneity
are equivalent, that is, one implies the other. Although shear contributes
positively to the focusing effect via the Raychaudhuri equation
1 1
Θ̇ + Θ2 = − ρ − 2σ 2 , (5.31)
3 2
its dynamical action can make the collapse incoherent and dispersive.
Depending on the rate of fall-off of shear near the singularity, its disper-
sive effect can play the critical role of delaying the formation of the apparent
horizon, without directly hampering the process of collapse. The dispersive
effect of shear always tends to delay the formation of the apparent horizon,
but it is able to expose the singularity only when the shear is strong enough
near the singularity.
Here, the effects of inhomogeneities and shear are analyzed for certain
collapse models. The implication is certainly not that it is only the shear
that can cause the changes in the trapped surfaces’ geometry. The important
point to emphasize is that modifications in the trapped surfaces’ geometry
can occur rather naturally in the process of a dynamical collapse, due to
physical agencies such as inhomogeneities and spacetime shear. In general,
the phenomena such as trapped surfaces formation and apparent horizons
are independent of any spacetime symmetries, and it would seem clear that
a naked singularity will not develop in a more general gravitational collapse
5.3 Non-spherical models 225
as well, unless there is a suitable delay of the apparent horizon. This suggests
that the forces such as shear will continue to be pivotal in determining the
final fate of physically realistic gravitational collapse scenarios, independent
of any spacetime symmetries. In any case, the main purpose here has been to
try to understand and find the physical mechanism that leads the collapse to
the development of a naked singularity rather than a blackhole, in some of the
well-known classes exhibiting such behavior. It is found that the shear and
inhomogeneity provide a covariant dynamical explanation of the phenomenon
of naked singularity formation in a spherical gravitational collapse.
The basic question is: what governs the geometry of the trapped surfaces,
or the formation, or otherwise, of the naked singularities in gravitational
collapse? In other words, what is it that causes the naked singularity rather
than a blackhole to develop as the final end product of the collapse? It turns
out from the above that physical agencies such as inhomogeneities in matter
profiles, as well as spacetime shear, play a key role in distorting the trapped
surface geometry, and could delay the trapped surface formation during the
collapse, thus giving rise to a naked singularity, rather than a blackhole, as
the collapse final state. This, in a way, provides the physical understanding of
the phenomena of blackhole and naked singularity final states in the collapse.
It is seen that in a spherical collapse, these phases are generic, and are seen
to be determined by the nature of the initial data from which the collapse
develops, and in terms of the allowed dynamical evolutions. Physical agencies
such as inhomogeneities and shear may cause them, and given the initial
data, there are non-zero measure classes of evolutions that evolve into either
of these outcomes.
C ≤ 2π (2GM/c2 ). (5.32)
5.3 Non-spherical models 227
Therefore, unlike the cosmic censorship conjecture, the hoop conjecture now
no longer rules out all the naked singularities, but only makes a definite
assertion on the occurrence of the event horizons in the gravitational collapse.
It is thus seen that the hoop conjecture allows for the occurrence of naked
singularities in general relativity, at least when the collapse is sufficiently
aspherical, and especially when one or two dimensions are sufficiently larger
than the others. A known example of this is given by Lin, Mestel, and Shu
(1965), which discusses gravitational collapse for uniform spheroidal objects
from the perspective of instability in Newtonian gravity (see also, Thorne,
1972; Shapiro and Teukolsky, 1991, 1992, for a discussion). Here, a non-
rotating homogeneous spheroid collapses, maintaining its homogeneity and
spheroidicity, but its deformations grow as the collapse progresses. If the
initial condition is that of a slightly oblate spheroid, the collapse results in a
pancake singularity through which the evolution could proceed and continue.
However, for a slightly prolate spheroidal configuration, the matter collapses
to a thin thread that ultimately results in a spindle singularity. This is more
serious in nature in that the gravitational potential and force and the tidal
forces blow up, as opposed to only the densities blowing up, in a mere shell-
crossing singularity. Even in the case of an oblate collapse, the passing of
matter through the pancake causes prolateness, and subsequently a spindle
singularity again results without the formation of any horizon.
It was indicated by the numerical calculations of Shapiro and Teukolsky
(1991) that a similar situation is maintained in general relativity, also in con-
formity with the hoop conjecture. They evolved collissionless gas spheroids
in full general relativity that collapse in all cases to singularities. When the
spheroid is sufficiently compact, a blackhole that contains the singularity
forms, but when the semi-major axis of the spheroid is sufficiently large, a
spindle singularity results without an apparent horizon forming. These have
to be treated as numerical results, as opposed to a full analytic treatment,
which need not be in contradiction to a suitably formulated version of the
cosmic censorship. The definition of naked singularity here was basically in
terms of the non-occurrence of trapped surfaces in a certain family of non-
spacelike surfaces. Sometimes this may, or may not, indicate a genuine naked
singularity, and a blackhole may still form. However, this gives rise to the
possibility of the occurrence of naked singularities in the collapse of finite
systems in asymptotically flat spacetimes, which could be in violation of the
weak cosmic censorship, but possibly in conformity with the hoop conjecture.
A somewhat broader statement in a similar spirit to the hoop conjecture
is the event horizon conjecture of Israel (1984, 1986a, 1986b). A general
statement for this conjecture is given by the requirement that an event hori-
zon must form whenever a matter distribution (satisfying appropriate energy
conditions) has passed a certain critical point in its gravitational collapse,
namely the formation of a closed trapped surface. A strong motivation for
228 Final fate of a massive star
The TBL dust collapse spacetimes, even though they generalize the homo-
geneity assumption, are special in some ways. They are spherically symmetric
and they have matter in the form of irrotational pressureless dust. It is
interesting to know whether naked strong curvature singularities occur in
more general situations than this. Here, it is discussed how naked strong
curvature singularities occur in the Szekeres spacetimes that do not have
any Killing vectors (Joshi and Królak, 1986). This result shows that naked
strong curvature singularities do not arise necessarily as a result of the spher-
ical symmetry. Nevertheless, the Szekeres spacetimes have the same special
form of matter as the TBL spacetimes, which is again irrotational pressure-
less dust. Moreover, Szekeres spacetimes are also special because they can
be matched to the Schwarzschild spacetime and they cannot contain any
gravitational radiation (Bonnor, 1976).
The Szekeres spacetime (Szekeres, 1975) is a solution of the Einstein
equations representing irrotational dust, where
and the units have been chosen so that c = 8πG = 1. The metric has the
diagonal form given by
where (r, x, y ) are comoving spatial coordinates. The solution is given below
for the case Y = ∂Y /∂r = 0,
f (0) = 0. (5.39)
The function R0 (r) = R(r, 0), which again indicates the physical radius of
the shells, is a monotonically increasing function of r. The freedom in the
choice of the radial coordinate r can be used to do the scaling as
R0 (r) = r. (5.40)
P F − 3F P
ρ= . (5.41)
P 2 R2 Y
Although for P > 0, the surfaces r = const., t = const. are spheres, the
solution is not spherically symmetric here because the spheres are no longer
concentric in this case, and their centers are given by (−a−1 b1 , −a−1 b2 ).
Szekeres has also analyzed the singularities and their causal structure in these
spacetimes. When R = 0, the singularity is of the first kind, and when Y = 0
the singularity is of the second kind. The singularities of the second kind are
familiar shell-crossing singularities that also occur in TBL spacetimes. As in
the TBL spacetimes, the shell-crossing singularities in Szekeres spaces can
also be both locally and globally naked (Szekeres, 1975). However, they are
generally believed to be mild and they will not be considered here. These
singularies can be eliminated by imposing a regularity condition,
Y > 0. (5.42)
The analysis here can again be similar to that in the TBL collapse models
case, as given by Joshi and Dwivedi (1993a). A set of new functions can be
introduced,
R F F
X= α
, η = r , Λ = α,
r F r
1 − 3η
1
P
Θ= , L=r , (5.45)
r3(α − 1)/2 P
where
1 Θ
H = ηX + √ , (5.47)
3 X
√
Λ 2
N = − 2 Θ − ηX 3/2
. (5.48)
2X 3
The tangents K a = dxa /dk for the outgoing radial (x = const., y = const.)
null geodesics can be written as
dt P
Kt = = , (5.49)
dk Y
dr P
Kr = = , (5.50)
dk PY Y
dx
Kx = = 0, (5.51)
dk
dy
Ky = = 0, (5.52)
dk
5.3 Non-spherical models 233
where u = rα . Note that when one writes the quantity U , one is essentially
working out the necessary condition for the null geodesics to emerge from
the singularity at R = 0, r = 0, with a well-defined tangent given by dR/dt.
Therefore, if the trajectories are emerging with a well-defined tangent from
the singularity, then the quantity U (x, u) is well-defined by definition. For
more details, see, for example, Joshi and Dwivedi (1993a). Consider the limit
X0 of the function X along the null geodesic terminating at the singularity
at R = 0, u = 0. Using the l’Hospital rule,
R dR
X0 = lim = lim = lim U (X, u) = U (X0 , 0).(5.56)
R→0, u→0 u R→0, u→0 du R→0, u→0
The necessary condition for the existence of the null geodesic outgoing from
the central singularity is the existence of the positive real root X0 of the
equation
V (X ) ≡ U (X, 0) − X = 0. (5.57)
where
η0 = η (0), Λ0 = Λ(0), Θ0 = Θ(0). (5.59)
234 Final fate of a massive star
This is exactly the same equation as in the marginally bound TBL case.
Consequently, the same analysis as by the TBL case given earlier by Joshi
and Dwivedi (1993a) applies here. Here, only a few key results will be sum-
marized. To show that the singularity is naked, there is a need to prove that
a solution of the geodesic equation exists such that the tangent X0 is real-
ized at the singularity. It can be proved that there is always at least a single
null geodesic outgoing from the central singularity. Therefore, the existence
of the real and positive root of the (5.58) is both a necessary and sufficient
condition for the existence of a naked singularity.
Some further remarks will now be made on the strength of the singularity,
which will be investigated in some detail. Let M be the spacetime manifold,
and let J (k ) be the quantity along a null geodesic in the spacetime as given
below, with λ : (k0 , 0] → M , being parametrized by an affine parameter k ,
0
J (k ) = Rab K a K b dk . (5.60)
k0
interval (k0 , 0] with respect to k . Therefore the limit limk→0 K t exists, and
Kt
Kr = . (5.63)
rα−1 (H − XL)
Therefore if α = 1 the limit limk→0 K r also exists. Suppose that α = 1, and
consider the limit limk→0 kRab K a K b . Applying the l’Hospital rule twice and
using (5.53),
η0 Λ0 H0
lim kRab K a K b = lim K r . (5.64)
k→0 X02 (α − N0 )2 k→0
Therefore when α = 1 the above limit is finite and the naked singularity just
satisfies the limiting focusing condition (Rab K a K b diverges logarithmically).
If α > 1 the above limit diverges and the naked singularity also satisfies the
limiting focusing condition (but not the strong focusing condition, unless
Λ0 = 0).
The above results show that under the regularity conditions, and for a
marginally bound collapse, any central naked singularity in the Szekeres
spacetime is always a strong curvature singularity in Królak’s sense. Fur-
thermore, as the considerations above reduce to the TBL case under various
appropriate limits, the above also implies that many classes of naked sin-
gularities occurring in the TBL dust collapse models are strong curvature
naked singularities necessarily. The significance of this is that these are
physically serious curvature singularities, which are not removable from the
spacetime. As discussed earlier in Chapter 4, all volume forms along non-
spacelike geodesics falling into such a singularity go to a vanishing value. In
other words, all physical objects are crushed to a zero volume as they fall
into a strong curvature singularity.
the loop quantum gravity formalisms. All the same, it would appear that the
goal of achieving a fully consistent and complete quantum theory of gravity
is yet to be reached. The possibility of singularity resolution through quan-
tum gravity is clearly linked to such efforts, and much work may be needed
before it can be fully realized.
The singularity theorems predicting the occurrence of the spacetime singu-
larities, however, contain three main assumptions under which the existence
of a singularity is predicted in the form of geodesic incompleteness in the
spacetime. These are in the form of a typical causality condition that ensures
a suitable and physically reasonable global structure of the spacetime, an
energy condition that requires the positivity of energy density at the classi-
cal level as seen by a local observer, and finally a condition demanding that
trapped surfaces must exist in the dynamical evolution of the universe, or
in the later stages of a continual gravitational collapse. A trapped region in
the spacetime consists of trapped surfaces, which are two-surfaces, such that
both in going as well as outgoing wavefronts normal to it must converge.
Such trapped surfaces then necessarily give rise to a spacetime singularity
either in a gravitational collapse or in cosmology.
For the same reason, the process of trapped surface formation in a gravita-
tional collapse is also central to blackhole physics. The role of such a trapping
of light and matter within the framework of Einstein’s theory of gravitation
was highlighted by Datt (1938) and Oppenheimer and Snyder (1939), within
the context of the continual collapse of a massive matter cloud. They studied
the collapse of pressureless dust clouds using general relativity, and showed
that it leads to the formation of an event horizon, and a blackhole as the
collapse endstate, assuming that the spatial density distribution within the
star was strictly homogeneous, that is, ρ = ρ(t) only. In the later stages
of the collapse, an apparent horizon develops that is the boundary of the
trapped region, thus giving rise to an event horizon. Once the collapsing star
has entered the event horizon, the causal structure of the spacetime there
would then imply that no non-spacelike curves from that region would escape
away, and the star is cut off from the faraway observers, thus giving rise to
a blackhole in the spacetime.
While blackhole physics has led to several interesting theoretical, as well
as observational, developments and discussions, it is necessary, however, to
study more realistic models of the gravitational collapse in order to put
blackhole physics on a sound footing. This is because the Oppenheimer–
Snyder–Datt scenario is rather idealized and pressures would play an
important role in the dynamics of any realistic collapsing star. Also, density
distributions within the star could not be completely homogeneous in any
physically realistic model, but would be higher at the center of the cloud,
with a typical negative gradient as one moved away from the center. The
5.4 Blackhole paradoxes 237
could be delayed or avoided during the collapse. This could then offer a some-
what natural resolution to the issues such as those above. The main issue
would then be that of finding a mechanism, either at a classical or quantum
level, so that much of the matter of the collapsing star can escape away and
would be thrown out during the final stages of the gravitational collapse,
thus resolving the problems such as those above, and also that of the infinite
density spacetime singularity. Towards realizing such a scenario, the usual
physical reasonability conditions could be imposed, and while regularity of
the initial data could be required, as well as the weak energy condition, the
pressure could be allowed to be negative.
Such a possibility at a quantum level is discussed in some detail in the
next section, where the role of quantum effects towards generating an out-
wards flow of matter and radiation from the collapsing cloud in the very
late stages of the collapse is explored. In any case, the occurrence of a
spacetime singularity indicates the breakdown of the Einstein gravitation
theory in these extreme regimes. The singularity problem and various black-
hole paradoxes at classical and quantum levels were mentioned above. Note
that quantum corrections could generate a strong negative pressure in the
interior of the cloud in the very late stages of the collapse, where the
classical theory should break down. The collapsing star could then radi-
ate away most of its matter as the process of the gravitational collapse
evolves, so as to avoid the formation of trapped surfaces and the spacetime
singularity.
where dΩ2 is the line element on a two-sphere. In this comoving frame, the
energy-momentum tensor of the scalar field is given as
1 2
T t
t = −ρ(t) = − Φ̇ + V (Φ) , (5.66)
2
1 2
T rr θ φ
= T θ = T φ = p(t) = Φ̇ − V (Φ) , (5.67)
2
F −Ḟ
ρ= p= , (5.68)
R2 R R2 Ṙ
F
Ṙ2 = . (5.69)
R
d 3
a Φ = −a3 V (Φ),Φ . (5.70)
dt
Since the aim here is to construct a continual collapse model, the class with
ȧ < 0, which is the collapse condition implying that the area radius of a
shell at a constant value of comoving radius r decreases monotonically, is
considered. In general, there may be classes of solutions where a scalar field
may disperse also (see for example, Choptuik, 1993). The objective, however,
is to examine whether the singularities forming in the scalar field collapse
could be naked, or would be necessarily covered within a blackhole, and if so
under what conditions. The singularity resulting from a continual collapse is
given by a = 0, that is, when the scale factor vanishes and the area radius
for all the collapsing shells becomes zero. At the singularity ρ → ∞.
The key factor that decides the visibility, or otherwise, of the singularity is
the geometry of the trapped surfaces that may form as the collapse evolves,
that is the two-surfaces in the spacetime from which both outgoing and in
going wavefronts necessarily converge. The boundary of the trapped region
in a spherically symmetric spacetime is given by the equation F = R, which
describes the apparent horizon for the spacetime. The spacetime region where
the mass function F satisfies F < R is not trapped, while F > R describes
a trapped region.
In terms of the scale variable a, the mass function can be written as
" #
3 1 3 1 2 2
F =r a Φ(a) ,a ȧ + V (Φ(a)) . (5.71)
3 2
242 Final fate of a massive star
This is because from the above Einstein equation for density, the mass
function can be solved as
1
F = ρ(t)R3 . (5.72)
3
From (5.72), it can be seen that
F 1
= ρ(t)r2 a2 . (5.73)
R 3
The above relation decides the trapping, or otherwise, in the spacetime as
the collapse develops. The classes of collapse solutions for a scalar field with
potential, and the trapping, or otherwise, as the collapse develops will be
constructed and investigated. This is relevant from the perspective of the
cosmic censorship conjecture in order to understand the development of a
blackhole or a naked singularity as the collapse outcomes.
Towards such a purpose, consider the class of models where, near the
singularity, the divergence of the density is given by ρ(t) ≈ 1/a(t). Then,
using (5.71), the above condition implies, near the singularity,
1 1
Φ(a)2,a ȧ2 + V (Φ(a)) = . (5.74)
2 a
Now, solving the equation of motion (5.69),
√
a
ȧ = − √ . (5.75)
3
The negative sign implies a collapse scenario where ȧ < 0. Using (5.75) and
(5.74) in the latter part of (5.68), Φ can be solved as
Φ(a) = − ln a. (5.76)
5
V (Φ) = eΦ . (5.77)
6
Therefore near the singularity,
1 2
ρ(t) ≈ , p(t) ≈ − . (5.78)
a(t) 3a(t)
In the model above the weak energy condition is satisfied as ρ > 0 and
ρ + p > 0, although the pressure would be negative. Also, it should be noted
that the pressure does not have to be negative from the initial epoch, because
the specific behavior of ρ in (5.78) has been required only near the singularity.
It is possible to choose a V (Φ) such that at the initial epoch 1/2Φ̇2 > V (Φ),
and then the pressure would be positive. But, near the singularity, V (Φ)
would behave according to (5.78), and hence the pressure should decrease
monotonically from the initial epoch and tend to −∞ at the singularity.
If, from an epoch t = t∗ (or equivalently for some a = a∗ ) the density
starts growing as a−1 , then integrating (5.75), the singular epoch can be
obtained as √
ts = t∗ + 2 3a∗ . (5.79)
Therefore, the collapse reaches the singularity in a finite comoving time,
where the matter energy density as well as the Kretschman scalar κ =
Rijkl Rijkl diverge. Note that from the equation of motion (3.31) it follows
that the metric function a is given by
2
√ 1 ∗
a(t) = a∗ − √ (t − t ) . (5.80)
2 3
This completes the interior solution within the collapsing cloud, thus giving
the required construction.
It can be seen from the above considerations that the absence, or otherwise,
of the trapped surfaces and the behavior of the pressure crucially depend on
the rate of divergence of the density ρ near the singularity. To examine this
more carefully, near the singularity set ρ = a−n with n > 0, as it is known
that ρ(t) must diverge as a(t) goes to zero in the limit of approach to the
singularity. In this case, solving the Einstein equations gives
(n − 3) −n
p= a . (5.81)
3
The corresponding values of Φ and V (Φ) are
√ n √nΦ
Φ = − n ln(a), V (Φ) = 1 − e . (5.82)
6
Again, calculating F/R in this general case,
F 1
= a2−n . (5.83)
R 3
Therefore, it can be seen that for low enough divergences (0 < n < 2) no
trapped surfaces form and there are negative pressures near the singularity.
For 2 ≤ n < 3, trapped surfaces do form, but the pressure still remains
244 Final fate of a massive star
where v is the retarded null coordinate and rv is the Vaidya radius. Matching
the area radius at the boundary,
rb a(t) = rv (v ). (5.86)
Then on the hypersurface Σ, the interior and exterior metrics are given by
and
2M (rv , v ) drv
ds2Σ+ =− 1− +2 dv 2 + rv2 dΩ2 . (5.88)
rv dv
Matching the first fundamental form on this hypersurface,
dv 1
= , (rv )Σ = rb a(t). (5.89)
dt Σ 1− 2M (rv ,v)
+ 2 dr v
rv dv
To match the second fundamental form (extrinsic curvature) for interior and
exterior metrics, note that the normal to the hypersurface Σ, as calculated
from the interior metric, is given by ni− = [0, a(t)−1 , 0, 0], and the non-
vanishing components of the normal derived from the generalized Vaidya
metric are
−1/2
2M (rv , v ) drv
n+ = − 1 −
v
+2 , (5.90)
rv dv
2M (rv , v ) drv
1− +
rv dv
nr+v = . (5.91)
2M (rv , v ) drv
1− +2
rv dv
Here the extrinsic curvature is defined as
1
Kab = Ln gab , (5.92)
2
that is, the second fundamental form is the Lie derivative of the metric with
respect to the normal vector n. The above expression is equivalent to
1
Kab = [gab,c nc + gcb nc ,a + gac nc ,b ]. (5.93)
2
−
Setting Kθθ − Kθθ
+
Σ
= 0 on the hypersurface Σ,
2M (rv , v ) drv
1− +
rv dv
rb a(t) = rv . (5.94)
2M (rv , v ) drv
1− +2
rv dv
Simplifying the above using (5.89) and (5.69) on the boundary, one gets
The above equations, together with (5.94), completely specify the matching
at the boundary of the collapsing scalar field.
It is known (see for example, Wang and Wu, 1999), that a generalized
Vaidya spacetime describes the matter that is a combination of matter fields
of Type I and Type II. Specifically, the energy–momentum tensor Tik for
(n)
the matter can be written as a linear superposition of two tensors T ik and
(m)
T ik given by
(n)
T ik = µli lk , (5.98)
(m)
T ik = (ρ + P )(li nk + ni lk ) + P gik . (5.99)
M = M (rv ), which are two different solutions of (3.102) (see for exam-
ple, Joshi and Dwivedi, 1999; Wang and Wu, 1999; Giambo, 2005). These
give two unique exterior spacetimes, both of which are subclasses of the
generalized Vaidya geometry described above.
Now it can be seen that at the singular epoch t = ts , 2M (rv , v )/rv →
0. Therefore the exterior metric around the singularity smoothly trans-
forms to
with the classical energy density and pressure of the scalar field of
Then it is easily seen that F/R = (κ/3)ln−4 a2−n r2 . Therefore, in the collaps-
ing phase as a → 0, whether or not the trapped surfaces form is determined
by the value of n. As discussed above, for 0 < n < 2, if no trapped surfaces
exist initially then no trapped surfaces would form until the epoch a(t) = 0
√ 2/n
with a(t) = 1 − n t/2 3 .
The absence of trapped surfaces is accompanied by a negative pressure,
implying that for a constant value of the comoving coordinate r, Ḟ is nega-
tive and so the mass contained in the cloud of that radius keeps decreasing.
This leads to a classical outward energy flux. As the collapse proceeds, the
scale factor vanishes in finite time and the physical densities blow up, leading
to a naked singularity. Since no trapped surfaces form during the collapse,
the outward energy flux will, in principle, be observable. However, near the
singularity, when the energy density is close to Planckian values, this clas-
sical picture has to be modified and the scenario incorporating quantum
gravity modifications into classical dynamics needs to be investigated. A
non-perturbative semi-classical modification, based on loop quantum grav-
ity for the interior, was considered by Goswami, Joshi, and Singh (2006),
250 Final fate of a massive star
Singularity
Classical epoch
t =0
which is discussed below in some detail. The underlying geometry for the
FRW spacetime in loop quantum cosmology is discrete and both the scale
factor and the inverse scale factor
operators have discrete eigenvalues. In
particular, a critical scale a∗ = jγ/3P exists below which the eigenval-
ues of the inverse scale factor become proportional to the positive powers
of the scale factor (see Fig. 5.2). Here γ ≈ 0.2375 is the Barbero–Immirzi
parameter, P is the Planck length, and j is a half-integer free parameter
that arises because the inverse scale factor operator is computed by tracing
over SU(2) holonomies in an irreducible spin j representation. The value
of this parameter is arbitrary and is constrained only by phenomenological
considerations.
The change in the behavior of the classical geometrical density (1/a3 ) for
scales a ∼ a∗ can be approximated by
dj (a) = D(q ) a−3 , q := a2 /a2∗ , a∗ := jγ/3 P , (5.106)
with
$
6
D(q ) = (8/77) q 3/2 7 (q + 1)11/4 − |q − 1|11/4
%6
− 11q (q + 1)7/4 − sgn (q − 1)|q − 1|7/4 . (5.107)
and
2 1 Ḋ(q ) Φ̇2 Ḋ(q )
peff = 1− − D(q ) V (Φ) − V (Φ) . (5.112)
3 (ȧ/a) D(q ) 2 3(ȧ/a)
F = (κ/3)(d−1 2 3 3
j Φ̇ /2 + a V (Φ)) r . (5.113)
252 Final fate of a massive star
For a < a∗ , as the scale factor decreases, the energy density and mass in the
interior decrease and the negative pressure increases strongly. This leads to a
strong burst of matter. The absence of trapped surfaces enables the quantum
gravity induced burst to propagate via the generalized Vaidya exterior to an
observer at infinity.
5.5 Resolution of a naked singularity 253
In the toy model considered above, it was shown that the classical outcome
and evolution of the collapse is radically altered by the non-perturbative
modifications to the dynamics, as induced by the quantum effects. These
considerations are, of course, within the mini-superspace setting, and the
general case of inhomogeneities and anisotropies remains open. However, the
possibility of such observable signatures in astrophysical bursts originating
from the quantum gravity regime near the spacetime singularity is indeed
intriguing, indicating that gravitational collapse scenarios can be used as
probes to test quantum gravity models.
References
255
256 References
Hagerdorn, R. (1968). ‘Hadronic matter near the boiling point’. Nuovo Cimento
A., 56, 1027.
Harada, T. (1998). ‘Final fate of the spherically symmetric collapse of a perfect
fluid’. Phys. Rev. D., 58, 104015.
Harada, T. and Maeda, H. (2001). ‘Convergence to a self-similar solution in
general relativistic gravitational collapse’. Phys. Rev. D., 63, 084022.
Harada, T., Iguchi, H. and Nakao, K. I. (2000). ‘Naked singularity explosion’.
Phys. Rev. D., 61, 101502.
Harada, T., Iguchi, H. and Nakao, K. (2002). ‘Physical processes in naked
singularity formation’. Prog. Theor. Phys., 107, 449.
Harrison, B. K., Thorne, K. S., Wakano, M. and Wheeler, J. (1965). Gravitation
Theory and Gravitational Collapse. Chicago: University of Chicago Press.
Hartle, J. (1978). ‘Relativistic stars, gravitational collapse, and black holes’. In
Relativity, Astrophysics and Cosmology, ed. W. Israel. Dordrecht: Reidel.
Hawking, S. W. (1971). ‘Gravitational radiation from colliding black holes’. Phys.
Rev. Lett., 26, 1344.
Hawking, S. W. (1975). ‘Particle creation by black holes’. Commun. Math. Phys.,
43, 199.
Hawking, S. W. and Ellis, G. F. R. (1973). The Large Scale Structure of Space-
time. Cambridge: Cambridge University Press.
Hawking, S. W. and Israel, W. (1979a). General Relativity: An Einstein
Centenary Survey. Cambridge: Cambridge University Press.
Hawking, S. W. and Israel, W. (1979b). ‘An introductory survey’. In General
Relativity: An Einstein Centenary Survey, eds S. W. Hawking and W. Israel.
Cambridge: Cambridge University Press.
Hawking, S. W. and Penrose, R. (1970). ‘The singularities of gravitational
collapse and cosmology’. Proc. Roy. Soc. Lond. A., 314, 529.
Hellaby, J. and Lake, K. (1985). ‘Shell crossings and the Tolman model’.
Astrophys. J., 290, 381.
Herrera, L. and Santos, N. O. (2005). ‘Cylindrical collapse and gravitational
waves’. Class. Quant. Grav., 22, 2407.
Hertog, T., Horowitz, G. and Maeda, K. (2004). ‘Generic cosmic-censorship
violation in anti-de Sitter space’. Phys. Rev. Lett., 92, 131101.
Hicks, N. J. (1965). Notes on Differential Geometry. Princeton, NJ: Van
Nostrand.
Hiscock, W. A., Williams, L. G. and Eardley, D. M. (1982). ‘Creation of particles
by shell-focusing singularities’. Phys. Rev. D., 26, 751.
Hubeny, V., Lio, X., Rangamani, M. and Shenker, S. (2004). ‘Comments on
cosmic censorship in AdS/CFT’, hep-th/0403198.
References 261
Yodzis, P., Seifert, H. -J. and Muller zum Hagen, H. (1974). ‘On the occur-
rence of naked singularities in general relativity II’. Commun. Math. Phys.,
37, 29.
Zehavi, I. and Dekel, A. (1999). ‘Constraints on the cosmological constant from
flows and supernovae’. Nature, 401, 252.
Index
269
270 Index