Novel Method for Incorporating Model Uncertainties into Gravitational Wave Parameter Estimates

Abstract

Posterior distributions on parameters computed from experimental data using Bayesian techniques are only as accurate as the models used to construct them. In many applications these models are incomplete, which both reduces the prospects of detection and leads to a systematic error in the parameter estimates. In the analysis of data from gravitational wave detectors, for example, accurate waveform templates can be computed using numerical methods, but the prohibitive cost of these simulations means this can only be done for a small handful of parameters. In this work a novel method to fold model uncertainties into data analysis is proposed; the waveform uncertainty is analytically marginalised over using with a prior distribution constructed by using Gaussian process regression to interpolate the waveform difference from a small training set of accurate templates. The method is well motivated, easy to implement, and no more computationally expensive than standard techniques. The new method is shown to perform extremely well when applied to a toy problem. While we use the application to gravitational wave data analysis to motivate and illustrate the technique, it can be applied in any context where model uncertainties exist.

6 pages, 3 figures, accepted for publication in Physical Review Letters

References in corpus (2)

• LISA detections of massive black hole inspirals: parameter extraction errors due to inaccurate template waveforms
• New advances in the Gaussian-process approach to pulsar-timing data analysis

Cited by in corpus (4)

• Eccentric, nonspinning, inspiral, Gaussian-process merger approximant for the detection and characterization of eccentric binary black hole mergers
• Surrogate model for an aligned-spin effective one body waveform model of binary neutron star inspirals using Gaussian process regression
• Statistical Gravitational Waveform Models: What to Simulate Next?
• 2 Fast 2 Fiducial: Gaussian processes for the interpolation and marginalization of waveform error in extreme-mass-ratio-inspiral parameter estimation