Gravitational collapse with tangential pressure

Gravitational collapse with tangential pressure

Daniele Malafarina Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India    Pankaj S. Joshi Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India

Using the general formalism for spherical gravitational collapse developed in PSJ (), we investigate here the final fate of a spherical distribution of a matter cloud, where radial pressures vanish but tangential pressures are non-zero. Within this framework, firstly we examine the effect of introducing a generic small pressure in a well-known black hole formation process, which is that of an otherwise pressure-free dust cloud. The intriguing result we find is that a dust collapse that was going to a black hole final state could now go to a naked singularity final configuration, when arbitrarily small tangential pressures are introduced. The implications of such a scenario are discussed in some detail. Secondly, the approach here allows us to generalize the earlier results obtained on gravitational collapse with non-zero tangential pressure, in the presence of a non-zero cosmological constant. Finally, we discuss the genericity of black hole and naked singularity formation in collapse with non-zero tangential pressure. The treatment here gives a unified and complete picture on collapse final states, in terms of black hole and naked singularity formation, generalizing the earlier results obtained for this class of collapse models. Thus the role of tangential stresses towards determining collapse endstates emerges in a straightforward and transparent manner in our treatment.

Gravitational collapse, black holes, naked singularity
04.20.Dw,04.20.Jb,04.70 Bw

I Introduction

Since the first proposal of the Cosmic Censorship Conjecture (CCC) Penrose (), a great effort has been devoted in understanding the mechanism by which black holes and naked singularities can form as the end state of the complete gravitational collapse of a massive object. However, no general proof, or even a suitable mathematically rigorous formulation of the CCC is available despite the efforts of many decades. Therefore it has become essential and inevitable that gravitational collapse within the framework of general relativity be studied in detail.

The current situation is that, despite the great amount of efforts devoted in the past years in studying the CCC and its implications, the issue of final fate of complete gravitational collapse of a massive body is far from being entirely resolved. It is clear now how naked singularities have to be considered as a general feature of general relativistic physics, and that they may appear as the end-state of collapse in a broad variety of situations (see e.g. Ref () and references therein). Still, it is not yet fully understood how different kinds of matter lead to the formation of either a black hole or a naked singularity. To give an example, how the presence of pressure in the form of tangential stresses and radial stresses might direct the collapse towards the formation of black holes or naked singularities is still a matter of debate.

The first gravitational collapse scenario studied in detail within the framework of the general theory of relativity was the so called Oppenheimer-Snyder-Datt model OSD (), that represents a collapsing ball made of homogeneous pressureless dust. It is seen here that an event horizon forms before the appearance of the final spacetime singularity at the end of collapse, therefore hiding the singularity inside a black hole. Further studies of inhomogeneous dust collapse though, such as those in the Lemaitre-Tolman-Bondi (LTB) models LTB (), showed that naked singularities can indeed form during collapse, thus violating the cosmic censorship dust (). However, the dust models are idealized scenarios that ignore pressures within the cloud. Therefore, collapse models with pressures, which are more physically realistic scenarios, were examined and it was soon discovered that for the collapse with non-zero pressures also, the end state is just as likely to be a naked singularity, as for models without pressures (see e.g. fluid ()).

In particular, collapse models with a vanishing radial pressure, but with non-vanishing tangential stresses, have been studied in some detail press (). These models are of particular interest, because they provide a somewhat clear insight into the role of pressure towards determining the collapse final states.

Using the formalism developed in PSJ () to study here the class of collapse models with a non-zero tangential pressures, three main results are obtained here. Firstly, we examine the effect of introducing small tangential pressure perturbations, on the black hole formation process, for a collapsing inhomogeneous dust cloud which is otherwise pressure-free. It is seen that a matter cloud that was collapsing to a black hole final state could change its course, and instead result in a naked singularity final configuration, when small tangential pressure perturbations are allowed in an otherwise pressure-free cloud. It is seen that the opposite scenario occurs as well. Secondly, we generalize here earlier results on gravitational collapse with pressure, in the presence of a non-zero cosmological constant. Finally, from the analysis of some collapse models, some interesting insights can be obtained on the genericity or otherwise of naked singularity formation in gravitational collapse, in the presence of a non-zero tangential pressure. Using the approach developed here, it is clearly seen how the black hole and naked singularity final states in collapse are distinguished, depending on the initial configurations of the matter density, the tangential stresses profiles and the velocities of the collapsing shells. Given this initial data, the Einstein equations then fully determine the allowed evolutions that will lead to the formation either of a black hole or a naked singularity.

Various subcases of this class of collapse models have been studied in the past by different authors, as we point out below, but our treatment here gives a unified and complete picture on collapse final states in terms of black hole and naked singularity formation, thus generalizing the earlier results. Thus the role of tangential stresses towards determining the final state of a collapsing matter cloud emerges transparently in our treatment. This helps us to see in a collective manner the different results in this class of collapse models, which are otherwise somewhat scattered in a variety of approaches that do not always make it straightforward to see the overall conclusions in totality. Within this framework, the known results are easily recovered, and new scenarios and insights are obtained.

In section II we review as needed the general formalism for gravitational collapse in the presence of pressures, as restricted to the case of non-vanishing tangential pressure and vanishing radial pressure. The set of Einstein equations that governs the dynamical evolution of the collapsing cloud and the procedure to solve them is outlined. The regularity and physical reasonableness conditions that one would like to impose upon the matter cloud, in order to avoid any unphysical behaviour of the matter fields are described. The initial and boundary conditions that determine the dynamical evolution of the matter cloud are also discussed. Section III discusses the evolution of collapse, a function related to the tangent to the outgoing null geodesics at the singularity is evaluated, and it is the sign of this quantity that governs the final outcome of collapse in terms of either a black hole or naked singularity. The occurrence of trapped surfaces during collapse is discussed, and the visibility of the singularity is studied through the analysis of outgoing radial null geodesics.

The formalism is then applied in section IV to examine the effect of introducing small tangential pressures on the black hole formation process, for an otherwise pressure-free inhomogeneous dust cloud collapse. In the section V, we discuss the genericity of black holes and naked singularities obtained as collapse endstates, for this class of collapse models. The homogeneous and inhomogeneous dust cases are also briefly discussed, giving a review of known results for dust, obtained as a special case of the above treatment. In section VI, we generalize some earlier results on collapse with tangential pressure in the presence of a cosmological constant, which is asymptotically de-Sitter or anti-de Sitter spacetime. We study how the presence of a cosmological term in the equations affects the dynamical behaviour of the cloud. The study is then applied to the classes studied in section IV, by adding a cosmological term and analyzing how the dynamics is altered. This constitutes a new class of solutions whose behaviour was not studied earlier. Finally, in section VII, we outline the interesting features of the above treatment in connection with previous results obtained by other authors.

Ii Einstein Equations and Regularity conditions

The general treatment to obtain the evolution and final fate of a spherical gravitational collapse with non-zero pressures was developed in PSJ (), and it constitutes the backbone upon which the present work is developed. One of the advantages of this approach is that it shows the final stages of collapse in a rigorous and yet transparent fashion, where all the crucial elements that enter towards the determination of the final outcome are enclosed in a single straightforward equation. It is clear that in order for the singularity to be visible, some outgoing nonspacelike geodesics must come out of it to reach the faraway observers in the universe. This visibility of outgoing geodesics is shown to be related to the tangent to the null geodesics near the singularity, which in turn is shown to be dependent on the mass distribution, the pressure and the velocity profiles for the collapsing cloud. Once these quantities are set at the initial time, the evolution, as governed by the Einstein equations, determines entirely the final outcome of collapse.

In spherical coordinates, the most general line element for a spherically symmetric collapsing cloud of matter depends upon three functions and of the comoving radial coordinate and the comoving time , and it can be written as,


We shall assume here through out that the radial pressure vanishes. The energy-momentum tensor, expressed in co-moving coordinates then takes the form,


The Einstein equations then relate the three functions to the energy density and pressures of the system, and they can then be written as,


where the function is the Misner-Sharp mass of the system and we have introduced the functions and , which are defined as,


Since there is a scaling gauge freedom available for the radial coordinate , we can define without any loss of generality a scale so that at the initial time we have . To this purpose, we introduce the scaling function such that,

The gravitational collapse condition is given here by . The area function can then be expressed as,


and the requirement of collapse is fulfilled once .

The singularity is then achieved by reaching (which corresponds to ), and it is realized along the singularity time curve . Note that as we approach the regular center of coordinates , when , the energy density remains bounded, although the physical radius goes to zero. Therefore, this indicates that the genuine curvature singularity, where the physical radii of all the matter shells reach a vanishing value and the energy density diverges, is reached only in the case when . This is reflected in the scaling function by the requirement that at the singularity. To avoid the presence of shell-crossing singularities (that can be shown to be gravitationally weak and thus possibly removable cross ()), we further assume that,


The shell focusing singularity at is a true curvature singularity, and therefore is not removable from the spacetime, and its visibility or otherwise is the main issue that concerns us here.

From equation (II.3), it follows that the Misner-Sharp mass in this case does not depend on time (), and therefore it must be a function of only. We are then left with four equations in the six unknown .

Collapse models with such a vanishing radial pressure, but non-zero tangential stresses were studied in the mass-area coordinates by Magli Magli (), where the analytic structure of these solutions was derived. The algebraic equation which governs the nature of the central singularity was written and its dependence upon the initial configuration and the equation of state for the matter field was analyzed. Still, given the non-trivial nature of the problem, the task to study the outcome of collapse in full generality proved to be non-attainable. Further work by Gonçalves, Jhingan and Magli Gonc-Jhin-Mag () studied the final stages of collapse in the limited case of marginally bound collapse while the analysis of one of us press () concentrated on models with some restrictions on the pressure profiles (namely, assuming for the metric function ). Also, the work by Harada, Nakao and Iguchi Nakao () analyzed the nakedness and strength of the central singularity obtaining some general conditions which were later specialized to the dust and ‘Einstein cluster’ scenarios. In this work the nakedness and curvature strength of the shell focusing singularity was studied in detail, providing a classification of the curvature strength of the singularity based on some conditions (namely, the strong curvature condition, gravity dominance condition, and the limiting focusing condition).

Another interesting model with vanishing radial stresses that has been studied is the so called ‘Einstein cluster’. It is a specific subclass of models with tangential stresses that is constituted by a collapsing cloud of counter rotating particles. Since the tangential stress is due to the particles’ angular momentum, the ‘Einstein cluster’ is a particularly interesting model, because it reproduces the effects of rotation without departing from spherical symmetry. The analysis of the singularity forming in the ‘Einstein cluster’ cluster () shows once more how departure from dust models to include pressures does not forbid the formation of naked singularities.

With the definitions (II.10), (II.8), we substitute the unknowns with . The definition of provides one of the two free functions and determines the energy density via equation (II.4). We can choose as the second free function and therefore evaluate through equation (II.5), which can be written in the form,


Without any loss of generality, we can obtain from (II.6) once we rewrite it as,


and we define a suitable function from


Then equation (II.13) can be integrated to obtain,


The arbitrary function above that results from the integration can be interpreted following the analogy with dust LTB models. Then turns out to be related to the velocity of the collapsing shells (more precisely to the velocity of the infalling particles at the boundary of the cloud), and once we write it as


we can see that the values of in the dust limit correspond to the bound (), unbound () or marginally bound () Lemaitre-Tolman-Bondi models.

Writing equation (II.7) in the form of a potential,


it is possible to analyze the dynamical behaviour of different pressure models with the usual phase space diagrams of classical mechanics. Some examples will be discussed in section VI. It turns out that depending on the velocity profile and the equation describing the tangential stresses, collapsing scenarios and bouncing scenarios are both allowed Magli (). This marks a considerable difference with the LTB dust case where the potential is always decreasing, thus allowing only collapse to occur.

Finally, from equation (II.7) we get the equation of motion:


where the negative sign has been considered in order to describe gravitational collapse. In some cases equation (II.18) can be integrated to obtain the function , thus solving completely the system of Einstein equations. In other situations, considerations on the behaviour of near the center of the collapsing cloud suffice to provide the analysis of the dynamical behaviour of the system approaching the singularity.

A viable physical description of the gravitational collapse of a massive body must be subdued to certain constraints and regularity conditions in order to ensure that the physical quantities that appear in the model are well defined and regular everywhere.

Requiring regularity of the energy density at the initial surface from which the collapse commences, and at the center of the cloud throughout collapse, imposes some constraints on the classes of mass profiles that can be considered. Specifically, the mass profile must have the form,


which immediately implies the regularity of the density function at the center of the cloud at all regular epochs. By this prescription, we ensure that the Misner-Sharp mass vanishes at the center of the cloud and that the energy density away from the singularity has a finite value at the center. A further condition that may be imposed sometimes is that the energy density does not present a cusp at , and that it is a smooth function of the radial coordinate . This can be achieved by requiring that the linear term of the energy density expansion near the center vanishes, which is reflected in the requirement that at . Of course, this is an additional requirement and in a general physically realistic scenario, this may or may not be realized. In any case, requiring smoothness brings in some additional mathematical simplification.

The requirement that there is no pressure gradient at implies that the tangential pressure must become equal to the radial pressure (i.e. zero in our case) at the center. This imposes the condition that goes as at least, near the center, and implies regularity for the metric function via equation (II.5). This condition then implies that we have in general,


with being a regular function and . Regularity of implies regularity of the function via equation (II.14), and we have,


with being again a regular function. The relation between and is given by the equation (II.14), which can be written as


As we have seen before, a specific behavior of is required in order to avoid shell crossing singularities. The choice of satisfying implies further regularity for the energy density. This in turn gives some constraints on the non-vanishing terms in the expansion of the energy density at a constant time surface. Namely, we find that in order to avoid the shell crossing singularities, we must require , which, from a physical point of view, is a reasonable condition since we expect the density to be higher at the center of the cloud, and decreasing as we move away from the center. In the case when , turns out to be a decreasing function from the center outwards if .

As noted, from the regularity of at the center of the cloud, we get the behaviour of as in equation (II.16), and this turns out to be related to the velocity profile, and thus the kinetic energy, of the infalling particles. Positivity of is thus reflected in the condition that


Furthermore, the requirement that implies that moving away from the center, the contribution of the kinetic energy to the total energy of the collapsing cloud decreases cross (). Therefore, the requirement of avoidance of shell crossing implies that some further condition must be imposed upon . In general, from the Misner-Sharp mass we get


which, in the limit of the dust case, becomes,


The marginally bound case () is obviously allowed, and we see that considering only the constant values for the function (which correspond in the dust case to the bound and unbound LTB models), the bound case, i.e. might imply the occurrence of shell crossing if . The bound case is therefore allowed under the restrictions imposed as above. We note from equation (II.23) that the permitted values for have a limited range, namely , where is the boundary of the cloud (an infinite cloud cannot therefore have a bound velocity profile). On the other hand, the unbound models, which are also allowed, describe a collapse scenario that is not uniquely gravitational. Some additional energy is impressed to the infalling particles from outside, causing the particles to have a non zero velocity at all times. We note that in the unbound case, solutions with at the initial time are not possible.

In the more general situation where , the fact that certain velocity profiles are allowed or not might depend on the boundary conditions imposed on the matter cloud. Therefore, for some specific cases, both positive and negative functions can be considered, without causing the presence of shell crossings.

At this point, we note that the requirement that the regularity conditions must be satisfied constrains the class of allowed models and the allowed equation of state. Regularity would be reasonable to demand from the point of view of the physical validity, since it guarantees that the metric functions and physical quantities are well behaved at any regular epochs of collapse, and at the center. We see that this imposes certain non-trivial restrictions on the possible choice of the matter models. In fact, we notice how regularity and the choice of a linear equation of state of the type , as was assumed in some previous works (see for example, Singh ()), are mutually exclusive. This is immediately seen from equation (II.5) which now becomes,


and once integrated and evaluated at the initial time it gives,


Thus the metric function diverges as tends to zero, while on the other hand regularity demands that at all times, near the center. We thus see that in comoving coordinates the assumption of a linear equation of state leads to a breakdown of the coordinate system at the center, or the metric functions not being regular. This would imply that the initial data set is not regular, or that a more suitable coordinate system is required. The inspection of the Kretschmann scalar for this model shows that for any , the central shell presents a curvature singularity (where diverges) at all times, and thus the initial data surface is not regular. This is confirmed also by the divergence at , of the energy density , as evaluated from equation (II.4). It is seen that at the initial epoch also, the energy density diverges necessarily at the center. For the tangential pressure models, the choice turns out to be consistent with the regularity condition only in the dust case, where .

In order for the matter model to be physically realistic, suitable energy conditions, in particular the weak energy conditions, must be satisfied. That is, the energy density as measured by local observers must be non-negative, and for any timelike vector we require,


which reduces to


Since the physical radius is always positive, to ensure the positivity of energy density we must require and or their ratio to be positive. As we have seen, positivity of ensures that there are no shell crossing singularities occurring, while the fact that is positive is granted by the choice of positive with a suitable behaviour. From equation (II.12) we see that positivity of is then given by


and is satisfied automatically (and not only) by any increasing function . The sign of is the only factor deciding positivity or negativity of the tangential pressure and it is worth noting that negative tangential stress profiles are also allowed and they satisfy energy conditions provided that the relation expressed in equation (II.30) is satisfied. Furthermore, for small values of , from equation (II.30) we see that regardless of the values taken by , there will always be a neighborhood of for which and therefore is automatically satisfied in a close neighborhood of the center.

Finally, specifying the initial conditions for the collapsing matter cloud consists in prescribing the values of the three metric functions and that of the density and stress profiles at the initial time , given as functions of Initial (). These can be written as,

As noted earlier, at the initial time we have taken the scale such that , which implies that the scaling function at the initial epoch. Furthermore, from we get .

Not all of these initial value functions can be chosen arbitrarily. In fact, the scaling function, the choice of the mass profile, and the Einstein equations impose some constraints. From equations (II.4) and (II.19), we see that


while from equation (II.12) we get,


Also, from equation (II.20) we can write,


In turn, can be related to the function , defined by equation (II.21), at the initial time via equation (II.14),


Finally, we can relate the initial condition for to the initial value of the function through equation (II.9) and equation (II.14),


If we want the initial collapse configuration to be not trapped, in order to see how different initial matter configurations evolve forming trapped surfaces, we must further require,


which will impose some restrictions on the possible choices of the radial boundary in accordance with . This condition reflects on the initial configuration for and since implies . Therefore to avoid trapped surfaces at the initial time the velocity of the infalling shells must satisfy


This means that the scenario where collapse commences from an equilibrium configuration where , can be taken only as a limiting case.

As we noted above, considering only the presence of the tangential pressures implies conservation of the Misner-Sharp mass during the collapse. Therefore, matching to an exterior spherically symmetric vacuum solution leads inevitably to consider the Schwarzschild metric matching (). The collapsing cloud has a compact support within the boundary taken at , and the pressure of the matter composing the collapsing star is assumed to vanish at the boundary. Matching conditions imply continuity of the metric and its first derivatives across this surface, although it should be noted that eventual discontinuities in the first derivatives of the metric across the boundary may be interpreted as an additional distribution of matter acting as a layer separating the interior from the exterior. In this case, the jump of the second fundamental form across the boundary is related to an additional stress-energy tensor layer confined to the boundary. Such a matching is in principle always possible, and together with the regularity conditions and the energy conditions, might impose some further restrictions on the initial configurations allowed.

Iii Collapse evolution and final state

As we show below, the possible future evolution of the collapsing matter cloud depends upon the choice made for the matter content, namely the two free functions that determine the mass profile and the tangential pressure profile, and the initial conditions.

In order to investigate the final outcome of collapse, we will study the time curves that lead the shell of matter labeled with a specific coordinate value to fall into the singularity. As already mentioned, the spacetime singularity is achieved for , and the divergence of the energy density (and eventually of the tangential pressure) indicates that it must be a true curvature singularity.

Given the fact that the function is monotonically decreasing in time, we can invert it and consider the time as a function of the variables and . This is equivalent to considering the area-radius coordinates as done in some earlier works (see e.g. Magli (); Nakao ()) with the relation between and given by equation (II.10). The function will identify the time at which the shell labeled by a specific value reaches the event . In this manner, the occurrence of the singularity is described by the time curve . The tangent of the singularity curve at , i.e. at the central singularity, is then seen to be related to the eventual existence of outgoing radial null geodesics escaping away from the singularity, and therefore to the possible visibility or otherwise of the singularity itself to faraway observers.

From equation (II.18) we can write,


with the minus sign chosen in accordance with , in order to describe the collapse. Integrating the above equation, while treating the radial coordinate as a constant, we get


It follows therefore that, the time taken by the shell at any given value of to reach the singularity, located at , is given by,


We note that is in general a differentiable function at least outside the singularity, which would be the case because of the regularity of all the functions involved. That is the case in many of the known collapse models. In any case, whenever it is at least a function, we can always write it as an expansion near the central shell in the form,




and , given by,


Note that the above quantities are defined only if . This may be called the reality condition for the gravitational collapse to take place. The choices of and , made from (II.16) and (II.21) imply that as , we must have , and therefore we see how cannot have constant terms or terms linear in in a close neighborhood of . This justifies the regularity conditions for and that were deducted earlier. Furthermore, we can see how the regularity condition for , that comes directly from that for , ensures that and do not diverge. Evaluating explicitly we get


By continuity, we can therefore take that the singularity curve , which corresponds to , is differentiable, or in any case even if it is only a function, we can take it that it has at least a well-defined tangent at the singularity. This is the case in important known collapse models such as the Tolman-Bondi-Lemaitre dust collapse, or the ‘Einstein cluster’ models, and also for the Oppenheimer-Snyder-Datt homogeneous dust collapse. The singularity curve is computed as the time it takes for the shell located at to reach the singularity and can then be written as,


where is the time at which the central shell becomes singular. It can be shown that for any generic initial configuration there exist classes of dynamical collapse evolutions that have a well defined singularity curve.

While the collapse evolves, the trapped surfaces may form in the collapsing cloud. If the trapped surfaces form at a time anteceding the formation of singularity then the latter will be covered, thus giving rise to a black hole at the end of the collapse. If, on the other hand, the trapped surfaces do not form, or form at a later time after the singularity, then the singularity might be visible to faraway observer Global aspects ().

The boundary of the trapped region is delimited by the apparent horizon. The equation for the apparent horizon in a spherically symmetric spacetime is given by , which for the metric (II.1) becomes . From equation (II.7) we can write it as,


The behavior of the apparent horizon for the collapsing cloud with vanishing radial pressure is analogous to that of the dust case, since in both these cases. If we note that as we must have on the apparent horizon, it is immediately clear that the only singularity that can eventually be visible in the present case is that at the center of the collapsing cloud. In fact, since it is not possible to satisfy near the singularity but away from the center, the region surrounding the singularity cannot be timelike, and therefore any singularity that might eventually form near the center with must not be visible.

The difference with the dust case is then given by the different possible behaviours of the function . It is this function, as given from equation (II.18), that determines the apparent horizon curve in dependence of the mass, velocity and the tangential stress profiles. We can see the apparent horizon as a curve given by


with . Inversely, we can also write it as a curve , which represents the time at which the shell labeled as becomes trapped.

In order to understand when the central singularity may be visible to faraway observers in the spacetime, we shall analyze the time curve of the apparent horizon which is given by,


where is the singularity time curve, whose initial point is . Near , equation (III.11) can be written in the form,


From this, it is easy to see how the presence of stresses affects the time of the formation of the apparent horizon. In fact, all the initial configurations that cause to be positive force the boundary of the trapped surfaces to form after the occurrence of the singularity, thus leaving the ‘door open’ for the null geodesics to possibly come out, away from the singularity.

We can verify the visibility of the central singularity through the behaviour of radial null geodesics. This means that in order to verify the visibility of the singularity we must check if there are any future directed null or non-spacelike geodesics that terminate in the past at the singularity (for a detailed analysis of outgoing geodesics in singular spacetimes originating from collapse see Giambo ()). It can be shown that if the singularity is radially censored then it must be completely censored. On the other hand, if there are any future directed non-spacelike curves coming out from the singularity and reaching external observers then the singularity is naked. Further, in PSJ () it was proven that is a necessary and sufficient condition for outgoing radial null geodesics to escape the singularity, at least locally. Outgoing radial null geodesics are given by the equation,


with the positive sign considered in order for to be increasing moving away from the singularity.

Now, introducing the variable , we can write equation (III.13) as


We thus look for a radial null geodesic terminating in the past at the singularity, whose equation near takes the form


with finite and positive. Such null geodesics terminate at the singularity with a definite tangent and correspond to the curve


For a specific choice of , we can solve equation (III.14) and evaluate the value of the tangent to the null geodesics at the central singularity which is given by , . Thus obtaining,


which is positive and finite at the central singularity if . It is now clear how is directly related to the tangent to the null geodesics at the singularity. In this case we conclude that the future directed radial null geodesics do come out of the singularity to reach far away observers.

Let us stress that it is the sign of the quantity that uniquely determines the visibility of the singularity, as in fact implies that the time curve of the apparent horizon is future directed increasing, and future directed null geodesics in that case do come out. Therefore the central singularity occurring at is not covered.

All those initial configurations for which , and cause to be positive will force the apparent horizon to appear after the formation of the singularity, and will allow the null geodesics to come out of the singularity that forms at the center of the collapsing cloud. We conclude that the condition not only implies that the trapped surfaces form at a later stage during collapse than the singularity, but it is also the sufficient condition for null geodesics from the singularity to be visible (at least locally) to external observers.

Iv Small pressure perturbations in dust collapse

We would now like to examine here, using the formalism above, how the final fate of the gravitational collapse is affected when pressure perturbations are introduced in a pressure-free inhomogeneous dust cloud, which was otherwise going to terminate in a black hole final state. This provides a useful insight into the role of pressure in gravitational collapse dynamics. Certain classes of collapse models with a non-zero tangential pressure are analyzed, their properties are discussed, and the collapse end-states are examined.

For the sake of definiteness, we consider the case when the Misner-Sharp mass is given by an expression near the center as,


Further expressing the function (and therefore the tangential pressure) in the form,


we have,


With these expressions, we are then able to write the quantity in dependence of the physically relevant profiles , , as,


where the function (which is related to via equation (II.14)) has been written as,


We see now immediately that the behaviour of the tangential pressure (be it positive, negative, increasing or decreasing) is here reflected in the terms and in , and can influence the final outcome of the gravitational collapse. In fact, it is the choice of the coefficients of the initial density, velocity, and stress profiles that determines the quantities , and , which are responsible for the behaviour of , which in turn determines the tangent to the singularity curve at the origin. As we noted earlier, the sign of uniquely determines the final fate of collapse in terms of either a black hole or a naked singularity. As it was shown in PSJ (), positivity of is a sufficient condition for outgoing null geodesics to come out from the singularity, thus making it visible to external observers. In some cases of physical interest it turns out that . This is likely to happen when marginally bound collapse with only quadratic terms in the expansion of energy density and pressure is considered, but it is not the only possibility. In general, whenever , the analysis of higher order terms must be carried out. In this case, we can write as


and evaluate as


We also note here that the negative tangential pressures are not forbidden by the energy conditions. Here we observe, how the presence of negative pressures, that may be due to the presence of exotic matter in the form of dark matter or some other repulsive effects, may influence the final stages of collapse, since it will affect the value of through the function given in (II.22).

The second order term in the equation for the singularity curve is given by equation (IV.7). It turns out that this term becomes relevant in many physically realistic scenarios which consider the mass profile and the density profile to be written with only quadratic terms in . From the regularity conditions analyzed earlier, we can further see that (that is reflected in the requirement that ) is a reasonable assumption. This is consistent with the analysis of energy density profiles, that, for physically realistic scenarios, have only quadratic terms in the expansion. The energy density , can thus be written as .

In the situation where mass and density profiles are given only in quadratic terms, we then have,


We can now examine a few physically relevant classes of collapse models with such a property, together with the prescription for marginally bound collapse (). This choice of the mass, pressure profile and velocity profile causes to vanish and the final outcome of collapse is decided by the term, which in this case reads,


The analysis is then carried out along the same lines described above.

Firstly, an important class of collapse models, also to be noted as the reference frame, is that of an inhomogeneous dust cloud. Inhomogeneous dust has been widely studied and the occurrence of naked singularities in such cases is well-known. From the framework given above, we recover the inhomogeneous dust when , and that together with the requirement that only quadratic terms appear in the density, yields


The collapse will end in a naked singularity for all negative values of . Further, in this case, the structure of the apparent horizon can be evaluated explicitly. Note that for the above case reduces to the Oppenheimer-Snyder-Datt collapse scenario, describing homogeneous dust, and all the terms s vanish as the singularity is simultaneous.

Let us now consider the effect of introducing small tangential pressure perturbations in an otherwise pressure-free inhomogeneous dust cloud collapse, that was going to a black hole final state. To this aim, we choose a class of collapse models restricted by the assumption that the second order term in the mass profile vanishes. We see that if then , and the final outcome of collapse for the inhomogeneous dust cloud is then decided by the next order term, namely . Exactly in the same manner as described above, all positive values of will cause the collapse to end in a black hole. If now we introduce an arbitrarily small tangential pressure of the form (IV.3), we see that near this can cause the collapse to end in a naked singularity, whenever the function is chosen in such a way that , as given by equation (IV.9) with , is positive. We therefore have constructed a class of models of small tangential pressure perturbations of the LTB collapse scenario that drastically changes the final outcome of collapse on introduction of small pressures.

A second illustrative class of models is that of a ‘quasi-Hookean’ pressures, which is given by a specific choice of the tangential stresses corresponding to


for . This choice leads to an energy function of the form


In this case, putting again , we get

Some more considerations on the dynamical behaviour of this particular case will be presented in section VI also. What we can see from the above is that the values of the constants and are responsible for the final outcome of collapse in terms of black hole or naked singularity, and that either of the outcomes are possible depending on the choices of these parameters.

Finally, it is interesting to study the behaviour of a subclass of the previous model, choosing the class of small tangential stresses in the form,


In this case we see that , and near we have and . Then the final outcome of collapse is decided by the sign of


We see that in this case it is possible that the introduction of a suitable tangential stress in the inhomogeneous dust cloud will cause it to end in a naked singularity, while the same model with a vanishing pressure resulted in a black hole, and viceversa. The condition to be satisfied by the tangential stress is .

To see this more explicitly, in analogy with the previous example, consider for instance


then the tangential stress near the center of the cloud will be given by,


We note that the tangential pressure as well as the density diverge at the singularity. The collapse here ends in a naked singularity whenever


as can easily be calculated from equation (IV.15).

Figure 1: The Naked singularity region and the black hole region in the parameter space for fixed values of , and .

In Fig. 1, the occurrence of black hole and naked singularity phases in the latest example are shown, depending on the possible values chosen in the parameters space of and , for fixed values of , and (taken at the initial time). The line corresponds to , and acts as a critical surface separating the two outcomes, in which case the evaluation of higher order is necessary in order to determine the final fate of collapse. The axis corresponds to the LTB inhomogeneous dust case, with the origin being the Oppenheimer-Snyder-Datt model.

We see clearly that to have a naked singularity in the dust case, we must require . This is not the case anymore in the more general situation with non-vanishing tangential stresses. In fact, such scenarios allow for positive values of to lead to a naked singularity, provided that the tangential stress (and hence in this example) is negative, as well as negative values of leading to black holes whenever a suitable positive tangential pressure is introduced. The interesting feature here is that the introduction of tangential pressures does not rule out the possibility of formation of naked singularity.

This shows that the conjecture suggested in Magli () does not hold (even in the more restrictive case of quadratic density and pressure profiles), and that naked singularities are a general feature of gravitational collapse with pressures. The conjecture stated that given a known Lemaitre-Tolman-Bondi collapse model ending in a black hole, all quadratic pressure perturbations to the model will necessarily originate a black hole. Now we see that adding a suitable generic pressure perturbation to the LTB model can in fact change the outcome of collapse.

¿From a physical point of view, the most relevant quadrant in the model shown in Fig. 1 is the first one (where and ), since it corresponds to positive pressures and radially decreasing energy density profile. We see that both the black hole and naked singularity outcomes are possible in this case.

We shall stress however that the model illustrated above, while presenting a simple and intriguing structure, does not depict the only possible structure of collapse with tangential pressures. Other examples, with different choices of the functions , and , might provide an entirely different picture for the final stages of collapse, thus indicating once more how rich and complex is the description of collapse resulting from general relativistic analysis.

This model shows also how the black hole formation process described in the well known homogeneous dust scenario is not ‘stable’ when arbitrarily small stress perturbations are introduced. Since, as shown in Fig. 1, it lies on the critical surface separating the ‘naked singularity configuration space’ from the ‘black hole configuration space’, any close neighborhood of the OSD model will contain initial data leading to both outcomes.

On the other hand, we see that in the context of this model, inhomogeneous dust collapse leading either to black hole or naked singularity (thus away from the critical surface where ), is ‘stable’ under small pressure perturbations, in the sense that it is always possible to find a small neighborhood of the configuration space which will develop the same final outcome.

The above considerations illustrate the role pressures can play towards determining the final fate of collapse. Clearly, depending on their nature, the non-zero tangential pressures can help create a naked singularity, or a black hole. It follows that a mere introduction of pressure in the collapsing cloud, by itself, cannot help us restore the cosmic censorship. Finally we note that the structure illustrated for the above model extends naturally to more general pressure profiles, still leaving unchanged the general behaviour where both naked singularities and black holes might generate at the end of collapse.

V Genericity of black hole and naked singularity final states

The consideration above helps us gain an important insight into the nature and genericity of the black hole and naked singularity formation in gravitational collapse. Let us take, for example, a scenario where the term is positive, thus leading the collapse to a naked singularity final fate. In that case, it is clear that there always exists a small enough perturbation in the initial data, which is in the form of small changes in the initial density, pressures, and velocities of the collapsing shells, which will still preserve by continuity the positive sign for the term, thus preserving the naked singularity final fate of the collapse. A similar conclusion will hold for the black hole final state of collapse as well.

In this sense, if a collapse evolution was to develop in a naked singularity, there is a small neighborhood around that specific initial set in the configuration space of initial data, which also continues to take the collapse to a naked singularity. In such a sense, the naked singularity formation (as well as the black hole formation) is generic in the gravitational collapse with non-zero tangential pressures case.

In the light of these results, we can see that the idea that the introduction of pressures must favor the black hole formation process is not correct, at least for non-vanishing tangential stresses. One may argue that this situation arises because we have neglected radial pressures, a choice that, while still not being the dust model, is not sufficiently ‘realistic’. Nevertheless the same arguments can be extended to the case of generic pressures (as it will be done in a future work), to show that collapse scenarios leading to a naked singularity have a small neighborhood in the configuration space of initial pressure profiles, that still leads to a naked singularity.

Since the parameters that determine the final fate of collapse cannot be chosen arbitrarily if one wishes to obtain a black hole, we can say that a specific final outcome of collapse requires some tuning in the initial data in contrast with the possible view that any initial configuration with pressures will necessary destroy the naked singularity formation picture and lead to a black hole.

The fact that clearly emerges from the present analysis is that the scenarios describing the formation of naked singularities as the endstate of collapse are not mere counterexamples to cosmic censorship whose relevance is limited and confined. We can therefore say that the presence of tangential stresses can ‘undress’ a singularity that was covered in the absence of stresses, or viceversa, a naked singularity dust model can be covered by the introduction of a suitable positive pressure.

It is therefore useful to relate the scenario presented here with the known collapse models without pressures, namely the dust collapse scenarios. The first collapse model within the framework of general relativity was studied by Oppenheimer, Snyder and Datt, and it described the collapse of a massive sphere made of homogeneous dust. The Lemaitre-Tolman-Bondi models then considered the collapse where inhomogeneous dust is allowed. It turned out that when inhomogeneities in the density distributions are allowed, then the appearance of naked singularities takes place as possible endstate of gravitational collapse. Dust collapse models are well known and widely studied and the structure of the formation of trapped surfaces and of the singularity is well established, also in the case of higher spacetime dimensions n-dim-dust ().

The dust collapse is obtained from the above framework, when we impose . In that case, from equation (II.5) it follows that , which together with the condition gives . This in turns implies that and . For inhomogeneous dust in general we have,


We see that inhomogeneous dust can lead to the formation of naked singularities or black holes, depending on the behaviour of the mass profile and the velocity profile JDM (). Furthermore, in the marginally bound case it is possible to evaluate explicitly the apparent horizon from equation (III.11) as


while in the bound and unbound cases, the apparent horizon curve can be read from equation (III.11).

Collapse is said to be homogeneous when . Homogeneity, together with the initial condition , implies some constraints on the possible matter models and velocity profiles. Now let us consider the following three conditions,

  • ,

  • ,

  • ,

where the second condition implies a simultaneous collapse, in which all shells of matter fall in the singularity at the same time. We can see that homogeneity, dust and the above conditions are not all mutually independent. In fact, in the case of collapse with vanishing radial stresses, which is under consideration here, we can prove the following statements:

  1. Homogeneity (i) and (ii),

  2. (ii) (i), (iii) and dust.

We can see from this that homogeneous collapse implies necessarily dust, while vice versa is not necessarily true (note, however, that when radial stresses are included it is possible to have homogeneous collapse also in the case of a perfect fluid with homogeneous pressures). Therefore, the Oppenheimer-Snyder-Datt model implies necessarily that the three conditions above are satisfied. We note that the simultaneous collapse (which implies homogeneous collapse in the case when ) implies that all shells fall into the singularity at the same time , thus excluding the possibility of its visibility due to the structure of the apparent horizon in that case. The quantity vanishes identically in this case for every . In order for the singularity arising in dust collapse scenarios to be visible, inhomogeneities in the density distribution must be allowed for and have to be taken into consideration, thus dropping the condition (i), which automatically implies the non-occurrence of a simultaneous collapse.

To see the first statement, it is sufficient to use equation (II.4) with the requirement of homogeneity. From that it follows that , which together with the requirement that implies (i) and (ii). To see the inverse implication, it is sufficient to put (i) and (ii) into equation (II.4). To prove the second statement, we make use of equation (II.7) written in the form,


Since does not depend on , it has the same value for any , which is also that obtained for , namely,


¿From this it follows immediately,


Vice versa, assuming dust, and conditions (i) and (iii) to hold, we see that does not depend on , which implies . That together with the requirement implies .

Vi Gravitational collapse with cosmological constant

It is interesting and indeed important to see how the presence of a non-zero cosmological constant can affect the dynamical evolution of a gravitationally collapsing massive matter cloud. The interest in models with cosmological constant arises primarily from the observation of the accelerating speed in the expansion of the universe, an observed phenomenon that indicates that the universe might be filled with some repulsive force known as ‘dark energy’. One of the possible ways by which dark energy can be described is through the introduction of a non-zero cosmological term in Einstein equations.

The earlier works on the subject were constrained by some simplifying assumptions. Namely the collapsing matter was taken to be a dust cloud or the function was chosen to have a specific dependence on and in order to simplify the integration of equation (II.13) (see Lambda ()).

In the present section we extend the previous results to the most general case and test how the dynamics and the final stages of collapse of the models described above are affected by the introduction of the cosmological constant. In doing so, we consider an example of dynamical collapse in the presence of tangential stress and cosmological term, that connects directly with the models described in section IV.

The reason for the interest in the introduction of only tangential stresses in the collapsing cloud with cosmological constant resides in the fact that, from equation (II.19), we see that the matching with an exterior vacuum background must be done with the well-known Schwarzschild-deSitter or anti-deSitter spacetimes. On the other hand, the presence of radial pressures would require the matching with a more complicated generalized Vaidya metric.

Einstein equations (II.3) and (II.4), in the presence of a cosmological constant term take the form,


while the other equations remain unchanged. The first one above leads to a Misner-Sharp mass of the form


while the energy density remains unchanged. The energy conditions here imply


however, all of the above need not be respected now, and we may have a weaker form of energy conditions holding, depending on what the sign of the cosmological constant is chosen to be. The matching with the exterior metric in this case needs to be done with a spacetime that is asymptotically either Schwarzschild-deSitter or Schwarzschild-anti-deSitter, depending on the sign of the term.

The integration of Einstein equations proceeds exactly in the same way as described above, and we therefore obtain,


where the quantity can now be interpreted as an effective potential. The regions of allowed motion are those for which (since ). The zeros of correspond to turning points of the dynamics where the collapse is halted and the shell bounces back. Therefore, if for some , the cloud will bounce back and not collapse indefinitely towards the center. Of course, this situation is not allowed in the dust cases when , is small and where gravity is the only force acting on the collapsing particles. Still, bouncing behaviours were found to be possible in scenarios with tangential stresses and bound velocity profiles. In the same manner, whenever , it is possible to observe such bouncing behaviours in the dynamics of the collapsing shells.

Typically, such a scenario will become relevant for collapsing clouds of rather large sizes in the universe, for example, for the case of a large cluster of galaxies collapsing under its own gravitational attraction, or possibly for a large collapsing star cluster. In any case, for the sizes of a collapsing star of an average size, the effects of a non-zero may not make any significant impact on the final outcome of its gravitational collapse.

Following a similar analysis of the time curves as that developed in section III, we can evaluate in this case as,


This leads the numerator of (which is the essential element in determining its sign) to have exactly the same value as in equation (IV.4):


We can thus see that the presence of the term affects collapse in the sense that it might introduce some positive roots in the effective potentials, and therefore would give rise to some bouncing behaviours. But it will not affect the terms that are responsible for the sign of , which determines the final outcome of collapse to be a black hole or a naked singularity. This appears reasonable, when we think of the presence of the cosmological constant as related to a dark energy in the universe, because the effect that the dark energy will have on collapse can be that of disrupting and preventing collapse from happening by causing a bounce. But it cannot influence the nature of the final fate of a collapse that is uniquely determined by the initial configurations for densities, pressures and velocities.

Furthermore, the presence of an added positive -term can lead to while , thus in fact allowing the evolution for some configurations that were forbidden in the case of vanishing .

The visibility of the singularity in the case of collapse with a cosmological constant differs from the case studied above also in the structure and formation of the apparent horizon. From equation (VI.3) and equation (III.9) we get


from which we see that in general the apparent horizon equation is a cubic equation in that can have zero, one or more than one positive roots, thus affecting considerably the formation of trapped surfaces. Still it is possible to show (see Lambda ()), that the presence of the cosmological constant does not prevent the possibility of formation of naked singularities, which continues to be a general feature of gravitational collapse.

It is interesting at this point to analyze some examples of dynamical behaviours based on the above considerations. We will thus make a specific, although reasonable, choice of the mass, pressure and velocity profiles and confront the effective potential given by equation (VI.5) written in the form,


for different values of .

Following the examples provided in section IV we will choose only quadratic terms in the mass and stress profiles and study a ‘quasi-Hookean’ equation for the pressure of the form