Abstract

The present works deals with gravitational collapse of cylindrical viscous heat conducting anisotropic fluid following the work of Misner and Sharp. Using Darmois matching conditions, the dynamical equations are derived and the effects of charge and dissipative quantities over the cylindrical collapse are analyzed. Finally, using the Miller-Israel-Steward causal thermodynamic theory, the transport equation for heat flux is derived and its influence on collapsing system has been studied.

1. Introduction

A challenging but curious issue in gravitational physics as well as in relativistic astrophysics is to know the final fate of a continual gravitational collapse. The stable configuration of a massive star persists as long as the inward pull of gravity is neutralized by the outward pressure of the nuclear fuel at the core of the star. Subsequently, when the star has exhausted its nuclear fuel there is no longer any thermonuclear burning and there will be endless gravitational collapse. However, depending on the mass of the collapsing star, the compact objects such as white dwarfs, neutron stars, and black holes are formed. White dwarf and neutron star gravity are counterbalanced by electron and neutron degeneracy pressure, respectively, while black hole is an example of the end state of collapse.

The study of gravitational collapse was initiated long back in 1939 by Oppenheimer and Snyder [1]. They have studied the collapse of a homogeneous spherical dust cloud in the frame work of general relativity. Then after a quarter century, a more realistic investigation was done by Misner and Sharp [2] with perfect fluid in the interior of a collapsing star. In both the studies, the exterior of the collapsing star was chosen as vacuum. Vaidya [3] formulated the nonvacuum exterior of a star having radiating fluid in the interior. An inhomogeneous spherically symmetric dust cloud was analytically studied by Joshi and Singh [4] and they have shown that the final fate of the collapsing star depends crucially on the initial density profile and the radius of the star. Debnath et al. [5] investigated collapse dynamics of the nonadiabatic fluid, considering quasi-spherical Szekeres space-time in the interior and plane symmetric Vaidya solution in the exterior region.

Although most of the works on collapse dynamics are related to spherical objects, still there is interesting information about self-gravitating fluids for collapsing object with different symmetries. The natural choice for nonspherical symmetry is axis-symmetric objects. The vacuum solution for Einstein field equations in cylindrically symmetric space-time was obtained first by Levi-Civita [6], but still it is a challenging issue of interpreting two independent parameters in the solution. Herrera and Santos [7] studied cylindrical collapse of nondissipative fluid with exterior Einstein-Rosen space-time and showed wrongly a nonvanishing radial pressure on the boundary surface and subsequently in collaboration with MacCallum [8] they corrected the result. Then Herrera and collaborators investigated cylindrical collapse of matter with [9] or without shear [10].

Further, the junction conditions due to Darmois [11] has a very active role in dealing collapsing problems. Sharif et al. [1214] showed the effect of positive cosmological constant on the collapsing process by using junction conditions between static exterior and nonstatic interior with a cosmological constant. Also Herrera et al. [15], using junction conditions, were able to prove that any conformally flat cylindrically symmetric static source cannot be matched to the Levi-Civita space-time. Then Kurita and Nakao [16] formulated naked singularity along the axis of symmetry, considering cylindrical collapse with null dust.

Moreover, from realistic point of view, it is desirable to consider dissipative matter in the context of collapse dynamics [1719]. Considering collapse of a radiating star with dissipation in the form of radial heat flow and shear viscosity, Chan [20] has showed that shear viscosity plays a significant role in the collapsing process. Collapse dynamics with dissipation of energy as heat flow and radiation have been studied by Herrera and Santos [18]. Subsequently, by considering of causal transport equations related to different dissipative components (heat flow, radiation, shear, and bulk viscosity), Herrera et al. [15, 21, 22] investigated the collapse dynamics. The same collapsing process with plane symmetric geometry or others has been examined by Sharif et al. [23, 24].

On the other hand, in the context of gravitational waves, the sources must have nonspherical symmetry. Further, cylindrical collapse of nondissipative fluid with exterior containing gravitational waves shows nonvanishing pressure on the boundary surface by using Darmois matching conditions. Recently, it has been verified [25, 26] in studying cylindrical collapse of anisotropic dissipative fluid with formation of gravitational waves outside the collapsing matter.

In the present work, following Misner and Sharp, collapse dynamics of viscous, heat conducting charged anisotropic fluid in cylindrically symmetric background will be studied. The paper is organized as follows. Section 2 deals with basic equations related to interior and exterior space-time. The junction conditions are evaluated and discussed in Section 3. The dynamical equations are derived and studied in Section 4. Finally, the process of mass, heat, and momentum transfer through transport equation is discussed in Section 5.

2. Interior and Exterior Space-Time: Basic Equations

Mathematically, the whole four-dimensional space-time manifold having a cylindrical collapsing process can be written as with . Here, , the collapsing cylindrical surface, is a time-like three-surface and is the boundary of the two four-dimensional submanifolds (interior) and (exterior).

In choosing comoving coordinates, the line element can be written as follows [25, 26]:where the metric coefficients are functions of and , that is, and so on. Also due to cylindrical symmetry, the coordinates are restricted as follows:

For compact notation, we write , .

The anisotropic fluid having dissipation in the form of shear viscosity and heat flow has the energy-momentum tensor of the following form [7, 9]:Here , , , and stand for energy density, the radial pressure, the tangential pressure, coefficient of shear viscosity, and radial heat flux vector, respectively. Also and are unit time-like and space-like vectors satisfying the following relations:Moreover, the shear tensor has the following expression:where is the acceleration vector and is the expansion scalar.

For the above metric, one may choose the unit time-like vector, space-like vector, and heat flux vector in a simple form asThe shear tensor has only nonzero diagonal components aswhere ,  ,  .

Also the acceleration vector and the expansion scalar have the following explicit expressions:In the above, by notation, we have used and .

If, in addition, we assume the above fluid distribution to be charged then the energy-momentum tensor for the electromagnetic field has the following form:where the Maxwell field tensor is related to the four potential as and the evolution of the field tensor corresponds to Maxwell equations where is the four-current vector.

As the charge per unit length of the cylinder is at rest with respect to comoving coordinates so the magnetic field will be zero in this local coordinate system [27, 28]. Hence the four potentials and the four currents take the following simple form:where is the scalar potential and is the charge density.

From the law of conservation of charge, , one obtains the total charge distribution interior to radius and per unit length of the cylinder as Now the explicit form of Maxwell’s equations (11) for the interior space-time is given by A first integral of (14) gives which satisfies identically the other Maxwell equation (15). Hence one obtains the electric field intensity as .

Further, in the interior space-time , the Einstein field equations have the following explicit form:

The gravitational energy per specific length in cylindrically symmetric space-time is defined as follows [2931]: In principle, is the charge associated with a general current which combines the energy-momentum of the matter and gravitational waves. It is usually referred to in the literature as -energy for the cylindrical symmetric space-time. For cylindrically symmetric model with killing vectors, the circumference radius and specific length are defined as follows [2931]:

,  ,  so that is termed as areal radius.

For the present model with the contribution of electromagnetic field in the interior region, the -energy takes the following form:

It should be noted that the above energy is also very similar to Tabu’s mass function in the plane symmetric space-time [32].

The exterior space-time manifold of the cylindrical surface is described by the metric in the retarded time coordinate as [33, 34] where is the usual retarded time, is the total mass inside , is the total charge bounded by , and is an arbitrary constant. Further, from the point of view of the interior manifold , the bounding three-surface (comoving surface) is described as and hence the interior metric on takes the following form: where

defines the time coordinate on and by notation implies the equality of both sides on the surface .

Similarly, from the perspective of the exterior manifold, the boundary three-surface is characterized by

so that the exterior metric on takes the form Here, by notation, we write .

3. Junction Conditions

In order to have a smooth matching of the interior and exterior manifolds over the bounding three surfaces (not a surface layer), the following conditions due to Darmois [11] are to be satisfied.

(i) The continuity of the first fundamental form is

(ii) The continuity of the second fundamental form is . This implies that the continuity of the extrinsic curvature over the hypersurface [11] iswhere is given by

In the above expression for extrinsic curvature, are the components of the outward unit normal to the hypersurface with respect to the manifolds (i.e., in the coordinates ) and have explicit expressions.

and with .

Also, in the above, the Christoffel symbols are evaluated for the metric in or accordingly and we choose , , as the intrinsic coordinates on for convenience.

The continuity of the 1st fundamental form gives

Now the nonvanishing components of extrinsic curvature are

Hence, continuity of the extrinsic curvature together with (33) gives the following relations over [34, 35]:

Thus (35) gives the total mass inside the boundary surface , while (36) shows the linear relationship between the energy for the cylindrically symmetric space-time with the bounding mass over . Further, (37) shows a linear relationship among the fluid parameters on the bounding surface . Hence radial pressure is in general nonzero on the bounding surface due to dissipative nature of the fluid and the charge on the bounding surface. But when dissipative components of the fluid are switched off, then the above result (uncharged) agrees with the results of Herrera et al. [8]. Also it should be noted that the radial pressure on the boundary does not depend on the charge bounded by ; it depends only on the charge on the surface .

4. Analysis of Dynamical Equations

From the conservation of energy-momentum, that is, , we can have two zero scalars, namely, and

Using (3) and (9), the explicit expressions for these two scalars are

Now following the formulation of Misner and Sharp [2], we introduce the proper time derivative and proper radial derivative as

so that the fluid velocity in the collapsing situation can be defined as follows [36]: Using (18)–(23) and (41), we can obtain the acceleration of a collapsing matter inside as Now combining (39) and (42), we obtain

Using (23) and the junction condition , we write [24]

Hence, using the field equations for the interior manifold, we obtain the time rate of change of -energy as Also the above equation can be interpreted as the variation of the total energy inside the collapsing cylinder. Note that, due to negativity of the fluid velocity , the first term on the r.h.s will contribute to the energy of the system, provided that the radial pressure is restricted as . Due to negativity, the second term indicates an outflow of energy in the form of radiation during the collapsing process. The third term is Coulomb-like force term and it will increase the energy of the system provided .

Further, using the Einstein field equations (17), (21), and the expression for - energy in (23), the radial derivative of the energy takes the following form:

This radial derivative can be interpreted as the energy variation between the adjacent cylindrical surfaces within the matter distribution. The first term on the r.h.s. is the usual energy density of the fluid element, while the second term and third term are the conditions due to the electromagnetic field. The fourth term represents contribution due to the dissipative heat flux and the last term will increase or decrease the energy of the system during the collapse of the cylinder provided .

Finally, the collapse dynamics is completely characterized by the equation of motion in (43). Normally, for collapsing situation, should be negative, that is, indicating an inward radial flow of the system. Consequently, terms on the r.h.s (of (43)) contributing negatively favour the collapse and positive terms oppose the collapsing process. In an extreme situation, the system will be in hydrostatic equilibrium if terms of both signs balance each other. Further, from dimensional analysis, the factor can be considered as an inertial mass density, independent of heat flux contribution. The first term on the r.h.s. of (43) can be identified as the gravitational force, indicating the effects of specific length and electric charge in the gravitational contribution. The second term has three contributing components: the pressure gradient (which is negative), local anisotropy of the fluid, and electromagnetic field term. The remaining terms represent the heat flux contribution and due to negativity they seem to leave the system along the radial outward streamlines.

5. Causal Thermodynamics: The Transport Equation

In causal thermodynamics, due to Miller-Israel-Stewart, the transport equation for heat flow is given by [21]where is the projection tensor of the -surface orthogonal to the unit time-like vector , represents the thermal conductivity, is the temperature, denotes the relaxation time, and is the inertial term due to Tolman. Now, due to cylindrical symmetry, the above transport equation (47) simplifies to with .

Note that although heat dissipation is described above through a physically reasonable transport equation (47) or (48), the shear viscosity is described according to the standard (noncausal) irreversible thermodynamics. Thus, there is no such relaxation time for shear viscosity and as such no causal evolution equation for it. Further, the thermodynamic viscous/heat coupling coefficients are not taken into consideration in the present work; that is, the present approach is only partially causal. However, a full causal approach of the dissipative collapse has been considered by Herrera et al. [22].

Now considering proper derivatives in (40) of the above equation and using the field velocity (in (41)) and equation of motion (i.e., (43)), one obtains the effects of heat flux or dissipation in the collapsing process as

with

The l.h.s. of (49) can be interpreted as Newtonian force with as the inertial mass density. So as , ; that is, there is no inertial force and collapse will be inevitable due to gravitational attraction. Further, the inertial mass density decreases as long as and it increases for . Moreover, due to equivalence principle, the gravitational mass also decreases or increases according to or and gives a clear distinction between the expanding and collapsing process due to dynamics of dissipative system. Note that although the gravitational force is affected by the same factor , the hydrodynamical force is free from it. Further, combination of all these terms on the r.h.s of (49) results in the l.h.s, that is, ; there will be gravitational collapse, while there will be expansion if the l.h.s is to be positive. It should be mentioned that such effects were first discussed by Herrera et al. [37] in context of thermal conduction in hydrostatic equilibrium. Interestingly, if continuously decreases from a value larger than unity to one less than unity, then there will be a phase transition (collapse to expansion) and bounce will occur. As a result, there will be loss of energy of the system and collapsing cylinder with nonadiabatic source causes emission of gravitational radiation. Therefore, there will be radiation outside the collapsing cylinder and hence the choice of the exterior metric (in (24)) is justified.

Conflicts of Interest

The authors declare that they have no conflicts of interest.