pith. sign in

arxiv: 2512.15877 · v2 · pith:2SEY32WPnew · submitted 2025-12-17 · 🌀 gr-qc · astro-ph.HE· hep-ph

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

classification 🌀 gr-qc astro-ph.HEhep-ph
keywords greybody factorringdowngravitational wavesblack hole mergersnumerical relativityfrequency domainSXS catalog
0
0 comments X

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.

The paper shows that the greybody factor of the remnant black hole modulates the post-merger ringdown signal in binary coalescences. A simple four-parameter model built on this factor matches the frequency-domain amplitude of many numerical relativity waveforms from the SXS catalog for comparable-mass aligned-spin systems. The approach reaches mismatches of order 10 to the minus five and improves on prior models by about two orders of magnitude. Analytical fits are given for the model parameters in terms of the initial masses and aligned spins. This setup supports new consistency tests of the ringdown phase that complement traditional black hole spectroscopy.

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

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

  • 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

Figures reproduced from arXiv: 2512.15877 by Emanuele Berti, Paolo Pani, Romeo Felice Rosato, Sophia Yi.

Figure 1
Figure 1. Figure 1: FIG. 1. Mismatch [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Illustration of the fitting procedure for the simulation [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Cross-validation of the polynomial model for the fitted [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: , where the dependence reduces to two dimensions and is therefore easier to visualize. The same procedure will then be applied in full generality to the complete three-dimensional parameter space (δ, χ+, χ−). We consider the training set and validation set defined above. We then compute the CV MSE for polynomials of increasing total degree N , and identify the optimal degree as the value that minimizes the… view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison between the values predicted by the [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Distribution of mismatches between the numerical [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Relative “importance” of different BBH parameters, [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Comparison between the FFT of the preprocessed [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Predictions for [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Mean mismatch for the [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

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)
  1. 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.
  2. 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

2 responses · 0 unresolved

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
  1. 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

  2. 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

1 steps flagged

Four-parameter greybody ringdown model achieves reported accuracy via per-waveform fitting to NR data followed by secondary regression on progenitor parameters

specific steps
  1. 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

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the domain assumption that greybody factors modulate ringdown amplitudes in a way capturable by four fitted parameters, plus the practical assumption that SXS waveforms are sufficiently accurate ground truth; no new entities are introduced.

free parameters (1)
  • four model parameters
    The model is explicitly described as four-parameter; these are fitted to the numerical waveforms and later given analytical expressions.
axioms (1)
  • domain assumption Greybody factor of the remnant black hole modulates the post-merger ringdown signal
    Invoked in the opening sentence as a recently shown fact that the present model builds upon.

pith-pipeline@v0.9.0 · 5433 in / 1426 out tokens · 22799 ms · 2026-05-16T21:17:52.970098+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Novel ringdown tests of general relativity with black hole greybody factors

    gr-qc 2026-04 unverdicted novelty 7.0

    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.

  2. Quasinormal mode/grey-body factor correspondence for Kerr black holes

    gr-qc 2025-12 conditional novelty 5.0

    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

53 extracted references · 53 canonical work pages · cited by 2 Pith papers · 16 internal anchors

  1. [1]

    The top-left panel of Fig

    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. [2]

    Using these values, we compute the corresponding reflectivity Rℓm(¯ω,χ)and construct the reduced spectrum Yℓm(ω) = Hℓm(ω)√ Rℓm(¯ω,χ)

    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. [3]

    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]

    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. [4]

    Within the validity region of the model, the values ofAℓmandpℓmobtained from different¯ωi are expected to be consistent within their uncertainties

    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. [5]

    importance

    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. [6]

    As clearly visible in Fig

    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. [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. [8]

    This prescription safely excludes the region where numerical noise affects the Fourier transform across theentirecatalog

    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. [9]

    Black hole spectroscopy: from theory to experiment

    E. Bertiet al., (2025), arXiv:2505.23895 [gr-qc]

  10. [10]

    A. G. Abacet al.(LIGO Scientific, Virgo, KAGRA), Phys. Rev. Lett.135, 111403 (2025), arXiv:2509.08054 [gr-qc]

  11. [11]

    (2025), arXiv:2509.08099 [gr-qc]

  12. [12]

    Crescimbeni, G

    F. Crescimbeni, G. Carullo, E. Berti, G. Caneva Santoro, M. H.-Y. Cheung, and P. Pani, (2025), arXiv:2511.02915 [gr-qc]

  13. [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]

  14. [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]

  15. [15]

    Cotesta, G

    R. Cotesta, G. Carullo, E. Berti, and V. Cardoso, Phys. Rev. Lett.129, 111102 (2022), arXiv:2201.00822 [gr-qc]

  16. [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. [17]

    S. Ma, L. Sun, and Y. Chen, Phys. Rev. D107, 084010 (2023), arXiv:2301.06639 [gr-qc]

  18. [18]

    Isi and W

    M. Isi and W. M. Farr, Phys. Rev. Lett.131, 169001 (2023), arXiv:2310.13869 [astro-ph.HE]

  19. [19]

    Carullo, R

    G. Carullo, R. Cotesta, E. Berti, and V. Cardoso, Phys. Rev. Lett.131, 169002 (2023), arXiv:2310.20625 [gr-qc]

  20. [20]

    C. D. Capano, M. Cabero, J. Westerweck, J. Abedi, S. Kastha, A. H. Nitz, Y.-F. Wang, A. B. Nielsen, and B. Krishnan, Phys. Rev. Lett.131, 221402 (2023), arXiv:2105.05238 [gr-qc]

  21. [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. [22]

    Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog

    R. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. D 103, 122002 (2021), arXiv:2010.14529 [gr-qc]

  23. [23]

    Tests of General Relativity with GWTC-3

    R. Abbottet al.(LIGO Scientific, VIRGO, KAGRA), (2021), arXiv:2112.06861 [gr-qc]

  24. [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. [25]

    Oshita, (2023), arXiv:2309.05725 [gr-qc]

    N. Oshita, Phys. Rev. D109, 104028 (2024), arXiv:2309.05725 [gr-qc]

  26. [26]

    R. F. Rosato, K. Destounis, and P. Pani, Phys. Rev. D 110, L121501 (2024), arXiv:2406.01692 [gr-qc]

  27. [27]

    R. F. Rosato, S. Biswas, S. Chakraborty, and P. Pani, Phys. Rev. D111, 084051 (2025), arXiv:2501.16433 [gr- qc]

  28. [28]

    Okabayashi and N

    K. Okabayashi and N. Oshita, Phys. Rev. D110, 064086 (2024), arXiv:2403.17487 [gr-qc]

  29. [29]

    Oshita, K

    N. Oshita, K. Takahashi, and S. Mukohyama, Phys. Rev. D110, 084070 (2024), arXiv:2406.04525 [gr-qc]

  30. [30]

    R. F. Rosato, S. Biswas, S. Chakraborty, and P. Pani, (2025), arXiv:2511.08692 [gr-qc]

  31. [31]

    Nollert, Phys

    H.-P. Nollert, Phys. Rev. D53, 4397 (1996), arXiv:gr- qc/9602032

  32. [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]

  33. [33]

    J. L. Jaramillo, R. Panosso Macedo, and L. Al Sheikh, Phys. Rev. X11, 031003 (2021), arXiv:2004.06434 [gr-qc]

  34. [34]

    M. H.-Y. Cheung, K. Destounis, R. P. Macedo, E. Berti, and V. Cardoso, Phys. Rev. Lett.128, 111103 (2022), arXiv:2111.05415 [gr-qc]

  35. [35]

    Berti, V

    E. Berti, V. Cardoso, M. H.-Y. Cheung, F. Di Filippo, F. Duque, P. Martens, and S. Mukohyama, Phys. Rev. D106, 084011 (2022), arXiv:2205.08547 [gr-qc]

  36. [36]

    Boyle et al

    M. Boyleet al., Class. Quant. Grav.36, 195006 (2019), arXiv:1904.04831 [gr-qc]

  37. [37]

    M. A. Scheelet al., Class. Quant. Grav.42, 195017 (2025), arXiv:2505.13378 [gr-qc]

  38. [38]

    S. A. Teukolsky, Phys. Rev. Lett.29, 1114 (1972)

  39. [39]

    W. H. Press and S. A. Teukolsky, Nature238, 211 (1972)

  40. [40]

    A. A. Starobinskii, Sov. Phys. JETP37, 28 (1973)

  41. [41]

    A. A. Starobinskil and S. M. Churilov, Sov. Phys. JETP 65, 1 (1974)

  42. [42]

    L. S. Finn and D. F. Chernoff, Phys. Rev. D47, 2198 (1993), arXiv:gr-qc/9301003

  43. [43]
  44. [44]

    Blanchet, T

    L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. Iyer, Phys. Rev. Lett.93, 091101 (2004), arXiv:gr- qc/0406012

  45. [45]

    M. H.-Y. Cheung, E. Berti, V. Baibhav, and R. Cotesta, Phys. Rev. D109, 044069 (2024), [Erratum: Phys.Rev.D 110, 049902 (2024), Erratum: Phys.Rev.D 112, 049901 (2025)], arXiv:2310.04489 [gr-qc]

  46. [46]

    Hastie, R

    T. Hastie, R. Tibshirani, and J. Friedman,The elements of statistical learning: data mining, inference and predic- tion, 2nd ed. (Springer, 2009)

  47. [47]

    Anselmo, C

    C. Anselmo, C. Pacilio, and D. Gerosa, (2025), arXiv:2512.05193 [gr-qc]

  48. [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]

  49. [49]

    Chenet al., Phys

    Y. Chenet al., Phys. Rev. D110, 064049 (2024), arXiv:2405.06197 [gr-qc]

  50. [50]

    Advanced Methods in Black-Hole Perturbation Theory

    P. Pani, Int. J. Mod. Phys. A28, 1340018 (2013), arXiv:1305.6759 [gr-qc]

  51. [51]

    W. H. Press and S. A. Teukolsky, Astrophys. J.185, 649 (1973)

  52. [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

  53. [53]

    Pedregosa, G

    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)