Pith. sign in

REVIEW 2 major objections 6 minor 29 references

Even with a perfect model, finite SNR systematically stretches the damping time of ringdown signals; stacking ~100 moderate events without care can fake a Kerr violation.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 06:55 UTC pith:H3VN4FJI

load-bearing objection Clean O(ρ^{-2}) calculation showing that τ is biased high without Q-suppression; the N≳100 catalog warning is real inside the stated model and worth acting on. the 2 major comments →

arxiv 2607.05486 v1 pith:H3VN4FJI submitted 2026-07-06 gr-qc astro-ph.HE

Finite signal-to-noise ratio bias in parameter estimation for damped oscillations: cautionary remark about catalog-level black-hole spectroscopy

Koutarou Kyutoku This is my paper
classification gr-qc astro-ph.HE
keywords black-hole spectroscopyfinite SNR biasdamping timeFisher matrixquasinormal modesKerr hypothesiscatalog combinationparameter estimation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Black-hole spectroscopy tests whether a remnant is a Kerr black hole by measuring the frequency and damping time of its ringdown. This paper shows that, for the simplest possible model—a single damped sinusoid with fixed start time in white noise—the damping time is biased high by an amount of order 1 over SNR squared that is not suppressed by the quality factor. Two independent sources produce the bias: the usual decreasing prior on amplitude, and a geometric correction to the likelihood surface that survives even with a flat prior. Both survive when the model and analysis are perfect. Because statistical errors shrink only as the square root of the number of events while the bias does not, a naïve catalog combination of roughly a hundred events at SNR around 10 can appear to reject the Kerr hypothesis. The paper supplies simple reparametrizations and prior adjustments that cancel the leading bias event by event.

Core claim

At finite signal-to-noise ratio the damping time of a pure damped sinusoid receives an O(ρ^{-2}) positive bias that is independent of the quality factor Q, arising both from the prior gradient on amplitude and from the curvature of the likelihood; the bias therefore accumulates under naïve catalog stacking and can produce a false rejection of the Kerr hypothesis with ≳100 events of ρ≈10.

What carries the argument

The beyond-Fisher expansion of the maximum-likelihood and Bayesian estimators to O(ρ^{-2}), evaluated on the four-parameter damped-sinusoid waveform, which isolates the Q-independent pieces of both the prior-gradient bias and the likelihood-geometry bias for the damping time.

Load-bearing premise

The quantitative size and sign of the bias are derived for a pure single-mode damped sinusoid with fixed start time in stationary white Gaussian noise; any real ringdown that contains overtones, higher harmonics or colored noise could change the result.

What would settle it

Inject pure single-mode damped sinusoids of known parameters at ρ=10 into white Gaussian noise, recover τ with both maximum-likelihood and Bayesian analyses under the stated priors, and check whether the ensemble-averaged bias matches the analytic O(ρ^{-2}) expressions; a null result would falsify the central claim.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Catalog-level black-hole spectroscopy that simply multiplies or averages marginalized posteriors on τ will systematically prefer longer damping times and can falsely exclude Kerr with ≳100 moderate events.
  • Adopting the damping rate γ=1/τ instead of τ removes the leading likelihood-geometry bias for maximum-likelihood estimators.
  • Choosing a power-law prior p(τ)∝τ^{-n_τ} with n_τ=(n_A+3)/2 cancels both the prior-gradient and likelihood-geometry biases for individual events.
  • Focusing only on high-SNR ringdowns avoids the finite-ρ bias but leaves the analysis more vulnerable to model incompleteness.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Any hierarchical or stacked analysis that treats point estimates or one-dimensional marginals as independent measurements of τ will inherit the same accumulation problem, not only black-hole spectroscopy.
  • The same Q-independent bias mechanism should appear in other damped-oscillator problems outside gravitational waves (e.g., free-induction-decay spectroscopy) whenever amplitude and lifetime are estimated jointly at moderate SNR.
  • If real catalogs already show a mild excess of long damping times, a first diagnostic is to re-analyze the same events under the bias-canceling prior and under the γ parametrization; disappearance of the excess would favor the finite-SNR explanation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. The paper computes O(ρ^{-2}) biases in matched-filter parameter estimation for a single-mode damped sinusoid (fixed start time, white Gaussian noise) as a minimal model of black-hole ringdown spectroscopy. Using the beyond-Fisher expansion, it shows that the damping time τ receives a positive bias that is not suppressed by the quality factor Q, arising from (i) the prior-gradient shift when the amplitude prior is decreasing (Eq. 16) and (ii) the likelihood-geometry bias of both the maximum-likelihood estimator (Eq. 21) and the Bayesian posterior mean (Eq. 30). Statistical errors and covariances are given explicitly (Eqs. 5–14), cross-checked under two φ-averaging schemes (main text and App. A). The paper argues that naive combination of N≳100 events with ρ=10 can produce a false catalog-level rejection of the Kerr hypothesis, and proposes simple mitigations (reparametrization to γ=1/τ for ML; adjusted power-law priors on τ).

Significance. If the analytic results hold under the stated idealizations, the paper supplies a concrete, previously under-emphasized systematic for catalog-level black-hole spectroscopy: a Q-independent positive bias in τ that can accumulate under naive joint inference even when the waveform model and single-event analysis are perfect. The derivations are transparent, parameter-free within the model, and algebraically cross-checked (two averaging procedures plus a referenced Mathematica notebook). The mitigation strategies are immediately usable. The work is therefore a useful caution for ongoing LVK-style remnant tests and for any pipeline that stacks ringdown posteriors or point estimates. The idealization (single mode, fixed start, white noise) is stated clearly and does not invalidate the formal result inside that model.

major comments (2)
  1. [Abstract; Secs. III–IV] Abstract and Secs. III–IV: the headline claim that catalog-level spectroscopy can report a false Kerr violation with ≳100 events at ρ=10 mixes two distinct mechanisms and combination assumptions. Prior-gradient bias (Fig. 2, n_A=1) yields N≳200 in the (f,τ) plane; ML likelihood-geometry bias alone yields N≳25 (text after Eq. 21). The abstract’s single ≳100 figure should be tied explicitly to a stated combination protocol (e.g., stacking of τ medians vs. full-likelihood hierarchical combination) so that the quantitative warning is falsifiable rather than illustrative.
  2. [Sec. IV] Sec. IV, paragraph after Eq. 26 and discussion of catalog combination: the paper asserts that appropriate Bayesian combination of full individual likelihoods suppresses the likelihood-geometry bias by N^{-1}, while “naive” joint inferences accumulate it. This distinction is load-bearing for the cautionary claim, yet “naive” is only exemplified (“combining only marginalized likelihood or point estimates”). A short, operational definition—or a one-paragraph sketch of how a hierarchical analysis of the form used in Refs. [27,28] would cancel the O(ρ^{-2}) term—would make the recommendation actionable for existing pipelines.
minor comments (6)
  1. [Sec. II] Sec. II, Eq. (3) and surrounding text: the phase-dependent factor in ρ^{2} is kept, then quantities are φ-averaged after normalizing by ρ^{2}. A one-sentence reminder that the leading-in-Q terms are robust under the alternative averaging of App. A would help readers who skip the appendix.
  2. [Fig. 1; Sec. II] Fig. 1 caption and text: Q^{2}≈10 for χ≈0.686 is used throughout the quantitative estimates. Stating once that the Q-independent leading bias terms dominate already for Q^{2}≳10 (and only strengthen for higher spin) would clarify that the caution is not restricted to equal-mass nonspinning remnants.
  3. [Sec. III] Sec. III, Eq. (16): the power-law indices n_A and n_τ are introduced without a brief statement of the prior support (e.g., that the power law is assumed only over the region where the likelihood is appreciable). A short clause would avoid confusion with improper priors.
  4. [Sec. IV.A] Sec. IV.A, reparametrization to γ=1/τ: the expansion (22)–(25) is clear, but it would help to note explicitly that the Jacobian of the reparametrization does not alter the Fisher information for the statistical error, only the O(ρ^{-2}) bias via the quadratic term.
  5. [App. A] References: the companion Mathematica notebook is mentioned only in App. A. Adding a data-availability or software note (or arXiv ancillary file pointer) would improve reproducibility.
  6. [Sec. II; Sec. IV] Typographical: “ST A TISTICAL” in the Sec. II heading appears to contain an extraneous space; “LIKELIHOOD-GEOMETR Y” in Sec. IV likewise.

Circularity Check

0 steps flagged

No circularity: O(ρ^{-2}) bias formulas are direct expansions of the likelihood and prior inside an explicitly idealized model; the Kerr comparison is an external benchmark.

full rationale

The paper derives the prior-gradient bias (Eq. 16) and the likelihood-geometry biases for the maximum-likelihood estimator (Eq. 21) and Bayesian mean (Eq. 30) by standard beyond-Fisher expansions of a single-mode damped sinusoid (Eq. 1) in white Gaussian noise. Every algebraic step is self-contained: the Fisher matrix, its inverse, the second-derivative contractions that produce Δ heta_ML and Δ heta_B, and the prior-gradient shift δ heta^a = C^{ab} abla_b ln p are computed explicitly and averaged over phase. None of these expressions is obtained by fitting a free parameter to data and then re-labeling the fit as a prediction, nor by importing a uniqueness theorem or ansatz from the author’s own prior work. The quantitative N ≳ 100 warning is simply the ratio of the derived O(ρ^{-2}) bias to the O(ρ^{-1}) statistical error under the stated idealizations; it does not close a definitional loop. Self-citations are limited to standard references (Vallisneri 2008, Cutler & Flanagan 1994, Berti et al. 2006) that supply the general expansion formulas or the Q(χ) fit; they are not load-bearing for the sign or Q-independence of the τ bias. The idealized model is openly declared as a limitation, not smuggled in as a hidden premise. Consequently the derivation chain contains no circular step.

Axiom & Free-Parameter Ledger

0 free parameters · 4 axioms · 0 invented entities

The calculation rests on standard matched-filter and Fisher-matrix technology plus a handful of domain idealizations that the author states explicitly. No free parameters are fitted to data; Q and ρ are free variables. No new physical entities are introduced.

axioms (4)
  • domain assumption Detector noise is stationary, Gaussian and white (constant Sn) over the frequency band of the damped sinusoid.
    Stated in Sec. II; used to define the inner product (Eq. 2) and all subsequent Fisher-matrix elements.
  • domain assumption The waveform is exactly a single-mode damped sinusoid with fixed start time t=0 (Eq. 1); no overtones, higher harmonics or nonlinear modes.
    Explicit modeling choice in Sec. II; the author notes that real ringdowns are more complicated.
  • standard math Beyond-Fisher expansion of the likelihood truncated at O(ρ^{-2}) is sufficient to capture the leading bias.
    Taken from Vallisneri (2008) and applied throughout Secs. III–IV.
  • domain assumption Priors on amplitude are power-law (p(A)∝A^{-n_A}) with n_A typically positive.
    Used to obtain the prior-gradient bias (Eq. 16); log-flat (n_A=1) is the common choice cited from LVK analyses.

pith-pipeline@v1.1.0-grok45 · 16139 in / 2467 out tokens · 22134 ms · 2026-07-11T06:55:31.877173+00:00 · methodology

0 comments
read the original abstract

We investigate biases in parameter estimation for damped oscillations motivated by applications to black-hole spectroscopy in gravitational-wave physics. Focusing on the simplest model of a single-mode damped sinusoid with a fixed start time in white noise, we show that, at finite signal-to-noise ratio \rho, the damping time is biased toward larger values without being suppressed by the quality factor for two reasons. One is the gradient of the prior, through which the damping time is affected by the typically decreasing prior on the amplitude. The other is a higher-order finite-\rho correction to the likelihood geometry. These biases arise even if the model and analysis are appropriate. Moreover, they could be exaggerated in naive joint inferences from catalog events. Quantitatively, if estimates from multiple events with \rho=10 are combined without due care, the catalog-level black-hole spectroscopy could report false violation of the Kerr hypothesis with >~100 events. We also propose simple strategies to mitigate these biases at the level of individual events.

Figures

Figures reproduced from arXiv: 2607.05486 by Koutarou Kyutoku.

Figure 1
Figure 1. Figure 1: FIG. 1. Square of the quality factor [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Contour enclosing 68 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

29 extracted references · 21 linked inside Pith

  1. [1]

    Dreyer, B

    O. Dreyer, B. Kelly, B. Krishnan, L. S. Finn, D. Garrison, and R. Lopez-Aleman, Black-hole spectroscopy: test- ing general relativity through gravitational-wave obser- vations, Classical and Quantum Gravity21, 787 (2004), arXiv:gr-qc/0309007 [gr-qc]

    Pith/arXiv arXiv 2004

  2. [2]

    Berti, V

    E. Berti, V. Cardoso, G. Carullo, J. Abedi, N. Af- shordi, S. Albanesi, V. Baibhav, S. Bhagwat, J. L. Bl´ azquez-Salcedo, B. Bonga, B. Bucciotti, G. Caneva Santoro, P. A. Cano, C. Capano, M. Ho-Yeuk Cheung, C. Chirenti, G. B. Cook, A. Ka-Wai Chung, M. De Ami- cis, K. Destounis, O. J. C. Dias, W. Del Pozzo, F. Duque, W. M. Farr, E. Finch, N. Franchini, K. ...

    Pith/arXiv arXiv 2025

  3. [3]

    Berti, V

    E. Berti, V. Cardoso, and A. O. Starinets, TOPI- CAL REVIEW: Quasinormal modes of black holes and black branes, Classical and Quantum Gravity26, 163001 (2009), arXiv:0905.2975 [gr-qc]

    Pith/arXiv arXiv 2009

  4. [4]

    The LIGO Scientific Collaboration and the Virgo Col- laboration, Tests of General Relativity with GW150914, Phys. Rev. Lett.116, 221101 (2016), arXiv:1602.03841 [gr-qc]

    Pith/arXiv arXiv 2016

  5. [5]

    The LIGO Scientific Collaboration and the Virgo Col- laboration, Tests of general relativity with binary black holes from the second LIGO-Virgo gravitational-wave transient catalog, Phys. Rev. D103, 122002 (2021), arXiv:2010.14529 [gr-qc]

    Pith/arXiv arXiv 2021

  6. [6]

    The LIGO Scientific Collaboration, the Virgo Collabora- tion, and the KAGRA Collaboration, Tests of general rel- ativity with GWTC-3, Phys. Rev. D112, 084080 (2025), arXiv:2112.06861 [gr-qc]

    Pith/arXiv arXiv 2025

  7. [7]

    The LIGO Scientific Collaboration, the Virgo Collabora- tion, and the KAGRA Collaboration, GWTC-4.0: Tests of General Relativity. III. Tests of the Remnants, arXiv e-prints , arXiv:2603.19021 (2026), arXiv:2603.19021 [gr- qc]

    Pith/arXiv arXiv 2026

  8. [8]

    Pacilio, D

    C. Pacilio, D. Gerosa, and S. Bhagwat, Catalog variance of testing general relativity with gravitational-wave data, Phys. Rev. D109, L081302 (2024), arXiv:2310.03811 [gr- qc]

    Pith/arXiv arXiv 2024

  9. [9]

    Thrane, P

    E. Thrane, P. D. Lasky, and Y. Levin, Challenges for testing the no-hair theorem with current and planned gravitational-wave detectors, Phys. Rev. D96, 102004 (2017), arXiv:1706.05152 [gr-qc]

    Pith/arXiv arXiv 2017

  10. [10]

    Baibhav, E

    V. Baibhav, E. Berti, V. Cardoso, and G. Khanna, Black hole spectroscopy: Systematic errors and ring- down energy estimates, Phys. Rev. D97, 044048 (2018), arXiv:1710.02156 [gr-qc]

    Pith/arXiv arXiv 2018

  11. [11]

    Isi and W

    M. Isi and W. M. Farr, Analyzing black-hole ring- downs, arXiv e-prints , arXiv:2107.05609 (2021), arXiv:2107.05609 [gr-qc]

    Pith/arXiv arXiv 2021

  12. [12]

    S. H. V¨ olkel and A. Dhani, Quantifying systematic biases in black hole spectroscopy, Phys. Rev. D112, 084076 (2025), arXiv:2507.22122 [gr-qc]

    arXiv 2025

  13. [13]

    Vallisneri, Use and abuse of the Fisher infor- mation matrix in the assessment of gravitational- wave parameter-estimation prospects, Phys

    M. Vallisneri, Use and abuse of the Fisher infor- mation matrix in the assessment of gravitational- wave parameter-estimation prospects, Phys. Rev. D77, 042001 (2008), arXiv:gr-qc/0703086 [gr-qc]

    Pith/arXiv arXiv 2008

  14. [14]

    Cutler and ´E

    C. Cutler and ´E. E. Flanagan, Gravitational waves from merging compact binaries: How accurately can one extract the binary’s parameters from the inspiral waveform?, Phys. Rev. D49, 2658 (1994), arXiv:gr- qc/9402014 [gr-qc]

    arXiv 1994

  15. [15]

    Djermoune and M

    E.-H. Djermoune and M. Tomczak, Statistical analysis of the Kumaresan-Tufts and Matrix Pencil methods in esti- mating a damped sinusoid, in2004 12th European Signal Processing Conference(IEEE, 2004) p. 73

    2004

  16. [16]

    Berti, V

    E. Berti, V. Cardoso, J. A. Gonz´ alez, and U. Sperhake, Mining information from binary black hole mergers: A comparison of estimation methods for complex exponen- tials in noise, Phys. Rev. D75, 124017 (2007), arXiv:gr- qc/0701086 [gr-qc]

    arXiv 2007

  17. [17]

    Berti, V

    E. Berti, V. Cardoso, and C. M. Will, Gravitational- wave spectroscopy of massive black holes with the space interferometer LISA, Phys. Rev. D73, 064030 (2006), arXiv:gr-qc/0512160 [gr-qc]

    Pith/arXiv arXiv 2006

  18. [18]

    Bhagwat, X

    S. Bhagwat, X. J. Forteza, P. Pani, and V. Ferrari, Ringdown overtones, black hole spectroscopy, and no- hair theorem tests, Phys. Rev. D101, 044033 (2020), arXiv:1910.08708 [gr-qc]

    Pith/arXiv arXiv 2020

  19. [19]

    Baibhav, M

    V. Baibhav, M. H.-Y. Cheung, E. Berti, V. Cardoso, G. Carullo, R. Cotesta, W. Del Pozzo, and F. Duque, Agnostic black hole spectroscopy: Quasinormal mode content of numerical relativity waveforms and limits of validity of linear perturbation theory, Phys. Rev. D108, 104020 (2023), arXiv:2302.03050 [gr-qc]

    Pith/arXiv arXiv 2023

  20. [20]

    Takahashi and H

    K. Takahashi and H. Motohashi, Iterative extraction of overtones from black hole ringdown, Classical and Quan- tum Gravity41, 195023 (2024), arXiv:2311.12762 [gr-qc]

    Pith/arXiv arXiv 2024

  21. [21]

    T. A. Clarke, M. Isi, P. D. Lasky, E. Thrane, M. Boyle, N. Deppe, L. E. Kidder, K. Mitman, J. Moxon, K. C. Nelli, W. Throwe, and N. L. Vu, Toward a self-consistent framework for measuring black hole ringdowns, Phys. Rev. D109, 124030 (2024), arXiv:2402.02819 [gr-qc]

    Pith/arXiv arXiv 2024

  22. [22]

    E. W. Leaver, Spectral decomposition of the perturbation response of the Schwarzschild geometry, Phys. Rev. D34, 384 (1986)

    1986

  23. [23]

    M. A. Scheel, M. Boyle, T. Chu, L. E. Kidder, K. D. Matthews, and H. P. Pfeiffer, High-accuracy waveforms for binary black hole inspiral, merger, and ringdown, Phys. Rev. D79, 024003 (2009), arXiv:0810.1767 [gr-qc]

    Pith/arXiv arXiv 2009

  24. [24]

    Echeverria, Gravitational-wave measurements of the mass and angular momentum of a black hole, Phys

    F. Echeverria, Gravitational-wave measurements of the mass and angular momentum of a black hole, Phys. Rev. D40, 3194 (1989)

    1989

  25. [25]

    L. S. Finn, Detection, measurement, and gravitational radiation, Phys. Rev. D46, 5236 (1992), arXiv:gr- qc/9209010 [gr-qc]

    arXiv 1992

  26. [26]

    H. Yang, K. Yagi, J. Blackman, L. Lehner, V. Paschalidis, F. Pretorius, and N. Yunes, Black Hole Spectroscopy with Coherent Mode Stacking, Phys. Rev. Lett.118, 161101 (2017), arXiv:1701.05808 [gr-qc]

    Pith/arXiv arXiv 2017

  27. [27]

    C. J. Moore, E. Finch, R. Buscicchio, and D. Gerosa, Testing general relativity with gravitational-wave cat- alogs: The insidious nature of waveform systematics, iScience24, 102577 (2021), arXiv:2103.16486 [gr-qc]

    Pith/arXiv arXiv 2021

  28. [28]

    Zhong, M

    H. Zhong, M. Isi, K. Chatziioannou, and W. M. Farr, Multidimensional hierarchical tests of general relativ- ity with gravitational waves, Phys. Rev. D110, 044053 (2024), arXiv:2405.19556 [gr-qc]

    Pith/arXiv arXiv 2024

  29. [29]

    Nakano, H

    H. Nakano, H. Takahashi, H. Tagoshi, and M. Sasaki, An effective search method for gravitational ringing of black holes, Phys. Rev. D68, 102003 (2003), arXiv:gr- qc/0306082 [gr-qc]

    arXiv 2003