Results from an injection and recovery study involving 35 numerical relativity simulations and three waveform models

We present Bayesian inference results from an extensive injection-recovery campaign to test the validity of three state of the art quasicircular gravitational waveform models: SEOBNRv5PHM, henom, henom, the latter with the SpinTaylorT4 implementation for its precession dynamics. We analyze 35 strongly precessing binary black hole numerical relativity simulations with all available harmonic content. Ten simulations have a mass ratio of <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mn>4</a:mn><a:mo>:</a:mo><a:mn>1</a:mn></a:math> and five, mass ratio of <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mn>8</c:mn><c:mo>:</c:mo><c:mn>1</c:mn></c:math>. Overall, we find that SEOBNRv5PHM is the most consistent model to numerical relativity, with the majority of true source properties lying within the inferred 90% credible interval. However, we find that none of the models can reliably infer the true source properties for binaries with mass ratio <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mn>8</e:mn><e:mo>:</e:mo><e:mn>1</e:mn></e:math> systems. We additionally conduct inspiral-merger-ringdown (IMR) consistency tests to determine if our chosen state of the art waveform models infer consistent properties when analysing only the inspiral (low frequency) and ringdown (high frequency) portions of the signal. For the simulations considered in this work, we find that the IMR consistency test depends on the frequency that separates the inspiral and ringdown regimes. For two sensible choices of the cutoff frequency, we report that henom can produce false GR deviations. Meanwhile, we find that henom is the most reliable model under the IMR consistency test. Finally, we reanalyze the same 35 simulations, but this time we incorporate model accuracy into our Bayesian inference. Consistent with the work in Hoy Incorporating model accuracy into gravitational-wave Bayesian inference, ., we find that this approach generally yields more accurate inferred properties for binary black holes with less biases compared to methods that combine model-dependent posterior distributions based on their evidence, or with equal weight.