Hawking radiation: Bogoliubov derivation, thermal spectrum, and black-hole evaporation

Intuition Beginner

Black holes are not completely black. In 1974 Stephen Hawking showed that when one quantises ordinary matter fields on the curved background of a black hole, the hole radiates particles at a precise temperature set by its surface gravity. The result, announced in Nature 248 (1974) 30 and developed in detail in Comm. Math. Phys. 43 (1975) 199, is one of the deepest predictions of theoretical physics and the only known quantum-mechanical leak from an object that classical general relativity forbids any signal to escape.

The physical picture often used in popular accounts goes like this. The quantum vacuum is restless: pairs of virtual particles flicker into existence and back out everywhere, even in empty space. Near a black-hole horizon these pairs can be torn apart. One member of the pair falls into the hole carrying negative energy as measured at infinity; the partner escapes outward, carrying positive energy. To a distant observer the hole appears to emit radiation, and to conserve total energy the hole loses mass at the rate set by the outgoing flux. The actual technical derivation is a calculation in quantum field theory on a collapsing-star background, but the heuristic captures the essential bookkeeping.

The temperature of the radiation is determined by the mass of the hole. In conventional units it is $T_H=\hslash c^3/(8\pi GM{k}_B)$. The Hawking temperature scales as one over the mass: small holes are hot, large holes are cold. A solar-mass black hole has $T_H\approx 6\times 10^{-8}$ K, vastly colder than the $2.7$ K cosmic microwave background — such a hole gains more mass from absorbing the CMB than it loses to Hawking radiation, and so is not observably evaporating today.

A primordial black hole formed in the very early universe with mass around $10^{14}$ to $10^{15}$ g would have a Hawking temperature of $10^{11}$ K and would be evaporating now in a final burst of gamma rays. The non-detection of such bursts in dedicated searches like Fermi GBM and HAWC places upper limits on the cosmological abundance of primordial black holes in this mass window. A hypothetical micro-black hole of mass $1$ kg would have Hawking temperature about $10^{23}$ K and would explode in roughly $10^{-22}$ seconds, releasing about $10^{17}$ J — about a megaton of TNT.

The lifetime of an evaporating black hole follows from the Stefan-Boltzmann law applied to a blackbody of horizon area $A=16\pi (GM/c^2{)}^2$ at temperature $T_H$. The result is $\tau =(5120\pi G^2/\hslash c^4)M^3$, scaling as the cube of the mass. A solar-mass black hole has Hawking lifetime around $10^{67}$ years — far longer than the current age of the universe ($1.4\times 10^{10}$ years) but finite. The eventual fate of every isolated black hole in our universe is to evaporate completely, given enough cosmic time.

Visual Beginner

The picture below sketches the Penrose-style causal diagram of a collapsing star together with the Hawking radiation streaming out to future null infinity. The collapsing matter is represented by a shaded region; the event horizon emerges as the boundary of the trapped region; the Hawking quanta are represented by wavy arrows fanning out from just outside the horizon toward future null infinity.

The picture encodes three pieces. First, the in-state of the quantised field, prepared on past null infinity before the star collapses, is the standard Minkowski vacuum: positive-frequency modes are those that look like plane waves with positive frequency at past null infinity. Second, the collapsing geometry has no time-translation symmetry during the collapse — the metric depends on time — so the notion of positive frequency mixes between past and future. Third, the asymptotic out-state at future null infinity, when expressed in terms of in-modes, is not the in-vacuum: it is a thermal state of out-particles at the Hawking temperature, with a Planck distribution governing the occupation of each out-mode.

Worked example Beginner

Take a primordial black hole with the mass that would make its Hawking lifetime equal to the current age of the universe. Compute its initial Hawking temperature, the typical energy of the emitted radiation, and the total energy that has already been released into the universe by such holes through their first cosmic-time evaporation.

Step 1. Mass for lifetime equal to age of universe. Solving $\tau =(5120\pi G^2/\hslash c^4)M^3=4.4\times 10^{17}$ s for $M$, with $G=6.67\times 10^{-11}$, $\hslash =1.05\times 10^{-34}$, $c=3\times 10^8$:

$$M^3=\frac{\hslash c^4\tau }{5120\pi G^2}=\frac{(1.05\times 10^{-34})(3\times 10^8{)}^4(4.4\times 10^{17})}{5120\pi (6.67\times 10^{-11}{)}^2}\approx 1.7\times 10^{33}{\text{ kg}}^3,$$

so $M\approx 1.2\times 10^{11}$ kg, equivalent to $1.2\times 10^{14}$ g.

Step 2. Initial Hawking temperature. Plug $M=1.2\times 10^{11}$ kg into $T_H=\hslash c^3/(8\pi GM{k}_B)$:

$$T_H=\frac{(1.05\times 10^{-34})(3\times 10^8{)}^3}{8\pi (6.67\times 10^{-11})(1.2\times 10^{11})(1.38\times 10^{-23})}\approx 10^{12}\text{ K}.$$

The typical energy per emitted quantum is $E\approx k_B{T}_H\approx 10^{-11}$ J $\approx 100$ MeV — gamma-ray photons, plus electron-positron pairs, plus higher-mass particles as the hole heats up toward the final burst.

Step 3. Final-flash energy. A black hole evaporating today has nearly its full initial mass-energy remaining in the late stages of the burn. With $M=1.2\times 10^{11}$ kg, the rest-mass energy is $M{c}^2=(1.2\times 10^{11})(9\times 10^{16})\approx 10^{28}$ J. Each primordial black hole at this mass scale completes its evaporation by releasing about $10^{28}$ J of energy as a localised burst of gamma rays. Searches for such bursts by the Fermi Gamma-ray Burst Monitor and the HAWC observatory bound the local density of evaporating primordial black holes to fewer than about $10^{-9}$ per cubic parsec per year.

What this tells us: a black hole of mass $10^{14}$ g evaporating now would emit gamma rays with energies in the 100 MeV range and release $10^{28}$ J in a final burst; non-observation of such bursts excludes primordial black holes from being a dominant component of dark matter in this mass window.

Check your understanding Beginner

Formal definition Intermediate+

We work in geometric units $G=c=\hslash =k_B=1$ where convenient, restoring physical units when computing observable numbers. Consider a free massless real scalar field $\varphi$ on a globally hyperbolic Lorentzian background $(M,g)$ obeying the Klein-Gordon equation

$${\square }_g\varphi =g^{ab}{\nabla }_a{\nabla }_b\varphi =0.$$

The space of classical solutions carries a natural symplectic form

$$\Omega ({\varphi }_1,{\varphi }_2)={\int }_{\Sigma }({\varphi }_1{\nabla }_a{\varphi }_2-{\varphi }_2{\nabla }_a{\varphi }_1)n^{a}d\Sigma ,$$

with $\Sigma$ any Cauchy surface and $n^a$ its future-directed unit normal. A positive-frequency decomposition of the solution space is a choice of complex structure $J$ on the real symplectic vector space, equivalent to a choice of complete orthonormal mode basis $\{u_k\}$ in the Klein-Gordon inner product $(u_k,u_{k^{\prime }})=-i\Omega ({\bar{u}}_k,u_{k^{\prime }})={\delta }_{k{k}^{\prime }}$, picking out which modes count as positive-frequency.

Definition (in/out vacua). Given a Cauchy surface in the asymptotic past on which the geometry is stationary (e.g., past null infinity ${\mathcal{I}}^{-}$ of an asymptotically flat spacetime), the in-modes $\{u_k^{\text{in}}\}$ are the positive-frequency solutions with respect to the past asymptotic time-translation Killing field, normalised as above. The in-vacuum $|0_{\text{in}}⟩$ is the state annihilated by the in-mode annihilators ${\hat{a}}_k=(u_k^{\text{in}},\hat{\varphi }{)}_{\text{KG}}$. The out-modes $\{u_k^{\text{out}}\}$ and out-vacuum $|0_{\text{out}}⟩$ are defined analogously with respect to the future asymptotic Killing field.

Definition (Bogoliubov coefficients). Since both $\{u_k^{\text{in}}\}$ and $\{u_k^{\text{out}}\}$ are complete bases of the same complexified solution space, the second is a linear combination of the first plus complex conjugates:

$$u_k^{\text{out}}=\sum _p\left({\alpha }_{kp}{u}_p^{\text{in}}+{\beta }_{kp}{\bar{u}}_p^{\text{in}}\right).$$

The coefficients $(\alpha ,\beta )$ are the Bogoliubov coefficients of the transformation. By Klein-Gordon orthonormality of both bases they satisfy

$$\sum _p({\alpha }_{kp}{\bar{\alpha }}_{k^{\prime }p}-{\beta }_{kp}{\bar{\beta }}_{k^{\prime }p})={\delta }_{k{k}^{\prime }},\sum _p({\alpha }_{kp}{\beta }_{k^{\prime }p}-{\beta }_{kp}{\alpha }_{k^{\prime }p})=0.$$

The induced relation between mode operators is ${\hat{b}}_k={\sum }_p({\bar{\alpha }}_{kp}{\hat{a}}_p-{\bar{\beta }}_{kp}{\hat{a}}_p^{†})$, where ${\hat{b}}_k=(u_k^{\text{out}},\hat{\varphi }{)}_{\text{KG}}$ are out-mode annihilators.

Definition (Hawking effect). Take the asymptotically flat Schwarzschild spacetime that arises by gravitational collapse of a regular spherically symmetric matter distribution. The Hawking effect is the statement that the expected number of out-particles in mode $k$ measured by an observer at future null infinity in the in-vacuum $|0_{\text{in}}⟩$ is

$$⟨0_{\text{in}}|{\hat{N}}_k^{\text{out}}|0_{\text{in}}⟩=\sum _p|{\beta }_{kp}{|}^2=\frac{{\Gamma }_{k\ell }}{e^{\omega /T_H}-1},$$

a Planck distribution at the Hawking temperature $T_H=\kappa /(2\pi )=1/(8\pi M)$ in geometric units (equivalently $T_H=\hslash c^3/(8\pi GM{k}_B)$ in conventional units), modulated by a greybody factor ${\Gamma }_{k\ell }$ that accounts for partial-wave backscattering off the Schwarzschild gravitational potential between the horizon and infinity.

Counterexamples to common slips

• The Hawking effect is not a particle-antiparticle pair production at the horizon. The popular heuristic of virtual pairs split by the horizon is not the actual derivation; in Hawking's calculation the radiation arises from the propagation of vacuum fluctuations through the gravitational-collapse geometry, with the relevant mathematical object being the Bogoliubov transformation between in-modes on ${\mathcal{I}}^{-}$ and out-modes on ${\mathcal{I}}^{+}$, not a localised pair-creation event near the horizon. The pair picture is a useful mnemonic only.
• The Hawking temperature does not equal the gravitational redshift. For an observer just outside the horizon, the locally measured temperature of the radiation diverges (because every photon emitted from the horizon is enormously blueshifted in the static frame); the Hawking temperature is this diverging local temperature redshifted to infinity, finite by definition.
• The thermal spectrum is not a Boltzmann distribution. It is a Planck distribution, $1/(e^{\omega /T_H}-1)$, because the field is bosonic. For a fermionic field one gets a Fermi-Dirac distribution $1/(e^{\omega /T_H}+1)$ with the same temperature.

Key theorem with proof Intermediate+

Theorem (Bogoliubov derivation of the Hawking spectrum on a collapsing-shell background). Let $\varphi$ be a free massless scalar field on the spherically symmetric spacetime obtained by gravitational collapse of a thin null shell of total ADM mass $M$. Prepare $\varphi$ in the Minkowski in-vacuum on ${\mathcal{I}}^{-}$. Then the asymptotic out-flux at ${\mathcal{I}}^{+}$ for each spherical-harmonic partial wave $\ell$ and outgoing frequency $\omega$ is a Planck distribution at temperature $T_H=\kappa /(2\pi )=1/(8\pi M)$, modulated by a greybody factor ${\Gamma }_{\omega \ell }$ that depends on the partial-wave transmission probability through the Schwarzschild potential:

$$\frac{d{N}_{\omega \ell }}{dtd\omega }=\frac{1}{2\pi }\frac{{\Gamma }_{\omega \ell }}{e^{\omega /T_H}-1}.$$

Proof. The geometry is Minkowski for $v<v_0$ and Schwarzschild for $v>v_0$ where $v$ is the advanced Eddington-Finkelstein coordinate, $v_0$ the moment the collapsing null shell crosses, and the two regions are matched along the shell. Outside the shell the metric is Schwarzschild $d{s}^2=-(1-2M/r)d{t}^2+(1-2M/r{)}^{-1}d{r}^2+r^{2}d{\Omega }^2$.

Take an outgoing s-wave mode $u_{\omega }^{\text{out}}$ at ${\mathcal{I}}^{+}$ of asymptotic form $u_{\omega }^{\text{out}}\sim e^{-i\omega u}/(r\sqrt{4\pi \omega })$, where $u=t-r^{\ast }$ is the retarded null coordinate and $r^{\ast }=r+2M\ln (r/2M-1)$ is the tortoise coordinate. Trace this mode backwards in time through the collapse geometry to find its form on ${\mathcal{I}}^{-}$.

Inside the collapsing shell the geometry is flat Minkowski, so the mode propagates as a free spherical wave; outside the shell it propagates on Schwarzschild. The mode equation in the Schwarzschild region reduces, after separation of variables, to the s-wave radial equation

$$\frac{d^2{\psi }_{\omega }}{d{r}^{\ast 2}}+{\omega }^2{\psi }_{\omega }=V(r){\psi }_{\omega },$$

with the Regge-Wheeler potential $V(r)=(1-2M/r)\cdot 2M/r^3$ for s-waves. A geometric-optics (high-frequency) treatment of the late-time outgoing modes is sufficient to derive the Hawking spectrum: high-frequency outgoing modes near ${\mathcal{I}}^{+}$ trace back to modes near the horizon with arbitrarily high local frequency, then propagate inward across the collapsing shell.

The key step is the relation between the affine parameter on a horizon-skimming null geodesic and the retarded time $u$ at ${\mathcal{I}}^{+}$. A null geodesic that originates at ${\mathcal{I}}^{-}$ with advanced-time coordinate $v$ just below the horizon-formation value $v_H$ propagates outward through the shell and arrives at ${\mathcal{I}}^{+}$ at retarded time

$$u\approx -\frac{1}{\kappa }\ln (v_H-v)+\text{const},$$

with $\kappa =1/(4M)$ the surface gravity. The logarithm comes from the exponential redshift of modes propagating outward near a horizon-generating null geodesic: a fixed change in retarded time $u$ corresponds to an exponentially small change in advanced time $v$.

Tracing the out-mode $u_{\omega }^{\text{out}}\sim e^{-i\omega u}$ back through this relation gives an in-mode of the form

$$u_{\omega }^{\text{in}}(v)\sim e^{-i\omega u(v)}=e^{(i\omega /\kappa )\ln (v_H-v)}\cdot \text{phase}=(v_H-v{)}^{i\omega /\kappa }\cdot \text{phase}$$

for $v<v_H$, and zero for $v>v_H$. This is the form of the late-time outgoing mode evaluated at ${\mathcal{I}}^{-}$ — a power-law-modulated wave with a sharp cutoff at $v=v_H$.

Fourier-decomposing this $v$-dependence into positive- and negative-frequency components with respect to the Minkowski past-asymptotic Killing field ${\partial }_v$ at ${\mathcal{I}}^{-}$ yields the Bogoliubov coefficients. The standard contour-integration argument (Hawking 1975 §2.5, also Birrell-Davies 1982 §8.1) evaluates

$${\alpha }_{\omega {\omega }^{\prime }}\sim {\int }_{-\infty }^{v_H}(v_H-v{)}^{i\omega /\kappa }{e}^{i{\omega }^{\prime }v}dv,{\beta }_{\omega {\omega }^{\prime }}\sim {\int }_{-\infty }^{v_H}(v_H-v{)}^{i\omega /\kappa }{e}^{-i{\omega }^{\prime }v}dv,$$

with the second integral obtained from the first by analytic continuation of ${\omega }^{\prime }$ through the upper half-plane, picking up the rotation phase $e^{-\pi \omega /\kappa }$ in the process. The ratio of the squared magnitudes is

$$\frac{|{\beta }_{\omega {\omega }^{\prime }}{|}^2}{|{\alpha }_{\omega {\omega }^{\prime }}{|}^2}=e^{-2\pi \omega /\kappa }.$$

Combining with the Bogoliubov normalisation ${\sum }_{{\omega }^{\prime }}(|{\alpha }_{\omega {\omega }^{\prime }}{|}^2-|{\beta }_{\omega {\omega }^{\prime }}{|}^2)=1$ gives

$$\sum _{{\omega }^{\prime }}|{\beta }_{\omega {\omega }^{\prime }}{|}^2=\frac{1}{e^{2\pi \omega /\kappa }-1}=\frac{1}{e^{\omega /T_H}-1},$$

with $T_H=\kappa /(2\pi )$. Including the partial-wave label $\ell$ and the greybody factor ${\Gamma }_{\omega \ell }$ that accounts for backscattering off the Regge-Wheeler potential gives the stated emission rate. $\square$

Bridge. The factor $e^{-2\pi \omega /\kappa }$ identifies the in-vacuum at ${\mathcal{I}}^{+}$ with a thermal state at temperature $T_H=\kappa /(2\pi )$, the foundational reason being the exponential redshift of horizon-skimming null modes; this is exactly the universal kinematic statement that any bifurcate Killing horizon carries a thermal state at temperature equal to the surface gravity divided by $2\pi$. The Bogoliubov calculation builds toward 13.06.05 pending (Page curve and the information paradox), where the time-dependent von Neumann entropy of the outgoing radiation is computed and the unitarity question becomes sharp. It also appears again in 13.09.08 (Bunch-Davies state on de Sitter), where the same surface-gravity-to-temperature dictionary gives the Gibbons-Hawking temperature $T_{dS}=H/(2\pi )$ of the de Sitter cosmological horizon; the bridge is between local horizon kinematics and global thermal physics, and the central insight is that the spectrum is fixed by the geometry alone, independent of the matter content of the collapse. The same pattern recurs in the Unruh effect, where the uniformly accelerated observer's Rindler horizon gives a thermal Rindler bath at $T_U=a/(2\pi )$, identifying surface gravity with proper acceleration via the equivalence principle.

Exercises Intermediate+

Advanced results Master

The 1974-75 Hawking calculations promoted black holes from classical gravitational solutions to thermodynamic objects with a precise quantum-mechanical temperature, and ignited the modern programme of black-hole information theory. We collect the central results.

Theorem 1 (Hawking effect; Hawking 1974, Nature 248, 30; 1975, Comm. Math. Phys. 43, 199). Let $\varphi$ be a free massless scalar field on the spacetime obtained by spherically symmetric gravitational collapse to a Schwarzschild black hole of mass $M$. Prepare $\varphi$ in the in-vacuum on ${\mathcal{I}}^{-}$. Then the asymptotic flux at ${\mathcal{I}}^{+}$ is a thermal flux at temperature $T_H=\hslash \kappa /(2\pi k_{B}c)=\hslash c^3/(8\pi GM{k}_B)$ in conventional units, modulated by partial-wave-dependent greybody factors. The result extends to fields of arbitrary spin, with Bose-Einstein occupation $1/(e^{\omega /T_H}-1)$ for integer-spin fields and Fermi-Dirac occupation $1/(e^{\omega /T_H}+1)$ for half-integer-spin fields, the same temperature in both cases.

The proof, sketched in the Key Theorem section above, traces high-frequency outgoing s-wave modes backward through the gravitational-collapse geometry; the logarithmic redshift relation between affine parameters on horizon-skimming null geodesics and retarded time at infinity gives the exponential redshift $e^{-\pi \omega /\kappa }$ in the ratio of Bogoliubov $\beta$-coefficients to $\alpha$-coefficients, which combined with normalisation yields the Planck distribution. The derivation is robust to the details of the collapsing matter: a thin null shell, a dust cloud, and a regular fluid collapse all give the same asymptotic out-spectrum. This universality is the first hint that the Hawking effect is a property of the horizon geometry alone, not of the collapse history.

Theorem 2 (Euclidean / conical-deficit derivation; Gibbons-Hawking 1977, Phys. Rev. D 15, 2752). The Wick-rotated Schwarzschild metric in imaginary time, $t\to -i\tau$, defines a positive-definite four-manifold whose geometry near $r=2M$ is regular only when the imaginary-time coordinate $\tau$ has period $\beta =8\pi GM/\hslash$. The Hawking temperature is then identified as $T_H=1/\beta$ by the standard thermofield-double / KMS argument: the gravitational path integral with this periodic Euclidean time computes the thermal partition function $Z=\mathrm{Tr}{e}^{-\beta H}$, and stationary observers in the Lorentzian section see a thermal state at this temperature.

The conical-deficit argument is the cleanest existing derivation of $T_H$ and makes manifest that the temperature is determined by the horizon kinematics alone. Restoring the calculation in physical units: the saddle-point evaluation of the Euclidean gravitational action gives $\beta F=-{\beta }^2/(16\pi G\hslash )+(\text{boundary terms})$, from which $F=-\beta /(16\pi G\hslash )+M=M/2$ (with $\beta =8\pi GM/\hslash$), giving the canonical-ensemble first law $dF=-SdT$ that reproduces the Bekenstein-Hawking entropy $S_{BH}=A/(4G\hslash )$. The Gibbons-Hawking 1977 paper is the foundational document of Euclidean quantum gravity in its black-hole application.

Theorem 3 (Parikh-Wilczek tunnelling derivation; Parikh-Wilczek 2000, Phys. Rev. Lett. 85, 5042). Hawking radiation can be derived as a semiclassical WKB tunnelling process: a Hawking quantum of energy $\omega$ tunnels through the horizon along a classically forbidden trajectory with action $\Delta S_{BH}/\hslash$, where $\Delta S_{BH}$ is the change in the Bekenstein-Hawking entropy. To leading order in $\omega /M$ the tunnelling probability $\exp (-\omega /T_H)$ reproduces the Boltzmann factor of the Planck distribution at $T_H$; the next-order correction $+4\pi {\omega }^2/\hslash$ encodes back-reaction of the emitted quantum on the black-hole mass.

The Parikh-Wilczek derivation makes the connection between Hawking radiation and the second law of black-hole thermodynamics maximally transparent: each emitted Hawking quantum is a unit of entropy flux through the horizon, with the emission probability set by the entropy cost of removing energy from the hole. This perspective is the conceptual seed of the modern programme that relates Hawking radiation to entanglement and to the gravitational path integral on replica geometries.

Theorem 4 (Greybody factors; Page 1976, Phys. Rev. D 13, 198; Page 1976, Phys. Rev. D 14, 3260). The asymptotic out-spectrum at ${\mathcal{I}}^{+}$ is modulated by partial-wave-dependent greybody factors ${\Gamma }_{s\ell \omega }$ that quantify the transmission probability for a wave of spin $s$, partial-wave $\ell$, and frequency $\omega$ to propagate from just outside the horizon through the angular-momentum barrier to infinity. For a Schwarzschild black hole, ${\Gamma }_{0\ell \omega }$ for s-waves at low frequency scales as $(2M\omega {)}^{2\ell +2}$, peaks near $2M\omega \sim 1$, and approaches unity at high frequency.

Page 1976 computed the spin-0, spin-1/2, spin-1, and spin-2 greybody factors and integrated to get the total luminosity and the species-dependent fractional emission rate. For an uncharged non-rotating black hole the radiation is dominated by photons and gravitons at high mass, by photons + electrons + positrons at intermediate mass, and by the full Standard Model spectrum (including all light hadrons and the QCD-confinement-induced jet of pions) once the Hawking temperature exceeds about $100$ MeV. The corrections to the Stefan-Boltzmann estimate of the lifetime from greybody factors and species multiplicities are factors of order $10$-$100$ but do not change the qualitative $M^3$ scaling.

Theorem 5 (Hartle-Hawking state; Hartle-Hawking 1976, Phys. Rev. D 13, 2188). On the maximally extended (Kruskal) Schwarzschild manifold, there is a unique Hadamard quasi-free state of the free Klein-Gordon field that is invariant under the bifurcate-Killing isometry generated by the time-translation Killing vector, regular on both the future and past horizons, and thermal at temperature $T_H$ to any static observer. This is the Hartle-Hawking state $|0_{HH}⟩$.

The Hartle-Hawking state is the natural state for the eternal black hole (a Schwarzschild geometry without a collapse history); it is the state in which the black hole is in thermal equilibrium with a bath of Hawking radiation at infinity. The Israel 1976 thermofield-double construction identifies the Hartle-Hawking state with the maximally entangled purification of the canonical thermal state on the exterior region, with the second exterior region (region IV of the Kruskal extension) playing the role of the thermofield partner. The Hartle-Hawking state on the Kruskal extension is the black-hole analogue of the Bunch-Davies state on de Sitter (see 13.09.08); both arise from the requirement that the Euclidean section be smooth at the horizon, and both are the unique Hadamard states respecting the maximal symmetry of the spacetime.

Theorem 6 (Unruh state; Unruh 1976, Phys. Rev. D 14, 870). For the spacetime of a black hole that forms by gravitational collapse (as opposed to the eternal-black-hole spacetime above), the appropriate Hadamard state is the Unruh state $|0_U⟩$, defined to be the in-vacuum on ${\mathcal{I}}^{-}$ and the Hawking state thermal at $T_H$ on the future horizon. The Unruh state is not thermal on the past horizon; it differs from the Hartle-Hawking state precisely on that asymptotic region.

The three canonical states on the Schwarzschild geometry — Boulware, Unruh, Hartle-Hawking — form a natural hierarchy. The Boulware state $|0_B⟩$ is the state regular at ${\mathcal{I}}^{-}$ and ${\mathcal{I}}^{+}$ but singular on both horizons; it is the state appropriate for a quasi-static observer outside a star whose surface is well above the would-be horizon radius. The Unruh state is appropriate for a black hole formed by collapse. The Hartle-Hawking state is appropriate for the eternal Kruskal-extended black hole. This taxonomy (Wald 1994, Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics, §5) is the canonical framework for QFT in a black-hole background.

Theorem 7 (Unruh effect and equivalence-principle identification; Davies 1975, J. Phys. A 8, 609; Unruh 1976, Phys. Rev. D 14, 870). A uniformly accelerated observer in Minkowski space with proper acceleration $a$ perceives the Minkowski vacuum as a thermal bath at the Unruh temperature $T_U=\hslash a/(2\pi c{k}_B)$. The Unruh temperature has the same surface-gravity-to-temperature dictionary as the Hawking temperature, $T={\kappa }_{\text{eff}}/(2\pi )$, with the proper acceleration playing the role of an effective surface gravity for the Rindler horizon of the accelerated observer.

The Unruh effect is the equivalence-principle counterpart of the Hawking effect: a static observer just outside a Schwarzschild horizon has a locally measured proper acceleration $a=\kappa /\sqrt{1-2M/r}$, and the locally measured temperature in the Unruh state agrees with the Tolman-redshifted Hawking temperature $T(r)=T_H/\sqrt{1-2M/r}$. Both effects are particular cases of the general thermal-state property of bifurcate Killing horizons established by Sewell 1982 Ann. Phys. 141, 201 via the Bisognano-Wichmann theorem in algebraic quantum field theory: the modular automorphism group of the von Neumann algebra of observables in a Rindler wedge of Minkowski space is the boost subgroup of the Poincaré group, and the resulting KMS condition at inverse temperature $\beta =2\pi /a$ identifies the Minkowski vacuum as a thermal state to the accelerated observer.

Theorem 8 (Page curve and the information paradox; Page 1993, Phys. Rev. Lett. 71, 3743; AMPS 2013, JHEP 1302, 062; Penington 2020, JHEP 09, 002; AHMST 2020, JHEP 1907, 063). Suppose a black hole forms from a pure quantum state and evaporates completely into Hawking radiation, and suppose the evaporation is unitary. Then the von Neumann entropy of the radiation as a function of evaporation time follows the Page curve: it rises while the radiation is small (the radiation is approximately thermal in this regime), peaks at the Page time $t_{\text{Page}}\sim M^3$ when about half the original Bekenstein-Hawking entropy has been emitted, and then decreases back to zero as the late-time radiation becomes purified by the earlier radiation.

The information paradox of Hawking 1976 Phys. Rev. D 14, 2460 was the observation that the naive Hawking spectrum is exactly thermal at all times, giving a von Neumann entropy that grows monotonically with time and reaches $S_{BH}$ at full evaporation — incompatible with unitary evolution of an initially pure state. The AMPS 2013 firewall argument sharpened the tension: monogamy of entanglement forbids the simultaneous assumption of smooth-horizon effective field theory, unitary evaporation, and the equivalence-principle prediction that an infalling observer sees nothing unusual at the horizon. The Penington 2020 and AHMST 2020 island-formula derivations resolve the paradox at the semiclassical level: the entanglement entropy of the radiation in the gravitational path integral is computed by a generalised-entropy minimisation $S({\rho }_R)={\min }_{\sigma }{\text{ext}}_{\sigma }[A(\sigma )/(4G\hslash )+S_{\text{matter}}({\Sigma }_R\cup \mathcal{I})]$ over candidate quantum-extremal surfaces $\sigma$ bounding a possibly non-empty island region $\mathcal{I}$ in the black-hole interior; before the Page time the empty-island extremum dominates and reproduces the Hawking spectrum, while after the Page time a non-empty-island extremum dominates and brings in the interior degrees of freedom that purify the late-time radiation.

Synthesis. Putting these together, Hawking radiation is the foundational reason that classical relativity, quantum field theory, and thermodynamics fit into a single coherent semiclassical framework, with the three derivations (Bogoliubov on the collapsing-shell geometry, Euclidean conical-deficit on the maximally extended manifold, and Parikh-Wilczek tunnelling across the horizon) giving the same Hawking temperature $T_H=\kappa /(2\pi )$ from independent starting points. The central insight is that the surface gravity $\kappa$ on a bifurcate Killing horizon plays the role of a temperature in the universal sense: the in-vacuum at past null infinity is identified, by the Bogoliubov transformation that propagates modes through the collapse geometry, with a Planck-distributed thermal flux at future null infinity, and the same result follows from Euclidean smoothness alone without reference to any collapse history. This is exactly the structure that identifies horizon kinematics with thermal physics; the bridge is the universal surface-gravity-to-temperature dictionary $T=\kappa /(2\pi )$ that recurs in the Unruh effect at $T_U=a/(2\pi )$ and in the Gibbons-Hawking de Sitter case at $T_{dS}=H/(2\pi )$, and the pattern generalises to any bifurcate Killing horizon via the Bisognano-Wichmann modular-automorphism argument.

The same chain of identifications puts the area theorem, the four laws of black-hole mechanics, and the Bekenstein-Hawking entropy formula on a single mechanical footing: the Hawking effect promotes $\kappa /(2\pi )\leftrightarrow T$ and $A/4\leftrightarrow S$ from analogies to identities, with the proportionality coefficient fixed by the Bogoliubov calculation. Putting this together with the Page-curve programme of the 1993 paper and the 2020 island-formula derivations identifies the late-time radiation's entanglement structure as the gravitational-path-integral analogue of the quantum-error-correcting code of holographic entanglement-wedge reconstruction; the bridge to AdS/CFT in Maldacena 1998 Adv. Theor. Math. Phys. 2, 231 is the explicit unitary boundary description that converts the thermal Hawking flux from an apparent unitarity violation into the standard thermalisation of a strongly coupled large-$N$ gauge theory. The information paradox, in this contemporary view, is not a contradiction but a consistency condition on quantum gravity: the gravitational path integral must reproduce the Page curve from semiclassical input alone, and the island-formula derivations are the most rigorous existing evidence that it does.

Full proof set Master

We give the central proof of the Hawking spectrum via the Bogoliubov transformation between in-modes on ${\mathcal{I}}^{-}$ and out-modes on ${\mathcal{I}}^{+}$ on a gravitationally collapsing geometry, following Hawking 1975 and the codification in Birrell-Davies 1982 §8 and Wald 1994 §6.

Proposition (Hawking spectrum from horizon-skimming geodesic redshift). Let $\varphi$ be a free massless scalar field on the spacetime of a spherically symmetric gravitational collapse to a Schwarzschild black hole of mass $M$, prepared in the in-vacuum $|0_{\text{in}}⟩$ on ${\mathcal{I}}^{-}$. Let $u_{\omega }^{\text{out}}$ be a late-time s-wave outgoing mode at ${\mathcal{I}}^{+}$ of asymptotic form $u_{\omega }^{\text{out}}\sim e^{-i\omega u}/(r\sqrt{4\pi \omega })$ with $u=t-r^{\ast }$ the retarded null coordinate. Then the in-vacuum expectation of the out-number operator is

$$⟨0_{\text{in}}|{\hat{N}}_{\omega }^{\text{out}}|0_{\text{in}}⟩=\frac{{\Gamma }_{\omega }}{e^{2\pi \omega /\kappa }-1},$$

with $\kappa =1/(4M)$ the surface gravity and ${\Gamma }_{\omega }$ the s-wave greybody factor.

Proof. The proof has four steps.

Step 1: geometric-optics propagation. The late-time outgoing s-wave mode $u_{\omega }^{\text{out}}$ at ${\mathcal{I}}^{+}$ has support concentrated on retarded times $u$ that grow without bound as the affine parameter on a horizon-skimming generator advances. Propagating the mode backward in time using the geometric-optics / WKB approximation (which is exact in the high-frequency limit and dominates the late-time spectrum), the mode follows a null geodesic from ${\mathcal{I}}^{+}$ back through the Schwarzschild exterior, through the collapsing matter, and out to ${\mathcal{I}}^{-}$. In the high-frequency limit the geometric-optics propagation is exact up to backscattering off the Regge-Wheeler potential, which is absorbed into the greybody factor ${\Gamma }_{\omega }$ and does not affect the thermal-spectrum derivation.

Step 2: redshift relation between affine parameter and retarded time. Let $v_H$ denote the advanced-time value at which the future event horizon forms on ${\mathcal{I}}^{-}$ (the value of the advanced-Eddington-Finkelstein coordinate $v$ at which the horizon-generating null geodesic reaches ${\mathcal{I}}^{-}$). A horizon-skimming outgoing null geodesic that originates on ${\mathcal{I}}^{-}$ with advanced time $v<v_H$ and propagates outward arrives at ${\mathcal{I}}^{+}$ at retarded time

$$u=u_0-\frac{1}{\kappa }\ln (v_H-v)+O(v_H-v),$$

with $u_0$ a constant determined by the matching across the collapsing shell. This relation is the central kinematic fact of the Hawking effect: a fixed change in retarded time $\delta u$ at ${\mathcal{I}}^{+}$ corresponds to an exponentially small change $\delta v\propto e^{-\kappa u}(v_H-v)$ in advanced time at ${\mathcal{I}}^{-}$. The exponential redshift of horizon-skimming modes is the geometric content of the Hawking effect.

Step 3: in-mode form on ${\mathcal{I}}^{-}$. Tracing the late-time out-mode $u_{\omega }^{\text{out}}\sim e^{-i\omega u}$ back through the geometric-optics propagation gives a mode on ${\mathcal{I}}^{-}$ of the form

$$f_{\omega }(v)=e^{-i\omega u(v)}\Theta (v_H-v)=e^{i(\omega /\kappa )\ln (v_H-v)}{e}^{-i\omega u_0}\Theta (v_H-v).$$

The mode is non-zero only for $v<v_H$ (modes with $v>v_H$ fall into the black hole and do not reach ${\mathcal{I}}^{+}$) and has a power-law-modulated oscillatory dependence on $v_H-v$. Fourier-decomposing this $v$-dependence into positive- and negative-frequency Minkowski components with respect to the past asymptotic time-translation Killing field ${\partial }_v$ at ${\mathcal{I}}^{-}$ gives the Bogoliubov coefficients.

Step 4: Bogoliubov ratio via analytic continuation. Define the Fourier components

$${\alpha }_{\omega {\omega }^{\prime }}=(u_{{\omega }^{\prime }}^{\text{in}},f_{\omega }{)}_{\text{KG}},{\beta }_{\omega {\omega }^{\prime }}=-({\bar{u}}_{{\omega }^{\prime }}^{\text{in}},f_{\omega }{)}_{\text{KG}},$$

with $u_{{\omega }^{\prime }}^{\text{in}}(v)=e^{-i{\omega }^{\prime }v}/\sqrt{4\pi {\omega }^{\prime }}$. The integrals evaluate to

$${\alpha }_{\omega {\omega }^{\prime }}\propto {\int }_{-\infty }^{v_H}(v_H-v{)}^{i\omega /\kappa }{e}^{i{\omega }^{\prime }v}dv,{\beta }_{\omega {\omega }^{\prime }}\propto {\int }_{-\infty }^{v_H}(v_H-v{)}^{i\omega /\kappa }{e}^{-i{\omega }^{\prime }v}dv.$$

Substituting $w=v_H-v$:

$${\alpha }_{\omega {\omega }^{\prime }}\propto e^{i{\omega }^{\prime }{v}_H}{\int }_0^{\infty }{w}^{i\omega /\kappa }{e}^{-i{\omega }^{\prime }w}dw,{\beta }_{\omega {\omega }^{\prime }}\propto e^{-i{\omega }^{\prime }{v}_H}{\int }_0^{\infty }{w}^{i\omega /\kappa }{e}^{i{\omega }^{\prime }w}dw.$$

The two integrals are related by analytic continuation ${\omega }^{\prime }\to -{\omega }^{\prime }$, which corresponds to rotating $w$ through the upper half-plane by $\pi$. Under this rotation $w^{i\omega /\kappa }\to e^{-\pi \omega /\kappa }{w}^{i\omega /\kappa }$ (taking the principal branch with branch cut on the negative real axis), so

$$|{\beta }_{\omega {\omega }^{\prime }}|=e^{-\pi \omega /\kappa }|{\alpha }_{\omega {\omega }^{\prime }}|,\frac{|{\beta }_{\omega {\omega }^{\prime }}{|}^2}{|{\alpha }_{\omega {\omega }^{\prime }}{|}^2}=e^{-2\pi \omega /\kappa }.$$

Combined with the Bogoliubov normalisation ${\sum }_{{\omega }^{\prime }}(|{\alpha }_{\omega {\omega }^{\prime }}{|}^2-|{\beta }_{\omega {\omega }^{\prime }}{|}^2)=1$ (taking the diagonal mode-counting interpretation), the squared $\beta$-sum is

$$⟨0_{\text{in}}|{\hat{N}}_{\omega }^{\text{out}}|0_{\text{in}}⟩=\sum _{{\omega }^{\prime }}|{\beta }_{\omega {\omega }^{\prime }}{|}^2=\frac{1}{e^{2\pi \omega /\kappa }-1},$$

a Planck distribution at temperature $T_H=\kappa /(2\pi )=1/(8\pi M)$. Including the partial-wave label and the greybody factor from backscattering gives ${\Gamma }_{\omega }/(e^{2\pi \omega /\kappa }-1)$, completing the proof. $\square$

Proposition (Euclidean conical-deficit derivation of $T_H$). The Wick-rotated Schwarzschild metric, $d{s}_E^2=(1-2M/r)d{\tau }^2+(1-2M/r{)}^{-1}d{r}^2+r^{2}d{\Omega }^2$ with $\tau =it$, is regular at $r=2M$ if and only if $\tau$ has period $\beta =8\pi GM/\hslash$ in physical units. The Hawking temperature is then $T_H=1/\beta =\hslash c^3/(8\pi GM{k}_B)$.

Proof. Introduce the proper-radial coordinate $\rho =2\sqrt{2M(r-2M)}$ near the horizon. Then $r-2M={\rho }^2/(8M)$ and $1-2M/r=({\rho }^2/(8M))/(2M+{\rho }^2/(8M))\approx {\rho }^2/(16M^2)$ to leading order in $\rho$. The Euclidean metric near $\rho =0$ becomes

$$d{s}_E^2\approx \frac{{\rho }^2}{16M^2}d{\tau }^2+d{\rho }^2+(2M{)}^{2}d{\Omega }^2.$$

The $(\rho ,\tau )$ block is the flat two-plane in polar coordinates $(\rho ,\theta )$ with $\theta =\tau /(4M)$ provided $\theta$ has period $2\pi$, i.e. $\tau$ has period $\beta =8\pi M$ in geometric units, or $\beta =8\pi GM/\hslash$ in conventional units. Any other period for $\tau$ produces a conical-deficit singularity at $\rho =0$, and the geometry fails to be smooth. By the standard thermofield-double argument, periodicity of imaginary time at period $\beta$ corresponds to a thermal state at temperature $T_H=1/\beta$. $\square$

Connections Master

• Black hole thermodynamics and area theorem 13.06.03. The direct predecessor unit develops the four laws of black-hole mechanics, the area theorem, and the Bekenstein-Hawking entropy formula in their classical-and-semiclassical forms, with the surface-gravity-to-temperature dictionary $T_H=\kappa /(2\pi )$ stated and the proportionality coefficient $S=A/4$ taken from Hawking 1974-75. The present unit supplies the missing derivation: the Bogoliubov calculation that fixes the temperature, the Euclidean-smoothness argument that re-derives it geometrically, and the Parikh-Wilczek tunnelling argument that connects the emission probability to the entropy change of the hole. With the present unit in place, the dictionary of 13.06.03 is no longer an analogy but a complete identity.

• Black holes 13.06.01. The Schwarzschild, Kerr, and Reissner-Nordström geometries developed in the predecessor unit provide the gravitational backgrounds on which the Klein-Gordon mode equation is solved and the Bogoliubov coefficients are computed. The present unit's Schwarzschild Hawking temperature $T_H=1/(8\pi M)$, Kerr generalisation $T_H={\kappa }_{Kerr}/(2\pi )$, and Reissner-Nordström analogues all follow from the surface-gravity formulas computed there.

• Bosonic Fock space and second quantisation 12.13.01. The in/out-state quantisation framework on which the Bogoliubov transformation acts is the Fock space construction of 12.13.01, with the in-vacuum $|0_{\text{in}}⟩$ defined as the cyclic vector for the in-mode creation operators and the Bogoliubov transformation acting as an automorphism of the canonical commutation relations. The present unit applies the Fock-space machinery in the curved-spacetime regime where the choice of positive-frequency decomposition is not canonical, leading to the in-vacuum / out-vacuum distinction.

• Bunch-Davies state on de Sitter 13.09.08. The de Sitter cosmological horizon carries a Gibbons-Hawking temperature $T_{dS}=H/(2\pi )$ derivable by exactly the same Euclidean-smoothness argument as the Schwarzschild case treated here, with the de Sitter Hubble rate $H$ playing the role of the surface gravity. The Bunch-Davies state on the Lorentzian section is the de Sitter analogue of the Hartle-Hawking state on the Kruskal extension of Schwarzschild. Both are the unique Hadamard states respecting the maximal symmetry of the spacetime, and both arise as thermal states of static-patch observers at the horizon temperature.

• Blackbody radiation: Planck, Stefan-Boltzmann, Wien 11.05.03. The asymptotic out-spectrum at ${\mathcal{I}}^{+}$ of the Hawking radiation, integrated over all partial waves and species, is a Planck spectrum at temperature $T_H$ (modulated by greybody factors of order unity). The total luminosity follows the Stefan-Boltzmann law $L=\sigma A{T}_H^4$, with the integration yielding the lifetime $\tau =(5120\pi G^2/\hslash c^4)M^3$ cubic in the initial mass. The present unit makes the blackbody pattern of 11.05.03 applicable to a gravitational system through the Bogoliubov mechanism that supplies the relevant thermal occupation.

Historical & philosophical context Master

Hawking's 1974 calculation is one of the great surprises of modern physics. Hawking had been attempting to refute Bekenstein's 1972-73 proposal ^{[Bekenstein1973]} that black holes carry an entropy proportional to their horizon area, on the grounds that an entropy without a temperature was thermodynamically inconsistent. The relevant temperature, $T_H=\kappa /(2\pi )$, appeared as a formal proportionality constant in the Bardeen-Carter-Hawking 1973 four-laws paper but was treated as an analogy rather than a physical temperature. Hawking computed the particle creation rate from a collapsing star using techniques he had developed with Hartle and others for cosmological particle creation in the late 1960s, and found, to his initial surprise, an exactly thermal flux at exactly the predicted temperature ^{[Hawking1974]}. The full Bogoliubov-transformation derivation appeared in Comm. Math. Phys. 43 (1975) 199 ^{[Hawking1975]}, and remains the canonical calculation; the Birrell-Davies 1982 codification follows it closely.

The Euclidean conical-deficit derivation due to Gibbons and Hawking 1977 ^{[GibbonsHawking1977]} provided a complementary, geometric perspective: the Hawking temperature is determined by the requirement that the Euclidean section of the Schwarzschild metric be smooth at the horizon, with no reference to a collapse history. The Hartle-Hawking 1976 path-integral derivation ^{[HartleHawking1976]} gave a third, the eternal-black-hole-thermofield-double construction that identifies the Hartle-Hawking state on the Kruskal extension with the maximally entangled purification of the canonical thermal density matrix on the exterior region. The Parikh-Wilczek 2000 tunnelling derivation ^{[ParikhWilczek2000]} gave a fourth, identifying the emission probability of a Hawking quantum with the WKB action equal to the change in the Bekenstein-Hawking entropy. The four derivations agree on the leading-order temperature and disagree only at the level of back-reaction corrections, which is where the unitarity question lives.

The information paradox surfaced almost immediately. Hawking's 1976 Phys. Rev. D paper argued that the exactly thermal nature of the asymptotic radiation, combined with the eventual complete evaporation of the hole, would convert an initially pure quantum state into a final mixed thermal state — violating unitarity. Don Page in 1993 ^[Page1993] turned this into a quantitative consistency condition: if the evaporation is unitary, the von Neumann entropy of the radiation must follow what is now called the Page curve, peaking at the Page time $t_{\text{Page}}\sim M^3$ when about half the original Bekenstein-Hawking entropy has been emitted. The Almheiri-Marolf-Polchinski-Sully 2013 firewall paper ^[AMPS2013] sharpened the tension: monogamy of entanglement forbids the simultaneous assumption of smooth horizons, unitary evaporation, and effective field theory outside the horizon. The 2019-2020 island-formula derivations of Penington ^{[Penington2020]} and Almheiri-Hartman-Maldacena-Shaghoulian-Tajdini ^[AHMST2020] provided what is currently the best semiclassical evidence that the Page curve is reproduced by the gravitational path integral via quantum-extremal-surface minimisation with non-empty interior islands after the Page time.

The conceptual significance of Hawking radiation extends beyond gravity. The Unruh effect (Unruh 1976 Phys. Rev. D 14, 870 ^[Unruh1976]) shows that an observer uniformly accelerated in Minkowski space measures a thermal bath at temperature $T_U=\hslash a/(2\pi c{k}_B)$, with the surface gravity replaced by the proper acceleration; the de Sitter Gibbons-Hawking temperature $T_{dS}=H/(2\pi )$ has the Hubble rate playing the analogous role. The three effects — Hawking, Unruh, Gibbons-Hawking — establish the universal thermal-state property of bifurcate Killing horizons, made rigorous by Sewell 1982 and Kay-Wald 1991 via the Bisognano-Wichmann theorem in algebraic quantum field theory. The structural unification has motivated thermodynamic-gravitational programmes including Jacobson 1995 Phys. Rev. Lett. 75, 1260 in which the Einstein equation is derived as a thermodynamic equation of state, and Verlinde 2011 JHEP 04, 029 in which gravity is interpreted as an entropic force.

The observational status of Hawking radiation remains indirect. No primordial black hole has been observed evaporating, and the gamma-ray and CMB-background constraints from Carr-Kohri-Sendouda-Yokoyama 2010 ^{[CarrKohriSendoudaYokoyama2010]} and subsequent re-analyses bound the cosmological abundance of primordial black holes in the $10^{14}$ to $10^{17}$ g evaporation window to a tiny fraction of the dark matter. Laboratory analogues — Unruh-DeWitt detectors, sonic-horizon black-hole analogues in flowing Bose-Einstein condensates due to Steinhauer 2016 Nature Phys. 12, 959 — have measured thermal Hawking-like spectra in the analogue-gravity setting, consistent with the kinematic universality of the effect. A direct astrophysical observation of Hawking radiation from a stellar-mass black hole is excluded by the eight orders of magnitude below the CMB temperature; the field remains primarily theoretical.

Bibliography Master

@article{Hawking1974,
  author = {Hawking, Stephen W.},
  title = {Black hole explosions?},
  journal = {Nature},
  volume = {248},
  year = {1974},
  pages = {30--31},
}

@article{Hawking1975,
  author = {Hawking, Stephen W.},
  title = {Particle creation by black holes},
  journal = {Comm. Math. Phys.},
  volume = {43},
  year = {1975},
  pages = {199--220},
}

@article{Davies1975,
  author = {Davies, Paul C. W.},
  title = {Scalar particle production in {S}chwarzschild and {R}indler metrics},
  journal = {J. Phys. A},
  volume = {8},
  year = {1975},
  pages = {609--616},
}

@article{Unruh1976,
  author = {Unruh, William G.},
  title = {Notes on black hole evaporation},
  journal = {Phys. Rev. D},
  volume = {14},
  year = {1976},
  pages = {870--892},
}

@article{GibbonsHawking1977,
  author = {Gibbons, Gary W. and Hawking, Stephen W.},
  title = {Action integrals and partition functions in quantum gravity},
  journal = {Phys. Rev. D},
  volume = {15},
  year = {1977},
  pages = {2752--2756},
}

@article{HartleHawking1976,
  author = {Hartle, James B. and Hawking, Stephen W.},
  title = {Path-integral derivation of black-hole radiance},
  journal = {Phys. Rev. D},
  volume = {13},
  year = {1976},
  pages = {2188--2203},
}

@article{Page1976,
  author = {Page, Don N.},
  title = {Particle emission rates from a black hole},
  journal = {Phys. Rev. D},
  volume = {13},
  year = {1976},
  pages = {198--206},
}

@article{ParikhWilczek2000,
  author = {Parikh, Maulik K. and Wilczek, Frank},
  title = {Hawking radiation as tunneling},
  journal = {Phys. Rev. Lett.},
  volume = {85},
  year = {2000},
  pages = {5042--5045},
}

@article{Page1993,
  author = {Page, Don N.},
  title = {Information in black hole radiation},
  journal = {Phys. Rev. Lett.},
  volume = {71},
  year = {1993},
  pages = {3743--3746},
}

@article{AMPS2013,
  author = {Almheiri, Ahmed and Marolf, Donald and Polchinski, Joseph and Sully, James},
  title = {Black holes: complementarity or firewalls?},
  journal = {JHEP},
  volume = {1302},
  year = {2013},
  pages = {062},
}

@article{Penington2020,
  author = {Penington, Geoffrey},
  title = {Entanglement wedge reconstruction and the information paradox},
  journal = {JHEP},
  volume = {09},
  year = {2020},
  pages = {002},
}

@article{AHMST2020,
  author = {Almheiri, Ahmed and Hartman, Thomas and Maldacena, Juan and Shaghoulian, Edgar and Tajdini, Amirhossein},
  title = {The entropy of {H}awking radiation},
  journal = {JHEP},
  volume = {1907},
  year = {2019},
  pages = {063},
}

@article{Bekenstein1973,
  author = {Bekenstein, Jacob D.},
  title = {Black holes and entropy},
  journal = {Phys. Rev. D},
  volume = {7},
  year = {1973},
  pages = {2333--2346},
}

@article{CarrKohriSendoudaYokoyama2010,
  author = {Carr, B. J. and Kohri, K. and Sendouda, Y. and Yokoyama, J.},
  title = {New cosmological constraints on primordial black holes},
  journal = {Phys. Rev. D},
  volume = {81},
  year = {2010},
  pages = {104019},
}

@book{Wald1994,
  author = {Wald, Robert M.},
  title = {Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics},
  publisher = {University of Chicago Press},
  year = {1994},
}

@book{BirrellDavies1982,
  author = {Birrell, N. D. and Davies, P. C. W.},
  title = {Quantum Fields in Curved Space},
  publisher = {Cambridge University Press},
  year = {1982},
}

@book{ParkerToms2009,
  author = {Parker, Leonard and Toms, David J.},
  title = {Quantum Field Theory in Curved Spacetime},
  publisher = {Cambridge University Press},
  year = {2009},
}