Dissipative Spherical Gravitational Collapse of Isotropic Fluid

We present a number of parametric class of exact solutions of a radiating star and the matching conditions required for the description of physically meaningful fluid. A number of previously known class of solutions have been rediscovered which describe well behaved nature of fluid distributions. The interior matter fluid is shear-free spherically symmetric isotropic and undergoing radial heat flow. The interior metric obeyed all the relevant physical and thermodynamic conditions and matched with Vaidya exterior metric over the boundary. Initially the interior solutions represent a static configuration of perfect 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 collapse begins and at the stage when collapsing configuration reaches the horizon of the black hole.


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.

Field Equations and Junction Conditions
Space-time of a radiating star is divided by its boundary into two significant regions, the interior space-time M − and the exterior space-time M + 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 ( ) ( ) 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 The metric in the interior of a shear-free spherically symmetric fluid distribution is given by The energy-momentum tensor for the matter distribution with radial heat flow is ( ) where  is the energy density of the fluid, p the isotropic pressure, w µ is the four-velocity and q µ the radial heat flow vector.Assuming comoving coordinates, we have The heat flow vector q µ is orthogonal to the velocity vector so that 0 q w µ µ = and takes the form 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 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 the coupling constant in geometrized units is 8π κ = ( ) . . 1 i e G c = = .The exterior space-time is described by Vaidya's metric [25] which represents an outgoing radial flow of radiation where v is the retarded time and

( )
M v 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]) ( ) ( ) ( ) ( ) where Σ m is the mass function calculated in the interior at Σ r r = (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 The total luminosity for an observer at rest at infinity is ( )

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 func-tions of r and t coordinates as In view of ( 19) and ( 20) the field Equations ( 8)-( 11) lead to the following system of equations where The isotropy of pressure would give the equation In the absence of dissipative force the Equation ( 14), ( ) ( ) , reduces to the condition ( ) where If we assume ( ) ( ) [20]), solution of ( 26) is 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 ( ) 1 g t = .
For collapsing configurations we must have ( ) 20) and (30) represents a static perfect fluid at t → −∞ and then the fluid gradually starts evolving into a non-adiabatic radiating collapse.Equations ( 29) and ( 30) become ( ) ( ) We observed that the function

( )
f t decreases monotonically from the value ( ) The new parametric class of solutions of Equation ( 26) obtained by Tewari [20] is given as ( ) where , , , , n l C C D and 2 D are constants and ( ) ( ) . For different values of n or l Equations ( 33) and (34) give a variety of solutions (in Table 1).For 0,1 n = ; 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], , horizon-free case studied by Tewari and Charan [21].We here present some more special solutions and a detailed study of a class of solutions.
) for a

Collapsing Radiating Star
In order to construct the new realistic model we assume 5 3 n = − , and from ( 33) and (34) we obtain ( ) ( ) ( ) In view of ( 36) and ( 37) we obtain from ( 24) and ( 25) The junction condition ( ) The central values of 0  and 0 p are given by ( ) 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.
The solution should have monotonically decreasing expressions for the pressure and density with respect to r .For this, using (38) and (39) respectively we get ( ) Thus the extrema of 0  occurs at the centre if 0 0 0 r ′ = ⇒ =  , ( ) holds for all the 2400 0, values of an 52 d 9 Thus the extrema of 0 p occurs at the centre if 0 0 0 p r ′ = ⇒ = , ( ) ( ) The above observed range of 2 D 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 By using ( 7), ( 21) -( 23), ( 31), ( 32), ( 36) and (37) the explicit expressions for  , p , q and Θ become ( ) ( ) ) We can see the physical parameters , , p q  are finite, positive, monotonically decreasing at any instant with respect to radial coordinate for Σ 0 r r ≤ ≤ .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 1 1 , C D and 2 C .
The total energy entrapped inside Σ is given by (15), which becomes, using (19), (20), (31), (32), (36) and (37) where ( ) ( ) 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 ( ) The Equations ( 55) and (56) show that L ∞ vanishes in the beginning when ( ) 1 f t → and at the stage when Σ z → ∞ .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 ( ) ( ) and using (32) in view of (57) we have

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]) ; , g w w w q q g w w T Tw where  (≥ 0) is the thermal conductivity and τ (≥ 0) is the relaxation time.To get a simple estimate of the temperature evolution, by setting 0 τ = in (59) we get If we assume thermal conductivity , where γ and Ω are positive constants, on integration (60) yields ( ) ( ) ( ) where ( ) T t is an arbitrary function of t.The effective surface temperature observed by external observer can be calculated from the expression The temperature inside the star is given by

Conclusion
We have proposed some special solutions of Tewari [20] and a new class of exact solutions corresponding to 5 3 n = − has been studied in detail.The interior metric is in separable form and the model seems to be physically 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.

Table 1 .
Different values of n or l.

Detailed Study of a New Class of Exact Solutions Corresponding to n
− + 4.