REVIEW 3 major objections 6 minor 1 cited by
Neutrino oscillations in gravitational and cosmological backgrounds
T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. 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 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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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).
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 (3)
- [Sec. 2, Eqs. (12)-(16)] 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.
- [Sec. 3, Eqs. (34)-(36) and the following discussion] 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.
- [Sec. 3, Eq. (40) and Fig. 1] 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.
minor comments (6)
- [Eq. (10)] The matrix display for g_ij is corrupted by a stray period after the last row; please fix the typesetting.
- [After Eq. (16)] 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.
- [After Eq. (28)] The inequality 'K >> m^2' should read 'K^2 >> m^2' to be consistent with the surrounding expansions.
- [Eq. (15) to Eq. (16)] 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.
- [Sec. 3, tensor discussion] 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.
- [General] There are numerous spelling and typographical errors ('FR W', 'Schwarchild', 'interestng', 'Wl', 'measuments'); the manuscript needs a careful proofreading pass.
Circularity Check
No significant circularity: the phase shifts are derived from the eikonal equation with stated assumptions, with no fitted inputs or load-bearing self-citations.
full rationale
The derivation is self-contained. The additional phase coefficient alpha_+ in Eqs. (12)-(15) is obtained by inserting the stated eikonal ansatz into the eikonal equation (11) and solving algebraically for alpha_+; no parameter is fitted to the target oscillation phase, and Eq. (16) is simply that coefficient evaluated with the paper's stated kinematic approximations (k2 = 0, k1 approximately E X/L, k1^2 << m^2). The FRW extension in Sec. 3 similarly solves Eq. (8) with the perturbed metric and displays the resulting coefficients in Eqs. (35)-(40); the cosmological inputs (H0, Omega_m, Omega_Lambda, z) and neutrino mass/energy values are external data, not outputs being predicted. The only author self-citation (Koutsoumbas et al. 2018 in ref. [16]) appears in a list of future-work directions and carries no weight in the derivation. Possible physical objections, such as whether the coefficient alpha_+ should be interpreted as an integrated path phase or whether the eikonal consistency condition holds for co-propagating neutrinos, concern the correctness of the perturbative calculation rather than circularity. Therefore no circular step is exhibited.
Assumptions & free parameters
assumptions (4)
- domain assumption The eikonal (WKB) approximation applies to the neutrino wavefunction, with slowly varying amplitude and phase, satisfying g^mu nu partial_mu S partial_nu S = m^2 - R/4.
- domain assumption The fermion field obeys the Lichnerowicz identity (gamma^mu nabla_mu)^2 = box - R/4, so the eikonal equation contains the curvature term; other treatments that use a covariantly constant spinor drop it.
- domain assumption Perturbations are small and first-order: h_+ and Phi are much less than 1, so only terms linear in the perturbation are kept.
- domain assumption The FRW background is flat with Omega_r + Omega_m + Omega_Lambda = 1 and standard expansion H^2 = H0^2 (Omega_r/a^4 + Omega_m/a^3 + Omega_Lambda).
Cite this review
Pith. "Pith review of Neutrino oscillations in gravitational and cosmological backgrounds." pith.science (2026). https://pith.science/paper/W64CRY7T
@misc{pith2026190902735,
author = {Pith},
title = {Pith review of: Neutrino oscillations in gravitational and cosmological backgrounds},
year = {2026},
howpublished = {\url{https://pith.science/paper/W64CRY7T}},
note = {Machine review of arXiv:1909.02735}
}
read the original abstract
We use the eikonal approximation in order to calculate the additional phase shift between two neutrino mass eigenstates during their propagation in a background of gravitational wave or scalar perturbations in the flat and the FRW spacetime metric. We comment on the dependence of the results on the characteristics of the perturbations, give some order-of-magnitude estimates, and find that, although small, the resulting phase difference persists for large redshifts, up to the validity of our approximations.
Figures
Forward citations
Cited by 1 Pith paper
-
Effects of gravitational lensing on neutrino oscillation in Hu-Sawicki f(R) gravity
Neutrino oscillation probabilities from lensed paths in Hu-Sawicki f(R) gravity are derived, showing sensitivity to λ and to neutrino mass parameters.
Reference graph
Works this paper leans on
-
[1]
S.Bilenky and B. Pontecorvo, Phys. Rep. 41, 225 (1978). L. Wolfenstein, Phys. Rev. D17, 2369 (1978). B. Kayser, Phys. Rev. D24, 110 (1981). C. Giunti, C. W. Kim and U. W. Lee, Phys. Rev. D44, 3635 (1991). C. Giunti, J. Phys. G 34, 93 (2007)
work page 1978
-
[2]
L. Stodolsky, Gen. Relativ. Gravit. 11, 391 (1979). D. V. Ahluwalia and C. Burgard, Gen. Relativ. Gravit. 28, 1161 (1996). S. I. Godunov and G. S. Pastukhov, Phys. At. Nucl. 74, 302 (2011). T. Bhattacharya, S. Habib and E. Mottola, Phys. Rev. D59, 067301 (1999). C. Y. Cardall and G. M. Fuller, Phys. Rev. D55, 7960 (1997). 12 R. M. Crocker, C. Giunti and D...
arXiv 1979
-
[3]
M. G. Aartsen et al. (IceCube Collaboration), Phys. Rev. Let. 113, 101101 (2014). S. Adrian-Martinez et al. (ANTARES and IceCube and LIGO Scientific and Virgo Collaborations), Phys. Rev D93, 122010 (2016)
work page 2014
-
[4]
E. Kh. Akhmedov, V. A. Rubakov and A. Yu. Smirnov, Phy. Rev. Lett. 81, 1359 (1998). G. C. Branco, R. Gonzalez Felipe and F. R. Joaquim, Rev. Mod. Phys. 84, 515 (2012). S. Blanchet and P. Di Bari, New. J. Phys. 14, 125012 (2012)
work page 1998
- [5]
-
[6]
L.Visinelli, Gen. Relativ. Gravit. 47, 62 (2015)
work page 2015
-
[7]
H. Lawson and M. Michelsohn, Spin Geometry , Princeton University Press (1989)
work page 1989
-
[8]
B. S. DeWitt, Rev. Mod. Phys. 29, 377 (1957). K. S. Cheng, J. Math. Phys. 13, 1723 (1972)
work page 1957
Show all 16 references
-
[9]
N. J. Poplawski, J. Math. Phys. 47, 072501 (2006)
2006
-
[10]
Dvornikov, Phys
M. Dvornikov, Phys. Rev. D99, 116021 (2019)
2019
-
[11]
Meszaros, D
P. Meszaros, D. B. Fox and K. Murase, Nature Rev. Phys. 1, 585 (2019). C. A. Arguelles, M. Bustamante, A. Kheirandish, S. Palomares- Ruiz, J. Salvado and A. C. Vincent, PoS ICRC2019 , 849 (2020). arXiv:1907.08690 [astro-ph.HE]
2019 arXiv
-
[12]
Aker et al
M. Aker et al. (KATRIN Collaboration), it Phys. Rev. Lett. 123, 221802 (2019). 13
2019
-
[13]
Ishihara (IceCube Collaboration), J
A. Ishihara (IceCube Collaboration), J. Phys. Conf. Ser. 718, 062027 (2016). I. Tamborra, J. Phys. Conf. Ser. 718, 052041 (2016). M. Ageron (ANTARES Collaboration), Nucl. Instrum. Meth. A 656, 11 (2011). A. Karle, J. Phys. Conf. Ser. 39, 379 (2006)
2016
-
[14]
Giunti and C
C. Giunti and C. W. Kim, Phys. Rev. D58, 017301 (1998). M. Beuthe, Phys. Rev. D66, 013003 (2002)
1998
-
[15]
Silk and L
J. Silk and L. Stodolsky, Phys. Lett. B639, 14 (2006)
2006
-
[16]
Buoninfante, G
L. Buoninfante, G. Gaetano, L. Petruzziello and L. Smaldone, Phys. Rev. D101, 024016 (2020). S. Abel, L. Buoninfante and A. Mazumdar, JHEP 01, 003 (2020). H. Swami, K. Lochan and K.M. Patel, Phys. Rev. D102, 024043 (2020). M. Dvornikov, Eur. Phys. J. C80, 474 (2020). G. Koutso...
2020 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.