A practical Bayesian method for gravitational-wave ringdown analysis with multiple modes
Pith reviewed 2026-05-23 03:51 UTC · model grok-4.3
The pith
FIREFLY analytically marginalizes amplitudes and phases of multiple quasi-normal modes to cut ringdown inference time from hours to minutes while keeping posteriors and evidence consistent.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
FIREFLY analytically marginalizes the amplitude and phase parameters of QNMs to reduce the computational cost and speed up the full-parameter inference from hours to minutes, while achieving consistent posterior and evidence. Rigorously based on the principle of Bayesian inference and importance sampling, the method is statistically interpretable, flexible in prior choice, and compatible with various advanced sampling techniques.
What carries the argument
The FIREFLY algorithm, which analytically marginalizes the amplitude and phase parameters of quasi-normal modes and applies importance sampling to the remaining parameters.
If this is right
- The computational speedup grows larger as the number of included quasi-normal modes increases.
- The method remains compatible with any advanced sampling technique used on the reduced parameter space.
- Prior choices on the remaining parameters stay fully flexible without altering the marginalization step.
- The approach supplies a statistically interpretable route to faster ringdown analyses that scales with detector sensitivity.
Where Pith is reading between the lines
- The same marginalization step could be adapted to other gravitational-wave signal models that separate linear amplitudes from nonlinear parameters.
- If the consistency holds under realistic noise, the method could support near-real-time ringdown parameter estimation during observing runs.
- Extending the importance-sampling weights to include detector calibration uncertainties would test whether the speedup survives additional systematic parameters.
Load-bearing premise
That analytical marginalization of amplitudes and phases followed by importance sampling produces posteriors and evidences that match full nested sampling or MCMC across the full parameter space and all prior choices.
What would settle it
A side-by-side run of FIREFLY and full nested sampling on the same simulated ringdown signal that contains three or more quasi-normal modes, checking whether the recovered log-evidence and marginal posteriors on the mode frequencies and damping times agree within sampling noise.
read the original abstract
Gravitational-wave (GW) ringdown signals from black holes (BHs) encode crucial information about the gravitational dynamics in the strong-field regime, which offers unique insights into BH properties. In the future, the improving sensitivity of GW detectors is to enable the extraction of multiple quasi-normal modes (QNMs) from ringdown signals. However, incorporating multiple modes drastically enlarges the parameter space, posing computational challenges to data analysis. Inspired by the $F$-statistic method in the continuous GW searches, we develope an algorithm, dubbed as FIREFLY, for accelerating the ringdown signal analysis. FIREFLY analytically marginalizes the amplitude and phase parameters of QNMs to reduce the computational cost and speed up the full-parameter inference from hours to minutes, while achieving consistent posterior and evidence. The acceleration becomes more significant when more QNMs are considered. Rigorously based on the principle of Bayesian inference and importance sampling, our method is statistically interpretable, flexible in prior choice, and compatible with various advanced sampling techniques, providing a new perspective for accelerating future GW data analysis.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces FIREFLY, an algorithm for Bayesian ringdown analysis of gravitational-wave signals that analytically marginalizes the linear amplitude and phase parameters of multiple quasi-normal modes (QNMs) and employs importance sampling over the nonlinear parameters (frequencies, damping times, etc.). It claims this reduces full-parameter inference from hours to minutes while producing posteriors and evidences consistent with standard nested sampling or MCMC, with the speedup increasing for larger numbers of modes. The method is presented as statistically interpretable, prior-flexible, and compatible with existing samplers, building on F-statistic ideas.
Significance. If the consistency of posteriors and evidences holds under the stated conditions, the approach offers a practical route to multi-mode ringdown inference at scale, which will become essential as detector sensitivity improves and higher-order QNMs become measurable. The analytic marginalization step is a standard consequence of the Gaussian likelihood for stationary noise, and the importance-sampling wrapper preserves the Bayesian interpretation without introducing new fitted quantities.
major comments (2)
- [Abstract] Abstract: the central claim that the method 'achieves consistent posterior and evidence' is load-bearing yet unsupported by any reported diagnostics (effective sample size, weight histograms, KL divergence to reference posteriors, or direct comparison tables) across mode counts, prior widths, or frequency separations. Without these, it is impossible to assess whether the importance-sampling approximation remains reliable when modes are closely spaced or priors are broad.
- [Abstract] The manuscript does not specify how the proposal distribution for importance sampling is constructed or adapted when the number of QNMs increases; if the proposal is derived from a lower-mode run or from a fixed approximation, the variance of the weights can grow rapidly, undermining the claimed consistency for the full multi-mode case.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive report. The two major comments identify important gaps in the validation and documentation of the importance-sampling step. We address each point below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that the method 'achieves consistent posterior and evidence' is load-bearing yet unsupported by any reported diagnostics (effective sample size, weight histograms, KL divergence to reference posteriors, or direct comparison tables) across mode counts, prior widths, or frequency separations. Without these, it is impossible to assess whether the importance-sampling approximation remains reliable when modes are closely spaced or priors are broad.
Authors: We agree that the current manuscript does not report the quantitative diagnostics listed. In the revised version we will add (i) effective sample sizes and weight histograms for representative cases, (ii) KL-divergence values between the importance-sampled and nested-sampling posteriors, and (iii) summary tables comparing posterior means, credible intervals and log-evidence differences across mode counts, prior widths and frequency separations. These additions will be placed in a new subsection of the results and referenced from the abstract. revision: yes
-
Referee: [Abstract] The manuscript does not specify how the proposal distribution for importance sampling is constructed or adapted when the number of QNMs increases; if the proposal is derived from a lower-mode run or from a fixed approximation, the variance of the weights can grow rapidly, undermining the claimed consistency for the full multi-mode case.
Authors: The referee correctly notes that the construction and adaptation of the proposal distribution are not described. We will revise the methods section to specify (a) how the proposal is obtained (typically from a preliminary run with fewer modes or from an analytic Gaussian approximation to the marginalized likelihood), (b) the procedure used to update the proposal when additional modes are introduced, and (c) any safeguards (e.g., weight truncation or adaptive resampling) employed to control weight variance. We will also add a brief discussion of the scaling of weight variance with mode number. revision: yes
Circularity Check
No circularity: standard analytic marginalization + importance sampling
full rationale
The paper derives FIREFLY from the standard Bayesian marginalization of linear amplitude/phase parameters under Gaussian noise (inspired by the external F-statistic) followed by importance sampling over nonlinear parameters. This is a direct application of established inference principles with no self-definitional reduction, no fitted input renamed as prediction, and no load-bearing self-citation chain. The consistency claim is presented as a consequence of the method rather than an input; the derivation remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Ringdown signals are linear superpositions of quasi-normal modes with independent amplitudes and phases.
Forward citations
Cited by 3 Pith papers
-
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.
-
FluxMC: Rapid and High-Fidelity Inference for Space-Based Gravitational-Wave Observations
FluxMC integrates flow matching with parallel tempering MCMC to converge in under five hours on high-fidelity IMRPhenomHM waveforms for massive black hole binaries, where standard methods fail after hundreds of hours ...
-
A Robust and Efficient F-statistic-based Framework for Consistent Bayesian Inference of Compact Binary Coalescences
F-statistic framework analytically maximizes over distance and polarization to enable faster Bayesian inference of compact binary coalescences with a new evidence formulation that matches full frequency-domain results...
Reference graph
Works this paper leans on
-
[1]
Observation of Gravitational Waves from a Binary Black Hole Merger
Abbott, B.P., et al. : Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett. 116(6), 061102 (2016) https://doi.org/10.1103/PhysRevLett.116.061102 arXiv:1602.03837 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.116.061102 2016
-
[2]
Tests of general relativity with GW150914
Abbott, B.P., et al. : Tests of general relativity with GW150914. Phys. Rev. Lett. 116(22), 221101 (2016) https://doi.org/10.1103/PhysRevLett.116.221101 arXiv:1602.03841 [gr-qc]. [Erratum: Phys.Rev.Lett. 121, 129902 (2018)]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.116.221101 2016
-
[3]
Abbott, R., et al.: Tests of General Relativity with GWTC-3 (2021) arXiv:2112.06861 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[4]
A gravitational-wave standard siren measurement of the Hubble constant
Abbott, B.P., et al. : A gravitational-wave standard siren measurement of the Hubble constant. Nature 551(7678), 85–88 (2017) https://doi.org/10.1038/nature24471 arXiv:1710.05835 [astro-ph.CO]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1038/nature24471 2017
-
[5]
GW170817: Measurements of Neutron Star Radii and Equation of State
Abbott, B.P., et al. : GW170817: Measurements of neutron star radii and equation of state. Phys. Rev. Lett. 121(16), 161101 (2018) https://doi.org/10.1103/PhysRevLett.121.161101 arXiv:1805.11581 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.121.161101 2018
-
[6]
The population of merging compact binaries inferred using gravitational waves through GWTC-3
Abbott, R., et al. : Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3. Phys. Rev. X 13(1), 011048 (2023) https://doi.org/10.1103/PhysRevX.13.011048 arXiv:2111.03634 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevx.13.011048 2023
-
[7]
Isi, M., Giesler, M., Farr, W.M., Scheel, M.A., Teukolsky, S.A.: Testing the no-hair theorem with GW150914. Phys. Rev. Lett. 123(11), 111102 (2019) https://doi.org/10.1103/PhysRevLett.123.111102 arXiv:1905.00869 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.123.111102 2019
-
[8]
Isi, M., Farr, W.M., Giesler, M., Scheel, M.A., Teukolsky, S.A.: Testing the Black-Hole Area Law with GW150914. Phys. Rev. Lett. 127(1), 011103 (2021) https://doi.org/10.1103/PhysRevLett.127.011103 arXiv:2012.04486 [gr-qc]
-
[9]
Cardoso, V., Franzin, E., Pani, P.: Is the gravitational-wave ringdown a probe of the event horizon? Phys. Rev. Lett. 116(17), 171101 (2016) https://doi.org/10.1103/PhysRevLett.116.171101 arXiv:1602.07309 [gr-qc]. [Erratum: Phys.Rev.Lett. 117, 089902 (2016)]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.116.171101 2016
-
[10]
Teukolsky, S.A.: Perturbations of a rotating black hole. 1. Fundamental equations for gravitational elec- tromagnetic and neutrino field perturbations. Astrophys. J. 185, 635–647 (1973) https://doi.org/10.1086/ 152444
work page 1973
-
[11]
Berti, E., Cardoso, J., Cardoso, V., Cavaglia, M.: Matched-filtering and parameter estimation of ringdown waveforms. Phys. Rev. D 76, 104044 (2007) https://doi.org/10.1103/PhysRevD.76.104044 arXiv:0707.1202 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.76.104044 2007
-
[12]
Dreyer, O., Kelly, B.J., Krishnan, B., Finn, L.S., Garrison, D., Lopez-Aleman, R.: Black hole spectroscopy: Testing general relativity through gravitational wave observations. Class. Quant. Grav.21, 787–804 (2004) https://doi.org/10.1088/0264-9381/21/4/003 arXiv:gr-qc/0309007
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/21/4/003 2004
-
[13]
Christensen, N., Meyer, R.: Markov chain Monte Carlo methods for Bayesian gravitational radiation data analysis. Phys. Rev. D 58, 082001 (1998) https://doi.org/10.1103/PhysRevD.58.082001
-
[14]
Sharma, S.: Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy. Ann. Rev. Astron. Astrophys. 55, 213–259 (2017) https://doi.org/10.1146/annurev-astro-082214-122339 arXiv:1706.01629 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1146/annurev-astro-082214-122339 2017
-
[15]
2004, in American Institute of Physics Conference Series, Vol
Skilling, J.: Nested Sampling. AIP Conf. Proc. 735(1), 395 (2004) https://doi.org/10.1063/1.1835238
-
[16]
Bayesian Analysis 1(4), 833–859 (2006) https://doi.org/10.1214/06-BA127
Skilling, J.: Nested sampling for general Bayesian computation. Bayesian Analysis 1(4), 833–859 (2006) https://doi.org/10.1214/06-BA127
-
[17]
Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO
Reitze, D., et al. : Cosmic Explorer: The U.S. Contribution to Gravitational-Wave Astronomy beyond LIGO. Bull. Am. Astron. Soc. 51(7), 035 (2019) arXiv:1907.04833 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[18]
The US Program in Ground-Based Gravitational Wave Science: Contribution from the LIGO Laboratory
Reitze, D., et al. : The US Program in Ground-Based Gravitational Wave Science: Contribution from the LIGO Laboratory. Bull. Am. Astron. Soc. 51, 141 (2019) arXiv:1903.04615 [astro-ph.IM] 19
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[19]
: The Einstein Telescope: A third-generation gravitational wave observatory
Punturo, M., et al. : The Einstein Telescope: A third-generation gravitational wave observatory. Class. Quant. Grav. 27, 194002 (2010) https://doi.org/10.1088/0264-9381/27/19/194002
-
[20]
Hild, S., et al.: Sensitivity Studies for Third-Generation Gravitational Wave Observatories. Class. Quant. Grav. 28, 094013 (2011) https://doi.org/10.1088/0264-9381/28/9/094013 arXiv:1012.0908 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/28/9/094013 2011
-
[21]
Laser Interferometer Space Antenna
Amaro-Seoane, P., et al. : Laser Interferometer Space Antenna. e-prints (2017) https://doi.org/10.48550/ arXiv.1702.00786 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[22]
Hu, W.-R., Wu, Y.-L.: The Taiji Program in Space for gravitational wave physics and the nature of gravity. Natl. Sci. Rev. 4(5), 685–686 (2017) https://doi.org/10.1093/nsr/nwx116
-
[23]
TianQin: a space-borne gravitational wave detector
Luo, J., et al. : TianQin: a space-borne gravitational wave detector. Class. Quant. Grav. 33(3), 035010 (2016) https://doi.org/10.1088/0264-9381/33/3/035010 arXiv:1512.02076 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/33/3/035010 2016
-
[24]
Bhagwat, S., Pacilio, C., Barausse, E., Pani, P.: Landscape of massive black-hole spectroscopy with LISA and the Einstein Telescope. Phys. Rev. D 105(12), 124063 (2022) https://doi.org/10.1103/PhysRevD. 105.124063 arXiv:2201.00023 [gr-qc]
- [25]
-
[26]
Pitte, C., Baghi, Q., Besan¸ con, M., Petiteau, A.: Exploring tests of the no-hair theorem with LISA. Phys. Rev. D 110(10), 104003 (2024) https://doi.org/10.1103/PhysRevD.110.104003 arXiv:2406.14552 [gr-qc]
-
[27]
Berti, E., Cardoso, V., Starinets, A.O.: Quasinormal modes of black holes and black branes. Class. Quant. Grav. 26, 163001 (2009) https://doi.org/10.1088/0264-9381/26/16/163001 arXiv:0905.2975 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/26/16/163001 2009
-
[28]
Jaranowski, P., Krolak, A., Schutz, B.F.: Data analysis of gravitational - wave signals from spinning neutron stars. 1. The Signal and its detection. Phys. Rev. D 58, 063001 (1998) https://doi.org/10.1103/ PhysRevD.58.063001 arXiv:gr-qc/9804014
work page internal anchor Pith review Pith/arXiv arXiv 1998
- [29]
-
[30]
Wang, H.-T., Yim, G., Chen, X., Shao, L.: Gravitational Wave Ringdown Analysis Using the F -statistic. Astrophys. J. 974(2), 230 (2024) https://doi.org/10.3847/1538-4357/ad7096 arXiv:2409.00970 [gr-qc]
-
[31]
Prix, R., Krishnan, B.: Targeted search for continuous gravitational waves: Bayesian versus maximum- likelihood statistics. Class. Quant. Grav. 26, 204013 (2009) https://doi.org/10.1088/0264-9381/26/20/ 204013 arXiv:0907.2569 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0264-9381/26/20/ 2009
-
[32]
LIGO Document, https://dcc.ligo.org/LIGO-T1500618/public
Prix, R.: Bayesian QNM search on black hole ringdown modes (applied to GW150914). LIGO Document, https://dcc.ligo.org/LIGO-T1500618/public
-
[33]
Gabbard, H., Messenger, C., Heng, I.S., Tonolini, F., Murray-Smith, R.: Bayesian parameter estimation using conditional variational autoencoders for gravitational-wave astronomy. Nature Phys.18(1), 112–117 (2022) https://doi.org/10.1038/s41567-021-01425-7 arXiv:1909.06296 [astro-ph.IM]
- [34]
-
[35]
Dax, M., Green, S.R., Gair, J., P¨ urrer, M., Wildberger, J., Macke, J.H., Buonanno, A., Sch¨ olkopf, B.: Neu- ral Importance Sampling for Rapid and Reliable Gravitational-Wave Inference. Phys. Rev. Lett. 130(17), 171403 (2023) https://doi.org/10.1103/PhysRevLett.130.171403 arXiv:2210.05686 [gr-qc]
-
[36]
Pacilio, C., Bhagwat, S., Cotesta, R.: Simulation-based inference of black hole ringdowns in the time domain. Phys. Rev. D 110(8), 083010 (2024) https://doi.org/10.1103/PhysRevD.110.083010 arXiv:2404.11373 [gr-qc] 20
-
[37]
Springer (2004), https://link.springer.com/ book/10.1007/978-1-4757-4145-2
Robert, C.P., Casella, G.: Monte carlo statistical methods. Springer (2004), https://link.springer.com/ book/10.1007/978-1-4757-4145-2
-
[38]
: The SXS Collaboration catalog of binary black hole simulations
Boyle, M., et al. : The SXS Collaboration catalog of binary black hole simulations. Class. Quant. Grav. 36(19), 195006 (2019) https://doi.org/10.1088/1361-6382/ab34e2 arXiv:1904.04831 [gr-qc]
-
[39]
Giesler, M., Isi, M., Scheel, M.A., Teukolsky, S.: Black Hole Ringdown: The Importance of Overtones. Phys. Rev. X 9(4), 041060 (2019) https://doi.org/10.1103/PhysRevX.9.041060 arXiv:1903.08284 [gr-qc]
-
[40]
Statistics and Computing 29(5), 891 (2019) https: //doi.org/10.1007/s11222-018-9844-0
Higson, E., Handley, W., Hobson, M., Lasenby, A.: Dynamic nested sampling: an improved algorithm for parameter estimation and evidence calculation. Statistics and Computing 29(5), 891 (2019) https: //doi.org/10.1007/s11222-018-9844-0
-
[41]
Computational Optimal Transport: With Ap- plications to Data Science
Peyr´ e, G., Cuturi, M.: Computational optimal transport: With applications to data science. Foundations and Trends in Machine Learning 11(5-6), 355–607 (2019) https://doi.org/10.1561/2200000073
-
[42]
Wang, Y., Shang, Y., Babak, S.: EMRI data analysis with a phenomenological waveform. Phys. Rev. D 86, 104050 (2012) https://doi.org/10.1103/PhysRevD.86.104050 arXiv:1207.4956 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.86.104050 2012
-
[43]
Hu, Q., Veitch, J.: Rapid Premerger Localization of Binary Neutron Stars in Third-generation Gravitational-wave Detectors. Astrophys. J. Lett. 958(2), 43 (2023) https://doi.org/10.3847/2041-8213/ ad0ed4 arXiv:2309.00970 [gr-qc]
-
[44]
Berti, E., Klein, A.: Mixing of spherical and spheroidal modes in perturbed Kerr black holes. Phys. Rev. D 90(6), 064012 (2014) https://doi.org/10.1103/PhysRevD.90.064012 arXiv:1408.1860 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.90.064012 2014
-
[45]
Finn, L.S.: Detection, measurement and gravitational radiation. Phys. Rev. D 46, 5236–5249 (1992) https://doi.org/10.1103/PhysRevD.46.5236 arXiv:gr-qc/9209010
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.46.5236 1992
-
[46]
Loredo, T.J., Wolpert, R.L.: Bayesian inference: more than Bayes’s theorem. Frontiers in Astronomy and Space Sciences 11, 1326926 (2024) https://doi.org/10.3389/fspas.2024.1326926 arXiv:2406.18905 [stat.ME]
-
[47]
Thrane, E., Talbot, C.: An introduction to Bayesian inference in gravitational-wave astronomy: parameter estimation, model selection, and hierarchical models. Publ. Astron. Soc. Austral. 36, 010 (2019) https: //doi.org/10.1017/pasa.2019.2 arXiv:1809.02293 [astro-ph.IM]. [Erratum: Publ.Astron.Soc.Austral. 37, e036 (2020)]
-
[48]
WIREs Computational Statistics 2(1), 54–60 (2010) https://doi.org/10.1002/wics.56
Tokdar, S.T., Kass, R.E.: Importance sampling: a review. WIREs Computational Statistics 2(1), 54–60 (2010) https://doi.org/10.1002/wics.56
-
[49]
Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy
Ashton, G., et al. : BILBY: A user-friendly Bayesian inference library for gravitational-wave astron- omy. Astrophys. J. Suppl. 241(2), 27 (2019) https://doi.org/10.3847/1538-4365/ab06fc arXiv:1811.02042 [astro-ph.IM] 21
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4365/ab06fc 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.