Dissipative Spherical Gravitational Collapse of Isotropic Fluid ()
1. Introduction
Modern astrophysics opens numerous challenges for researchers to study various aspects of gravitational collapse. It is one of the most interesting phenomena in relativistic astrophysics. When a body does not have substantially strong pressure gradient force, it may continue collapsing because of its own gravity. Therefore, a detailed description of gravitational collapse of massive stars and the modelling of the structure of compact objects under various conditions is the problem that attracts significant attention of researchers in relativistic astrophysics. With reference to cosmic censorship conjecture gravitational collapse must terminate into a black hole (Penrose [1] ), but there are several counter examples where a naked singularity has more chances to be formed (Joshi and Malafarnia [2] and references therein). Still there is no established theory which can determine whether there will be the formation of a black hole or a naked singularity.
To understand the nature of collapse, it is necessary to form a realistic model, which requires solving non linear differential equations. Solving such equations is a very difficult task; various efforts have been made in this direction. The pioneering work in this area was started by Oppenheimer and Snyder [3] , in which they assumed a spherically symmetric distribution of matter, adiabatic flow and the equation of state in the form of dust with Schwarzschild exterior. Later on, taking into account the outgoing radiation from collapsing spherical fluid, Vaidya [4] initiated the problem and the modified equations proposed by Misner [5] and Lindquist et al. [6] for an adiabatic distribution of matter.
It is an established fact that gravitational collapse is a highly dissipating energy process (Herrera and Santos [7] , Herrera et al. [8] , Mitra [9] and references therein) which plays a dominant role in the formation and evolution of stars. In diffusion approximation the dissipative model is described by the heat flow type vector. Santos [10] studied the junction conditions of collapsing spherically symmetric shear-free non-adiabatic fluid with radial heat flow which was based on relativistic models suggested by Glass [11] . On the similar ground a number of stellar models (Maiti [12] , de Oliveira et al. [13] , Bonnor et al. [14] , Banerjee et al. [15] , Herrera et al. [16] , Tewari [17] - [20] , Tewari and Charan [21] , Sharma and Tikekar [22] , Ivanov [23] , Pinheiro and Chan [24] and also references therein) have been reported with the impact of various factors such as shear, inhomogeneity, anisotropy, electromagnetic field and various dissipative processes on the evolution.
Shear-free fluid distributions with radial heat flow are often studied in order to simplify the calculations and allow realistic analytic solutions. In view of the above arguments we present some special solutions of Tewari [20] and demonstrate a detailed study of one such solution in order to construct a realistic model of collapsing radiating star. The interior space-time metric is matched with Vaidya exterior metric (Vaidya [25] ) over the boundary, and the final fate of our model is formation of a black hole. The paper is organised as follows. In Section 2 the field equations and the junction conditions which match the interior metric of the collapsing fluid with the exterior metric are given. In Section 3 a new class of exact solutions of the field equations and a table of some special solutions are presented. In Section 4 a detailed study of a class of solutions for a collapsing radiating star is given. Section 5 describes temperature profile of the solution and finally in Section 6 some concluding remarks have been made.
2. Field Equations and Junction Conditions
Space-time of a radiating star is divided by its boundary into two significant regions, the interior space-time
and the exterior space-time
for a stellar model. Each of these specific regions is described by a peculiar smooth time-like three dimensional hyper surface containing
as its boundary. When approaching
from the exterior or the interior space time, we require
(1)
with this the line elements match on the boundary
.
Spherically symmetric collapsing distribution of dissipative fluid in the form of heat flow with a boundary of a time like spherical surface
with metric is considered as
(2)
The metric in the interior of a shear-free spherically symmetric fluid distribution is given by
(3)
The energy-momentum tensor for the matter distribution with radial heat flow is
(4)
where
is the energy density of the fluid,
the isotropic pressure,
is the four-velocity and
the radial heat flow vector. Assuming comoving coordinates, we have
(5)
The heat flow vector
is orthogonal to the velocity vector so that
and takes the form
(6)
The line element (3) corresponds to shear-free spherically symmetric fluid (Glass [26] ) as the shear tensor vanishes identically. The fluid collapse rate
of the fluid distribution (3) is given as
(7)
here and hereafter the dots and primes stand respectively for differentiation with respect to t and r.
Non-trivial Einstein’s field equations in view of (3) and (4) are given by following system of equations
(8)
(9)
(10)
(11)
the coupling constant in geometrized units is
.
The exterior space-time is described by Vaidya’s metric [25] which represents an outgoing radial flow of radiation
(12)
where
is the retarded time and
is the exterior Vaidya mass.
The junction conditions for matching two line elements (3) and (12) continuously across a spherically symmetric time-like hyper surface
are well known and obtained by (Santos [10]
(13)
(14)
(15)
where
is the mass function calculated in the interior at
(Cahill et al. [27] and Misner and Sharp [28] ).
Some other characteristics of the model such as the surface luminosity and the boundary redshift observed on
are
(16)
(17)
The total luminosity for an observer at rest at infinity is
(18)
3. Solution of the Field Equations
In order to solve the field equations we choose a particular form of the metric coefficients given in (3) into functions of
and
coordinates as
(19)
(20)
In view of (19) and (20) the field Equations (8)-(11) lead to the following system of equations
(21)
(22)
(23)
where
(24)
(25)
The isotropy of pressure would give the equation
(26)
In the absence of dissipative force the Equation (14),
, reduces to the condition
and yields at ![]()
(27)
where
(28)
If we assume
(Tewari [20] ), solution of (26) is
(29)
(30)
where β is an arbitrary constant and the constant of integration in (30) has been eliminated by the means of transformation in t. The solution (29), (30) is identical to the solution presented by de Oliveira et al. [13] , Bonnor et al. [14] with
. For collapsing configurations we must have
. From Equation (29) we have
as
is positive. We choose
in order to have
as f → 1 the solution (19), (20) and (30) represents a static perfect fluid at
and then the fluid gradually starts evolving into a non-adiabatic radiating collapse.
Equations (29) and (30) become
(31)
(32)
We observed that the function
decreases monotonically from the value
at
to
at
.
The new parametric class of solutions of Equation (26) obtained by Tewari [20] is given as
(33)
(34)
where
and
are constants and
(35)
n is real if
or
.
For different values of n or l Equations (33) and (34) give a variety of solutions (in Table 1). For
; we rediscover the Schwarzschild interior solution and the collapsing radiating star model in this case has been studied by de Oliveira et al. [13] and Bonnor et al. [14] and for n = −1, the solution reduces to Banerjee et al. [15] ,
for
, it reduces to Tewari [19] [20] , for
, horizon-free case studied by Tewari and Cha-
ran [21] . We here present some more special solutions and a detailed study of a class of solutions.
![]()
Table 1. Different values of n or l.
4. Detailed Study of a New Class of Exact Solutions Corresponding to
for a
Collapsing Radiating Star
In order to construct the new realistic model we assume
, and from (33) and (34) we obtain
(36)
(37)
In view of (36) and (37) we obtain from (24) and (25)
(38)
(39)
The junction condition
gives
(40)
The central values of
and
are given by
(41)
(42)
A physically reasonable solution should satisfy certain conditions. They are
(i) The central values of pressure, density and metric potential component should be non-zero positive definite.
This condition is satisfied if
and subjecting to the condition
yields D2 > 0, which holds good if
.
(ii) The solution should have monotonically decreasing expressions for the pressure and density with respect to
. For this, using (38) and (39) respectively we get
(43)
Thus the extrema of
occurs at the centre if
,
(44)
![]()
(45)
Thus the extrema of
occurs at the centre if
,
(46)
The above observed range of
is covered under energy condition (i).
Thus the expressions of right hand side of (44) and (46) are hold good for all the already mentioned ranges of parameters satisfying physical viability of the model, showing thereby that the density and pressure are maximum at the centre and monotonically decreasing towards the boundary surface.
Further, it is mentioned here that the boundary of the collapsing radiating star is established only when
and.
By using (7), (21) - (23), (31), (32), (36) and (37) the explicit expressions for
,
,
and
become
(47)
(48)
(49)
(50)
where
(51)
We can see the physical parameters
are finite, positive, monotonically decreasing at any instant with respect to radial coordinate for
. Initially collapse is zero and it becomes infinite at the final phase of the configuration. The constant
is positive for the given range of the parameters
and
.
The total energy entrapped inside
is given by (15), which becomes, using (19), (20), (31), (32), (36) and (37)
(52)
where
(53)
Using (16) - (20), (31), (32), (36) and (37) the luminosity and the red shift observed on
and luminosity observed by a distant observer are given by
(54)
(55)
(56)
The Equations (55) and (56) show that
vanishes in the beginning when
and at the stage when
.
The black hole formation time when the collapse reaches the horizon of the black hole occurs when the surface redshift goes to infinity, for this the term in the parentheses in Equation (56) goes to zero and we obtain
(57)
and using (32) in view of (57) we have
(58)
5. Temperature Profile
Now we will investigate the evolution of temperature of the collapsing star. In view of extended irreversible thermodynamics, the relativistic Maxwell-Cattaneo relation for temperature governing the heat transport within the collapsing matter in the truncated Israel-Stewart theory has the form (Israel et al. [29] , Maartens [30] , and Martinez [31] )
(59)
where
(≥ 0) is the thermal conductivity and
(≥ 0) is the relaxation time. To get a simple estimate of the temperature evolution, by setting
in (59) we get
(60)
If we assume thermal conductivity
, where
and
are positive constants, on integration (60) yields
(61)
where
is an arbitrary function of t.
The effective surface temperature observed by external observer can be calculated from the expression (Schwarzschild [32] )
(62)
The effective surface temperature is zero in the beginning
and at the stage when collapsing configuration reaches the horizon of the black hole
.
where the constant
in Photon is given by
(63)
where
and
denoting respectively Boltzmann and Plank constants.
Choosing
which represents radiation interaction with matter through the diffusive approximation (Misner and Sharp [33] ). The arbitrary function
is determined by using (61) and (62) as
(64)
The temperature inside the star is given by
(65)
6. Conclusion
We have proposed some special solutions of Tewari [20] and a new class of exact solutions corresponding to
has been studied in detail. The interior metric is in separable form and the model seems to be physical-
ly and thermodynamically sound as it corresponds to well-behaved nature for fluid density, isotropic pressure and radiation flux density throughout the fluid sphere. Initially the interior solutions represent a static configuration of non-dissipative fluid which then gradually starts evolving into radiating collapse. The apparent luminosity as observed by the distant observer at rest at infinity and the effective surface temperature are zero in remote past at the instant when the collapse begins and at the stage when collapsing configuration reaches the horizon of the black hole.
7. Acknowledgements
We thank the anonymous referee for valuable suggestions.