Journal of GravityVolume 2013 (2013), Article ID 659605, 8 pageshttp://dx.doi.org/10.1155/2013/659605

Research Article

Collapse of a Relativistic Self-Gravitating Star with Radial Heat Flux: Impact of Anisotropic Stresses

Ranjan Sharma and Shyam Das

Department of Physics, P.D. Women’s College, Jalpaiguri, West Bengal 735101, India

Received 28 February 2013; Accepted 22 April 2013

Academic Editor: Sergei Odintsov

Copyright © 2013 Ranjan Sharma and Shyam Das. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We develop a simple model for a self-gravitating spherically symmetric relativistic star which begins to collapse from an initially static configuration by dissipating energy in the form of radial heat flow. We utilize the model to show how local anisotropy affects the collapse rate and thermal behavior of gravitationally evolving systems.

1. Introduction

In cosmology and astrophysics, there exist many outstanding issues relating to a dynamical system collapsing under the influence of its own gravity. In view of Cosmic Censorship Conjecture, the general relativistic prediction is that such a collapse must terminate into a space-time singularity covered under its event horizon though there are several counter examples where it has been shown that a naked singularity is more likely to be formed (see [1] and references therein). In astrophysics, the end stage of a massive collapsing star has long been very much speculative in nature [1, 2]. From classical gravity perspective, to get a proper understanding of the nature of collapse and physical behavior of a collapsing system, construction of a realistic model of the collapsing system is necessary. This, however, turns out to be a difficult task because of the highly nonlinear nature of the governing field equations. To reduce the complexity, various simplifying methods are often adopted and the pioneering work of Oppenheimer and Snyder [3] was a first step in this direction when collapse of a highly idealized spherically symmetric dust cloud was studied. Since then, various attempts have been made to develop realistic models of gravitationally collapsing systems to understand the nature and properties of collapsing objects. It got a tremendous impetus when Vaidya [4] presented a solution describing the exterior gravitational field of a stellar body with outgoing radiation and Santos [5] formulated the junction conditions joining the interior space time of the collapsing object to the Vaidya exterior metric [4]. These developments have enabled many investigators to construct realistic models of gravitationally evolving systems and also to analyze critically relevance of various factors such as shear, density inhomogeneity, local anisotropy, electromagnetic field, viscosity, and so forth, on the physical behaviour of collapsing bodies [6–52]. In the absence of any established theory governing gravitational collapse, such investigations have been found to be very useful to get a proper understanding about systems undergoing gravitational collapse.

The aim of the present work is to develop a simple model of a collapsing star and investigate the impact of pressure anisotropy on the overall behaviour of the collapsing body. Anisotropic stresses may occur in astrophysical objects for various reasons which include phase transition, density inhomogeneity, shear, and electromagnetic field [10, 53, 54]. In [53], it has been shown that influences of shear, electromagnetic field, and so forth on self-bound systems can be absorbed if the system is considered to be anisotropic, in general. Local anisotropy has been a well-motivated factor in the studies of astrophysical objects and its role on the gross features of static stellar configurations have been investigated by many authors (see, e.g., [10, 11, 55–59] and references therein). For dynamical systems, though pressure anisotropy is, in general, incorporated in the construction, very few works have been reported till date where impacts of anisotropy have been discussed explicitly [12–14, 26, 60]. Appropriate junction conditions for an anisotropic fluid collapsing on the background space time described by the Vaidya metric have been obtained in [27]. Considering a spheroidal geometry of Vaidya and Tikekar [61], Sarwe and Tikekar [50] have analyzed the impact of geometry vis-a-vis matter composition on the collapse of stellar bodies which begin to collapse from initial static configurations possessing equal compactness. Sharma and Tikekar [51, 52] have investigated the evolution of nonadiabatic collapse of a shear-free spherically symmetric object with anisotropic stresses on the background of space-time obtained by introducing an inhomogeneous perturbation in the Robertson-Walker space-time.

In the present work, we have developed a model describing a shear-free spherically symmetric fluid distribution radiating away its energy in the form of radial heat flux. The star begins its collapse from an initially static configuration whose energy-momentum tensor describing the material composition has been assumed to be anisotropic, in general. To develop the model of the initial static star, we have utilized the Finch and Skea [62] ansatz which has earlier been found to be useful to develop physically acceptable models capable of describing realistic stars in equilibrium [63–66]. The back ground space-time for a static configuration for the given ansatz has a clear geometrical interpretation as may be found in [65]. By assuming a particular form of the anisotropic parameter, we have solved the relevant field equations and constructed a model for the initial static stellar configuration which could either be isotropic or anisotropic in nature. The solution provided by Finch and Skea [62] is a subclass of the solution provided here. Since, the solution presented here provides a wider range of values of the anisotropic parameter, it enables us to examine the impact of anisotropic stresses on the evolution of a large class of initial static configurations.

Our paper has been organized as follows. In Section 2, we have presented the basic equations governing the system undergoing non-adiabatic radiative collapse. In Section 3, by assuming a particular anisotropic profile, we have solved the relevant field equations to develop a model for the initial static star. In Section 4, by stipulating the boundary conditions across the surface separating the stellar configuration from the Vaidya [4] space-time, we have solved the surface equation which governs the evolution of the initial static star that begins to collapse when the equilibrium is lost. In Section 5, we have analyzed the impact of anisotropy on gravitationally collapsing systems by considering evolution of initial static stars which could either be isotropic or anisotropic. Impact of anisotropy on the evolution of temperature has been analyzed in Section 6. Finally, some concluding remarks have been made in Section 7.

2. Equations Governing the Collapsing System

We write the line element describing the interior space-time of a spherically symmetric star collapsing under the influence of self-gravity (in standard coordinates , , , and ) as where, , , and are yet to be determined. Note that in (1), if we set , the metric corresponds to a static spherically symmetric configuration.

We assume that the matter distribution of the collapsing object is an imperfect fluid described by an energy-momentum tensor of the form where represents the energy density and , respectively, denote fluid pressures along the radial and transverse directions; is the -velocity of the fluid; is a unit space-like -vector along the radial direction; and is the heat flux vector which is orthogonal to the velocity vector so that and .

The Einstein’s field equations governing the evolution of the system are then obtained as (we set )

In (3)–(6), a “prime” denotes differentiation with respect to and a “dot” denotes differentiation with respect to . Making use of (4) and (5), we define the anisotropic parameter of the collapsing object as

We rewrite (3)–(5) as where , and denote the energy-density, radial pressure, and tangential pressure, respectively, of the initial static star.

3. Interior Space-Time of the Initially Static Configuration

In our construction, we assume that an initially static star (with in (1)), described by metric potentials , , and anisotropy , starts collapsing if, for some reasons, it loses its equilibrium. To develop a model of the initially static configuration, we first assume that the anisotropic parameter is separable in its variables so that . Equation (7) then reduces to which is independent of . Equation (11) can only be solved if any two of the unknown functions (, , and ) are specified. To develop a physically reasonable model of the initial static configuration, we first utilize the Finch and Skea [62] ansatz for the metric potential given by where is the curvature parameter describing the geometry of the configuration. In the static case, the ansatz (12) has a clear geometric characterization and has been found to yield realistic models for compact stellar objects [64–66]. Substituting (12) in (11), we obtain where we have used the following transformation:

To solve (13), we assume a particular anisotropic profile where is the measure of anisotropy. The motivation for choosing the particular form of the anisotropy parameter is the following: (1) it is physically reasonable as the anisotropy vanishes at the center (, i.e., ) as expected, and (2) it provides a solution of (13) in closed form. Note that corresponds to an initial static star which is isotropic in nature. Substituting (15) in (13), we get whose solution is obtained as where and are integration constants, is the Bessel function of first kind of order , and is the Bessel function of second kind of order . It is obvious that the solution is valid for . At , the Bessel function encounters a singularity and, therefore, we shall deal with the case separately.

Special Cases

Case 1 ( (). That is, initial static configuration is isotropic in nature). Equation (13) reduces to whose solution is found to be where, and are integration constants. This particular solution has been found previously in [62]. However, the solution (19) can be obtained directly from the solution (17) by substituting and setting and .

Case 2 ( ()). Equation (16) reduces to whose solution is given by where and are integration constants [64]. Unfortunately, the above solution can not be regained from the general solution (17) due to the properties of Bessel functions and should be treated separately.

We thus have a model for an initially static stellar configuration which could either be isotropic or anisotropic in nature. After loss of equilibrium, the initially static star starts collapsing, and to generate a solution of the subsequent evolving system, we need to determine . This will be taken up in the following section.

4. Exterior Space-Time and Junction Conditions

In our construction, evolution of the collapsing object is governed by the function which can be determined from the boundary conditions across the boundary surface joining the interior space-time and the exterior space-time described by the Vaidya [4] metric

In (22), denotes the retarded time and represents the total mass of the collapsing star. The junction conditions are determined by assuming a time-like -surface which separates the interior and the exterior manifolds [5]. Continuity of the metric space times () and the extrinsic curvatures () across the surface then yield the following matching conditions linking smoothly the interior () and the exterior () space-times across the boundary:

From (25) and (27), the mass enclosed within the boundary surface at any instant within a boundary surface can be written as

Combining (6), (9), and (26), together within the condition , we deduce the surface equation governing the collapsing matter in the form where we have defined

Note that for a given initial static configuration appears as a constant in (29). Following Bonnor et al. [6], we write (29) as a first order differential equation which admits a solution

Note that at , that is, at the onset of collapse, and as .

We, therefore, have a complete description of the interior and exterior space-times of the collapsing body. In the following section, we shall analyze the impact of anisotropy by making use of the solutions thus obtained.

5. Physical Analysis

To understand the nature of collapse, we assume that a star starts its collapse at (i.e., ) and the initial static star is described by the parameters , , and . We assume that the collapsing object has an initial radius and mass satisfying the condition . If for some reason, instability develops inside the star, it begins to collapse. The comoving boundary surface then starts shrinking until it reaches its Schwarzschild horizon value , where denotes the time of formation of the black hole corresponding to the value of .

From (28), the mass of the evolving star at any instant within the boundary radius may be written as where mass of the initial static star has the form

Consequently, the condition yields and the time of black hole formation is obtained as

The model parameters (namely, , , and ) involving the initial static star are determined using the following boundary conditions: where we have matched the static interior space-time to the Schwarzschild exterior and imposed the condition that the radial pressure () must vanish at the boundary.

Now, to get an insight about the effects of anisotropy, we consider different initial static configurations, both isotropic and anisotropic (). We consider different initially static stellar configurations of identical initial masses and radii (we have assumed and radius  km) with different values of the anisotropic parameter . Using numerical procedures, we have calculated the corresponding model parameters and also evaluated the time of formation of black holes and radius () and mass () of the black hole formed. Our results have been compiled in Table 1. Variations of for different choices of the anisotropic parameter have been shown in Figure 1. We have also calculated the rate of collapse in our model for different anisotropic parameters. The collapse rate, in our model, is obtained as

Table 1: Data for collapsing systems having different values of the anisotropic parameter . All the configurations start collapsing with initial masses and radii  km.

(km)

(km)
()

0.9
0.2848
−0.6341
20.8427
0.0244
0.2298
−2.3914
4.596
2.2980

0.6
0.3800
−0.4994
20.8427
0.0163
0.2298
−3.5816
4.596
2.2980

0.3
0.4065
−0.4619
20.8427
0.0142
0.2298
−4.1074
4.596
2.2980

0
0.42067
−0.4418
20.8427
0.0131
0.2298
−4.4578
4.596
2.2980

−0.3
0.4300
−0.4287
20.8427
0.0124
0.2298
−4.7239
4.596
2.2980

−0.6
0.4367
−0.4192
20.8427
0.0118
0.2298
−4.9401
4.596
2.2980

−0.9
0.4418
−0.4120
20.8427
0.0114
0.2298
−5.1230
4.596
2.2980

−1.5
0.4493
−0.4012
20.8427
0.0108
0.2298
−5.4231
4.596
2.2980

Figure 1: Variation of for different anisotropic parameters.

The collapse rate turns out to be for , respectively. Since, the collapsing object contracts in size as time evolves , in our construction, turns out to be negative. However, from the absolute values of , we note that for a positive anisotropy () the collapse rate increases as compared to an isotropic star while for a negative value of (), the rate decreases. From Table 1, we note the following.

Case  1. When (implying ), the horizon is formed at a faster rate as compared to , that is, isotropic case.

Case  2. When (implying ), the horizon is formed at a slower rate as compared to , that is, isotropic case.

Similar observations may be found in [12]. However, we note that mass and radius of the collapsed configuration do not depend on anisotropy in our formulation.

6. Thermal Behaviour

To analyze the impact of anisotropy on the evolution of temperature of the collapsing system, we use the relativistic Maxwell-Cattaneo relation for temperature governing the heat transport [60, 67, 68] given by where is the thermal conductivity and is the relaxation time. For the line element (1), (41) reduces to

Following an earlier work [51], we write the relativistic Fourier heat transport equation by setting in (42). For , combining (6) and (42), we get

Let us now assume that the thermal conductivity varies as , where and are constants. Equation (43) then yields where we have used (31). Integrating the above equation, we get

To get a simple estimate of temperature evolution, we set , , and and evaluate the surface temperature at any instant as

Time evolution of the surface temperature for different anisotropic parameter has been shown in Figure 2, where we have used the data from Table 1.

Figure 2: Evolution of surface temperature for different anisotropic parameters.

7. Discussions

We have developed a minimalistic and basic framework of a gravitationally collapsing system which has allowed us to examine the impact of pressure anisotropy explicitly. In our construction, we have ignored the effects of various other factors relevant to collapsing systems such as shear, viscosity, and charge. Moreover, the line element describing the interior space-time has been assumed to be spherically symmetric where the metric functions have been chosen to be separable in its variables. Though formulation of a more general framework is always preferred to describe a realistic situation, the simple model developed here provides a mechanism to capture the impact of anisotropy on gravitational collapse successfully. Though, as pointed out in [53], effects of factors like shear, charge, and so forth on self-bound systems can be absorbed by considering the system to be anisotropic, in general, we intend to formulate a more general framework so as to examine the combined impacts of various factor relevant to gravitationally collapsing systems. These issues will be taken up elsewhere.

Acknowledgment

R. Sharma is thankful to the Inter University Center for Astronomy and Astrophysics (IUCAA), Pune, India, where a part of this work was carried out under its Visiting Research Associateship Program.

References

P. S. Joshi and D. Malafarina, “Recent developments in gravitational collapse and spacetime singularities,” International Journal of Modern Physics D, vol. 20, Article ID 2641, 2011.

S. Chandrasekhar, “Stellar configurations with degenerate cores,” Observatory, vol. 57, pp. 373–377, 1934.

J. R. Oppenheimer and H. Snyder, “On continued gravitational contraction,” Physical Review, vol. 56, no. 5, pp. 455–459, 1939.

P. C. Vaidya, “The gravitational field of a radiating star,” Proceedings of the Indian Academy of Science A, vol. 33, no. 5, pp. 264–276, 1951.

N. O. Santos, “Non-adiabatic radiating collapse,” Monthly Notices of the Royal Astronomical Society, vol. 216, pp. 403–410, 1985.

W. B. Bonnor, A. K. G. de Oliveira, and N. O. Santos, “Radiating spherical collapse,” Physics Reports, vol. 181, no. 5, pp. 269–326, 1989.

A. K. G. de Oliveira and N. O. Santos, “Nonadiabatic gravitational collapse,” Astrophysical Journal, vol. 312, pp. 640–645, 1987.

A. K. G. de Oliveira, N. O. Santos, and C. A. Kolassis, “Collapse of a radiating star,” Monthly Notices of the Royal Astronomical Society, vol. 216, pp. 1001–1011, 1985.

A. K. G. de Oliveira, C. A. Kolassis, and N. O. Santos, “Collapse of a radiating star revisited,” Monthly Notices of the Royal Astronomical Society, vol. 231, pp. 1011–1018, 1988.

L. Herrera, J. Ospino, and A. di Prisco, “All static spherically symmetric anisotropic solutions of Einstein’s equations,” Physical Review D, vol. 77, no. 2, Article ID 027502, 2008.

L. Herrera, N. O. Santos, and A. Wang, “Shearing expansion-free spherical anisotropic fluid evolution,” Physical Review D, vol. 78, no. 8, Article ID 084026, 10 pages, 2008.

L. Herrera, A. di Prisco, J. L. Hernandez-Pastora, and N. O. Santos, “On the role of density inhomogeneity and local anisotropy in the fate of spherical collapse,” Physics Letters A, vol. 237, no. 3, pp. 113–118, 1998.

L. Herrera and N. O. Santos, “Thermal evolution of compact objects and relaxation time,” Monthly Notices of the Royal Astronomical Society, vol. 287, no. 1, pp. 161–164, 1997.

L. Herrera, A. di Prisco, J. Martin, J. Ospino, N. O. Santos, and O. Troconis, “Spherically symmetric dissipative anisotropic fluids: a general study,” Physical Review D, vol. 69, no. 8, Article ID 084026, 2004.

L. Herrera and W. Barreto, “Relativistic gravitational collapse in comoving coordinates: the post-quasistatic approximation,” International Journal of Modern Physics D, vol. 20, no. 07, pp. 1265–1288, 2011.

L. Herrera and N. O. Santos, “Local anisotropy in self-gravitating systems,” Physics Reports, vol. 286, no. 2, pp. 53–130, 1997.

L. Herrera, J. Martin, and J. Ospino, “Anisotropic geodesic fluid spheres in general relativity,” Journal of Mathematical Physics, vol. 43, no. 10, pp. 4889–4897, 2002.

L. Herrera, “Physical causes of energy density inhomogenization and stability of energy density homogeneity in relativistic self-gravitating fluids,” International Journal of Modern Physics D, vol. 20, no. 9, Article ID 1689, 2011.

A. di Prisco, L. Herrera, G. L. Denmat, M. A. H. MacCallum, and N. O. Santos, “Nonadiabatic charged spherical gravitational collapse,” Physical Review D, vol. 76, no. 6, Article ID 064017, 2007.

A. DiPrisco, L. Herrera, and V. Varela, “Cracking of homogeneous self-gravitating compact objects induced by uctuations of local anisotropy,” General Relativity and Gravitation, vol. 29, no. 10, pp. 1239–1256, 1997.

G. Pinheiro and R. Chan, “Radiating gravitational collapse with shear viscosity revisited,” General Relativity and Gravitation, vol. 40, no. 10, pp. 2149–2175, 2008.

G. Pinheiro and R. Chan, “Radiating gravitational collapse with an initial inhomogeneous energy density distribution,” General Relativity and Gravitation, vol. 43, no. 5, pp. 1451–1467, 2011.

G. Pinheiro and R. Chan, “Radiating shear-free gravitational collapse with charge,” General Relativity and Gravitation, vol. 45, no. 1, pp. 243–261, 2013.

R. Chan, “Radiating gravitational collapse with shear viscosity,” Monthly Notices of the Royal Astronomical Society, vol. 316, no. 3, pp. 588–604, 2000.

R. Chan, “Radiating gravitational collapse with shearing motion and bulk viscosity,” Astronomy and Astrophysics, vol. 368, no. 1, p. 325, 2001.

R. Chan, M. F. A. da Silva, and J. F. V. da Rocha, “Gravitational collapse of self-similar and shear-free fluid with heat flow,” International Journal of Modern Physics D, vol. 12, article 347, no. 3, 2003.

R. Chan, L. Herrera, and N. O. Santos, “Dynamical instability for shearing viscous collapse,” Monthly Notices of the Royal Astronomical Society, vol. 267, p. 637, 1994.

R. Tikekar and L. K. Patel, “Non-adiabatic gravitational collapse of charged radiating fluid spheres,” Pramana, vol. 39, no. 1, pp. 17–25, 1992.

S. D. Maharaj and M. Govender, “Collapse of a charged radiating star with shear,” Pramana, vol. 54, no. 5, pp. 715–727, 2000.

M. Sharif and K. Iqbal, “Spherically symmetric gravitational collapse,” Modern Physics Letters A, vol. 24, no. 19, Article ID 1533, 2009.

M. Sharif and A. Siddiqe, “Singularity in gravitational collapse of plane symmetric charged vaidya spacetime,” Modern Physics Letters A, vol. 25, no. 33, Article ID 2831, 2010.

M. Sharif and A. Siddiqe, “Dynamics of charged plane symmetric gravitational collapse,” General Relativity and Gravitation, vol. 43, no. 1, pp. 73–91, 2011.

M. Sharif and G. Abbas, “Charged perfect fluid cylindrical gravitational collapse,” Journal of the Physical Society of Japan, vol. 80, Article ID 104002, 2011.

M. Sharif and G. Abbas, “Dynamics of non-adiabatic charged cylindrical gravitational collapse,” Astrophysics and Space Science, vol. 335, no. 2, pp. 515–521, 2011.

M. Sharif and S. Fatima, “Charged cylindrical collapse of anisotropic fluid,” General Relativity and Gravitation, vol. 43, no. 1, pp. 127–142, 2011.

W. Barreto, B. Rodríguez, L. Rosales, and O. Serrano, “Self-similar and charged radiating spheres: an anisotropic approach,” General Relativity and Gravitation, vol. 39, no. 1, pp. 23–39, 2007.

C. Ghezzi, “Relativistic structure, stability, and gravitational collapse of charged neutron stars,” Physical Review D, vol. 72, Article ID 104017, 2005.

R. Goswami and P. S. Joshi, “Gravitational collapse of an isentropic perfect fluid with a linear equation of state,” Classical and Quantum Gravity, vol. 21, no. 15, Article ID 3645, 2004.

R. Goswami and P. S. Joshi, “What role do pressures play in determining the final end state of gravitational collapse?” Classical and Quantum Gravity, vol. 19, article 129, no. 20, 2002.

R. Goswami and P. S. Joshi, “Black hole formation in perfect fluid collapse,” Physical Review D, vol. 69, Article ID 027502, 2004.

S. Chakraborty and T. Bandyopadhyay, “Collapsing inhomogeneous dust in the background of perfect (or anisotropic) fluid,” International Journal of Modern Physics D, vol. 17, no. 08, pp. 1271–1281, 2008.

D. M. Eardley and L. Smarr, “Time functions in numerical relativity: marginally bound dust collapse,” Physical Review D, vol. 19, no. 8, pp. 2239–2259, 1979.

F. C. Mena, B. C. Nolan, and R. Tavakol, “Role of anisotropy and inhomogeneity in Lemaitre-Tolman-Bondi collapse,” Physical Review D, vol. 70, no. 8, Article ID 084030, 2004.

S. Thirukkanesh and S. D. Maharaj, “Charged relativistic spheres with generalized potentials,” Mathematical Methods in the Applied Sciences, vol. 32, no. 6, pp. 684–701, 2009.

K. S. Govinder, M. Govender, and R. Maartens, “On radiating stellar collapse with shear,” Monthly Notices of the Royal Astronomical Society, vol. 299, no. 3, pp. 809–810, 1998.

D. Schäfer and H. F. Goenner, “The gravitational field of a radiating and contracting spherically-symmetric body with heat flow,” General Relativity and Gravitation, vol. 42, no. 11, pp. 2119–2130, 2000.

B. V. Ivanov, “Collapsing shear-free perfect fluid spheres with heat flow,” General Relativity and Gravitation, vol. 44, no. 7, pp. 1835–1855, 2012.

B. V. Ivanov, “Self-gravitating spheres of anisotropic fluid in geodesic flow,” International Journal of Modern Physics D, vol. 20, no. 3, article 319, 2011.

B. V. Ivanov, “Evolving spheres of shear-free anisotropic fluid,” International Journal of Modern Physics A, vol. 25, no. 20, pp. 3975–3991, 2010.

S. Sarwe and R. Tikekar, “Non-adiabatic gravitational collapse of a superdense star,” International Journal of Modern Physics D, vol. 19, no. 12, pp. 1889–1904, 2010.

R. Sharma and R. Tikekar, “Space-time inhomogeneity, anisotropy and gravitational collapse,” General Relativity and Gravitation, vol. 44, no. 10, pp. 2503–2520, 2012.

R. Sharma and R. Tikekar, “Non-adiabatic radiative collapse of a relativistic star under different initial conditions,” Pramana, vol. 79, no. 3, pp. 501–509, 2012.

B. V. Ivanov, “The importance of anisotropy for relativistic fluids with spherical symmetry,” International Journal of Theoretical Physics, vol. 49, no. 6, pp. 1236–1243, 2010.

R. L. Bowers and E. P. T. Liang, “Anisotropic spheres in general relativity,” Astrophysical Journal, vol. 188, pp. 1262–1274, 1974.

B. V. Ivanov, “Maximum bounds on the surface redshift of anisotropic stars,” Physical Review D, vol. 65, no. 10, Article ID 104011, 2002.

M. K. Mak and T. Harko, “Anisotropic stars in general relativity,” Proceedings of the Royal Society A, vol. 459, no. 2030, pp. 393–408, 2003.

T. Harko and M. K. Mak, “Anisotropic relativistic stellar models,” Annalen Der Physik, vol. 11, no. 1, pp. 3–13, 2002.

K. Dev and M. Gleiser, “Dark energy and the statistical study of the observed image separations of the multiply-imaged systems in the class statistical sample,” International Journal of Modern Physics D, vol. 13, no. 06, 2004.

S. Karmakar, S. Mukherjee, R. Sharma, and S. D. Maharaj, “The role of pressure anisotropy on the maximum mass of cold compact stars,” Pramana, vol. 68, no. 6, pp. 881–889, 2007.

J. Martinez, “Transport processes in the gravitational collapse of an anisotropic fluid,” Physical Review D, vol. 53, no. 12, pp. 6921–6940, 1996.

P. C. Vaidya and R. Tikekar, “Exact relativistic model for a superdense star,” Journal of Astrophysics and Astronomy, vol. 3, no. 3, pp. 325–334, 1982.

M. R. Finch and J. E. F. Skea, “A realistic stellar model based on an ansatz of Duorah and Ray,” Classical and Quantum Gravity, vol. 6, no. 4, article 467, 1989.

S. Hansraj and S. D. Maharaj, “Charged analogue of finch-skea stars,” International Journal of Modern Physics D, vol. 15, no. 08, Article ID 1311, 2006.

R. Tikekar and K. Jotania, “A relativistic two-parameter core-envelope model of compact stars,” Gravitation and Cosmology, vol. 15, no. 2, pp. 129–133, 2009.

R. Tikekar and K. Jotania, “On relativistic models of strange stars,” Pramana, vol. 68, no. 3, pp. 397–406, 2007.

K. Jotania and R. Tikekar, “Paraboloidal space-times and relativistic models of strange stars,” International Journal of Modern Physics D, vol. 15, no. 8, pp. 1175–1182, 2006.

W. Israel and J. Stewart, “Transient relativistic thermodynamics and kinetic theory,” Annals of Physics, vol. 118, no. 2, pp. 341–372, 1979.

R. Maartens, “Dissipative cosmology,” Classical and Quantum Gravity, vol. 12, no. 6, Article ID 1455, 1995.