Recognition: 2 theorem links
· Lean TheoremBilby: A user-friendly Bayesian inference library for gravitational-wave astronomy
Pith reviewed 2026-05-13 13:49 UTC · model grok-4.3
The pith
Bilby provides a Python library for expert Bayesian parameter estimation of gravitational-wave sources with simple syntax.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Bilby is a Python code that provides expert-level parameter estimation infrastructure with straightforward syntax and tools that facilitate use by beginners. It allows accurate and reliable gravitational-wave parameter estimation on both real, freely-available data from LIGO/Virgo, and simulated data. Examples cover compact binary mergers and other signals including supernovae and binary neutron star merger remnants, while additional functionality supports population studies via hierarchical Bayesian modelling to infer distributions such as the black hole mass function from ensembles of observations.
What carries the argument
Bilby, the Python library that supplies the syntax and tools for setting up Bayesian sampling, signal models, and data analysis for gravitational-wave signals.
If this is right
- Users can switch signal models or implement new likelihood functions without rewriting core sampling infrastructure.
- Analyses of real LIGO/Virgo data and simulated injections can be performed with the same interface.
- Hierarchical modelling allows inference of population-level parameters, such as the shape of the black hole mass distribution, from multiple events.
- New detectors can be added to the analysis framework through modular extensions.
Where Pith is reading between the lines
- New researchers could test custom signal models on public data without first building an entire inference pipeline.
- Faster iteration on analysis methods might become possible once the library handles the underlying sampling and data interfaces.
- Integration with other open tools could expand the set of signals and detectors that non-specialists can study routinely.
Load-bearing premise
The library's code is correctly implemented and well-tested so that the provided examples produce reliable results without hidden errors in the sampling or data processing.
What would settle it
Running one of the compact binary merger examples on real LIGO/Virgo data and obtaining parameter posteriors that match independent published results from other codes to within statistical expectations.
read the original abstract
Bayesian parameter estimation is fast becoming the language of gravitational-wave astronomy. It is the method by which gravitational-wave data is used to infer the sources' astrophysical properties. We introduce a user-friendly Bayesian inference library for gravitational-wave astronomy, Bilby. This python code provides expert-level parameter estimation infrastructure with straightforward syntax and tools that facilitate use by beginners. It allows users to perform accurate and reliable gravitational-wave parameter estimation on both real, freely-available data from LIGO/Virgo, and simulated data. We provide a suite of examples for the analysis of compact binary mergers and other types of signal model including supernovae and the remnants of binary neutron star mergers. These examples illustrate how to change the signal model, how to implement new likelihood functions, and how to add new detectors. Bilby has additional functionality to do population studies using hierarchical Bayesian modelling. We provide an example in which we infer the shape of the black hole mass distribution from an ensemble of observations of binary black hole mergers.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces Bilby, a Python library for Bayesian parameter estimation in gravitational-wave astronomy. It claims to deliver expert-level infrastructure with straightforward syntax suitable for both beginners and experts, enabling accurate analysis of real LIGO/Virgo data and simulations for compact binary mergers, supernovae, binary neutron star remnants, and other signals. Additional functionality supports hierarchical Bayesian population studies, illustrated by an example inferring the black hole mass distribution from binary black hole merger observations. The manuscript provides usage examples demonstrating how to modify signal models, implement new likelihoods, and add detectors.
Significance. If the implementation and validation claims hold, Bilby would be a significant contribution to the field by lowering the barrier to reliable gravitational-wave parameter estimation and population inference. The emphasis on extensibility, open examples, and support for both real and simulated data positions it as a practical tool that could standardize workflows and facilitate broader community participation in LIGO/Virgo analyses.
minor comments (2)
- [Abstract and Examples] The abstract states that results are 'accurate and reliable' on real and simulated data, but the manuscript should include at least one quantitative comparison (e.g., posterior overlap or credible-interval agreement) against an established code such as LALInference for a standard binary black hole injection to support this claim.
- [Population studies section] In the population-inference example, clarify the prior choices and the exact form of the hierarchical likelihood used for the black-hole mass distribution; this would help readers reproduce the reported shape inference.
Simulated Author's Rebuttal
We thank the referee for their positive review of our manuscript and their recommendation to accept. The referee's summary accurately reflects the scope and goals of Bilby as a user-friendly library for Bayesian inference in gravitational-wave astronomy, including its support for both real and simulated data as well as hierarchical population studies.
Circularity Check
No circularity: software library description with no derivations or fitted predictions
full rationale
The paper introduces the Bilby Python library for Bayesian parameter estimation in gravitational-wave astronomy. It describes the API syntax, usage examples for compact binary mergers and other signals, hierarchical population inference, and states that results are reliable on real LIGO/Virgo and simulated data. No mathematical derivations, equations, parameter fittings, or predictions are presented that could reduce to inputs by construction. The work is a software announcement whose soundness depends on the released code and external testing rather than any self-referential chain. No self-citations or ansatzes are load-bearing in a circular manner.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
Cost.FunctionalEquationwashburn_uniqueness_aczel unclearBayesian parameter estimation is fast becoming the language of gravitational-wave astronomy... We introduce a user-friendly Bayesian inference library for gravitational-wave astronomy, Bilby.
-
Foundation.DAlembert.Inevitabilitybilinear_family_forced unclearBilby provides expert-level parameter estimation infrastructure with straightforward syntax... for the analysis of compact binary mergers and other types of signal model
Forward citations
Cited by 24 Pith papers
-
GW240925 and GW250207: Astrophysical Calibration of Gravitational-wave Detectors
The first informative astrophysical calibration of gravitational-wave detectors is reported using GW240925 and GW250207.
-
Testing the Kerr hypothesis beyond the quadrupole with GW241011
GW241011 data shows consistency with Kerr black holes for both quadrupole and octupole moments and delivers the first observational bounds on spin-induced octupole deviations.
-
End-to-End Population Inference from Gravitational-Wave Strain using Transformers
Dingo-Pop uses a transformer to perform amortized, end-to-end population inference from GW strain data in seconds, bypassing per-event Monte Carlo sampling.
-
Lessons from binary dynamics of inspiralling equal-mass boson-star mergers
Numerical simulations of equal-mass boson-star mergers reveal larger waveform deviations from black-hole binaries in late inspiral and merger, plus odd multipole excitations for certain scalar-field phases, with some ...
-
Novel ringdown tests of general relativity with black hole greybody factors
GreyRing model based on greybody factors reproduces numerical relativity ringdown signals with mismatches of order 10^{-6} and enables a new post-merger consistency test of general relativity applied to GW250114.
-
GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
LIGO and Virgo detected 39 compact binary coalescence events in O3a, including 13 new ones, with black hole binaries up to 150 solar masses and the first significantly asymmetric mass ratios.
-
Cracking Gravitational Wave Multiple Ringdown Modes in Space
FIREFLY algorithm enables 200-fold faster multi-mode ringdown analysis for space-borne gravitational wave detectors while remaining compatible with time-delay interferometry.
-
Fast neural network surrogate for multimodal effective-one-body gravitational waveforms from generically precessing compact binaries
Neural network surrogate approximates precessing compact binary gravitational waveforms up to 1000x faster than the base EOB model with validated accuracy.
-
Post-Newtonian inspiral waveform model for eccentric precessing binaries with higher-order modes and matter effects
pyEFPEHM extends prior PN models to include higher-order quasi-circular phasing, generalized precession solutions, and eccentric corrections up to 1PN in selected multipoles for eccentric precessing binaries with matt...
-
GW231123: False Massive Graviton Signatures from Unmodeled Point-Mass Lensing
Unmodeled point-mass lensing produces a spurious nonzero graviton mass posterior in GW231123 that vanishes when lensing is included in the analysis.
-
Computationally efficient models for the dominant and sub-dominant harmonic modes of precessing binary black holes
IMRPhenomXPHM is a new computationally efficient phenomenological model for precessing binary black hole gravitational-wave signals that incorporates higher-order modes via twisting-up maps from non-precessing waveforms.
-
The Impact of Spin Priors on Parameterized Tests of General Relativity
Spin prior choices propagate into tests of GR via the 1.5PN deviation parameter δφ̂3 in a non-trivial, event-dependent way, with stronger effects for short-inspiral events and partial degeneracy with χ_eff when the de...
-
Bayesian Analysis of Gravitational Wave Microlensing Effects from Galactic Double White Dwarfs
Bayesian analysis of simulated Taiji observations shows microlensing from lenses above 10^5 solar masses can be distinguished from unlensed DWD signals when separation is below 3 Einstein radii, while lower masses or ...
-
GW250114: testing Hawking's area law and the Kerr nature of black holes
GW250114 data confirm the remnant black hole ringdown frequencies lie within 30% of Kerr predictions and that the final horizon area is larger than the sum of the progenitors' areas to high credibility.
-
Mitigating Systematic Errors in Parameter Estimation of Binary Black Hole Mergers in O1-O3 LIGO-Virgo Data
Parametric models incorporating waveform phase and amplitude uncertainties mitigate systematic errors in gravitational wave parameter estimation, producing consistent results across models and raw/deglitched data for ...
-
Gravitational-wave astronomy requires population-informed parameter estimation
Population-informed hierarchical parameter estimation is required for unbiased astrophysical interpretation of gravitational-wave events rather than using standard individual posteriors with reference priors.
-
Basilic: An end-to-end pipeline for Bayesian burst inference and model classification in gravitational-wave data
Basilic is an end-to-end Bayesian pipeline for gravitational-wave burst inference and model classification, with a case study showing signal degeneracies between binary black hole mergers and cosmic strings.
-
Investigating the effect of sensitivity of KAGRA on sky localization of gravitational-wave sources from compact binary coalescences
KAGRA enhances sky localization of binary neutron star mergers in the LVK network via added baselines, with measurable gains at current sensitivity and larger improvements as range reaches ~30 Mpc.
-
GW190711_030756 and GW200114_020818: astrophysical interpretation of two asymmetric binary black hole mergers in the IAS catalog
Two asymmetric BBH mergers are characterized with mass ratios 0.35 and ≤0.20; one shows high spins, negative χ_eff, and strong precession, suggesting an emerging population of massive rapidly spinning systems.
-
Improved Constraints on Non-Kerr Deviations from Binary Black Hole Inspirals Using GWTC-4 Data
Bayesian constraints from GWTC-4 binary black hole inspirals show Johannsen metric deformation parameters α13 and ε3 consistent with zero, supporting the Kerr hypothesis.
-
Tests of General Relativity with GWTC-3
No evidence for physics beyond general relativity is found in the analysis of 15 GW events from GWTC-3, with consistency in residuals, PN parameters, and remnant properties.
- Signatures of a subpopulation of hierarchical mergers in the GWTC-4 gravitational-wave dataset
- The Science of the Einstein Telescope
- Inferring the properties of a population of compact binaries in presence of selection effects
Reference graph
Works this paper leans on
-
[1]
The core package The core package provides all of the code required for general problems of inference. It provides a unified inter- face to several different samplers listed below, standard sets of priors including arbitrary user-defined options, and a universal result object that stores all important information from a given simulation. Prior and likelihood...
-
[2]
The gw package The gw package provides the core functionality for parameter estimation specific to transient gravitational waves. Building on the core package, this provides prior specifications unique to such problems, e.g., a prior that is uniform in co-moving volume distance, as well as the standard likelihood used when studying gravitational- wave trans...
work page 2015
-
[3]
B. P. Abbott et al., Phys. Rev. Lett.116, 241102 (2016)
work page 2016
-
[4]
B. P. Abbott et al. , Phys. Rev. X 6, 041014 (2016)
work page 2016
-
[5]
B. P. Abbott et al. , Phys. Rev. X 6, 041015 (2016)
work page 2016
-
[6]
B. P. Abbott et al. , (2018), arXiv:1805.11581
work page Pith review arXiv 2018
- [7]
-
[8]
B. P. Abbott et al., Phys. Rev. Lett.119, 161101 (2017)
work page 2017
-
[9]
B. P. Abbott et al. , Astrophys. J. 848, L13 (2017)
work page 2017
-
[10]
B. P. Abbott et al. , Nature 551, 85 (2017)
work page 2017
- [11]
- [12]
-
[13]
Optimal Search for an As- trophysical Gravitational-Wave Background,
R. Smith and E. Thrane, “Optimal Search for an As- trophysical Gravitational-Wave Background,” (2018)
work page 2018
-
[14]
B. Farr, D. E. Holz, and W. M. Farr, Astrophys. J. 854, L9 (2018)
work page 2018
- [15]
- [16]
-
[17]
B. P. Abbott et al., Phys. Rev. Lett.116, 221101 (2016)
work page 2016
-
[18]
B. P. Abbott et al., Phys. Rev. Lett.119, 141101 (2017)
work page 2017
-
[19]
B. P. Abbott et al., Phys. Rev. Lett.120, 031104 (2018)
work page 2018
-
[20]
B. P. Abbott et al., Phys. Rev. Lett.120, 201102 (2018)
work page 2018
-
[21]
M. Vallisneri, J. Kanner, R. Williams, A. Weinstein, and B. Stephens, in Journal of Physics Conference Se- ries, Journal of Physics Conference Series, Vol. 610 (2015) p. 012021
work page 2015
- [22]
- [23]
-
[24]
gwastro/pycbc: 1.12.3 release,
A. Nitz, I. Harry, D. Brown, C. M. Biwer, et al. , “gwastro/pycbc: 1.12.3 release,” (2018)
work page 2018
-
[25]
B. P. Abbott et al., Phys. Rev. Lett.116, 061102 (2016)
work page 2016
-
[26]
B. P. Abbott et al., Phys. Rev. Lett.116, 241103 (2016)
work page 2016
-
[27]
B. P. Abbott et al., Phys. Rev. Lett.118, 221101 (2017)
work page 2017
-
[28]
B. P. Abbott et al. , Astrophys. J. 851, L35 (2017)
work page 2017
- [29]
-
[30]
E. Thrane and C. Talbot, “An introduction to Bayesian inference in gravitational-wave astronomy: parameter estimation, model selection, and hierarchical models,” (2018), arXiv:1809.02293
-
[31]
Pierce, Types and Programming Languages (The MIT Press, 2002)
B. Pierce, Types and Programming Languages (The MIT Press, 2002)
work page 2002
- [32]
-
[33]
D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Goodman, Publ. Astron. Soc. Pac. 125, 306 (2013)
work page 2013
-
[34]
https://www.python.org/dev/peps/pep-0008/
-
[35]
https://www.python.org/dev/peps/pep-0020/
-
[36]
https://pypi.org/project/BILBY/
-
[37]
https://git.ligo.org/lscsoft/bilby/ ()
-
[38]
https://lscsoft.docs.ligo.org/bilby/ ()
-
[39]
W. D. Vousden, W. M. Farr, and I. Mandel, Mon. Not. R. Astron. Soc. 455, 1919 (2016)
work page 1919
-
[40]
J. Salvatier, T. Wiecki, and C. Fonnesbeck, PeerJ Com- puter Science 2, 55 (2016)
work page 2016
- [41]
- [42]
- [43]
- [44]
-
[45]
http://kylebarbary.com/nestle/
-
[46]
https://github.com/joshspeagle/dynesty
- [47]
-
[48]
Hierarchical Data Format, version 5,
The HDF Group, “Hierarchical Data Format, version 5,” (1997-NNNN), http://www.hdfgroup.org/HDF5/
work page 1997
-
[49]
LSC, “LIGO Scientific Collaboration Algorithm Li- brary software packages: LALInference, LALSimu- lation, and LALInspiral,” https://wiki.ligo.org/ DASWG/LALSuite
-
[50]
D. Macleod, S. Coughlin, A. L. Urban, T. Massinger, et al. , (2018), 10.5281/zenodo.1346349
-
[51]
M. van der Sluys, V. Raymond, I. Mandel, C. R¨ over, N. Christensen, V. Kalogera, R. Meyer, and A. Vecchio, Classical and Quantum Gravity 25, 184011 (2008)
work page 2008
- [52]
- [53]
- [54]
- [55]
- [56]
-
[57]
B. Farr, E. Ochsner, W. M. Farr, and R. O’Shaughnessy, Phys. Rev. D 90, 024018 (2014)
work page 2014
-
[58]
Marginalisation of the time and phase pa- rameters in cbc parameter estimation,
W. M. Farr, “Marginalisation of the time and phase pa- rameters in cbc parameter estimation,” (2014), https: //dcc.ligo.org/T1400460-v2/public
work page 2014
-
[59]
L. P. Singer and L. R. Price, Phys. Rev. D 93, 024013 (2016)
work page 2016
-
[60]
L. P. Singer et al. , Astrophys. J. 829, L15 (2016)
work page 2016
-
[61]
Modelling calibration errors in cbc waveforms,
W. M. Farr, B. Farr, and T. Littenberg, “Modelling calibration errors in cbc waveforms,” (2014), https: //dcc.ligo.org/LIGO-T1400682/public
work page 2014
-
[62]
This example is found in the Bilby git repository at https://git.ligo.org/lscsoft/bilby/blob/master/ examples/injection_examples/basic_tutorial.py
- [63]
-
[64]
B. P. Abbott et al. , Astrophys. J. 848, L12 (2017)
work page 2017
-
[65]
S. De, D. Finstad, J. M. Lattimer, D. A. Brown, E. Berger, and C. M. Biwer, (2018), arXiv:1804.08583
work page Pith review arXiv 2018
-
[66]
This example is found in the Bilby git repository at https://git.ligo.org/lscsoft/bilby/blob/master/ examples/injection_examples/binary_neutron_ star_example.py
-
[67]
´E. ´E. Flanagan and T. Hinderer, Phys. Rev. D 77, 021502 (2008). 14
work page 2008
-
[68]
B. D. Lackey and L. Wade, Phys. Rev. D 91, 043002 (2015)
work page 2015
-
[69]
M. E. Lower, E. Thrane, P. D. Lasky, and R. Smith, Phys. Rev. D 98, 083028 (2018)
work page 2018
-
[70]
B. R. Iyer, T. Souradeep, C. S. Unnikrishnan, S. Dhurandhar, S. Raja, A. Sengupta, and I. Con- sortium), “Ligo-india, proposal of the consortium for indian initiative in gravitational-wave observa- tions (indigo),” (2011), https://dcc.ligo.org/ LIGO-M1100296/public
work page 2011
- [71]
- [72]
-
[73]
B. P. Abbott et al. , Classical Quantum Gravity 34, 44001 (2017)
work page 2017
- [74]
-
[75]
https://git.ligo.org/gwinc/pygwinc
-
[76]
This example is found in the Bilby git repository at https://git.ligo.org/lscsoft/bilby/blob/ master/examples/injection_examples/Australian_ detector.py
- [77]
- [78]
- [79]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.