Gravitational-wave Parameter Estimation in Non-Gaussian Noise Using Score-based Likelihood Characterization

Gravitational-wave (GW) parameter estimation typically assumes that instrumental noise is Gaussian and stationary. Obvious departures from this idealization are typically handled on a case-by-case basis, e.g., through bespoke procedures to “clean” non-Gaussian noise transients (glitches), as was famously the case for the GW170817 neutron-star binary. Although effective, this data manipulation can bias key astrophysical inferences, such as binary precession, and compound unpredictably when combining multiple observations. Alternative bias-free methods, like joint noise-signal inference, remain too computationally expensive for large-scale execution. Here we take a different approach: rather than explicitly modeling individual non-Gaussianities to then apply the traditional GW likelihood, we seek to learn the true distribution of instrumental noise without presuming Gaussianity and stationarity in the first place. Assuming only noise additivity, we employ score-based diffusion models to learn an empirical noise distribution directly from detector data and then combine it with a deterministic waveform model to provide an unbiased estimate of the likelihood function. We validate the method by performing inference on a subset of GW parameters from 400 mock observations, containing real LIGO noise from either the Livingston or Hanford detectors. We show that the proposed method can recover the true parameters even in the presence of loud glitches, and that the inference is unbiased over a population of signals without applying any cleaning to the data. This work provides a promising avenue for extracting unbiased source properties in future GW observations over the coming decade.