Statistical Gravitational Waveform Models: What to Simulate Next?

Abstract

Models of gravitational waveforms play a critical role in detecting and characterizing the gravitational waves (GWs) from compact binary coalescences. Waveforms from numerical relativity (NR), while highly accurate, are too computationally expensive to produce to be directly used with Bayesian parameter estimation tools like Markov-chain-Monte-Carlo and nested sampling. We propose a Gaussian process regression (GPR) method to generate accurate reduced-order-model waveforms based only on existing accurate (e.g. NR) simulations. Using a training set of simulated waveforms, our GPR approach produces interpolated waveforms along with uncertainties across the parameter space. As a proof of concept, we use a training set of IMRPhenomD waveforms to build a GPR model in the 2-d parameter space of mass ratio $q$ and equal-and-aligned spin ${\chi }_1={\chi }_2$. Using a regular, equally-spaced grid of 120 IMRPhenomD training waveforms in $q\in [1,3]$ and ${\chi }_1\in [-0.5,0.5]$, the GPR mean approximates IMRPhenomD in this space to mismatches below $4.3\times 10^{-5}$. Our approach can alternatively use training waveforms directly from numerical relativity. Beyond interpolation of waveforms, we also present a greedy algorithm that utilizes the errors provided by our GPR model to optimize the placement of future simulations. In a fiducial test case we find that using the greedy algorithm to iteratively add simulations achieves GPR errors that are $\sim 1$ order of magnitude lower than the errors from using Latin-hypercube or square training grids.

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
• Validating the effective-one-body model of spinning, precessing binary black holes against numerical relativity
• Modeling the source of GW150914 with targeted numerical-relativity simulations
• A Parameter Estimation Method that Directly Compares Gravitational Wave Observations to Numerical Relativity
• Novel Method for Incorporating Model Uncertainties into Gravitational Wave Parameter Estimates

Cited by in corpus (1)

• Eccentric, nonspinning, inspiral, Gaussian-process merger approximant for the detection and characterization of eccentric binary black hole mergers