Thunderbolt Over the Horizon: The Collapse of Semiclassical Gravity in Black Hole Evaporation ✦
We study the formation and evaporation of a black hole in a four-dimensional semiclassical model with diffeomorphism invariance, taking into account a one-loop trace anomaly. The Einstein equations with backreaction are solved for spherically symmetric collapse of a null shell. The results demonstrate the emergence of a spacelike singularity of the thunderbolt type after the apparent horizon retreats; the singularity extends far beyond the hole, into a region with parametrically small semiclassical curvature. The effect is caused by a nonlinear instability of equations with higher-derivative terms and is typical for models with anomalous quantum corrections. The appearance of a thunderbolt indicates a breakdown of the semiclassical effective field theory at macroscopic scales and undermines the standard formulation of the information paradox.
Context
The standard black hole evaporation scenario, proposed by Stephen Hawking, relies on the semiclassical gravity approximation, where quantum fields propagate on a fixed curved background. However, including the backreaction of radiation on the geometry—key to resolving the information paradox—has remained technically out of reach in four dimensions. Understanding how a black hole evaporates and what's left behind is essential for unifying quantum mechanics and general relativity. The work by Lowe and Thorlacius fills this gap, using a model that captures the quantum conformal anomaly—a fundamental effect linking geometry to quantum fields.
Methods
The researchers considered a spherically symmetric collapse of a null shell that forms a black hole, and solved the full system of semiclassical equations of motion with backreaction. They used the Riegert action, which localizes the anomalous contributions via an auxiliary scalar field. The equations, written in double-null coordinates, were integrated using the method of characteristics with a fourth-order Runge-Kutta algorithm on a rectangular grid—a prime example of numerical modeling in general relativity. This approach allowed them to track the formation of the horizon, Hawking radiation, and the subsequent dynamics up to late times, including the region beyond the evaporated hole. Crucially, they preserved the theory's diffeomorphism invariance, ensuring correct treatment of the speed of light and causal structure.
Results
The main result is the emergence of a spacelike 'thunderbolt' singularity after the apparent horizon recedes and the black hole effectively disappears. This singularity develops in a region where classical spacetime curvature is negligibly small, and spreads over macroscopic distances, reaching spacelike infinity. The Kretschmann scalar, an invariant measure of curvature, shows a sharp surge and sign change, pointing to a catastrophic divergence. Numerical tests confirm this is not a grid artifact but a consequence of a nonlinear instability inherent to the fourth-order semiclassical equations. Interestingly, the metric itself near the singularity behaves universally: the exponential factor tends to zero while the radial function grows. This is reminiscent of redshift behavior near a horizon, but here it becomes extreme.
Implications
These results mean that semiclassical gravity—a standard tool of theoretical physics—predicts its own failure at macroscopic scales, not just in the Planck regime. This strikes at the formulation of the information paradox, which assumed the semiclassical description holds arbitrarily far from a black hole. Consequently, any arguments about the fate of quantum information based solely on semiclassics lose their force. It may be necessary to introduce new physical principles into effective field theory to avoid such catastrophic singularities. Curiously, a connection emerges with black hole entropy, first described by Jacob Bekenstein.
Future development
Future research could go in several directions. First, testing the hypothesis that fine-tuning the initial data could eliminate the 'thunderbolt'—though the authors think it unlikely without modifying the theory. Second, searching for generalized models that smoothly continue the evolution beyond the semiclassical breakdown. Third, exploring observational consequences: if such singularities arise in real black holes, their evaporation could be accompanied by macroscopic catastrophes, placing constraints on primordial black holes. More fundamentally, this opens a new window into quantum gravity: any viable candidates (string theory, loop gravity) must explain how to avoid the 'thunderbolt'. It's not ruled out that the uncertainty principle, which underpins Hawking radiation, could play a key role.
Impact
The discovery will impact quantum gravity, astrophysics, and quantum information. In particular, it could change our understanding of how black holes interact with their surroundings and challenge scenarios of their time evolution in cosmology.
Next steps
In the near term, it's necessary to analyze the behavior of the solution for larger black hole masses and to check whether the instability actually develops on a light-crossing timescale rather than the Hawking evaporation timescale. It is also important to investigate whether higher-order quantum corrections could alter the dynamics, perhaps building on ideas developed by William Unruh.
Key open problems
The work is directly connected to the Hawking information paradox and the problem of quantizing gravity. It demonstrates that the semiclassical approach is not self-consistent even in a weak field, raising the question of the limits of effective field theories in curved spacetime. This failure resonates with unresolved questions about quantum information and entropy in the presence of gravity.
🎯 The term 'thunderbolt' was coined by Hawking and Stewart in 1993, but back then it was thought not to appear in two-dimensional models. Now it turns out that four-dimensional gravity with quantum corrections inevitably leads to this dramatic finale.
Key numbers
- mass of the black hole in the model: 2.2 Planck units
- anomalous coefficients a, b, c: 10
- singularity development time: ~10 light-crossing times
- numerical grid size: up to 20,000 × 4,000 nodes
- convergence accuracy: <1%