Fast and faithful interpolation of numerical relativity surrogate waveforms using meshfree approximation

Abstract

Several theoretical waveform models have been developed over the years to capture the gravitational wave emission from the dynamical evolution of compact binary systems of neutron stars and black holes. As ground-based detectors improve their sensitivity at low frequencies, the real-time computation of these waveforms can become computationally expensive, exacerbating the steep cost of rapidly reconstructing source parameters using Bayesian methods. This paper describes an efficient numerical algorithm for generating high-fidelity interpolated compact binary waveforms at an arbitrary point in the signal manifold by leveraging computational linear algebra techniques such as singular value decomposition and meshfree approximation. The results are presented for the time-domain \texttt{NRHybSur3dq8} inspiral-merger-ringdown (IMR) waveform model that is fine tuned to numerical relativity simulations and parameterized by the two component-masses and two aligned spins. For demonstration, we target a specific region of the intrinsic parameter space inspired by the previously inferred parameters of the \texttt{GW200311\_115853} event -- a binary black hole system whose merger was recorded by the network of advanced-LIGO and Virgo detectors during the third observation run. We show that the meshfree interpolated waveforms can be evaluated in $\sim 2.3$ ms, which is about $\times 38$ faster than its brute-force (frequency-domain tapered) implementation in the \textsc{PyCBC} software package at a median accuracy of $\sim \mathcal{O}(10^{-5})$. The algorithm is computationally efficient and scales favourably with an increasing number of dimensions of the parameter space. This technique may find use in rapid parameter estimation and source reconstruction studies.

15 pages, 9 figures

References in corpus (5)

• An improved effective-one-body model of spinning, nonprecessing binary black holes for the era of gravitational-wave astrophysics with advanced detectors
• The third and a half post-Newtonian gravitational wave quadrupole mode for quasi-circular inspiralling compact binaries
• Accelerated gravitational-wave parameter estimation with reduced order modeling
• Fast likelihood evaluation using meshfree approximations for reconstructing compact binary sources
• Efficient Reduced Order Quadrature Construction Algorithms for Fast Gravitational Wave Inference

Cited by in corpus (2)

• A Deep Learning Powered Numerical Relativity Surrogate for Binary Black Hole Waveforms
• Auto-encoder model for faster generation of effective one-body gravitational waveform approximations