REVIEW 4 major objections 4 minor 37 references
The imprint of cosmic expansion history on the propagation of gravitational waves
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper derives a next-to-leading-order correction to gravitational-wave amplitudes from propagation through an expanding Robertson-Walker universe and argues it breaks the chirp-mass/redshift degeneracy.
desk verdict Clever formal setup undercut by a sign error in the central waveform and a Lambda-CDM conversion that doesn't follow from the derived equations — claims too big for the math as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the amplitude correction factor $F^-(\omega,t)$, which emerges from separating the metric perturbation into angular modes with tensor spherical harmonics. The radial part gives the oscillatory $1/r$ wave with comoving frequency $\omega$, while the time part obeys $(-\frac{11}{4}\dot{a}^2 - \frac{11}{2}a\ddot{a} + \omega^2)T + a^2\ddot{T}=0$; solving this under the short-wavelength assumption $\dot{a}\ll\omega$ and $\ddot{a}\ll\omega^2$ produces $F^\pm$. This factor is what converts the expansion history into a concrete waveform correction, and its $1/\omega$ form is what makes the effect grow at low frequencies.
What would settle it
Compute the exact numerical solution of the separated time equation for a Λ-CDM background at $f=10^{-12}$ Hz and $z=5000$ without assuming $\dot{a}\ll\omega$ and $\ddot{a}\ll\omega^2$, and compare the resulting amplitude at detection with Eqs. (20)-(23); if the exact solution disagrees with the WKB correction by more than the claimed effect, the central claim fails in that regime.
Extended reading notes
Core claim
The central claim is that the outgoing gravitational wave in a Robertson-Walker background is, to next-to-leading order, $h_+ - ih_\times \sim \frac{F^-(\omega,t)}{a(t)r} e^{i\omega[\beta(r)-\eta(t)]}$, with $F^\pm = 1 \pm \frac{1}{2i\omega}\int_{t_{\rm emis}}^{t}(3\dot{a}^2/a + 5\ddot{a})\,dt'$. For a flat Λ-CDM universe this becomes Eq. (23), an explicit integral over the density parameters $\Omega_r$, $\Omega_m$, $\Omega_k$, $\Omega_\Lambda$ and the dark-energy equation-of-state parameter $d$. Because the factor depends on the whole expansion history along the propagation path, the waveform carries cosmological information that is absent in the standard flat-space result. The authors further show that the corrected waveform is not invariant under the scaling $(M_c,L,t)\to(M_c\lambda,L\lambda,t\lambda)$, so the degeneracy between chirp mass and redshift is lifted in principle.
Load-bearing premise
The correction factor stands on the assumption that the expansion rate and acceleration of the universe are much smaller than the gravitational-wave frequency at every point along the propagation path.
Editorial extensions
If this is right
- Gravitational waves at frequencies $\lesssim 10^{-12}$ Hz would be noticeably amplified by the correction factor, making the effect relevant for proposed ultra-low-frequency observations.
- With the degeneracy broken, a single binary inspiral could in principle yield both chirp mass and redshift, rather than only the redshifted combination, if the expansion-history term is known.
- The integral over $\sigma(x)$ offers a direct probe of the cosmological density parameters and the dark-energy equation of state from gravitational-wave propagation alone.
- Neglecting $F^-$ in parameter estimation at high redshift and low frequency would introduce a systematic error in the inferred source parameters.
- The radial equation suggests a cutoff frequency $\omega^2+3k<0$, below which GW modes would not oscillate; this is a testable prediction of the same formalism.
Reading between the lines
- A natural next step would be to solve the time equation exactly for a Λ-CDM background at $f\sim10^{-12}$ Hz, without the short-wavelength approximation, to check whether the amplitude correction survives where it becomes large.
- If the correction holds, the same mechanism should show up in gravitational-wave memory or in stochastic backgrounds; the frequency dependence $1/f$ could distinguish this cosmological effect from astrophysical or detector noise.
- Because the correction depends on the equation-of-state parameter $d$, a sufficiently loud high-redshift event could, in principle, constrain dark energy; whether any realistic detector has the required sensitivity remains open.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript treats odd-parity tensor perturbations of the Robertson-Walker metric using a Regge-Wheeler decomposition, separates the radial and time equations, solves the time equation in a slow-expansion WKB approximation, and obtains a propagation factor F^-(ω,t) given by Eq. (20). It then specializes this factor to Lambda-CDM, Eq. (23)-(24), and claims that it amplifies gravitational waves at frequencies below about 10^-12 Hz and that the corrected waveform Eq. (22) breaks the usual chirp-mass/redshift degeneracy. The derivation uses no fitted parameters, and the flat-spacetime limit is stated as a benchmark. The central quantitative claims, however, rest on an internally inconsistent exponent in Eq. (20) and on a Lambda-CDM reduction of Eq. (18) that does not produce Eq. (23).
Significance. If the central result were correct, the paper would offer a genuinely interesting next-to-leading-order correction to gravitational-wave propagation in a cosmological background, with possible consequences for standard-siren cosmology and for the chirp-mass/redshift degeneracy. The approach is standard and the paper is not circular: the derivation starts from Einstein equations and no parameter is fitted to force the claimed correction. These strengths are, however, outweighed by load-bearing errors. The sign error in Eq. (20) converts a phase factor into a spurious exponential growth, and Eq. (23) does not follow from Eq. (18). Consequently, the abstract's amplification claim and the degeneracy-breaking claim are not supported by the present derivation.
major comments (4)
- [Eq. (20)] The phase in Eq. (20) is inconsistent with the definition of η(t) given just before Eq. (17). With η(t)=i∫_{t_emis}^t dt'/a(t'), the exponential exp{iω[β(r)-η(t)]} equals exp{iωβ(r)} exp{ω∫_{t_emis}^t dt'/a(t')}, which grows exponentially with time rather than describing an outgoing wave. In the flat limit a=1 and β(r)=r this reduces to exp(iωr) exp(ω(t-t_emis)), not the exp[iω(r-t)]/r stated in the text. The outgoing solution should carry the factor e^{-ωη}=e^{-iω∫dt/a} that appears in Eq. (17); as written, Eq. (20) is the source of the claimed low-frequency amplification and cannot be used to obtain Eq. (22).
- [Eqs. (18), (23)-(24)] Equation (23) does not follow from Eq. (18). Substituting x=a(t), H^2/H0^2 = Ω_r x^-4 + Ω_m x^-3 + Ω_k x^-2 + Ω_Λ x^{-3(1+w)}, and dt=dx/(xH) into the integral of Eq. (18) gives F^- = 1 - (H0/(2iω)) ∫_{1/(1+z)}^1 [-2Ω_r x^-4 + (1/2)Ω_m x^-3 + 3Ω_k x^-2 + (8 - (15/2)(1+w))Ω_Λ x^{-3(1+w)}] / sqrt(D) dx, where D is the Friedmann sum. This integrand is different from Eq. (24): the numerator has x^-3 and x^-4 terms that are absent from the quoted σ(x), and the denominator is not the one used there. In addition, with ω=2πf the prefactor should be 1/(4π i f), not 1/(4 i f); the factor π is missing. The quantitative claims for f≲10^-12 Hz and the degeneracy-breaking claim therefore rest on an equation that is not derived.
- [Discussion] The numerical check of the slow-expansion approximation is inconsistent with Lambda-CDM parameters. At z=5000, a≈2×10^-4, Ω_r≈9×10^-5, and H0≈2.3×10^-18 s^-1, one has dot a=aH≈√Ω_r H0(1+z)≈1.1×10^-16 and |ddot a|≈aH^2≈6×10^-29 s^-2. For f=10^-15 Hz, using ω=2πf gives (dot a/ω)^2≈3×10^-4 and |ddot a|/ω^2≈1, not the quoted 10^-5 and 10^-9; if the text literally means ω=10^-15 s^-1, both ratios are even larger. Thus the WKB assumption ddot a/ω^2≪1 is not satisfied at the claimed parameters, and the stated validity estimate is not supported.
- [Abstract and Eq. (18)] Even after correcting the exponent issue, the first-order correction F^- in Eq. (18) is purely imaginary: F^- = 1 - (1/(2iω))∫(3dot a^2/a + 5ddot a)dt = 1 + (i/(2ω))∫(3dot a^2/a + 5ddot a)dt. Hence |F^-| = 1 + O(1/ω^2), and there is no first-order amplitude amplification. The abstract's statement that the effect 'will obviously amplify GWs' at low frequencies therefore does not follow from the derived F^-; it follows only from the growing exponential produced by the sign error in Eq. (20).
minor comments (4)
- [Eqs. (14)-(18)] The separation constant is denoted ω but is never identified as an angular frequency, and Eq. (23) later uses f without a 2π conversion; the missing factor π in the prefactor appears related to this ambiguity.
- [Eq. (16)] The replacement α→ω uses the condition k/ω^2≪1, but for the low frequencies discussed here this inequality can fail; the domain of validity of the radial solution and of the O(r^-2) dropping should be stated explicitly.
- [Final paragraph] The text says 'axis perturbation' where 'axial perturbation' is intended; this typo should be corrected.
- [Eq. (22) and Ref. [24]] The statement that the first-order phase term agrees with Ref. [24] is not checkable as written because it relies on Eq. (20), whose exponent has the sign error; the comparison should be re-verified after correcting the propagation factor.
Circularity Check
No significant circularity: the derivation is self-contained and not reduced to its inputs.
full rationale
The derivation is not circular. The correction factor F^- is obtained from an explicit perturbative solution of the Robertson-Walker background equations (Eqs. 2-20), with the asymptotic assumptions ȧ << ω and ä << ω² stated before Eq. (17). No parameter is fitted to target data, and the ΛCDM specialization in Eqs. (23)-(24) is a substitution of the Friedmann expansion law into the derived integral, not an input used to define the result. The flat-spacetime limit (k=0, a=1) and the phase limit attributed to Ref. [24] are external benchmarks, so the central waveform Eq. (22) is not validated by its own assumptions. A possible algebraic mismatch between Eq. (18) and Eq. (23) would be a derivation error, not circularity; circularity would require the claimed prediction to be identical by construction to a fitted input or to a self-citation, and no such identity appears in the paper.
Assumptions & free parameters
assumptions (5)
- domain assumption The universe on large scales is described by the Robertson-Walker metric with homogeneity and isotropy, and the Friedmann equations hold.
- domain assumption The GW is a small metric perturbation with delta T = 0; matter is not perturbed by the wave.
- ad hoc to paper The expansion is slow compared to the GW frequency: dot-a/omega << 1 and double-dot-a/omega^2 << 1.
- ad hoc to paper The radial equation can be solved assuming 1-k r^2 ~ O(1) and dropping O(r^-2) terms.
- ad hoc to paper Axial (odd-parity) perturbations are sufficient to describe GW propagation from any source.
Cite this review
Pith. "Pith review of The imprint of cosmic expansion history on the propagation of gravitational waves." pith.science (2026). https://pith.science/paper/WJUVG42M
@misc{pith2026241215578,
author = {Pith},
title = {Pith review of: The imprint of cosmic expansion history on the propagation of gravitational waves},
year = {2026},
howpublished = {\url{https://pith.science/paper/WJUVG42M}},
note = {Machine review of arXiv:2412.15578}
}
abstract
Gravitational waves (GWs) are regarded as standard sirens for Cosmology. GWs from compact binary coalescence (CBC) can directly determine the luminosity distance but usually can not obtain information about the redshift. However, if the universe is not flat but accelerating, GWs should carry this cosmological effect. In this Letter, for the first time, we explore how the expansion of the Universe affects GW propagation by perturbing the Robertson-Walker metric. We achieve a comprehensive and rigorous formalism at the next-leading order to describe the cosmological acceleration in GWs from any kind of sources. Theoretically, this cosmological effect will obviously amplify GWs at frequencies as low as $10^{-12}$ Hz.
Reference graph
Works this paper leans on
- [24]
-
[1]
B. P. Abbott et al. (LIGO Scientific Collabora- tion and Virgo Collaboration), Observation of grav- itational waves from a binary black hole merger, Phys. Rev. Lett. 116, 061102 (2016)
work page 2016
-
[2]
( 3), we have δSµν = 3 8π ( ¨a a + k + ˙a2 a2 ) hµν
and Eq. ( 3), we have δSµν = 3 8π ( ¨a a + k + ˙a2 a2 ) hµν . (7) Following Regge and Wheeler, we decompose h into tensor functions with the correct symmetry behavior un- der rotation on the 2-sphere dσ 2. For the axial pertur- bation (with the parity ( −)l+1) for given l and m, h has the same form as that of Regge and Wheeler [27], where the Ylm(θ, φ ) i...
-
[3]
F. Acernese et al. , Advanced virgo: a second- generation interferometric gravitational wave detector, Classical and Quantum Gravity 32, 024001 (2014)
work page 2014
-
[4]
K. Somiya and (for the KAGRA Collabo- ration), Detector configuration of kagra–the japanese cryogenic gravitational-wave detector, Classical and Quantum Gravity 29, 124007 (2012)
work page 2012
-
[5]
are indepen- dent; therefore, we have three differential equations for only two unknown functions, h0 and h1. After some al- gebraic simplification, we obtain the differential equation satisfied by h1, r ( −ra ˙a ∂h 1 ∂t + ra2 ∂ 2h1 ∂t 2 + kr3 ∂ 2h1 ∂r 2 +kr2 ∂h 1 ∂r − r ∂ 2h1 ∂r 2 + 2 ∂h 1 ∂r ) + h1 ( 3r2a¨a + 7r2 ˙a2 + 3kr2 + l2 + l − 2 ) = 0, (12) and the ...
-
[6]
R. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Gwtc-2: Compact bi- nary coalescences observed by ligo and virgo dur- ing the first half of the third observing run, Phys. Rev. X 11, 021053 (2021)
work page 2021
-
[7]
will ef- fectively amplify the GWs’ amplitude
Show all 37 references
-
[8]
Agazie et al
G. Agazie et al. , The nanograv 15 yr data set: Observations and timing of 68 millisecond pulsars, The Astrophysical Journal Letters 951, L9 (2023)
2023
-
[9]
T. L. S. Collaboration, J. Aasi, et al. , Advanced ligo, Classical and Quantum Gravity 32, 074001 (2015)
2015
-
[10]
D. J. Reardon et al., Search for an isotropic gravitational- wave background with the parkes pulsar timing array, The Astrophysical Journal Letters 951, L6 (2023)
2023
-
[11]
H. Xu, S. Chen, et al. , Searching for the nano- hertz stochastic gravitational wave background with the chinese pulsar timing array data release i, Research in Astronomy and Astrophysics 23, 075024 (2023)
2023
-
[12]
B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Gwtc-1: A gravitational-wave transient catalog of compact binary mergers observed by ligo and virgo during the first and second observing runs, Phys. Rev. X 9, 031040 (2019)
2019
-
[13]
Gouttenoire, First-order phase transition interpreta- tion of pulsar timing array signal is consistent with solar- mass black holes, Phys
Y. Gouttenoire, First-order phase transition interpreta- tion of pulsar timing array signal is consistent with solar- mass black holes, Phys. Rev. Lett. 131, 171404 (2023)
2023
-
[14]
( 15) are the radial part and time part, respectively
and Eq. ( 15) are the radial part and time part, respectively. The solutions to these equa- tions describe the propagation behavior of GWs. From our example for CBC, the correction factor Eq. (
-
[15]
Abbott et al
R. Abbott et al. (LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration), Gwtc-3: Compact binary coalescences observed by ligo and virgo during the second part of the third observing run, Phys. Rev. X 13, 041039 (2023)
2023
-
[16]
, The second data release from the european pulsar timing array - i
EPTA Collaboration, Antoniadis, J., et al. , The second data release from the european pulsar timing array - i. the dataset and timing analysis, A&A 678, A48 (2023)
2023
-
[17]
Xue, Z.-Q
X. Xue, Z.-Q. Xia, et al. (PPTA Collabora- tion), High-precision search for dark photon dark matter with the parkes pulsar timing array, Phys. Rev. Res. 4, L012022 (2022)
2022
-
[18]
D. J. D’Orazio and A. Loeb, Using gravitational wave parallax to measure the hubble parameter with pulsar timing arrays, Phys. Rev. D 104, 063015 (2021)
2021
-
[19]
Roshan and G
R. Roshan and G. White, Using gravitational waves to see the first second of the universe (2024), arXiv:2401.04388 [hep-ph]
2024 arXiv
-
[20]
Discussion— In this Letter, we study the perturbation of the R W metric and separate variables using spheri- cal tensors
will become exp [iω (r − t)] /r which exactly is a monochromatic grav- itational wave in flat spacetime. Discussion— In this Letter, we study the perturbation of the R W metric and separate variables using spheri- cal tensors. Eq. (
-
[21]
X. Xue, L. Bian, J. Shu, Q. Yuan, X. Zhu, N. D. R. Bhat, S. Dai, Y. Feng, B. Goncharov, G. Hobbs, E. Howard, R. N. Manchester, C. J. Rus- sell, D. J. Reardon, R. M. Shannon, R. Spiewak, N. Thyagarajan, and J. Wang, Constraining cosmologi- cal phase transitions with the parkes ...
2021
-
[22]
( 18) are not invariant under this transformation, which indi- cates its potential to break the degeneracy between the chirp mass Mc and the redshift z
with amplitude correction Eq. ( 18) are not invariant under this transformation, which indi- cates its potential to break the degeneracy between the chirp mass Mc and the redshift z. There are still some interesting points for future work. In Eq. ( 16), if ω 2 + 3k < 0, the ra...
-
[23]
If one neglects its contribution, it may result in a non-negligible systemic error in parameter estimation
will become considerable at high redshift and low frequency, e.g., z > 1000 and f < 10−12 Hz. If one neglects its contribution, it may result in a non-negligible systemic error in parameter estimation. To obtain these equations, we have assumed that ˙a2/ω 2 → 0 and ¨a/ω 2 → 0....
-
[25]
Alonso- ´Alvarez, J
G. Alonso- ´Alvarez, J. M. Cline, and C. Dewar, Self-interacting dark matter solves the final par- sec problem of supermassive black hole mergers, Phys. Rev. Lett. 133, 021401 (2024)
2024
-
[26]
B. P. Abbott et al. , A gravitational-wave stan- dard siren measurement of the hubble constant, Nature 551, 85 (2017)
2017
-
[27]
Cutler and E
C. Cutler and E. E. Flanagan, Gravitational waves from merging compact binaries: How accurately can one ex- tract the binary’s parameters from the inspiral wave- form?, Phys. Rev. D 49, 2658 (1994)
1994
-
[28]
Riazuelo and J.-P
A. Riazuelo and J.-P. Uzan, Quintessence and gravita- tional waves, Phys. Rev. D 62, 083506 (2000)
2000
-
[29]
Nishizawa, Generalized framework for testing grav- ity with gravitational-wave propagation
A. Nishizawa, Generalized framework for testing grav- ity with gravitational-wave propagation. i. formulation, Phys. Rev. D 97, 104037 (2018)
2018
-
[30]
I. D. Saltas, I. Sawicki, L. Amendola, and M. Kunz, Anisotropic stress as a signature of nonstandard propagation of gravitational waves, Phys. Rev. Lett. 113, 191101 (2014)
2014
-
[31]
N. Seto, S. Kawamura, and T. Nakamura, Possibility of direct measurement of the acceleration of the universe using 0.1 hz band laser interferometer gravitational wave antenna in space, Phys. Rev. Lett. 87, 221103 (2001)
2001
-
[32]
J. M. Bardeen, Gauge-invariant cosmological perturba- tions, Phys. Rev. D 22, 1882 (1980)
1980
-
[33]
D. H. Lyth and D. Wands, Conserved cosmological per- turbations, Phys. Rev. D 68, 103515 (2003)
2003
-
[34]
Regge and J
T. Regge and J. A. Wheeler, Stability of a schwarzschild singularity, Phys. Rev. 108, 1063 (1957)
1957
-
[35]
Amaro-Seoane et al
P. Amaro-Seoane et al. (LISA Collaboration), Laser interferometer space antenna (2017), arXiv:1702.00786 [astro-ph.IM]
2017 arXiv
-
[36]
Hu and Y.-L
W.-R. Hu and Y.-L. Wu, The Taiji Program in Space for gravitational wave physics and the na- ture of gravity, National Science Review 4, 685 (2017) , https://academic.oup.com/nsr/article-pdf/4/5/685/31566708/nwx116
2017
-
[37]
Luo, L.-S
J. Luo, L.-S. Chen, H.-Z. Duan, Y.-G. Gong, S. Hu, J. Ji, Q. Liu, J. Mei, V. Milyukov, M. Sazhin, C.-G. Shao, V. T. Toth, H.-B. Tu, Y. Wang, Y. Wang, H.- C. Yeh, M.-S. Zhan, Y. Zhang, V. Zharov, and Z.-B. Zhou, Tianqin: a space-borne gravitational wave detec- tor, Classical an...
2016
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.