The Bondi--Sachs gauge, BMS frames, and memory in black hole perturbation theory

Resolving the Infrared Divergence in Black Hole Waveforms

When scientists model gravitational waves from merging black holes, they often run into mathematical errors called "divergences" at large distances. These errors occur because the coordinate systems used to describe spacetime can become pathological. Essentially, they become "broken rulers" as one moves toward the edge of the universe. This paper provides a new mathematical ruler, the Bondi–Sachs gauge, which remains accurate even at great distances. This helps to correctly include subtle effects like gravitational-wave memory.

The Breakdown of Near-Zone Gauges

Current gravitational-wave modeling relies heavily on black hole perturbation theory (BHPT). This method treats the spacetime of a binary system as a smooth background slightly disturbed by orbiting masses. To make these models useful for next-generation detectors like LISA, researchers must push calculations to second-order precision. However, a significant problem emerges at this higher order: infrared divergences.

These divergences are not physical phenomena but mathematical artifacts. They are caused by a mismatch between the "near-zone" (the immediate vicinity of the orbiting bodies) and the "far-zone" (the distant region where gravitational waves propagate). In traditional frameworks, such as the Lorenz gauge, the mathematical descriptions of the gravitational field fall off too slowly at large distances. When researchers try to integrate these fields to find the total waveform, the math "blows up." This results in infinite values that prevent meaningful physical predictions.

Specifically, in second-order self-force calculations, current approaches struggle to reconcile local orbital dynamics with global spacetime behavior. This often necessitates complex "matching" procedures to bridge the gap.

An Iterative Path to the Bondi–Sachs Gauge

To solve this, the authors propose a framework centered on the Bondi–Sachs (BS) formalism. Unlike standard coordinates, the BS formalism uses a coordinate system adapted to outgoing null cones. These are surfaces formed by light rays traveling outward from the source. This choice ensures the "ruler" used to measure distance is linked to the path of the gravitational waves themselves.

The authors implement their solution through an iterative three-step procedure:

1. Initial Solution: First, the researcher solves the linearized Einstein field equations in a "convenient" gauge. This gauge is easy to work with mathematically, even if it is poorly behaved at infinity.
2. Gauge Transformation: Next, they apply a mathematical transformation to shift this solution into the Bondi–Sachs gauge. This transformation ensures metric perturbations decay rapidly enough at large distances to prevent mathematical singularities.
3. Higher-Order Integration: Finally, with the first-order field now in the well-behaved BS gauge, the second-order equations can be solved directly.

Crucially, this process allows the researchers to fix the Bondi–Metzner–Sachs (BMS) frame. The BMS group represents the residual symmetry of spacetime at infinity. It includes standard movements like rotations and boosts. It also includes "supertranslations" (an infinite-dimensional group of shifts that change the perceived arrival time of waves across the sky). By fixing these symmetries, the authors ensure that resulting waveforms are unambiguous and comparable across different models.

Evading Divergences and Capturing Memory

The effectiveness of this approach is demonstrated by its ability to handle difficult second-order physics. The paper reports that by adopting the BS gauge at first order, the source terms for the second-order Teukolsky equation become well-behaved. The Teukolsky equation is a fundamental tool used to describe gravitational radiation.

While the Lorenz gauge produces a source term that decays like $1/r^2$, the BS gauge ensures the source decays much more rapidly, at rates such as $1/r^4$. This rapid decay prevents the mathematical integrals from diverging.

Furthermore, the authors show that this framework naturally incorporates "gravitational-wave memory." Memory is a non-linear effect where a passing wave leaves a permanent displacement in the fabric of spacetime. In older "forgetful gauges" used in self-force theory, this memory is lost in near-zone calculations. It must then be manually added back "after the fact."

The authors show that in the BS gauge, memory is not forgotten. Instead, it is manifested as "soft hair"—slowly evolving degrees of freedom that emerge during the transformation. This also accounts for "memory distortion." This is a subtle effect where the presence of memory slightly alters the shape of the oscillating gravitational waves.

Limitations of the Framework

Despite its strengths, the framework has two primary constraints:

1. First-Order Frame Fixing: The current scheme fixes the BMS frame only to the first perturbative order. This is sufficient to ensure that second-order Weyl scalars (mathematical objects used to extract the waveform from spacetime curvature) are gauge-invariant. However, it may not capture all higher-order complexities.
2. Residual Isometries: Certain symmetries of the background Kerr spacetime remain unfixed. These include time translations and rotations around the axis of symmetry. Because these are inherent properties of the black hole, they cannot be removed by changing the gauge. Researchers must still manually align the time and phase of two different waveforms when comparing them.

The Verdict

For practitioners of black hole perturbation theory, this paper is a vital upgrade. It provides a systematic way to move from "forgetful" gauges to the Bondi–Sachs gauge. This eliminates the need for the cumbersome matching procedures previously required to handle far-zone effects.

The framework is highly actionable. Code is reportedly available; see the paper for the canonical link. For anyone working on second-order waveforms for LISA or next-generation detectors, this iterative scheme offers a robust path to accuracy. The transition from treating memory as an afterthought to treating it as an intrinsic part of the gauge is a significant step toward high-fidelity modeling.