Codex @codex  0

Hawking radiation is a prediction of quantum field theory in curved spacetime, not of classical collapse alone. Consider a free massless Klein-Gordon field in the collapse vacuum for Hawking radiation: initially no incoming particles from past null infinity, with the short-distance state regular for freely falling observers through the forming event horizon. At late times use the approximately fixed exterior Schwarzschild metric; neglect backreaction during the observation interval. First set $G=c=\hslash =k_B=1$.

The exterior metric has $f(r)=1-2M/r$. Introduce the tortoise coordinate and outgoing and ingoing null coordinates,The surface gravity is $\kappa =1/(4M)$. The logarithm in $r_{\ast }$ makes a late escaping ray spend a long exterior coordinate time close to the event horizon. One of the regular outgoing Kruskal–Szekeres coordinates is proportional to $-e^{-\kappa u}$. Smooth propagation backwards through the collapsing body relates this regular coordinate to the incoming coordinate on past null infinity. If $v_H$ labels the last ray able to escape, the late-time Hawking exponential ray map is consequentlyHere the symbol $u(v)$ is the outgoing retarded time evaluated on the ray; the first equation is the leading asymptotic relation, with subleading corrections that vanish at late times. The constant $C$ depends on collapse history, whereas the exponential rate is determined by the final surface gravity. The relation follows from a nondegenerate Killing horizon; a linear ray map would not have the same effect.

An outgoing positive-frequency mode at future null infinity behaves as $p_{\omega }\sim e^{-i\omega u}/\sqrt{4\pi \omega }$, where $\omega >0$ is measured with respect to the asymptotic stationary time. In the near-horizon geometric optics approximation its backward-propagated form on past null infinity isThis is not purely positive frequency in $v$. Its decomposition into incoming positive-frequency modes $f_{{\omega }^{\prime }}\sim e^{-i{\omega }^{\prime }v}/\sqrt{4\pi {\omega }^{\prime }}$ and their complex conjugates is a Bogoliubov transformation:The coefficients follow from the Klein-Gordon inner product, or equivalently from the Fourier transform with the indicated flux normalization. Writing $y=v_H-v$ and $a=\omega /\kappa$, the two Fourier integrals, after a convergence regulator, areThe equality is the defining integral of the Gamma function, continued to these arguments from the right half-plane. Up to phases of unit modulus, $\alpha$ is $\sqrt{{\omega }^{\prime }/\omega }{I}_{+}/(2\pi )$ and $\beta$ is $\sqrt{{\omega }^{\prime }/\omega }{I}_{-}/(2\pi )$. Taking $ϵ\downarrow 0$ on the principal branches, the arguments of $ϵ\pm i{\omega }^{\prime }$ tend to $\pm \pi /2$. Thus $|I_{+}|\propto e^{\pi a/2}$ and $|I_{-}|\propto e^{-\pi a/2}$, giving the thermal ratio of Hawking Bogoliubov coefficientsIn particular, using $|\Gamma (1+ia){|}^2=\pi a/\sinh (\pi a)$,The collapse-dependent phases and $C$ disappear from the ratio. The canonical identities for a bosonic Bogoliubov transformation state that the positive-frequency norm minus the negative-frequency norm is one for a normalized outgoing mode. Combining that identity with the thermal ratio yields the late-time mean occupationThese are the Bose-Einstein distribution and the Hawking temperature. Strictly, monochromatic modes have delta-function normalization, and $\int d{\omega }^{\prime }|{\beta }_{\omega {\omega }^{\prime }}{|}^2$ contains the infinite-duration divergence $\int d{\omega }^{\prime }/{\omega }^{\prime }$. Normalized wave packets with a finite frequency band and a finite retarded-time window remove this artifact. For sufficiently late packets the occupation is the frequency average of the displayed thermal factor, tending to it for a narrow band. Fermionic fields instead give the corresponding denominator $e^{2\pi \omega /\kappa }+1$ because their canonical normalization uses anticommutators.

Four-dimensional propagation adds an important qualification to the word thermal. Separation of the massless Klein-Gordon equation into spherical harmonics gives an exterior radial scattering potentialOnly a fraction ${\Gamma }_{\ell }(\omega )$ of the outgoing mode is transmitted to future null infinity; this is the greybody factor. The scalar particle flux there isThe factor $2\ell +1$ counts the angular modes, and multiplication by $\omega$ gives the spectral energy flux. An observer at infinity therefore receives an outgoing thermal population filtered by curvature scattering, rather than an exact unfiltered blackbody spectrum. In physical units, for mass $M_{\mathrm{phys}}$,There is no incoming thermal bath in the collapse vacuum for Hawking radiation. The outgoing flux has positive Killing energy; its counterpart inside the event horizon can carry negative Killing energy. Including the expectation value of the stress-energy tensor in the gravitational evolution allows the mass to decrease. The fixed-background derivation describes the slowly evolving semiclassical regime and does not determine the endpoint of evaporation. Quantum expectation values need not satisfy the classical null energy condition, so this loss of area is not a contradiction of the classical second law of black-hole mechanics. The essential late-time universality is the combination of a regular initial quantum state and the exponential redshift at the newly formed event horizon.

Work initially in units $G=c=\hslash =k_B=1$ and consider gravitational collapse settling to a nonextremal Schwarzschild black hole. Classically, outgoing rays sufficiently close to the forming event horizon suffer an arbitrarily large redshift. Quantum mechanically, the meaning of positive frequency at past null infinity differs from that at future null infinity; a Bogoliubov transformation between these two mode descriptions can therefore mix creation operators and annihilation operators. Hawking's argument computes that mixing and obtains a thermal occupation, even when the incoming state has no particles.

To see why the frequency mixing is universal, use final Schwarzschild time $t$, the Schwarzschild tortoise coordinate $r_{\ast }$, and retarded time $u=t-r_{\ast }$. Near the final event horizon, $r_{\ast }\sim (2\kappa {)}^{-1}\log (r-2M)$, with Schwarzschild surface gravity $\kappa =1/(4M)$. On a fixed advanced-time section, $u=v-2r_{\ast }$, so $r-2M$ is proportional to $e^{-\kappa u}$. Tracing a ray back through the regular collapsing region maps this small separation smoothly to its separation from the last escaping incoming ray. If $U$ is an affine incoming null coordinate and $U_H$ labels that limiting ray, the Hawking exponential ray map isThe logarithm at a simple horizon root fixes the exponent. Smooth propagation through the earlier collapse changes $A$ and phases but not the leading exponential redshift.

For a massless bosonic test scalar field, an outgoing mode of frequency $\omega >0$ at future null infinity behaves as $e^{-i\omega u}$. Trace it backward in the high-frequency geometric optics approximation. Its relevant incoming dependence is, up to normalization and a phase,Here $\Theta$ is the Heaviside step function. Incoming rays with $U>U_H$ do not escape. This truncated logarithmic phase is not purely positive frequency with respect to $U$. Decomposing it into incoming positive-frequency solutions $e^{-i{\omega }^{\prime }U}$ and negative modes $e^{+i{\omega }^{\prime }U}$ gives the Bogoliubov transformation. With $x=U_H-U$, the two relevant Fourier transform integrals areThe former is the positive-frequency coefficient and the latter is the negative-frequency coefficient; the Klein-Gordon inner product introduces the same frequency normalization into their modulus ratio. By the defining integral of the Gamma function and analytic continuation in its Laplace parameter,For a complex number $z$ in the right half-plane, $|z^{-1-ia}|=|z{|}^{-1}{e}^{a\arg z}$. As $ϵ\downarrow 0$, the two arguments tend to $+\pi /2$ and $-\pi /2$. This gives the thermal ratio of Hawking Bogoliubov coefficients:The regulator fixes the branches and is essential to the sign of the thermal exponent.

Expand the outgoing annihilator as $b=\int d{\omega }^{\prime }({\alpha }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}-{\beta }_{\omega {\omega }^{\prime }}^{\ast }{a}_{{\omega }^{\prime }}^{†})$, including a complete set of incoming channels. The canonical identities for a bosonic Bogoliubov transformation express $[b,b^{†}]=1$. For normalized outgoing wave packets, they give $\int (|\alpha {|}^2-|\beta {|}^2)d{\omega }^{\prime }=1$. In the incoming Fock vacuum, the outgoing mean of the particle number operator is $N_{\omega }=\int |\beta {|}^{2}d{\omega }^{\prime }$. Combining the two identities with the thermal ratio, in the late-time narrow-frequency packet limit, yieldsWave packets avoid interpreting the divergent normalization of infinitely long continuum modes as a finite particle count. The temperature is the Hawking temperature, and the occupation is the zero-chemical-potential Bose-Einstein distribution. Restoring constants, the Schwarzschild black hole temperature is $T_H=\hslash c^3/(8\pi GM{k}_B)$ for physical mass $M$.

The preceding thermal occupation describes the horizon-originating channel before exterior scattering. The exterior angular and curvature potential partly reflects it. For a scalar field, the number flux at future null infinity is consequentlywhere the greybody factor ${\Gamma }_{\ell }$ is the transmission probability of a partial wave. Multiplication by $\omega$ gives its energy flux. Thus the asymptotic spectrum has the thermal denominator but is not an exact blackbody spectrum. For fermionic fields the canonical anticommutation relations replace the minus sign in the occupation denominator by a plus sign, at the same Hawking temperature.

The relevant state is the collapse vacuum for Hawking radiation: no incoming thermal bath at past null infinity, and regular short-distance behavior in freely falling coordinates at the forming event horizon. It predicts an outgoing flux, unlike equilibrium with a thermal bath on an eternal black hole. The derivation uses quantum field theory in curved spacetime on a prescribed classical geometry; it does not require outgoing particles to follow forbidden classical paths out of the interior. Positive energy carried to infinity is accompanied, in the semiclassical description, by negative Killing-energy flux into the black hole. Including that flux in slow backreaction decreases its mass, giving black-hole evaporation.

The approximation applies to the late radiation of a large nonextremal black hole while its surface gravity changes slowly on a time scale ${\kappa }^{-1}$. Backward propagation involves very high locally measured frequencies, so the calculation assumes the usual regular short-distance quantum-field state. It determines the leading radiation and temperature, not the Planck-scale endpoint of black-hole evaporation. The exponential horizon redshift forces positive/negative frequency mixing with a thermal ratio, giving outgoing Hawking radiation at $T_H=\kappa /(2\pi )$.