{"id":"ad7eb2c2-c3b7-4772-bff8-468fe40711f6","arxiv_id":"1909.02735","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Co-propagating neutrinos and gravitational waves generate an energy-dependent phase shift that persists to high redshift.","lead":"This paper calculates how spacetime ripples from gravitational waves or density fluctuations nudge the oscillation phase of neutrinos as they travel. The effect is small, but it grows with energy and does not fade for very distant sources.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (16) drops the source-term phase: for co-propagating neutrinos the GW-induced oscillation phase cancels to O(ω m^2 L/E^2), so the claimed E^3 enhancement is an artifact.","rationale":"The reader's weakest assumption focused on the eikonal validity when α+ becomes large in the co-propagating limit. The concern raised here is different and more fundamental: even when the eikonal solution is valid, the quantity computed in Eq. (16) is not the oscillation phase. The oscillation phase is the difference of the action evaluated along the neutrino worldline from source to detector. In the standard plane-wave treatment, the phase at the source is subtracted implicitly by setting the source at the coordinate origin. For the gravitational-wave solution (12), the oscillatory term α+ e^{iω(x3-t)} has a nonzero value at the source, and the mass-dependent part of that value does not cancel from the observable phase. Along a co-propagating trajectory the exponent x3-t is nearly constant, so the source and detector terms almost cancel; the residual is of order ωL(E-K)/E = ωL(k⊥²+m²)/(2E²). For the numbers used in the paper this residual is about 10⁻¹⁶, meaning the claimed O(1) phase is suppressed by sixteen orders of magnitude. This directly invalidates the central claim of an energy-growing gravitational contribution to the oscillation phase. It also affects the FRW persistence claim, since the local value of α̃ at the detector cannot be the accumulated path phase. The paper's own statement that the phase difference is zero for strictly co-propagating same-source neutrinos is consistent with this cancellation but is not applied to the off-axis case, where k⊥≠0 only changes the suppression factor, not the cancellation. Because the main novel result depends on this missing subtraction, the central claim does not stand as presented. A single endpoint-subtraction calculation settles the issue cleanly.","tokens_in":9092,"tokens_out":38453,"duration_ms":408295,"concrete_test":"Re-evaluate the phase difference for the two-detector setup using ΔS = [S1(x2,t2)-S1(x1,t1)] - [S2(x2,t2)-S2(x1,t1)] with S as in Eq. (12), source at (t1,x1), detector at (t2,x2). Confirm the result is δα(e^{iω(L-T)}-1), not δα; then insert the paper's values (X=10³ km, L=10⁸ pc, ω=100 Hz, E=10⁵ TeV, m≈1 eV) and compute |1-e^{iωL(K/E-1)}| ≈ ωL m²/(2E²). If this factor is ~10⁻¹⁶, the claimed phase of order unity is suppressed to <10⁻¹³, falsifying Eq. (16) as the observable oscillation phase.","verdict_should_be":"REJECT","load_bearing_attack":"Equation (16) identifies the mass-dependent coefficient α+ of the oscillating term in S (Eq. (12)) with the observable oscillation phase. In the WKB/plane-wave treatment the oscillation phase is the action difference over the path, ΔS = [S1(x2,t2)-S1(x1,t1)] - [S2(x2,t2)-S2(x1,t1)], not S1(det)-S2(det). Including the α terms at both endpoints multiplies the amplitude by (e^{iω(x3,2-t2)} - e^{iω(x3,1-t1)}). Along the co-propagating neutrino trajectory with K≈E-(k⊥²+m²)/(2E), this factor has magnitude |1-e^{iωL(K/E-1)}| ≈ ωL(k⊥²+m²)/(2E²). With the paper's quoted parameters (L≈10⁸ pc, ω≈100 Hz, E≈10⁵ TeV, m²≈1 eV², k⊥≈EX/L), this factor is ~10⁻¹⁶. Hence the O(1) phase claimed below Eq. (16) becomes ≲10⁻¹³, and the E³/ω enhancement is an artifact of neglecting the phase at the source. The paper itself notes the same-source phase difference is zero; the same cancellation applies to the off-axis detector because the GW phase is nearly constant along the trajectory. The FRW persistence claim is affected likewise: the local value of α̃ at the detector is not the accumulated path phase. The reader's concern about eikonal breakdown is secondary: even if α+ is small, the missing endpoint subtraction removes the enhancement.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the eikonal (WKB) approximation to the Klein-Gordon/Dirac equation in a curved background to compute corrections to the neutrino oscillation phase induced by tensor (gravitational-wave) and scalar metric perturbations, in both flat spacetime and an FRW background. In the flat-spacetime tensor case, for neutrinos propagating nearly parallel to a plane gravitational wave with a small transverse momentum, the authors obtain a mass-dependent coefficient alpha_+ (Eq. (15)) and estimate an additional oscillation phase delta-phi ~ (E^3/omega)(X^2/L^2)(delta-m^2/m^4) A_gw (Eq. (16)), which they claim can reach order unity for E ~ 10^5 TeV. In the FRW background, analogous coefficients are derived for tensor perturbations (Eq. (36)) and scalar perturbations (Eq. (40)), and the authors claim that, unlike the unperturbed contributions, the perturbation-induced phase differences persist at large redshifts. The paper also discusses the phenomenological difficulties of observing such effects, including the need for coincident gravitational-wave and neutrino sources and the coherence bound of Eq. (17).","tokens_in":9355,"tokens_out":22207,"duration_ms":231668,"significance":"The eikonal formalism in Sec. 2 is standard and the algebraic derivation of the coefficients alpha_+ and alpha~_+ is straightforward and internally consistent; the paper is clearly written and openly states its assumptions. If the central claim were correct, the E^3/omega growth would be a novel and phenomenologically interesting signature, and the high-redshift persistence would distinguish the perturbation contributions from the constant FRW terms. However, the main quantitative claim is vitiated by an incorrect identification of the eikonal coefficient with the observable oscillation phase: the source-to-detector action difference, not the value of the phase at the detector, determines the oscillation phase. This error removes the E^3/omega enhancement in the tensor case and undercuts the FRW tensor persistence claim. The algebraic results in Eqs. (15) and (36) may still be of technical interest, but the central physical conclusion does not follow as presented.","major_comments":[{"comment":"The identification of the coefficient alpha_+ in Eq. (15) with the observable oscillation phase is inconsistent with the standard WKB definition of the phase difference. The oscillation phase is the difference of actions accumulated from a common source to the detector, delta-phi = -[(S_1(det)-S_1(src))-(S_2(det)-S_2(src))], not the value of S at the detector with the source phase set to zero. With the ansatz Eq. (12) and the source at (t_1,x_31)=(0,0), the gravitational-wave contribution to each eigenstate is alpha_i(e^{i omega (L-T)}-1). Along the trajectory defined by Eq. (14), L-T roughly equals (k_perp^2+m^2)L/(2E^2), so this factor is of order omega (k_perp^2+m^2)L/(2E^2); for the quoted parameters (omega ~ 100 Hz, L ~ 10^8 pc, E ~ 10^5 TeV, m^2 ~ 1 eV^2, k_perp ~ EX/L) this is ~10^-16. The O(1) phase claimed below Eq. (16) is therefore an artifact of dropping the source-term phase, and the E^3/omega enhancement in Eq. (16) does not follow. The paper's own remark that the phase difference is zero when the neutrino beam and the gravitational wave come from the same source is the leading instance of this cancellation; the cancellation persists for the off-axis detector because the gravitational-wave phase is nearly constant along the ultra-relativistic trajectory.","section":"Sec. 2, Eqs. (12)-(16)"},{"comment":"The claim that the tensor-perturbation contribution to the oscillation phase 'does not tend to a constant for large z' is based on the observation that the amplitude alpha~_+ in Eq. (36) does not vary with the scale factor. This is not the relevant quantity: the accumulated phase difference again requires subtraction of the source term, alpha~_i(e^{i omega (x_3-eta)}|_{det} - e^{i omega (x_3-eta)}|_{src}). Along the highly relativistic neutrino trajectory in a flat FRW background, eta is approximately equal to x_3, so the same near-cancellation as in the flat-space case occurs, and the FRW tensor correction is suppressed by a similarly small factor of order omega (m~^2+k_perp^2)L/(2E_0^2). No calculation of the integrated phase difference is given, so the persistence claim for tensor perturbations is unsupported.","section":"Sec. 3, Eqs. (34)-(36) and the following discussion"},{"comment":"The scalar-perturbation phase difference in Eq. (40) is complex-valued, but the oscillation probability in Eq. (5) requires a real phase. As written, the expression contains an explicit factor 1/i, and the second panel of Fig. 1 adds this complex number to the real function F_2(z). The authors need to specify that the real part is taken, and to justify discarding the imaginary part, which would correspond to an amplitude modulation rather than a phase modification. Without this specification, the displayed result is not well defined.","section":"Sec. 3, Eq. (40) and Fig. 1"}],"minor_comments":[{"comment":"The matrix display for g_ij is corrupted by a stray period after the last row; please fix the typesetting.","section":"Eq. (10)"},{"comment":"The dimensional statement 'delta-m^2 ~ m^2 ~ 10^{-w} eV' is inconsistent, since m^2 has units of energy squared; the text should specify the mass-squared units clearly.","section":"After Eq. (16)"},{"comment":"The inequality 'K >> m^2' should read 'K^2 >> m^2' to be consistent with the surrounding expansions.","section":"After Eq. (28)"},{"comment":"The step leading from Eq. (15) to Eq. (16) is not shown; the relation between alpha_+ and the oscillation phase, and the approximations used for k_1, k_+, and the detector geometry, should be written out explicitly.","section":"Eq. (15) to Eq. (16)"},{"comment":"The statement that the tensor-perturbation contribution 'does not tend to a constant for large z' is asserted without a direct calculation of the z-dependence of the endpoint-subtracted expression; the discussion should be tied to an explicit integral or estimate.","section":"Sec. 3, tensor discussion"},{"comment":"There are numerous spelling and typographical errors ('FR W', 'Schwarchild', 'interestng', 'Wl', 'measuments'); the manuscript needs a careful proofreading pass.","section":"General"}],"recommendation":"reject","confidential_remarks":"The central error is mathematical, not a matter of taste or consensus: the observable oscillation phase is the action difference between the two mass eigenstates over the path, and the missing endpoint subtraction in the tensor case suppresses the claimed E^3/omega enhancement by a factor of order 10^-16 for the quoted parameters. This is a load-bearing flaw that invalidates the main claim of the paper and the FRW tensor persistence result. A corrected tensor calculation would likely reduce the effect to a negligible level, leaving only the scalar-perturbation result, which the manuscript itself describes as too small for phenomenological relevance. I therefore recommend rejection, though a future revision that properly computes the integrated phase difference and focuses on the scalar case could be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nI read this paper with the stress-test note in hand, and I think the note is right: the claimed E^3/omega enhancement in Eq. (16) is an artifact of dropping the source term. The authors compute the mass-dependent coefficient alpha_+ of the oscillating phase and treat its value at the detector as the oscillation phase. But the physical phase is the difference in action along the path, and the alpha contribution enters as (alpha_1 - alpha_2)(e^{i phi(L)} - e^{i phi(0)}). Along a co-propagating trajectory e^{i phi(L)} - e^{i phi(0)} is roughly omega (k_perp^2 + m^2) L / (2 E^2), which for their benchmark numbers is about 10^-16. Multiply that into delta_alpha and the order-unity phase becomes 10^-13 or smaller. The E^3/omega growth is canceled by the same E^2 in the bracket. The same issue affects the FRW persistence claim: the local alpha_tilde at the detector is not the accumulated path phase.\n\nWhat the paper does well: the eikonal setup with the Lichnerowicz curvature term is careful, the tensor-perturbation algebra through Eqs. (15) and (36) checks out, and the authors are transparent about the two-source coincidence problem. The FRW extension is a natural step and the separation between the F1/F2 terms and the perturbative corrections is clearly organized.\n\nSoft spots beyond the main one: the scalar-perturbation phase in Eq. (40) is complex and never real-averaged; the figure inflates c_s to 10^-3 and still shows only a modest effect; and the eikonal consistency for alpha_+ is not quantified, though this is secondary -- omega alpha is small even when alpha itself is large.\n\nBottom line: the algebraic machinery is useful as a formal exercise, but the phenomenological claim does not survive. I would not cite the numerical result. A referee should catch the endpoint issue quickly. The paper deserves a serious referee rather than a desk reject because the derivation is clean and the mistake is instructive. I would send it out, expecting either a major revision or removal of the central claim.\n\nBest,\n[Your name]","headline":"Eq. (16)'s E^3 enhancement evaporates once you subtract the source-phase term; the eikonal algebra is sound but the central numerical claim is an artifact.","tokens_in":9917,"tokens_out":12205,"would_cite":false,"duration_ms":116705,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that neutrino flavor oscillations acquire an additional phase from gravitational wave and scalar perturbation backgrounds, with the extra phase growing like $E^3/\\omega$ and persisting to high redshift.","keywords":["neutrino oscillations","eikonal approximation","gravitational waves","FRW spacetime","scalar perturbations","curved spacetime neutrino propagation","ultra-high-energy neutrinos","phase difference"],"falsifier":"A direct numerical solution of the Klein-Gordon or Dirac equation in the plane gravitational wave background, across the parameter range used for Eq. (16), would settle whether the exact oscillation phase follows the predicted $E^3/\\omega$ growth or departs from it once $A_{gw}E^2/(\\omega m^2)$ approaches unity.","tokens_in":8839,"feed_emoji":"🌊","tokens_out":8636,"duration_ms":83225,"temperature":0.7,"pith_summary":"The paper sets out to show that neutrino flavor oscillations pick up an additional phase when the neutrinos propagate through gravitational wave or scalar perturbation backgrounds, both in flat spacetime and in an expanding FRW universe. Using an eikonal approximation, it finds that the extra phase difference between two mass eigenstates scales as $E^3/\\omega$ for neutrinos traveling with a gravitational wave, along with geometric and mass factors, rather than the usual $1/E$ dependence. In cosmological backgrounds, the perturbation-induced phase does not tend to a constant at large redshift, unlike the standard FRW contributions, so the effect would persist for sources at high $z$. The authors estimate that for typical astrophysical parameters the phase could reach order unity for ultra-high-energy neutrinos, although a gravitational wave burst and a neutrino signal would have to coincide in a very small time window.","feed_headline":"Gravitational waves add an energy-cubed phase to neutrino oscillations","feed_subtitle":"The extra phase survives to high redshift, unlike the standard cosmological phase, and can reach order one for extreme neutrino energies.","key_machinery":"The central object is the eikonal ansatz with an oscillating correction, $S=-Et+k_1x_1+k_2x_2+Kx_3+\\alpha_+e^{i\\omega(x_3-t)}$ (and its FRW counterpart), inserted into the generalized eikonal equation $g^{\\mu\\nu}\\partial_\\mu S\\,\\partial_\\nu S=m^2-R/4$. The decisive step is that for neutrinos moving in the same direction as the gravitational wave, $E-K\\approx(k_\\perp^2+m^2)/(2E)$ is very small, so the denominator in the equation for $\\alpha_+$ is tiny and the correction grows with neutrino energy. In FRW backgrounds the perturbation amplitudes carry a factor $1/a$, which cancels the scale-factor dependence and keeps the extra phase from saturating at large redshift. This ansatz carries the argument because it turns a metric perturbation directly into a definite phase coefficient.","core_discovery":"The central claim is that the phase $S$ of a neutrino propagating in a gravitational wave background can be written as the flat-space phase plus a small oscillating term $\\alpha_+ e^{i\\omega(x_3-t)}$, and the eikonal equation fixes $\\alpha_+ \\propto c_+ k_+^2 E/[\\omega(k_\\perp^2+m^2)]$. For two detectors separated by a transverse distance $X$, this yields an extra oscillation phase difference $\\delta\\tilde{\\phi}_{12} \\sim (E^3/\\omega)(X^2/L^2)(\\delta m^2/m^4)A_{gw}$. In FRW spacetime, with tensor perturbation $h_+=(c_+/a)e^{i\\omega(x_3-\\eta)}$ and scalar perturbation $\\Phi=c_s e^{ikx_3}$, the oscillating phase amplitude does not scale with the scale factor, so the additional phase difference does not tend to a constant at large redshift, unlike the standard $F_1$ and $F_2$ contributions. The result is therefore largest for high-energy neutrinos co-propagating with the wave and persists to large redshifts within the validity of the approximation.","pith_inferences":["If the $E^3/\\omega$ scaling survives a full non-perturbative treatment, gravitational wave backgrounds should imprint an energy-dependent oscillation signature distinct from standard matter effects; this testable distinction is not developed in the paper.","The same eikonal machinery applied to a stochastic gravitational wave background would have to average over propagation directions, which should erase the co-propagating enhancement; the paper notes an earlier average (reference [10]) but does not analyze the directional dependence itself.","Because the paper does not check whether $A_{gw}E^2/(\\omega m^2)$ remains small in the co-propagating regime, the natural next step is a numerical or higher-order eikonal solution to see whether the phase saturates or changes form where the perturbative estimate would break.","If the persistence at high redshift is correct, primordial gravitational wave or scalar perturbation backgrounds could in principle leave a redshift-dependent oscillation phase that distinguishes them from standard FRW contributions."],"forward_implications":["The gravitational contribution to the oscillation phase grows with neutrino energy rather than falling like the standard $1/E$ term, so the highest-energy neutrinos receive the largest correction.","In FRW backgrounds the perturbation-induced phase does not approach a constant at large redshift, so high-redshift sources remain distinguishable by the perturbed oscillation phase.","The enhancement appears only for neutrinos co-propagating with the gravitational wave; neutrinos moving opposite to the wave do not get the $E^3/\\omega$ amplification.","Order-of-magnitude estimates place the effect within the reach of ultra-high-energy neutrino telescopes for $E\\sim 10^5$ TeV, although a gravitational wave burst and a neutrino burst must coincide in a small time window."],"supporting_citations":[{"why":"Supplies the standard neutrino oscillation probability and the baseline phase difference $\\delta m^2 L/(2E)$ that the new terms modify.","marker":"[1]"},{"why":"Supplies the curved-spacetime phase-difference results, including the Schwarzschild correction to which the new flat-space estimate is compared.","marker":"[2]"},{"why":"Treats a similar gravitational-wave-induced neutrino phase problem with a different geometry, providing the contrast for the enhancement found here.","marker":"[5]"},{"why":"Gives the unperturbed FRW phase-difference functions $F_1,F_2$ that tend to constants at large redshift, the behavior the perturbations are shown to change.","marker":"[6]"},{"why":"Provides the earlier solution of the eikonal equation in a gravitational wave background that the present ansatz extends.","marker":"[9]"},{"why":"Averages the gravitational wave effect over a stochastic background and therefore misses the co-propagating enhancement, serving as the comparison case.","marker":"[10]"},{"why":"Supplies the neutrino mass-squared values used in the order-of-magnitude estimates for the extra phase.","marker":"[12]"},{"why":"Documents the energy reach of ultra-high-energy neutrino telescopes used to assess observability of the phase.","marker":"[13]"},{"why":"Gives the coherence-length bound used to check that the astrophysical baseline does not decohere the beam.","marker":"[14]"}],"fun_headline_variants":["Gravitational waves add energy-cubed neutrino phase","GWs shift neutrino oscillations with energy-cubed phase","High-energy neutrinos get GW phase at every redshift","Gravitational wave phase persists to high redshift for neutrinos","Neutrino oscillation phase from GWs survives redshift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the gravitational wave modifies the neutrino phase by only a small oscillating term, so the first-order eikonal equation applies; in the co-propagating geometry that produces the enhancement, this requires $A_{gw}E^2/(\\omega m^2)$ to remain small, and if it does not, the predicted phase scaling fails.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational waves add energy-cubed neutrino phase","GWs shift neutrino oscillations with energy-cubed phase","High-energy neutrinos get GW phase at every redshift","Gravitational wave phase persists to high redshift for neutrinos","Neutrino oscillation phase from GWs survives redshift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000705,"raw_usage":{"total_tokens":3134,"prompt_tokens":855,"completion_tokens":2279,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":2199}},"tokens_in":471,"tokens_out":2279,"duration_ms":17357,"temperature":1.0,"reasoning_tokens":2199,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:40:58.599964+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct numerical solution of the Klein-Gordon or Dirac equation in the plane gravitational wave background, across the parameter range used for Eq. (16), would settle whether the exact oscillation phase follows the predicted $E^3/\\omega$ growth or departs from it once $A_{gw}E^2/(\\omega m^2)$ approaches unity.","supporting_citations":[{"cited_title":"Pontecorvo, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the standard neutrino oscillation probability and the baseline phase difference $\\delta m^2 L/(2E)$ that the new terms modify."},{"cited_title":"Dvornikov, Phys","cited_arxiv_id":null,"evidence_quote":"Treats a similar gravitational-wave-induced neutrino phase problem with a different geometry, providing the contrast for the enhancement found here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the unperturbed FRW phase-difference functions $F_1,F_2$ that tend to constants at large redshift, the behavior the perturbations are shown to change."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier solution of the eikonal equation in a gravitational wave background that the present ansatz extends."},{"cited_title":"Dvornikov, Phys","cited_arxiv_id":null,"evidence_quote":"Averages the gravitational wave effect over a stochastic background and therefore misses the co-propagating enhancement, serving as the comparison case."},{"cited_title":"Aker et al","cited_arxiv_id":null,"evidence_quote":"Supplies the neutrino mass-squared values used in the order-of-magnitude estimates for the extra phase."},{"cited_title":"Ishihara (IceCube Collaboration), J","cited_arxiv_id":null,"evidence_quote":"Documents the energy reach of ultra-high-energy neutrino telescopes used to assess observability of the phase."},{"cited_title":"Giunti and C","cited_arxiv_id":null,"evidence_quote":"Gives the coherence-length bound used to check that the astrophysical baseline does not decohere the beam."}],"review_version":1}