IJAA  Vol.6 No.2 , June 2016
Spherical Gravitational Collapse of Anisotropic Radiating Fluid Sphere
Abstract: We here present a relativistic model for a spherically symmetric anisotropic fluid to study the various factors of physical and thermal phenomenon during the evolution of a collapsing star dissipating energy in the form of radial heat flow. We also proposed a table of some new parametric class of solutions which will be useful for constructing the new compact star models. The constructed algorithm obeys all the relevant requirements of a realistic model and matched with Vaidya exterior metric over the boundary. At the initial stage 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.

Received 22 March 2016; accepted 15 May 2016; published 18 May 2016

1. Introduction

In order to construct the new realistic models, it is desirable to solve the Einstein's field equations but due to non linear character of the equations it is a very difficult task; various efforts have been made in this direction. The maiden exact solution of spherical gravitational collapse was due to Oppenheimer and Snyder [4] and which apparently suggested that such a collapse results in formation of Black Holes. Taking into account the outgoing radiation from collapsing spherical fluid Vaidya [5] [6] initiated the problem describing the exterior field and Santos [7] presented the junction conditions of collapsing spherically symmetric shear-free non-adiabatic fluid with radial heat flow.

Herrera and Santos [8] and Mitra [9] established the fact that gravitational collapse is a high energy dissipating process which plays a dominant role in the formation and evolution of stars. Historically the dissi- pation of energy from collapsing fluid distribution is described in two limiting cases. The first case describes the free streaming approximation (pure radiation) while second one is diffusion approximation (heat flow). The prominent work in pure radiation is due to Tewari and he solved the Einstein’s field equations with a new approach and developed the Quasar models [10] - [13] . While a number of realistic models in diffusion approximation with the impact of various factors such as inhomogeneity, anisotropy, viscosity, electromagnetic field and various dissipative processes on the evolution are critically discussed by de Oliveira et al. [14] ; Bonnor et al. [15] ; Herrera et al. [16] ; Maharaj and Govender [17] ; Ivanov [18] ; Tewari [19] [20] ; Pinheiro and Chan [21] ; Tewari and Charan [22] [23] ; Sharif and Iftikhar [24] .

Bowers and Liang [25] did the pioneering work in the field of anisotropic fluid models which enabled researchers to study the effect of anisotropy on the physical behavior of a star undergoing gravitational collapse. Herrera and Santos [26] explored the properties of anisotropic self-gravitating spheres using the perturbation method. Herrera et al. [27] investigated that the local pressure anisotropy is one of the responsible factors for inhomogeneities in energy density and in this series some interesting noticeable work is reported by Ivanov [28] [29] ; Reddy et al. [30] ; Tewari and Charan [31] [32] ; Govender et al. [33] and a different approach with a new concept Horizon function is very recently introduced for anisotropic fluids by Ivanov [34] [35] .

The main objective of this work is to present a simple anisotropic collapsing radiating fluid model and discuss all the relevant thermal and physical conditions by taking Tewari and Charan [31] solution as seed solution and we propose some new exact solutions for this parametric class of solution which are useful in construction of the compact stars models. The interior metric is spherically symmetric, shear-free, anisotropic radiating away its energy in the form of radial heat and contracting in size during the process of collapse. The interior metric matched with the Vaidya exterior metric [6] over the boundary. Final fate of our model is the formation of a black hole. The paper is organised as: In Section 2 the space-times and the junction conditions for collapsing system are given. In Section 3 solution of Einstein’s field equations is presented. Section 4 describes a parametric class of solution and some new solutions. In Section 5 detailed study of the model in which we have obtained expressions of various physical parameters is described. Section 6 contains the temperature profile of the collapsing body. Finally in Section 7 some concluding remarks have been made.

2. The Space-Times and Junction Conditions

The interior space-time of a shear-free spherically symmetric fluid with the coordinates is given by


The energy-momentum tensor for the matter distribution with anisotropy in pressure is


where is the energy density of the fluid, the radial pressure, the tangential pressure, is the four- velocity, the radial heat flow vector and is a unit space like four vector along the radial direction.

Assuming comoving coordinates, we have. The heat flow vector is orthogonal to the velocity vector so that and takes the form.

Since the interior fluid is radiating energy in the form of heat therefore the exterior space-time of a collapsing radiating star is described by Vaidya’s outgoing exterior metric [6]


where v is the retarded time and is the Vaidya mass.

The junction conditions for radiating star matching two line elements (1) and (3) at the boundary continuously across a spherically symmetric time-like hyper surface are very well known and given by Santos [7]




where is the mass function calculated in the interior at (Cahill et al. [36] ; Misner and Sharp [37] ). The primes and dots stand for differentiation with respect to r and t respectively.

Some other physical and thermal features of the collapsing matter are the surface luminosity and the boundary redshift observed on are



The total luminosity for an observer at rest at infinity is


3. Solution of Einstein’s Field Equations

In order to solve the Non-trivial Einstein’s field equations which are generated by (1) and (2), we choose a separable form of the metric coefficients given in (1) into functions of r and t coordinates as and. The coupling constant in geometrized units is taken as and in view of (1) and (2) with the help of separable metric coefficients we get the following expressions for field equations









Here the quantities with the suffix 0 corresponds to the static star model with metric components.

In the absence of non-adiabatic dissipative forces the Equation (5), , reduces to the condition and yields at




Tewari [20] by assuming solved the Equation (17) and the solution so obtained is identical to the solution presented by de Oliveira et al. [14] and Bonnor et al. [15] with.



Here we observed that the function decreases monotonically from the value at to at. It interpret that the collapse begins in the remote past and gradually starts evolve into radiating collapse.

4. The New Parametric Class of Solution and Some Proposed Solutions

To find a new parametric class of exact solutions of pressure anisotropy equation which is created by the Equa- tions (11) and (12), Tewari and Charan [31] assumed that the anisotropy evolves and got the following time independent differential equation


Making an adhoc relationship between the variables in (21), the above mentioned authors obtained the following solution




where and are constants and


where n is real if or.

For different values of n or l Equations (22) and (23) give a variety of solutions and they are categorized as isotropic pressure and homogeneous density, isotropic pressure and inhomogeneous density while some inhomogeneous density and anisotropic pressure. For we get the homogeneous density and anisotropic pressure, for we get the isotropic pressure and homogeneous density and in this case solution reduces to Banerjee et al. [38] and for we get isotropic pressure and inhomogeneous density. In order to maintain the inhomogeniety and anisotropy for a collapsing radiating star, a Horizon-free case has been studied by Tewari and Charan [31] for and in an other study the same authors presented a Black hole model for [32] . Since a number of solutions can be obtained with the above mentioned parametric class of solution, so keeping this point in mind we here find some exact solutions and listed them in a table and they will be fruitful for further study to construct the various models of radiating and static stars.

5. Physical and Thermal Analysis of a Specific Model

To construct a new relativistic model, in the present study we assume, to maintain the anisotropy for a collapsing radiating star and using (22)-(24) we get







The junction condition gives


Here from (30) and (31), we are seeing that at the centre radial and transverse pressures are equal and anisotropy vanishes there.

A physically reasonable solution should satisfy certain energy conditions and they are:

The central values of both the pressures, energy density and gravitational potential should be non-zero positive definite and the solution should have monotonically decreasing expressions for the pressures and density with the increase of r. Thus in view of these conditions, we write the bounds of model parameters and.

Now using (10)-(13), (19), (26) and (27) the expressions for and q reduce the following





We can see the physical parameters are finite, positive, monotonically decreasing at any instant with respect to radial coordinate for.

The fluid collapse rate with the help of (1), (19), (26) and (27) is


where by using (18), (26) and (27) we have


The constant is positive for the given range of the parameters and.

The total energy entrapped inside the surface is given by using (6), (18), (26), (27) and (38)




Using (4) and (27), we get the physical radius of the collapsing radiating star as


Using (7)-(9), (18), (26), (27) and (38) the expressions for the surface luminosity, the boundary redshift on and the luminosity for distant observer at rest at infinity respectively are




To obtain the black hole formation time, the surface redshift goes to infinity, for this the term in the parentheses of (44) goes to zero and we have



6. Temperature Profile for Collapsing Radiating Star

To investigate the temperature inside and on the boundary surface, we utilize temperature gradient law (Israel et al. [39] ; Maartens [40] ; and Martinez [41] )


where is the thermal conductivity and is the relaxation time. To get a simple estimate of the temperature evolution, we set relaxation time as zero in above expression and obtained


The effective surface temperature observed by external observer can be calculated from the expression (Schwarzschild [42] )


where for photons the constant is given by, here k and denoting respectively Boltzmann and Plank constants. Choosing which represents radiation interaction with matter through the diffusive approximation (Misner and Sharp [43] .

The arbitrary function is determined by using (48) and (49) as


Temperature distribution throughout the interior of the collapsing radiating star is given by


It follows that the surface temperature of the collapsing star tends to zero at the beginning of the collapse and the stage of formation of black hole.

7. Conclusion

We here presented a new radiating fluid model collapsing in the influence of its own gravity using Tewari and Charan [31] solution as seed solution. The interior fluid is spherically symmetric shear-free anisotropic and radiating away its energy in the form of radial heat flow. We have obtained a class of exact solutions by assigning different values to the parameter n and proposed solutions are listed in the Table 1. Keeping in mind pressure anisotropy a simple radiating star model for studied in detail. Physical and thermal features of the model are significantly sound as it corresponds to well-behaved nature for the fluid density, both the radial

Table 1. Parametric class of solutions.

and tangential pressures and the flux density throughout the fluid sphere. We observed that the function decreases monotonically from the value at to at. The apparent luminosity

as observed by the distant observer at rest at infinity is zero in remote past at the instance when the collapse begins and at the stage of black hole formation.

Particularly we have constructed a radiating star model by taking suitable parameters, , ,. Initially when both the mass and physical radius of a core of pre Supernovae were and 15.8 km and they remain and 6.30 km respectively at the time of black hole formation.

The time of black hole formation is observed as 0.62996 S. We observed that the model is well behaved for the chosen constraints and there are a number of such sets for which the solution is well behaved. Thus the parametric class of solution and constructed table is very fruitful for further study of radiating and static compact star’s modeling. The surface temperature of the collapsing radiating star tends to zero at the beginning of the

collapse and the stage of formation of black hole.


Authors express their gratitude and thanks to the anonymous referee for his rigorous review and valuable suggestions.

Cite this paper: Tewari, B. , Charan, K. and Rani, J. (2016) Spherical Gravitational Collapse of Anisotropic Radiating Fluid Sphere. International Journal of Astronomy and Astrophysics, 6, 155-165. doi: 10.4236/ijaa.2016.62013.

[1]   Joshi, P.S. and Malafarina, D. (2011) Recent Developments in Gravitational Collapse and Spacetime Singularities. International Journal of Modern Physics D, 20, 2641-2729.

[2]   Penrose, R. (1969) Gravitational Collapse: The Role of General Relativity. La Rivista del Nuovo Cimento, 1, 252-276.

[3]   Virbhadra, K.S. (2009) Relativistic Images of Schwarzschild Black Hole Lensing. Physical Review D, 79, 083004.

[4]   Oppenheimer, J.R. and Snyder, H. (1939) Continued Gravitational Contraction. Physical Review, 56, 455.

[5]   Vaidya, P.C. (1951) The Gravitational Field of a Radiating Star. Proceedings of the National Academy of Sciences, India, Section A: Physical Sciences, 33, 264-276.

[6]   Vaidya, P.C. (1953) Newtonian Time in General Relativity. Nature, 171, 260-261.

[7]   Santos, N.O. (1985) Non-Adiabatic Radiating Collapse. Monthly Notices of the Royal Astronomical Society, 216, 403-410.

[8]   Herrera, L. and Santos, N.O. (2004) Dynamics of Dissipative Gravitational Collapse. Physical Review D, 70, 084004.

[9]   Mitra, A. (2006) Why Gravitational Contraction Must Be Accompanied by Emission of Radiation Both in Newtonian and Einstein Gravity. Physical Review D, 74, 024010.

[10]   Tewari, B.C. (1988) Radiating Fluid Spheres in General Relativity. Astrophysics and Space Science, 149, 233-239.

[11]   Tewari, B.C. (1994) Relativistic Radiating Fluid Distribution. Indian Journal of Pure and Applied Physics, 32, 504.

[12]   Tewari, B.C. (2006) Relativistic Model for Radiating Star. Astrophysics and Space Science, 306, 273-277.

[13]   Tewari, B.C. (2010) Radiating Fluid Balls in General Relativity. VDM Verlag, Saarbrucken.

[14]   de Oliveira, A.K.G., Santos, N.O. and Kolassis, C.A. (1985) Collapse of a Radiating Star. Monthly Notices of the Royal Astronomical Society, 216, 1001-1011.

[15]   Bonnor W. B., de Oliveira, A.K.G. and Santos, N.O. (1989) Radiating Spherical Collapse. Physical Review Letters, 181, 269-326.

[16]   Herrera, L., Di Prisco, A., Ospino, J., Fuenmayor, E. and Triconis, O. (2009) Structure and Evolution of Self-Gravi-tating Objects and the Orthogonal Splitting of the Riemann Tensor. Physical Review D, 79, Article ID: 064025.

[17]   Maharaj, S.D. and Govender, M. (2005) Radiating Collapse with Vanishing Weyl Stresses. International Journal of Modern Physics D, 14, 667-676.

[18]   Ivanov, B.V. (2012) Collapsing Shear-Free Perfect Fluid Spheres with Heat Flow. General Relativity and Gravitation, 44, 1835-1855.

[19]   Tewari, B.C. (2012) Relativistic Collapsing Radiating Stars. Astrophysics and Space Science, 342, 73-77.

[20]   Tewari, B.C. (2013) Collapsing Shear-Free Radiating Fluid Spheres. General Relativity and Gravitation, 45, 1547-1558.

[21]   Pinheiro, G. and Chan, R. (2013) Radiating Shear-Free Gravitational Collapse with Charge. General Relativity and Gravitation, 45, 243-261.

[22]   Tewari, B.C. and Charan, K. (2014) Radiating Star, Shear-Free Gravitational Collapse without Horizon. Astrophysics and Space Science, 351, 613-617.

[23]   Tewari, B.C. and Charan, K. (2015) Dissipative Spherical Gravitational Collapse of Isotropic Fluid. Journal of Modern Physics, 6, 453-462.

[24]   Sharif, M. and Iftikhar, S. (2015) Charged Dissipative Collapse of Shearing Viscous Star. Astrophysics and Space Science, 357, 79.

[25]   Bowers, R.L. and Liang, E.P.T. (1974) Anisotropic Spheres in General Relativity. Astrophysical Journal, 188, 657-665.

[26]   Herrera, L. and Santos, N.O. (1997) Local Anisotropy in Self-Gravitating Systems. Physics Reports, 286, 53-130.

[27]   Herrera, L., Di Prisco, A., Hernandez-Pastoraand, J.L. and Santos, N.O. (1998) On the Role of Density Inhomogeneity and Local Anisotropy in the Fate of Spherical Collapse. Physics Letters A, 237, 113-118.

[28]   Ivanov, B.V. (2010) Evolving Spheres of Shear-Free Anisotropic Fluid. International Journal of Modern Physics A, 25, 3975-3991.

[29]   Ivanov, B.V. (2011) Self-Gravitating Spheres of Anisotropic Fluid in Geodesic Flow. International Journal of Modern Physics D, 20, 319-334.

[30]   Reddy, K.P., Govender, M. and Maharaj, S.D. (2015) Impact of Anisotropic Stresses during Dissipative Gravitational Collapse. General Relativity and Gravitation, 47, 35.

[31]   Tewari, B.C. and Charan, K. (2015) Horizon Free Eternally Collapsing Anisotropic Radiating Star. Astrophysics and Space Science, 357, 107.

[32]   Tewari, B.C. and Charan, K. (2015) Gravitational Collapse, Shear-Free Anisotropic Radiating Star. arxiv:1503.02165

[33]   Govender, M., Bogadi, R.S., Lortan, D.B. and Maharaj, S.D. (2016) Radiating Collapse in the Presence of Anisotropic Stresses. International Journal of Modern Physics D, 25, Article ID: 1650037.

[34]   Ivanov, B.V. (2016) All Solutions for Geodesic Anisotropic Soherical Collapse with Shear and Heat Flow. Astrophysics and Space Science, 361, 18.

[35]   Ivanov, B.V. (2016) A Different Approach to Anisotropic Spherical Collapse with Shear and Heat Radiation. International Journal of Modern Physics D, 25, Article ID: 1650049.

[36]   Cahill, M.E. and McVittie, G.C. (1970) Spherical Symmetry and Mass Energy in General Relativity I. General Theory. Journal of Mathematical Physics, 11, 1382-1391.

[37]   Misner, C.W. and Sharp, D.H. (1964) Relativistic Equations for Adiabatic Spherically Symmetric Gravitational Collapse. Physical Review B, 136, 571-576.

[38]   Banerjee, A., Chatterjee, S. and Dadhich, N. (2002) Spherical Collapse with Heat Flow and without Horizon. Modern Physics Letters A, 17, 2335-2339.

[39]   Israel, W. and Stewart, J. (1979) Transient Relativistic Thermodynamics and Kinetic Theory. Annals of Physics, 118, 341-372.

[40]   Maartens, R. (1995) Dissipative Cosmology. Classical and Quantum Gravity, 12, 1455-1465.

[41]   Martinez, J. (1996) Transport Processes in the Gravitational Collapse of an Anisotropic Fluid. Physical Review D, 53, 6921-6940.

[42]   Schwarzschild, M. (1958) Structure and Evolution of stars. Dover, New York.

[43]   Misner, C.W. and Sharp, D.H. (1965) Spherical Gravitational Collapse with Energy Transport by Radiative Diffusion. Physics Letters, 15, 279-281.