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 (24)

• Scikit-learn: Machine Learning in Python
• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW151226: Observation of Gravitational Waves from a 22-Solar-Mass Binary Black Hole Coalescence
• GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2
• Frequency-domain gravitational waves from non-precessing black-hole binaries. II. A phenomenological model for the advanced detector era
• Properties of the Binary Black Hole Merger GW150914
• Astrophysical Implications of the Binary Black-Hole Merger GW150914
• Frequency-domain gravitational waves from non-precessing black-hole binaries. I. New numerical waveforms and anatomy of the signal
• A simple model of complete precessing black-hole-binary gravitational waveforms
• An improved effective-one-body model of spinning, nonprecessing binary black holes for the era of gravitational-wave astrophysics with advanced detectors
• Effective-one-body model for black-hole binaries with generic mass ratios and spins
• Inspiral-merger-ringdown waveforms of spinning, precessing black-hole binaries in the effective-one-body formalism
• Validating the effective-one-body model of spinning, precessing binary black holes against numerical relativity
• Fast prediction and evaluation of gravitational waveforms using surrogate models
• Fast and accurate prediction of numerical relativity waveforms from binary black hole coalescences using surrogate models
• Frequency domain reduced order models for gravitational waves from aligned-spin compact binaries
• Frequency domain reduced order model of aligned-spin effective-one-body waveforms with generic mass-ratios and spins
• Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence
• Effects of waveform model systematics on the interpretation of GW150914
• Modeling the source of GW150914 with targeted numerical-relativity simulations
• A Parameter Estimation Method that Directly Compares Gravitational Wave Observations to Numerical Relativity
• Improving gravitational-wave parameter estimation using Gaussian process regression
• Novel Method for Incorporating Model Uncertainties into Gravitational Wave Parameter Estimates
• An architecture for efficient gravitational wave parameter estimation with multimodal linear surrogate models

Cited by in corpus (45)

• Surrogate model of hybridized numerical relativity binary black hole waveforms
• Multi-messenger Bayesian parameter inference of a binary neutron-star merger
• Constraints on the neutron star equation of state from AT2017gfo using radiative transfer simulations
• Non-parametric inference of the neutron star equation of state from gravitational wave observations
• Including higher order multipoles in gravitational-wave models for precessing binary black holes
• Enhancing Gravitational-Wave Science with Machine Learning
• Laying the foundation of the effective-one-body waveform models SEOBNRv5: improved accuracy and efficiency for spinning non-precessing binary black holes
• Eccentric, nonspinning, inspiral, Gaussian-process merger approximant for the detection and characterization of eccentric binary black hole mergers
• High-accuracy mass, spin, and recoil predictions of generic black-hole merger remnants
• Frequency domain reduced order model of aligned-spin effective-one-body waveforms with higher-order modes
• GRANDMA Observations of Advanced LIGO's and Advanced Virgo's Third Observational Campaign
• The Science of the Einstein Telescope
• Surrogate model for an aligned-spin effective one body waveform model of binary neutron star inspirals using Gaussian process regression
• Implications of the search for optical counterparts during the first six months of the Advanced LIGO's and Advanced Virgo's third observing run: possible limits on the ejecta mass and binary properties
• Cover Your Basis: Comprehensive Data-Driven Characterization of the Binary Black Hole Population
• A Surrogate Model for Gravitational Wave Signals from Comparable- to Large- Mass-Ratio Black Hole Binaries
• The first six months of the Advanced LIGO's and Advanced Virgo's third observing run with GRANDMA
• A hybrid post-Newtonian -- effective-one-body scheme for spin-precessing compact-binary waveforms
• Ain't No Mountain High Enough: Semi-Parametric Modeling of LIGO-Virgos Binary Black Hole Mass Distribution
• On standardizing kilonovae and their use as standard candles to measure the Hubble constant
• Waveform accuracy and systematic uncertainties in current gravitational wave observations
• Possible Causes of False General Relativity Violations in Gravitational Wave Observations
• Implications of the search for optical counterparts during the second part of the Advanced LIGO's and Advanced Virgo's third observing run: lessons learned for future follow-up observations
• A Precessing Numerical Relativity Waveform Surrogate Model for Binary Black Holes: A Gaussian Process Regression Approach
• Physics of eccentric binary black hole mergers: A numerical relativity perspective
• Regression methods in waveform modeling: a comparative study
• Fast post-adiabatic waveforms in the time domain: Applications to compact binary coalescences in LIGO and Virgo
• Reduced Order and Surrogate Models for Gravitational Waves
• A flexible Gaussian process reconstruction method and the mass function of the coalescing binary black hole systems
• Gaussian Processes for Glitch-robust Gravitational-wave Astronomy
• Rapid and accurate parameter inference for coalescing, precessing compact binaries
• Applications of machine learning in gravitational wave research with current interferometric detectors
• KilonovaNet: Surrogate Models of Kilonova Spectra with Conditional Variational Autoencoders
• Probabilistic Model for the Gravitational Wave Signal from Merging Black Holes
• Incorporation of model accuracy in gravitational wave Bayesian inference
• Accounting for numerical-relativity calibration uncertainty in gravitational-wave modeling and inference
• A fully precessing higher-mode surrogate model of effective-one-body waveforms
• Active Learning for Computationally Efficient Distribution of Binary Evolution Simulations
• Enhancing gravitational waveform models through dynamic calibration
• Optimizing the Placement of Numerical Relativity Simulations using a Mismatch Predicting Neural Network
• Predicting the motion of a high-Q pendulum subject to seismic perturbations using machine learning
• Machine Learning for Quantum-Enhanced Gravitational-Wave Observatories
• Kilonova-Targeting Lightcurve Classification for Wide Field Survey Telescope
• Standing Wave Decomposition Gaussian Process
• Conditional Autoencoder for Generating Binary Neutron Star Waveforms with Tidal and Precession Effects