{"id":"8fb59a33-fbfc-407d-847a-955d053305ab","arxiv_id":"1908.08472","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A gravitational wave in an FLRW universe has an effective wave number different from its redshifted frequency, producing an angle-dependent enhancement in pulsar timing residuals that could measure H0.","lead":"This paper predicts that gravitational waves moving through an expanding universe can create a narrow, ring-like zone in the sky where pulsar timing signals get strongly amplified, and that the ring's angular size depends on the local value of the Hubble constant. The authors derive the effect, map out its dependence on source distance, wave frequency, and pulsar distance, and propose using it as a new, local probe of cosmic expansion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (15), the source of the (1+H0T) amplitude and the k_eff differing from w_eff phase, is asserted without derivation and is never reconciled with the standard FLRW tensor-mode equation, so the PTA observable in Eq. (21) may require a different h variable.","rationale":"The reader identifies Eq. (15) and the coordinate-transformation mapping as the weakest assumption; I agree with that identification. The whole paper is internally organized around Eq. (15): its unusual -H∂_T term produces the growing amplitude and the phase that determines the enhancement angle, and the same equation is used to claim consistency with the transformed SdS wave in Eq. (18). But the derivation of Eq. (15) is compressed into a single \"one gets\", and the paper does not specify whether h_ij is the dimensionless comoving strain D or the full metric perturbation δg_ij = a²D. In standard FLRW perturbation theory the two obey different equations, with opposite signs in the Hubble-friction term, and they enter the timing residual differently. This is not a mere disagreement with a convention; it is an internal consistency gap in the argument that connects the theoretical wave solution to the observable residual. The proposed check settles it: either an independent derivation reproduces Eq. (15) in the variables used in Eq. (21), or the claim that the enhancement measures H0 is not supported. Because the reader's REJECT verdict already rests on the same gap, I do not propose a different verdict; the paper should be strengthened by an explicit derivation and consistency check before it can be accepted.","tokens_in":19637,"tokens_out":22308,"duration_ms":249605,"concrete_test":"Derive the linearized Einstein equation for the exact TT perturbation h_μν whose spatial components are those of Eq. (18) in the FLRW ΛCDM background, without assuming Eq. (15), and check whether Eq. (18) satisfies the full field equations to O(H0). If the correct equation in these variables is D'' + 3H D' - a^{-2}∇²D = 0 for the strain D entering Eq. (21), recompute Eq. (23) with h replaced by h/a² (equivalently, replace the amplitude factor 1+H0(...) by 1-H0(...) in the integrand) and compare the resulting peak angle with Eq. (25). A shift of the peak beyond the accuracy claimed in Table 1 would overturn the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every observable prediction (Eqs. 18, 23, 25) flows from the O(H0) wave equation (15): the sign and coefficient of the H0∂_T term turn a flat-space harmonic into the growing amplitude (1+H0T) and into the phase Θ in Eq. (24) whose stationary point yields Eq. (25). But Eq. (15) is not derived. Eq. (9) is introduced with \"one gets\", and Eq. (15) is obtained by dropping the source and the H^2, ä/a terms; the surviving -H0∂_T term is never traced to the linearized Ricci or Einstein equations. This matters because the standard FLRW tensor-mode equation for the dimensionless comoving strain D, the quantity that enters the timing-residual formula (21), is D'' + 3H D' - a^{-2}∇²D = 0, with +3H, not -H. The paper's h could be either D or δg = a²D; only the latter obeys Eq. (15), but then it is not the strain that should be inserted in Eq. (21) without an additional a^{-2} factor. Eq. (23) uses the (1+H0T) growth directly. The paper never states which convention h obeys, so the chain Eq. (15) → Eq. (18) → Eq. (23) is unsupported at its most load-bearing point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript claims that gravitational waves propagating in a Friedmann-Lemaitre-Robertson-Walker (FLRW) background have, at first order in the Hubble constant H0, an effective wavenumber that differs from the effective frequency, and that this difference produces a sharp enhancement in pulsar timing residuals for a specific angle between the source and the pulsar as seen from Earth. The paper rederives this effect by linearizing the Einstein equations around FLRW, obtains an approximate wave equation in Eq. (15), and gives a wave solution in Eq. (18) that is also obtained by transforming a Schwarzschild-de Sitter wave to comoving coordinates. It then computes the timing residual in Eq. (23), verifies the approximate enhancement relation Eq. (25) numerically in Table 1, studies the dependence of the signal on frequency, amplitude, source distance, and pulsar distance, and proposes a pulsar timing array search strategy based on 'enhancement rings.' The central physical claim is that the angular position of the enhancement can be used as a local probe of H0.","tokens_in":19888,"tokens_out":21746,"duration_ms":217666,"significance":"If the central wave solution and the resulting phase formula were correct, the paper would present a genuinely interesting possibility: a single-source pulsar timing array measurement of H0 from a local (z much less than 1) geometry, with a concrete observational signature that is distinctive and difficult to mimic by changing the strain amplitude or frequency. The manuscript is clearly written, uses realistic pulsar distances, and provides a useful parametric study of the proposed signal. However, the significance is conditional: the entire chain of predictions depends on the linearized wave equation Eq. (15), whose derivation is not shown and which appears inconsistent with the standard FLRW tensor-mode equation. The numerical checks in Section 5.1 are internal consistency checks rather than independent validations. Thus the paper is not yet in a publishable state; the central derivation must be established or corrected before the observational proposal can be assessed.","major_comments":[{"comment":"The wave equation is not derived. Equation (15) is introduced after 'one gets' and after dropping higher-order terms, but the paper never states what metric variable h_ij denotes. In the standard FLRW transverse-traceless decomposition, the dimensionless comoving strain obeys h_ddot + 3H h_dot - a^{-2} grad^2 h = 0, whereas Eq. (15) has a single term -H0 h_T. The sign and coefficient of this Hubble-friction term control the (1+H0 T) amplitude growth and the phase Theta in Eq. (24), so the chain Eq. (15) -> Eq. (18) -> Eq. (23) -> Eq. (25) is unsupported at its most load-bearing point.","section":"Section 2, Eq. (15)"},{"comment":"The passage from Eq. (9) to Eq. (15) is also internally unclear. The radial operator in Eq. (9) as printed is -R^2/(2a^2) d/dR(1/R^2 d/dR), which expands to -1/(2a^2) partial_R^2 + 1/(a^2 R) partial_R, not the operator -1/a^2(partial_R^2 + 2/R partial_R) that appears in Eq. (15). Because Eq. (15) is the equation that the wave solution Eq. (18) is claimed to satisfy, this mismatch needs to be displayed and resolved explicitly.","section":"Section 2, Eqs. (9) and (15)"},{"comment":"The numerical verification compares the integral in Eq. (23), which is built from the assumed wave form Eq. (18), with the approximate relation Eq. (25) taken from Ref. [21]. Both expressions derive from the same assumed phase structure, so Table 1 demonstrates self-consistency between an approximate stationary-phase formula and a numerical integration of the same model, but it does not validate Eq. (15) against independent physics. The Conclusion's statement that the effect is 'firmly established' therefore overstates what has been shown in this manuscript.","section":"Section 5.1 and Table 1"}],"minor_comments":[{"comment":"The symbol Delta in Eq. (17) is not defined; it should be specified (presumably Delta = 1) or removed, since the coordinate transformation is central to the derivation of Eq. (18).","section":"Section 3, Eq. (17)"},{"comment":"There are typographical slips in the frequency discussion, e.g. 'w = 1, 3 nHz' should presumably read 'w = 1.3 nHz'.","section":"Section 5.2"},{"comment":"The proposed observational protocol does not include a quantitative signal-to-noise estimate or an expected source-rate calculation; the paper appropriately notes that integration times are not addressed, but the feasibility claims in the conclusions should be softened accordingly.","section":"Sections 5 and 6"}],"recommendation":"reject","confidential_remarks":"The manuscript is part of a series (Refs. [13-15,20,21]) and depends heavily on prior work for the coordinate transformation and for Eq. (25). The refereeing outcome hinges on whether Eq. (15) can be derived within the standard linearized FLRW tensor formalism; if the conventional +3H friction term is the correct one, the claimed enhancement is likely an artifact of the chosen variable. The paper should not be considered for publication until this is resolved. I therefore recommend rejection of the present version, while noting that a corrected derivation, if it can be supplied, could change that assessment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a re-derivation of an effect already claimed by the same group, not a new discovery. What's new is a cleaner wave-equation route, a full ΛCDM background, and a systematic parameter study. The internal arithmetic is mostly sound: (18) does satisfy (15) to O(H0), and Table 1's match between the peak formula (25) and the numerical integration is reassuring. The observational protocol—three pulsars in an enhancement ring, using real IPTA pulsars—is concrete and useful.\n\nThe soft spot is exactly where the reader placed it. Eq. (15) is load-bearing: it produces the (1+H0T) amplitude and the k_eff ≠ w_eff phase. But it is introduced with \"one gets\" after dropping terms, and it is never reconciled with the standard FLRW tensor-mode equation. If h is the dimensionless comoving strain, the standard equation has +3H∂_T, not −H∂_T. If h is δg = a²D, then you can't insert it into (21) without an extra a^{-2} factor. The paper doesn't say which h it is, and every observable prediction flows from that choice. This is a genuine gap, not a manufactured one.\n\nA related weakness: the verification in Sec. 5.1 is a self-consistency check, not an independent test. The approximate formula (25) and the wave solution (18) come from the same coordinate transformation and the same prior papers [13,15,21]. No external derivation or simulation confirms the effect. The paper also gives no noise budget for the proposed three-pulsar detection; it says integration time was not considered. For a proposal whose selling point is observability, that is a real omission, though minor relative to the equation issue.\n\nThe citation pattern is insular but not dishonest. The central prior results are all from the same collaboration, which would be fine if the derivation were solid; the problem is that it isn't independently checkable from this paper.\n\nWho is this for? Someone working on PTA single-source searches who wants a concrete, falsifiable target could use this as a pointer. But nobody should build a search on it until the wave-equation issue is settled. My recommendation: send it to a referee, not desk reject. The paper is clear, internally consistent, and potentially important, but the referee should press hard on the definition of h and the derivation of Eq. (15), or ask for an explicit comparison with the standard FLRW tensor-mode calculation.","headline":"An internally consistent but self-referential proposal for an H0-dependent PTA enhancement; the central wave equation is asserted, not derived, and may not be the right equation for the timing residual.","tokens_in":20501,"tokens_out":1890,"would_cite":false,"duration_ms":20543,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83F05"],"pacs":["04.30.-w","98.80.-k"],"model":"deepseek-v4-flash","headline":"In an expanding universe, the effective wavenumber of a gravitational wave differs from its redshifted frequency, and the paper argues that this mismatch makes pulsar timing residuals peak at an angle determined by H0.","keywords":["pulsar timing arrays","Hubble constant","gravitational waves","FLRW metric","effective wavenumber","timing residual enhancement","local H0 measurement","cosmological constant"],"falsifier":"Substitute the proposed solution (19) into the standard linearized tensor-mode equation on an FLRW background, $h'' + 3H h' - a^{-2}\\nabla^2 h = 0$: if it is not a solution at $\\mathcal{O}(H_0)$, the claimed wavenumber mismatch is an artifact of Eq. (15). Observationally, for a known source at distance $Z_E$, the enhanced-residual ring must appear at the angle given by $H_0 \\simeq (2c/Z_E)\\sin^2(\\alpha/2)$; a deep pulsar-timing search around that angle for a known single source would settle the claim.","tokens_in":19301,"feed_emoji":"📡","tokens_out":8545,"duration_ms":76283,"temperature":0.7,"pith_summary":"This paper claims that the standard treatment of gravitational waves in cosmology, which only red-shifts the frequency, misses a first-order effect of the Hubble constant. In FLRW coordinates the physical wave has an effective wavenumber $k_{\\rm eff} = w(1 - H_0 R/2)$ that differs from the effective frequency $w_{\\rm eff} = w(1 - H_0 R)$. When a pulsar timing residual integrates the wave along the line of sight, the mismatch produces a sharp enhancement at an angle satisfying $H_0 \\simeq (2c/Z_E) \\sin^2(\\alpha/2)$, so the angular position of the peak is a local measure of $H_0$. The paper rederives this result from the linearized wave equation, extends it to a universe with dust, radiation, and a cosmological constant, and shows the peak survives changes in frequency, amplitude, and pulsar distance. If correct, this could give pulsar timing arrays both an easier route to single-source detection and a new local measurement of the Hubble constant.","feed_headline":"Pulsar timing may read H0 from a single sky angle","feed_subtitle":"A gravitational-wave effect concentrates the timing signal at an angle set by the Hubble constant and the source distance.","key_machinery":"The load-bearing object is the first-order wave solution $h_{ij} = (e_{ij}/R)(1+H_0 T)\\cos(w_{\\rm eff} T - k_{\\rm eff} R)$, together with the $\\mathcal{O}(H_0)$ coordinate transformation $t = T + R^2 H_0/2$, $r = R(1+\\Delta T H_0)$ that converts the source-frame harmonic (16) into this comoving form. The physical content is the difference between $w_{\\rm eff}$ and $k_{\\rm eff}$: the wave's phase velocity is no longer unity in comoving coordinates. The same solution is also obtained from the linearized FLRW wave equation (15), which at $\\mathcal{O}(H_0)$ contains the single extra term $-H_0 \\partial_T$. The timing residual integral (22)--(24) carries that phase mismatch along the pulsar--Earth path, and the enhancement relation (25) follows from the stationary-phase condition of that integral.","core_discovery":"The central claim, stated in the paper's own terms, is that the metric perturbation of a gravitational wave in FLRW comoving coordinates, to first order in $H_0$, is $h_{ij} = (e_{ij}/R)(1+H_0 T)\\cos(w_{\\rm eff} T - k_{\\rm eff} R)$ with $w_{\\rm eff} = w(1-H_0 R)$ and $k_{\\rm eff} = w(1-H_0 R/2)$. A naive harmonic wave with only a redshifted frequency is not a solution of the linearized equation; the anharmonicity induced by the coordinate transformation between the Schwarzschild--de Sitter source frame and the FLRW observer frame is needed. Because $k_{\\rm eff} \\neq w_{\\rm eff}$, local experiments see only the usual frequency shift, but a pulsar timing residual, which accumulates over a long null geodesic, does not cancel, and its magnitude peaks sharply when the angle between source and pulsar satisfies $H_0 \\simeq (2c/Z_E)\\sin^2(\\alpha/2)$. This is the mechanism by which the paper proposes to measure $H_0$ locally.","pith_inferences":["A reader should not assume Eq. (15) is the standard FLRW tensor-mode equation: the usual comoving perturbation equation has a $+3H\\partial_T$ friction term, while Eq. (15) has $-H_0 \\partial_T$. Whether the proposed solution survives the standard equation is a check the paper leaves for follow-up work.","If the wavenumber mismatch is a gauge artifact of the particular coordinate mapping, the angular peak would not persist in a fully gauge-invariant treatment; testing Eq. (19) in a gauge-invariant formulation would settle whether the predicted ring is physical.","A practical cross-check with existing IPTA data would be to look for simultaneous anomalous residuals in pulsars whose mutual separations match a single enhancement ring for a known candidate supermassive-black-hole merger; the absence of such a ring would count against the effect.","The relation between ring width and pulsar distance (narrower peaks for farther pulsars) could be used to confirm the effect statistically: all pulsars in the ring should show peaks centered at the same $\\alpha_{\\rm max}$ but with widths that follow the paper's FWHM dependence."],"forward_implications":["A single supermassive-black-hole merger observed by a pulsar timing array would produce a ring-shaped enhancement region on the sky; any pulsar inside the ring should show a residual many times larger than the rest of the signal.","The angular radius of that ring, with the source distance known, yields a local value of $H_0$ at redshift well below 1, independent of the usual distance-ladder and CMB calibrations.","Because the peak position is insensitive to the wave's frequency, amplitude, and polarization, the effect cannot be mimicked by changing source parameters.","At first order all cosmological components enter only through $H_0$, so a single measurement gives the Hubble constant rather than the separate densities; the paper notes next-order corrections could separate the components.","Stochastic-background searches that average over many pulsars and assume decorrelated signals would wash out this single-source effect, so PTA detection strategies may need to look for individual events in triplets of pulsars."],"supporting_citations":[{"why":"Original derivation of the enhancement in pulsar timing residuals from propagation in a curved metric; this work rederives and extends it.","marker":"[13]"},{"why":"Provides the coordinate transformation between Schwarzschild--de Sitter and FLRW coordinates at first order, including all cosmological components.","marker":"[15]"},{"why":"Earlier treatment of gravitational waves in the presence of a cosmological constant, which motivates the source-frame wave.","marker":"[20]"},{"why":"Independent derivation of the relation between H0 and the enhancement angle, used here as a cross-check on Eq. (25).","marker":"[21]"},{"why":"Supplies the timing-residual integral and null-geodesic parametrization used to compute the gravitational-wave timing residual.","marker":"[22]"},{"why":"First IPTA data release providing the real pulsar distances used in the numerical signal characterization.","marker":"[27]"}],"fun_headline_variants":["Gravitational waves morph with H0, pulsar timing sees angle","Pulsar timing boost reveals H0 from source angle","Measure H0 with a single pulsar-source angle","GW signal amplifies at H0-dependent angle in PTAs","H0 from pulsar timing: the angle that amplifies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the $\\mathcal{O}(H_0)$ coordinate transformation between the Schwarzschild--de Sitter source frame and FLRW comoving coordinates, and the linearized wave equation (15) with its single $-H_0\\partial_T$ term, correctly describe the physical gravitational-wave metric perturbation; the paper does not reconcile this equation with the standard FLRW tensor-mode equation, which contains a $+3H\\partial_T$ friction term.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational waves morph with H0, pulsar timing sees angle","Pulsar timing boost reveals H0 from source angle","Measure H0 with a single pulsar-source angle","GW signal amplifies at H0-dependent angle in PTAs","H0 from pulsar timing: the angle that amplifies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1714,"prompt_tokens":992,"completion_tokens":722,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":637}},"tokens_in":608,"tokens_out":722,"duration_ms":6404,"temperature":1.0,"reasoning_tokens":637,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:41:44.668164+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Substitute the proposed solution (19) into the standard linearized tensor-mode equation on an FLRW background, $h'' + 3H h' - a^{-2}\\nabla^2 h = 0$: if it is not a solution at $\\mathcal{O}(H_0)$, the claimed wavenumber mismatch is an artifact of Eq. (15). Observationally, for a known source at distance $Z_E$, the enhanced-residual ring must appear at the angle given by $H_0 \\simeq (2c/Z_E)\\sin^2(\\alpha/2)$; a deep pulsar-timing search around that angle for a known single source would settle the claim.","supporting_citations":[{"cited_title":"Local measurement of {\\Lambda} using pulsar timing arrays","cited_arxiv_id":"1209.3724","evidence_quote":"Original derivation of the enhancement in pulsar timing residuals from propagation in a curved metric; this work rederives and extends it."},{"cited_title":"On the propagation of gravitational waves in a $\\Lambda$CDM universe","cited_arxiv_id":"1711.08315","evidence_quote":"Provides the coordinate transformation between Schwarzschild--de Sitter and FLRW coordinates at first order, including all cosmological components."},{"cited_title":"Gravitational waves in the presence of a cosmological constant","cited_arxiv_id":"1106.4511","evidence_quote":"Earlier treatment of gravitational waves in the presence of a cosmological constant, which motivates the source-frame wave."}],"review_version":1}