⚛️ general relativity

Analytical Time-Domain Waveforms for Black-Hole Scattering in the Post-Minkowskian Expansion

This paper presents the first analytic computation of time-domain gravitational waveform multipoles for hyperbolic black-hole scattering at leading post-Minkowskian order to all orders in velocity, utilizing cut loop integrals and dimensional regularization to derive logarithmic expressions that facilitate the integration of post-Minkowskian fluxes into binary merger waveform models.

Original authors: Stefano De Angelis, Vincent F. He, Julio Parra-Martinez

Published 2026-10-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Stefano De Angelis, Vincent F. He, Julio Parra-Martinez

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Gravity is the invisible architecture of the universe, the force that holds planets in orbit and bends the path of light. For decades, scientists have listened to the ripples in this fabric, known as gravitational waves, to understand the violent collisions of black holes. To interpret these signals, researchers rely on mathematical models that predict what the waves should look like. These models are built on two main pillars: one that works well when objects move slowly and gently, and another that handles the extreme speeds and intense gravity found near black holes. While the slow-motion version is very precise, the high-speed version has historically been difficult to translate into the specific, time-based shapes needed to match real-world observations. Without a clear picture of how these waves evolve moment by moment at high speeds, our ability to decode the most energetic events in the cosmos remains incomplete.

A team of physicists has now bridged this gap by calculating the exact shape of gravitational waves produced when two black holes fly past each other at incredible speeds without colliding. This scenario, known as hyperbolic scattering, is the high-speed cousin of the merging black holes that detectors like LIGO and Virgo usually observe. The researchers focused on the very first stage of this interaction, where the gravitational pull is weak enough to be treated as a small disturbance, but the speed of the black holes is arbitrary, meaning it can be anything from a slow drift to nearly the speed of light. Their goal was to produce a clean, analytical description of the waves in the time domain, which is the way a signal appears as it changes over time, rather than as a collection of frequencies.

The team achieved this by developing a new mathematical technique to break down the complex, three-dimensional spray of gravitational radiation into simpler, distinct components. Imagine trying to understand the sound of a complex chord by isolating each individual note; the researchers did something similar for gravitational waves, separating the signal into specific patterns based on their shape and orientation. By treating the calculation as a problem of cutting and integrating loops in a specific way, they were able to use powerful computational tools to solve equations that had previously been too difficult to handle. This approach allowed them to bypass the complicated, oscillating functions that usually appear in frequency-based calculations and instead find solutions made of simple, familiar mathematical building blocks like fractions and logarithms.

The result is a set of precise formulas that describe the gravitational wave signal for any speed, from slow to near-light speed, in the frame of reference of one of the black holes. These formulas depend on only two variables: the relative speed of the black holes and the time elapsed since the closest approach. The researchers validated their findings by comparing them against existing numerical simulations and independent calculations, finding perfect agreement in every case. They provided the first closed-form expressions for the dominant wave patterns, known as quadrupolar and octupolar modes, which are the primary shapes of the radiation emitted during such an encounter.

This work is significant because it provides a missing piece of the puzzle for understanding binary black hole systems. While most gravitational wave detections come from black holes that spiral together and merge, the physics of how they scatter at high speeds contains the same fundamental information needed to model the merger process. By capturing the full speed dependence without approximation, these new formulas offer a rigorous benchmark for future models. They allow scientists to test their theories against a known, exact solution rather than relying on truncated approximations that might miss subtle relativistic effects. The team has also made the detailed coefficients for these wave patterns available for other researchers to use, paving the way for more accurate models of black hole collisions and a deeper understanding of how gravity behaves in the most extreme environments.

Technical Summary: Analytical Time-Domain Waveforms for Black-Hole Scattering in the Post-Minkowskian Expansion

Problem Statement
Gravitational-wave observations require waveform models that accurately describe the motion of binary bodies and the radiation they emit. While Effective-One-Body (EOB) models and other waveform families successfully combine numerical relativity with analytic approximations, they traditionally rely on the post-Newtonian (PN) expansion, which assumes small velocities. In contrast, the post-Minkowskian (PM) expansion assumes weak fields but allows for arbitrary velocities, making it suitable for high-energy black-hole scattering. Although PM results for scattering observables (such as the scattering angle and impulse) have been computed to high precision, the gravitational waveform itself has not yet been directly incorporated into EOB models. A primary obstacle is the lack of analytic time-domain multipole moments derived from PM scattering amplitudes. Existing frequency-domain PM waveforms involve Bessel-like functions and require complex angular projections and Fourier transforms to obtain time-domain multipoles, a process that often necessitates subtle frame choices and lacks closed-form analytic expressions valid to all orders in velocity.

Methodology
The authors present a method to compute the time-domain multipoles of the gravitational waveform from hyperbolic black-hole scattering at leading PM order, retaining full dependence on velocity. The core of their approach involves a covariant partial-wave projection performed as a cut loop integral.

Covariant Projection: The waveform is decomposed into spin-weighted spherical harmonics. Instead of standard angular integration, the authors express these harmonics as symmetric-trace-free (STF) tensors constructed from covariant kinematic vectors. This allows the projection onto linearly polarized modes ( $f^+_{\ell m}$ and $f^\times_{\ell m}$ ) to be formulated as an integral over on-shell graviton momentum.
Cut Loop Integral Formulation: By combining the partial-wave projection with the Fourier transform to impact-parameter space and the integration over graviton frequency, the problem is recast as a cut Feynman integral with non-integer propagator powers. This formulation utilizes generalized unitarity within Worldline Quantum Field Theory (WQFT).
Integration Techniques: The resulting integrals are reduced to a basis of master integrals using Integration-by-Parts (IBP) identities and the reverse unitarity method. The authors derive canonical differential equations for these master integrals with respect to two dimensionless variables, $x$ (related to the relative Lorentz factor $\gamma$ ) and $t$ (related to retarded time $u$ ). These equations are solved via path-ordered iterated integration, with boundary conditions determined in the low-velocity limit.

Key Contributions and Results
The paper provides the first closed analytic expressions for the time-domain quadrupolar ( $\ell=2$ ) and octupolar ( $\ell=3$ ) modes of the leading-order PM waveform in the rest frame of one of the black holes.

Functional Form: Unlike frequency-domain results that involve Bessel-like functions, the derived time-domain multipoles are composed entirely of elementary functions: rational functions, logarithms, and arctangents. The general form of the multipole moments is given by:
$f^p_{\ell m}(x, t) \propto \frac{1}{(\gamma^2 - 1)^{\ell/2}} \left[ {}_0R^p_{\ell m} + {}_1R^p_{\ell m} \log x + {}_2R^p_{\ell m} \arctan(\dots) + {}_3R^p_{\ell m} \log(\dots) + {}_4R^p_{\ell m} \log(\dots) \right]$
where the coefficients ${}_iR^p_{\ell m}$ are rational functions of the variables $x$ and $t$ .
Velocity Dependence: The expressions retain the full dependence on velocity (encoded in $x$ ) and retarded time. Consequently, they encode an infinite sequence of PN corrections at a fixed PM order, providing analytic benchmarks beyond any finite velocity truncation.
Validation: The results are validated against three independent methods:
1. Direct numerical projection of the leading-order PM waveform.
2. The PN expansion of the multipolar post-Minkowskian (MPM) waveform.
3. An independent self-force calculation of the $\ell=2$ multipole up to 20PN order.
  The authors report complete agreement across all validation checks.
Symmetry Properties: The paper details the symmetry of the modes under reflection ( $y \to -y$ ), establishing relations between modes with opposite $m$ values and between the two polarizations.

Significance
The authors claim this work establishes explicit relativistic counterparts to the PN modes currently used in waveform models. By providing analytic time-domain multipoles valid to all orders in velocity, the results pave the way for incorporating post-Minkowskian fluxes directly into waveform models for binary mergers. The methodology offers a systematic route to computing higher multipoles and higher PM orders, potentially resolving the mismatch between amplitude-based PM calculations and traditional PN waveform modeling. The work also facilitates the investigation of analytic continuation from scattering to bound waveforms and provides a common basis for comparing PM, PN, numerical relativity, and self-force calculations.

Technical Summary: Analytical Time-Domain Waveforms for Black-Hole Scattering in the Post-Minkowskian Expansion

More like this