Horizons, Charges, and Gravitational-Wave Extraction
Scrollable formula. Use Left and Right arrow keys to reveal hidden content.
Scrollable table. Use Left and Right arrow keys to reveal additional columns.
Scrollable figure. Use Left and Right arrow keys to inspect hidden content.
This chapter follows numerical information from a strong-field spatial slice to radiation at future null infinity. It treats horizon location, quasi-local charges, curvature modes, radial extrapolation, characteristic extraction, fluxes, and memory as one connected measurement problem. The organizing principle is that every reported quantity carries a surface, frame, normalization, resolution, and independent consistency check.
Helpful background. Formulations of General Relativity supplies Cauchy and characteristic decompositions; Differential Equations for Relativity supplies nonlinear solves and evolution methods; Black-Hole Geometry and Mechanics defines the horizon structures; Gravitational Radiation and Waveforms defines asymptotic observables; and Computational Relativity supplies the volume-wide numerical framework.
Enter this chapter
Section titled “Enter this chapter”Begin with Apparent-Horizon Finding, where the vanishing outgoing expansion becomes a nonlinear surface equation and multiple branches must be searched deliberately. Continue to Quasi-Local Mass, Spin, and Multipoles, which makes the approximate-symmetry and normalization dependence of horizon charges explicit.
The radiation path starts with Curvature-Based Wave Extraction, including the tetrad, harmonic, and finite-sphere conventions behind . Finite-Radius Extrapolation asks when a stable inverse-radius fit can approximate null infinity. Cauchy–Characteristic Extraction replaces that fit with a null evolution and its own worldtube, initial-data, and asymptotic-frame errors. Finally, Fluxes, Memory, and Consistency Checks tests whether radiation and remnant quantities close the appropriate charge balances.
The continuum definitions remain with the geometry and radiation volumes. Here the task is computational: translate those definitions into finite algorithms and show how an implementation can be wrong even when a residual is small.
What each measurement tells you
Section titled “What each measurement tells you”An apparent horizon answers a question about a particular spatial slice: where does the outgoing light-front area stop increasing? Its area and spin characterize a two-surface near the black hole. A waveform answers a different question: how does the asymptotic transverse geometry change with retarded time? Comparing their charges through a balance law connects the near and far measurements, but does not make their gauges or errors identical.
A useful first calculation is the time-symmetric Schwarzschild sphere: derive its expansion, find the isotropic radius , and recover area . A useful first radiative calculation is one sinusoidal mode: differentiate strain to news and curvature, then compute its positive energy flux. The leaf pages work through both before discussing production data.
One chain with distinct uncertainties
Section titled “One chain with distinct uncertainties”The chapter has four coupled numerical layers:
| Layer | Primary result | Refinement variable | Discriminating check |
|---|---|---|---|
| Surface | Marginal surface and area | Volume and angular resolution | Kerr geometry and independent expansion residual |
| Charge | Spin, mass, multipoles | Surface representation and eigensolver tolerance | Coordinate rotation and late Kerr limit |
| Waveform | Curvature, news, strain modes | Cauchy, extraction, and characteristic grids | Linear wave and gauge Schwarzschild tests |
| Balance | Fluxes, recoil, memory | Mode cutoff and time window | ADM–radiation–remnant and BMS relations |
Do not collapse these layers into one “waveform error.” A horizon solver can converge while its parent volume grid remains under-resolved. A finite-radius fit can be stable while the tetrad mixes Coulomb curvature into the radiative scalar. A characteristic evolution can converge while worldtube reconstruction remains biased. A flux budget can close accidentally through correlated errors or a tuned integration constant.
Bishop and Rezzolla organize the major extraction methods and their relationships in Bishop and Rezzolla 2016, §§2–8. Thornburg supplies the corresponding algorithmic taxonomy for horizons in Thornburg 2007, §§3.1–3.3. Together they motivate a chapter-wide rule: change both numerical scale and mathematical route whenever feasible.
Recurring benchmarks
Section titled “Recurring benchmarks”Three benchmark families connect the pages. First, Schwarzschild and Kerr slices have known horizon areas, stationary charges, and—under an asymptotically appropriate tetrad—no outgoing radiation. Coordinate shifts and rotations test whether an implementation confuses coordinate geometry with an invariant.
Second, a linear Teukolsky or perturbed-Schwarzschild wave supplies a known radiative solution. Extract it at several radii and resolutions, then compare finite-radius fitting with characteristic propagation. In a nonlinear characteristic code compared against a linear solution, changing wave amplitude distinguishes expected quadratic modeling error from truncation error.
Third, a resolved binary-black-hole evolution connects initial ADM charges, individual and common horizons, outgoing modes, and final remnant charges. It has no exact full solution, so it is a consistency benchmark rather than an analytic validation case. Its value comes from crossing methods: independent surface residuals, multiple charge definitions, radius and worldtube variation, and flux closure.
For every family, pin the input checksum, coordinate and tetrad conventions, precision, solver revisions, tolerances, extraction radii, harmonic cutoff, and time interval. A future rerun should emit machine-readable residuals and refinement tables. The pages specify such interfaces but do not claim that the calculations have been executed.
Reading a claimed result
Section titled “Reading a claimed result”Ask five questions in order. What geometric surface or asymptotic cut supports the quantity? Which normal, tetrad, harmonic, and Bondi frame define it? Which resolution or tolerance sequence demonstrates a stable regime? Which systematic variation changes the mathematical route rather than merely redrawing the same output? Which failure criterion would cause the claim to be withdrawn?
The most important limitation is slicing and frame dependence. Apparent horizons are slice dependent; approximate spins depend on a chosen vector field; finite-radius modes depend on extraction spheres and tetrads; angular momentum and memory depend on asymptotic frame and integration choices. These dependencies do not make the quantities useless. They specify what must accompany them.
Review the chapter
Section titled “Review the chapter”A binary calculation reports a convergent common-horizon mass and a convergent CCE energy flux. Their energy balance misses closure by , much larger than either quoted truncation error. The discrepancy is unchanged under characteristic refinement but falls when the Cauchy grid is refined. Where should the investigation begin?
Solution
Begin upstream with the shared Cauchy solution and its interpolation to both horizon and worldtube. Characteristic refinement cannot remove an error already present in its boundary data, and separate convergence of two derived quantities does not expose their correlated Cauchy bias. Refine the evolution, repeat the horizon and worldtube measurements together, inspect constraints and outer-boundary effects, and retain the balance residual as an unresolved systematic until the common refinement explains it.