Modeling the frequency-domain ringdown amplitude of comparable-mass mergers with greybody factors
Pith reviewed 2026-05-16 21:17 UTC · model grok-4.3
The pith
A four-parameter greybody model reproduces the frequency-domain ringdown amplitude of comparable-mass aligned-spin mergers with mismatches around 10^{-5}.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The greybody factor of the remnant black hole modulates the post-merger ringdown signal, and a simple four-parameter model based on it accurately reproduces the frequency-domain amplitude across a large set of comparable-mass aligned-spin numerical relativity waveforms. The model achieves mismatches of order 10^{-5} on SXS catalog data, improves existing models by roughly two orders of magnitude, supplies analytical fits for its parameters in terms of progenitor masses and spins, and opens the door to new ringdown consistency tests.
What carries the argument
The greybody factor of the remnant black hole, which modulates the post-merger ringdown signal and serves as the basis for a four-parameter frequency-domain amplitude model.
If this is right
- The model supplies ready-to-use analytical expressions for ringdown amplitudes in terms of initial binary parameters.
- It enables direct consistency tests between the ringdown phase and the remnant black-hole properties.
- The approach improves waveform modeling accuracy by two orders of magnitude for the targeted class of mergers.
- It provides an alternative route to black-hole spectroscopy that focuses on amplitude rather than frequency alone.
Where Pith is reading between the lines
- The same greybody-based construction could be tested on precessing or unequal-mass systems to check whether the four-parameter form remains sufficient.
- If the model holds more broadly, it could tighten constraints on deviations from general relativity in the ringdown regime of gravitational-wave events.
- Combining the amplitude fits with existing quasinormal-mode frequency measurements might yield joint tests of both amplitude and frequency content.
Load-bearing premise
The greybody factor of the remnant black hole modulates the post-merger ringdown amplitude in a way that a simple four-parameter model can capture for comparable-mass aligned-spin mergers.
What would settle it
Applying the four-parameter model to new numerical relativity waveforms outside the SXS comparable-mass aligned-spin range and finding mismatches well above 10^{-5} would show the model does not generalize as claimed.
Figures
read the original abstract
It was recently shown that, in a binary coalescence, the greybody factor of the remnant black hole modulates the post-merger ringdown signal. In this work, we demonstrate that a simple four-parameter model based on the greybody factor accurately reproduces the frequency-domain amplitude of a large set of comparable-mass, aligned-spin numerical relativity waveforms from the SXS catalog, achieving mismatches of order ${\cal O}(10^{-5})$ and improving existing models by roughly two orders of magnitude. We also identify the optimal initial frequency for applying the model in the frequency domain and provide analytical fits of the model parameters in terms of the progenitor masses and aligned spins. Our results pave the way for new consistency tests of the ringdown phase, complementary to traditional black hole spectroscopy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that a simple four-parameter model based on the greybody factor of the remnant black hole accurately reproduces the frequency-domain amplitude of the ringdown signal for comparable-mass, aligned-spin binary black hole mergers. Using a large set of SXS numerical-relativity waveforms, the model achieves mismatches of order O(10^{-5}), improves on existing models by roughly two orders of magnitude, identifies an optimal initial frequency for application in the frequency domain, and supplies analytical fits for the four parameters in terms of progenitor masses and aligned spins.
Significance. If the empirical accuracy holds under full scrutiny, the work provides a compact, analytically tractable model for ringdown amplitudes that could support new consistency tests of general relativity complementary to black-hole spectroscopy. The use of a broad SXS catalog, the reported mismatch improvement, and the regression of parameters onto progenitor quantities are concrete strengths that aid reproducibility.
major comments (2)
- The central claim that the greybody factor modulates the ringdown amplitude rests on an empirical four-parameter fit rather than a derivation from the transmission probability through the effective potential. No section derives the functional form from first principles; parameters are fitted per waveform and then regressed, leaving open whether the O(10^{-5}) mismatches reflect a genuine physical mechanism or a flexible phenomenological ansatz tuned to the catalog.
- Details on data-selection criteria for the SXS set, error-bar treatment, and full validation plots (including dependence on start time and overtone content) are required to substantiate the cross-catalog accuracy claim; the abstract-only presentation limits verification of the load-bearing result.
Simulated Author's Rebuttal
We thank the referee for their positive assessment and recommendation for minor revision. We address each major comment below and will update the manuscript accordingly to improve clarity and reproducibility.
read point-by-point responses
-
Referee: The central claim that the greybody factor modulates the ringdown amplitude rests on an empirical four-parameter fit rather than a derivation from the transmission probability through the effective potential. No section derives the functional form from first principles; parameters are fitted per waveform and then regressed, leaving open whether the O(10^{-5}) mismatches reflect a genuine physical mechanism or a flexible phenomenological ansatz tuned to the catalog.
Authors: We agree that the four-parameter model is phenomenological: its functional form is chosen to reproduce the shape of the greybody transmission probability, but the parameters themselves are determined by fitting to numerical-relativity waveforms and then regressed onto progenitor quantities. This construction is motivated by the prior result (cited in the introduction) that the remnant greybody factor modulates the ringdown amplitude, yet we do not claim a first-principles derivation from the effective potential. The low mismatches across the SXS catalog provide empirical support for the physical relevance of the ansatz, but we acknowledge it remains a tuned model. In the revised manuscript we will expand the discussion in Section II to explicitly state the phenomenological character, the motivation from the greybody factor, and the limitations of the approach. revision: partial
-
Referee: Details on data-selection criteria for the SXS set, error-bar treatment, and full validation plots (including dependence on start time and overtone content) are required to substantiate the cross-catalog accuracy claim; the abstract-only presentation limits verification of the load-bearing result.
Authors: We will add the requested information. The revised manuscript will include a new subsection in Section III that specifies the SXS waveform selection criteria (mass-ratio and spin ranges, resolution cuts, and exclusion of cases with significant eccentricity), describes how numerical truncation errors are estimated and propagated into the mismatch calculation, and presents additional validation figures showing mismatch versus start time and versus the number of overtones retained. These plots and the associated text will be placed in the main body and an appendix to allow full verification of the reported accuracy. revision: yes
Circularity Check
Four-parameter greybody ringdown model achieves reported accuracy via per-waveform fitting to NR data followed by secondary regression on progenitor parameters
specific steps
-
fitted input called prediction
[Abstract]
"we demonstrate that a simple four-parameter model based on the greybody factor accurately reproduces the frequency-domain amplitude of a large set of comparable-mass, aligned-spin numerical relativity waveforms from the SXS catalog, achieving mismatches of order O(10^{-5}) and improving existing models by roughly two orders of magnitude. [...] We also identify the optimal initial frequency for applying the model in the frequency domain and provide analytical fits of the model parameters in terms of the progenitor masses and aligned spins."
The four parameters are fitted to each individual NR waveform to produce the quoted low mismatches; the 'analytical fits' are then obtained by regressing those fitted parameter values against progenitor masses and spins. The reproduction accuracy on the catalog is therefore achieved by construction of the per-waveform fit rather than by an independent derivation from the greybody transmission probability.
full rationale
The paper's central demonstration is that a four-parameter model based on the remnant greybody factor reproduces the frequency-domain ringdown amplitude of SXS waveforms to O(10^{-5}) mismatches. This accuracy is obtained by fitting the four parameters directly to each numerical-relativity waveform in the catalog; the subsequent provision of analytical expressions for those parameters in terms of progenitor masses and spins is a secondary fit to the already-fitted values. The initial premise that the greybody factor modulates the ringdown signal is adopted from prior work without independent derivation of the specific functional form in this manuscript. Consequently the low mismatches validate the flexibility of the chosen ansatz on the fitted data rather than constituting an a-priori prediction, producing partial circularity of the fitted-input-called-prediction type.
Axiom & Free-Parameter Ledger
free parameters (1)
- four model parameters
axioms (1)
- domain assumption Greybody factor of the remnant black hole modulates the post-merger ringdown signal
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
a simple four-parameter model based on the greybody factor accurately reproduces the frequency-domain amplitude... H(2)ℓm(ω)=MAℓm√Rℓm(ω¯,χ)/ω¯pℓm
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the reflectivity Rℓm(ω¯,χ) encodes the reflection probability... computed by integrating the Teukolsky equation
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 2 Pith papers
-
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.
-
Quasinormal mode/grey-body factor correspondence for Kerr black holes
WKB analysis of the Teukolsky equation establishes a quasinormal-mode to greybody-factor correspondence for Kerr black holes that holds in the eikonal limit for gravitational perturbations and matches numerics at high...
Reference graph
Works this paper leans on
-
[1]
We compute the Fourier transform of one of the numerical relativity simulations from the SXS cata- log, as described in Appendix A, and we discard it above the cutoff frequency¯ωcut. The top-left panel of Fig. 2 shows the resulting spectrum for the sim- ulation SXS:BBH:3982, focusing on the dominant (ℓ,m) = (2,2)mode
-
[2]
The remnant mass M and dimensionless spin χ are available in the SXS metadata. Using these values, we compute the corresponding reflectivity Rℓm(¯ω,χ)and construct the reduced spectrum Yℓm(ω) = Hℓm(ω)√ Rℓm(¯ω,χ)
-
[3]
We fitYℓm(ω)with our model MA ℓm/¯ωpℓm. Since the minimum frequency of validity of the model is not known a priori, we consider a set of starting frequencies ¯ωi such that ¯ωi∈[0,¯ωcut]. The lowest frequency shown in Fig. 2 is¯ω= 0.2, which still lies within the inspiral region for the example simulation SXS:BBH:3982. In order to improve the stability of ...
-
[4]
We then study the behavior of the fitted parameters Aℓmand pℓmas functions of the starting frequency ¯ωi. Within the validity region of the model, the values ofAℓmandpℓmobtained from different¯ωi are expected to be consistent within their uncertainties. Outside this region,Aℓmand pℓmvary significantly with ¯ωi, as the model no longer reproduces the data a...
-
[5]
Finally, we compare the fitted model with the data. In the top-left panel of Fig. 2, the dashed pink line shows the best-fit model for the simulation SXS:BBH:3982, while the solid black line corre- sponds to the numerical spectrum. The bottom-left panel shows the mismatch, computed as in Eq.(4), as a function of the starting frequency¯ωi. As ex- pected, t...
-
[6]
How to apply the hyper-fits When applying the model to numerical data, it is im- portant to select the frequency window in which the power-law behavior of the spectrum is reliably captured and the data are not compromised by the noise. As clearly visible in Fig. 5 (and consistently across all simulations), the magnitude of the Fourier transformHℓm(ω)exhib...
-
[7]
A safe, universal choice is ¯ωin ≳0.75 ¯ω∗, which ensures that the model is applied well within the power-law regime for the vast majority of simu- lations
-
[8]
The upper frequency¯ωf can be defined as the fre- quency at which the magnitude of the Fourier trans- form decreases to one twentieth of its value at¯ω∗. This prescription safely excludes the region where numerical noise affects the Fourier transform across theentirecatalog. However, forsimulationsinwhich the Fourier-transform amplitude is intrinsically v...
-
[9]
Black hole spectroscopy: from theory to experiment
E. Bertiet al., (2025), arXiv:2505.23895 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[10]
A. G. Abacet al.(LIGO Scientific, Virgo, KAGRA), Phys. Rev. Lett.135, 111403 (2025), arXiv:2509.08054 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[11]
(2025), arXiv:2509.08099 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[12]
F. Crescimbeni, G. Carullo, E. Berti, G. Caneva Santoro, M. H.-Y. Cheung, and P. Pani, (2025), arXiv:2511.02915 [gr-qc]
-
[13]
Observational Black Hole Spectroscopy: A time-domain multimode analysis of GW150914
G. Carullo, W. Del Pozzo, and J. Veitch, Phys. Rev. D 99, 123029 (2019), [Erratum: Phys.Rev.D 100, 089903 (2019)], arXiv:1902.07527 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[14]
M. Isi, M. Giesler, W. M. Farr, M. A. Scheel, and S. A. Teukolsky, Phys. Rev. Lett.123, 111102 (2019), arXiv:1905.00869 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[15]
R. Cotesta, G. Carullo, E. Berti, and V. Cardoso, Phys. Rev. Lett.129, 111102 (2022), arXiv:2201.00822 [gr-qc]
-
[16]
Searching for a ringdown overtone in GW150914,
E. Finch and C. J. Moore, Phys. Rev. D106, 043005 (2022), arXiv:2205.07809 [gr-qc]
- [17]
- [18]
-
[19]
G. Carullo, R. Cotesta, E. Berti, and V. Cardoso, Phys. Rev. Lett.131, 169002 (2023), arXiv:2310.20625 [gr-qc]
- [20]
-
[21]
Ringdown of GW190521: Hints of multiple quasinormal modes with a precessional interpretation,
H. Siegel, M. Isi, and W. M. Farr, Phys. Rev. D108, 064008 (2023), arXiv:2307.11975 [gr-qc]
-
[22]
R. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. D 103, 122002 (2021), arXiv:2010.14529 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[23]
Tests of General Relativity with GWTC-3
R. Abbottet al.(LIGO Scientific, VIRGO, KAGRA), (2021), arXiv:2112.06861 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[24]
Guptaet al., (2024), 10.21468/SciPostPhysComm- Rep.5, arXiv:2405.02197 [gr-qc]
A. Guptaet al., (2024), 10.21468/SciPostPhysComm- Rep.5, arXiv:2405.02197 [gr-qc]
-
[25]
Oshita, (2023), arXiv:2309.05725 [gr-qc]
N. Oshita, Phys. Rev. D109, 104028 (2024), arXiv:2309.05725 [gr-qc]
- [26]
-
[27]
R. F. Rosato, S. Biswas, S. Chakraborty, and P. Pani, Phys. Rev. D111, 084051 (2025), arXiv:2501.16433 [gr- qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[28]
K. Okabayashi and N. Oshita, Phys. Rev. D110, 064086 (2024), arXiv:2403.17487 [gr-qc]
- [29]
-
[30]
R. F. Rosato, S. Biswas, S. Chakraborty, and P. Pani, (2025), arXiv:2511.08692 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [31]
-
[32]
Can environmental effects spoil precision gravitational-wave astrophysics?
E. Barausse, V. Cardoso, and P. Pani, Phys. Rev. D89, 104059 (2014), arXiv:1404.7149 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[33]
J. L. Jaramillo, R. Panosso Macedo, and L. Al Sheikh, Phys. Rev. X11, 031003 (2021), arXiv:2004.06434 [gr-qc]
work page internal anchor Pith review arXiv 2021
- [34]
- [35]
-
[36]
M. Boyleet al., Class. Quant. Grav.36, 195006 (2019), arXiv:1904.04831 [gr-qc]
- [37]
-
[38]
S. A. Teukolsky, Phys. Rev. Lett.29, 1114 (1972)
work page 1972
-
[39]
W. H. Press and S. A. Teukolsky, Nature238, 211 (1972)
work page 1972
-
[40]
A. A. Starobinskii, Sov. Phys. JETP37, 28 (1973)
work page 1973
-
[41]
A. A. Starobinskil and S. M. Churilov, Sov. Phys. JETP 65, 1 (1974)
work page 1974
-
[42]
L. S. Finn and D. F. Chernoff, Phys. Rev. D47, 2198 (1993), arXiv:gr-qc/9301003
work page internal anchor Pith review Pith/arXiv arXiv 1993
-
[43]
C. Cutler and E. E. Flanagan, Phys. Rev. D49, 2658 (1994), arXiv:gr-qc/9402014
work page internal anchor Pith review Pith/arXiv arXiv 1994
-
[44]
L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. Iyer, Phys. Rev. Lett.93, 091101 (2004), arXiv:gr- qc/0406012
- [45]
- [46]
-
[47]
C. Anselmo, C. Pacilio, and D. Gerosa, (2025), arXiv:2512.05193 [gr-qc]
-
[48]
Mixing of spherical and spheroidal modes in perturbed Kerr black holes
E. Berti and A. Klein, Phys. Rev. D90, 064012 (2014), arXiv:1408.1860 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[49]
Y. Chenet al., Phys. Rev. D110, 064049 (2024), arXiv:2405.06197 [gr-qc]
-
[50]
Advanced Methods in Black-Hole Perturbation Theory
P. Pani, Int. J. Mod. Phys. A28, 1340018 (2013), arXiv:1305.6759 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[51]
W. H. Press and S. A. Teukolsky, Astrophys. J.185, 649 (1973)
work page 1973
-
[52]
Eigenvalues and eigenfunctions of spin-weighted spheroidal harmonics in four and higher dimensions
E. Berti, V. Cardoso, and M. Casals, Phys. Rev. D73, 024013 (2006), [Erratum: Phys.Rev.D 73, 109902 (2006)], arXiv:gr-qc/0511111
work page internal anchor Pith review Pith/arXiv arXiv 2006
-
[53]
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cour- napeau, M. Brucher, M. Perrot, and E. Duchesnay, J. Mach. Learn. Res.12, 2825 (2011)
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.