Simulation-based Inference towards Gravitational-wave waveform systematics in Intermediate-Mass Binary Black Holes

TL;DR

This paper introduces a simulation-based inference framework using neural networks to rapidly estimate gravitational-wave source parameters while accounting for waveform model uncertainties, significantly reducing computation time.

Contribution

The authors develop a neural posterior estimation method that marginalizes over waveform model discrepancies, enabling fast, accurate, and systematic-aware gravitational-wave inference.

Findings

Achieves near-equivalent accuracy to traditional methods in less than a second per event.

Successfully marginalizes over waveform model uncertainties using a single trained neural network.

Reduces inference time by several orders of magnitude compared to nested sampling.

Abstract

Parameter estimation for gravitational-wave signals is computationally demanding due to the high dimensionality of the parameter space and the cost of repeated waveform generation in traditional Bayesian inference. These analyses require on the order of 10^8 likelihood evaluations and waveform generations, resulting in inference times of hours to days per event. Furthermore, discrepancies between waveform models introduce systematic uncertainties that can bias inferred source properties. To address these challenges, we propose a novel framework based on Simulation-based Inference (SBI) and Neural Posterior Estimation (NPE) and apply it to signals from Intermediate-Mass Black Holes (IMBH). In this framework, we train a single amortised neural posterior estimator on a large simulated dataset generated using two state-of-the-art waveform approximants, IMRPhenomXPHM and SEOBNRv5PHM. By treating the waveform model index as a latent variable, the network learns to produce posterior distributions that are naturally marginalized over the discrepancies of the two waveform models. Once trained, the model enables direct posterior sampling in milliseconds per event, eliminating the need for likelihood evaluations while simultaneously accounting for model systematics. We demonstrate that this approach recovers accurate posterior distributions for IMBH signals injected into Gaussian noise, achieving close agreement with traditional nested-sampling results while reducing inference time by several orders of magnitude. Our results show that NPE can robustly incorporate waveform-model systematics within a unified framework, offering a scalable path toward rapid, systematics-aware gravitational-wave inference. Establishing these methods as promising alternatives to classical likelihood-based pipelines for current and future high-mass gravitational-wave observations.