Self-Consistent Black-Hole Evaporation

A self-consistent evaporation calculation must do more than place Hawking radiation on a prescribed metric: it must determine the quantum state, renormalized stress tensor, and evolving geometry together, while showing where the approximation remains reliable. At the 2026-08-20 evidence cutoff, this programme is successful in controlled two-dimensional and adiabatic settings but remains incomplete for generic four-dimensional collapse through the late endpoint. The central deficit is a controlled global solution, not evidence that semiclassical backreaction is internally inconsistent.

Helpful background. Hawking Flux, Backreaction, and Evaporation Dynamics derives the fixed-background flux, mass-loss estimate, and semiclassical Einstein equation used below.

The coupled mean-field problem

The target equation is

$$G_{ab}[g]+\Lambda g_{ab}+H_{ab}^{\mathrm{ren}}[g] =8\pi G_{\mathrm{ren}}\,\langle T_{ab}[g]\rangle_{\omega,\mathrm{ren}},$$

supplemented by a state prescription for $\omega$, initial and boundary data, and the constraints. The local tensor $H_{ab}^{\mathrm{ren}}$ represents the curvature counterterms required by renormalization. A solution is therefore not specified by a luminosity law $dM/du$ alone. One must show that the stress tensor is conserved on the evolving metric, that the geometry used to compute it is the geometry produced by it, and that neglected fluctuations and higher-curvature terms remain small.

Assessment at 2026-08-20	Meaning
Question status	Open
Evidence status	Insufficient evidence for a generic four-dimensional completion; supported in reduced and adiabatic regimes
Confidence	High in the scope distinction; lower for the fate of the newest construction
Domain	Semiclassical gravity with specified matter, state, symmetry, and renormalization prescription
What could change the conclusion	An independent four-dimensional verification or endpoint evolution

A starting project: reproduce a coupled two-dimensional solution

This is a suitable direction for someone who can vary a gravity–matter action, handle characteristic initial data, and distinguish a renormalized stress tensor from a prescribed flux. Read the linked evaporation background and the CGHS action first. Then use Russo, Susskind, and Thorlacius 1992, §2, Eqs. (2.1)–(2.5), open PDF as a bounded analytic benchmark: the transformed fields, constraints, and incoming-shock solution are specified together.

Reproduce the field redefinition and substitute the exact solution into both the evolution equations and constraints, including the distributional shock matching. Translate the paper’s normalization before comparing geometric quantities. Check the linear-dilaton vacuum and the region where the inverse field redefinition becomes singular. This tests the coupled model much more directly than integrating a four-dimensional mass-loss formula. It remains a two-dimensional semiclassical calculation with a special counterterm.

The output is a derivation and executable residual check, with regular-region residuals separated from matching conditions and strong-coupling boundaries. A laptop and symbolic algebra suffice for this baseline. Extending to state-dependent stress tensors in four dimensions needs QFT-in-curved-spacetime expertise, renormalization control, and numerical-relativity collaboration; it should not be the first implementation milestone.

One candidate extension is a numerical characteristic evolution for a specified smooth incoming pulse in the same model, compared with the analytic vacuum and shock limits. Check prior work and available codes before treating that exercise as novel. Stop or pivot if constraint error does not converge, a boundary condition is being chosen to force the desired answer, or the approximation reaches its validity boundary. A successful solver demonstrates control of the stated model and regime; it does not determine real black-hole endpoints.

What has actually been solved

Two-dimensional dilaton gravity supplies the cleanest analytic model solutions. The CGHS model couples gravity to conformal matter and makes Hawking flux and backreaction analytically tractable Callan et al. 1992, Eqs. (1)–(16). The RST modification adds a specific local term that permits a closed semiclassical treatment and exposes how the endpoint depends on boundary conditions Russo, Susskind, and Thorlacius 1992, §§II–IV. These models establish that a coupled quantum-state/geometry evolution can be coherent. They do not establish that spherical reduction preserves the decisive four-dimensional degrees of freedom or fixes the endpoint.

In four dimensions, the usual description is adiabatic: a large hole passes through a sequence of slowly changing exterior geometries, and a renormalized Unruh-state stress tensor drives a mass loss of order $dM/du\sim-\hbar/M^2$. This is reliable while curvature, fractional mass change per light-crossing time, and stress-tensor fluctuations are controlled. It is a local-in-regime result, not a collapse-to-endpoint spacetime.

Results and open questions at 2026-08-20

Evidence	Result	Domain	Limitation
Renormalized semiclassical and stochastic-gravity framework	Gives the mean equation, noise kernel, and validity diagnostics	Globally hyperbolic backgrounds with a specified Hadamard state	A framework, not a four-dimensional evaporating solution
CGHS and RST models	Solve radiation plus backreaction in $1+1$ dimensions	Large matter number and chosen dilaton model	Dimensional reduction and endpoint prescription are structural assumptions
Fixed-background spectra and mass-loss calculations	Robustly support an outward flux during the large-hole phase	Curvature well below the cutoff and slowly varying exterior	Geometry is not solved self-consistently to the endpoint
2026 quasi-stationary construction	Proposes a closed-form $3+1$ solution with an Unruh-like state	Spherical, quasi-stationary assumptions	New preprint; global validity and independent reproduction remain unsettled

Hu and Verdaguer formulate the key distinction between a small expectation value and small fluctuations: the semiclassical equation is a mean-field equation, and the noise kernel tests whether a sharply defined mean geometry is adequate Hu and Verdaguer 2008, §§2–4, 6. This matters near rapidly evolving or high-curvature regions. A fixed-background phenomenological review, useful for spectra and lifetimes, correspondingly should not be counted as a self-consistent solution Auffinger 2023, §§2–3.

The recent proposal by del Río and Marañón-González raises the frontier by constructing a quasi-stationary four-dimensional geometry and stress tensor together del Río and Marañón-González 2026, pp. 1–5, especially p. 4 and Eqs. (1)–(14), PDF. Its scientific weight depends on checks of renormalization, the quantum state, constraints, approximation hierarchy, and continuation beyond the quasi-stationary domain. Until those checks are independently reproduced, it is important positive evidence rather than a settled global solution.

Validity is a quantitative condition

A persuasive calculation should monitor at least four dimensionless ratios: curvature invariants relative to the ultraviolet scale, the adiabatic time divided by the light-crossing time, the rms stress fluctuation relative to its mean at a stated smearing scale, and the constraint residual relative to the physical source. It should also close the energy balance between Bondi mass loss and flux at infinity. None of these checks can be replaced by regularity of a coordinate chart.

Horizon language also matters. The event horizon is teleological; a trapping or apparent horizon is foliation dependent; and a shrinking marginal surface need not encode the same causal statement. A result about one should not be silently promoted to another. Likewise, an Unruh-like state is selected by collapse-motivated boundary conditions, whereas Hartle–Hawking and Boulware states describe different physics.

Adverse evidence and failed routes

The principal unresolved gap is the absence of a generally accepted $3+1$ evolution from regular collapse through a controlled Planckian endpoint. Substituting a Hawking luminosity into a Vaidya metric does not fill this gap because it prescribes rather than derives the stress tensor. Conversely, a regular stress tensor near one horizon cut does not prove global regularity or small fluctuations.

Dimensional models are sometimes overread in opposite directions. Their solvability neither proves a universal remnant nor proves complete disappearance. The RST endpoint reflects its action and boundary prescription. A four-dimensional proposal can also fail without invalidating the mean-field programme: loss of adiabatic control, state dependence, or a large noise kernel would delimit that solution rather than refute semiclassical gravity everywhere.

What would change the assessment

The status would move toward partially resolved if two independent implementations produced the same four-dimensional causal structure from the same explicit initial state, while satisfying constraints, stress conservation, renormalization conditions, energy balance, and a published error budget. Evolution into a regime where curvature approaches the effective-theory cutoff would still require a controlled matching rule or ultraviolet completion.

A decisive negative result could instead be a theorem showing that no Hadamard state and renormalization prescription can satisfy the coupled equations under clearly physical collapse data, or an unavoidable fluctuation instability before substantial evaporation. Either result must state symmetry, matter content, horizon definition, and domain of dependence.

Research checklist

• Specify dimension, symmetry, matter fields, quantum state, and renormalized couplings.
• Distinguish an imposed luminosity from a computed conserved $\langle T_{ab}\rangle$.
• Evolve the gravitational constraints and quote normalized residuals.
• Compare curvature, adiabatic, fluctuation, and higher-derivative scales.
• Track ADM or Bondi energy together with horizon and asymptotic fluxes.
• State the final time reached and why the approximation is still controlled there.
• Test sensitivity to state, regulator, boundary data, and resolution.

Source and data map

The two-dimensional claims are analytic and trace to the CGHS and RST actions and equations. The validity criterion comes from the stochastic-gravity framework. The four-dimensional frontier claim is a preprint-level construction and should be cited with its version. Fixed-background evaporation tables are comparison inputs, not substitute backreaction data. No independent reproduction was completed for this assessment; a useful rerun would publish the stress-tensor implementation, boundary data, constraint norms, curvature hierarchy, and energy budget.

References

• Auffinger, J. “Primordial Black Hole Constraints with Hawking Radiation—A Review.” Progress in Particle and Nuclear Physics 131, 104040 (2023). doi:10.1016/j.ppnp.2023.104040.
• Callan, C. G., Jr., S. B. Giddings, J. A. Harvey, and A. Strominger. “Evanescent Black Holes.” Physical Review D 45, R1005–R1009 (1992). doi:10.1103/PhysRevD.45.R1005.
• del Río, A., and F. J. Marañón-González. “Black-hole evaporation from the semiclassical Einstein equations in 3+1 dimensions.” arXiv:2607.11737v1 (2026). arXiv.
• Hu, B. L., and E. Verdaguer. “Stochastic Gravity: Theory and Applications.” Living Reviews in Relativity 11, 3 (2008). doi:10.12942/lrr-2008-3.
• Russo, J. G., L. Susskind, and L. Thorlacius. “End Point of Hawking Radiation.” Physical Review D 46, 3444–3449 (1992). doi:10.1103/PhysRevD.46.3444.