# NRHJSur3dq8: a spectral numerical-relativity surrogate family for non-eccentric aligned-spin binary-black-hole waveforms

Preprint · Vaishak Prasad

- Archive id: arxiv:2609.09088
- Published: 2026-09-08
- arXiv: https://arxiv.org/abs/2609.09088
- Source URL: https://arxiv.org/abs/2609.09088
- PDF: https://arxiv.org/pdf/2609.09088
- Provenance: arxiv

## Abstract

We present NRHJSur3dq8, a family of numerical-relativity surrogates for non-eccentric aligned-spin binary-black-hole waveforms inspired by an action-angle parameterization. The model provides both the dynamics: the Hamiltonian values, azimuthal action $J\_φ(Φ)$, the orbital clock $τ(Φ)$ as functions of its phasing, and the radiative pieces: the waveforms, their near-analytic derivatives, and modeling errors associated with each. The surrogated variables are conditioned into pieces that vary slowly on the radiation-reaction timescale, represented on shared $h$-adaptive Chebyshev elements, and regressed across the intrinsic parameters with a Gaussian-process model. A time-parameterized merger-ringdown set of elements is attached where the adiabatic description breaks down. Fixed-hyperparameter leave-one-out validation of the central NRHJSur3dq8\_AA model variant against numerical relativity gives a $(2,2)$ unweighted mismatch of $1.6^{+62}\_{-1.5}\\times10^{-6}$ and an all-mode mismatch of $2.8^{+84}\_{-2.7}\\times10^{-6}$ over the full inspiral-merger-ringdown span. Sky-averaged and signal-to-noise-ratio-weighted over $32$ orientations per simulation, the model's full inspiral-merger-ringdown in-sample mismatch against the SXS numerical-relativity simulations used for training is $8.9\\times10^{-7}$ to $2.8\\times10^{-5}$ across total masses $40$-$120 M\_\\odot$. Waveforms are hybridized using post-Newtonian waveforms. Beyond accuracy, NRHJSur3dq8 provides two capabilities uncommon in numerical-relativity surrogates: nearly analytic derivatives of the strain with respect to the physical parameters, and a calibrated Gaussian-process predictive uncertainty of the model. A full non-hybridized inspiral-merger-ringdown waveform evaluates in $12.5$ ms on a single CPU call. Batched evaluation takes $0.43$ ms per waveform at a batch of $4096$, measured at a starting frequency of $25$ Hz.

_Abstract reproduced under CC0-1.0 (arXiv descriptive metadata)._

---
Space Kitty Archive · https://space-kitty.com · API: https://space-kitty.com/api/v1 · MCP: https://space-kitty.com/mcp · Premium: $0.20/request via x402 (Base Sepolia (testnet)). Metadata: arXiv (CC0), OpenAlex (CC0), Europe PMC (CC BY/CC0 works), ClinicalTrials.gov.