A Parameter Estimation Method that Directly Compares Gravitational Wave Observations to Numerical Relativity

Abstract

We present and assess a Bayesian method to interpret gravitational wave signals from binary black holes. Our method directly compares gravitational wave data to numerical relativity simulations. This procedure bypasses approximations used in semi-analytical models for compact binary coalescence. In this work, we use only the full posterior parameter distribution for generic nonprecessing binaries, drawing inferences away from the set of NR simulations used, via interpolation of a single scalar quantity (the marginalized log-likelihood, $\ln \mathcal{L}$) evaluated by comparing data to nonprecessing binary black hole simulations. We also compare the data to generic simulations, and discuss the effectiveness of this procedure for generic sources. We specifically assess the impact of higher order modes, repeating our interpretation with both $l\le 2$ as well as $l\le 3$ harmonic modes. Using the $l\le 3$ higher modes, we gain more information from the signal and can better constrain the parameters of the gravitational wave signal. We assess and quantify several sources of systematic error that our procedure could introduce, including simulation resolution and duration; most are negligible. We show through examples that our method can recover the parameters for equal mass, zero spin; GW150914-like; and unequal mass, precessing spin sources. Our study of this new parameter estimation method demonstrates we can quantify and understand the systematic and statistical error. This method allows us to use higher order modes from numerical relativity simulations to better constrain the black hole binary parameters.

30 pages, 22 figures; submitted to PRD

References in corpus (14)

• Inspiral, merger and ringdown of unequal mass black hole binaries: a multipolar analysis
• Model Waveform Accuracy Standards for Gravitational Wave Data Analysis
• Analysis of spin precession in binary black hole systems including quadrupole-monopole interaction
• LISA detections of massive black hole inspirals: parameter extraction errors due to inaccurate template waveforms
• Where post-Newtonian and numerical-relativity waveforms meet
• Anatomy of the binary black hole recoil: A multipolar analysis
• Extraction of Gravitational Waves in Numerical Relativity
• Gravitational-wave observations of binary black holes: Effect of non-quadrupole modes
• Modeling the source of GW150914 with targeted numerical-relativity simulations
• Gravitational waves from BH-NS binaries: Effective Fisher matrices and parameter estimation using higher harmonics
• Perturbative extraction of gravitational waveforms generated with Numerical Relativity
• Comparison of subdominant gravitational wave harmonics between post-Newtonian and numerical relativity calculations and construction of multi-mode hybrids
• Intrinsic selection biases of ground-based gravitational wave searches for high-mass BH-BH mergers
• Numerical Relativity meets Data Analysis: Spinning Binary Black Hole Case

Cited by in corpus (2)

• Statistical Gravitational Waveform Models: What to Simulate Next?
• Accurate inspiral-merger-ringdown gravitational waveforms for non-spinning black-hole binaries including the effect of subdominant modes