Three Neutrino Oscillations in Uniform Matter
Ara Ioannisian, Stefan Pokorski

TL;DR
This paper analytically solves the three-neutrino oscillation problem in uniform matter, providing explicit formulas for oscillation probabilities with modified mixing angles and mass eigenvalues, enhancing understanding of neutrino behavior in matter.
Contribution
It introduces an analytical solution for three-neutrino oscillations in constant matter density using successive diagonalizations, offering a clear parametric form similar to vacuum oscillations.
Findings
Explicit formulas for oscillation probabilities in matter.
Modified mixing angles and mass eigenvalues in matter.
Simplified analytical approach for neutrino oscillations.
Abstract
Following similar approaches in the past, the Schrodinger equation for three neutrino propagation in matter of constant density is solved analytically by two successive diagonalizations of 2x2 matrices. The final result for the oscillation probabilities is obtained directly in the conventional parametric form as in the vacuum but with explicit simple modification of two mixing angles ( and ) and mass eigenvalues.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsNeutrino Physics Research · Quantum, superfluid, helium dynamics · Dark Matter and Cosmic Phenomena
Three Neutrino Oscillations in Uniform Matter
Yerevan Physics Institute, Alikhanian Br. 2, 375036 Yerevan, Armenia
Institute for Theoretical Physics and Modeling, 375036 Yerevan, Armenia
Stefan Pokorski
Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, ul. Pasteura 5, PL-02-093 Warsaw, Poland
Abstract:
Following similar approaches in the past, the Schrodinger equation for three neutrino propagation in matter of constant density is solved analytically by two successive diagonalizations of 2x2 matrices. The final result for the oscillation probabilities is obtained directly in the conventional parametric form as in the vacuum but with explicit simple modification of two mixing angles ( and ) and mass eigenvalues.
The MSW effect for the neutrino propagation in matter attracts a lot of experimental and theoretical attention.
On the theoretical side, a large number of numerical simulations of the MSW effect in matter with a constant or varying density has been performed. Although, in principle, sufficient for comparing the theory predictions with experimental data, they do not provide a transparent physical interpretation of the experimental results. Therefore, several authors have also published analytical or semi-analytical solutions to the Schroedinger equation for three neutrino propagation in matter of constant density, in various perturbative expansions [1, 2, 3].
We solve the Schroedinger equation in matter with constant density by two successive diagonalizations of 2x2 matrices (similar approaches have been used in the past, in particular in ref. [2] and [3]). The final result for the oscillation probabilities is obtained directly in the conventional parametric form as in the vacuum but with modified two mixing angles and mass eigenvalues111The results of this paper have been presented as private communication by one of us (A.I) to the members of the T2HKK collaboration in December 2017., similarly to the well known results for the two-neutrino propagation in matter. The three neutrino oscillation probabilities in matter have been presented in the same form as here in the recent ref. [4], where the earlier results obtained in ref. [3] are rewritten in this form. The form of our final results can also be obtained after some simplifications from ref.[2]. Our approach can be easily generalized to non-constant matter density by dividing the path of the neutrino trajectory in the matter to layers and assuming constant density in each layer.
In the electroweak basis the neutrino Hamiltonian is
[TABLE]
The matrix () is the neutrino mixing matrix in the vacuum (matter). The mass squared differences are defined as () and (, positive sign is for normal mass ordering and negative sign for inverted one). and are eigenvalues of the neutrino Hamiltonian (we always can add proportional to unity diagonal matrix to the neutrino Hamiltonian). is the neutrino weak interaction potential energy ( is electron number density) and we take it in this section to be x-independent.
We work in the auxiliary basis [5, 6, 7] and do two rotations for diagonalization of the neutrino Hamiltonian in eq 1.
For the mixing angles and in matter we get
[TABLE]
[TABLE]
where
[TABLE]
And for differences between eigenvalues of the neutrino Hamiltonian, , we have
[TABLE]
In our approximation 23 angle and the CP phase remain unchanged: and [8].
The oscillation probabilities () have the same forms as for the vacuum oscillations with mass eigenstates as above and with replacements and . For the transition we have
[TABLE]
This approximate solution is valid for all energies. For anti-neutrino oscillations , V -V and . For normal mass hierarchy is positive and for inverted mass hierarchy it is negative.
Our solutions are illustrated in Fig. 1 for oscillation at DUNE distance for several values of and compared with the oscillation probabilities in the vacuum, shown by the dotted curves.
The most important effect is the dependence of the oscillation probability on the angle which has larger (smaller) values in matter than in the vacuum for normal (inverted) neutrino mass hierarchies (and opposite for antineutrinos). Thus the oscillation probabilities have larger(lower) oscillation amplitudes for normal (inverted) neutrino mass hierarchies (and opposite for antineutrinos). In oder words the matter of the Earth is amplifying the effect of the mass ordering on neutrino oscillations. The dependence on the angle enters multiplicatively in the first three terms of eq (7), whereas the fourth term is small in the region of the first maximum. Therefore the matter effects relative to the oscillations in the vacuum do not depend on the value of , as it is seen in Fig. 1 . Moving to the next resonances (lower energies) the difference between oscillations in matter and in the vaccuum remain qualitatively similar, although some small differences can be seen due to the fact that the change in the angle is smaller.
In Fig. 2 we show the accuracy of the analytical solutions comparing them with numerical/exact results.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S. T. Petcov, Phys. Lett. B 214 (1988) 259. E. K. Akhmedov, R. Johansson, M. Lindner, T. Ohlsson and T. Schwetz, JHEP 0404 (2004) 078. A. Cervera, A. Donini, M. B. Gavela, J. J. Gomez Cadenas, P. Hernandez, O. Mena and S. Rigolin, Nucl. Phys. B 579 (2000) 17 Erratum: [Nucl. Phys. B 593 (2001) 731]. H. Nunokawa, S. J. Parke and J. W. F. Valle, Prog. Part. Nucl. Phys. 60 (2008) 338
- 2[2] M. Blennow and A. Y. Smirnov, Adv. High Energy Phys. 2013 (2013) 972485
- 3[3] P. B. Denton, H. Minakata and S. J. Parke, JHEP 1606 (2016) 051
- 4[4] P. B. Denton, H. Minakata and S. J. Parke, ar Xiv:1801.06514 [hep-ph].
- 5[5] P. I. Krastev and S. T. Petcov, Phys. Lett. B 205 (1988) 84.
- 6[6] O. L. G. Peres and A. Y. Smirnov, Nucl. Phys. B 680 (2004) 479
- 7[7] A. Ioannisian and S. Pokorski, Phys. Lett. B 782 (2018) 641
- 8[8] S. T. Petcov and Y. L. Zhou, Phys. Lett. B 785 (2018) 95
