Precision Ringdown Measurements of Binary Black Hole Remnants

Achal Kumar; Marek J. Szczepa\'nczyk; Poulami Dutta Roy; Sergey Klimenko

Precision Ringdown Measurements of Binary Black Hole Remnants

T0 review · 3 major / 5 minor · reviewed 2026-08-02 · deepseek-v4-flash

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

Pith's one-line read A denoised ringdown analysis of 16 binary black hole mergers finds the remnant's dominant mode frequency and damping time consistent with general relativity, with combined fractional deviations within ±2.8% (frequency) and about ±8% (dampin

desk verdict A genuinely new ringdown pipeline with careful injection work, but the headline precision rests on correction terms calibrated only in a narrow subspace; needs a direct comparison to prior constraints and a stated treatment of IMR circularity. read the letter →

arxiv 2607.13275 v2 pith:TCRWY7MX submitted 2026-07-14 gr-qc

Achal Kumar , Poulami Dutta Roy , Marek J. Szczepańczyk , Sergey Klimenko This is my paper · ORCID

classification gr-qc

keywords ringdownquasinormalmodesblackholespectroscopygravitationalwavesCoherentWaveBurstgeneralrelativitytestsGWTC-3binaryremnants

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish that the ringdown of a binary black hole remnant can be measured with higher precision than previous analyses by fitting quasinormal modes to denoised, reconstructed signals rather than raw strain data. Using a cumulative-energy reference to start the fit closer to merger, it constrains the dominant (2,2,0) mode's frequency and damping time for 16 events from the third gravitational-wave transient catalog. All 16 measurements are consistent with the predictions of general relativity, and the combined deviations are δf = 0.003 ± 0.028 and δτ = 0.050 (+0.081/−0.086) at 90% confidence. If correct, this demonstrates a path to sharper strong-field tests of gravity as detector sensitivity grows.

What carries the argument

The key mechanism is the denoised signal reconstruction by cWB—a coherent time-frequency pixel selection that estimates the signal without assuming a waveform model—followed by a two-mode damped-sinusoid fit for the (2,2,0) and (2,2,1) modes. The ringdown window is anchored not to the noisy signal peak but to the normalized cumulative energy e(t) = E(t)/E(T_end), set at 0.82, which reduces reference-time jitter. Correction factors calibrated on simulated injections account for reconstruction bias and allow the window to start at ~3.5 t_Mf after merger, with the systematic uncertainty budget validated through injection-based coverage studies.

What would settle it

Calibrate the same correction pipeline on an injection set that mirrors the catalog's parameter space—unequal masses, aligned spins up to χ≈0.9, and nonzero inclination—then check whether the recovered δf220 and δτ220 for the injected signals remain within the 90% intervals; any systematic offset exceeding the quoted uncertainty would show the result is not universal.

Watch

Extended reading notes

Core claim

The central claim is that for all 16 binary black hole mergers analyzed, the measured frequency and damping time of the dominant (2,2,0) quasi-normal mode of the remnant are consistent with general relativity. The combined fractional deviations are δf220 = 0.003 ± 0.028 and δτ220 = 0.050 (+0.081/−0.086) at 90% confidence, both consistent with zero. The paper further claims that its cWB-reconstruction-based method yields tighter constraints than previous ringdown analyses because the cumulative-energy reference time lets the fit start at approximately 3.5 t_Mf after the merger, and because the denoising reduces the impact of non-Gaussian detector noise.

Load-bearing premise

The correction that makes the early ringdown usable is calibrated only on equal-mass, non-spinning, face-on binary injections, yet is applied as a single ~5% shift to every real event, some of which have remnant spins up to 0.94 and unequal masses; if that mapping is off by a few percent, the quoted constraints shift by more than their intervals.

Editorial extensions

If this is right

If the claimed precision is real, GR tests from ringdown alone can reach the sub-3% level in frequency with a modest number of events.
The cumulative-energy reference could be adopted by other ringdown analyses to reduce systematic uncertainty in the start time.
The approach applies directly to future, more sensitive detector networks, where the same method should produce even tighter bounds.
The overtone amplitude estimates provide a way to check the validity of the assumed linear perturbation theory.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The correction functions are derived from a narrow injection set (equal mass, non-spinning, face-on), yet are applied universally; a test with spinning, unequal-mass injections would reveal whether the claimed precision holds for all events in the catalog.
The combined constraint on the remnant spin, χf = 0.709 (+0.050/−0.061), is inferred from the f·τ product; this offers a population-level test of formation scenarios, though it depends on the same correction universality.
One could extend this pipeline to subdominant modes such as (3,3) or (2,1) to attempt genuine black-hole spectroscopy; the denoising approach might make those modes accessible.
The frequency systematic uncertainty (inflated by 62%) suggests that the pixel selection partially cancels the Gaussian noise reduction; a refined pixel-selection rule could recover the full gain.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

You should know two things about this one. First, the method is a real advance: cWB denoising plus a cumulative-energy ringdown window that starts ~3.5 Mf after peak, with injection-calibrated corrections and coverage validation. That combination is new and the systematic-error work is unusually careful. Second, the central precision claim—"tighter constraints than previous measurements"—is not actually demonstrated against any published number, and the universality of the correction terms is the load-bearing assumption.

What the paper does well: the injection study is honest and thorough. They inject SEOBNRv4HM signals into O4 data, measure reconstruction biases, correct the frequency for SNR-dependent pixel loss, and validate the error budget with coverage plots (Fig. 4). The frequency uncertainty inflation of 62% and the damping-time inflation of 5% are backed by those studies. That is exactly the kind of systematic accounting you want in a ringdown paper. The GR-consistency result for 16 events is plausible, and the combined constraints are reported with asymmetric uncertainties and a joint spin estimate that avoids an obvious bias.

Now the soft spots, in proportion. The correction terms (Eqs. 4–7) are calibrated entirely on q=1, nonspinning, zero-inclination injections, and the paper declares them universal with negligible uncertainty. But the corrections are ~5% in both frequency and damping time at the chosen analysis point e(Tw)=0.82—larger than the quoted combined frequency uncertainty of ±2.8%. Real events include remnant spins up to 0.94 and overtone amplitudes up to 0.94, far outside the calibration subspace. The coverage studies cannot validate extrapolation to those parameters because the injections live in the same narrow subspace. If the true correction for a high-spin, unequal-mass event differs by 2–3% from the adopted zero-spin curve, the reported central values shift by more than the error bars. That is not a fatal flaw, but it is a genuine weakness in the headline precision claim, and the paper's own admission that corrections vary 4–6% over the mass range makes the "negligible" statement hard to accept.

Two smaller issues: the abstract claims tighter constraints without any direct comparison to previous ringdown measurements, and the IMR-predicted f220/tau220 used as the GR baseline are themselves uncertain—particularly for the high-mass events that dominate this sample, where the IMR fit is substantially informed by the same post-merger data being tested. The paper does not discuss that circularity or state how IMR prediction uncertainty is propagated into delta f and delta tau.

For a ringdown practitioner this paper deserves a serious referee, not a desk reject. The method is substantive and the core result is probably right, but the precision claim needs to be supported with a comparison table, broader calibration (or conservative error inflation), and a clear caveat on the IMR circularity. I would send it to review and ask for those revisions.

Referee Report

3 major / 5 minor

Summary. The paper presents a new ringdown analysis pipeline for binary black hole remnants. The method uses cWB-reconstructed waveforms rather than raw strain, defines the ringdown window via the normalized cumulative energy at e(t)=0.82 (~3.5 t_Mf after peak), fits a two-mode damped-sinusoid model with BILBY, and applies parameter-dependent correction functions calibrated on SEOBNRv4HM injections. The analysis is applied to 16 GWTC-3 events with SNR>12 and M>50 M_sun. The reported combined fractional deviations from GR are δf220 = 0.003^{+0.028}_{-0.028} and δτ220 = 0.050^{+0.081}_{-0.086}, consistent with GR, and the paper claims this yields tighter QNM constraints than previous measurements.

Significance. If the method is sound, it represents a useful step toward precision ringdown tests: the cumulative-energy reference time is more robust than the usual peak-time reference, the cWB denoising is a practical way to mitigate detector noise, and the systematic-error accounting via injection-based coverage studies (Fig. 4) is a strength. However, the central precision claim rests on the assertion that the correction functions are universal. The calibration set is narrow (q=1, nonspinning, zero inclination, SNR 12-50), while the real events include remnant spins up to χf≈0.94 and overtone amplitudes up to ε≈0.94. The corrections at the chosen window are ~5% in both frequency and damping time, which is larger than the quoted combined frequency uncertainty (±2.8%). The paper's own text reports 4-6% variation over the remnant-mass range and dependence on aligned spins, so the claimed 'negligible' uncertainty of the universal correction is load-bearing and not validated by the coverage studies. The claim of tighter constraints than previous measurements is also not supported by any explicit numerical comparison.

major comments (3)

[Secs. II B 3, III A; Eqs. (4)-(7), Fig. 3, Table I] The universal correction functions are calibrated exclusively on SEOBNRv4HM injections with q=1, nonspinning components, and zero inclination. Yet Table I includes events with remnant spin as high as χf≈0.94 and overtone amplitude ε≈0.94 (e.g., GW190521_074359). The correction at e=0.82 is approximately 5% in both frequency and damping time, larger than the quoted combined frequency uncertainty of ±2.8%. The paper itself states that the corrections vary by 4-6% over the mass range and depend on aligned spins producing rapidly rotating remnants. The coverage validation in Fig. 4 uses injections in the same q=1, nonspinning, face-on subspace and therefore cannot certify the extrapolation. A 2-3% miscalibration for a subset of events would shift δf220 and δτ220 by more than the quoted intervals. This is the central load-bearing assumption and must be supported by injection studies spanning
[Abstract and Sec. V] The abstract and conclusion claim the method 'yields tighter constraints on the QNM frequency and damping time than previous measurements' and specifically refer to the GWTC-4 ringdown analysis [47]. No numerical comparison to [47] (or to other published ringdown constraints, e.g., the frequency-domain SEOBNRv5PHM results or direct QNM fits) is provided anywhere in the manuscript. Without such a comparison, the central precision claim is not demonstrated. The authors should tabulate the published 90% intervals and show quantitatively how their combined δf220 and δτ220 improve on them, including the impact of any differences in event selection and prior choices.
[Eq. (11) and surrounding text] The fractional deviations δf220 and δτ220 are defined relative to f_IMR and τ_IMR, the values predicted by the IMR analysis. The IMR parameter estimation fits the full signal with GR waveforms that include a ringdown model; therefore the 'GR prediction' is not independent of the same post-merger data that the ringdown analysis uses. If a real deviation from GR were present in the ringdown, the IMR fit could partially absorb it into the inferred mass and spin, biasing δ toward zero. The manuscript does not quantify this effect or discuss how it affects the interpretation of the reported consistency with GR. This is particularly relevant because the claimed precision (e.g., ±2.8% in combined δf220) is smaller than the systematic effects being neglected.

minor comments (5)

[Eq. (1)] The exponentials are missing the time variable: the model should read h(t) = A exp[i(2π f0 t + φ0)] + ε A exp[i(2π f1 t + φ1)], or an equivalent convention should be stated.
[Eq. (2)] The quadrature H(t) is not defined in the text. Presumably it is the Hilbert transform of h(t), but this should be stated explicitly.
[Sec. II B 2] The 'three cycles' duration of the ringdown window is not precisely defined: is it three cycles of the fundamental mode, the dominant instantaneous frequency, or the waveform envelope? Please specify.
[Sec. III B] The statement that the systematic uncertainty in damping time is increased by only 5% while frequency is increased by 62% is clear, but the derivation of the total uncertainty (Table I caption) would benefit from an explicit formula linking σ_stat and σ_sys to the quoted 90% intervals.
[General] There are minor typos and formatting issues, e.g., 'SEOBNRv4HM' is not consistently typeset, and the reference list contains a preprint number (2603.19021) that should be updated if published.

Circularity Check

0 steps flagged · score 2.0 of 10

No construction-level circularity; the analysis is calibrated on GR simulations and tested against IMR predictions, though the 'universal' correction is extrapolated beyond its calibration subspace.

full rationale

The paper's central chain is: cWB reconstructs a denoised waveform; a cumulative-energy reference defines an early ringdown window; a correction calibrated on SEOBNRv4HM injections accounts for merger/bias; a two-mode damped-sinusoid fit returns (f, τ); deviations from IMR-based GR predictions are formed via Eq. (11). The correction terms (Eqs. 4–7) are fitted to simulated signals with known GR answers, so applying them to real events is calibration, not a self-referential prediction. The comparison baseline in Eq. (11) uses f_IMR, τ_IMR from a separate IMR parameter estimation; although the IMR waveform includes ringdown information, the ringdown fit is an independent estimator, so δf220 and δτ220 are not identically zero by construction. No equation in the paper reduces the measured quantity to its input. The self-citations for cWB (refs [40–43]) cite an established search/reconstruction algorithm and are not load-bearing in a circular sense. The main weakness is an extrapolation: the universal correction is calibrated only on q=1, non-spinning, zero-inclination injections, while reported events have high remnant spins (χf up to 0.94) and unequal masses. This is a systematic-uncertainty concern, not a circularity. The 5% correction magnitude exceeds the quoted combined frequency uncertainty, but this affects accuracy, not logical circularity. Score 2 reflects minor self-citation and the partially data-derived GR baseline, with no identified construction-level circular step.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central claim rests on a chain of calibration steps fitted to simulated GR injections plus standard GR assumptions, rather than on a first-principles derivation. The free parameters are the empirical constants defining the ringdown window and the cWB reconstruction-bias corrections; the axioms are the GR mapping between remnant (M, χ) and QNM parameters, the assumption that the merger-ringdown transition is captured by SEOBNRv4HM, the assumption that cWB reconstruction errors are correctable, and the use of IMR-based parameters as an independent GR baseline. No new physical entities are introduced.

free parameters (5)

Cumulative-energy window threshold = 0.82 (modified by (1 − N/(3·SNR²))·0.82 at low SNR)
Chosen by hand/calibration in Secs. II B 2-3 to balance merger contamination vs SNR; sets the window start at ~3.5 t_Mf after peak, the point where ~5% corrections are applied.
SNR-dependent frequency-bias correction = 1/(1 − e^{−SNR/5.9})
Empirical fit to cWB reconstruction bias in Sec. III A; at SNR=12 the correction is ~15% in frequency — the largest single correction applied.
Window-offset vs remnant-mass relation = T_w − T_p ≈ 0.028·M_f + 0.94 (in t_Mf)
Linear fit from injections (Eq. 7, Sec. III A) used to map the estimated M_f back to a cumulative-energy value for the corrections.
Systematic-error inflation factors = 62% (frequency), 5% (damping time) added in quadrature to statistical errors
Chosen in Sec. III B so that injection coverage matches nominal (Fig. 4); the 62% combines a 27% jitter inflation with a 100%-of-statistical pixel-selection systematic.
Error-model scalings from injections = σ_e≈129/SNR, σ_p≈389/SNR, σ_sys,f≈87/SNR %, σ_sys,τ≈197/SNR %, σ_stat,f≈110/SNR %, σ_stat,τ≈584/SNR %
Numerical coefficients (Eqs. 8-10) fitted to injection recovery, used to derive the jitter and systematic inflation factors.

assumptions (6)

domain assumption Kerr QNM relations: f220 and τ220 are uniquely determined by remnant (M, a) via GR perturbation theory
Used to build the two-mode model and to map measurements to Mf and χf (Sec. II B 1); standard GR result invoked via refs [32-35].
ad hoc to paper The linear QNM regime is already reached at e(t)=0.82 (~3.5 t_Mf after peak), with remaining merger contamination removable by GR-calibrated corrections
The whole 'move closer to the merger' strategy rests on corrections fit to SEOBNRv4HM (Secs. II B 3, III); this is not established from first principles.
domain assumption cWB reconstruction faithfully recovers the ringdown signal, with TF pixel-selection biases that are small and empirically correctable
Analysis is performed on reconstructed waveforms rather than raw strain; biases are handled by SNR-dependent corrections (Secs. II A, III A).
domain assumption Noise in O4 injection data is representative of the noise in the GWTC-3 events analyzed (O1/O2/O3a/O3b)
Systematic errors and coverage are calibrated with injections into O4 data (Sec. III), then applied to events from earlier observing runs.
domain assumption IMR parameter estimates (f_IMR, τ_IMR, M, χ) are unbiased GR predictions suitable as the comparison baseline
Used in Eq. (11) and in event selection (M>50 M⊙); for high-mass events the IMR fit is dominated by merger-ringdown data, making the baseline partially non-independent (Sec. IV).
ad hoc to paper The correction curves depend only weakly on source parameters beyond remnant mass, so a universal calibration applies
Eqs. (6)-(7) assume universal SNR- and mass-dependent window corrections across the event population; the paper asserts rather than bounds the exceptions (aligned spins, unequal masses, waveform approximant).

how reviews work

0 comments

Cite this review

Pith. "Pith review of Precision Ringdown Measurements of Binary Black Hole Remnants." pith.science (2026). https://pith.science/paper/TCRWY7MX

@misc{pith2026260713275,
  author       = {Pith},
  title        = {Pith review of: Precision Ringdown Measurements of Binary Black Hole Remnants},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TCRWY7MX}},
  note         = {Machine review of arXiv:2607.13275}
}

abstract

The ringdown gravitational wave from a binary black hole (BBH) merger is a superposition of quasi-normal modes (QNMs) of the remnant black hole. In general relativity (GR), QNMs are damped harmonic oscillations with frequencies and damping times uniquely determined by the remnant's mass and spin. The measurement of the ringdown modes and performing black hole spectroscopy provides a tool to test the validity of GR. In this work, we present RingCWB, a ringdown analysis method based on coherent WaveBurst (cWB), an unmodeled pipeline for the detection and reconstruction of gravitational-wave signals. This method yields tighter constraints on the QNM frequency and damping time than previous measurements. The improved precision results from the noise reduction achieved by the cWB reconstruction and the enhanced ringdown analysis, which probes the remnant properties at earlier times, closer to the merger. We have analysed publicly available binary black hole (BBH) detections from the third Gravitational-Wave Transient Catalog (GWTC-3). For all events considered, the measured frequency and damping time of the dominant $(l,m)=(2,2)$ mode are found to be consistent with the predictions of GR. A combined analysis further strengthens these constraints, yielding fractional deviations in frequency $\delta f_{220} = -0.005_{-0.028}^{+0.028}$ and damping time $\delta\tau_{220} = 0.032_{-0.090}^{+0.108}$, consistent with zero within the quoted uncertainties.

Figures

Figures reproduced from arXiv: 2607.13275 by the authors.

**Figure 2.** FIG. 2: The reconstructed time-domain signals (left [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** FIG. 3: Correction terms for the (2 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: FIG. 4: Coverage of the marginalized posteriors for the [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]

Discussion (0). Continue with ORCID to comment.

Pith tools

Reviewed August 2, 2026 · model on record in the stance chip above.