Skip to content

Hawking Radiation from Gravitational Collapse

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

Hawking radiation in a collapse spacetime follows from an exponential relation between affine null coordinates. Outgoing rays received at increasingly late retarded time trace backward ever closer to the forming event horizon, so equal late-time intervals correspond to exponentially compressed early-time intervals. Fourier transforming that redshift mixes positive and negative early frequencies and yields TH=κ/(2π)T_H=\kappa/(2\pi) on a fixed, asymptotically stationary background.

Helpful background. Flat-Spacetime Quantum Field Theory provides Bogoliubov transformations. Partial Differential Equations and Green Functions supplies geometric-optics propagation. Causality and Global Structure develops event horizons and null infinity. Quantum Fields in Curved Spacetime provides Hadamard in-states. Black-Hole Geometry and Mechanics supplies surface gravity and horizon generators. Black-Hole Perturbation Theory provides scattering modes and flux normalization.

Consider a free massless scalar on a collapse geometry with a late stationary, nonextremal black-hole exterior. Let uu be retarded time on future null infinity, normalized to asymptotic stationary time; it labels outgoing rays and is not an affine parameter along each ray. Trace an escaping radial wave packet backward through the exterior, through the regular center before collapse, and to past null infinity, whose affine advanced coordinate UU labels the incoming rays. Massive fields require a different asymptotic scattering description because their particles do not reach null infinity.

There is a critical early ray U=UHU=U_H whose outgoing continuation becomes the event horizon. In a collapse model with exponential late peeling, the leading relation is

U(u)=UH−Ae−κu+o(e−κu),A>0.U(u)=U_H-Ae^{-\kappa u} +o(e^{-\kappa u}), \qquad A>0.

For a differentiably controlled remainder, dU/du∼Aκe−κudU/du\sim A\kappa e^{-\kappa u} is positive and tends to zero. The constant AA and subleading terms depend on collapse details, while the leading exponent is fixed by the surface gravity κ\kappa normalized to the asymptotic time generator. A stronger remainder estimate requires more assumptions about how the exterior settles; it is not implied by smoothness alone.

The figure isolates the causal structure in an idealized spherical null-shell collapse. Follow each incoming ray from past null infinity I−\mathscr I^- to the regular center, then outward to future null infinity I+\mathscr I^+. The last escaping rays approach the dashed event horizon. Radial light rays are at 45∘45^\circ in this conformal drawing; neither their drawn spacing nor a coordinate distance measures the exponential redshift.

Incoming rays reflect at the regular center and escape to future null infinity. Rays approaching the critical incoming ray approach the event horizon; the null shell ends at a spacelike singularity and the horizon at future timelike infinity.

Escaping rays in null-shell collapse reflect before the critical incoming ray and remain outside the event horizon. The exponential relation concerns the early affine label UU and late retarded time uu, not compactified distances. I±\mathscr I^\pm denote null infinity, i±i^\pm timelike infinity, and i0i^0 spatial infinity. Fixed-background radial conformal schematic, not to scale; no evaporation is modeled.

A late outgoing mode has phase e−iωue^{-i\omega u} with ω>0\omega>0. Using the leading map and tracing it to the early region gives

pω(U)∝Θ(UH−U)(UH−UA)iω/κ.p_\omega(U)\propto \Theta(U_H-U) \left(\frac{U_H-U}{A}\right)^{i\omega/\kappa}.

The logarithmic phase has a branch point at UHU_H. It is not purely positive frequency with respect to the regular affine coordinate UU. Expand it in early modes e−iω′Ue^{-i\omega'U} and their conjugates:

pω=∫0∞dω′(αωω′fω′+βωω′fω′∗).p_\omega =\int_0^\infty d\omega' \left( \alpha_{\omega\omega'}f_{\omega'} +\beta_{\omega\omega'}f_{\omega'}^* \right).

To see the thermal factor explicitly, set x=UH−U>0x=U_H-U>0, b=ω/κb=\omega/\kappa, and regulate the two Fourier integrals by ϵ>0\epsilon>0:

I±(ω′)=∫0∞dx xibe±iω′x−ϵx=Γ(1+ib)(ϵ∓iω′)−1−ib.I_\pm(\omega') =\int_0^\infty dx\,x^{ib}e^{\pm i\omega'x-\epsilon x} =\Gamma(1+ib)(\epsilon\mp i\omega')^{-1-ib}.

The harmless phase A−ibe±iω′UHA^{-ib}e^{\pm i\omega'U_H} and common flux-normalization factors cancel from the ratio. Projection onto positive early frequency gives α∝I−\alpha\propto I_-; projection onto negative early frequency gives β∝I+\beta\propto I_+. On the principal branch, arg⁡(ϵ±iω′)→±π/2\arg(\epsilon\pm i\omega')\to\pm\pi/2 for ω′>0\omega'>0. Thus ∣I−∣2|I_-|^2 contains e+πbe^{+\pi b} and ∣I+∣2|I_+|^2 contains e−πbe^{-\pi b}, yielding

∣βωω′∣2=e−2πω/κ∣αωω′∣2|\beta_{\omega\omega'}|^2 =e^{-2\pi\omega/\kappa} |\alpha_{\omega\omega'}|^2

up to the common continuum normalization. Combining this ratio with the Bogoliubov norm identity produces the Planck factor

⟨Nω⟩=1e2πω/κ−1\langle N_\omega\rangle =\frac{1}{e^{2\pi\omega/\kappa}-1}

before the exterior scattering potential is included. The original collapse calculation and its wave-packet normalization are in Hawking 1975, §2, Eqs. (2.16)–(2.29), pp. 208–211, with the published correction in Hawking 1976, erratum p. 206.

Comparison with (eω/T−1)−1(e^{\omega/T}-1)^{-1} gives

TH=κ2π.T_H=\frac{\kappa}{2\pi}.

For Schwarzschild with physical mass MM, κ=1/(4GM)\kappa=1/(4GM) and TH=1/(8πGM)T_H=1/(8\pi GM) in the site’s units. In SI units the temperature is TH=ℏc3/(8πGMkB)T_H=\hbar c^3/(8\pi G M k_B). Its factor of ℏ\hbar shows that the radiation is a quantum effect even though the exponential ray map comes from classical geometry.

Exact monochromatic modes run for infinite time and carry delta-function normalization. Hawking wave packets localize both frequency and retarded time, making the late-time limit and particle count finite per packet. The calculation must use the same affine normalizations in the Klein–Gordon product, the ray-tracing map, and the definition of κ\kappa.

The exterior angular-momentum barrier multiplies the occupation by a transmission probability. It does not change the near-horizon temperature, but it makes the spectrum at infinity nonblackbody. That step is treated on Greybody Factors and Asymptotic Flux.

The derivation assumes a regular in-vacuum of Hadamard short-distance form, smooth collapse, a nonextremal horizon, a late approximately stationary exterior, and geometric optics for the large backward blueshift. A dimension check requires ω/κ\omega/\kappa to be dimensionless. A sign check requires dU/du>0dU/du>0 and U→UHU\to U_H from the escaping side as u→∞u\to\infty. Reversing the exponential sign would spread rather than compress late rays and destroy the thermal continuation.

The large backward frequency

ωearly∼ωdudU=ωκ(UH−U)\omega_{\mathrm{early}} \sim\omega\frac{du}{dU} =\frac{\omega}{\kappa(U_H-U)}

is the trans-Planckian precursor of a fixed late frequency. The leading result is robust under a class of ultraviolet modifications only when additional adiabaticity and state assumptions hold; it is not an unconditional proof of ultraviolet independence.

Backreaction is absent. The calculation is reliable over intervals short compared with order-one change of MM and does not derive an evaporation endpoint or information recovery. It describes field propagation on the specified collapse geometry.

Two rays arrive separated by Δu\Delta u at late time. Using the exponential map, find the ratio of their early affine separations when the second pair is shifted later by an additional interval δu\delta u.

Solution

For a small fixed arrival separation, the corresponding early separation is approximately ΔU=(dU/du)Δu=Aκe−κuΔu\Delta U=(dU/du)\Delta u=A\kappa e^{-\kappa u}\Delta u. Shifting both arrivals later by δu\delta u multiplies this separation by e−κδue^{-\kappa\delta u}. Therefore equal detector-time bins sample exponentially finer early affine intervals. This calculation checks the redshift directly and explains why a logarithmic phase, rather than the detailed trajectory of the matter surface, controls the late spectrum.

  • Hawking, Stephen W. “Black Hole Explosions?” Nature 248 (1974): 30–31. DOI.
  • Hawking, Stephen W. “Erratum: Particle Creation by Black Holes.” Communications in Mathematical Physics 46 (1976): 206. DOI.
  • Hawking, Stephen W. “Particle Creation by Black Holes.” Communications in Mathematical Physics 43 (1975): 199–220. DOI.