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.

Waveform systematics in gravitational-wave inference of signals from binary neutron star merger models incorporating higher-order modes information

Waveform systematics in gravitational-wave inference of signals from binary neutron star merger models incorporating higher order modes information

Enriching the Symphony of Gravitational Waves from Binary Black Holes by Tuning Higher Harmonics

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

Effects of waveform systematics on inferences of neutron star population properties and the nuclear equation of state

Waveform accuracy and systematic uncertainties in current gravitational wave observations

Insights into Binary Neutron Star Merger Simulations: A Multi-Code Comparison

New and robust gravitational-waveform model for high-mass-ratio binary neutron star systems with dynamical tidal effects

Incorporating waveform calibration error in gravitational-wave modeling and inference for SEOBNRv4

Systematic bias from waveform modeling for binary black hole populations in next-generation gravitational wave detectors

Identifying thermal effects in neutron star merger remnants with model-agnostic waveform reconstructions and third-generation detectors

A phenomenological gravitational waveform model of binary black holes incorporating horizon fluxes

Accuracy limitations of existing numerical relativity waveforms on the data analysis of current and future ground-based detectors

Kilohertz gravitational waves from binary neutron star mergers: Numerical-relativity informed postmerger model

Unraveling information about supranuclear-dense matter from the complete binary neutron star coalescence process using future gravitational-wave detector networks

Numerical-Relativity-Informed Effective-One-Body model for Black-Hole-Neutron-Star Mergers with Higher Modes and Spin Precession

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

Analytically improved and numerical-relativity informed effective-one-body model for coalescing binary neutron stars