Black holes

Notes on Black holes (Draft)

“The black holes of nature are the most perfect macroscopic objects there are in the universe: the only elements in their construction are our concepts of space and time. And since the general theory of relativity provides only a single unique family of solutions for their descriptions, they are the simplest objects as well.”

— S. Chandrasekhar, The Mathematical Theory of Black Holes (1983)

Abstract

Black holes, are considered here with short notes. Sections include a sampling of the literature concerning: 1. The general theory of black holes, 2. Generic behavior of the black hole event horizon, 3. Scattering from black hole sources, and 4. Gravitational radiation from black hole collisions.

Introduction

These notes collect and organize the literature on black hole physics as it connects to computation. The subject sits at the intersection of several mathematical traditions: the global geometric methods of Penrose and Hawking, the perturbation theory developed around the Kerr solution by Teukolsky and Chandrasekhar, the membrane paradigm of Thorne and collaborators, and the numerical relativity program that culminated in the 2005 breakthrough simulations of binary black hole coalescence. What unifies these traditions is the event horizon—a null hypersurface whose geometry, dynamics, and response to perturbation encode most of what we mean by “black hole physics.”

The notes are organized into four sections reflecting the logical structure of the subject. Section 1 treats the foundational theorems: thermodynamic laws, exact solutions, uniqueness results, and the formation problem. Section 2 addresses the event horizon itself—its geometric characterization, the membrane paradigm, null hypersurface formalism, and the dynamics of horizon mergers. Section 3 covers wave propagation and scattering in black hole backgrounds, from Matzner’s pioneering 1968 analysis through modern exact methods. Section 4 treats gravitational radiation: binary dynamics, the self-force problem, data analysis methods, and observational advances.

The bibliography comprises 31 entries (labeled S1–S31) spanning 1968 to 2024. The collection emphasizes mathematical rigor and geometric insight, reflecting the character of black hole physics as one of the most precisely formulated areas of theoretical physics. Several entries serve double duty—Wald’s General Relativity (S3), for instance, provides both the causal-structure machinery needed for precise horizon definitions and the variational framework underlying the first law of black hole mechanics.

A note on scope: these notes focus on classical and semiclassical black hole physics as it bears on computation. Quantum gravity proper—the information paradox, firewall debates, holographic entanglement—lies outside the present scope, though the thermodynamic foundations treated in Section 1 connect directly to those questions. The emphasis throughout is on the mathematical structures that underlie numerical implementation: the 3+1 decomposition, null surface geometry, perturbation theory, and geodesic focusing.

The intellectual lineage of these notes deserves brief comment. The scattering theory of Section 3 traces to Matzner’s 1968 analysis, which predates the golden age of black hole thermodynamics (1971–1974) and even the Penrose process (1969). The perturbation theory underlying gravitational wave template construction was developed by Regge and Wheeler (1957), Zerilli (1970), and Teukolsky (1973), decades before the first detection. The membrane paradigm of Section 2 emerged from astrophysical motivations in the 1980s. The numerical relativity breakthrough came in 2005. And the observational era began on September 14, 2015, with the detection of GW150914. The bibliography thus spans five decades of theoretical development that, taken together, constitutes one of the most complete correspondences between mathematical prediction and experimental confirmation in the history of physics.

The organization into four sections reflects a logical rather than historical ordering. One reads the general theorems first because they define the objects of study; the event horizon section second because the horizon is the defining feature of a black hole; scattering theory third because it provides the link between geometry and observable quantities; and gravitational wave physics last because it is the observational payoff. Within each section, the sources are arranged to support a progressive development of understanding, from foundational to specialized.

1. General Theorems

The mathematical theory of black holes rests on a small number of foundational results: the laws of black hole mechanics, the uniqueness theorems for stationary solutions, the singularity theorems, and the formation problem. These results establish what black holes are (the exact solutions), what they must do (the thermodynamic laws), and when they form (trapped surface criteria). The literature in this section provides the conceptual and mathematical infrastructure on which everything else depends.

1.1 Thermodynamic Foundations

The analogy between black hole mechanics and thermodynamics, established in its definitive form by Bardeen, Carter, and Hawking (1973) [S1], remains one of the most remarkable structural correspondences in theoretical physics. The four laws are:

• Zeroth law. The surface gravity $\kappa$ of a stationary black hole is constant over the event horizon. This is the analog of thermodynamic equilibrium: just as temperature is uniform throughout a body in thermal equilibrium, surface gravity is uniform across the horizon of a stationary black hole.

• First law. For perturbations of a Kerr–Newman black hole, the changes in mass $M$, horizon area $A$, angular momentum $J$, and charge $Q$ satisfy \(dM = \frac{\kappa}{8\pi}\,dA + \Omega\,dJ + \Phi\,dQ,\) where $\Omega$ is the angular velocity and $\Phi$ the electrostatic potential at the horizon. The identification $\kappa/8\pi \leftrightarrow T$ and $A \leftrightarrow S$ gives the Bekenstein–Hawking entropy formula $S = A/4$ (in natural units).

• Second law. The area of the event horizon never decreases in any classical process: $\delta A \geq 0$. Hawking’s area theorem, proved using the Raychaudhuri equation and the null energy condition, is the gravitational analog of the second law of thermodynamics.

• Third law. It is impossible to reduce the surface gravity to zero by any finite sequence of operations—the analog of the unattainability formulation of the third law.

The Hawking (1973) paper [S1] treated these laws as a formal analogy. Bekenstein had already suggested that the analogy was more than formal—that black holes carry genuine entropy proportional to their area. Hawking’s 1974 discovery that quantum effects cause black holes to radiate thermally at temperature $T = \kappa/2\pi$ confirmed Bekenstein’s conjecture and elevated the laws of black hole mechanics to laws of black hole thermodynamics. The present notes do not cover Hawking radiation in detail, but the thermodynamic laws inform several computational diagnostics: monitoring horizon area growth in numerical simulations provides a check on the second law, while the first law relates the mass and spin extracted from gravitational waveforms to horizon geometry.

The biographical overview by Carr et al. (2020) [S9] provides historical context for Hawking’s contributions, including the area theorem, the singularity theorems (with Penrose), and the information paradox. As a survey of one person’s trajectory through the field, it serves as a useful map of how these foundational results developed and interconnected.

1.2 Mathematical Framework

Three texts provide the mathematical infrastructure. They occupy different positions in the difficulty landscape and serve complementary purposes.

Wald (1984) [S3] is the standard graduate-level treatment of general relativity, distinguished from other texts (Misner–Thorne–Wheeler, Carroll, Schutz) by its emphasis on global methods. The later chapters develop Penrose diagrams, causal structure, the Penrose and Hawking singularity theorems, and the uniqueness theorems for stationary black holes. For readers approaching black hole physics through numerical relativity, Wald’s treatment of the initial value formulation (Chapter 10) and his careful definitions of event horizons, apparent horizons, and trapped surfaces (Chapter 12) are essential. The mathematical precision of the definitions matters computationally: the event horizon is a global construct requiring knowledge of the entire future development of the spacetime, while the apparent horizon is quasi-local and can be found on individual time slices—a distinction with direct algorithmic consequences.

Chandrasekhar (1983) [S2] is a different sort of text entirely. Where Wald develops the general framework, Chandrasekhar performs exhaustive calculations for the specific exact solutions. The book treats the Schwarzschild, Reissner–Nordström, Kerr, and Kerr–Newman solutions in turn, developing perturbation theory for each and analyzing their stability properties and quasi-normal mode spectra. The level of calculational detail is extreme by modern standards; Chandrasekhar works through every separation of variables and every radial equation explicitly. For the reader willing to follow the calculations, the reward is complete mastery of the perturbation theory that underlies gravitational wave template construction. The quasi-normal mode spectra computed here are the theoretical predictions tested by LIGO/Virgo ringdown analyses.

Chruściel (2020) [S8] represents the modern mathematical approach. Beginning with global Lorentzian geometry (causal structure, topology, geodesic completeness), the text develops the theory of trapped surfaces, apparent horizons, and uniqueness theorems with full mathematical rigor. The treatment of “numerical black holes”—relating the mathematical definitions to objects that can be located on computational grids—is particularly relevant for the computational focus of these notes. This text occupies the highest position in the difficulty hierarchy: it assumes familiarity with differential geometry at the level of a mathematics (not physics) graduate program.

1.3 Formation of Black Holes

Christodoulou (2008) [S4] addresses a question that had remained open since the singularity theorems of the 1960s: can black holes form dynamically from regular initial data, without symmetry assumptions? The Penrose and Hawking theorems established that singularities form under generic conditions, but singularity formation and black hole formation are logically distinct—a singularity could in principle be naked (visible to distant observers) rather than hidden behind a horizon.

Christodoulou’s monograph proves that trapped surfaces—and hence black holes, given cosmic censorship—can form from the focusing of gravitational waves in vacuum spacetimes. The proof introduces the “short pulse method,” constructing initial data as a thin shell of gravitational radiation and showing that the nonlinear focusing of this pulse generates a trapped surface. The result establishes that Minkowski spacetime, while stable against small perturbations (the Christodoulou–Klainerman theorem), is dynamically unstable to black hole formation under sufficiently strong gravitational wave focusing.

The mathematical machinery is formidable: the proof uses the full Einstein equations in double null coordinates, derives sharp estimates on the Riemann curvature tensor, and tracks the evolution of null geodesic congruences through the focusing region. For the computationally minded reader, the relevant insight is structural: the formation mechanism is nonlinear focusing of null geodesics, governed by the Raychaudhuri equation (treated in Section 2.3 below). The short pulse method has since been extended and simplified by Klainerman, Luk, and Rodnianski, among others.

1.4 The Kerr Solution

Astrophysical black holes rotate. The Kerr solution (1963), describing the unique stationary, axisymmetric, vacuum black hole, is therefore the physically relevant exact solution. Three entries provide complementary perspectives.

Kerr (2008) [S5] is a historical and technical account by Roy Kerr of his discovery of the rotating black hole metric. The paper describes the algebraic approach that led to the solution—the use of Kerr–Schild coordinates, null tetrads, and algebraic classification of the Weyl tensor—and situates the discovery in the context of 1960s general relativity. Kerr–Schild coordinates have found renewed importance in numerical relativity: they provide horizon-penetrating coordinates that avoid the coordinate singularity of Boyer–Lindquist coordinates at the event horizon, making them natural for computational work.

Wiltshire, Visser, and Scott (2009) [S6] edit a comprehensive monograph on Kerr spacetime geometry, treating geodesic structure, the Penrose process (energy extraction from the ergosphere), frame dragging, and the multiple horizons (inner and outer) of the Kerr solution. The edited volume format provides multiple perspectives on the same geometric object, which is useful given the richness of Kerr phenomenology.

Teukolsky (2014) [S7] reviews the perturbation theory of the Kerr metric, centered on the Teukolsky master equation. The remarkable mathematical fact underlying this work is the separability of the wave equation in Kerr spacetime—a property that Chandrasekhar called “miraculous” and that depends on the existence of a hidden symmetry encoded in the Killing tensor. The Teukolsky equation governs perturbations of all spin weights (scalar, electromagnetic, gravitational) in a single formalism. Its solutions yield quasi-normal mode frequencies (the characteristic “ringing” of a perturbed black hole), scattering amplitudes, and stability results. For gravitational wave physics, the Teukolsky formalism provides the theoretical backbone for ringdown waveform templates and black hole spectroscopy.

The separability deserves further comment because it is both mathematically deep and computationally consequential. The Kerr metric admits four constants of motion for geodesics: energy $E$, axial angular momentum $L_z$, rest mass $\mu$, and the Carter constant $Q$—the last being associated with the Killing tensor rather than a Killing vector, and hence a “hidden” symmetry with no simple spacetime interpretation. This four-constant integrability makes the geodesic equations fully separable, reducing the problem to quadratures. For wave equations, the analogous property is that the angular and radial dependences separate completely, reducing partial differential equations to ordinary differential equations (ODEs). This reduction from PDEs to ODEs is what makes analytic progress possible: QNM frequencies can be computed to high precision, scattering cross-sections can be evaluated by numerical integration of ODEs, and the stability of the Kerr solution can be analyzed mode by mode. Without separability, one would be forced to solve the full PDE system numerically—a vastly harder computational problem.

The perturbative stability of the Kerr solution remains, after fifty years, an open mathematical problem. Whiting (1989) proved mode stability (no exponentially growing modes), and recent work by Dafermos, Holzegel, Rodnianski, and others has established quantitative decay estimates for scalar fields on Kerr backgrounds. The full nonlinear stability of Kerr, analogous to the Christodoulou–Klainerman stability proof for Minkowski space, remains unproven—though significant progress has been made by Klainerman and Szeftel.

2. The Event Horizon

The event horizon is the central geometric object in black hole physics. Defined globally as the boundary of the causal past of future null infinity, it is a null hypersurface whose generators are null geodesics that, once they join the horizon, never leave. The geometry of this surface—its expansion, shear, and topology—encodes the dynamical state of the black hole. Computing the event horizon in numerical simulations requires tracking null geodesics backward in time from late in the evolution, a fundamentally global construction that distinguishes it from the quasi-local apparent horizon (the outermost marginally trapped surface on a given time slice).

2.1 The Membrane Paradigm

Thorne, Price, and Macdonald (1986) [S10–S11] introduced the membrane paradigm, a reformulation of black hole physics that treats the event horizon as a two-dimensional membrane endowed with physical properties. The horizon membrane has:

• An electrical surface resistivity of $377\,\Omega$ (the impedance of free space—a coincidence that reflects the universality of the coupling between electromagnetic fields and the horizon).
• A shear viscosity and bulk viscosity, governing the dissipative response of the horizon to tidal perturbations.
• A surface temperature proportional to the surface gravity (the Hawking temperature).

The power of the paradigm is computational: it replaces the global causal structure of the black hole spacetime with boundary conditions on a two-dimensional surface. Electromagnetic fields threading the horizon satisfy equations formally identical to those of a resistive membrane in flat spacetime. Tidal deformations of the horizon are governed by a viscous fluid equation. This reformulation provides physical intuition for phenomena—the Blandford–Znajek mechanism for jet production, tidal heating of the horizon during inspiral—that are difficult to visualize in the full four-dimensional geometry.

For numerical relativity, the membrane paradigm motivates the use of horizon quantities (area, angular momentum, multipole moments) as simulation diagnostics. The horizon’s viscous response also connects to modern developments in the fluid/gravity correspondence, where the long-wavelength dynamics of certain black branes is exactly dual to Navier–Stokes hydrodynamics.

2.2 Null Hypersurface Geometry

The mathematical characterization of the event horizon as a null hypersurface requires machinery beyond the standard 3+1 decomposition of spacelike surfaces. Two references develop this formalism.

Poisson (2004) [S12] provides a modern toolkit for relativists, with detailed treatments of the geometry of null surfaces, junction conditions (Israel formalism), and the boundary conditions appropriate for horizons in numerical simulations. The text is pitched at the working relativist rather than the mathematician, making it more immediately applicable to computational work than the more rigorous treatments of Chruściel.

Gourgoulhon and Jaramillo (2005) [S13] develop a 3+1 perspective on null hypersurface geometry, connecting the isolated horizon formalism of Ashtekar and collaborators with the numerical relativity 3+1 decomposition. This paper bridges two communities: the mathematical relativists who developed quasi-local horizon definitions (isolated horizons, dynamical horizons, trapping horizons) and the numerical relativists who need to locate and characterize horizons on computational grids. The 3+1 reformulation of Damour’s black hole mechanics provides a unified framework for extracting physical quantities (mass, angular momentum, multipole moments) from numerically computed horizons.

The distinction between quasi-local and global horizon definitions is computationally significant. The event horizon requires knowledge of the full future evolution (one must determine which null geodesics eventually escape to infinity), while the apparent horizon—the outermost marginally trapped surface—can be found on a single time slice by solving an elliptic equation. In practice, numerical simulations locate apparent horizons during the evolution (using, e.g., the methods of Thornburg) and construct the event horizon in post-processing by integrating null geodesics backward from late times.

The hierarchy of horizon definitions merits explicit enumeration, as the terminology is a source of persistent confusion. An event horizon is the boundary of the causal past of future null infinity—a teleological definition requiring knowledge of the entire future spacetime. An apparent horizon is the outermost marginally trapped surface on a given spacelike slice, defined locally in time but depending on the choice of foliation. An isolated horizon (Ashtekar et al.) is a null surface satisfying quasi-equilibrium conditions (vanishing expansion, time-independent induced geometry), appropriate for black holes that are not actively accreting or merging. A dynamical horizon (Ashtekar and Krishnan) generalizes the isolated horizon to non-equilibrium situations, defined as a spacelike surface foliated by marginally trapped 2-spheres.

For numerical work, the apparent horizon is the primary real-time diagnostic because it can be found algorithmically on each time slice. The standard method solves the equation $\theta_{(\ell)} = 0$ for the outgoing null expansion on trial surfaces, typically using spectral methods (expanding the surface in spherical harmonics) with Newton’s method iteration. The event horizon, computed in post-processing, provides the true boundary and the correct area for thermodynamic considerations, but at the cost of requiring the full spacetime history.

2.3 The Raychaudhuri Equation

Kar and SenGupta (2007) [S14] provide a comprehensive review of the Raychaudhuri equation, the fundamental evolution equation for geodesic congruences in curved spacetime. For a congruence of null geodesics with tangent vector $k^a$, the Raychaudhuri equation governs the evolution of the expansion $\theta$:

\[\frac{d\theta}{d\lambda} = -\frac{1}{2}\theta^2 - \sigma_{ab}\sigma^{ab} + \omega_{ab}\omega^{ab} - R_{ab}k^a k^b,\]

where $\sigma_{ab}$ is the shear, $\omega_{ab}$ is the twist, and $R_{ab}$ is the Ricci tensor. For a hypersurface-orthogonal congruence (vanishing twist), the null energy condition ($R_{ab}k^a k^b \geq 0$) ensures that $\theta$ decreases monotonically, driving geodesics to focus and form caustics.

This equation is the workhorse of classical black hole physics. Hawking’s area theorem follows from it: the generators of the event horizon form a twist-free null congruence, so the null energy condition guarantees non-decreasing expansion, which in turn guarantees non-decreasing area. The singularity theorems of Penrose and Hawking use the Raychaudhuri equation to show that geodesic focusing leads to geodesic incompleteness. Christodoulou’s formation proof (Section 1.3) is ultimately an analysis of when the nonlinear terms in the Raychaudhuri equation dominate.

Computationally, the Raychaudhuri equation governs the behavior of null geodesic generators that constitute the event horizon. When generators join the horizon (at caustics or “crease” points), the expansion undergoes characteristic behavior that must be tracked accurately in numerical event-horizon finders.

2.4 Merger Dynamics

The geometry of the event horizon during a binary black hole merger is topologically nontrivial. Before merger, there are two disjoint horizons; after merger, there is one. The transition involves a topology change on the horizon surface, occurring through the formation and annihilation of “crease sets”—singularities in the null generators where new generators join the horizon.

Emparan and Martínez (2016) [S15] study event horizon dynamics during mergers, particularly in higher-dimensional spacetimes where the richer topology of horizon cross-sections allows novel phenomena (e.g., “pinch-off” transitions from a single connected horizon to multiple components, the reverse of merger). While the higher-dimensional aspects go beyond the scope of astrophysical computation, the analysis of topology change mechanisms applies equally to the four-dimensional case.

Emparan et al. (2018) [S16] investigate the extreme mass ratio limit, where a small black hole spirals into a much larger one. In this limit the deformation of the large horizon can be treated perturbatively, but the merger itself remains a genuinely nonlinear event. The extreme mass ratio regime is the target of the LISA space-based gravitational wave detector, making horizon dynamics in this regime both mathematically tractable and observationally relevant.

Gadioux (2024) [S17] examines the formation and evolution of crease sets—the loci where new null generators join the event horizon through caustic singularities. These geometric features encode information about the merger history and provide a diagnostic for the dynamical state of the horizon. This recent work connects to the broader program of understanding event horizon topology change, building on earlier work by Siino, Husa and Winicour, and others on the structure of the event horizon during gravitational collapse and merger.

3. Scattering from Black Holes

The scattering of waves from black holes connects the geometric theory of Sections 1–2 to observable quantities: absorption cross-sections, scattering amplitudes, quasi-normal mode frequencies, and (ultimately) gravitational waveforms. The mathematical framework is the theory of wave equations on curved backgrounds—the Klein–Gordon, Maxwell, and linearized Einstein equations in Schwarzschild and Kerr spacetimes. The remarkable separability properties of these equations in the Kerr background (encoded in the Teukolsky equation) make analytic progress possible, while the physical phenomena—glory scattering, superradiance, quasi-normal ringing—provide both intuition and computational benchmarks.

3.1 Foundations

Matzner (1968) [S18] established the scattering cross-section for massless scalar waves by a Schwarzschild black hole, pioneering the analysis of wave-optical phenomena in black hole spacetimes. The key physical insight is that the photon sphere at $r = 3M$ acts as a critical impact parameter: waves with angular momentum near the critical value orbit multiple times before scattering, producing interference effects analogous to glory scattering in optics. The forward scattering amplitude diverges logarithmically (the gravitational analog of Coulomb scattering), and the backward glory produces a characteristic diffraction pattern whose angular scale is set by the photon sphere radius.

This early work is foundational for several reasons. First, it established the partial wave decomposition as the primary computational tool for black hole scattering, expanding the scattered field in spin-weighted spherical harmonics and solving ordinary differential equations for the radial functions. Second, it identified the physically important length scales: the Schwarzschild radius $2M$, the photon sphere $3M$, and the wavelength $\lambda$ of the incident radiation. The dimensionless ratio $M/\lambda$ governs the relative importance of geometric optics ($M/\lambda \gg 1$) versus wave-optical ($M/\lambda \sim 1$) effects. Third, it motivated the development of more sophisticated approaches—WKB methods, complex angular momentum techniques, exact analytic solutions—that the subsequent literature pursues.

3.2 Spectral Theory

Persides (1973) [S19–S20] investigated the spectral properties of the Laplace operator in Schwarzschild spacetime and derived the radial wave equation governing scattering. These two papers establish the mathematical framework—the function spaces, boundary conditions, and spectral properties—for the scattering problem as a well-posed mathematical question. The radial equation, after separation of angular variables, takes the form of a one-dimensional Schrödinger equation with an effective potential that depends on the black hole mass, the angular momentum quantum number $\ell$, and the spin of the field. The effective potential has a single maximum near the photon sphere and vanishes at both the horizon ($r \to 2M$) and spatial infinity ($r \to \infty$), giving a standard scattering problem with well-defined transmission and reflection coefficients.

The spectral theory perspective is important for understanding quasi-normal modes (QNMs)—the damped oscillations of a perturbed black hole. QNMs are not normal modes in the usual sense (the operator is not self-adjoint due to radiation boundary conditions at the horizon and infinity), but rather resonances—poles of the scattering matrix in the complex frequency plane. The machinery Persides develops is the starting point for the QNM calculations that are central to gravitational wave ringdown analysis.

3.3 Comprehensive Treatments

Futterman, Handler, and Matzner (1988) [S21] is the comprehensive monograph on scattering from black holes, covering scalar, electromagnetic, and gravitational waves in Schwarzschild and Kerr backgrounds. The text develops the partial wave formalism systematically, treats absorption cross-sections (how much energy a black hole captures), glory scattering (the backward-directed interference pattern), and superradiant amplification (the wave analog of the Penrose process, where waves scattered from a rotating black hole can extract rotational energy). This book is the primary reference for anyone implementing scattering calculations.

Bezerra et al. (2013) [S22] extend the scattering formalism to charged (Reissner–Nordström) and rotating (Kerr) backgrounds. The separation of the Klein–Gordon equation in Kerr–Newman spacetime introduces the angular eigenvalue problem (spheroidal harmonics) as an additional computational element.

Li et al. (2022) [S23] present exact analytic solutions for certain black hole scattering problems, advancing beyond the perturbative and numerical treatments of earlier work. Exact solutions provide benchmarks for numerical codes and can reveal qualitative features—such as the analytic structure of the scattering amplitude in the complex angular momentum plane—that approximate methods obscure.

4. Gravitational Waves

The detection of gravitational waves from binary black hole mergers by LIGO in 2015 transformed black hole physics from a theoretical discipline into an observational science. The theoretical infrastructure underlying gravitational wave astronomy includes post-Newtonian approximations for the inspiral phase, numerical relativity for the merger, and perturbation theory for the ringdown. The self-force program provides a systematic framework for extreme mass ratio systems. The entries in this section cover the theoretical, computational, and observational aspects of gravitational radiation from black hole systems.

4.1 Binary Dynamics

Trimble and Thorne (2018) [S24] review the physics of binary systems and their gravitational wave signatures. The conceptual framework for understanding gravitational waves from binaries involves three regimes, each requiring different analytical and computational tools:

The inspiral phase, where the binary components orbit at separations much larger than the Schwarzschild radii, is well described by post-Newtonian (PN) theory—a systematic expansion in powers of $v/c$ beyond Newtonian gravity. The leading gravitational wave emission is quadrupolar, with the Einstein quadrupole formula giving the luminosity $L = (32/5)(M\omega)^{10/3}$ for a circular binary with total mass $M$ and orbital frequency $\omega$ (in geometric units).

The merger phase, where the two bodies plunge together and their horizons coalesce, can only be computed by solving the full Einstein equations numerically. This phase produces the peak gravitational wave amplitude and encodes information about the nonlinear dynamics of general relativity that is inaccessible to perturbative methods. The 2005 breakthrough by Pretorius, using generalized harmonic coordinates, and the independent 2006 successes by Campanelli et al. (moving punctures with BSSN formalism) and Baker et al. (similar approach), demonstrated that the merger phase is computationally tractable and produces waveforms in remarkable agreement with the post-Newtonian and perturbation theory predictions in the appropriate limits. The merger waveform is dominated by the $(\ell,m) = (2,2)$ mode and lasts only a few cycles, but it carries the most luminous gravitational wave emission—during the peak, GW150914 radiated energy at a rate exceeding the combined luminosity of all stars in the observable universe.

The ringdown phase, where the newly formed black hole settles to a Kerr solution by radiating quasi-normal modes, is described by black hole perturbation theory. The QNM frequencies depend only on the mass and spin of the final black hole (by the no-hair theorem), making ringdown observations a test of general relativity. The dominant mode is the fundamental $(\ell,m) = (2,2)$ mode, with frequency $f_{220} \approx 1.525\,(100\,M_\odot/M)\,\text{kHz}$ and damping time $\tau_{220} \approx 0.856\,(M/100\,M_\odot)\,\text{ms}$ for a Kerr black hole with spin $a/M = 0.67$. Detecting multiple QNM overtones would enable black hole spectroscopy—a direct test of the no-hair theorem and a probe of the strong-field regime of general relativity.

4.2 Self-Force and Radiation Reaction

The self-force program addresses the problem of motion and radiation reaction for a compact object orbiting a much larger black hole—the extreme mass ratio inspiral (EMRI) scenario. The physical question is: how does the gravitational field of a small body affect its own trajectory? The answer involves regularizing the divergent self-field of a point particle in curved spacetime and extracting the finite “radiation reaction” force that drives the inspiral.

Quinn (2000) [S25] analyzes the self-force for a scalar charge in curved spacetime, which serves as a technically simpler model problem for the gravitational case. The scalar self-force problem shares the essential mathematical structure—a Green’s function that must be decomposed into a singular “direct” part (to be subtracted) and a regular “tail” part (which contributes to the self-force)—without the gauge complications of the gravitational case.

de Felice and Bini (2010) [S26] develop the formalism for classical measurements in curved spacetimes, addressing how local observers characterize tidal fields, precession effects, and radiative forces. This work provides the conceptual framework for interpreting self-force results in terms of physically measurable quantities.

Poisson, Pound, and Vega (2011) [S27] is the comprehensive Living Reviews article on the motion of point particles in curved spacetime. The article develops the gravitational self-force through matched asymptotic expansions—treating the small body’s neighborhood with its own strong-field solution and matching to the external spacetime treated as a perturbation of the large black hole’s geometry. The resulting “MiSaTaQuWa” equation (named for Mino, Sasaki, Tanaka, Quinn, and Wald) gives the regularized equation of motion. This formalism underlies the waveform models being developed for LISA, the space-based gravitational wave detector that will observe EMRIs in the millihertz band.

The self-force problem illustrates a recurring theme in these notes: the interplay between exact mathematical structure and computational implementation. The formal solution involves the retarded Green’s function of the linearized Einstein equation, which can be decomposed into a “direct” part (supported on the past light cone) and a “tail” part (supported inside the light cone, arising from backscattering off the spacetime curvature). The direct part diverges at the particle’s location and must be regularized; the tail part is finite and provides the physical self-force. Several regularization schemes have been developed—mode-sum regularization (Barack and Ori), effective source methods (Vega and Detweiler), and puncture schemes—each with distinct computational trade-offs. The first-order self-force is now computed routinely for generic orbits in Kerr spacetime, and second-order calculations (needed for LISA accuracy requirements) are an active area of development.

The gravitational self-force connects to the broader theme of radiation reaction in field theories. The electromagnetic analog—the Abraham–Lorentz–Dirac equation for the self-force on a charged particle—has been known since the early twentieth century and exhibits pathologies (runaway solutions, preacceleration) that the gravitational case largely avoids due to the equivalence principle. The scalar model problem treated by Quinn [S25] shares the mathematical structure of the gravitational case without the gauge complications, making it a valuable testing ground for regularization methods and numerical algorithms.

4.3 Detection and Data Analysis

Jaranowski and Królak (2023) [S28] survey the statistical and computational methods for analyzing gravitational wave signals in detector data. The dominant technique is matched filtering: cross-correlating the detector output with a bank of theoretical waveform templates and searching for statistically significant matches. The method requires accurate waveform models spanning the parameter space (masses, spins, sky location, orientation), which in turn requires the inspiral-merger-ringdown waveform models produced by combining post-Newtonian theory, numerical relativity, and perturbation theory. The review covers the Gaussian noise case in detail, treating signal detection as a hypothesis testing problem within the Neyman–Pearson framework, and discusses parameter estimation using Bayesian methods.

Romano and Cornish (2017) [S29] review a complementary detection problem: the stochastic gravitational wave background—the superposition of unresolved signals from many sources. The stochastic background has multiple astrophysical components (unresolved compact binaries, supernovae, neutron star mergers) and potentially cosmological components (inflation, phase transitions, cosmic strings). Detection requires cross-correlating the outputs of multiple detectors, exploiting the fact that the gravitational wave signal is correlated between detectors while the noise is not. The overlap reduction function, which encodes the detector geometry and the antenna pattern, determines the sensitivity of a given detector pair to the stochastic background.

4.4 Observational Breakthroughs

LIGO Scientific Collaboration et al. (2016) [S30] announced the first direct detection of gravitational waves: the event GW150914, observed on September 14, 2015. The signal matched the waveform predicted by general relativity for the inspiral, merger, and ringdown of two black holes with masses approximately $36\,M_\odot$ and $29\,M_\odot$ at a luminosity distance of roughly 410 Mpc. The final black hole had mass approximately $62\,M_\odot$ with spin parameter $a/M \approx 0.67$, meaning roughly $3\,M_\odot c^2$ of energy was radiated as gravitational waves in less than a second.

The detection confirmed several theoretical predictions simultaneously: the existence of gravitational waves (predicted by Einstein in 1916), the existence of stellar-mass binary black holes (whose formation channels had been debated), and the accuracy of numerical relativity waveform predictions for the merger phase. Subsequent observing runs have detected dozens of binary black hole mergers, enabling population studies and increasingly precise tests of general relativity.

The significance of GW150914 for the program of these notes cannot be overstated. The matched filter analysis that identified the signal required a bank of waveform templates produced by combining post-Newtonian theory (inspiral), numerical relativity (merger), and black hole perturbation theory (ringdown)—precisely the three theoretical frameworks treated in Sections 1, 3, and 4 of these notes. The parameter estimation that determined the source masses and spin used Bayesian methods with the same template families. And the test of general relativity—comparing the observed ringdown frequency with the prediction from the no-hair theorem—relied on the QNM calculations rooted in the Teukolsky formalism of Section 1.4 and the scattering theory of Section 3.

The gravitational wave catalog has grown rapidly since GW150914. The first three observing runs (O1, O2, O3) of LIGO and Virgo detected approximately 90 compact binary coalescences, predominantly binary black hole mergers but also binary neutron star mergers (GW170817, with its electromagnetic counterpart) and neutron star–black hole mergers. The fourth observing run (O4), begun in 2023, has further expanded the catalog. The population of observed binary black holes has yielded constraints on formation channels (isolated binary evolution versus dynamical formation in dense stellar environments), the black hole mass spectrum (including evidence for features at masses where pair-instability supernovae are expected to suppress black hole formation), and the distribution of black hole spins.

Pretorius (2023) [S31] provides a comprehensive survey of gravitational wave physics by one of the pioneers of numerical relativity. Pretorius achieved the first stable simulation of a binary black hole merger and gravitational wave emission in 2005—a breakthrough that, together with the independent work of Campanelli et al. and Baker et al. in 2006, solved a problem that had resisted numerical attack for decades. The review covers the current state of gravitational wave observations, theoretical developments, and prospects for next-generation detectors (Einstein Telescope, Cosmic Explorer, LISA).

6 Reading Paths

The collection supports several reading paths depending on the reader’s goals:

Mathematical foundations. S3 (Wald) → S2 (Chandrasekhar) → S8 (Chruściel) → S4 (Christodoulou). This path develops the rigorous framework, progressing from the standard graduate treatment through exhaustive calculations of exact solutions to modern geometric methods and the formation problem.

Computational black hole physics. S3 (Wald, Chapters 10, 12) → S12 (Poisson toolkit) → S13 (Gourgoulhon–Jaramillo 3+1 null surfaces) → S14 (Raychaudhuri equation) → S10–11 (membrane paradigm). This path develops the formalism needed for numerical horizon finding and characterization.

Gravitational wave physics. S7 (Teukolsky) → S21 (Futterman–Handler–Matzner scattering) → S27 (Poisson–Pound–Vega self-force) → S28 (Jaranowski–Królak data analysis) → S30 (LIGO detection). This path follows the chain from perturbation theory through waveform modeling to observation.

The Kerr solution. S5 (Kerr discovery) → S6 (Wiltshire et al. monograph) → S7 (Teukolsky perturbations) → S2 (Chandrasekhar, Kerr chapters). This path provides a complete picture of the astrophysically relevant black hole solution.

Bibliography

General theorems

ID	Link	Notes
S1	Hawking “Blackhole mechanics” 1973	Landmark work establishing the four laws of black hole mechanics and their deep analogy to thermodynamics.
S2	Chandrasekhar “Mathematical Theory of Blackholes” 1983	Rigorous mathematical treatment of exact black-hole solutions and their stability properties in general relativity.
S3	Wald “General Relativity” 1984	Authoritative graduate-level text covering the foundations of general relativity and the global theorems governing black holes.
S4	Christodoulou “Formation of Black Holes in General Relativity” 2008	Detailed proof demonstrating gravitational collapse in vacuum spacetimes leads to black-hole formation.
S5	Kerr “Discovering the solution” 2008	Historical and technical overview of Roy Kerr’s derivation of the rotating black-hole metric.
S6	Wiltshire et al “The Kerr Spacetime” 2009	Comprehensive monograph analyzing the geometry and physical properties of the Kerr solution.
S7	Teukolsky “The Kerr Metric” 2014	Survey of perturbation theory and separability in Kerr spacetime for wave propagation analyses.
S8	Crusciel “Geometry of Black holes” 2020	Mathematical exposition of black-hole geometry emphasizing horizon uniqueness theorems.
S9	Carr et al “Stephen William Hawking CH CBE” 2020	Biographical and scientific overview summarizing Hawking’s key contributions to black-hole physics.

The event horizon

ID	Link	Notes
S10	Thorne et al “Membrane paradigm” 1986	Introduces the membrane paradigm, modeling the event horizon as a two-dimensional viscous surface with physical properties.
S11	Thorne et al “Membrane paradigm” 1986	Extended discussion of horizon electrodynamics and fluid analogies within the membrane framework.
S12	Poisson “A Relativist’s Toolkit” 2004	A modern text including detailed consideration of null surfaces, with a detailed derivation of horizon boundary conditions.
S13	Gourgoulhon and Jaramillo “3+1 perspective on null surfaces” 2005	Applies a 3+1 decomposition to characterize the geometry and evolution of null hypersurfaces such as event horizons.
S14	Kar and SenGupta “Raychaudhuri equations” 2007	Derives Raychaudhuri’s equation in various geometric settings to analyze horizon focusing and caustic formation.
S15	Emparan and Martínez “Black hole fusion” 2016	Studies the dynamics of event horizons during black-hole mergers in higher-dimensional spacetimes.
S16	Emparan et al “Black hole fusion in the extreme mass ratio” 2018	Numerical investigation of horizon deformation in extreme-mass-ratio inspirals.
S17	Gadioux “Evolution of creases on the event horizon” 2024	Examines the formation and evolution of crease sets on dynamical horizons during non-linear interactions.

Scattering from black holes

ID	Link	Notes
S18	Matzner “Scattering from Blackholes” 1968	Early analysis of massless scalar wave scattering by a Schwarzschild black hole, establishing cross-sections.
S19	Persides “Laplacian on Schwarzschild” 1973	Investigates spectral properties of the Laplace operator on the Schwarzschild background.
S20	Persides “Radial Schwarzschild” 1973	Derives and analyzes the radial wave equation governing scattering in Schwarzschild spacetime.
S21	Futterman et al “Scattering from black holes” 1988	Comprehensive monograph covering quantum and classical wave scattering by black holes across spin fields.
S22	Bezzerra “Scattering from charge and rotating black holes” 2013	Extends scattering theory to charged (Reissner–Nordström) and rotating (Kerr) black-hole backgrounds.
S23	Li et al “Exact Scattering” 2022	Presents exact analytic solutions for wave scattering in certain black-hole geometries.

Gravitational waves

ID	Link	Notes
S24	Trimble and Thorne “Binaries” 2018	Review of binary black-hole dynamics and their astrophysical and gravitational-wave implications.
S25	Quinn “Radiating scalar particle” 2000	Analyzes the self-force and radiation reaction of a scalar charge in curved spacetime.
S26	de Felice and Bini “Classical Measurement in Curved Spacetime” 2010	Develops a formalism for how classical observers measure tidal and radiative effects in curved backgrounds.
S27	Posisson “The radiating particle” 2011	Reviews motion and self-force of point particles emitting gravitational radiation in general relativity.
S28	Jaranowski and Królak “Gravitational-wave data” 2023	Survey of statistical and computational methods for analyzing and interpreting gravitational-wave signals.
S29	Romano and Cornish “Stochastic gravitational wave backgrouns” 2017	Reviews theoretical models and detection strategies for stochastic gravitational-wave backgrounds.
S30	LIGO et al “Gravitational Waves” 2020	Landmark LIGO/Virgo collaboration report presenting the first direct detection of gravitational waves from a binary merger.
S31	Pretorius “Survey of Gravitational Waves” 2023	Comprehensive overview of numerical relativity simulations and observational prospects in gravitational-wave astronomy.