{"id":"70603ee4-c0b3-4c02-a10f-431bf627db30","arxiv_id":"2412.15578","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A perturbation-theory treatment of gravitational waves in a Robertson-Walker universe yields a next-to-leading-order correction factor to the waveform, but the quantitative expression appears to contain sign and factor errors.","lead":"This paper derives a next-order correction to how gravitational waves propagate through an expanding universe, claiming the expansion imprints a small frequency- and redshift-dependent factor on the waveform. The result could in principle help separate redshift from chirp mass in gravitational-wave observations, but the effect is tiny and the derived formulas contain errors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (23) does not follow from Eq. (18): converting the derived correction to Lambda-CDM yields a different integrand, so the quantitative claims (amplification, degeneracy breaking) lack a valid derivation.","rationale":"The paper's central claim is that a next-to-leading-order WKB correction to gravitational-wave propagation in RW spacetime, Eq. (20), is given quantitatively by Eq. (23) in Lambda-CDM, and that it produces amplification at low frequencies and breaks the chirp-mass/redshift degeneracy. For that claim to hold, Eq. (23) must be the Lambda-CDM specialization of Eq. (18). I checked the WKB step: Eq. (18) is consistent with the time ODE Eq. (15), so the derivation up to the general F^- is credible. The load-bearing weak point is the Lambda-CDM conversion. A standard change of variable x=a(t) with the Friedmann expressions yields a structurally different integrand, with powers x^{-4}, x^{-3}, x^{-2} and a different dark-energy term, rather than the sigma(x) in Eq. (24). The factor 1/(2 i omega)=1/(4 pi i f) also exposes a missing pi in Eq. (23). Because all numerical claims and the degeneracy-breaking argument use Eq. (23), the central result is not reliably established as written. I also note the paper's own validity estimate at z=5000, f=10^{-15} Hz is numerically off by about an order of magnitude; this does not by itself invalidate the WKB derivation but reinforces the need for a corrected conversion. No code or data is provided. The reader's REJECT verdict stands; I do not see a need to alter it, though the reason is more specific: Eq. (23) is not the integral of Eq. (18).","tokens_in":7062,"tokens_out":26510,"duration_ms":213534,"concrete_test":"Independently convert Eq. (18) to redshift using the standard Lambda-CDM Friedmann equation: substitute x=a(t), use dt=dx/(xH), and compare the resulting integrand term by term with sigma(x) in Eq. (24). Then evaluate both the corrected F^- and the paper's Eq. (23) at z=5000 for f=10^{-12} Hz and f=10^{-15} Hz. If the integrand does not match Eq. (24) in the x-powers, or if the numerical corrections disagree beyond the pi/coefficient factor, the central quantitative formula is invalid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (18) is internally consistent with Eq. (15): substituting T=sqrt(a) F e^{±i omega ∫ dt/a} and keeping first order gives 2 i omega a F_dot + a^2 F_ddot - (3 a_dot^2 + 5 a a_ddot) F = 0, hence Eq. (18) at order 1/omega. The failure is the Lambda-CDM specialization. With x=a(t), H^2/H0^2 = Omega_r x^{-4} + Omega_m x^{-3} + Omega_k x^{-2} + Omega_DE x^{-3(1+w)} and dt = dx/(x H), Eq. (18) integrates to F^- = 1 - (H0/(2 i omega)) ∫_{1/(1+z)}^1 [ -2 Omega_r x^{-4} + (1/2) Omega_m x^{-3} + 3 Omega_k x^{-2} + (8 - 7.5(1+w)) Omega_DE x^{-3(1+w)} ] / sqrt(D) dx, with D the Friedmann sum. This is not Eq. (24): sigma(x) contains Omega_m x and Omega_r x^2, no x^{-3} or x^{-4} terms, and a different denominator. The coefficient is also 1/(2 i omega)=1/(4 pi i f), not 1/(4 i f), so Eq. (23) is missing a factor pi. Since Eq. (23) is the only quantitative bridge to the claimed amplification at 10^{-12} Hz and to breaking the chirp-mass/redshift degeneracy, those claims are unsupported. The Discussion validity check is also inconsistent with Lambda-CDM: at z=5000 and f=10^{-15} Hz, a_dot/H0 ≈ sqrt(Omega_r)(1+z) ≈ 45 while omega/H0 ≈ 2.7e3, so (a_dot/omega)^2 ≈ 2.7e-4, not ~1e-5.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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).","tokens_in":7494,"tokens_out":18225,"duration_ms":148001,"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":[{"comment":"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).","section":"Eq. (20)"},{"comment":"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.","section":"Eqs. (18), (23)-(24)"},{"comment":"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.","section":"Discussion"},{"comment":"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).","section":"Abstract and Eq. (18)"}],"minor_comments":[{"comment":"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.","section":"Eqs. (14)-(18)"},{"comment":"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.","section":"Eq. (16)"},{"comment":"The text says 'axis perturbation' where 'axial perturbation' is intended; this typo should be corrected.","section":"Final paragraph"},{"comment":"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.","section":"Eq. (22) and Ref. [24]"}],"recommendation":"reject","confidential_remarks":"The two load-bearing errors—the sign inconsistency in Eq. (20) and the incorrect Lambda-CDM reduction in Eq. (23)—are decisive for rejection. A revision that merely fixes typographical issues will not suffice; the derivation of the propagation factor and its cosmological specialization would need to be redone, and the claimed low-frequency amplification would need to be re-evaluated with the corrected expression."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Clever formal machinery, but the central equations don't hold together. The paper's real contribution is using tensor spherical harmonics on a Robertson-Walker background and writing down the next-to-leading-order amplitude correction F^-; that is new relative to the cited geometric-optics and flat-space papers, and the phase limit matching Ref. [24] is a sensible consistency check. The derivation to Eqs. (15)-(18) is plausible, and the internal consistency of Eq. (18) with the time equation checks out.\n\nThe problems start at Eq. (20). With eta(t) defined as i integral dt'/a, the exponent i omega [beta(r) - eta(t)] becomes i omega beta(r) + omega integral dt'/a, a growing real exponential, not a phase; it also fails to reduce to exp[i omega (r - t)] in the flat limit. That is a sign error or a bad convention, and it undercuts the waveform Eq. (20) as written.\n\nEq. (23) is worse. The Lambda-CDM specialization does not follow from Eq. (18) via the Friedmann equations. Converting the integral properly gives an integrand with Omega_r x^{-4} and Omega_m x^{-3} terms and a different denominator from the sigma(x) in Eq. (24); the prefactor is 1/(4 pi i f), not the printed 1/(4 i f). So the low-frequency form of F^-, the claimed amplification near 10^-12 Hz, and the chirp-mass/redshift degeneracy argument all rest on equations that are not derived. The order-of-magnitude check at z=5000, f=10^-15 Hz also comes out around 10 times larger than the paper's estimate (still small, but the numbers as stated don't hold).\n\nNone of this means the idea is dead. The formal approach is worth pursuing, and the flaws look fixable rather than fundamental. But the version on arXiv overstates both the novelty ('first time', 'obvious amplification') and the quantitative support. I would send it to a referee — the method and the claim are substantive enough to deserve expert scrutiny — but the referee's job is clear: verify Eqs. (18)-(24), and don't let the Letter proceed until the Lambda-CDM conversion is done correctly. As printed, I would not accept it.","headline":"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.","tokens_in":8011,"tokens_out":7040,"would_cite":false,"duration_ms":55205,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83F05","83C25"],"pacs":["04.30.-w","98.80.-k"],"model":"deepseek-v4-flash","headline":"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.","keywords":["gravitational waves","Robertson-Walker metric","cosmological perturbation","chirp mass-redshift degeneracy","standard sirens","tensor spherical harmonics","cosmic acceleration","low-frequency gravitational waves"],"falsifier":"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.","tokens_in":6799,"feed_emoji":"🌌","tokens_out":8096,"duration_ms":69443,"temperature":0.7,"pith_summary":"This paper tries to establish that gravitational waves propagating through an expanding universe acquire a small but calculable correction to their amplitude, on top of the usual redshift and phase shift. The authors derive the correction by perturbing the Robertson-Walker metric with tensor spherical harmonics and solving the wave equation to next-to-leading order in the expansion rate. In a Λ-CDM universe the correction factor is $F^- = 1 - \\frac{H_0}{4if}\\int_{1/(1+z)}^{1}\\sigma(x)\\,dx$, and it becomes significant for high-redshift, ultra-low-frequency waves. The result matters because it would break the chirp-mass/redshift degeneracy that otherwise prevents gravitational-wave signals from giving redshifts directly.","feed_headline":"Cosmic expansion imprints GWs, breaking mass-redshift degeneracy","feed_subtitle":"A next-order waveform correction could let standard sirens measure redshift directly.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the tensor spherical-harmonic decomposition method used to reduce the Robertson-Walker perturbation equations into radial and time parts.","marker":"[27]"},{"why":"Supplies the geometric-optics phase correction for gravitational waves in a flat Λ-CDM universe that the next-to-leading-order amplitude correction extends and compares against.","marker":"[24]"},{"why":"Establishes the chirp-mass/redshift degeneracy as an exact invariance under the scaling transformation that the new waveform is claimed to break.","marker":"[20]"}],"fun_headline_variants":["Expansion history imprinted in gravitational waves","Cosmic expansion breaks GW mass-redshift degeneracy","Next-order GW effect from cosmic acceleration","Gravitational waves carry expansion history, lift degeneracy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Expansion history imprinted in gravitational waves","Cosmic expansion breaks GW mass-redshift degeneracy","Next-order GW effect from cosmic acceleration","Gravitational waves carry expansion history, lift degeneracy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":1859,"prompt_tokens":876,"completion_tokens":983,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":492,"completion_tokens_details":{"reasoning_tokens":924}},"tokens_in":492,"tokens_out":983,"duration_ms":9338,"temperature":1.0,"reasoning_tokens":924,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:18:16.052913+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Cutler and E","cited_arxiv_id":null,"evidence_quote":"Supplies the tensor spherical-harmonic decomposition method used to reduce the Robertson-Walker perturbation equations into radial and time parts."},{"cited_title":"Smarra, B","cited_arxiv_id":null,"evidence_quote":"Supplies the geometric-optics phase correction for gravitational waves in a flat Λ-CDM universe that the next-to-leading-order amplitude correction extends and compares against."},{"cited_title":"Discussion— In this Letter, we study the perturbation of the R W metric and separate variables using spheri- cal tensors","cited_arxiv_id":null,"evidence_quote":"Establishes the chirp-mass/redshift degeneracy as an exact invariance under the scaling transformation that the new waveform is claimed to break."}],"review_version":1}