Pith. sign in

REVIEW 3 major objections 4 minor 25 references

Fibonacci Fast Convergence for Neutrino Oscillations in Matter

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A simple rotation sequence diagonalizes neutrino oscillations in matter with errors that shrink at a Fibonacci rate, beating perturbation theory.

desk verdict Useful, honest paper: the Fibonacci convergence of the rotation method is real and well demonstrated, but the solar-resonance rescue is numerically shown, not proven. read the letter →

arxiv 1909.02009 v2 pith:UE6X4D5R submitted 2019-09-04 hep-ph hep-th

Peter B. Denton , Stephen J. Parke , Xining Zhang This is my paper · ORCID
classification hep-phhep-th PACS 14.60.Pq
keywords neutrinooscillationsmattereffectsHamiltoniandiagonalizationrotationmethodJacobiFibonacciconvergencesolarresonanceperturbationtheory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper claims that repeated 2x2 rotations of the three-flavor neutrino Hamiltonian in matter converge exponentially fast, with the error order following the Fibonacci sequence. This matters because the exact diagonalization of the matter Hamiltonian is analytically unwieldy and ordinary perturbative expansions degrade near level crossings. The authors identify a concrete optimal recipe: first apply the vacuum (2-3) rotation, then at each step rotate away the largest remaining off-diagonal element. With this recipe, a few rotations suppress the off-diagonal Hamiltonian to order epsilon to a Fibonacci power, yielding high-precision oscillation probabilities at constant per-step cost.

What carries the argument

The load-bearing mechanism is the Fibonacci recursion of rotation orders. At each step one selects a 2x2 submatrix of the Hamiltonian, diagonalizes it with a unitary rotation, and repeats; the pivot is chosen as the largest off-diagonal element (the LODE strategy). The recursion works because a rotation with sine of order epsilon^a converts the existing off-diagonal scales epsilon^a and epsilon^b into new scales epsilon^b and $epsilon^{{a+b}}$, so the leading order after consecutive rotations follows Fibonacci numbers. Near the solar resonance, the argument is saved by a cancellation between the numerator and denominator of the correction metric, which the paper emphasizes occurs only for the special sequence of rotations it identifies.

What would settle it

At a matter potential just above the solar resonance, for example Ye rho E = 0.250 GeV g/cm3, implement the LODE rotation sequence listed in the paper and compute max_{j>k} |2E (H1)_jk / $\Delta$ lambda_jk| after each rotation; the paper predicts a suppression to about $10^{-21}$ after seven rotations, so a failure to reach comparably fast suppression would falsify the Fibonacci convergence claim.

Watch

Extended reading notes

Core claim

The central discovery is that, after an initial vacuum (2-3) rotation, the off-diagonal entries of the neutrino Hamiltonian are hierarchically ordered, and a pivot rule that always selects the largest off-diagonal element produces a Fibonacci recursion in the orders of smallness. If the leading off-diagonal scales are epsilon^a and epsilon^b, a rotation with angle of order epsilon^a creates new off-diagonal terms of order epsilon^b and $epsilon^{{a+b}}$; hence after each rotation the leading order is the sum of the two previous orders, which is the Fibonacci rule. For the fully populated case with a=b=1, the off-diagonal Hamiltonian after N rotations scales as $epsilon^{{F_{N+1}}$}, the eigenvector corrections scale as $epsilon^{{F_{N+1}}$}, and the eigenvalue corrections scale as $epsilon^{{2F_{N+1}}$}. The paper also finds that near the solar resonance, where the diagonal gap is small, a specific rotation sequence produces a numerator-denominator cancellation in the correction metric, and numerical checks confirm that the largest-off-diagonal strategy remains superior to the largest-rotation-angle strategy for all matter potentials.

Load-bearing premise

The Fibonacci rate assumes that the two diagonal entries selected for each rotation differ by an order-one amount in units of the small parameter, so the rotation angle inherits the size of the off-diagonal element it removes; at the solar resonance this gap is itself small, and the claimed rate survives only through a numerator-denominator cancellation that is verified numerically, not derived.

Editorial extensions

If this is right

  • After a handful of rotations, the off-diagonal Hamiltonian is suppressed to order epsilon^{F_{N+1}}, so the eigenvector errors reach precision such as 10^-8 with only four matter rotations at typical energies near the solar resonance.
  • Because each additional rotation is a single 2x2 diagonalization with constant complexity, the rotation method becomes increasingly more efficient than perturbation theory as the desired precision grows.
  • The optimal strategy begins with the vacuum (2-3) rotation, which makes the Hamiltonian real and leaves the matter potential unchanged; the largest-off-diagonal pivot rule then performs at least as well as the largest-rotation-angle rule for all energies.
  • The Fibonacci convergence rate is expected to extend beyond three-flavor oscillations, including scenarios with non-standard interactions or additional sterile neutrinos.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The Fibonacci counting argument is not specific to neutrinos: any nearly diagonal Hermitian matrix with hierarchical off-diagonal entries and order-one diagonal gaps should exhibit a similar rotation-sequence convergence, which may be useful in other coupled-oscillator or quantum-mechanical contexts.
  • The special cancellation at the solar resonance suggests that a fully analytic proof of Fibonacci convergence at all matter potentials may exist, but the paper only verifies the cancellation numerically, so establishing such a proof would strengthen the general claim.
  • The sharp contrast between the largest-off-diagonal and largest-rotation-angle pivot rules hints that pivot choice in iterative diagonalization can matter much more than typical matrix-analysis folklore suggests, which could inform numerical algorithms outside neutrino physics.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies a Jacobi-like rotation method for diagonalizing the 3-flavor neutrino Hamiltonian in matter, in which one successively applies 2x2 unitary rotations to reduce off-diagonal elements. It claims that, with a suitable pivot-selection strategy (largest off-diagonal element, LODE) and after an initial vacuum (2-3) rotation, the size of the off-diagonal part decreases with the number of rotations following the Fibonacci sequence, leading to exponential convergence. The authors also compare LODE with the largest-rotation-angle strategy (LROT), concluding that LODE is superior or equal at all matter potentials, and they provide numerical figures and tables over a range of neutrino energies. A generic order-counting argument in Sec. 2.2 is used to explain the Fibonacci rate, and the special behavior near the solar resonance is handled by a claimed numerator-denominator cancellation in the eigenvector-correction metric.

Significance. If the result holds, the paper offers a practical high-precision approximation scheme for neutrino oscillations with a simple, constant-cost iterative step and a striking exponential convergence rate. The order-counting derivation of the Fibonacci rate is elegant and, away from resonance, convincing; the LODE-versus-LROT comparison is a genuine practical insight. The paper does not fit any free parameter to the observed convergence rate: the rotation angles are algebraic functions of the Hamiltonian and the numerical inputs are standard global-fit neutrino parameters. The central weakness is that the all-matter-potential claim, including the solar-resonance region, rests on a cancellation that is numerically demonstrated for one benchmark rather than derived, and the exact-resonance point is not explicitly treated.

major comments (3)
  1. [Sec. 2.2, Eqs. (8)-(10) and Sec. 3.1, Eq. (22)] The generic Fibonacci order-counting argument assumes that the diagonal gap of the pivot pair is O(1) in units of the small parameter epsilon, so that sin(alpha) inherits the order of the off-diagonal element. At the solar resonance this hypothesis fails, as the paper acknowledges. The rescue is the statement, after Eq. (10) and again after Eq. (22), that a special sequence produces a cancellation between numerator and denominator in the eigenvector corrections. This cancellation is not derived: Eq. (22) is an asymptotic expansion at a approximately Delta_m21^2, and no explicit expression for the eigenvector-correction metric is given to show that the cancellation persists through all subsequent rotations. Since the abstract claims the Fibonacci rate without restricting the matter-potential range, this is a load-bearing gap in the proof of the central claim.
  2. [Sec. 3.1, Eqs. (19)-(22) and Sec. 3.3, Eq. (27)] At the exact solar resonance, a = Delta_m21^2 cos(2 theta12)/c13^2, the diagonal elements lambda_- and lambda_0 coincide, so the metric in Eq. (27), max|2E (H1)_jk / Delta_lambda_jk|, has a zero denominator for the (1-2) pair after the vacuum (2-3) and matter (1-3) rotations. The first-order perturbation theory used to define this metric is not valid at this point, and the paper does not specify how the degenerate subspace is handled before the next rotation is applied. The numerical figures may avoid the exact resonance point, but the text claims convergence 'at the solar resonance' and therefore needs an explicit treatment of the degenerate case.
  3. [Sec. 3.3, Table 2 and Fig. 2] The numerical evidence for the resonance-region cancellation is limited to one benchmark (Ye rho E = 0.250 GeV g/cm^3) in Table 2 and to curves in Fig. 2 without numerical values exactly at the resonance or scans over the oscillation parameters that enter the cancellation (theta13, theta12, the mass ordering, and delta_CP). No code is provided. Given that the claimed all-matter-potential Fibonacci convergence depends on the stability of this cancellation, the paper should supply either an analytic proof or a more systematic numerical demonstration that includes the exact resonance point and variations of the relevant parameters.
minor comments (4)
  1. [Sec. 3.3] The line 'sin 2theta12 = 0.31, sin 2theta13 = 0.022, sin 2theta23 = 0.58' is presumably a typo for sin^2(theta12), sin^2(theta13), and sin^2(theta23); the numerical values are the global-fit values for the squared sines, not for sin(2theta).
  2. [Introduction] The expression 'cos(1 3 cos^{-1}(...))' appears to be a formatting error and should read 'cos(1/3 cos^{-1}(...))'.
  3. [Eqs. (18) and (B.1)] The matrices displayed after Eq. (18) and in Eq. (B.1) are hard to parse because the 2x2 block structure is not typeset clearly; please add explicit row and column indices.
  4. [Sec. 3.3, Eq. (27)] The metric in Eq. (27) is called the size of the first-order corrections to the eigenvectors, but the relationship between that metric and the usual first-order perturbation expression should be stated more explicitly, in particular the definition of Delta_lambda_jk and the role of the factor 2E.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the Fibonacci convergence claim is derived from the rotation order-counting and checked numerically with standard inputs; no fitted parameter is renamed as a prediction.

full rationale

The central claim of Sec. 2.2 is an order-counting derivation from the explicit rotation formulas (Eqs. 2-5 and Eq. 9). After a rotation with sin α ~ O(ε^a), Eq. 10 shows the remaining off-diagonal orders become b and a+b, and the next rotation then gives orders a+b and a+b+b, i.e. the Fibonacci recurrence. This is a mathematical consequence of the algebraic rotation formulas, not an input fitted to the claimed convergence rate. The neutrino application in Sec. 3 specifies a definite strategy—vacuum (2-3) rotation followed by the LODE pivot—and verifies the resulting corrections numerically in Fig. 2 and Table 2 using standard global-fit neutrino parameters. No free parameter is tuned to make the Fibonacci rate emerge. The solar-resonance caveat is explicitly acknowledged: Sec. 2.2 states that the generic argument fails when the diagonal gap is O(ε), and Sec. 3.1 supplies the analytic estimate sin(θ~13 - θ13) = s13 c13 (a/Δm2_ee) {1 + O(...)} (Eq. 22) that produces the numerator-denominator cancellation. The paper also notes that "this cancellation occurs only for a special choice of the sequence of rotations." The stability of this cancellation is verified numerically rather than proven for all parameters, but that is a completeness or correctness concern, not circularity. Self-citations (e.g. refs. [11], [13], [20], [21]) describe the rotation method or related eigenvalue results; they do not supply or pre-empt the Fibonacci-rate result, so they are not load-bearing for the paper's central claim. The derivation is self-contained: no equation is defined in terms of the result it is said to predict, and no fitted quantity is repackaged as a prediction.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

All numbers used in the numerical section are standard neutrino parameters from a global fit; nothing is fitted in this paper. The algorithm's rotation angles are deterministic functions of the Hamiltonian, so no free parameters are introduced. The central derivation relies on the standard matter Hamiltonian for constant density, on the measured hierarchy of mixing and mass parameters, and on elementary properties of 2x2 unitary transformations. One additional ad hoc input is the asserted cancellation at the solar resonance, which is not proven in full generality. No new physical entities are postulated.

assumptions (4)
  • domain assumption The Wolfenstein matter Hamiltonian in Eq. 12 correctly describes three-flavor neutrino propagation in constant matter density.
    Invoked as the starting Hamiltonian; the matter potential term diag(a,0,0) is treated as exact and constant. The paper's whole convergence analysis is conducted on this matrix.
  • domain assumption The measured neutrino parameter hierarchy holds: ϵ = Δm21^2/Δm_ee^2 ≈ 0.03, s13 ≈ 0.15, giving off-diagonal orders 0.15 ≫ 0.015 ≫ 0.0021 (Eq. 17).
    The Fibonacci exponent bookkeeping and the ordering of LODE versus LROT performance depend on this hierarchy. The authors note if s13 were an order of magnitude smaller the initial strategy ranking changes, although LODE still wins.
  • standard math The 2x2 complex rotation in Eqs. 2-5 diagonalizes the selected submatrix and leaves the rest of the Hermitian matrix with bounded elements.
    Standard Jacobi-type rotation property used throughout Sec. 2 and Appendix B.
  • ad hoc to paper At the solar resonance the matrix elements and differences of eigenvalues satisfy a special cancellation for the LODE sequence so that the metric in Eq. 27 remains small.
    The generic Fibonacci argument assumes diagonal gaps of order O(1), which fails at the solar resonance. The paper asserts the cancellation (Sec. 2.2 'we emphasize that this cancellation occurs only for a special choice of the sequence of rotations') and verifies it numerically, but does not prove it in full generality.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Fibonacci Fast Convergence for Neutrino Oscillations in Matter." pith.science (2026). https://pith.science/paper/UE6X4D5R

@misc{pith2026190902009,
  author       = {Pith},
  title        = {Pith review of: Fibonacci Fast Convergence for Neutrino Oscillations in Matter},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UE6X4D5R}},
  note         = {Machine review of arXiv:1909.02009}
}
read the original abstract

Understanding neutrino oscillations in matter requires a non-trivial diagonalization of the Hamiltonian. As the exact solution is very complicated, many approximation schemes have been pursued. Here we show that one scheme, systematically applying rotations to change to a better basis, converges exponentially fast wherein the rate of convergence follows the Fibonacci sequence. We find that the convergence rate of this procedure depends very sensitively on the initial choices of the rotations as well as the mechanism of selecting the pivots. We then apply this scheme for neutrino oscillations in matter and discover that the optimal convergence rate is found using the following simple strategy: first apply the vacuum (2-3) rotation and then use the largest off-diagonal element as the pivot for each of the following rotations. The Fibonacci convergence rate presented here may be extendable to systems beyond neutrino oscillations.

Figures

Figures reproduced from arXiv: 1909.02009 by the authors.

Figure 1
Figure 1. Here we show the relative growth in precision by using rotations or following per [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. The size of the corrections to the eigenvectors, max [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 7 canonical work pages

  1. [1]

    Wolfenstein, Neutrino Oscillations in Matter, Phys

    L. Wolfenstein, Neutrino Oscillations in Matter, Phys. Rev. D17 (1978) 2369–2374, [,294(1977)]. doi:10.1103/PhysRevD.17.2369. 13

    doi:10.1103/physrevd.17.2369 1978

  2. [2]

    Neutrino oscillation probabilities through the looking glass

    G. Barenboim, P. B. Denton, S. J. Parke, C. A. Ternes, Neutrino Oscillation Probabilities through the Looking Glass, Phys. Lett. B791 (2019) 351–360. arXiv:1902.00517, doi:10.1016/j.physletb.2019.03.002

    work page Pith review arXiv 2019

  3. [3]

    C. G. J. Jacobi, ¨Uber ein leichtes verfahren, die in der theorie der s¨ akularst¨ orungen vorkommenden gleichungen numerisch aufzul¨ osen, Crelle’s Journaldoi:10.1515/crll.1846.30.51

    doi:10.1515/crll.1846.30.51

  4. [4]

    S. T. Petcov, S. Toshev, Three Neutrino Oscillations in Matter: Analytical Results in the Adiabatic Approximation, Phys. Lett. B187 (1987) 120–126. doi:10.1016/0370-2693(87)90083-9

    work page doi:10.1016/0370-2693(87)90083-9 1987

  5. [5]

    A Simple Parameterization of Matter Effects on Neutrino Oscillations

    M. Honda, Y. Kao, N. Okamura, T. Takeuchi, A Simple parameterization of matter effects on neutrino oscillations arXiv:hep-ph/0602115

    work page Pith review arXiv

  6. [6]

    Matter Effect on Neutrino Oscillations from the violation of Universality in Neutrino Neutral Current Interactions

    M. Honda, N. Okamura, T. Takeuchi, Matter Effect on Neutrino Oscilla- tions from the violation of Universality in Neutrino Neutral Current Inter- actionsarXiv:hep-ph/0603268

    work page Pith review arXiv

  7. [7]

    Kopp, Efficient numerical diagonalization of hermitian 3 x 3 matrices, Int

    J. Kopp, Efficient numerical diagonalization of hermitian 3 x 3 matrices, Int. J. Mod. Phys. C19 (2008) 523–548. arXiv:physics/0610206, doi: 10.1142/S0129183108012303

    arXiv 2008

  8. [8]

    S. K. Agarwalla, Y. Kao, T. Takeuchi, Analytical approximation of the neutrino oscillation matter effects at largeθ13, JHEP 04 (2014) 047. arXiv: 1302.6773, doi:10.1007/JHEP04(2014)047

    arXiv 2014

Show all 25 references
  1. [9]

    Blennow, A

    M. Blennow, A. Yu. Smirnov, Neutrino propagation in matter, Adv. High Energy Phys. 2013 (2013) 972485. arXiv:1306.2903, doi:10.1155/2013/ 972485

    2013 arXiv

  2. [10]

    Minakata, S

    H. Minakata, S. J. Parke, Simple and Compact Expressions for Neutrino Oscillation Probabilities in Matter, JHEP 01 (2016) 180. arXiv:1505. 01826, doi:10.1007/JHEP01(2016)180

    2016 doi

  3. [11]

    P. B. Denton, H. Minakata, S. J. Parke, Compact Perturbative Expressions For Neutrino Oscillations in Matter, JHEP 06 (2016) 051. arXiv:1604. 08167, doi:10.1007/JHEP06(2016)051

    2016 doi

  4. [12]

    Compact perturbative expressions for neutrino oscillations in matter

    P. B. Denton, S. J. Parke, Addendum to “Compact perturbative expressions for neutrino oscillations in matter”, JHEP 06 (2018) 109. arXiv:1801. 06514, doi:10.1007/JHEP06(2018)109

    2018 doi

  5. [13]

    P. B. Denton, S. J. Parke, X. Zhang, Rotations Versus Perturbative Ex- pansions for Calculating Neutrino Oscillation Probabilities in Matter, Phys. Rev. D98 (3) (2018) 033001. arXiv:1806.01277, doi:10.1103/PhysRevD. 98.033001

    2018 arXiv

  6. [14]

    Martinez-Soler, H

    I. Martinez-Soler, H. Minakata, Perturbing Neutrino Oscillations Around the Solar Resonance, PTEP 2019 (7) (2019) 073B07. arXiv:1904.07853, doi:10.1093/ptep/ptz067. 14

    2019 arXiv

  7. [15]

    S. J. Parke, X. Zhang, Compact Perturbative Expressions for Oscillations with Sterile Neutrinos in Matter arXiv:1905.01356

    1905 arXiv

  8. [16]

    B. Yue, W. Li, J. Ling, F. Xu, A new analytical approximation for a light sterile neutrino oscillation in matter arXiv:1906.03781

    1906 arXiv

  9. [17]

    Yokomakura, K

    H. Yokomakura, K. Kimura, A. Takamura, Matter enhancement of T vi- olation in neutrino oscillation, Phys. Lett. B496 (2000) 175–184. arXiv: hep-ph/0009141, doi:10.1016/S0370-2693(00)01288-0

    2000 arXiv

  10. [18]

    Kimura, A

    K. Kimura, A. Takamura, H. Yokomakura, Exact formula of probability and CP violation for neutrino oscillations in matter, Phys. Lett. B537 (2002) 86–94. arXiv:hep-ph/0203099, doi:10.1016/S0370-2693(02) 01907-X

    2002 arXiv

  11. [19]

    Kimura, A

    K. Kimura, A. Takamura, H. Yokomakura, Exact formulas and simple CP dependence of neutrino oscillation probabilities in matter with constant density, Phys. Rev. D66 (2002) 073005. arXiv:hep-ph/0205295, doi: 10.1103/PhysRevD.66.073005

    2002 arXiv

  12. [20]

    P. B. Denton, S. J. Parke, X. Zhang, Eigenvalues: the Rosetta Stone for Neutrino Oscillations in Matter arXiv:1907.02534

    1907 arXiv

  13. [21]

    P. B. Denton, S. J. Parke, T. Tao, X. Zhang, Eigenvectors from Eigenval- uesarXiv:1908.03795

    1908 arXiv

  14. [22]

    K. J. Kelly, S. J. Parke, Matter Density Profile Shape Effects at DUNE, Phys. Rev. D98 (1) (2018) 015025. arXiv:1802.06784, doi:10.1103/ PhysRevD.98.015025

    2018 arXiv

  15. [23]

    S. F. King, S. Molina Sedgwick, S. J. Parke, N. W. Prouse, Effects of Matter Density Profiles on Neutrino Oscillations for T2HK and T2HKK. arXiv:2001.05505

    2001 arXiv

  16. [24]

    Nunokawa, S

    H. Nunokawa, S. J. Parke, R. Zukanovich Funchal, Another possible way to determine the neutrino mass hierarchy, Phys. Rev. D72 (2005) 013009. arXiv:hep-ph/0503283, doi:10.1103/PhysRevD.72.013009

    2005 arXiv

  17. [25]

    Esteban, M

    I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo, M. Maltoni, T. Schwetz, Global analysis of three-flavour neutrino oscillations: synergies and tensions in the determination ofθ23,δCP , and the mass ordering, JHEP 01 (2019) 106. arXiv:1811.05487, doi:10.1007/JHEP01(2019)106. 15

    2019 arXiv

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.