Abstract
We model the gravitational collapse of heavy massive shells including its main quantum corrections. Among these corrections, quantum improvements coming from Quantum Einstein Gravity are taken into account, which provides us with an effective quantum spacetime. Likewise, we consider dynamical Hawking radiation by modeling its backreaction once the horizons have been generated. Our results point towards a picture of gravitational collapse in which the collapsing shell reaches a minimum nonzero radius (whose value depends on the shell initial conditions) with its mass only slightly reduced. Then, there is always a rebound after which most (or all) of the mass evaporates in the form of Hawking radiation. Since the mass never concentrates in a single point, no singularity appears.
Singularity free gravitational collapse in an effective dynamical quantum spacetime
R. Torres
Department of Applied Physics, UPC, Barcelona, Spain.
[.1cm]
KEYWORDS: Gravitational Collapse, Black Holes, Hawking radiation, Quantum Gravity.
1 Introduction
It is expected that a large enough object would collapse classically until a horizon forms. Then Hawking radiation would appear and the mass of the object should be reduced. Less is known on the details of the later evolution. In fact, a complete investigation of the process would require a complete consistent theory of quantum gravity together with the calculational tools to achieve a description of the scenario. Since such apparatus is not currently available, for the moment one can only resort to the study of toy models in which the known main quantum contributions are taken into account. By this means, one can try to probe some of the features that one could expect from a full theory of quantum gravity.
In this letter we will work in this direction. In our toy model two main simplifications will be carried out. First, we assume the existence of a spherically symmetric spacetime in which the collapse takes place. Second, we choose as our collapsing object a thin shell. In other words, we assume that the spacetime is split in two different regions with a common spherically symmetric timelike boundary corresponding to the thin shell.
This second simplification deserves some comments. Clearly it means that we would be able to probe gravitational collapse only whenever the approximation in which one can neglect the shell thickness remains valid. An investigation of the conditions under which this is possible was carried out in [1] (see also [2]). The authors considered a shell composed of a number of (swave) scalar particles with mass bound together by gravitational interaction. The particles form a radially localized bound state corresponding to a finite thickness shell, whose mean position approximately follows a classical collapse. Moreover, during the collapse the average shell thickness decreases according to [1]
where is Newton’s gravitational constant. On the other hand, the fluctuations associated with the quantum nature of matter become dominant for a radius of the shell of the order of the Compton wavelength of the constituent quanta. So that they are negligible as long as [1]
In other words, if we want to probe the last stages of the collapsing phase by using the thinshell approximation we should use a shell composed of a high number of very heavy particles (large ’s) and, thus, we would be using a heavy massive shell. Only under these conditions the results obtained using the thinshell approximation are likely to be similar to the results that one would obtain in a real collapsing situation
With regard to the shell exterior region , we will describe it with a portion of an improved Schwarzschild solution with mass equal to the shell mass. Specifically, we choose for the exterior region an effective improved solution coming from Quantum Einstein Gravity that incorporates quantum corrections to the classical solution. It does so by taking into account the effect of virtual gravitons. I.e., just as in quantum electrodynamics the virtual pairs imply the existence of a screening effect leading to a running electric charge, when one considers the existence of virtual gravitons one obtains an antiscreening effect leading to a running gravitational constant, which is used to get the improved Schwarzschild solution ([5] and references therein). A summary of this effective solution will be carried out in section 2. On the other hand and following this approach, the shell massless interior region will have to be described by a portion of Minkowski’s spacetime (what is equivalent to a massless improved Schwarzschild solution).
Since the improved exterior solution possesses horizons, the tunneling of virtual particles through them is expected to produce Hawking radiation. Thus, in section 3 we will also model the effect of the backreaction to the radiation in the effective solution. Then, in section 4 we will consider the matching of the interior and exterior solutions through the spherically symmetric thin shell by using Israel’s formalism [6]. This will provide us with the shell evolution equation which will allow us to analyze its different attractive and repulsive contributions. Finally in section 5 the numerical integration of the evolution equation will be carried out and the results will be interpreted.
2 Exterior: Improved Schwarzschild solution
As explained in the introduction, in order to model a collapsing shell we should first establish the exterior to that shell (). In this work we want the exterior to incorporate the main quantum corrections to the classical solution. This can be done by using a renormalization group improved Schwarzschild solution found by Bonanno and Reuter [5] that can be written as
(2.1) 
where
(2.2) 
is Newton’s universal gravitational constant, is the mass measured by an observer at infinity and and are constants coming from the nonperturbative renormalization group theory and from an appropriate “cutoff identification”, respectively. The qualitative properties of this solution are fairly insensitive to the precise value of . In this way, in order to simplify the calculations, it is usual to choose [5][7]. On the other hand, can be found by comparison with the standard perturbative quantization of Einstein’s gravity (see [8] and references therein). It can be deduced that its precise value is , but again the properties of the solution do not rely on its precise value as long as it is strictly positive.
If we define
the horizons of the improved solution can be found by solving . Then, it is easy to see that the horizons correspond to the number of positive real solutions of a cubic equation and depend on the sign of its discriminant or, equivalently, on whether the mass is bigger, equal or smaller than a critical value . In particular, the value implies
where is the Planck mass. If then the equation has two positive real solutions satisfying . The existence of an inner solution represents a novelty with regard to the classical spacetime. However, it is interesting to remark that it is a result common to different approaches to Quantum Gravity. (See, for example, [5][9][10][11]). The outer solution can be considered as the improved Schwarzschild horizon, i.e., the Schwarzschild horizon with quantum corrections taken into account. The ‘improvement’ in this horizon can be made apparent for masses much bigger than Planck’s mass if one expands in terms of obtaining
The global structure of the improved solution for resembles the global structure of the ReissnerNordström spacetime with mass bigger than its charge (). A Penrose diagram corresponding to the improved solution for the case is shown in fig. 1.
On the other hand, if then there is only one positive real solution to the cubic equation and the global structure resembles that of a extremal ReissnerNordström solution (), whereas if the equation has not positive real solutions.
3 Hawking radiation from the horizons
We will now summarize the results on Hawking radiation in the quantum improved solution. A more complete description can be found in [12][13] which, in turn, are based on the tunneling approach by Parikh and Wilczek [14]. We consider Hawking radiation coming out from an improved black hole satisfying thanks to the tunneling process occurring both in the outer and in the inner horizons.
Since the coefficients of the metric do not depend on there is a killing vector which is straightforwardly found to be timelike for , lightlike for , spacelike for , lightlike for and timelike for .
The possibility of tunneling is based on the fact that the killing vector is spacelike for (region II in fig.1), what allows the existence of negative energy states.
Let us consider that a pair of photons is created in region II, where the 2spheres are closed trapped surfaces.
A pair of test photons would be classically forced to move inwards until reaching . However, for nontest photons energy conservation modifies this picture near the outer (inner) horizon since the positive energy photon produced in the pair could ‘tunnel’ the outer (inner) horizon and move outwards in region I (III ,respectively
In order to compute the tunneling rate we will rewrite the improved Schwarzschild’s solution in Painlevélike coordinates [15] so as to have coordinates which are not singular at the horizons. In order to do this it suffices to introduce a new coordinate replacing the Schwarzschildlike time such that and fix by demanding the constant time slices to be flat. In this way one gets:
(3.1) 
where can now take the values .
Let us consider pair production occurring close to one of the horizons with the positive energy particle tunneling from region II outwards. The standard results of the WKB method for the tunneling through a potential barrier that would be classically forbidden can be directly applied due to the infinite redshift near the horizon [14]. In particular, the semiclassical emission rate will be given by , where is the particle action. Therefore, we have to compute the imaginary part of the action for an outgoing positive energy particle which crosses one of the horizons outwards from to ,
(3.2) 
where we have used Hamilton’s equation , the equation for null geodesics and the fact that the BH loses mass after the emission of a shell (i.e., ) and, thus, becomes , which stands for with replaced by (see [12][13] for details).
If we define
and
where is the position of the outer (‘+’) or inner (‘’) horizons when is replaced by , then, by deforming the contour of integration so as to ensure that positive energy solutions decay in time, one can then write (3.2) as
(3.3) 
where the subindex ’’ or ’’ corresponds to the tunneling through the outer or inner horizon, respectively.
The semiclassical rate through every horizon will be
(3.4) 
When quadratic terms are neglected we can develop Im up to first order in as
obtaining a thermal radiation for the quantum black hole () with (positive) temperature at every horizon
(3.5) 
On the other hand, if we consider the full consequences of energy conservation, the distribution function for the emission of photons can be written as (see [16] –correcting the result of [17]–)
What for our quantum corrected solution becomes, at every horizon,
(3.6) 
with the additional requirement that, according to the properties of (see [12][13]), the energy of the emitted particles must satisfy .
Using (3.6), the flow of positive energy due to the tunneled particles at the horizons can be written approximately as [18][19]
(3.7) 
where we have taken into account the possible backscattered region in the greybody factor (more details in [12][13]) and we are compelled to take into account in the integration limits that the maximum energy of a radiated particle could be .
3.1 Modeling the backreaction
In order to modelize the evaporation, let us first write the improved Schwarzschild’s metric (2.1) in terms of ingoing EddingtonFinkelsteinlike coordinates , where
as
(3.8) 
This solution does not reflect the backreaction associated to the lost of mass due to the tunneling effect. However, we can model the mass lost taking into account that, whenever a pair of virtual particles is created, when the particle with positive energy escapes to infinity its companion, with negative energy, falls inwards and reduces the total mass. In this way, if we consider negative energy massless particles following ingoing null geodesics constant, the mass becomes a decreasing function . The metric which incorporates the effect of the decreasing mass due to the ingoing null radiation is (3.8) with replaced by , i.e., it corresponds to an improved ingoing Vaidya solution [7] that for the case that we will treat in our specific model takes the form
(3.9) 
where, from (2.2), the running gravitational constant takes the form
(3.10) 
Let us comment that when one considers the backreaction the solutions to become dynamical (marginally trapped) horizons .
On the other hand, the flux of negative energy particles directed inwards equals the flux of outgoing radiated particles and, therefore
(3.11) 
where comes from the combination of fluxes and coming from the outer and inner horizons, respectively.
4 Collapsing model
In order to model the collapsing thin shell we will now use Israel’s formalism [6][20]. We, therefore, assume that the spherically symmetric spacetime is split in two different regions with a common spherically symmetric timelike boundary : The thin shell. We choose that the shell is described by coordinates such that is the shell proper time and the parametric equations from can be locally written in EddingtonFinkelsteinlike coordinates as (and similarly for ). In order to have a welldefined geometry at the first fundamental forms (or induced line elements) of the boundary
where , must agree when computed from or , i.e., . A first important consequence of this fact in the case of spherically symmetric scenarios is that the areal coordinates () for the interior and exterior regions must agree on the shell. In this way, if we want to describe the evolution of the shell, the same function can be used from either the point of view of or .
Let be the unit spacelike vector () pointing from to . If we now define as a coordinate such that with at the hypersurface , the energymomentum tensor of the spacetime would have the form
(4.1) 
where is the Heaviside stepfunction and , which is tangent to , is the energymomentum tensor of the hypersurface. The extrinsic curvature or second fundamental form of is defined as
and it is related to the energymomentum tensor of the hypersurface through the Lanczos equations
(4.2) 
where and we are using the shell energymomentum threetensor and .
In our model we will consider that is composed of radially moving noninteracting particles, so that
(4.3) 
where is the massenergy density of the layer and is the shell 4velocity. As seen from the exterior region we will have
(4.4) 
where the dot stands for derivative with respect to . On the other hand, using the metric (3.9), the normalization condition implies
(4.5) 
where and . From this last definition we see that a necessary condition for to be timelike is
(4.6) 
In physical terms this can be interpreted as requiring the collapsing shell to have a minimum speed while it traverses the region where the 2spheres are closed trapped surfaces (). On the other hand, the extrinsic curvature of when computed from the exterior region using (3.9) is
(4.7)  
(4.8)  
(4.9) 
Since the interior of the shell has no mass, one can get the results for the interior region from the previous ones with . In particular, (3.9) becomes Minkowski’s solution and
(4.10)  
(4.11)  
(4.12)  
(4.13) 
4.1 The evolution of the shell
The Lanczos equations (4.2) provide us with
(4.14)  
(4.15) 
which can be easily combined using (4.7, 4.8, 4.11, 4.12) in order to obtain an equation for the massenergy density
(4.16) 
and an evolution equation
(4.17) 
In order to analyze the meaning of this equation, first note that , and that . Then we can differentiate three force terms in the right hand side of the equation. First, the usual term coming from the selfgravitating shell
that satisfies which implies that the shell, in the regions where the quantum corrections are small, is accelerated inwards. Second, the Hawkingradiation term
that satisfies . Specifically, in this approach this term does not contribute previous to the formation of a horizon since there is not Hawking radiation. However, once a horizon appears one has due to the negative energy coming from the horizon which is absorbed by the shell. Clearly, implies that the absorbtion of negative energy contributes with a force opposed to the collapse. Finally, the third term,
(4.18) 
satisfies since . The term can be interpreted as a force of quantum origin due to the antiscreening effect of virtual gravitons. This force is opposed to the collapse, however, it is easy to see that its effect is only relevant for .
5 Results from the numerical integration
In order to get the evolution of the shell we should numerically integrate (4.17) using that, from (4.5),
(5.1) 
We should also take into account that is a piecewise defined function satisfying

if the areal radius of the shell satisfies . I.e., according to the tunneling picture, when the collapse has not still reached the stage in which a dynamical outer horizon appears, no Hawking radiation exists.

if , where we should use (3.7) in order to take into account the radiation coming from the horizon which is absorbed by the shell. I.e., once the collapse of the shell generates an outer horizon, negative energy Hawking radiation diminishes the shell mass. The absorption coefficient informs us about the fraction of energy that is absorbed by the shell. In our model we will consider .

if , where we should use (3.7) in order to take into account the radiation coming from both the horizons and . I.e. when the collapse of the shell also generates an inner horizon then the negative energy coming from both the outer and the inner horizon contribute to the mass loss of the shell.
The numerical integration of the system of differential equations composed of the evolution equation (4.17) and the mass equation (5.1) during the collapse of the shell from its initial areal radius until is exemplified in fig. 2. In this figure we can see how the shell collapses creating an outer horizon . However, the figure does not show neither the generation of the inner horizon nor the shell’s later behaviour since they are indistinguishable at the chosen drawing scale. In the same figure we have also plotted the evolution of the mass. It is important to note that, even if negative energy Hawking radiation is being absorbed once the shell generates its outer horizon, the initial mass is so big and the collapse so fast that the mass remains practically constant at this stage. In this way, one can conclude that the behaviour of the heavy massive shell previous to the generation of the inner horizon is practically identical to the behaviour of a collapsing classical shell.
A continuation of fig. 2 is shown in figure 3. It describes the situation when the antiscreening effect of the virtual gravitons have already forced the generation of an inner horizon . We specifically show the evolution of the shell beyond this horizon. As has been exemplified in the graphic, the numerical integration of the system of differential equations provides us with a generic rebound of the shell which can be interpreted as due to the repulsive antiscreening force (4.18).
After the rebound takes place the shell again approaches , what implies that tends to zero. Then, since , equation (4.5) tells us that will diverge and that the shell radius will reach the inner horizon (see fig. 1).
If we just consider the results obtained so far
From the point of view of the shell things look quite different. We have seen that the shell collapses with negligible mass decrease while . However, in its short travel in region III the shell reaches the inner horizon after the rebound while managing to reduce its mass from a large value down to the Planckian value (in agreement with the external observer). It is not hard to understand this fast decrease from the point of view of the shell since we have with diverging at . In figure 4 we show the results of the mass evolution as a function of the shell proper time. There one can visualize how it is only when the shell is very close to the inner horizon that all the negative energy from the horizons is received in the form of an sudden implosion of negative energy radiation that causes the shell to lose most of its mass.
By collecting our results we have drawn the complete picture of the gravitational collapse of the shell in a Penrose diagram (fig. 5). Since the improved solution has an endogenous instability at [13] the actual behaviour of the spacetime around and beyond this horizon is not clear. This has been represented by a dashed line in the figure.
6 Conclusions
In this letter we have described the collapse of a heavy massive thin shell. We have commented that only for large enough masses the behaviour is expected to be similar to the strictly classical one during the collapsing phase. In this way, the generation of an (outer) horizon will be compulsory. We have seen that, once the horizon has been generated, a negligible amount of Hawking radiation will leave the horizon surroundings towards the future null infinity (in the form of positive energy particles –mostly photons) and towards the collapsing shell (in the form of negative energy particles –again, mostly photons) during this phase. There are two reasons why Hawking radiation is negligible at this stage. First, the numerical computations show that the time taken for the shell between the generation of the outer and the inner horizon is (and should be (4.6)) small. Second, the initial mass of the modeled shell is large and, therefore, Hawking radiation from the outer horizon (that approximately behaves as at this stage) is negligible.
We have seen that there is a repulsive force acting on the shell due to the antiscreening effect of the virtual gravitons. However, this force only becomes relevant for , i.e., its effects appear while the shell is gravitationally trapped inside the outer horizon (). Therefore, from the dynamical point of view modeled by the collapsing shell, it can be interpreted that this quantum repulsive force is the ultimate responsible that allows the shell to stop being gravitationally trapped or, equivalently, that causes the later behaviour after generating an inner horizon.
Moreover, since and for , the repulsive force (4.18) eventually creates an impassable barrier when due to the divergence of . In this way, there will always be a rebound point for the shell underneath the inner horizon . This rebound has been explicitly displayed with the help of some numerical computations in fig. 3. Therefore, our results point towards the conclusion that the total collapse of the shell and the subsequent creation of a singularity could be avoided.
Once the rebound point has been reached, the evolution continues with an expanding areal radius that approaches the inner horizon , but this time with a diverging exterior null time . We have seen that, from the point of view of the shell, this means that it is only when it is very close to this horizon that it suddenly receives a huge amount of negative energy coming from the outer and inner horizons, what would be interpreted by the shell as an implosion of negative energy radiation. Due to this fact, the shell should lose most (or all) of its mass almost instantaneously according to its own proper time. Strictly speaking, the model indicates that the mass should suddenly reach the Planckian value . However, this cannot be guaranteed since the exterior region has an endogenous instability [13] at the inner horizon () whose resolution is not at all clear. Thus, the ultimate result of the collapse (remnant, total evaporation, etc.) remains unclear. In our opinion, the resolution of the inner horizon instability and the knowledge of the ultimate result of the collapse would probably require the help of a full Theory of Quantum Gravity.
Footnotes
 Email: ramon.torresherrera@upc.edu
 Email: f.fayos@upc.edu
 On the contrary, a light shell with would possess a markedly quantum nature [3][4].
 Of course, one can also have the positive energy photon tunneling outwards in region I’ (III’, respectively).
 It is important to remark that in this section we are dealing only with the improved Schwarzschild solution. In the complete collapsing model, Hawking radiation will be directed towards the shell and the expression equivalent to (3.11) will have to be reconsidered. See section 5.
 I.e., here we are not taken into account possible instabilities and quantum gravity effects that will be commented later.
References
 Alberghi G L, Casadio R, Vacca G P and Venturi G 1999 Class. Quantum Grav. 16 131
 Ortíz L and Ryan M P 2007 Gen. Rel. Grav. 39 1087
 Hajicek P 1992 Commun. Math. Phys. 150 545
 Corichi A, CruzPacheco G, Minzoni A, Padilla P, Rosenbaum M, Ryan M P, Smyth N F and Vukasinac T 2002 Phys. Rev. D 65 064006
 Bonanno A and Reuter M 2000 Phys. Rev. D 62 043008
 Israel W 1966 Il Nuovo Cimento 44 B 1
 Bonanno A and Reuter M 2006 Phys. Rev. D 73 083005
 BjerrumBohr N E J, Donoghue J F and Holstein B R 2003 Phys. Rev. D 68 084005
 Modesto L 2010 Int. J. Theor. Phys. 49 1649
 AmelinoCamelia G 2006 Class. Quantum Grav. 23 2585
 Nicolini P 2009 Int. J. Phys. A 24 1229
 Torres R, Fayos F and LorenteEspín O 2013 Phys. Lett. B 720 198
 Torres R 2013 Phys. Lett. B 724 338
 Parikh M K and Wilczek F 2001 Phys. Rev. Lett. 85 5042 (arXiv:hepth/9907001v3)
 Painlevé P 1921 C. R. Acad. Sci. (Paris) 173 677
 KeskiVakkuri E and Kraus P 1997 Nucl. Phys. B491 249 (arXiv:hepth/9610045)
 Kraus P and Wilczeck F 1994 Nucl. Phys. B433 665 (arXiv:hepth/9408003); Nucl. Phys. B437 231 (arXiv:hepth/9411219)
 Birrell N D and Davies P C W 1982 Quantum Fields in Curved Space Cambridge: University Press
 Fabbri A and NavarroSalas J 2005 Modeling Black Hole Evaporation London: Imperial College Press
 Poisson E 2004 A relativist’s toolkit Cambridge: Cambridge University Press (UK)