Complete Quantum Stress Tensor Inside a Four Dimensional Schwarzschild Black Hole: A Divergent Focusing Source

Technical Summary: Complete Quantum Stress Tensor Inside a Four-Dimensional Schwarzschild Black Hole

Problem Statement
The renormalized stress-energy tensor (RSET), $\langle T_{ab} \rangle_{\text{ren}}$, serves as the local source term in the semiclassical Einstein equation ($G_{ab} = 8\pi G \langle T_{ab} \rangle_{\text{ren}}$). While the RSET outside black holes has been extensively studied, a controlled, complete RSET throughout the entire interior of a four-dimensional black hole has remained unavailable. Previous interior calculations were limited to specific observables (e.g., vacuum polarization $\langle \Phi^2 \rangle_{\text{ren}}$ in restricted regions, trace information, or selected flux components at Cauchy horizons) and could not provide the full tensorial structure required to determine local energy density, principal pressures, and the complete null projection entering the Raychaudhuri equation. Furthermore, no controlled renormalized local observable had been obtained in the deep ultraviolet (UV) neighborhood of a four-dimensional spacelike black-hole singularity, leaving the nature of the leading quantum source (transported Hawking flux vs. local vacuum polarization) and its effect on singularity focusing unknown.

Methodology
The authors compute the complete RSET for a massless minimally coupled scalar field in both the Unruh and Hartle–Hawking states throughout the Schwarzschild interior ($0 < r < 2M$). The calculation employs an angular-splitting renormalization scheme combined with high-order large-$\ell$ asymptotic subtraction.

1. Renormalization: The RSET is defined via point splitting, subtracting the local Christensen–DeWitt–Schwinger term. The authors integrate mode functions over frequency $\omega$ for fixed angular momentum $\ell$, subtracting large-frequency WKB terms at the integrand level.
2. Tail Subtraction: After frequency integration, residual large-$\ell$ self-cancellation terms remain. The authors numerically extract coefficients at sufficiently large $\ell$, subtract them, and verify convergence over a smaller $\ell$ range. This procedure yields results with residuals at the percent level.
3. Validation: The results are validated through three independent checks:
   • Covariant Conservation: $\nabla_a \langle T^{ab} \rangle_{\text{ren}} = 0$.
   • Trace Identity: $\langle T^a_a \rangle_{\text{ren}} = -(1/2)\Box \langle \Phi^2 \rangle_{\text{ren}} + A_{\text{geom}}$.
   • Global Integration: To avoid amplifying numerical noise from differentiation, the conservation and trace equations are rewritten in integral form and compared against the direct numerical results.
   • Exterior Matching: The results smoothly match independent exterior calculations at the horizon.

Key Results
The study determines all independent components of $\langle T^a_b \rangle_{\text{ren}}$ and $\langle \Phi^2 \rangle_{\text{ren}}$ from the event horizon down to $r/M \simeq 10^{-4}$.

1. Scaling Behavior near the Singularity:
   As $r \to 0$, both the Unruh and Hartle–Hawking states approach the same conserved scaling solution:
   $$M^4 \langle T^a_b \rangle_{\text{ren}} \simeq \left(\frac{r}{M}\right)^{-6} \tau^a_b$$
   where $\tau^a_b = \text{diag}(\lambda_t, \lambda_r, \lambda_\perp, \lambda_\perp)$. The numerical coefficients are found to be $\lambda_t \simeq 2.6357 \times 10^{-3}$, $\lambda_r \simeq -1.1926 \times 10^{-3}$, and $\lambda_\perp \simeq 3.3494 \times 10^{-3}$. These coefficients satisfy the constraints imposed by covariant conservation and the trace identity to within 0.3% relative deviation.

2. Vacuum Polarization vs. Hawking Flux:
   The state-dependent Unruh flux component scales as $\langle T^r_t \rangle_{\text{ren}} \sim r^{-2}$, which is suppressed by a factor of $(r/M)^4$ relative to the diagonal mixed components ($\sim r^{-6}$). Consequently, the leading UV source is local vacuum polarization induced by the high-curvature geometry, not transported Hawking radiation. This explains why physically distinct states (Unruh and Hartle–Hawking) converge to the same leading tensor near the singularity.

3. Energy Conditions and Focusing:
   The limiting tensor is strongly anisotropic and violates the dominant energy condition (principal pressures exceed energy density). However, it satisfies the null energy condition. The leading null projection is positive:
   $$T_{kk} \equiv \langle T_{ab} \rangle_{\text{ren}} k^a k^b > 0$$
   This implies that the RSET enters the Raychaudhuri equation with a focusing sign ($R_{ab}k^a k^b > 0$).

4. Renormalization Freedom:
   The leading $r^{-6}$ behavior and the sign of the null projection are invariant under standard finite renormalization freedom. On the Ricci-flat Schwarzschild background, conserved local curvature ambiguities either vanish or are subleading compared to the computed $r^{-6}$ term.

Significance and Claims
The paper claims to close two major gaps in the literature: providing the first complete RSET for a 4D black hole interior and obtaining the first controlled renormalized local observable in the deep UV near a spacelike singularity.

The primary physical conclusion is that, at the level of the complete fixed-background Schwarzschild RSET, the leading semiclassical source does not support the intuition that quantum defocusing smooths the singularity. Instead, it supplies a divergent local focusing source in the local Raychaudhuri equation. The authors emphasize that this is a local, fixed-background result; it does not constitute a proof of a global semiclassical singularity theorem, nor does it preclude singularity resolution via backreacted geometry or effective quantum gravity models. However, it establishes that any proposed defocusing mechanism in the Schwarzschild interior must explicitly overcome the leading $r^{-6}$ focusing source identified here. The work provides a benchmark for future backreaction studies and a concrete target for effective models of singularity resolution.