Abstract:The current and upcoming generations of gravitational wave experiments represent an exciting step forward in terms of detector sensitivity and performance. For example, key upgrades at the LIGO, Virgo and KAGRA facilities will see the next observing run (O4) probe a spatial volume around four times larger than the previous run (O3), and design implementations for e.g. the Einstein Telescope, Cosmic Explorer and LISA experiments are taking shape to explore a wider frequency range and probe cosmic distances. In this context, however, a number of very real data analysis problems face the gravitational wave community. For example, it will be crucial to develop tools and strategies to analyse (amongst other scenarios) signals that arrive coincidentally in detectors, longer signals that are in the presence of non-stationary noise or other shorter transients, as well as noisy, potentially correlated, coherent stochastic backgrounds. With these challenges in mind, we develop peregrine, a new sequential simulation-based inference approach designed to study broad classes of gravitational wave signal. In this work, we describe the method and implementation, before demonstrating its accuracy and robustness through direct comparison with established likelihood-based methods. Specifically, we show that we are able to fully reconstruct the posterior distributions for every parameter of a spinning, precessing compact binary coalescence using one of the most physically detailed and computationally expensive waveform approximants (SEOBNRv4PHM). Crucially, we are able to do this using only 2% of the waveform evaluations that are required in e.g. nested sampling approaches. Finally, we provide some outlook as to how this level of simulation efficiency and flexibility in the statistical analysis could allow peregrine to tackle these current and future gravitational wave data analysis problems.

EMRI_MC: A GPU-based Python code for Bayesian inference of EMRI waveforms

Fast and Fourier: Extreme Mass Ratio Inspiral Waveforms in the Frequency Domain

Rapid Parameter Estimation for Extreme Mass Ratio Inspirals Using Machine Learning

High-speed reconstruction of long-duration gravitational waves from extreme mass ratio inspirals using sparse dictionary learning

Parallelized Inference for Gravitational-Wave Astronomy

BHPWAVE: An adiabatic gravitational waveform model for compact objects undergoing quasi-circular inspirals into rotating massive black holes

Massively parallel Bayesian inference for transient gravitational-wave astronomy

Probing fundamental physics with Extreme Mass Ratio Inspirals: a full Bayesian inference for scalar charge

Simulation-based inference for stochastic gravitational wave background data analysis

Impact of relativistic waveforms in LISA's science objectives with extreme-mass-ratio inspirals

Rapid determination of LISA sensitivity to extreme mass ratio inspirals with machine learning

GERry: A Code to Optimise the Hunt for the Electromagnetic Counter-parts to Gravitational Wave Events

Peregrine: Sequential simulation-based inference for gravitational wave signals

Bilby: A User-friendly Bayesian Inference Library for Gravitational-wave Astronomy

Adiabatic gravitational waveform model for compact objects undergoing quasicircular inspirals into rotating massive black holes

Bayesian Method for Waveform Analysis with GPU Acceleration

Fast and faithful interpolation of numerical relativity surrogate waveforms using meshfree approximation

A test for LISA foreground Gaussianity and stationarity. II. Extreme mass-ratio inspirals

PyROQ: a Python-based Reduced Order Quadrature Building Code for Fast Gravitational Wave Inference

Tidally-induced nonlinear resonances in EMRIs with an analogue model

Eryn : A multi-purpose sampler for Bayesian inference