The book / chapter 22 · optional deeper trail
CHAPTER 22 · OPTIONAL DEEPER TRAIL

Singularities and black-hole thermodynamics

Follow a cloud of geodesics toward the point where geometry meets thermodynamics.

3 worked examples in this chapter
Before you begin
THE QUESTION

How can a statement about converging paths lead to black-hole entropy?

BRING WITH YOU

By the end: Read the focusing argument and distinguish classical horizon laws from semiclassical radiation.

22.1 From one falling observer to a cloud of them#

A geodesic describes one freely falling observer. Geodesic deviation describes the changing separation of nearby observers. A congruence is a smooth family of worldlines filling a region without crossing there. Think of an enormous cloud of tiny spacecraft, each with its engines off, carrying rulers to monitor its neighbors.

The cloud can expand, distort, and rotate. Those are different effects. A sphere becoming a larger sphere expands. A sphere becoming a same-volume ellipsoid shears. An ellipsoid turning without changing shape rotates.

Use c=1c=1 through the focusing discussion and let the congruence have unit tangent uμu^\mu, with uμuμ=1u^\mu u_\mu=-1. Its local rest-space metric is

hμν=gμν+uμuν.h_{\mu\nu}=g_{\mu\nu}+u_\mu u_\nu.

Define the projected velocity-gradient tensor

Bμν=hμαhνββuα.B_{\mu\nu}=h_\mu{}^\alpha h_\nu{}^\beta\nabla_\beta u_\alpha.

It tells neighboring observers how their relative velocity depends, to first order, on their spatial separation. Split this three-dimensional linear map into its trace, symmetric traceless part, and antisymmetric part:

Bμν=13θhμν+σμν+ωμν.\boxed{B_{\mu\nu}=\frac13\theta h_{\mu\nu} +\sigma_{\mu\nu}+\omega_{\mu\nu}.}

Explicitly,

θ=μuμ,σμν=B(μν)13θhμν,ωμν=B[μν].\theta=\nabla_\mu u^\mu, \qquad \sigma_{\mu\nu}=B_{(\mu\nu)}-\frac13\theta h_{\mu\nu}, \qquad \omega_{\mu\nu}=B_{[\mu\nu]}.

The equality between the projected trace and μuμ\nabla_\mu u^\mu follows from differentiating the fixed normalization of uu.

θ\theta is expansion, σ\sigma is shear, and ω\omega is vorticity or twist. For an infinitesimal comoving volume V\mathcal V,

θ=1VdVdτ=ddτlnV.\theta=\frac1{\mathcal V}\frac{d\mathcal V}{d\tau} =\frac{d}{d\tau}\ln\mathcal V.

The trace is the fractional volume-growth rate because the first-order fractional change of a determinant is the trace of the underlying linear deformation. That is the same determinant identity that appeared when varying g\sqrt{-g}.

There is also a direct bridge to geodesic deviation. For a commuting family of geodesics with separation vector ξ\xi, [u,ξ]=0[u,\xi]=0 implies

uξ=ξu.\nabla_u\xi=\nabla_\xi u.

In a transported rest frame, this is “separation velocity equals BB times separation.” Differentiating again and using the geodesic-deviation equation gives a matrix evolution law with a term B2-B^2 and a curvature term. Raychaudhuri’s equation is its trace.

22.2 Deriving the Raychaudhuri equation#

Start with the expansion and differentiate along the flow:

dθdτ=uρρ(μuμ).\frac{d\theta}{d\tau} =u^\rho\nabla_\rho(\nabla_\mu u^\mu).

Commute the derivatives, using the Riemann convention of this book. The contraction contributes a minus Ricci term. Then apply the product rule:

dθdτ=μ(uρρuμ)(μuρ)(ρuμ)Rμνuμuν.\frac{d\theta}{d\tau} =\nabla_\mu(u^\rho\nabla_\rho u^\mu) -(\nabla_\mu u^\rho)(\nabla_\rho u^\mu) -R_{\mu\nu}u^\mu u^\nu.

The first term is μaμ\nabla_\mu a^\mu, where aμ=uuμa^\mu=\nabla_u u^\mu. It vanishes for a geodesic congruence. The second is the trace of the square of the velocity-gradient map. For geodesic flow its rest-space decomposition gives

tr(B2)=13θ2+σμνσμνωμνωμν.\operatorname{tr}(B^2) =\frac13\theta^2+\sigma_{\mu\nu}\sigma^{\mu\nu} -\omega_{\mu\nu}\omega^{\mu\nu}.

Why the vorticity minus sign? A real antisymmetric matrix squares to a matrix with nonpositive trace. In two dimensions, the rotation generator

(0ww0)\begin{pmatrix}0&-w\\w&0\end{pmatrix}

squares to w2I-w^2I. The trace-free symmetric matrix has a positive sum of squared eigenvalues. Cross terms vanish because a symmetric tensor contracted with an antisymmetric tensor is zero, and shear has zero trace.

Therefore, for a timelike geodesic congruence in four spacetime dimensions,

dθdτ=13θ2σμνσμν+ωμνωμνRμνuμuν.\boxed{ \frac{d\theta}{d\tau} =-\frac13\theta^2 -\sigma_{\mu\nu}\sigma^{\mu\nu} +\omega_{\mu\nu}\omega^{\mu\nu} -R_{\mu\nu}u^\mu u^\nu. }

For a normalized accelerated congruence, add +μaμ+\nabla_\mu a^\mu to this expression, with the same projected definitions. Rockets can alter the cloud’s behavior; geodesic focusing theorems do not silently include thrust.

The terms describe distinct contributions:

  • θ2/3-\theta^2/3: convergence can reinforce itself even without curvature.
  • σ2-\sigma^2: stretching along some directions can accelerate volume focusing.
  • +ω2+\omega^2: rotation opposes the simple focusing tendency.
  • Rμνuμuν-R_{\mu\nu}u^\mu u^\nu: curvature directly changes the trace of relative acceleration.

The Weyl tensor is absent from the explicit last term, but it can generate shear, which then affects expansion through σ2-\sigma^2. Vacuum curvature can matter enormously even when Rμν=0R_{\mu\nu}=0.

22.3 How an inequality becomes a finite-time prediction#

Suppose the geodesics are orthogonal to a family of spacelike slices. Locally their covector is then uμ=Nμtu_\mu=-N\partial_\mu t, where tt labels the slices and NN normalizes uu. Antisymmetrizing its derivative cancels the second derivatives of tt. The terms left contain a factor μt\partial_\mu t; projecting both indices into a slice removes them. Thus ωμν=0\omega_{\mu\nu}=0 for this flow.

The converse local statement, that zero twist allows such orthogonal slices, is the codimension-one Frobenius integrability theorem. The focusing calculation below needs only the forward implication just checked.

Also suppose the timelike convergence condition holds:

Rμνuμuν0.R_{\mu\nu}u^\mu u^\nu\ge0.

Raychaudhuri then implies

dθdτ13θ2.\frac{d\theta}{d\tau}\le-\frac13\theta^2.

If the cloud starts converging, θ0<0\theta_0<0, divide by θ2>0\theta^2>0 and differentiate its reciprocal:

ddτ(1θ)13.\frac{d}{d\tau}\left(\frac1\theta\right)\ge\frac13.

If the geodesics extend through the required interval, this inequality forces the initially negative reciprocal toward zero within proper time no larger than 3/θ03/|\theta_0|: θ\theta becomes unboundedly negative, and the smooth congruence develops a focal point or caustic by that bound. An earlier incomplete geodesic endpoint is another possibility. The focusing argument does not by itself guarantee that the geodesics exist for the entire interval.

This is not yet a spacetime singularity. Aim a family of straight worldlines at the same event in Minkowski spacetime. Their congruence focuses; each worldline continues perfectly well. What breaks is the single-valued smooth velocity field used to describe that overlapping family.

The volume also makes the focusing bound explicit. Put v=(V/V0)1/3v=(\mathcal V/\mathcal V_0)^{1/3}, with v(0)=1v(0)=1. While the congruence is regular,

vv=θ˙3+θ290.\frac{v''}{v}=\frac{\dot\theta}{3}+\frac{\theta^2}{9}\le0.

Its initial slope is v(0)=θ0/3<0v'(0)=\theta_0/3<0. Since that slope cannot increase, v(τ)1+θ0τ/3v(\tau)\le1+\theta_0\tau/3. A positive smooth volume therefore cannot persist beyond 3/θ03/|\theta_0|. This establishes focusing if the geodesics extend that far; inferring spacetime incompleteness still requires the additional global argument in a singularity theorem.

22.4 Energy conditions and focusing#

Einstein’s equation converts curvature conditions into matter conditions, but it does not itself require ordinary matter to satisfy those conditions. They are additional hypotheses.

For a perfect fluid, useful pointwise conditions are:

Condition General idea Perfect-fluid inequalities
Null energy condition, NEC Tμνkμkν0T_{\mu\nu}k^\mu k^\nu\ge0 for every null kk ϵ+p0\epsilon+p\ge0
Weak energy condition, WEC Every timelike observer measures nonnegative local energy density ϵ0\epsilon\ge0, ϵ+p0\epsilon+p\ge0
Dominant energy condition, DEC Energy density is nonnegative and its flux is causal ϵp\epsilon\ge\lvert p\rvert
Strong energy condition, SEC (Tμν12Tgμν)vμvν0(T_{\mu\nu}-\tfrac12Tg_{\mu\nu})v^\mu v^\nu\ge0 for all timelike vv ϵ+p0\epsilon+p\ge0, ϵ+3p0\epsilon+3p\ge0

To verify the fluid inequalities, use its orthonormal rest frame and a unit timelike observer vμ^=γ(1,β)v^{\hat\mu}=\gamma(1,\boldsymbol\beta), with β<1|\boldsymbol\beta|<1. The measured density is

Tμνvμvν=(ϵ+p)γ2p.T_{\mu\nu}v^\mu v^\nu=(\epsilon+p)\gamma^2-p.

At rest this is ϵ\epsilon; as γ\gamma becomes large its sign is controlled by ϵ+p\epsilon+p. Nonnegativity for every observer therefore requires both WEC inequalities. A null vector has the form kμ^=q(1,n)k^{\hat\mu}=q(1,\mathbf n) with n=1|\mathbf n|=1, giving Tμνkμkν=q2(ϵ+p)T_{\mu\nu}k^\mu k^\nu=q^2(\epsilon+p) and the NEC.

For the SEC, add T/2=(ϵ+3p)/2T/2=(-\epsilon+3p)/2 to the timelike contraction. Its value at rest is (ϵ+3p)/2(\epsilon+3p)/2 and its large-γ\gamma coefficient is again ϵ+p\epsilon+p. For the DEC, the energy-flux vector seen by vv is Jμ^=Tμ^ν^vν^=(γϵ,γpβ)J^{\hat\mu}=-T^{\hat\mu}{}_{\hat\nu}v^{\hat\nu}=(\gamma\epsilon,-\gamma p\boldsymbol\beta). Requiring it to be future causal for every β<1|\boldsymbol\beta|<1 gives ϵ0\epsilon\ge0 and ϵ2p2\epsilon^2\ge p^2, or ϵp\epsilon\ge|p|.

For Λ=0\Lambda=0, the SEC implies timelike convergence. If Λ\Lambda remains on the geometric side, however,

Rμνuμuν=8πGN(Tμνuμuν+12T)Λ.R_{\mu\nu}u^\mu u^\nu =8\pi G_N\left(T_{\mu\nu}u^\mu u^\nu+\frac12T\right)-\Lambda.

For a comoving perfect-fluid observer this is 4πGN(ϵ+3p)Λ4\pi G_N(\epsilon+3p)-\Lambda. Positive Λ\Lambda can defeat timelike focusing. Equivalently, move it into an effective vacuum stress tensor with pΛ=ϵΛp_\Lambda=-\epsilon_\Lambda: it violates the SEC when its density is positive, while saturating the NEC.

For null vectors, the trace and cosmological terms vanish because gμνkμkν=0g_{\mu\nu}k^\mu k^\nu=0. The NEC therefore implies null convergence in GR even with Λ\Lambda.

Classical scalar potentials can violate the SEC, and quantum fields can violate classical pointwise energy conditions more broadly. These are reasons to inspect a theorem’s hypotheses carefully, not to call the theorem mistaken.

22.5 Trapped surfaces and what Penrose actually proved#

A null ray has no unit rest frame: its tangent kk is perpendicular to itself. To isolate the two transverse directions, choose another null vector ll with kl=1k\cdot l=-1 and define

qμν=gμν+kμlν+lμkν.q_{\mu\nu}=g_{\mu\nu}+k_\mu l_\nu+l_\mu k_\nu.

This tensor removes both kk and ll components and supplies a positive metric on the remaining two-dimensional screen. Project the separation-velocity map onto that screen and split it into trace, shear, and twist as before. The trace is now an area-growth rate, so the trace part of the two-dimensional map is θqμν/2\theta q_{\mu\nu}/2.

For affinely parametrized geodesic rays, taking the trace of the same evolution calculation gives

dθdλ=12θ2σμνσμν+ωμνωμνRμνkμkν.\frac{d\theta}{d\lambda} =-\frac12\theta^2-\sigma_{\mu\nu}\sigma^{\mu\nu} +\omega_{\mu\nu}\omega^{\mu\nu} -R_{\mu\nu}k^\mu k^\nu.

The shear and twist here live on the positive-definite two-dimensional screen transverse to the rays. An affine parameter is essential; a nonaffine parameter introduces an additional term proportional to θ\theta.

Take a compact spacelike two-surface without a boundary, such as a sphere, and send future light rays orthogonally away from it in both null-normal directions. For an ordinary sphere in flat space, the outward bundle grows in area and the inward bundle shrinks. A future trapped surface has negative expansion in both directions. Even the outward-directed light bundle initially loses cross-sectional area.

This is a local geometric condition on the surface and its null normals. It does not say that a photon locally travels more slowly than light, nor does it require a coordinate speed to become negative.

One standard form of Penrose’s theorem says that a sufficiently regular spacetime with a noncompact Cauchy hypersurface, null convergence, and a closed future trapped surface must be future null geodesically incomplete. The theorem combines focusing with global causal geometry; it does not assume spherical symmetry. The original result is Penrose’s 1965 paper, “Gravitational Collapse and Space-Time Singularities”.

Geodesic incompleteness means at least one inextendible geodesic has finite affine length in the relevant direction; for timelike geodesics, proper time is the physical parameter. It does not universally mean that a curvature scalar tends to infinity. The theorem does not supply a location, a topology, or a detailed microscopic description of “the singularity.”

Even incompleteness must be interpreted carefully: deleting one point from otherwise regular Minkowski spacetime creates incomplete geodesics artificially. One must consider extendibility and which spacetime has actually been specified. Conversely, a coordinate singularity that disappears in a larger smooth chart is not evidence that physics has ended.

The theorem applies without assuming a spherical source. Its conclusion is incompleteness under the stated global and convergence conditions, even when the matter distribution is less symmetric than the examples solved earlier.

Focusing is not yet a singularityFive straight timelike paths converge to a caustic, alongside the focusing time bound. This flat-spacetime schematic intentionally has no curvature. Raychaudhuri can force a congruence to focus, but singularity theorems require additional causal and global assumptions to conclude geodesic incompleteness.36 / FOCUSING IS NOT YET A SINGULARITYpositionaffine timecausticConvergence can have a deadlineBut straight lines can cross in flat spacetime.Global hypotheses do the extra work.
36 /
Focusing is not yet a singularity. This flat-spacetime schematic intentionally has no curvature. Raychaudhuri can force a congruence to focus, but singularity theorems require additional causal and global assumptions to conclude geodesic incompleteness.

22.6 Horizons are about causal access; predictability needs another definition#

In an asymptotically flat spacetime, the black-hole region consists of events unable to send a causal signal to future null infinity I+\mathscr I^+, the ideal destination of escaping light. Its boundary is the event horizon:

H+=J(I+).\mathcal H^+=\partial J^-(\mathscr I^+).

J(I+)J^-(\mathscr I^+) denotes the causal past of that destination. This definition depends on the entire future spacetime. An event horizon is not generally locatable using measurements made at one point and one instant.

A marginally outer trapped surface instead has vanishing outward null expansion, usually with negative inward expansion in the black-hole setting. Under suitable conditions an apparent horizon is the outer boundary of the trapped region on a chosen slice. These are valuable tools for dynamical calculations, but they depend on slicing and need not coincide with the event horizon.

For predictability, define the future domain of dependence D+(Σ)D^+(\Sigma): an event belongs to it if every past-inextendible causal curve through that event intersects Σ\Sigma. No causal influence can arrive there without passing through the supplied initial data.

A Cauchy surface meets every inextendible causal curve exactly once. A spacetime admitting such a surface is globally hyperbolic, in the usual boundary-free setting. This is the natural arena for the initial-value description of Chapter 20.

A Cauchy horizon bounds a domain of dependence. Beyond it, initial data on the chosen surface no longer control every possible incoming influence. It is different from an event horizon: crossing a black-hole event horizon need not destroy local predictability. Inner horizons in idealized charged or rotating solutions motivate the question of whether such extendible boundaries survive generic perturbations. Strong cosmic censorship studies this issue, and its precise claims depend on the matter model and the allowed regularity of extensions.

A final geometric tool helps organize these questions. Under a smooth positive conformal rescaling gμνΩ2gμνg_{\mu\nu}\mapsto\Omega^2g_{\mu\nu}, null cones are unchanged. Proper times and distances change. Conformal diagrams exploit this separation to compress enormous regions while preserving causal relations; they are maps of who can signal whom, not faithful scale drawings.

22.7 A horizon gets a temperature#

PREPARATION FOR THIS SECTION

What temperature adds to a theory

What does a thermal probability distribution know that a spacetime metric alone does not?

See the idea

A macrostate specifies coarse measurements such as energy and volume. A microstate specifies the finer physical configuration. A probability distribution assigns nonnegative numbers pip_i to the allowed microstates with ipi=1\sum_i p_i=1. An average measurement of a quantity AiA_i is A=ipiAi\langle A\rangle=\sum_i p_iA_i. The distribution is extra physical information.

Work it out
  1. Measure uncertainty with entropy

    For probabilities pip_i, define S=kBipilnpiS=-k_B\sum_i p_i\ln p_i, with 0ln00\ln0 interpreted by continuity as zero. kBk_B converts dimensionless logarithmic uncertainty into entropy units. For N equally likely states, substitution gives S=kBlnNS=k_B\ln N.

    pi=1/NS=kBlnN.p_i=1/N\quad\Longrightarrow\quad S=k_B\ln N.

    Why this step works Entropy depends on the probabilities and the chosen description of states.

  2. Find the canonical distribution

    Maximize S/kBS/k_B with fixed normalization and mean energy. Multipliers give lnpi1αβEi=0-\ln p_i-1-\alpha-\beta E_i=0, hence pi=eβEi/Zp_i=e^{-\beta E_i}/Z, where Z=jeβEjZ=\sum_j e^{-\beta E_j} normalizes the probabilities. For an equilibrium heat reservoir, β=1/(kBT)\beta=1/(k_BT).

    pi=eEi/(kBT)Z,pjpi=e(EjEi)/(kBT).p_i=\frac{e^{-E_i/(k_BT)}}Z,\qquad \frac{p_j}{p_i}=e^{-(E_j-E_i)/(k_BT)}.

    Why this step works The energy difference relative to the thermal energy kBT sets the population ratio.

  3. Relate entropy to energy transfer

    Reversible heat exchange obeys dS=δQrev/TdS=\delta Q_{\rm rev}/T. With fixed particle number, pressure work gives dE=TdSpdVdE=T\,dS-p\,dV. Irreversible heat transfer does not justify replacing the full entropy balance by δQ/T\delta Q/T with an arbitrary system temperature.

    (SE)V,N=1T.\left(\frac{\partial S}{\partial E}\right)_{V,N}=\frac1T.

    Why this step works Thermodynamic derivatives specify what is held fixed.

Go deeper

For two states with energies 0 and Δ>0\Delta>0, Z=1+eΔ/(kBT)Z=1+e^{-\Delta/(k_BT)}. At high positive temperature their probabilities approach one half; at low positive temperature the lower state dominates. Classical thermodynamics and probability do not by themselves imply that a black hole emits particles. The Hawking calculation additionally uses quantum fields, a choice of state, and a regularity or collapse condition. Euclidean periodicity supplies a consistency relation between a candidate thermal state and the near-horizon geometry; it is not the entire flux calculation.

Test the idea

FIRST, PREDICT

Does a metric by itself specify a quantum field’s thermal population?

Compare the reasoning

Temperature is simply another coordinate time.

Temperature controls energy populations; a coordinate labels events.

No. The field theory and state are additional information.

The same geometry can support distinct physical states with different observations.

Yes. Geometry replaces the probability distribution.

Geometrical regularity can constrain a state, but does not erase the need to specify one.

A hint

Separate the arena, dynamical law, and state.

NOW CHANGE THE EXAMPLE

A system has four equally likely states. Find S/(kBln2)S/(k_B\ln2).

A hint

Use S=kBln4S=k_B\ln4.

Work through the solution

ln4/ln2=2\ln4/\ln2=2.

Temperature is a property of a specified physical state, and entropy depends on the information retained.

PREPARATION FOR THIS SECTION

A pure whole can have an uncertain part

How can losing access to one subsystem create entropy without destroying the whole quantum state?

See the idea

A finite quantum system is described by a normalized vector of complex amplitudes. In an orthonormal basis i|i\rangle, write ψ=iaii|\psi\rangle=\sum_i a_i|i\rangle with iai2=1\sum_i|a_i|^2=1. The probability of basis outcome i is ai2|a_i|^2. This is a postulate of quantum mechanics, separate from GR. The notation ψ\langle\psi| means the conjugate transpose of ψ|\psi\rangle.

Work it out
  1. Turn a vector into a density matrix

    For a known pure state, ϱ=ψψ\varrho=|\psi\rangle\langle\psi|. For a mixture of preparations ψj|\psi_j\rangle with probabilities pjp_j, use ϱ=jpjψjψj\varrho=\sum_jp_j|\psi_j\rangle\langle\psi_j|. This positive matrix has trace one. An observable is a Hermitian matrix A; its mean is tr(ϱA)\operatorname{tr}(\varrho A).

    S=kBtr(ϱlnϱ).S=-k_B\operatorname{tr}(\varrho\ln\varrho).

    Why this step works Diagonalizing the density matrix reduces the entropy formula to the classical probability formula for its eigenvalues.

  2. Trace over what is inaccessible

    For two two-state systems, take Ψ=(00+11)/2|\Psi\rangle=(|00\rangle+|11\rangle)/\sqrt2. The full density matrix is pure and has entropy zero. If only the first system is accessible, define ϱA=b=0,1bϱABb\varrho_A=\sum_{b=0,1}\langle b|\varrho_{AB}|b\rangle, where the inner products act on the second system. The cross terms vanish because 01=0\langle0|1\rangle=0.

    ϱA=12(1001),SA=kBln2.\varrho_A=\frac12\begin{pmatrix}1&0\\0&1\end{pmatrix},\qquad S_A=k_B\ln2.

    Why this step works A partial trace keeps exactly the information needed for measurements on the accessible subsystem.

  3. Keep evolution separate from restriction

    A closed system evolves by a unitary matrix UU, with UU=IU^\dagger U=I and ϱUϱU\varrho\mapsto U\varrho U^\dagger. This preserves its eigenvalues and therefore its entropy. Restricting to a subsystem is a different operation and can change the entropy you assign to that subsystem.

    tr[(UϱU)2]=tr(ϱ2).\operatorname{tr}[(U\varrho U^\dagger)^2]=\operatorname{tr}(\varrho^2).

    Why this step works Unitary similarity preserves spectra; a partial trace does not preserve the full state’s spectrum.

Go deeper

A free field can be decomposed into normal modes; in a chosen stationary background each mode behaves as a quantum oscillator with energies En=ω(n+1/2)E_n=\hbar\omega(n+1/2). Summing the canonical geometric series gives mean occupation n=1/(eω/(kBT)1)\langle n\rangle=1/(e^{\hbar\omega/(k_BT)}-1) for bosonic modes with zero chemical potential. The oscillator spectrum and field quantization are additional quantum inputs; this short preparation does not derive them from classical geometry. In curved spacetime different choices of positive-frequency modes can define different vacuum states. A full Hawking-flux derivation must specify those modes, the state, and transmission through the exterior curvature potential. The finite example above explains why an exterior thermal description alone is not proof that global evolution destroyed information.

Test the idea

FIRST, PREDICT

A subsystem has a mixed density matrix. Must the full system also be mixed?

Compare the reasoning

No. An entangled pure state can have a mixed reduced state.

The two-system example has zero total entropy and a nonzero reduced entropy.

Yes. Pure states have pure parts.

Entanglement is precisely why this inference fails.

A mixed reduced state proves nonunitary evolution.

Partial tracing is not the same operation as evolving the full closed system.

A hint

Calculate the partial trace of the explicitly supplied two-system state.

NOW CHANGE THE EXAMPLE

For ϱA=I2/2\varrho_A=I_2/2, calculate the purity tr(ϱA2)\operatorname{tr}(\varrho_A^2).

A hint

Square the two diagonal entries and add.

Work through the solution

1/4+1/4=1/21/4+1/4=1/2. A normalized pure state instead has purity 1.

Mixed information about a part and loss of information in the whole are distinct claims.

Restore cc and \hbar. For a Schwarzschild black hole,

rs=2GNMc2,A=4πrs2.r_s=\frac{2G_NM}{c^2}, \qquad A=4\pi r_s^2.

A short derivation reveals why a temperature appears. Near the horizon write r=rs+xr=r_s+x, with xrsx\ll r_s, so 1rs/rx/rs1-r_s/r\approx x/r_s. Continue the stationary time coordinate to imaginary time t=iτEt=-i\tau_E. The near-horizon radial-time metric becomes

dsE2xrsc2dτE2+rsxdx2.ds_E^2\approx\frac{x}{r_s}c^2d\tau_E^2+\frac{r_s}{x}dx^2.

Introduce a radial proper-distance coordinate ϱ=2rsx\varrho=2\sqrt{r_sx}. Direct substitution gives

dsE2dϱ2+ϱ2(cdτE2rs)2.ds_E^2\approx d\varrho^2 +\varrho^2\left(\frac{c\,d\tau_E}{2r_s}\right)^2.

This is a flat plane in polar coordinates. Its angular coordinate is cτE/(2rs)c\tau_E/(2r_s). Smoothness at the origin requires that angle to have period 2π2\pi, so imaginary time has period

ΔτE=4πrsc=8πGNMc3.\Delta\tau_E=\frac{4\pi r_s}{c} =\frac{8\pi G_NM}{c^3}.

Quantum statistical mechanics identifies thermal equilibrium with imaginary-time period /(kBT)\hbar/(k_BT). Equating these periods yields

TH=c38πGNMkB.\boxed{T_H=\frac{\hbar c^3}{8\pi G_NMk_B}.}

This Euclidean argument characterizes the regular stationary thermal construction. Deriving the outgoing flux in a collapse spacetime involves a quantum-field calculation with the appropriate state and boundary conditions; it yields the same Hawking temperature at infinity, with frequency-dependent transmission factors modifying an ideal blackbody spectrum. The calculation is semiclassical, not a complete quantization of spacetime. See Wald’s research review of black-hole thermodynamics.

The commonly illustrated story of a particle pair appearing exactly on the horizon is a heuristic, not this derivation. Hawking radiation depends on quantum field modes, the state, and the global relation between early and late notions of positive frequency. It is not adequately explained by assigning ordinary local particle trajectories to vacuum fluctuations.

WORKED EXAMPLE

Separate the geometry from the quantum assumption

Which part of black-hole temperature can we derive with calculus, and what new physics must we import?

See the idea

A formula involving temperature needs more than a metric. We need to say what temperature means, how quantum states evolve, and why thermal equilibrium can be represented using imaginary time. We can check the algebra of that connection here while keeping its physical assumptions visible.

Work it out
  1. Give temperature and entropy operational jobs

    For a reversible change with the other thermodynamic work variables fixed, the first law reads dE=TdSdE=T\,dS. Energy EE is measured in joules, temperature TT in kelvin, and entropy SS in joules per kelvin. Thus 1/T=dS/dE1/T=dS/dE along this family of equilibrium states. In quantum statistical mechanics, a state of energy EE has thermal weight proportional to eE/(kBT)e^{-E/(k_BT)}, where Boltzmann’s constant kBk_B converts temperature to an energy scale. This thermal-weight rule is new physical input.

    dE=TdS,w(E)eE/(kBT).dE=T\,dS,\qquad w(E)\propto e^{-E/(k_BT)}.

    Why this step works A thermal weight is not derived from Einstein’s equation. Specifying the thermodynamic ensemble and its constraints is part of using it.

  2. Compare a quantum phase with a thermal weight

    A stationary quantum state of energy EE evolves with phase eiEt/e^{-iEt/\hbar}, where \hbar is the reduced Planck constant. Formally substitute t=iτEt=-i\tau_E; this is an analytic continuation, not a clock traveling into an imaginary physical direction. Since (i)(i)=1(-i)(-i)=-1, the phase becomes eEτE/e^{-E\tau_E/\hbar}. Matching it to the thermal weight gives the imaginary-time interval ΔτE=/(kBT)\Delta\tau_E=\hbar/(k_BT).

    eiEt/  t=iτE  eEτE/,ΔτE=kBT.e^{-iEt/\hbar}\ \xrightarrow{\ t=-i\tau_E\ }\ e^{-E\tau_E/\hbar},\qquad \Delta\tau_E=\frac{\hbar}{k_BT}.

    Why this step works The algebra explains the length of the thermal interval. The full periodicity statement comes from the thermal trace and quantum correlation functions; it is not established by this substitution alone.

  3. Let regular geometry determine that interval

    The near-horizon calculation in this section puts the Euclidean radial-time metric into polar form dϱ2+ϱ2dθE2d\varrho^2+\varrho^2d\theta_E^2, with θE=cτE/(2rs)\theta_E=c\tau_E/(2r_s). A smooth plane requires angular period 2π2\pi; a different period produces a conical tip. Therefore ΔτE=4πrs/c\Delta\tau_E=4\pi r_s/c. Equating this geometric period to the thermal interval gives the Schwarzschild temperature at infinity.

    TH=c4πrskB=c38πGNMkB.T_H=\frac{\hbar c}{4\pi r_sk_B}=\frac{\hbar c^3}{8\pi G_NMk_B}.

    Why this step works Smooth Euclidean geometry fixes a period. Quantum statistical mechanics supplies its interpretation as a temperature.

  4. Check what the first law then predicts

    For the nonrotating, uncharged equilibrium family, use E=Mc2E=Mc^2. Substituting THT_H in dS=dE/THdS=dE/T_H gives a term proportional to MdMM\,dM. Integrating produces an entropy proportional to M2M^2. Because A=4πrs2=16πGN2M2/c4A=4\pi r_s^2=16\pi G_N^2M^2/c^4, it is also proportional to horizon area. The integration leaves an arbitrary additive constant S0S_0.

    dS=8πkBGNMcdM,S=kBAc34GN+S0.dS=\frac{8\pi k_BG_NM}{\hbar c}\,dM,\qquad S=\frac{k_BAc^3}{4\hbar G_N}+S_0.

    Why this step works The first law fixes entropy differences. The conventional additive normalization does not emerge from this integration.

Go deeper

The stationary Euclidean construction is an equilibrium argument. For bosonic fields, the corresponding thermal correlation functions are periodic in imaginary time; fermionic fields have the associated antiperiodic condition. A collapse spacetime’s outgoing Hawking flux requires a quantum-field calculation with a specified state and boundary conditions. It yields the same asymptotic temperature scale, while scattering by the spacetime modifies the observed spectrum. This is semiclassical physics: quantum fields on a classical background, with controlled treatment of backreaction where applicable. It does not derive a microscopic count of black-hole states, explain the final stage of evaporation, or establish a complete quantum theory of gravity. Those are open tasks, not consequences hidden in the elementary algebra.

Test the idea

FIRST, PREDICT

Which statement correctly identifies the extra physical input in the Euclidean argument?

Compare the reasoning

Einstein’s classical equation alone derives the thermal Boltzmann weights.

The Boltzmann weights and quantum phase evolution are additional physical principles, not results of classical GR.

Imaginary time is the proper time recorded by an observer inside the horizon.

Imaginary time is an analytic-continuation tool. It is not the real proper time of a traveler.

Quantum statistical mechanics relates an imaginary-time interval to temperature.

Yes. Geometry determines the interval required for regularity; quantum thermal physics provides the relation to temperature.

A hint

Locate the first appearance of both \hbar and kBk_B in the chain of reasoning.

NOW CHANGE THE EXAMPLE

Compare two nonrotating, uncharged black holes with masses MM and 3M3M. Using the conventional S0=0S_0=0, what is S(3M)/S(M)S(3M)/S(M)?

A hint

The radius is proportional to mass, and the entropy is proportional to the horizon area.

Work through the solution

rs(3M)=3rs(M)r_s(3M)=3r_s(M), so the area and entropy each increase by 32=93^2=9. The temperature instead falls to one third.

The temperature calculation joins regular geometry to explicitly stated quantum and thermodynamic input.

22.8 Entropy is written in area#

For the nonrotating, uncharged family, the first law is

d(Mc2)=THdSBH.d(Mc^2)=T_H\,dS_{\mathrm{BH}}.

Substitute THT_H and solve for the entropy change:

dSBH=8πkBGNMcdM.dS_{\mathrm{BH}}=\frac{8\pi k_BG_NM}{\hbar c}\,dM.

Integrating gives, with the conventional additive normalization,

SBH=4πkBGNM2c=kBAc34GN=kBA4P2,P=GNc3.S_{\mathrm{BH}}=\frac{4\pi k_BG_NM^2}{\hbar c} =\boxed{\frac{k_BAc^3}{4\hbar G_N} =\frac{k_BA}{4\ell_P^2}}, \qquad \ell_P=\sqrt{\frac{\hbar G_N}{c^3}}.

For ordinary extensive matter at fixed local conditions, doubling volume doubles entropy. Here the entropy scales with horizon area: doubling the Schwarzschild mass multiplies both area and entropy by four. The formula determines that scaling without specifying microscopic constituents of the horizon.

Since THM1T_H\propto M^{-1}, a Schwarzschild black hole becomes hotter as it loses mass. Its heat capacity, C=d(Mc2)/dTH=Mc2/THC=d(Mc^2)/dT_H=-Mc^2/T_H, is negative. A slightly hotter hole loses energy to a bath and becomes hotter still, rather than relaxing back by the usual positive-heat-capacity mechanism. Equilibrium therefore requires an analysis of the whole system and its boundary conditions.

The classical horizon-area theorem requires the relevant convergence and global regularity assumptions. Hawking evaporation does not contradict it: the quantum stress tensor need not satisfy the classical energy hypothesis, and the horizon area can decrease. The thermodynamic quantity then involves generalized entropy,

Sgen=kBA4P2+Soutside,S_{\mathrm{gen}}=\frac{k_BA}{4\ell_P^2}+S_{\mathrm{outside}},

with the quantum-field entropy and gravitational parameters treated consistently under renormalization. This means that when short-distance field contributions are regulated, the theory’s parameters must be adjusted consistently so the physical prediction does not depend on the arbitrary regulator. Chapter 23 develops that distinction between a regulated intermediate expression and a prediction. The generalized second law has substantial support and proofs in specified settings; it should not be promoted without qualifications to a theorem covering every unknown quantum-gravitational process.

22.9 Classical area increase, with the assumptions visible#

For horizon-generating null geodesics, the vorticity vanishes and the screen has two dimensions. In an affine parameter λ\lambda, Raychaudhuri becomes

dθdλ=12θ2σabσabRμνkμkν.\frac{d\theta}{d\lambda} =-\frac12\theta^2-\sigma_{ab}\sigma^{ab} -R_{\mu\nu}k^\mu k^\nu.

The null energy condition, together with Einstein’s equation, makes the last contraction nonnegative; the cosmological term drops out because gμνkμkν=0g_{\mu\nu}k^\mu k^\nu=0. If the horizon expansion were negative, the inequality would force a future caustic. Under the global regularity and predictability assumptions of the classical area theorem, horizon generators cannot end in that way on the future horizon. This yields nonnegative expansion and nondecreasing horizon area. The global step is essential; a local differential equation alone does not prove the theorem.

For a simple numerical illustration, imagine two initially well-separated, nonspinning holes of equal mass MM, ending in a nonspinning hole of mass MfM_f. Since a Schwarzschild area is 16πGN2M2/c416\pi G_N^2M^2/c^4, area increase requires

Mf22M2,Mf2M.M_f^2\ge2M^2,\qquad M_f\ge\sqrt2\,M.

With initial total energy approximately 2Mc22Mc^2, the radiated fraction is consequently at most 11/229.3%1-1/\sqrt2\simeq29.3\% under these idealizations. This is an upper bound, not the predicted emission efficiency. A spinning remnant requires the Kerr area formula, so applying AM2A\propto M^2 blindly to a measured merger would be wrong.

For an uncharged stationary rotating black hole, the first law can be written in SI units as

d(Mc2)=κsgc28πGNdA+ΩHdJ.d(Mc^2)=\frac{\kappa_{\rm sg}c^2}{8\pi G_N}\,dA+\Omega_H\,dJ.

Here κsg\kappa_{\rm sg} is surface gravity with acceleration units, ΩH\Omega_H is horizon angular velocity, and JJ is angular momentum. For Schwarzschild, surface gravity is the limiting hovering acceleration after correcting by its redshift factor:

κsg=limrrs+N(r)aproper(r)=GNMrs2=c22rs.\kappa_{\rm sg}=\lim_{r\to r_s^+}N(r)a_{\rm proper}(r) =\frac{G_NM}{r_s^2}=\frac{c^2}{2r_s}.

The local hovering acceleration diverges; this redshifted limit remains finite and uses the clock normalization at infinity.

For the nonextremal Kerr family, aK<m|a_K|<m, the corresponding quantities are

A=4π(r+2+aK2),ΩH=caKr+2+aK2,κsg=c2(r+r)2(r+2+aK2).\begin{aligned} A&=4\pi(r_+^2+a_K^2),\\ \Omega_H&=\frac{ca_K}{r_+^2+a_K^2},\\ \kappa_{\rm sg}&=\frac{c^2(r_+-r_-)}{2(r_+^2+a_K^2)}. \end{aligned}

These formulas use the horizon radii defined in Section 17.7. The area follows by integrating the horizon cross-section: gθθgϕϕ=(r+2+aK2)sinθ\sqrt{g_{\theta\theta}g_{\phi\phi}}=(r_+^2+a_K^2)\sin\theta. The angular velocity is the horizon limit of gtϕ/gϕϕ-g_{t\phi}/g_{\phi\phi}. To check the first law within this family, differentiate r+22mr++aK2=0r_+^2-2mr_++a_K^2=0, obtaining (r+m)dr+=r+dmaKdaK(r_+-m)dr_+=r_+dm-a_Kda_K. Substitute that relation into dAdA, with M=c2m/GNM=c^2m/G_N and J=c3maK/GNJ=c^3ma_K/G_N. The two terms on the right reduce to c4dm/GN=d(Mc2)c^4dm/G_N=d(Mc^2). This family calculation verifies the coefficients; the general horizon mechanics law has broader hypotheses. With TH=κsg/(2πckB)T_H=\hbar\kappa_{\rm sg}/(2\pi c k_B) and SBH=kBAc3/(4GN)S_{\rm BH}=k_BAc^3/(4\hbar G_N), the area term equals THdSBHT_HdS_{\rm BH}. Surface gravity is constant on an equilibrium horizon under the zeroth law’s assumptions. The classical second law is area increase; quantum evaporation calls for generalized entropy instead. The various third-law formulations need additional qualifications and are not needed for this derivation. Wald’s account of the laws and their assumptions.

Two areas constrain one remnantTwo equal initial Schwarzschild horizons are compared with a limiting final circle whose area equals their sum. The equality illustration is the bound’s limiting case, not an achievable merger prediction. Initial binding energy is neglected, and all holes are assumed nonspinning; Kerr horizons require a different area formula.37 / TWO AREAS CONSTRAIN ONE REMNANTArea increase is a bound, not an efficiency predictionAn ideal comparison of nonspinning black holes.Minimum final mass from area alone
37 /
Two areas constrain one remnant. The equality illustration is the bound’s limiting case, not an achievable merger prediction. Initial binding energy is neglected, and all holes are assumed nonspinning; Kerr horizons require a different area formula.

22.10 The black-hole information question#

Classical uncertainty means we do not know which state a system has. Quantum theory also has entanglement: two subsystems can have a definite joint state even when neither has a definite pure state on its own. A pure state describes the complete quantum state; a mixed state describes uncertainty or the reduced description of a subsystem. Unitary evolution is the reversible state evolution of an isolated quantum system in ordinary quantum mechanics. These definitions are enough to state the puzzle, though not to reproduce a quantum-field calculation.

In the leading semiclassical account of a collapsing black hole, outgoing radiation is entangled with degrees of freedom behind the horizon. An observer with access only to the exterior describes approximately thermal radiation, modified by propagation through the surrounding geometry. A thermal-looking spectrum alone does not prove that all correlations are absent.

The tension appears when we combine several claims: a pure initial state, complete evaporation with no remaining hidden system, a final radiation state with irretrievably lost correlations, and unitary evolution of the entire isolated process. Those claims cannot all hold. The classical no-hair description of a stationary exterior is not, by itself, a proof that a quantum state has no microscopic information.

The entropy of radiation expected in a unitary evaporation rises while the radiation is entangled with the remaining hole, then eventually falls to zero if all that remains is the final pure radiation state. This qualitative rise-and-fall behavior is called the Page curve. In specified semiclassical models, calculations of radiation entropy include an island: an interior region whose field correlations contribute to the entropy assigned to the radiation. The calculation varies candidate island boundaries to make the generalized entropy stationary and selects the smallest admissible value. The no-island candidate can dominate early and an island candidate later, producing a Page-shaped curve. The underlying gravitational integral sums over candidate field and geometry histories weighted by their action. Its stationary contributions are called saddle points; additional contributions of this kind produce the island prescription in these models. These results concern specified quantum-gravitational calculations; they do not constitute direct measurements of astrophysical evaporation or a microscopic account valid for every black hole. Almheiri and collaborators’ review.

The boundary of the book is visible here. The classical Einstein equation alone cannot decide how quantum information is recovered. It supplies the geometry in which the question becomes sharp.

The information question has a shapeA schematic entropy curve rises and returns to zero, while a dashed comparison curve continues to rise. Axes are qualitative and the curves are not a quantitative evaporation solution. A final pure radiation state has zero fine-grained entropy for the whole radiation system; individual portions can remain mixed.38 / THE INFORMATION QUESTION HAS A SHAPEfraction of evaporation elapsedradiation entropy (schematic)Unitary expectationRise, then fall as information is recovered.Uncorrected semiclassical trendKeeps growing in the leading approximation.A thermal-looking spectrum can stillcontain correlations between quanta.
38 /
The information question has a shape. Axes are qualitative and the curves are not a quantitative evaporation solution. A final pure radiation state has zero fine-grained entropy for the whole radiation system; individual portions can remain mixed.

The idea to keep

Classical geometry plus energy and global assumptions yields focusing theorems. Hawking temperature needs quantum-field input in addition to GR.

Why is a finite-time caustic not automatically a spacetime singularity?

Even converging straight paths in flat spacetime form a caustic. Singularity theorems need additional global hypotheses to infer geodesic incompleteness.

Figure detail

Scroll to explore at full resolution. Colors follow your reading theme.