Radiating-collapsing models satisfying Karmarkar condition

This paper presents a class of exact spherical symmetric solutions of the Einstein equations admitting heat-conducting anisotropic fluid as a collapsing matter. The exterior spacetime is assumed to be the Vaidya metric. This class of solutions is shown to satisfy all the energy conditions throughout the interior of the star, and the luminosity is time independent, radiating uniformly throughout the collapse.


Introduction
There has been extensive research in the field of gravitational collapse. Since the pioneering work on the gravitational collapse of homogeneous dust, [1] [2], it is now accepted that for a gravitational collapse of homogeneous pressureless matter, the central singularity remains hidden behind the horizon, implying that the end state of the continual gravitational collapse of homogeneous dust cloud must be a black hole [3]. Further studies have examined various aspects of gravitationally collapsing stellar systems for different kinds of matter distributions, and details may be found in [4,5] [6]- [10]. These studies have thrown light to many interesting facts which must hold for the collapse processes to be physically realistic. For example, for the continuous and smooth matching of the interior collapsed spacetime to the exterior Vaidya spacetime over the timelike hypersurface Σ, the radial pressure must not vanish at the boundary of the collapsing radiant star, but instead be proportional to the heat flux [12] [13].
In order to have a physically well behaved model of a gravitation collapse with generic energy momentum tensors, one not only needs to find physically consistent analytical solutions of the Einstein field equations, but also must ensure validity of the energy conditions as well. The practice usually followed are: to specify the spacetime symmetry, or the gravitational potentials, imposing an equation of state, or restricting the matter content to find the solutions of the gravitational collapse. However, there exists a class I condition, which is useful to obtain classes of solutions. This condition arise from some well known geometric theorems as follows: First, an (n + 1)-dimensional space V n+1 can be embedded into a pseudo Euclidean space E n+2 of dimension (n+2) [14], and that all the spherically symmetric spacetime are in general of class II. Next, the necessary and sufficient condition for any Riemannian space to be embedding class I is that it satisfies the Karmarkar condition [15] [16]. Thus the Karmarkar condition is a useful condition which gives new solutions.
Recently, there has been a renewed interest in obtaining solutions using these conditions, [17]- [19]. The study of the non-static radiating metric with timelike Karmarkar condition, when the temporal dependence of the model is linear, shows that the collapse proceeds without the formation of the horizon [20]. The assumption is that the metric coefficients (which are in general, functions of r and t), be separable, for example g rr = b(r)f (t) (see the equation (2.1) for example). In their study, they assume one of gravitational potential, which we denote by b(r), to be constant. Next, using the results of the horizon-free collapse model [21] and a linear form of time function f (t) into the Karmarkar condition the other metric function g tt is obtained. It must be noted that this choice of b(r) is special in [20], and one may generalise. We introduce a form of pressure anisotropy ∆ = p t − p r , between the radial and the tangential pressures, and using the Karmarkar condition, show that here too the gravitational potentials may be determined. Clearly, these gravitational potentials generalise those in [20]. For these models, we study the thermodynamical evolution and their temperature profiles, since these quantities play prominent role during dissipative gravitational collapse. In fact, they decide the departure from the thermodynamical equilibrium. Quite naturally, the study of these transport processes in the context of irreversible thermodynamics in radiative gravitational collapse has gained considerable attention and [26]- [34] give further details.
The paper has been organized as follows. In section 2 we present a description of field equations of spherical symmetric anisotropic system and the junction conditions for smoothly matching of the interior spacetime with the exterior Vaidya spacetime across the timelike surface Σ. Also, we present the solutions to the Einstein field equation and the explicit expressions for physical quantities. In the light of these solutions we analyze the physical evidences of our model by verifying the energy conditions and it has been shown that all these energy conditions are well satisfied throughout the collapse. In section 3 we study the thermodynamical evolution and temperature profiles of the radiative gravitational collapse. Finally, discussion of the result accompanied with concluding remarks are given in section 4.

Fluid distribution and field equations
In comoving-synchronous coordinates, the general spherically symmetric metric is given by (2.1) We consider the interior energy-momentum tensor for the anisotropic fluid distribution with radial heat flow of the form where ρ, p r and p t are the energy density, radial pressure and tangential pressure respectively, u µ , X µ and q µ are unit time-like 4-velocity vector, unit space-like vector along radial 4-vector and radial heat flow vector respectively. These satisfy u µ u µ = −X µ X µ = −1 and u µ X µ = u µ q µ = 0.
In the comoving co-ordinates the 4-velocity and unit space-like vector and radial heat flow vector of the fluid are given by 3) The magnitude of the expansion scalar Θ for the metric (2.1) is given by The non vanishing components of Einstein-Maxwell field equations for the metric (2.1), energy momentum tensor (2.2) and (2.3) are (using units with c = 1 = 8πG) where dot and prime are the derivatives with respect to t and r respectively. From equations (2.5)-(2.8) we can see that the number of field equations are less than that of the number of unknowns. Also, the form of three metric potentials fix all the unknown physical quantities of the system. So, in order to study the collapsing phenomena we need to fix forms for these metric potentials.
Let us begin with the Israel-Darmois junction conditions [35]- [37]. The interior manifold is M − and exterior manifold is M + to be matched across the bounding timelike three space Σ, at r = r b . The exterior spacetime M + , is described by the Vaidya spacetime having outgoing radial flow of the radiation given by [11] The junction conditions require the matching of metric as well as the extrinsic curvatures where K ± ij is the extrinsic curvature to Σ given by [41] Here x l are the coordinates of interior and exterior spacetimes, z i are the coordinates that defines the hypersurface Σ and n i are the unit normal vector to Σ. The junction condition on metric functions, given by (2.10) at the hypersurface Σ gives where τ is the time coordinate defined only on the hypersurface Σ. To match the extrinsic curvatures, we need the normal vector fields. The unit normal vectors on the hypersurface Σ for the interior and exterior spacetime are given by The non vanishing components of the extrinsic curvature for metrics (2.1) and (2.9) are given by Now, from the second junction condition (2.11), one must have the equality for the θθ components since the spherical part matches for the interior and the exterior. Thus, the equation K + θθ = K − θθ at hypersurface Σ, and along with equations (2.13), (2.14) and (2.15) give where 2m is the total energy entrapped inside the hypersurface Σ [39,40]. Now, again from the matching condition (2.11), the matching of the Substituting equations (2.13), (2.14) and (2.23) into the equation (2.22) we have Now, differentiating (2.25) with respect to the τ and using equations (2.23) and (2.25), we can write the equation (2.24) and comparing with equations (2.6) and (2.8) we have the useful equation We also require the total Luminosity for an observer at rest at infinity is given by [9] Differentiating the equation (2.23) with respect to t and using equation (2.6), (2.13) and (2.24), above equation (2.27) becomes (2.28) The boundary redshift can be used to determine the time of formation of the horizon. The boundary redshift Z Σ is given by [9] dv dτ Σ = 1 + Z Σ .
Since the number of field equations are more than that of the number of unknowns in the system. These unknowns depend upon the form of metric potentials. To find the metric potentials, we take the pressure anisotropy ∆ to be: The full expression of anisotropy is given in equation (2.30). However, this form is too complicated to be solved in full detail. So, to simplify, we separate a set of terms which are significant in the following sense: The anisotopy in pressure ∆ is so assumed that the model is physically significant, such that ∆ vanishes at the center r = 0 of the cloud and is regular towards the boundary. Furthermore, this choice makes the original pressure anisotropic equation (2.30) as differential equation of only one function, which one can easily integrate to find solution. The vanishing part of the equation (2.30) which gives the following solution for the function b(r): where C 3 and C 4 are constant of integration.
One may now obtain time dependent solutions. It is known that an (n + 1)-dimensional space V n+1 can be embedded into a pseudo Euclidean space E n+2 of dimension (n + 2) if there exists a symmetric tensor b µν which satisfies the Gauss-Codazzi equations [14].
where e = ±1 (+ or −, when the normal to the manifold is spacelike or timelike respectively) and b µ ν are the coefficient of the second differential form. The necessary and sufficient condition for any Riemannian space to be an embedding class I is the Karmarkar condition [15,16], and for our case, the condition reduces to: The non vanishing components of the Riemann tensor for the shear free spherically symmetric metric (2.1) are Although we have found the form of one of the metric potential b(r) from the choice of anisotropy factor, however, the class I condition equation (2.44) is still nonlinear in its temporal and radial behavior. For a collapsing model both the conditions (2.26) and (2.44) must be simultaneously satisfied. It has been found that one of the solutions of (2.44) is a linear solution [21] f where C Z > 0. Using equation (2.34) and (2.45) into the above equation (2.44) we obtain the form of a(r) as where C 1 and C 2 are integration constants.
The dependency of gravitational potentials a(r) and b(r) f (t) in perfect fluid and the Karmarkar condition leads to two exact solutions: the Schwarzschild [22] or the Kohler and Chao [23] solution. The Kohler and Chao solution is only physical for unbound configuration such as cosmological model as the radial pressure vanishes at r → ∞. The pressure isotropy condition ∆ = 0 and the Karmarkar condition (2.44) are satisfied when either C 1 , C 3 or C 4 vanish. However, for the metric potential a(r) to be positive and greater than zero, both C 3 and C 4 should not vanish throughout the collapse. This implies that in the static case, when f (t) = 1, from pressure isotropy and the Karmarkar condition we must have C 1 = 0. Thus, the for static case, Karmarkar condition together with the pressure isotropy yields the Schwarzschild like form of the gravitational potentials given by Plot of the radial heat flux q, equation (2.52), with respect to t and r. This plot shows that the heat flux is positive, implying heat is radiating throughout the collapse. As the star deviates from the equilibrium, it will start radiating heat and as the contraction of the cloud increases it will start radiating more and more heat flux. It can be seen from the graph that the heat flux starts to increase as the collapse starts at t = −100 and heat flux is positive and maximum as it reaches the end state of the collapse at t = 0. Now physical quantities (2.5)-(2.8) interms of these exact solutions of the radiating star becomes The boundary condition (2.26) in the view of (2.50)-(2.52) becomes 2ff +ḟ 2 − 2xḟ = y, (2.55) It can be seen from the Figs. 1a, 1b and 2a that density, radial pressure and tangential pressure are positive throughout the collapse. Also, Fig. 2b shows that the radial heat flux is also positive throughout collapse. Also the expansion scalar (2.4) and the Misner sharp mass function (2.23) have the form (2.57) The metric potentials a(r) and b(r) should be greater than zero throughout the collapsing phenomena. This implies for b(r) > 0, both C 3 and C 4 should be negative. For second metric potential to be positive i.e. a(r) > 0 we must have C 2 Z < C 2 . This implies at the center r = 0, C Z < C 1 − 2C 2 C 3 C 4 . Also, the condition of positivity of the radial pressure p r throughout the collapsing phenomena imposes constraints on the range of C 1 as C 2 C 3 C 4 < C 1 < 2C 2 C 3 C 4 .
Using equations (2.34), (2.45), (2.46) and (2.50) into the (2.28), we found that the Luminosity of the radiating collapse becomes time independent and vanishes if R = (b + rb ) = 0. Which implies that the red shift diverges at the time of formation of the blackhole i.e. at R = 0. As we know that R implies the shell crossing singularity and for a smooth collapse there should be no shell crossing during the collapse, which means shell crossing singularity should occur either at the same or later epoch than that of shell focusing singularity( i.e at R = 0). As the red shift diverges at R = 0, this means the only place where blackhole can form is at the central singularity where both the shells crossing and shell focusing singularity forms simultaneously. This is because for a physically realistic smooth model there should be R > 0 throughout the collapse. It is clear from the Fig. 3a that the expansion scalar Θ is negative throughout implying the contracting behavior of the system. Also, Fig. 3b shows that the mass is linear function of time t.
Let us now ask if the solutions also satisfy the energy conditions as well. Energy conditions plays important role in the study of the astronomical phenomena like collapsing stellar models. In this section we will analyze the physical evidences of our model by verifying the energy conditions. The energy conditions namely weak energy condition (WEC), null energy condition (NEC), dominant energy condition (DEC) and strong energy conditions (SEC) will be satisfied at all points in the stellar model if the following inequalities are satisfied simultaneously [9] [22] E1 : (ρ + p r ) 2 − 4q 2 ≥ 0 (SEC/DEC/WEC) E6 : 2p t + (ρ + p r ) 2 − 4q 2 ≥ 0 (SEC) Beside these energy conditions, a physically reasonable stellar model should also satisfy E7 : ρ > 0, p r > 0, p t > 0, and ρ < 0, p r < 0, p t < 0.
Here, we can see that the validity of the E1 and E2 inequalities implies that the E4 inequality is satisfied. In the same fashion, the validity of the E1, E2 and E7 inequalities ensures that the E5 and E6 inequalities are satisfied. So, in general we only need to check the validity of the E1, E2, E3 and E7. For our radiating stellar model, it can be seen from the figures 4a, 4b and 4c that all these energy conditions are well satisfied throughout the interior of the collapsing star.
We would like to shed some light on the stability of our model by investigating the shear free condition. It has been found that the evolution of the shear depends upon the scalar Y T F which is given by [24] where m is Misner Sharp mass function given by (2.23) or interms of the density and heat flux it can be written as [24] m = 1 2 Here, R = r b(r) f (t) is the radius of the collapsing cloud. Equation (2.60) shows that the scalar Y T F depends upon the density of the cloud ρ, pressure anisotropy term ∆ and mass function m. It has been shown that for a geodesic fluid, the presence of the inhomogeneous density ρ, anisotopy term ∆ and heat flux q can generates the shearing effects within the fluid. To check the stability of the shear free condition for our model, we use the forms of the pressure anisotropy (2.31), density (2.49) and mass function into the equation (2.60), which gives Y T F = ∆. This shows that although we have consider the initial shear free collapse, at later time, the presence of the anisotropy term ∆ generates the shearing effects in the fluid distribution.

Thermal properties
The study of the thermodynamical evolution and the temperature profiles of the radiating stars plays prominent role during dissipative gravitational collapse as they decide the departure from the thermodynamical equilibrium. Previous studies of the shear free and shearing dissipative gravitational collapses shows that the relaxation effects plays significant role in the temperature profiles towards the end state of the dissipative gravitational collapse see [27,28,29,32,42] and references therein. To study the thermodynamical evolution and the temperature profiles of inside the collapsing star, we will use the causal transport equation for the metric (2.1) given by [26,27,43] where, α > 0, β > 0, γ > 0 and σ > 0 are constants and h µν = g µν + u µ u ν . Also, are physically reasonable choices of the mean collision time, between massive and massless particles τ c , thermal conductivity k and the relaxation time τ respectively [27] [34]. τ represents the causality index, measures the strength of relaxational effects and τ = 0 or β = 0 represents the noncausal case. Using these forms of the physical quantities given in equation (3.3), the form of the causal heat transport equation (3.2) becomes The noncausal solution of heat equations are obtained by setting β = 0 i.e. τ = 0 in the above transport equation (3.4) [34] The causal solution of the above transport equation (3.4) are given by [34] (a T ) 4 where F (t) is the function of integration. The function F (t) is determined by invoking boundary conditions where L ∞ is the total luminosity for an observer at infinity and δ > 0 is constant. As with previous Figure 5: Plot of the causal and noncausal temperature T of the star with respect to r for σ = 0.
investigations [31,34], Fig.5 shows that both the causal and noncausal temperature are same at the boundary of the star. However, at later stages of the collapse, relaxation effects plays significant role and they differ from the noncausal case. This behavior can be seen from the Fig. 5 that with β > 0 the relaxations effects grows and the causal temperature remains greater than that of noncausal temperature throughout the interior of the star. These results are in agreement with the earlier results obtained for the shear free collapse [30,32].

Conclusion
In this paper, we presented new exact solutions of the Einstein field equations for spherical symmetric systems with heat-conducting anisotropic fluid as a collapsing matter. The interior spacetime has been smoothly matched with the exterior Vaidya metric of radiant star across the timelike hypersurface Σ. For this purpose, we assume the gravitational potentials to be separated into their radial and temporal coordinates. Now to find their exact forms, we consider a special form of the pressure anisotropy as given by the equation (2.31) and find one of the metric potential. To find the other gravitational potential, we employed the Karmarkar condition (2.44), which makes the interior spacetime to be class I.
We have investigated the physical quantities like density (2.49), radial pressure (2.50) and tangential pressure (2.51) and it can be seen from the Figs. 1a, 1b and 2a that they are positive throughout the collapse. From Fig. 2b it is clear that the radial heat flux (2.52) is finite and positive throughout collapse. From Fig. 3a, the negative form of the expansion scalar (2.56) implies the collapsing behavior of the system. We have also investigated the mass function and total luminosity of the collapsing star. It is clear from the Fig. 3b that the mass (2.57) depends linearly on the temporal coordinate. It has been found that the luminosity is time independent for this class of solutions and it radiates uniformly throughout the collapse. Thus the solution is a physically permissible solution of the Einstein theory. We also checked the physically viability of the collapsing model by studying the energy conditions E1, E2 and E3 together with the positivity of the density, radial and tangential pressure profiles and heat flux. It is clear from the figures 4a, 4b and 4c that all three energy conditions E1, E2 and E3 are satisfied. It shows that the present class of solutions seams to be representing the physically viable collapse model. Also, we have found the surface temperature of the collapsing star at a large past time. It is clear from the Fig. 5 that both the causal and noncausal temperatures are same at the boundary Σ, however, differs at all interior points of collapsing star. At later stages of the collapse, relaxation effects (3.7) plays significant role and remains greater than that of the noncausal case (3.5). From Fig. 5 shows that with β > 0 the relaxations effects grows and the causal temperature remains greater than that of noncausal temperature throughout the interior of the star. These results are in agreement with the earlier results obtained for the shear free collapse.