Self Similar Collapse and the Raychaudhuri equation

The role of the Raychaudhuri equation in studying gravitational collapse is discussed. A self-similar distribution of a scalar field along with an imperfect fluid in a conformally flat spacetime is considered for the purpose. The general focusing condition is found out and verified against the available exact solutions. The connection between the Raychaudhuri equation and the critical phenomena is also explored.


Introduction
It is widely accepted that at the end of its life-cycle, a massive astronomical body undergoes a phase of gravitational collapse. Followed by the pioneering work on an idealized model star collapsing to zero volume by Datt [1] and Oppenheimer and Snyder [2], many attempts have been made to find a more complete and increasingly generalized description of gravitational collapse. For a detailed review of the relevance of gravitational collapse in physics and the open problems, we refer to the references [3,4].
It is not a trivial task to simplify the Einstein field equations enough to guarantee the extraction of an exact solution that can describe the collapsing evolution. This often compels one to resort to the study of simplified systems which can describe the essential physics. Christodoulou studied the global initial value problem for Einstein's equations in the spherically symmetric case with a massless scalar field and a e-mail: sgc14ip003@iiserkol.ac.in b e-mail: soumya.chakrabarti@saha.ac.in c e-mail: adg@iiserkol.ac.in d e-mail: narayan@iiserkol.ac.in showed that the scalar field, depending on the initial collapsing profile, can either converge towards a zero proper volume or bounce, dispersing towards infinity [5][6][7]. Choptuik studied a similar problem numerically and pointed out the now well known critical phenomena in massless scalar field collapse [8]. Similar critical phenomena were studied by Brady, Chambers and Goncalves [9] and Gundlach [10]. Throughout the last decade, many attempts have been made to generalize the ideas of critical phenomena for less simplified scenarios. For example, critical phenomena in a collapsing system with angular momentum was studied by Olabarrieta, Ventrella, Choptuik and Unruh [11]. The critical behavior of a spherically symmetric collision-less matter at the threshold of black hole formation was studied by Olabarrieta and Choptuik [12]. Analytical investigations of the critical collapse problem, and a search for a theoretical explanation for the behavior discovered by Choptuik were carried out by Brady [13], under the assumption that the collapse is selfsimilar. For a detailed review on the critical phenomena in Gravitational Collapse we refer to the review by Gundlach and Martin-Garcia [14].
The current work stems from a different motivation -we intend to investigate the role of the Raychaudhuri equation [15,16], in the context of the study of gravitational collapse and critical phenomena. This equation has been used to study relativistic charged collapse by Kouretsis and Tsagas [17]. In the present work we study the evolution of a conformally flat geometry, minimally coupled with a scalar field, along with the presence of a fluid. This specific case of spacetime metric is well-studied for radiating and/or shear-free stars [18][19][20][21][22][23][24][25] and is receiving increasing interest in the context of gravitational collapse quite recently [26][27][28]. We make the additional assumption that the evolution is self-similar in nature. Gravitational collapse of a self-similar fluid with heat flow in a conformally flat spacetime was studied by arXiv:1912.05824v1 [gr-qc] 12 Dec 2019 Chan, Silva and Rocha [29]. Collapsing models of fluids with the assumption of self-similarity in more general spacetimes were studied in [30][31][32][33] and references therein. In this work, we first discuss the usefulness of the Raychaudhuri equation in studying the dynamics of the system. Generic conditions regarding the evolution of the spacetime can be provided using the Raychaudhuri equation. These conditions are useful in the absence of an exact solution. We then work out exact solutions of field equations for some special cases making use of some simplifying assumptions. We show the consistency of the conclusions regarding the dynamics from exact solutions, if available, with those arrived at from the Raychaudhuri equation. The role of the Raychaudhuri equation in the context of the critical phenomena is also studied.
The paper is organized as follows. In section 2 we introduce our system and discuss the assumptions involved. The field equations and conservation equations for the system are also presented in this section. The Raychaudhuri equation and focusing condition are briefly discussed in the section 3. In this section, we focus on the role of the Raychaudhuri equation in extracting information about the evolution of spacetime. This section also includes a discussion on the consistency of the results with exact solutions of the field equations, where available, and a connection between focusing condition and the critical phenomena. The final section 4 includes some concluding remarks.

The System
The spacetime which we consider is a conformally flat spherically symmetric spacetime for which the metric can be written as, where 1 A 2 (r,t) is the conformal factor which governs the evolution of the 2-sphere.
The contribution to the energy-momentum tensor comes from a scalar field and a fluid, where and T fluid µν = (ρ + p t ) u µ u ν + p t g µν + (p r − p t ) χ µ χ ν + q u µ χ ν + u ν χ µ .
ρ, p t , p r and q are the energy density, tangential pressure, radial pressure, and radial heat flux of the fluid respectively; u µ = Aδ µ 0 is the velocity of the fluid and χ µ = Aδ µ 1 is a unit spacelike vector along the radial direction. We assume the system to be self-similar of the first kind in nature i.e. the metric admits a homothetic Killing vector [34][35][36]. In the present case, we write, where z = t r . With this choice the field equations become self similar (i.e. the only independent variable in these equations is z).

Field Equations
The Einstein field equations (in the units 8πG = 1) for the metric (1) with energy momentum tensor (2) can be written as, and where dot and prime denote differentiation with respect to t and r respectively. In terms of the self-similar variable z, equation (9) can be written as, Equation (10) suggests that if φ is constant and q = 0, B will be proportional to z. Similarly, we can write down equations (6)- (8) in terms of z as, where equation (10) has been used to replace terms containing derivatives of the scalar field. Using equations (11) and (12), one can easily show that, while (12) and (13) will give,

Conservation Equations
We shall now write down the conservation equations for this system under the assumption that the energy momentum tensors corresponding to the fluid and the scalar field are conserved independently. The conservation equation for the scalar field yields the wave equation, For the present metric (1), this equation (16) translates into, In terms of the self-similar variable z, we have, The conservation equations for the fluid is given by, from which we will get two non-trivial equations, and For the self similar case, if we use equation (14), both equations (20) and (21) yield the same equation as,

Focusing Condition
We have chosen a comoving observer, so that u α = Aδ α 0 . As the metric is conformally flat, the shear term and rotation term will vanish. Now, we know that the congruence will focus [37] within a finite affine parameter value if, which, with equation (23), leads to, Here the left hand side of the expression is related to the matter sector via the Einstein equations and the right hand side is the divergence of acceleration. Focusing may be forbidden when the divergence of acceleration opposes the gravitational attraction and its contribution dominates over that due to the matter part. If the divergence of acceleration term is sufficiently high, the evolution may lead to a complete dispersal.
For the metric (1), the condition (25) in terms of the conformal factor A(r,t) is given by, For the self similar case, the above condition (26) becomes,

Raychaudhuri Equation and the dynamics of spacetime
It is difficult to find exact solutions of the field equations for our system without any simplifying assumptions. The Raychaudhuri equation can be applied to provide generic conditions, regarding the dynamics of the spacetime. These conditions may lead us to some useful information about the evolution. We will now try to find such conditions. We have the condition for focusing (27) from the Raychaudhuri equation as, Using (15) this condition can be written as, These generic conditions like (28) or (29), may appear useful. These are generic conditions if satisfied by the conformal factor or the energy momentum tensor components, will lead to the formation of a singularity. Thus, there is an important role of these conditions in dictating the dynamics of the spacetime. Let us illustrate this point with a simple but very important example. If the fluid is a perfect isotropic fluid having no radial heat flux, the left hand side of the condition (29) is equal to zero. Thus, the condition is satisfied only when z → ∞ which means that there will be a central singularity at r = 0 or a singularity forms as t → ∞ which is inconsequential.
Thus, the possibility of the formation of a singularity at a finite future or the avoidance of a central singularity will depend on pressure anisotropy or heat flux.
We will discuss a few more examples. The condition (28) or (29) can be recast into different forms using the field equations. For example, if we use equation (10) and (15), the condition (28) can be written as, Now, let us consider the case where the scalar field is absent or a constant (equivalent to a cosmological constant) and q = 0. In this case the above condition (30) will be satisfied only when which implies dB dz = 0 (using equation (15)) and there will not be any evolution of the spacetime. Therefore, we can conclude that formation of a singularity in this case can be avoided. However, q = 0 or dφ dz = 0 or both may lead to the formation of a singularity. If we have a large rate of change of the scalar field, compared to the other terms, present in the condition (30), singularity formation is inevitable.
Let us discuss another example where dφ dz = 0 and p r = p t . The condition (30) then translates into, which gives, Thus, there will be singularity formation only when the heat flux of the fluid satisfies this constraint (33).

Exact solutions and the Raychaudhuri equation
We have not used exact solutions of the field equations in our discussion so far. It is worthwhile to check the consistency of the conclusions arrived at using the Raychaudhuri equation with those using exact solutions whenever available. For many special cases, i.e., with various sources, exact solutions for the metric can be found out. Some of them are listed in the Appendix A. With an exact solution, one can explicitly find whether there is a collapse or an expansion and whether the collapse, if there is any, results in a singularity. We have checked that in all such cases the conclusions are consistent with those obtained from the focusing condition (27) found out from the Raychaudhuri equation. One such nontrivial example is described in detail below.
We will discuss the case where a scalar field and a fluid both are present and we make an assumption that p r = p t but q = 0. For details, see Appendix A.5. The solution for B and A in this case are respectively given by, and where β is a constant given by p t = β p r . For β > 0, there may be a zero proper volume singularity when t → ∞ and/or r → ∞, which can be excluded from the discussion. There will be a dispersal in this case at z = k F 1 2β . When β < 0, formation of a singularity at z = k F Using the solution for B (equation (34)), one can show, , the first term in the right hand side of the above expression dominates. Only when β < 0, the condition (28), obtained using the Raychaudhuri equation, is fulfilled in this region and a singularity within finite z is unavoidable. On the other hand if β > 0, a dispersal will occur in this region. Singularity formation as z → ∞ is a possibility in this case. For β < 0, it is easy to show that the right hand side of the equation (36) must be negative as z → 0 which makes a dispersal possible here.

Critical Phenomena and the Raychaudhuri equation
From the exact solutions in Appendix A, we have found that collapse or dispersal correspond to the the quantity d dz 1 B dB dz being positive or negative, respectively. Therefore, a transition from collapse to dispersal or vice versa will be accompanied with a change of the sign of this quantity. Thus, a relation between focusing condition and the critical phenomena seems to be indicated. We will illustrate this point using the example discussed in the previous section.
In this example, let us take the z → k F 1 2β into consideration. When β is negative a singularity forms at z = k F . This shows that the critical parameter β determines the signature of d dz 1 B dB dz which in turn determines the focusing (collapse) or dispersal of the spherically symmetric distribution and the critical value of β is zero.

Conclusion
We have found the focusing condition (28) for a self similar matter distribution in a conformally flat spherically symmetric spacetime. Although the spacetime, we have considered, has a stringent symmetry requirement, the matter distribution is quite general and includes anisotropic fluid pressure, heat flux, as well as a minimally coupled scalar field. The nature of the conclusions drawn from this condition is verified against quite a few exact solutions that are available, and one of them is discussed in detail.
Thanks to the investigations as in references [8][9][10] and others, it is now well known that a scalar field collapse may have a critical phenomenon associated with it. The condition developed from the Raychaudhuri equation helps visualizing the critical phenomena in general. This remains valid even when the matter distribution includes a fluid. We have discussed one example, where exact solution is available, in detail. In this example, we could see the existence of a critical parameter, and could also determine the critical value (β = 0). It is intriguing to note that the parameter has nothing to do with the scalar field in this case and is determined only by the fluid, as β is simply a parameter connecting the radial and transverse fluid pressure. With the help of this example we have explicitly discussed the relation between the focusing condition and the critical phenomena.
which is inconsistent for finite z values. Thus, a consistent solution of the field equations with a massless scalar field as matter source is not possible under the mentioned assumptions.
Appendix A.2: Scalar field with a non-zero potential Even if we include a potential in the energy momentum tensor of the scalar field, we have the same solutions for B and A as in the previous case (equations (A.2) and (A.3) respectively). Here with equation (A.2), equations (11), (12) and (13) all yield, Therefore, the potential must remain constant. From equation (10) we have, which is consistent with the equation (18). Thus, the net energy momentum tensor effectively behaves like that of a cosmological constant. From equation (A.3), we can conclude that there is a singularity only when t approaches infinity where the scale factor (inverse of A) becomes zero which signifies an ever collapsing solution. At t = 0, A = 0 i.e. the scale factor becomes infinite which signifies a dispersal. Appendix A.3: Scalar field along with a perfect fluid If we include perfect isotropic fluid (p r = p t and q = 0) along with the scalar field, the solution for B will not change as dictated by equation (15). From equation (14) we have, ρ, p, φ and V (φ ) will remain constant in this case. Thus the matter turns out to be the cosmological constant for a consistent solution. Clearly, the conclusions in this case remain the same as in the previous case.
In the examples, discussed so far we have not observed any possibilities of singularity formation within finite time and the matter turns out to be the cosmological constant.
Appendix A.4: Fluid with isotropic pressure and radial heat flux This particular case, with p r = p t and q = 0, has already been studied by Chan, Silva and Rocha [29]. They found that there is a singularity formation due to collapse as t → 0. In this case from equations (15) and (10)  where C and D are constants of integration. From the above equation (A.9) we have, The expression for A confirms that there is a zero proper volume singularity at t = 0. The scale factor also becomes zero when t → ∞ and/or r → ∞. These possibilities are not worth considering. There may be a dispersal at z 2 = t 2 r 2 = − 2D C only if C and D are of opposite signs.
Appendix A.5: Scalar field along with an anisotropic fluid In this case, we have p r = p t but q = 0. Now, from (14) and (15)  The same equation of state as in (A.11) was obtained by Brandt et al [32] from the requirement of self-similarity even with a relaxation of the conformally flat condition. In this case, the conservation equation for the fluid (equation (22)) takes the form, It is difficult to find a solution for the conformal factor in this case without further simplifications. Thus, we will assume that the tangential pressure is proportional to radial pressure i.e. p t = β p r . Brandt et al worked with a similar assumption for a spherically symmetric spacetime which is not conformally flat [31]. With this assumption, solving equation (A.13) we have, where n is constant of integration. If we replace this in equation (A.12), the solution for B comes out as, where F is constant of integration and k = n(1 − β ) 2 . Consequently, the solution for A is, (A. 16)