REVIEW 2 major objections 5 minor 79 references
Analytic waveform derivatives make Fisher-matrix tests of gravity stable, fast, and able to map how well space- and ground-based detectors constrain non-GR effects by PN order.
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 →
T0 review · grok-4.5
2026-07-11 20:42 UTC pith:SMQGRHC6
load-bearing objection Solid, usable analytic Fisher pipeline for ppE tests; the trends and complementarity claims hold under the stated assumptions. the 2 major comments →
Tests of general relativity using analytic derivatives of parametrized post-Einsteinian gravitational waveforms within the Fisher-matrix framework
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Fully analytic expressions for the derivatives of frequency-domain TaylorF2 waveforms with respect to both binary parameters and parametrized post-Einsteinian deformation parameters make Fisher-matrix calculations stable and efficient, and those calculations reveal systematic trends in the constraints on non-GR effects as functions of post-Newtonian order, detector type, and source population, together with clear complementarity between space- and ground-based detectors for low-frequency-accumulating effects.
What carries the argument
The amplitude-phase decomposition of the waveform derivative, ∂_θ h̃ = h̃ (X_θ + i Y_θ), with all X_θ and Y_θ written in closed form for the GR phase coefficients and for every ppE correction function Ξ and Υ; this identity turns every Fisher integral into a real-valued frequency integral free of finite-difference step-size choices.
Load-bearing premise
The calculation rests on a 2PN-truncated GR baseline and on the claim that neglected eccentricity and spin-precession effects stay smaller than the statistical errors for the parameters being studied; if those systematics are larger, the reported uncertainties and complementarity statements no longer hold.
What would settle it
Recompute a representative set of Fisher matrices with a higher-order or precessing/eccentric waveform model for the same detector bands and populations; if the recovered σ_δ values and the space-versus-ground ranking shift by more than the statistical width claimed in the paper, the central forecasts fail.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives fully analytic waveform derivatives for the frequency-domain TaylorF2 inspiral within the parametrized post-Einsteinian (ppE) framework and uses them to construct Fisher-matrix forecasts for non-GR effects. After writing the waveform as amplitude and phase factors (Eqs. 9–10), the authors obtain closed-form expressions for the logarithmic amplitude and phase derivatives X_θ and Y_θ with respect to the full parameter set (Eqs. 16–27 and Appendices A–B). These derivatives are applied to space-based (LISA, Taiji, TianQin), ground-based (LVK, ET) and multiband configurations, combined with several BBH population models, to map how constraints on a deformation parameter δ scale with PN order, chirp mass, observation time and detector band. The same machinery is then specialized to a suite of concrete non-GR and environmental effects (noncommutative gravity, dCS, EdGB, scalar–tensor, charge, graviton mass, Ġ, Ṁ, dynamical friction). The central claims are that the analytic approach is numerically stable and substantially faster than finite differences (Fig. 1), that the resulting constraints exhibit systematic PN-order and detector-type trends (Figs. 3–5), and that space- and ground-based detectors are complementary for low-frequency-accumulating corrections (Fig. 6 and Sec. 6).
Significance. If the forecasts hold under the stated assumptions, the paper supplies a transparent, reusable analytic toolkit that removes step-size systematics from large-scale ppE Fisher studies and makes the PN-order and detector-complementarity trends quantitatively accessible. The explicit derivative catalogue in Appendices A–B, the runtime comparison in Fig. 1, and the population-averaged constraints in Fig. 6 are concrete, reusable contributions that future forecast and MCMC pipelines can adopt directly. The work therefore has clear methodological value for precision tests of gravity with next-generation detector networks, even though the underlying Fisher and 2PN approximations are standard.
major comments (2)
- Sec. 2 (paragraphs after Eq. 4) asserts that systematic biases from neglecting eccentricity and spin-precession remain sub-dominant to statistical Fisher errors for the parameters considered, citing only the authors’ earlier work [43]. Because the reported complementarity statements and the layered PN-order trends rest on the absolute scale of σ_δ, a short quantitative check (or an explicit statement of the mass/spin/eccentricity regime in which the claim holds) is needed before those trends can be regarded as robust forecasts rather than purely statistical scalings.
- The GR baseline is truncated at 2PN while ppE corrections are varied from –5.5PN to +2PN (Sec. 2 and Sec. 5.1). For positive-PN corrections the missing higher-order GR phase terms can correlate with the ppE parameters; the paper should either recompute a representative subset of the Fig. 4/6 constraints with a 3.5PN (or higher) baseline or demonstrate that the relative detector rankings and the 0PN degeneracy feature are insensitive to that truncation.
minor comments (5)
- Fig. 1 caption and text: clarify whether the reported speed-up factors include the cost of assembling the full multi-parameter Fisher matrix or only the derivative evaluation itself.
- Eq. (35) and the subsequent scaling discussion introduce an index n that is easily confused with the PN order n_PN; a brief notational remark would help.
- Table 1 lists f_high = 1 Hz for all space-based detectors; a short justification relative to the actual transfer-function cut-offs would be useful.
- In Appendix B.7–B.9 the simplified constants F_α, F_β and the sound-speed dependence of γ_DF are stated without numerical values or ranges; adding them would improve reproducibility of Fig. 6.
- A few typographical inconsistencies appear (e.g., “Liuet al.”, “O’Beirnet al.” in the introduction; missing spaces before citations).
Circularity Check
No significant circularity: analytic ppE–TaylorF2 derivatives and Fisher forecasts are direct consequences of the stated waveform model; self-citations supply only background assumptions and noise curves.
specific steps
-
self citation load bearing
[Sec. 2 (paragraphs after Eq. 4) and citation [43]]
"A related assessment examined how different PN waveform models affect parameter estimation for SBBHs observed by space-based GW detectors [43]. … In particular, the 2PN inspiral waveform is sufficient for the Fisher-matrix forecast and scaling analysis performed here. … A related assessment is presented in our previous work [43]."
The justification that the 2PN truncation (and the neglect of eccentricity/spin-precession systematics) remains adequate for the present forecasts rests solely on the authors’ own prior paper. The assumption is not re-derived or independently validated here; it is imported by self-citation. The step is minor because it affects only the domain of validity of the forecasts, not the algebraic construction of the analytic derivatives themselves.
full rationale
The derivation chain is self-contained. Section 2 writes the standard TaylorF2 GR phase (Eqs. 1–4) plus the conventional ppE amplitude/phase corrections (Eqs. 5–10). Section 3 then obtains the Fisher integrands by elementary logarithmic differentiation (Eqs. 16–27), with all GR phase derivatives listed explicitly in Appendix A and the theory-specific Ξ, ϒ derivatives in Appendix B; none of these steps is defined in terms of the final constraints. The subsequent numerical results are ordinary Fisher forecasts under stated detector curves, population models and the 2PN truncation. Self-citations ([43], [53], [55]) justify the 2PN baseline, noise models and the neglect of TDI, but they are not used to force the analytic expressions or the reported PN-order/detector trends. There is no parameter fitting followed by “prediction,” no uniqueness theorem imported from the authors, and no renaming of a known empirical pattern. Circularity burden is therefore minimal.
Axiom & Free-Parameter Ledger
free parameters (2)
- PN truncation order of GR baseline
- Observation-time cut-offs T_obs
axioms (4)
- domain assumption Stationary-phase approximation yields the frequency-domain TaylorF2 inspiral waveform truncated at 2PN
- domain assumption Fisher-matrix covariance equals the inverse of the noise-weighted inner product of waveform derivatives (high-SNR Gaussian noise)
- domain assumption Sky- and polarization-averaged detector response can be absorbed into an effective sensitivity curve with C=1
- domain assumption ppE amplitude and phase corrections of the form α u^a and β u^b capture the leading inspiral deviations of the theories considered
read the original abstract
Testing gravity beyond general relativity (GR) is essential for probing fundamental physics in the strong-field and highly dynamical regime accessed by gravitational-wave (GW) observations. In this work, we derive analytic expressions for waveform derivatives in the Fisher-matrix formalism within the parametrized post-Einsteinian framework, using the frequency-domain inspiral waveform. These analytic derivatives enable stable and efficient Fisher-matrix calculations without relying on finite-difference schemes. We apply this method to a wide range of detector configurations, including space-based, ground-based, and multiband observations, and combine it with different binary black hole population models. Our results reveal clear and systematic trends in the constraints on non-GR effects as functions of post-Newtonian order, detector type, and source population. They also demonstrate the complementarity between space- and ground-based detectors, particularly for effects that accumulate during the low-frequency inspiral. The analytic approach substantially reduces computational cost and avoids numerical systematics associated with step-size choices, making it well suited for large-scale parameter studies. These results provide robust forecasts for the capability of future GW observations to constrain a broad class of non-GR effects and environmental influences, highlighting the scientific potential of upcoming detector networks for precision tests of gravity.
Reference graph
Works this paper leans on
-
[1]
T. Baker, D. Psaltis, C. Skordis, Astrophys. J.802, 63 (2015). DOI 10.1088/0004-637X/802/1/63
-
[2]
Will, Living Rev
C.M. Will, Living Rev. Rel.17, 4 (2014). DOI 10. 12942/lrr-2014-4
2014
-
[3]
M. Kramer, et al., Phys. Rev. X11(4), 041050 (2021). DOI 10.1103/PhysRevX.11.041050
-
[4]
B. Bertotti, L. Iess, P. Tortora, Nature425, 374 (2003). DOI 10.1038/nature01997
-
[5]
Stairs, Living Rev
I.H. Stairs, Living Rev. Rel.6, 5 (2003). DOI 10.12942/ lrr-2003-5
2003
-
[6]
T. Clifton, P.G. Ferreira, A. Padilla, C. Skordis, Phys. Rept.513, 1 (2012). DOI 10.1016/j.physrep.2012.01. 001
-
[7]
Joyce, L
A. Joyce, L. Lombriser, F. Schmidt, Ann. Rev. Nucl. Part. Sci.66, 95 (2016). DOI 10.1146/ annurev-nucl-102115-044553
2016
-
[8]
E. Berti, et al., Class. Quant. Grav.32, 243001 (2015). DOI 10.1088/0264-9381/32/24/243001
-
[9]
T. Damour, G. Esposito-Farese, Phys. Rev. Lett.70, 2220 (1993). DOI 10.1103/PhysRevLett.70.2220
-
[10]
B.P. Abbott, et al., Phys. Rev. Lett.116(6), 061102 (2016). DOI 10.1103/PhysRevLett.116.061102
-
[11]
K.A. Postnov, L.R. Yungelson, Living Rev. Rel.17, 3 (2014). DOI 10.12942/lrr-2014-3
-
[12]
A. Ghosh, et al., Phys. Rev. D94(2), 021101 (2016). DOI 10.1103/PhysRevD.94.021101
-
[13]
A. Ghosh, N.K. Johnson-Mcdaniel, A. Ghosh, C.K. Mishra, P. Ajith, W. Del Pozzo, C.P.L. Berry, A.B. Nielsen, L. London, Class. Quant. Grav.35(1), 014002 (2018). DOI 10.1088/1361-6382/aa972e
-
[14]
B.P. Abbott, et al., Phys. Rev. Lett.116(22), 221101 (2016). DOI 10.1103/PhysRevLett.116.221101. [Er- ratum: Phys.Rev.Lett. 121, 129902 (2018)]
-
[15]
S. Bhagwat, X.J. Forteza, P. Pani, V . Ferrari, Phys. Rev. D101(4), 044033 (2020). DOI 10.1103/PhysRevD.101. 044033
-
[16]
A.G. Abac, et al., Phys. Rev. Lett.135(11), 111403 (2025). DOI 10.1103/kw5g-d732
-
[17]
R. Abbott, et al., Phys. Rev. D112(8), 084080 (2025). DOI 10.1103/PhysRevD.112.084080
-
[18]
Collaboration, the Virgo Collaboration, the KAGRA Collaboration
T.L.S. Collaboration, the Virgo Collaboration, the KAGRA Collaboration. Gwtc-4.0: Updating the gravitational-wave transient catalog with observations from the first part of the fourth ligo-virgo-kagra ob- serving run (2025). URLhttps://arxiv.org/abs/ 2508.18082
Pith/arXiv arXiv 2025
-
[19]
M. Punturo, et al., Class. Quant. Grav.27, 194002 (2010). DOI 10.1088/0264-9381/27/19/194002
-
[20]
Amaro-Seoane, et al
P. Amaro-Seoane, et al. Laser interferometer space an- tenna (2017). URLhttps://arxiv.org/abs/1702. 00786
2017
-
[21]
W.R. Hu, Y .L. Wu, Natl. Sci. Rev.4(5), 685 (2017). DOI 10.1093/nsr/nwx116
-
[22]
J. Luo, et al., Class. Quant. Grav.33(3), 035010 (2016). DOI 10.1088/0264-9381/33/3/035010
-
[23]
R.G. Cai, Z.K. Guo, B. Hu, C. Liu, Y . Lu, W.T. Ni, W.H. Ruan, N. Seto, G. Wang, Y .L. Wu, Fund. Res.4, 1072 (2024). DOI 10.1016/j.fmre.2023.10.007
-
[24]
Mishra, K.G
C.K. Mishra, K.G. Arun, B.R. Iyer, B.S. Sathyaprakash, Phys. Rev. D82, 064010 (2010). DOI 10.1103/ PhysRevD.82.064010
2010
-
[25]
N. Yunes, F. Pretorius, Phys. Rev. D80, 122003 (2009). DOI 10.1103/PhysRevD.80.122003
-
[26]
N. Cornish, L. Sampson, N. Yunes, F. Pretorius, Phys. Rev. D84, 062003 (2011). DOI 10.1103/PhysRevD.84. 062003
-
[27]
K. Chatziioannou, N. Yunes, N. Cornish, Phys. Rev. D86, 022004 (2012). DOI 10.1103/PhysRevD.86. 022004. [Erratum: Phys.Rev.D 95, 129901 (2017)]
-
[28]
T. Narikawa, H. Tagoshi, PTEP2016(9), 093E02 (2016). DOI 10.1093/ptep/ptw126
-
[29]
N. Yunes, X. Siemens, K. Yagi, Living Rev. Rel.28(1), 3 (2025). DOI 10.1007/s41114-024-00054-9
-
[30]
S. Tahura, K. Yagi, Phys. Rev. D98(8), 084042 (2018). DOI 10.1103/PhysRevD.98.084042. [Erratum: Phys.Rev.D 101, 109902 (2020)]
-
[31]
Z. Carson, K. Yagi, Phys. Rev. D101(4), 044047 (2020). DOI 10.1103/PhysRevD.101.044047
-
[32]
C. Liu, L. Shao, J. Zhao, Y . Gao, Mon. Not. Roy. As- tron. Soc.496(1), 182 (2020). DOI 10.1093/mnras/ staa1512
doi:10.1093/mnras/ 2020
-
[33]
O’Beirne, N.J
L. O’Beirne, N.J. Cornish, S.J. Vigeland, S.R. Taylor, Phys. Rev. D99(12), 124039 (2019). DOI 10.1103/ 21 PhysRevD.99.124039
2019
-
[34]
C. Shi, M. Ji, J.d. Zhang, J. Mei, Phys. Rev. D108(2), 024030 (2023). DOI 10.1103/PhysRevD.108.024030
-
[35]
C. Liu, W.H. Ruan, Z.K. Guo, Phys. Rev. D102(12), 124050 (2020). DOI 10.1103/PhysRevD.102.124050
-
[36]
G. Wang, W.B. Han, Phys. Rev. D103(6), 064021 (2021). DOI 10.1103/PhysRevD.103.064021
-
[37]
J. Wu, J. Li, Phys. Rev. D110(8), 084057 (2024). DOI 10.1103/PhysRevD.110.084057
-
[38]
Blanchet, Living Rev
L. Blanchet, Living Rev. Rel.17, 2 (2014). DOI 10. 12942/lrr-2014-2
2014
-
[39]
A. Buonanno, B. Iyer, E. Ochsner, Y . Pan, B.S. Sathyaprakash, Phys. Rev. D80, 084043 (2009). DOI 10.1103/PhysRevD.80.084043
-
[40]
B.S. Sathyaprakash, B.F. Schutz, Living Rev. Rel.12, 2 (2009). DOI 10.12942/lrr-2009-2
-
[41]
T. Robson, N.J. Cornish, C. Liu, Class. Quant. Grav. 36(10), 105011 (2019). DOI 10.1088/1361-6382/ ab1101
-
[42]
S. Khan, S. Husa, M. Hannam, F. Ohme, M. Pür- rer, X. Jiménez Forteza, A. Bohé, Phys. Rev. D93(4), 044007 (2016). DOI 10.1103/PhysRevD.93.044007
-
[43]
J. Wu, J. Li, X. Liu, Z. Cao, Phys. Rev. D109(10), 104014 (2024). DOI 10.1103/PhysRevD.109.104014
-
[44]
H. Phurailatpam, G. Achamveedu, M. Haney, T. Li, S. Tiwari. On constraining initial orbital eccentricity of purely inspiral events (2025). URLhttps://arxiv. org/abs/2508.12697
Pith/arXiv arXiv 2025
-
[45]
X. Liu, Z. Cao, Z.H. Zhu, Class. Quant. Grav.41(19), 195019 (2024). DOI 10.1088/1361-6382/ad72ca
-
[46]
C.B. Owen, C.J. Haster, S. Perkins, N.J. Cornish, N. Yunes, Phys. Rev. D108(4), 044018 (2023). DOI 10.1103/PhysRevD.108.044018
-
[47]
K. Kremer, et al., Phys. Rev. D99(6), 063003 (2019). DOI 10.1103/PhysRevD.99.063003
-
[49]
M. Vallisneri, Phys. Rev. D77, 042001 (2008). DOI 10.1103/PhysRevD.77.042001
-
[50]
C. Cutler, E.E. Flanagan, Phys. Rev. D49, 2658 (1994). DOI 10.1103/PhysRevD.49.2658
-
[51]
Z. Luo, Z. Guo, G. Jin, Y . Wu, W. Hu, Results Phys.16, 102918 (2020). DOI 10.1016/j.rinp.2019.102918
-
[52]
S.J. Huang, Y .M. Hu, V . Korol, P.C. Li, Z.C. Liang, Y . Lu, H.T. Wang, S. Yu, J. Mei, Phys. Rev. D102(6), 063021 (2020). DOI 10.1103/PhysRevD.102.063021
-
[53]
J. Wu, J. Li, Phys. Rev. D108(12), 124047 (2023). DOI 10.1103/PhysRevD.108.124047
-
[54]
M. Vallisneri, J. Crowder, M. Tinto, Class. Quant. Grav. 25, 065005 (2008). DOI 10.1088/0264-9381/25/6/ 065005
-
[55]
J. Wu, M. Sun, J. Li, Class. Quant. Grav.43(6), 065017 (2026). DOI 10.1088/1361-6382/ae5144
-
[56]
S. Hild, et al., Class. Quant. Grav.28, 094013 (2011). DOI 10.1088/0264-9381/28/9/094013
-
[57]
E. Poisson, Phys. Rev. D52, 5719 (1995). DOI 10.1103/PhysRevD.52.5719. [Addendum: Phys.Rev.D 55, 7980–7981 (1997)]
-
[58]
M. Colpi, A. Sesana,Gravitational Wave Sources in the Era of Multi-Band Gravitational Wave Astronomy (2017), pp. 43–140. DOI 10.1142/9789813141766_ 0002
-
[59]
N. Aghanim, et al., Astron. Astrophys.641, A6 (2020). DOI 10.1051/0004-6361/201833910. [Erratum: As- tron.Astrophys. 652, C4 (2021)]
-
[60]
R. Abbott, et al., Phys. Rev. X13(1), 011048 (2023). DOI 10.1103/PhysRevX.13.011048
-
[61]
R. Abbott, et al., Astrophys. J. Lett.913(1), L7 (2021). DOI 10.3847/2041-8213/abe949
-
[62]
H.T. Wang, et al., Phys. Rev. D100(4), 043003 (2019). DOI 10.1103/PhysRevD.100.043003
-
[63]
S.X. Yi, G. Nelemans, C. Brinkerink, Z. Kostrzewa- Rutkowska, S.T. Timmer, F. Stoppa, E.M. Rossi, S.F. Portegies Zwart, Astron. Astrophys.663, A155 (2022). DOI 10.1051/0004-6361/202141634
-
[64]
S.X. Yi, F. Stoppa, G. Nelemans, E. Cator, Astron. As- trophys.663, A156 (2022). DOI 10.1051/0004-6361/ 202243127
-
[65]
Z. Zhao, S. Shashank, D. Das, C. Bambi. Testing black hole metrics with binary black hole inspirals (2025). URLhttps://arxiv.org/abs/2510.04703
Pith/arXiv arXiv 2025
-
[66]
E. Harikumar, V .O. Rivelles, Class. Quant. Grav.23, 7551 (2006). DOI 10.1088/0264-9381/23/24/024
-
[67]
H. Song, H. Li, Z. Lyu, J. Zhu, J.C. Wang, P. Ji. Con- straining noncommutative spacetime with gw150914 and gw190814 (2025). URLhttps://arxiv.org/ abs/2507.23442
arXiv 2025
-
[68]
R. Jackiw, S.Y . Pi, Phys. Rev. D68, 104012 (2003). DOI 10.1103/PhysRevD.68.104012
-
[69]
T.L. Smith, A.L. Erickcek, R.R. Caldwell, M. Kamionkowski, Phys. Rev. D77, 024015 (2008). DOI 10.1103/PhysRevD.77.024015
-
[70]
Wang, P.C
H.T. Wang, P.C. Li, J.L. Jiang, G.W. Yuan, Y .M. Hu, Y .Z. Fan, Eur. Phys. J. C81(8), 769 (2021). DOI 10. 1140/epjc/s10052-021-09555-1
2021
-
[71]
T. Jacobson, Phys. Rev. Lett.83, 2699 (1999). DOI 10.1103/PhysRevLett.83.2699
-
[72]
M.W. Horbatsch, C.P. Burgess, JCAP05, 010 (2012). DOI 10.1088/1475-7516/2012/05/010
-
[73]
K.i. Maeda, N. Ohta, Y . Sasagawa, Phys. Rev. D80, 104032 (2009). DOI 10.1103/PhysRevD.80.104032
-
[74]
Gao, S.P
B. Gao, S.P. Tang, H.T. Wang, J. Yan, Y .Z. Fan, Phys. Rev. D110(4), 044022 (2024). DOI 10.1103/ PhysRevD.110.044022 22
2024
-
[75]
J. Zhang, S.Y . Zhou, Phys. Rev. D97(8), 081501 (2018). DOI 10.1103/PhysRevD.97.081501
-
[76]
N. Yunes, F. Pretorius, D. Spergel, Phys. Rev. D81, 064018 (2010). DOI 10.1103/PhysRevD.81.064018
-
[77]
K. Yagi, N. Tanahashi, T. Tanaka, Phys. Rev. D83, 084036 (2011). DOI 10.1103/PhysRevD.83.084036
-
[78]
A. Genova, E. Mazarico, S. Goossens, F.G. Lemoine, G.A. Neumann, D.E. Smith, M.T. Zu- ber, Nature Commun.9(1), 289 (2018). DOI 10.1038/s41467-017-02558-1
-
[79]
A. Toubiana, et al., Phys. Rev. Lett.126(10), 101105 (2021). DOI 10.1103/PhysRevLett.126.101105
-
[80]
G. Caneva Santoro, S. Roy, R. Vicente, M. Haney, O.J. Piccinni, W. Del Pozzo, M. Martinez, Phys. Rev. Lett. 132(25), 251401 (2024). DOI 10.1103/PhysRevLett. 132.251401
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.