Abstract:Gravitational-wave (GW) observations of binary black-hole (BBH) coalescences are expected to address outstanding questions in astrophysics, cosmology, and fundamental physics. Realizing the full discovery potential of upcoming LIGO-Virgo-KAGRA (LVK) observing runs and new ground-based facilities hinges on accurate waveform models. Using linear-signal approximation methods and Bayesian analysis, we start to assess our readiness for what lies ahead using two state-of-the-art quasi-circular, spin-precessing models: \texttt{SEOBNRv5PHM} and \texttt{IMRPhenomXPHM}. We ascertain that current waveforms can accurately recover the distribution of masses in the LVK astrophysical population, but not spins. We find that systematic biases increase with detector-frame total mass, binary asymmetry, and spin-precession, with most such binaries incurring parameter biases, extending up to redshifts $\sim3$ in future detectors. Furthermore, we examine three ``golden'' events characterized by large mass ratios, significant spin magnitudes, and high precession, evaluating how systematic biases may affect their scientific outcomes. Our findings reveal that current waveforms fail to enable the unbiased measurement of the Hubble-Lemaître parameter from loud signals, even for current detectors. Moreover, highly asymmetric systems within the lower BH mass-gap exhibit biased measurements of the secondary-companion mass, which impacts the physics of both neutron stars and formation channels. Similarly, we deduce that the primary mass of massive binaries ($ > 60 M_\odot$) will also be biased, affecting supernova physics. Future progress in analytical calculations and numerical-relativity simulations, crucial for calibrating the models, must target regions of the parameter space with significant biases to develop more accurate models. Only then can precision GW astronomy fulfill the promise it holds.

Are all models wrong? Falsifying binary formation models in gravitational-wave astronomy

A Trifecta of Modelling Tools: A Bayesian Binary Black Hole Model Selection combining Population Synthesis and Galaxy Formation Models

What You Don't Know Can Hurt You: Use and Abuse of Astrophysical Models in Gravitational-wave Population Analyses

Discriminating between Neutron Stars and Black Holes with Imperfect Knowledge of the Maximum Neutron Star Mass

Is there an excess of black holes around $20 M_{\odot}$? Optimising the complexity of population models with the use of reversible jump MCMC

Incorporating model accuracy into gravitational-wave Bayesian inference

Detecting the Undetected: Overcoming Biases in Gravitational-Wave Population Studies

Systematic Biases in Estimating the Properties of Black Holes Due to Inaccurate Gravitational-Wave Models

Gaussian mixture models of the total mass distribution of stellar black holes from LIGO-Virgo GWTC-2: Implications on the origin of GW190521

One to many: comparing single gravitational-wave events to astrophysical populations

Comparing gravitational waveform models for binary black hole mergers through a hypermodels approach

Poking Holes: Looking for Gaps in LIGO/Virgo's Black Hole Population

Constraining Formation Models of Binary Black Holes with Gravitational-Wave Observations

Probabilistic Model for the Gravitational Wave Signal from Merging Black Holes

Physical Models for the Astrophysical Population of Black Holes: Application to the Bump in the Mass Distribution of Gravitational Wave Sources

Searching for binary black hole sub-populations in gravitational wave data using binned Gaussian processes

Do gravitational wave observations in the lower mass gap favor a hierarchical triple origin?

Distinguish the environmental effects and modified theory of gravity with multiple massive black-hole binaries

Growing Pains: Understanding the Impact of Likelihood Uncertainty on Hierarchical Bayesian Inference for Gravitational-Wave Astronomy

Things that might go bump in the night: Assessing structure in the binary black hole mass spectrum