Parametrized spin-precessing inspiral-merger-ringdown waveform model for tests of general relativity

The coalescence of binary black holes (BBHs) provides a unique arena to test general relativity (GR) in the dynamical, strong-field regime. To this end, we present , a parametrized, multipolar, spin-precessing waveform model for BBHs in quasicircular orbits, built within the effective-one-body formalism. Compared to its predecessor, , our model introduces parametrized deviations from GR not only in the plunge-merger-ringdown stages, but also in the inspiral phase through modifications to the conservative dynamics. Additionally, it incorporates, for the first time, spin-precession effects. The free deviation parameters can be used to perform null tests of GR using current and future gravitational-wave observations. We validate through Bayesian parameter estimation, focusing on the quasinormal-mode frequency and damping time of the <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mo stretchy="false">(</a:mo><a:mo>ℓ</a:mo><a:mo>,</a:mo><a:mi>m</a:mi><a:mo>,</a:mo><a:mi>n</a:mi><a:mo stretchy="false">)</a:mo><a:mo>=</a:mo><a:mo stretchy="false">(</a:mo><a:mn>2</a:mn><a:mo>,</a:mo><a:mn>2</a:mn><a:mo>,</a:mo><a:mn>0</a:mn><a:mo stretchy="false">)</a:mo></a:math> mode. Our analysis of synthetic signals from numerical-relativity (NR) simulations of highly precessing BH mergers shows that, while correctly recovers consistency with GR, neglecting spin precession can lead to false detections of deviations from GR even at current detector sensitivity. Conversely, when analyzing a synthetic signal from a NR simulation of a binary boson-star merger, the model successfully identifies a deviation from a GR BBH signal. Finally, we reanalyze 12 events from the third Gravitational Wave Transient Catalog. Using a hierarchical combination of these events, we constrain fractional deviations in the frequency and damping time of the (2,2,0) quasinormal-mode to <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mi>δ</g:mi><g:msub><g:mi>f</g:mi><g:mn>220</g:mn></g:msub><g:mo>=</g:mo><g:mn>0.0</g:mn><g:msubsup><g:mn>0</g:mn><g:mrow><g:mo>−</g:mo><g:mn>0.06</g:mn></g:mrow><g:mrow><g:mo>+</g:mo><g:mn>0.06</g:mn></g:mrow></g:msubsup></g:math> and <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:mi>δ</i:mi><i:msub><i:mi>τ</i:mi><i:mn>220</i:mn></i:msub><i:mo>=</i:mo><i:mn>0.1</i:mn><i:msubsup><i:mn>5</i:mn><i:mrow><i:mo>−</i:mo><i:mn>0.24</i:mn></i:mrow><i:mrow><i:mo>+</i:mo><i:mn>0.26</i:mn></i:mrow></i:msubsup></i:math> at 90% credibility. These results are consistent with those from the LIGO-Virgo-KAGRA Collaboration, which did not account for spin-precession effects.