Validating Sequential Monte Carlo for Gravitational-Wave Inference
arXiv:2506.18977 · doi:10.1093/mnras/staf1458
Abstract
Nested sampling (NS) is the preferred stochastic sampling algorithm for gravitational-wave inference for compact binary coalenscences (CBCs). It can handle the complex nature of the gravitational-wave likelihood surface and provides an estimate of the Bayesian model evidence. However, there is another class of algorithms that meets the same requirements but has not been used for gravitational-wave analyses: Sequential Monte Carlo (SMC), an extension of importance sampling that maps samples from an initial density to a target density via a series of intermediate densities. In this work, we validate a type of SMC algorithm, called persistent sampling (PS), for gravitational-wave inference. We consider a range of different scenarios including binary black holes (BBHs) and binary neutron stars (BNSs) and real and simulated data and show that PS produces results that are consistent with NS whilst being, on average, 2 times more efficient and 2.74 times faster. This demonstrates that PS is a viable alternative to NS that should be considered for future gravitational-wave analyses.
arXiv comment14 pages, 17 figures, 6 tables, documentation for the accompanying data release is available at http://gw-smc.michaeljwilliams.me/
References in corpus (61)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- Advanced Virgo: a 2nd generation interferometric gravitational wave detector
- GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs
- Advanced LIGO
- MCMC using Hamiltonian dynamics
- GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run
- Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
- Gravitational Waves from Mergin Compact Binaries: How Accurately Can One Extract the Binary's Parameters from the Inspiral Waveform?
- Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy
- Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
- A simple model of complete precessing black-hole-binary gravitational waveforms
- Computationally efficient models for the dominant and sub-dominant harmonic modes of precessing binary black holes
- Open data from the first and second observing runs of Advanced LIGO and Advanced Virgo
- PolyChord: next-generation nested sampling
- GW150914: First results from the search for binary black hole coalescence with Advanced LIGO
- Bayesian inference for compact binary coalescences with BILBY: Validation and application to the first LIGO--Virgo gravitational-wave transient catalogue
- An introduction to Bayesian inference in gravitational-wave astronomy: parameter estimation, model selection, and hierarchical models
- BayesWave: Bayesian Inference for Gravitational Wave Bursts and Instrument Glitches
- Surrogate models for precessing binary black hole simulations with unequal masses
- On Particle Methods for Parameter Estimation in State-Space Models
- Open data from the third observing run of LIGO, Virgo, KAGRA and GEO
- PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals
- Bayesian coherent analysis of in-spiral gravitational wave signals with a detector network
- Real-time gravitational-wave science with neural posterior estimation
- Improving the NRTidal model for binary neutron star systems
- 4-OGC: Catalog of gravitational waves from compact-binary mergers
- Massively parallel Bayesian inference for transient gravitational-wave astronomy
- Observation of Gravitational Waves from the Coalescence of a Compact Object and a Neutron Star
- Fast and Accurate Inference on Gravitational Waves from Precessing Compact Binaries
- SEOBNRv5PHM: Next generation of accurate and efficient multipolar precessing-spin effective-one-body waveforms for binary black holes
- Nested Sampling with Normalising Flows for Gravitational-Wave Inference
- Calibration Uncertainty for Advanced LIGO's First and Second Observing Runs
- Efficient Bayesian inference for multimodal problems in cosmology
- Observing and measuring the neutron-star equation-of-state in spinning binary neutron star systems
- General-relativistic precession in a black-hole binary
- Interpreting Binary Neutron Star Mergers: Describing the Binary Neutron Star Dynamics, Modelling Gravitational Waveforms, and Analyzing Detections
- PESummary: the code agnostic Parameter Estimation Summary page builder
- Nested Sampling Methods
- Parameter estimation with gravitational waves
- Reconstructing the calibrated strain signal in the Advanced LIGO detectors
- Bilby-MCMC: An MCMC sampler for gravitational-wave inference
- SWIGLAL: Python and Octave interfaces to the LALSuite gravitational-wave data analysis libraries
- Accelerating astronomical and cosmological inference with Preconditioned Monte Carlo
- Rapid Parameter Estimation of Gravitational Waves from Binary Neutron Star Coalescence using Focused Reduced Order Quadrature
- Eryn : A multi-purpose sampler for Bayesian inference
- Sampling Errors in Nested Sampling Parameter Estimation
- Gravitational-wave astronomy with a physical calibration model
- Importance nested sampling with normalising flows
- Removing degeneracy and multimodality in gravitational wave source parameters
- : Bayesian inference of multimessenger astrophysical data, methods and application to gravitational-waves
- Revisiting the evidence for precession in GW200129 with machine learning noise mitigation
- VARAHA: A Fast Non-Markovian sampler for estimating Gravitational-Wave posteriors
- Fast marginalization algorithm for optimizing gravitational wave detection, parameter estimation and sky localization
- Adaptive Kernel Density Estimation proposal in gravitational wave data analysis
- Validating Sequential Monte Carlo for Gravitational-Wave Inference
- Persistent Sampling: Enhancing the Efficiency of Sequential Monte Carlo
- MCMC for multi-modal distributions
- Sequential Kalman Tuning of the -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems