REVIEW 3 major objections 6 minor 24 references
Analytical Expressions for Neutrino Oscillation
T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper compiles analytic expressions for neutrino oscillation in vacuum and matter and adds a new result: an analytic formula for the density dependence of the mixing angle θ23 at high matter densities.
desk verdict The standard compilation is fine for teaching, but the claimed new θ23 matter result does not follow from the stated derivation. 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 matter Hamiltonian $H_{\text{mat}} = U_{23}U_{13}U_{12}\,\mathrm{diag}(m_1^2, m_2^2, m_3^2)/(2E)\,U_{12}^\dagger U_{13}^\dagger U_{23}^\dagger + V_{CC}\,\mathrm{diag}(1,0,0)$ with $\delta_{CP}=0$. The derivation of the new $\tilde{\theta}_{23}$ formula uses an expansion of the mixing matrix in powers of the small parameter $\epsilon = \cos\tilde{\theta}_{13}$ around the high-density point $\tilde{\theta}_{13}\simeq \pi/2$, keeping the first order in $\epsilon$ and assuming the heaviest mass eigenstate decouples ($\tilde{\lambda}_3 \gg \tilde{\lambda}_{1,2}$). Matching the off-diagonal entries of the expanded Hamiltonian yields an explicit ratio for $\tan\tilde{\theta}_{23}$ in terms of the vacuum angles and $\Delta m^2_{31}$, $\Delta m^2_{21}$, together with the consistency relation $\tilde{\theta}_{12} = \alpha - \pi/2 - \tilde{\theta}_{23}$.
What would settle it
Take the benchmark parameters of the matter-effects section and diagonalize the exact three-family Hamiltonian (Eq. 15) for densities from $10^3$ to $10^6$ g/cm$^3$; extract $\tilde{\theta}_{23}$ numerically and compare with the analytic formula. If the analytic curve disagrees where $\cos\tilde{\theta}_{13}$ is not yet small, or fails to converge asymptotically at the highest densities, the claimed expression is falsified.
Extended reading notes
Core claim
The paper claims that, in the high-density regime where the matter potential dominates, the three-family neutrino Hamiltonian can be diagonalized analytically to first order in the small parameter $\epsilon = \cos\tilde{\theta}_{13}$, yielding an explicit ratio for $\tan\tilde{\theta}_{23}$ in terms of the vacuum angles and mass splittings, together with the consistency relation $\tilde{\theta}_{12} = \alpha - \pi/2 - \tilde{\theta}_{23}$. The authors state that the matter dependence of $\theta_{23}$ is a result not found in the main reviews, and they verify in a figure that the asymptotic values agree with numerical diagonalization. The rest of the paper is a compilation of standard vacuum and matter results, with the two-family and three-family approximations organized by density regime.
Load-bearing premise
The derivation of the analytic $\tilde{\theta}_{23}$ formula assumes the matter density is so high that $\cos\tilde{\theta}_{13}$ is very small ($\tilde{\theta}_{13}\simeq\pi/2$) and the heaviest mass eigenstate dominates; if those conditions are not met, the formula does not give the mixing angle.
Editorial extensions
If this is right
- For any density regime, the paper provides closed-form expressions for the matter-modified angles: $\tilde{\theta}_{12}$ at low densities, $\tilde{\theta}_{13}$ claimed valid for all densities, and $\tilde{\theta}_{23}$ plus $\tilde{\theta}_{12}$ at high densities, so oscillation probabilities can be assembled analytically.
- The new $\tilde{\theta}_{23}$ formula closes a gap in the standard reviews, meaning the matter dependence of the atmospheric mixing angle no longer requires a numerical calculation in the high-density limit.
- In the high-density regime the two high-density angles are linked by $\tilde{\theta}_{12} = \alpha - \pi/2 - \tilde{\theta}_{23}$, so fixing one determines the other.
- The paper's figures show that the analytic asymptotic values for $\tilde{\theta}_{23}$ and $\tilde{\theta}_{12}$ converge to the numerical calculation, supporting the claim that the formulas are correct in their stated limit.
Reading between the lines
- If the $\tilde{\theta}_{23}$ formula is numerically reliable beyond the asymptotic region shown, it could replace interpolated diagonalization in supernova neutrino transport codes, where matter potentials are very large and the standard review formulas stop.
- The same small-$\epsilon$ expansion could be applied with $\delta_{CP}\neq 0$, potentially yielding analytic matter-dependent CP-violating phases; the paper explicitly restricts itself to $\delta_{CP}=0$.
- A fuller analytic coverage of the three-family matter Hamiltonian would require an intermediate-density expression for $\tilde{\theta}_{12}$ bridging the low-density formula and the high-density asymptotic value; the paper leaves that bridge implicit.
- Because the paper's parameters are fixed benchmarks, a straightforward extension is to map the validity range of the $\tilde{\theta}_{23}$ formula across the physically allowed values of $\theta_{13}$ and $\Delta m^2_{31}$.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript presents a pedagogical derivation of neutrino oscillation probabilities in vacuum and in constant matter, together with asymptotic expansions for low and high densities. It covers standard two- and three-flavour results, and its claimed new contribution is an analytic expression for the matter dependence of the mixing angle θ23 at high densities (Section III.B). The paper includes numerical comparisons with exact diagonalization and points to a GitHub repository for the supporting programs.
Significance. If correct, the paper would provide a useful self-contained introduction to neutrino oscillation phenomenology, and the high-density θ23 formula would fill a gap that the authors claim is absent from standard reviews. The derivations of the standard results are mostly standard and the numerical support in Figs. 1–3 is helpful. However, the new θ23 result rests on a single algebraic step that does not survive direct evaluation, so the significance of the paper cannot be assessed until that step is corrected.
major comments (3)
- [III.B, after Eq. (22)] The derivation of tan θ̃23 does not follow from the stated relation. Evaluating (Hmat)12/(Hmat)13 with Eq. (15), δCP = 0, and the PMNS entries of Eq. (3) gives [Δ s13 s23 + δ(s12 c12 c23 − s12² s23 s13)] / [Δ s13 c23 − δ(s12 c12 s23 + s12² c23 s13)], where Δ = Δm31²/(4E) and δ = Δm21²/(4E). The printed formula, after factoring 2c13 from numerator and denominator, is [(Δ + δ cos2θ12) s13 s23 + δ sin2θ12 c23] / [(Δ + δ cos2θ12) s13 c23 − δ sin2θ12 s23]. These expressions differ by terms such as δ c12² s13 s23 + δ s12 c12 c23 in the numerator, and with the Section III parameters they give θ̃23 ≈ 48° and ≈ 54°, respectively. The central new result is therefore not supported by the derivation as written.
- [III.B, Eq. (22)] The high-density expansion is not rigorously controlled. The derivation assumes cos θ̃13 ≈ ε with ε ≪ 1 and λ̃3 ≫ λ̃1,2, but the final formula for tan θ̃23 is written in terms of vacuum parameters and no estimate of the neglected O(ε²) terms or of the required density range is provided. Since this expansion is the basis for the paper's only new result, the domain of validity should be quantified explicitly.
- [III.B, Eq. (22)] The matrix equation defining Ũ is not readable as printed: the 2×2 block O is inserted as a single entry, and the displayed matrix does not have consistent dimensions. In addition, the relation between α and θ12, θ23, and the formula θ̃12 = α − π/2 − θ̃23 are stated without derivation, so the reader cannot verify the angle identifications. Please rewrite Eq. (22) as an explicit 3×3 matrix and derive the relation involving α.
minor comments (6)
- [Eq. (19)] The symbol Δm² should carry the subscript 21 (Δm²₂₁); as printed it is ambiguous with the two-family Δm² of Eq. (11).
- [Section III] Throughout Section III and in Eq. (17), the notation 'm2' is used for the squared mass (e.g., 'm2_1/(2E)'); please use m1² consistently to avoid confusion with the mass m2.
- [General] There are several typos: 'familes' in Section III.A, 'eingenstates' before Eq. (7), and 'more date is needed' near the end of the Introduction.
- [Fig. 4] The figure appears to show only the asymptotic horizontal lines; since the text claims agreement with the numerical calculation, the numerical curves should be shown or the caption should explain their absence.
- [Conclusions] The paper mentions supporting numerical programs in a GitHub repository but gives no URL or reference; please add it.
- [Conclusions] The claim that the θ23 matter dependence is 'not found in the main reviews' should be accompanied by explicit citations to those reviews.
Circularity Check
No significant circularity; the derivations are self-contained from the stated Hamiltonians, with external measured parameters used only as inputs.
full rationale
The paper's central claimed new item, the matter-dependent θ23, is derived by expanding the diagonalizing matrix in powers of ε = cos θ̃13 (Eq. 22) and matching the first-order off-diagonal entries of Ũ diag(λ̃)Ũ† to the known matter Hamiltonian of Eq. (15). The step 'tan θ̃23 = (Hmat)12/(Hmat)13' is the matching equation itself, not a restatement of an input: the right-hand side is computed from the original PMNS matrix and mass splittings, and the left-hand side is the unknown matter angle being solved for. The final closed form is presented as the result of 'some algebraic manipulation' of those entries. No parameter is fitted to the target quantity, no uniqueness theorem from prior work by the authors is invoked, and the only self-reference is a GitHub repository mentioned as supporting numerical programs, which is not load-bearing. The low-density θ̃12, the two-family results, and the eigenvalues all follow from direct diagonalization of the stated Hamiltonians, while the measured oscillation parameters are boundary inputs for the illustrative plots rather than fitted outputs. A possible algebraic discrepancy between the printed tan θ̃23 expression and a direct evaluation of H12/H13 is a correctness concern, not evidence of circularity, because the expression does not reduce to its own input by construction. Therefore no circular step is identified.
Assumptions & free parameters
assumptions (5)
- standard math Standard quantum mechanical evolution of neutrino flavor states via the Schrödinger equation with the PMNS mixing matrix.
- domain assumption Wolfenstein matter potential VCC is the only matter effect, with constant density in the analyzed regimes.
- domain assumption The CP-violating phase δCP is set to zero in the three-family analysis.
- domain assumption Neutrino mass ordering is normal, with Δm2_21 and Δm2_31 positive, and the measured values are taken from experiments.
- ad hoc to paper At high densities, cos θ̃13 is small (ϵ) and the heavy mass eigenstate decouples (λ̃3 ≫ λ̃1,2).
Cite this review
Pith. "Pith review of Analytical Expressions for Neutrino Oscillation." pith.science (2026). https://pith.science/paper/EYZHR7XA
@misc{pith2026241216424,
author = {Pith},
title = {Pith review of: Analytical Expressions for Neutrino Oscillation},
year = {2026},
howpublished = {\url{https://pith.science/paper/EYZHR7XA}},
note = {Machine review of arXiv:2412.16424}
}
read the original abstract
Research in neutrino physics has been very active, both in experimental advances, with a new generation of detectors in operation and planning, and in theoretical discussions regarding the fundamental nature of the neutrino. This scientific dynamism has attracted many new students to the field. One of the first topics studied in neutrino physics by newcomers is the formalism of neutrino flavor oscillations and its associated phenomenology. We present this work as a compilation of this basic knowledge, through a step-by-step approach that facilitates an efficient understanding of this vast theoretical and experimental landscape.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
Y. Fukuda et al. (Super-Kamiokande), Evidence for oscillation of atmospheric neutrinos, Phys. Rev. Lett. 81, 1562 (1998), arXiv:hep-ex/9807003
arXiv 1998
-
[2]
Q. R. Ahmad et al. (SNO), Measurement of the rate of νe + d → p + p + e− interactions produced by 8B solar neutrinos at the Sudbury Neutrino Observatory, Phys. Rev. Lett. 87, 071301 (2001), arXiv:nucl-ex/0106015
arXiv 2001
-
[3]
Y. Cai, J. Herrero-Garc ´ ıa, M. A. Schmidt, A. Vicente, and R. R. Volkas, From the trees to the forest: a review of radiative neutrino mass models, Front. in Phys. 5, 63 (2017), arXiv:1706.08524 [hep-ph]
arXiv 2017
-
[4]
(17) The remaining problem is solved by the diagonalization process of the non-diagonal subspace, which results in λ1,2 = m2 1 + m2 2 4E + VCC c2 13 2 ± s δ21c2θ12 − VCC c2 13 2 2 + (δ21s2θ12 )2. (18) The reduced 2-dimensional subsystem has an analytical solution for the mixing angle, as already shown in Eq. (13), which gives tan(2˜θ12) = ∆m2 4Eν s2θ12 ∆m...
-
[5]
M. J. Dolinski, A. W. P. Poon, and W. Rodejohann, Neutrinoless Double-Beta Decay: Status and Prospects, Ann. Rev. Nucl. Part. Sci. 69, 219 (2019), arXiv:1902.04097 [nucl-ex]
arXiv 2019
-
[6]
C. Giunti and C. W. Kim, Fundamentals of Neutrino Physics and Astrophysics (2007)
work page 2007
-
[7]
Navas et al
S. Navas et al. (Particle Data Group), Review of particle physics, Phys. Rev. D 110, 030001 (2024)
2024
-
[8]
S. Dell’Oro, S. Marcocci, M. Viel, and F. Vissani, Neutrinoless double beta decay: 2015 review, Adv. High Energy Phys. 2016, 2162659 (2016), arXiv:1601.07512 [hep-ph]
arXiv 2016
Show all 24 references
-
[9]
Suekane and T
F. Suekane and T. J. de Castro Bezerra, Double chooz and a history of reactor θ13 experiments, Nuclear Physics B 908, 74 (2016)
2016
-
[10]
Nath and N
A. Nath and N. K. Francis, Detection techniques and investigation of different neutrino experiments, International Journal of Modern Physics A 36, 2130008 (2021)
2021
-
[11]
Adamson et al
P. Adamson et al. (MINOS+), Precision Constraints for Three-Flavor Neutrino Oscillations from the Full MINOS+ and MINOS Dataset, Phys. Rev. Lett. 125, 131802 (2020), arXiv:2006.15208 [hep-ex]
2020
-
[12]
de Gouvˆ ea, G
A. de Gouvˆ ea, G. Jusino S´ anchez, and K. J. Kelly, Very light sterile neutrinos at NOvA and T2K, Phys. Rev. D 106, 055025 (2022), arXiv:2204.09130 [hep-ph]
2022 arXiv
-
[13]
(9) Using the unitarity property of the PMNS matrix, the rows (or columns) of the matrix must satisfy the orthogonality conditions X k UµkU ∗ ek = 0, which can be replaced in Eq
sin2 ∆M 2L 4E . (9) Using the unitarity property of the PMNS matrix, the rows (or columns) of the matrix must satisfy the orthogonality conditions X k UµkU ∗ ek = 0, which can be replaced in Eq. (9). This makes it possible to rewrite the appearance probabilities as Pµe = 4 |Uµ...
-
[14]
Abbasi et al
R. Abbasi et al. ((IceCube Collaboration)*, IceCube), Measurement of atmospheric neutrino mixing with improved IceCube DeepCore calibration and data processing, Phys. Rev. D 108, 012014 (2023), arXiv:2304.12236 [hep-ex]
2023
-
[15]
Abe et al
K. Abe et al. (Super-Kamiokande), Atmospheric neutrino oscillation analysis with external constraints in Super- Kamiokande I-IV, Phys. Rev. D 97, 072001 (2018), arXiv:1710.09126 [hep-ex]
2018 arXiv
-
[16]
Abe et al
K. Abe et al. (T2K), Updated T2K measurements of muon neutrino and antineutrino disappearance using 3.6 ×1021 protons on target, Phys. Rev. D 108, 072011 (2023), arXiv:2305.09916 [hep-ex]
2023
-
[17]
Kajita, Atmospheric neutrinos and discovery of neutrino oscillations, Proceedings of the Japan Academy, Series B 86, 303 (2010)
T. Kajita, Atmospheric neutrinos and discovery of neutrino oscillations, Proceedings of the Japan Academy, Series B 86, 303 (2010)
2010
-
[18]
Gando et al
A. Gando et al. (KamLAND), Constraints on θ13 from A Three-Flavor Oscillation Analysis of Reactor Antineutrinos at KamLAND, Phys. Rev. D 83, 052002 (2011), arXiv:1009.4771 [hep-ex]
2011 arXiv
-
[19]
Wolfenstein, Neutrino Oscillations in Matter, Phys
L. Wolfenstein, Neutrino Oscillations in Matter, Phys. Rev. D 17, 2369 (1978)
1978
-
[20]
Abe et al
K. Abe et al. (T2K, Super-Kamiokande), First joint oscillation analysis of Super-Kamiokande atmospheric and T2K accelerator neutrino data, (2024), arXiv:2405.12488 [hep-ex]
2024
-
[21]
The eigenvalues of the non-diagonal matrix can be calculated through the usual diagonalization process
In addition, the notation c2θ = cos 2θ and s2θ = sin 2θ was employed. The eigenvalues of the non-diagonal matrix can be calculated through the usual diagonalization process. Furthermore, the terms proportional to identity in Eq. (11) should be added to the full eigenvalues: λ1...
-
[22]
Li et al
Z. Li et al. (Super-Kamiokande), Measurement of the tau neutrino cross section in atmospheric neutrino oscillations with Super-Kamiokande, Phys. Rev. D 98, 052006 (2018), arXiv:1711.09436 [hep-ex]. 14
2018 arXiv
-
[23]
sin2 ∆M 2L 4E = 1 − (s2 23 sin2(2θ13) + sin2(2θ23)c4
-
[24]
S. P. Mikheyev and A. Y. Smirnov, Resonance Amplification of Oscillations in Matter and Spectroscopy of Solar Neutrinos, Sov. J. Nucl. Phys. 42, 913 (1985)
1985
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.