Intermediate-mass-ratio black hole binaries: intertwining numerical and perturbative techniques

We describe in detail full numerical and perturbative techniques to compute the gravitational radiation from intermediate-mass-ratio black-hole-binary inspirals and mergers. We perform a series of full numerical simulations of nonspinning black holes with mass ratios q=1/10 and q=1/15 from different initial separations and for different finite-difference resolutions. In order to perform those full numerical runs, we adapt the gauge of the moving punctures approach with a variable damping term for the shift. We also derive an extrapolation (to infinite radius) formula for the waveform extracted at finite radius. For the perturbative evolutions we use the full numerical tracks, transformed into the Schwarzschild gauge, in the source terms of the Regge-Wheller-Zerilli Schwarzschild perturbations formalism. We then extend this perturbative formalism to take into account small intrinsic spins of the large black hole, and validate it by computing the quasinormal mode frequencies, where we find good agreement for spins |a/M|<0.3. Including the final spins improves the overlap functions when comparing full numerical and perturbative waveforms, reaching 99.5% for the leading (l,m)=(2,2) and (3,3) modes, and 98.3% for the nonleading (2,1) mode in the q=1/10 case, which includes 8 orbits before merger. For the q=1/15 case, we obtain overlaps near 99.7% for all three modes. We discuss the modeling of the full inspiral and merger based on a combined matching of post-Newtonian, full numerical, and geodesic trajectories.

💡 Research Summary

The paper presents a comprehensive methodology that combines full numerical relativity (NR) simulations with black‑hole perturbation theory to generate highly accurate gravitational‑wave (GW) signals from intermediate‑mass‑ratio black‑hole binaries (IMRBBs). The authors focus on non‑spinning binaries with mass ratios q = 1/10 and q = 1/15, performing a suite of NR runs at several initial separations and finite‑difference resolutions. To keep the moving‑puncture gauge stable over many orbits, they modify the standard shift condition by introducing a variable damping term. This adaptive gauge reduces coordinate drift and allows the simulations to cover up to eight inspiral cycles before merger even for the more extreme q = 1/15 case.

Because NR waveforms are extracted at finite radii, the authors derive an extrapolation formula that corrects both the 1/r fall‑off and the associated time‑delay, effectively mapping the finite‑radius ψ₄ data to an asymptotic waveform with errors below 10⁻⁴. The extrapolated NR waveforms serve as the benchmark for all subsequent comparisons.

For the perturbative side, the NR trajectories (binary separation and orbital frequency) are transformed from the NR gauge to Schwarzschild coordinates. These trajectories become source terms in the Regge‑Wheeler‑Zerilli (RWZ) equations governing metric perturbations on a Schwarzschild background. The authors extend the RWZ formalism to include first‑order corrections due to a small spin of the larger black hole, parameterised by a/M. They validate this spin‑augmented perturbative model by computing quasinormal‑mode (QNM) frequencies and find agreement within 0.5 % for |a/M| < 0.3.

Waveforms generated from the spin‑corrected RWZ equations are then compared to the NR results using overlap integrals. For the dominant (ℓ,m) = (2,2) and (3,3) modes, the overlaps reach 99.5 % (q = 1/10) and 99.7 % (q = 1/15). The sub‑dominant (2,1) mode also shows excellent agreement, with an overlap of 98.3 % for q = 1/10. Including the final spin of the remnant improves the overlaps by roughly 0.2–0.5 %, demonstrating the importance of even modest spin effects in the perturbative source.

Beyond the technical implementation, the paper outlines a unified modeling pipeline for the full inspiral‑merger‑ringdown signal. At large separations, post‑Newtonian (PN) expansions provide accurate orbital dynamics. In the intermediate regime, NR simulations supply the precise strong‑field evolution. Finally, from the late inspiral through merger and ringdown, the spin‑augmented perturbative approach, combined with geodesic trajectories, yields a seamless continuation of the waveform. This hybrid strategy leverages the strengths of each method while mitigating their individual limitations, offering a practical route to generate IMR waveforms for gravitational‑wave data analysis.

The results have immediate relevance for current detectors (LIGO, Virgo, KAGRA) and future observatories (Einstein Telescope, Cosmic Explorer), where IMR events are expected to be observed with sufficient signal‑to‑noise to demand sub‑percent waveform accuracy. By demonstrating that a modest spin correction and careful gauge handling can push NR‑perturbation overlaps above 99 %, the study provides a robust template‑generation framework that can be extended to more extreme mass ratios, modestly spinning primaries, and eventually to eccentric or precessing systems.