Power-Law Approach of the Stress-Energy Tensor to the Unruh State after Gravitational Collapse

Michael Wilson

arxiv: 2604.26840 · v1 · submitted 2026-04-29 · 🌀 gr-qc · hep-th· math-ph· math.MP

Power-Law Approach of the Stress-Energy Tensor to the Unruh State after Gravitational Collapse

Michael Wilson This is my paper

Pith reviewed 2026-05-07 12:34 UTC · model grok-4.3

classification 🌀 gr-qc hep-thmath-phmath.MP

keywords stress-energy tensorUnruh statenull shell collapselate-time asymptoticsbranch cutPrice's lawscalar field

0 comments

The pith

The renormalized stress-energy tensor approaches the Unruh state as a power law t_s^{-3} after gravitational collapse, with nonzero coefficient at null infinity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes the late-time decay rate of the difference between the in-vacuum and Unruh-state stress-energy tensor for a massless scalar field in a null-shell collapse spacetime. At finite exterior radius an upper bound of order t_s^{-3} is derived from the Hadamard difference and the branch-cut structure of the retarded Green function. At future null infinity the leading coefficient C_uu in the u_s^{-3} term is shown to be nonzero and negative by direct evaluation of the small-frequency branch-cut residue combined with Planck suppression at high frequencies. This power law originates from the same ω² ln ω branch point in the ℓ=0 Wronskian that underlies Price's law. The result confirms earlier numerical conjectures and provides an analytic understanding of the approach to the Unruh state.

Core claim

We show that Δ⟨T_uu⟩ at future null infinity behaves as C_uu u_s^{-3} with C_uu ≠ 0 as u_s → ∞. The coefficient is fixed by computing the branch-cut residue at small frequency explicitly and demonstrating that high-frequency contributions are Planck-suppressed, leaving a dominant term of definite sign. The decay exponent is set by the ω² ln ω singularity in the Wronskian of the zero-mode radial wave equation.

What carries the argument

The branch-cut residue arising from the ω² ln ω singularity in the Wronskian of the ℓ=0 radial wave equation, which determines the leading late-time power-law correction to the Unruh state.

If this is right

The difference |Δ⟨T_μν⟩| is bounded by C(r) t_s^{-3} at any finite radius.
The uu-component at null infinity has a nonzero negative leading coefficient C_uu.
The same singularity structure is conjectured to produce a t_s^{-7} decay for gravitational perturbations.
The result rules out exponential or faster decay for the scalar field case.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The persistent 1/u^3 tail may contribute to the semiclassical backreaction on the geometry during the final stages of evaporation.
Extending the calculation to include massive fields could show how mass terms modify the power-law tail.
The negative sign of C_uu implies a specific energy flux direction consistent with outgoing Hawking radiation.
Resolving the gauge issues would allow a similar analysis for spin-2 perturbations and test the conjectured t^{-7} bound.

Load-bearing premise

The small-frequency branch-cut residue provides the dominant contribution to the late-time integral after high-frequency Planck suppression is applied.

What would settle it

Numerical computation of the full frequency integral for the stress-energy tensor at asymptotically late retarded times that finds C_uu = 0 or C_uu > 0 would falsify the claim.

read the original abstract

We establish the rate at which the renormalized stress--energy tensor of a massless minimally coupled scalar field in the in-vacuum state of a collapsing null-shell spacetime approaches the corresponding Unruh-state value. At finite exterior radius, we establish the upper bound \[ |\Delta\langle T_{\mu\nu}\rangle|\leq C(r)\,t_s^{-3} \] from the Cauchy-surface decomposition of the Hadamard difference and the branch-cut structure of the retarded Green function. At future null infinity, we show that the leading coefficient in the late-time expansion \[ \Delta\langle T_{uu}\rangle\sim C_{uu}\,u_s^{-3} \] is nonzero, by computing the branch-cut residue explicitly at small frequency and using the Planck suppression of the thermal spectrum at large frequency to show that the dominant contribution to $C_{uu}$ has a definite sign. The result gives \[ \Delta\langle T_{uu}\rangle\big|_{\Iscr^+}(u_s) \sim C_{uu}\,u_s^{-3}, \qquad u_s\to\infty, \] with $C_{uu}\neq 0$. The exponent is determined by the $\omega^2\ln\omega$ branch-point singularity in the Wronskian of the $\ell=0$ radial wave equation, the same structure responsible for Price's law. The sign $C_{uu}<0$ is supported by a physical argument and by the numerical mode data of Gholizadeh Siahmazgi, Anderson, and Fabbri. The result confirms their conjecture that the approach is a power law. We conjecture that the same mechanism gives an analogous $t_s^{-7}$ bound for gravitational perturbations ($\ell_{\min}=2$), though the extension to the spin-2 case involves gauge issues not addressed here.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper computes an explicit nonzero value for the leading u^{-3} coefficient in the late-time approach of the stress-energy tensor to the Unruh state at null infinity.

read the letter

This paper pins down that the leading coefficient C_uu in Delta ~ C_uu u_s^{-3} at future null infinity is nonzero. They do it by evaluating the branch-cut residue of the small-frequency omega^2 ln omega singularity in the l=0 Wronskian, then arguing that Planck suppression kills the high-frequency tail so the residue sets both the power and the sign of the leading term. At finite radius they also give a clean t_s^{-3} upper bound from the Cauchy-surface split of the Hadamard difference and the retarded Green function structure. The exponent matches the one that produces Price's law, which is a nice consistency check. The sign is backed by the cited numerical mode data rather than derived from scratch, but the nonzero character comes from the direct residue calculation. That turns the earlier conjecture into a concrete result with an explicit nonzero C_uu. The work is careful about the in-vacuum state in the collapsing null-shell geometry and stays within the massless minimally coupled scalar case. The main soft spot is the step that claims the high-frequency integral is negligible at order 1/u^3 once the exponential cutoff is applied. If the mode-sum prefactors grow with omega before the Planck factor cuts in, or if contour deformation picks up other contributions of comparable size, the residue would not automatically dominate. The paper notes this relies on the specific integral representation, so that part needs a close look. The gravitational perturbation conjecture is left open because of gauge issues, which is reasonable. This is for people who need controlled late-time expansions when computing back-reaction or refining Hawking radiation calculations in evaporating black holes. It supplies a missing explicit coefficient that several groups have been waiting for. The math is grounded enough and the new result is sharp enough that it deserves a serious referee to verify the frequency integral and contour steps.

Referee Report

1 major / 1 minor

Summary. The paper derives the late-time approach of the renormalized stress-energy tensor for a massless minimally coupled scalar field in the in-vacuum state of a collapsing null-shell spacetime to the Unruh state. It proves an upper bound |Δ⟨T_μν⟩| ≤ C(r) t_s^{-3} at finite exterior radius via Cauchy-surface decomposition of the Hadamard difference and the branch-cut structure of the retarded Green function. At future null infinity it shows that the leading coefficient C_uu in Δ⟨T_uu⟩ ∼ C_uu u_s^{-3} is nonzero and negative by explicit evaluation of the ω² ln ω branch-cut residue in the ℓ=0 Wronskian at small frequency, combined with Planck suppression of the high-frequency tail in the Unruh thermal spectrum; the sign is further supported by numerical mode data. The exponent is traced to the same singularity that produces Price's law, and a conjecture is offered for an analogous t_s^{-7} bound in the spin-2 case.

Significance. If the central claims hold, the work supplies a parameter-free analytic confirmation of the power-law decay to the Unruh state, directly linking the decay rate and the nonzero coefficient to the explicit small-frequency residue of the radial Wronskian. The combination of Cauchy-surface decomposition, contour-residue evaluation, and external numerical data for the sign constitutes a concrete advance over purely numerical conjectures. The result strengthens the connection between classical tail behavior (Price's law) and quantum stress-tensor asymptotics in dynamical black-hole backgrounds.

major comments (1)

[late-time expansion at future null infinity] The argument that the small-frequency branch-cut residue determines the leading 1/u_s^3 coefficient at I^+ rests on the claim that Planck suppression renders the high-frequency tail negligible. The frequency integral representation of Δ⟨T_uu⟩(u_s) obtained from the Cauchy-surface decomposition must be shown explicitly to admit a splitting or contour deformation in which any polynomial or sub-exponential growth in ω is overpowered by the exponential Planck factor, without additional poles or branch cuts contributing at the same order. Without this control the dominance of the explicit ℓ=0 residue (and therefore the nonzero value and sign of C_uu) is not yet secured.

minor comments (1)

[final paragraph] The conjecture for gravitational perturbations (t_s^{-7} bound) is stated without derivation; a brief outline of the expected ℓ_min=2 Wronskian singularity and the gauge issues would clarify the scope of the claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comment on the late-time expansion at future null infinity. We address the point below and will strengthen the relevant section in the revised version.

read point-by-point responses

Referee: The argument that the small-frequency branch-cut residue determines the leading 1/u_s^3 coefficient at I^+ rests on the claim that Planck suppression renders the high-frequency tail negligible. The frequency integral representation of Δ⟨T_uu⟩(u_s) obtained from the Cauchy-surface decomposition must be shown explicitly to admit a splitting or contour deformation in which any polynomial or sub-exponential growth in ω is overpowered by the exponential Planck factor, without additional poles or branch cuts contributing at the same order. Without this control the dominance of the explicit ℓ=0 residue (and therefore the nonzero value and sign of C_uu) is not yet secured.

Authors: We agree that an explicit splitting strengthens the argument. In the revised manuscript we will insert a dedicated paragraph after the frequency-integral representation (currently Eq. (4.12) or equivalent) that splits the ω-integral at a fixed cutoff ω_0 > 0 independent of u_s. The low-frequency piece (0 < ω < ω_0) is controlled by the explicit ω² ln ω branch-point residue of the ℓ=0 Wronskian; its contribution is exactly C_uu u_s^{-3} plus higher-order terms that vanish faster. The high-frequency piece (ω > ω_0) is bounded using the exponential decay of the Unruh factor |1 - e^{-2πω/κ}|^{-1} ≲ e^{-2πω/κ} together with the known polynomial growth (at most ω^2) of the radial mode functions and the retarded Green function for fixed r. The resulting integral is then O(u_s^{-N}) for any N or exponentially small, hence o(u_s^{-3}). No other poles or branch cuts enter at order u_s^{-3} because the only singularity producing a 1/u^3 tail is the ℓ=0 branch point (higher-ℓ modes yield faster decay by Price’s law). This decomposition is justified by the absolute convergence of the integral for each fixed u_s and by the analytic continuation properties already established in Sec. 3. We will also note that the sign of C_uu remains negative once the high-frequency remainder is shown to be negligible, consistent with the numerical mode data cited in the paper. revision: yes

Circularity Check

0 steps flagged

No significant circularity; explicit residue computation and external numerical data establish nonzero C_uu

full rationale

The paper derives the leading late-time coefficient C_uu by direct evaluation of the ω² ln ω branch-point residue in the ℓ=0 Wronskian of the radial wave equation, combined with the standard Planck factor from the Unruh thermal spectrum to suppress the high-frequency tail. This is an explicit analytic step, not a fit or redefinition. The sign is corroborated by a physical argument plus external numerical mode data from Gholizadeh Siahmazgi et al. (distinct authors). No self-citation is load-bearing for the central claim, no ansatz is smuggled, and the result does not rename a known empirical pattern by construction. The derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The derivation rests on standard properties of retarded Green functions and Hadamard states in QFT on curved spacetime; no new free parameters, ad-hoc constants, or postulated entities are introduced.

axioms (2)

standard math The retarded Green function for the massless scalar field in the collapsing null-shell spacetime possesses a branch-cut structure in the complex frequency plane whose small-frequency behavior is governed by the Wronskian of the radial wave equation.
Invoked to extract the late-time power-law decay from the contour integral representation of the Hadamard difference.
domain assumption The in-vacuum state on the collapsing null-shell background and the definition of the Unruh state are well-defined Hadamard states whose difference can be expressed via the retarded Green function.
Standard setup in quantum field theory on curved spacetime for this geometry.

pith-pipeline@v0.9.0 · 5637 in / 1526 out tokens · 67467 ms · 2026-05-07T12:34:20.379305+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

[1]

S. W. Hawking, Commun. Math. Phys.43, 199 (1975); Erratum:46, 206 (1976)

work page 1975
[2]

W. G. Unruh, Phys. Rev. D14, 870 (1976)

work page 1976
[3]

Fabbri and J

A. Fabbri and J. Navarro-Salas,Modeling Black Hole Evaporation(Imperial College Press, London, 2005)

work page 2005
[4]

Late-time behavior of scalar field modes for a collapsing null shell spacetime and for the Unruh state in Schwarzschild spacetime,

S. Gholizadeh Siahmazgi, P. R. Anderson, and A. Fabbri, “Late-time behavior of scalar field modes for a collapsing null shell spacetime and for the Unruh state in Schwarzschild spacetime,” Phys. Rev. D112, 045013 (2025)

work page 2025
[5]

Gholizadeh Siahmazgi, P

S. Gholizadeh Siahmazgi, P. R. Anderson, and A. Fabbri, WakeSpace Scholarship (2025),http://hdl.handle.net/10339/111099

work page 2025
[6]

R. H. Price, Phys. Rev. D5, 2419 (1972)

work page 1972
[7]

Wilson, (2025)

M. Wilson, “Threshold resolvent singularities and the infrared structure of linearized gravity,” arXiv:2511.05345 [math-ph] (2025) [preprint; under review]. 11

work page arXiv 2025
[8]

E. W. Leaver, Phys. Rev. D34, 384 (1986)

work page 1986
[9]

E. S. C. Ching, P. T. Leung, W. M. Suen, and K. Young, Phys. Rev. Lett.74, 4588 (1995)

work page 1995
[10]

E. S. C. Ching, P. T. Leung, W. M. Suen, and K. Young, Phys. Rev. D52, 2118 (1995)

work page 1995
[11]

Andersson, Phys

N. Andersson, Phys. Rev. D55, 468 (1997)

work page 1997
[12]

Dafermos and I

M. Dafermos and I. Rodnianski, Invent. Math.162, 381 (2005)

work page 2005
[13]

Regge and J

T. Regge and J. A. Wheeler, Phys. Rev.108, 1063 (1957)

work page 1957
[14]

Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, Oxford, 1983)

S. Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, Oxford, 1983)

work page 1983
[15]

W. A. Hiscock, Phys. Rev. D23, 2813 (1981)

work page 1981
[16]

M. R. R. Good, P. R. Anderson, and C. R. Evans, Phys. Rev. D94, 065010 (2016)

work page 2016
[17]

P. R. Anderson, S. Gholizadeh Siahmazgi, R. D. Clark, and A. Fabbri, Phys. Rev. D 102, 125035 (2020)

work page 2020
[18]

P. R. Anderson, R. D. Clark, A. Fabbri, and M. R. R. Good, Phys. Rev. D100, 061703(R) (2019)

work page 2019
[19]

N. D. Birrell and P. C. W. Davies,Quantum Fields in Curved Space(Cambridge Uni- versity Press, Cambridge, 1982)

work page 1982
[20]

R. M. Wald,Quantum Field Theory in Curved Spacetime and Black Hole Thermody- namics(University of Chicago Press, Chicago, 1994)

work page 1994
[21]

Theorems on the uniqueness and thermal properties of stationary, nonsingular, quasifree states on spacetimes with a bifurcate Killing horizon,

B. S. Kay and R. M. Wald, “Theorems on the uniqueness and thermal properties of stationary, nonsingular, quasifree states on spacetimes with a bifurcate Killing horizon,” Phys. Rep.207, 49 (1991)

work page 1991
[22]

Vacuum polarization in Schwarzschild spacetime,

P. Candelas, “Vacuum polarization in Schwarzschild spacetime,” Phys. Rev. D21, 2185 (1980)

work page 1980
[23]

On the derivation of Hawking radiation associated with the formation of a black hole,

K. Fredenhagen and R. Haag, “On the derivation of Hawking radiation associated with the formation of a black hole,” Commun. Math. Phys.127, 273 (1990)

work page 1990
[24]

Analogue quantum simulation of the Hawking effect in a polariton superfluid,

M. J. Jacquet, L. Giacomelli, Q. Valnais, M. Joly, F. Claude, E. Giacobino, Q. Glorieux, I. Carusotto, and A. Bramati, “Analogue quantum simulation of the Hawking effect in a polariton superfluid,” Eur. Phys. J. D76, 152 (2022). 12

work page 2022

[1] [1]

S. W. Hawking, Commun. Math. Phys.43, 199 (1975); Erratum:46, 206 (1976)

work page 1975

[2] [2]

W. G. Unruh, Phys. Rev. D14, 870 (1976)

work page 1976

[3] [3]

Fabbri and J

A. Fabbri and J. Navarro-Salas,Modeling Black Hole Evaporation(Imperial College Press, London, 2005)

work page 2005

[4] [4]

Late-time behavior of scalar field modes for a collapsing null shell spacetime and for the Unruh state in Schwarzschild spacetime,

S. Gholizadeh Siahmazgi, P. R. Anderson, and A. Fabbri, “Late-time behavior of scalar field modes for a collapsing null shell spacetime and for the Unruh state in Schwarzschild spacetime,” Phys. Rev. D112, 045013 (2025)

work page 2025

[5] [5]

Gholizadeh Siahmazgi, P

S. Gholizadeh Siahmazgi, P. R. Anderson, and A. Fabbri, WakeSpace Scholarship (2025),http://hdl.handle.net/10339/111099

work page 2025

[6] [6]

R. H. Price, Phys. Rev. D5, 2419 (1972)

work page 1972

[7] [7]

Wilson, (2025)

M. Wilson, “Threshold resolvent singularities and the infrared structure of linearized gravity,” arXiv:2511.05345 [math-ph] (2025) [preprint; under review]. 11

work page arXiv 2025

[8] [8]

E. W. Leaver, Phys. Rev. D34, 384 (1986)

work page 1986

[9] [9]

E. S. C. Ching, P. T. Leung, W. M. Suen, and K. Young, Phys. Rev. Lett.74, 4588 (1995)

work page 1995

[10] [10]

E. S. C. Ching, P. T. Leung, W. M. Suen, and K. Young, Phys. Rev. D52, 2118 (1995)

work page 1995

[11] [11]

Andersson, Phys

N. Andersson, Phys. Rev. D55, 468 (1997)

work page 1997

[12] [12]

Dafermos and I

M. Dafermos and I. Rodnianski, Invent. Math.162, 381 (2005)

work page 2005

[13] [13]

Regge and J

T. Regge and J. A. Wheeler, Phys. Rev.108, 1063 (1957)

work page 1957

[14] [14]

Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, Oxford, 1983)

S. Chandrasekhar,The Mathematical Theory of Black Holes(Oxford University Press, Oxford, 1983)

work page 1983

[15] [15]

W. A. Hiscock, Phys. Rev. D23, 2813 (1981)

work page 1981

[16] [16]

M. R. R. Good, P. R. Anderson, and C. R. Evans, Phys. Rev. D94, 065010 (2016)

work page 2016

[17] [17]

P. R. Anderson, S. Gholizadeh Siahmazgi, R. D. Clark, and A. Fabbri, Phys. Rev. D 102, 125035 (2020)

work page 2020

[18] [18]

P. R. Anderson, R. D. Clark, A. Fabbri, and M. R. R. Good, Phys. Rev. D100, 061703(R) (2019)

work page 2019

[19] [19]

N. D. Birrell and P. C. W. Davies,Quantum Fields in Curved Space(Cambridge Uni- versity Press, Cambridge, 1982)

work page 1982

[20] [20]

R. M. Wald,Quantum Field Theory in Curved Spacetime and Black Hole Thermody- namics(University of Chicago Press, Chicago, 1994)

work page 1994

[21] [21]

Theorems on the uniqueness and thermal properties of stationary, nonsingular, quasifree states on spacetimes with a bifurcate Killing horizon,

B. S. Kay and R. M. Wald, “Theorems on the uniqueness and thermal properties of stationary, nonsingular, quasifree states on spacetimes with a bifurcate Killing horizon,” Phys. Rep.207, 49 (1991)

work page 1991

[22] [22]

Vacuum polarization in Schwarzschild spacetime,

P. Candelas, “Vacuum polarization in Schwarzschild spacetime,” Phys. Rev. D21, 2185 (1980)

work page 1980

[23] [23]

On the derivation of Hawking radiation associated with the formation of a black hole,

K. Fredenhagen and R. Haag, “On the derivation of Hawking radiation associated with the formation of a black hole,” Commun. Math. Phys.127, 273 (1990)

work page 1990

[24] [24]

Analogue quantum simulation of the Hawking effect in a polariton superfluid,

M. J. Jacquet, L. Giacomelli, Q. Valnais, M. Joly, F. Claude, E. Giacobino, Q. Glorieux, I. Carusotto, and A. Bramati, “Analogue quantum simulation of the Hawking effect in a polariton superfluid,” Eur. Phys. J. D76, 152 (2022). 12

work page 2022