A compact analytical approximation for a light sterile neutrino oscillation in matter

The existence of light sterile neutrinos is a long-standing question in particle physics. Several experimental “anomalies” might be explained by introducing eV mass scaled light sterile neutrinos. Many experiments are actively searching for such light sterile neutrinos through neutrino oscillation. For long baseline experiments, the matter effect should be treated carefully for precise calculation of the neutrino oscillation probabilities. However, this is usually time-consuming or analytically complex. In this manuscript, we adopt a Jacobi-like method to diagonalize the Hermitian Hamiltonian matrix and derive analytically simplified neutrino oscillation probabilities for 3 (active) + 1 (sterile)-neutrino mixing for a constant matter density. These approximations can reach a considerably high numerical accuracy while retaining their analytical simplicity and fast computing speed. This would be useful for current and future long baseline neutrino oscillation experiments.


Introduction θ 13
± Neutrino oscillation has been indisputably established by atmospheric, solar, reactor, and accelerator experimental results [1]. After recent reactor experiments [2][3][4] discovered the last unknown mixing angle in the 3neutrino mixing framework, neutrino oscillation measurements entered a precision era. The nonzero neutrino mass provides convincing evidence of new physics beyond the standard model. Introducing right-handed neutrinos is a natural way to introduce the neutrino mass. In the standard electro-weak V-A theory, right-handed neutrinos cannot couple with W and Z 0 bosons. While the electron collider experimental data [5] constrain the number of active light neutrino flavors to three, other new types of light neutrinos must be sterile. Currently, there is no theoretical constraint on the sterile neutrino mass. They could be very massive (10 15 GeV), as suggested by the see-saw mechanism; they could also be dark matter in the keV mass range; in addition, they might also be as light as sub-eV, which would explain the CMB measurement. If light sterile neutrinos mix with the active neutrinos, their signature might be observed with neutrino oscillation experiments. The LSND observed signal events from the source from decay at rest, which suggests a sterile neutrino with mass greater than 0.4 eV [6,7]. Recently, the MiniBooNE experiments reported a excess of electron-like events when combining both the and beam configurations. The significance of the combined LSND and MiniBooNE excesses can even reach 6 [8], although the source of the low energy excess from MiniBooNE remains unclear. Experimental hints of the existence of eV mass scaled sterile neutrinos have also been derived from short baseline reactor neutrino experiments [9][10][11]. However, the uncertainties associated with theoretical reactor antineutrino flux calculations might be underestimated, leading to an observed excess of antineutrino events at 4 -6 MeV relative to the predictions [12][13][14][15][16]. Therefore, whether the reactor anomaly can be completely attributed to the theoretical modeling or sterile neutrinos remains uncertain.
It is worth mentioning that although eV-scale sterile neutrinos can help explain several experimental anomalies, they are not quite theoretically motivated. In addition, they are also in tension with the muon neutrino disappearance results, especially with respect to recent results from IceCube [17] and MINOS/MINOS+ [18]. The most recent combined analysis with the MINOS+, Bugey, and Daya Bay experiments set a very strong limit on sterile neutrino mixing [19], which can almost completely exclude the LSND and MiniBooNE sterile neutrino hypothesis in the eV-scale region. However, a more convincing and direct testing would come from muon decay at rest experiments, such as the proposal of JSNS2 [20]. Certainly, the existence of eV mass scale sterile neutrinos therefore requires further evidence. Many reactor and accelerator neutrino experiments have been actively searching for sterile neutrinos at various mass scales [21][22][23][24][25].
For long baseline accelerator neutrino experiments [26,27], the neutrino matter effect plays an important role in neutrino mass hierarchy [28,29] and CP violation [30] measurements. As first indicated by Wolfenstein, neutrinos propagating in matter will oscillate differently from those in a vacuum [31]. The presence of electrons in matter changes the energy levels of the propagation eigenstates of neutrinos due to charged current coherent forward scattering of the electron neutrinos. Later, Mikheyev and Smirnov [32] further observed that the matter effect can produce resonant maximal flavor transitions when neutrinos propagate through matter at certain electron densities. Super-Kamiokande observes an indication of different solar neutrino fluxes during the night and day for solar neutrinos passing through additional terrestrial matter in the earth at different periods [33]. For sterile neutrino and other new physics explorations, the matter effect has to be calculated carefully and precisely, especially for long baseline neutrino oscillation experiments. Neutrino oscillation in matter can be solved accurately using numerical or analytical calculations [34] with a complex matrix diagonalization algorithm. In practice, analytic approximations are more commonly used in neutrino experiments and are useful for understanding the oscillation features. High-precision analytical expressions for 3-neutrino oscillation in matter have been thoroughly studied [35][36][37][38][39][40][41][42][43][44][45][46][47][48]. Some of them employ the perturbation theory and rely on expansions in parameter . Given the large observed, higher order corrections associated with are needed to achieve numerical accur-acy. Thus, the oscillation expressions usually become quite complicated. In Ref. [49], the Jacobi method was introduced to diagonalize the real Hermitian matrix. It maintains the same analytical expressions for the neutrinos propagating in matter as for those in vacuum in terms of the effective neutrino mixing angles and mass-squared differences in matter. For sterile neutrinos, the oscillation expressions will be very complicated if additional light sterile neutrinos exist [22]. Compared with standard 3-neutrino mixing, the simplest 3 (active) + 1 (sterile)-neutrino mixing has 3 additional mixing angles (i.e., , , and ) and 2 additional CP phases (i.e., and ). Furthermore, as sterile neutrinos do not interact with matter, the neutral current potential for active neutrinos also needs to be considered. N. Klop et al. [50] developed a method to convert 3+1-neutrino mixing with matter effects into a Non-Standard Interaction (NSI) problem in the 3-neutrino mixing case. Here, we follow the rotation strategy introduced in Ref. [49] and adopt the Jacobi-like method [51,52], which can diagonalize the Hermitian complex matrix, to derive analytical approximations for the 3+1-neutrino oscillation in matter. While retaining the simplicity of the formula, these expressions can also achieve very good numerical accuracy and fast calculation speed. This could be very useful for the current and forthcoming neutrino oscillation experiments. This paper is structured as follows. Section 2 introduces the fundamental theory of neutrino mixing and oscillation, including sterile neutrinos and the matter effect. The basic idea behind the Jacobi-like method and the derivation of analytical approximations for sterile neutrino oscillation probabilities are presented in section 3. Lastly, the accuracy of the proposed work is demonstrated in section 4 with two long baseline accelerator neutrino experiments. Further details on the Jacobi-like method and formula derivation are presented in the appendix.
In practice, when sterile neutrinos are much heavier than active neutrinos ( ), due to the finite detector space and energy resolution, the rapid oscillation frequency associated with the large mass-squared differences between the 4th and other mass eigenstates will be averaged out, leading to . The neutrino oscillation equation can then be simpli- with . In this paper, we prefer to use the full oscillation formula to preserve the rapid oscillations induced by sterile neutrinos.

Matter effect
When active neutrinos propagate through matter, the evolution equation is modified by coherent interaction potentials, which are generated through coherent forward elastic weak charged-current (CC) and neutral-current (NC) scattering in a medium. All active neutrinos can interact with the electrons, neutrons, and protons in matter through the exchange of a boson in the NC process. However, only electron neutrinos participate in the CC process with electrons through the exchange of .
For electron neutrinos, the CC potential is proportional to the electron number density.
, where is the Fermi coupling constant, and is the electron number density. The NC potentials generated by electrons and protons will cancel each other because they have opposite signs and the number densities of electrons and protons are basically the same in the Earth. The net NC potential, , is only sensitive to the neutron number density, . Both and need to swap signs for antineutrinos.
For 3-flavor neutrino oscillation, only the CC potential needs to be considered for the electron neutrino eigenstate, while the NC potential is a common term for all neutrino flavors and has no net effect on neutrino oscillation. However, the NC potential cannot be neglected in the 3+1-flavor neutrino case, as sterile neutrinos do not interact with matter. The effective Hamiltonian in the flavor eigenstate representation for 3+1-flavor neutrino mixing is given by where is the neutrino Hamiltonian in vacuum and V is the matter effect potential. and for neutrinos are given by respectively, where is the mass density. Similar to and , both and have to swap signs for antineutrinos. In this work, we assume a constant . If there is no special declaration, will be set to g/cm by default.
The evolution of the neutrino flavor state can be calculated using the Schr dinger equation . After diagonalizing the effective Hamiltonian matrix , we can calculate the neutrino oscillation probability in matter using the equation with the effective mixing matrix and effective masssquared differences . In the following approximations, we will rotate the Hamiltonian from the mass eigenstate. For simplicity, the effective Hamiltonian in the mass eigenstate can be written as where the Hermitian matrix element yields

U U = UR
In this case, the effective mixing yields , where R is the diagonalization matrix on H.

Analytical approximation
As shown in Ref. [34], the exact solution for the effective mixing matrix and effective mass-squared differences can be obtained analytically. However, obtaining higher-precision analytical approximations for neutrino oscillation in matter would be more convenient and time-saving. Here, we introduce a Jacobilike method, which is a unitary transformation operation method, to diagonalize the complex Hermitian matrix. Then, we present the effective mixing matrix and effective mass-squared differences of the 3+1-flavor neutrino mixing framework for both neutrinos and antineutrinos. Consequently, high accuracy can be obtained for the calculation of neutrino oscillation probabilities in matter. The Jacobi-like method, which originates from the × Jacobi eigenvalue algorithm, is an effective matrix rotation approach to diagonalize a complex Hermitian matrix. Here, we start with an example of solving a 2 2 Hermitian matrix. A Hermitian matrix can be diagonalized as with a rotation matrix where , and The choice of a ± sign for A is optional. For simplicity, we select it to be the same sign as that of and in Eq. (8) for the matter effect in the 3+1 framework. The sign in the denominator of is correlated with the exchange of the values of and in Eq. (13). In this work, we adopt ( ) for the i-j submatrix diagonalization if ( ). After rotation, the eigenvalues of M can be obtained as ϕ In summary, this method can be easily used to diagonalize a complex Hermitian matrix through rotation, in which the complex factor is used to address the complex diagonalization.

Application of the Jacobi-like method to 3+1-flavor neutrino mixing
To accurately diagonalize the neutrino Hamiltonian Hermitian matrix using the Jacobi-like method, in principle, we need to perform infinite iterations of a submatrix rotation. However, in practice, with only two continuous rotations on the effective Hamiltonian, we can already obtain the analytical approximations for neutrino oscillation in matter with very high accuracy. The diagonalized Hamiltonian yieldŝ where , and and are the rotation matrices. After some mathematical simplifications, can be expressed as , which has the same form as standard neutrino mixing U. For simplicity, we only show the major results of and in this section. The complete derivations are presented in Chinese Physics C Vol. 44 This is very similar to that in vacuum (i.e., Eq. (5)), except for one additional effective phase in the submatrix . , , , and are the effective angles and phases as functions of E in and . In addition, , , , and are the same as those in vacuum. In the diagonalization process, it is always better to first apply a rotation to the submatrix that has the largest absolute ratio of the off-diagonal element to the difference of the diagonal elements. As is the smallest masssquared difference compared with the others, we can start with the submatrix rotation first.
After the first rotation with the submatrix, we can obtain the effective angle and effective phase represented as functions of and through the combination of : . The and signs in correspond to the neutrino and antineutrino cases, respectively. After the first rotation ( (30) and (64)), we can obtain the eigenvalues of the effective Hamiltonian submatrix as θ 14 After partial diagonalization on the 1-2 submatrix, the off-diagonal elements of the 1-3 and 2-3 submatrices become the relatively largest of the rest of the submatrices for both the neutrino and antineutrino cases due to the smallness of the sterile neutrino mixing angles (i.e., , ). In the second rotation, we adopt the ( ) 1) rotation matrix for the neutrino (antineutrino) case. After the second rotation, we can obtain and as functions of and : In the equation for , the upper (lower) sign in front of corresponds to the NH (IH) (i.e., normal hierarchy (inverted hierarchy)) case, and ( ) corresponds to the neutrino (antineutrino) case. In the above equations, and have different expressions for the neutrinos and antineutrinos. For the neutrino case, , While for the antineutrino case, In this rotation, we can diagonalize the 2-3 (1-3) submatrix for neutrinos (antineutrinos) in Eq. (B11) (Eq. (B45)), resulting in two eigenvalues . The formulas for are for the neutrino case, and for the antineutrino case.
When the mixing between sterile neutrinos and active neutrinos is relatively small and the neutrino beam en- GeV, the off-diagonal elements in the effective Hamiltonian will be very small compared with the diagonal elements after two of the above rotations are performed. Namely, the effective Hamiltonian is approximately diagonalized. So far, all the effective parameters (i.e., , , , and ) in have been presented. The diagonal terms in the effective Hamiltonian in the new representation can be treated as . After subtracting the smallest neutrino (antineutrino) mass ( ), we can obtain the effective neutrino (antineutrino) mass-squared difference as for the neutrino case, and for the antineutrino case.
Till this point, all the effective parameters in 4-flavor neutrino oscillation have been presented, and hence the neutrino oscillation probabilities can be easily calculated using Eq. (9). As both the CC and NC potentials in matter are proportional to the neutrino energy, the values of these effective parameters in and are also energy dependent, as shown in Figs. 1 and 2. The effective matrix in matter introduces two effective mixing angles and , two effective CP phases and , and effective mass-squared differ-

Discussion
, where is an additional parameter introduced from the Jacobi-like method. These effective parameters are clearly energy dependent, as shown in Figs. 1 and 2.
In Fig. 1, when MeV, and are very close to the and values in vacuum. The value of increases (decreases) rapidly up to the maximum (the minimum ) in the neutrino (antineutrino) energy range from MeV to GeV, leading to ( ). In contrast, begins to change after GeV. It can go up to assuming NH for neutrinos and IH for antineutrinos when GeV; it tends to go down to for the other two combinations. When GeV, both the effective CP phases are close to their corresponding vacuum oscillation values ( and ). When the energy increases above GeV, the influence of the matter effect on and is not negligible.
In Fig. 2, the effect of matter also changes the values of the effective neutrino mass-squared differences . When MeV, , , and are close to their vacuum values.
begins to vary when MeV, while starts to change after GeV. In the case of eV with the current sterile neutrino limits, is insensitive to the matter effect when GeV. As the neutrino energy increases, the matter effect shifts the values of the effective more than and when GeV. It should be noted that has a dip structure around GeV for the antineutrino IH case. This feature also shows up in the 3-fla- , and with respect to neutrino energy. In this figure, we assume eV , , , [56], [57] , and . The solid and dashed lines represent the effective angles and phases, respectively. The shift of with a factor of is caused by the transmission of the sign "-" from and , where we set within , in Eq. (18). In the right-hand plot, as the energy rises, tends to be negative. At that point of time, we shift its negative sign to the phase to ensure that is in the interval , resulting in a shift.
In general, as shown in Figs. 1 and 2, the matter effect is negligible on both and 1-2 neutrino mixing when ( ) (or equivalently MeV). When the energy increases, it is clear that the 1-2 neutrino mixing submatrix is affected by matter considerably more than the other submatrices. So does . However, and 1-3 mixing neutrino mixing still hold stable when ( ) (or equivalently GeV). Furthermore, mixing between active and sterile neutrinos has little impact. From a mathematical point of view, in the function of rotation angles yielding A is proportional to the values of , , and , and is inversely proportional to . Hence, the smallness of these mixing angles and large can effectively suppress the values of the corresponding rotation angles to a negligible level in the submatrices. Therefore, after the rotations on the 1-2 and 2-3 (1-3) submatrices of the neutrino Hamiltonian, the effective Hamiltonian matrix is approximately diagonal.
The discussion above pertains to the general feature of our derived oscillation formula. In some particular cases, the oscillation formula can be simplified: In such cases, the neutrino mixing matrix is real and it is not necessary to introduce an extra phase in Eq. (17) for the matrix diagonalization. Thus, this reverts to the original Jacobi method. The neutrino oscillation forms are identical to those in vacuum with , , and .
• No active-sterile neutrino mixing ( The analytical approximations will reduce to 3-flavor neutrino oscillations.

Accuracy of the approximations
All the neutrino oscillation probabilities can be expressed with Eq. (9) based on the effective and calculated in section 3.2. In this section, we first check the accuracy of these approximations. Then, we highlight the accuracy of this work for two specific long-baseline accelerator neutrino experiments, namely the Tokai to Hyper-Kamiokande (T2HK) and Deep Underground Neutrino Experiment (DUNE).

General accuracy
The accuracy of our approximations can be quantified with , which is defined as the numerical difference between the approximations and exact solutions for neutrinos with flavor type converting to type.
To check the validity of our approximations, Figure 3 presents the results, as a function of neutrino energy and travel baseline, on four major neutrino oscillation channels, including , disappearance, and and appearance. The input values for the oscillation parameters are listed in Table 1. Here, we conservatively assume the unknown sterile neutrino associated mixing angles to be as large as and unknown phases to be maximal, i.e, . 20 10 −3 10 −4 As shown in Fig. 3, when the neutrino energy is below GeV, the accuracy is better than and for the neutrino disappearance and appearance channels, re- spectively. The accuracy will improve by one or two orders when the neutrino energy is smaller. The numerical accuracy of our approximations was relatively good for most of the long baseline neutrino oscillation experiments, as the difference is approximately one order smaller than the oscillation probabilities. The inaccuracy of our approach can be mainly attributed to the limited rota-tion iterations that were applied. One approach to further improve the accuracy is to apply some corrections based on the perturbative method after matrix rotation, as shown in Ref. [58]. Naturally, this will increase the complexity of the analytical expression. Our current approximations achieve a good balance between the numerical accuracy and simplicity of the approximation expres-  sions. It should be noted that the exact accuracy of the approximations depends on the exact neutrino mixing parameter inputs. Once the sterile neutrino mixing angles are smaller than and is larger than eV 2 , the accuracy of our expressions will further impove compared with the current evaluation.

T2HK and DUNE
∼ As shown in Eqs. (7) and (8), the matter effect is proportional to the neutrino beam energy and its propagation distance. Hence, the matter effect can significantly modify the neutrino oscillation features for long-baseline neutrino oscillation experiments, such as T2HK and DUNE. These experiments have relatively high energy beams at GeV and baselines of hundreds and thousands of kilometers. Here, we use them as examples to demonstrate 3+1-neutrino oscillation and check the accuracy of our approximations. 1300 θ 23 DUNE is a next generation on-axis long-baseline accelerator neutrino experiment. It proposes the use of Liquid Argon (LAr) detectors located deep underground, km away from the beam source. Its main physics goals are to solve three challenging issues in the neutrino sector, namely neutrino mass hierarchy, CP asymmetry, and the octant of . It can search for electron and tau neutrino (anti-neutrino) appearance and muon neutrino (antineutrino) disappearance channels from both the neutrino and antineutrino beam modes. 295 2.5 • T2HK is a proposed long-baseline experiment that has the primary objective of measuring CP asymmetry. The far detector is km away and off-axis from the J-PARC beam in Japan, when using the water Cher- Considering the existence of sterile neutrinos with a relatively large mass-squared difference eV , the high-frequency oscillation feature is clearly shown in the muon neutrino disappearance and electron neutrino appearance modes in Fig. 4. Given the matter effect and CP-violation phases, the electron neutrino and antineutrino appearance probabilities and are very different. Compared with the numerical calculations, the accuracy of the analytical approximations can reach in the case of Table 1 for neutrinos and antineutrinos, respectively, for the appearance mode. For the disappear- ance mode, the accuracies of and can reach in the case of Table 1. In Fig. 5, we compare the accuracy of this work with that of two previous studies for the T2HK and DUNE experiments. Our work clearly shows about an order of better accuracy compared with the other two [50,60], especially for the DUNE experiment. To achieve a higher accuracy of the approximation, we can continue introducing a perturbation correction later on the effective neutrino mixing and mass-squared differences, as adopted in [58]. However, given that the accuracy of this work is already substantially good for current and near-future experiments, we do not believe this is necessary.  Table 1 except the phases. We adopt and according to the convention used by N. Klop et al. [50]. Similar to the examples shown in their paper, we combine this with a Neutral Current NSI solution [59] to produce the oscillation probability.

Summary
The search for light sterile neutrinos is an area of great interest in the neutrino field. Many long baseline neutrino experiments continue to actively search for light sterile neutrinos in various mass regions. Both CC and NC induced matter effects are quite important for those experiments. The analytical approximations of neutrino oscillation are preferred in experimental neutrino research because they save considerable time and are helpful for understanding the neutrino oscillation features.
In this manuscript, we introduced a Jacobi-like method to derive simplified analytical expressions with good U ∆ m 2 i j (i, j = 1, 2, 3, 4) accuracy for neutrino oscillation in matter. Compact expressions of the effective mixing matrix and effective mass-squared differences were presented. The accuracy of this work is sufficient for a majority of long baseline neutrino experiments.
In addition, the Jacobi-like method is a general method for diagonalizing complex Hermitian matrices. It can also be extended to other physics topics, such as 3 (active) + N (sterile)-neutrino mixing and neutrino nonstandard interactions.
We would like to extend our thanks to Yu-Feng Li for the valuable discussions and to Neill Raper for helping with the editing of this manuscript. In the 3 + 1 framework, neutrino mixing can be written as a matrix (5). 6 rotation angles with 3 Dirac phases 1) are found in this matrix. All the elements of the mixing matrix are listed in 4 × 4 Table A1. Indeed, if we set the angles and Dirac phases introduced by sterile neutrinos to 0, the matrix will reduce to 3-flavor neutrino mixing.
Appendix B: Jacobi-like method 4 × 4 The matter effect for the 3+1 framework is more difficult to calculate than for the 3 neutrino framework because additional parameters are involved in the neutrino Hamiltonian, which is a complex Hermitian matrix and difficult to diagonalize. In this work, we adopt a rotation technique known as the Jacobi-like method to solve the diagonalization for complex Hermitian matrices. In the following subsections, we provide all the technical details for the neutrino and antineutrino cases separately.

B.1 Neutrino case
In this subsection, we show the diagonalization of the effective Hamiltonian with the matter effect and simplify the expressions of the effective mixing and mass-squared differences for the neutrino case. In the case of neutrinos, we find that the absolute values of the elements and in Eq. (10) are the relatively largest off-diagonal values because of the smallness of . Hence, we should di-agonalize the 1-2 submatrix first.

First rotation:
The rotation matrix can be written as We employ the Jacobi-like method to derive yielding with the eigenvalues of the 1-2 submatrix The corresponding off-diagonal terms in are

(B8)
After the first rotation, we can diagonalize the 2-3 submatrix later because it has a relatively large off-diagonal element and is useful for simplifying the effective matrix U using Eq. (B22).
Second rotation: The second rotation matrix yields The diagonal angle is compatible with where the upper sign is for NH , and the lower sign is for IH . is the amplitude of in Eq. (B3) yielding . The additional complex factor is given by . After two rotations, we obtain where the diagonal terms and obey The corresponding off-diagonal elements obey (B16) After two continuous rotations, the matrix is approximately diagonalized. We can evaluate the matrix diagonalization using the ratio defined as which compares the relative size of the off-diagonal elements in with the difference between the two corresponding diagonal terms. Figure B1 shows the ratios for various submatrices with respect to the neutrino energy. All of these are considerably smaller than , which means that the corresponding rotation angles on the submatrix of will be considerably smaller than based on section 3.1. Thus, the matrix can be treated as an approximately diagonalized matrix in our application range.
After diagonalization, the effective neutrino mixing matrix and mass-squared differences are given by , and will be 0. To establish the beauty of the mathematical form (it is convenient for to completely have the same form with U) and a good understanding of the matter effect on the oscillation parameters, we continue to simplify the effective mixing matrix below. First, we can observe that there are two 1-2 submatrix rotations and next to each other in Eq. (B18b). Therefore, we combine them into one submatrix. The corresponding processes are given by  Table 1. Here, we set the same limit on in with .

B.1.2 Simplification
To maintain a similar oscillation expression to the extent possible as its form in vacuum, we replace with where can be expressed as and yields To check the validity of Eq. (B22), we compare all the corresponding elements between and . As shown in Fig. B2, the differences between those two matrices are quite small in our application range.
Subsequently, we obtain  ) cos θ 13 sin θ 13 (1 + tan 2 Equally, we set the limit on within . We find (B28) For the same purpose of simplification, we replace by where D 13 (e iΘ 13 , 1, e −iΘ 13 , 1) except for an additional . We conclude that if , we can obtain , , , and . Meanwhile, the effective Hamiltonian becomes a real Hermitian matrix. In that case, our method reduces to the Jacobi method, which is a way to address the diagonalization of real Hermitian matrices.

B.2 Antineutrino case
In this subsection, we diagonalize the effective Hamiltonian in matter and simplify the expressions of the effective mixing and mass-squared differences for the antineutrino case.

B.2.1 Diagonalization process
In the case of antineutrinos, we rotate the 1-2 submatrix first, similar to the neutrino case.

First rotation: The rotation is given by
with and the additional complex factor . After the first rotation, we obtain the new effective Hamiltonian

(B42)
Second rotation: The rotation matrix is given by (c The rotation angle yields where the upper sign is for NH , and the lower sign is for IH with and a corresponding complex factor . After the two operations, we obtain the new effective Hamiltonian in the new representation:  (color online) Differences between and . The solid and dashed lines correspond to NH and IH, respectively. The input oscillation parameters are listed in Table 1. Figure B4 shows that the off-diagonal terms of are negligible as a good approximation. So far, we have obtained the effective mixing matrix and effective mass-squared . Using the above rotations, we obtain Similar to the neutrino case, we combine and into one submatrix below:   Figure B5 quantifies the differences between all the elements of and . It shows that the differences are allowable for exchange in the application range of our approxima-H ′′  Table 1.
R 13 (ω 2 , ϕ 2 ) R 12 Fig. B5. (color online) Differences between and . The solid and dashed lines correspond to NH and IH, respectively. The input oscillation parameters are listed in Table 1. It can be observed that can be simplified as We also set a bound on within . Subsequently, we obtain U ≈ R 34 R 24 R 14 R 13 R 13 ( θ 13 , δ 13 )D 13 R 12 .

U
In this subsection, we summarize the effective mixing matrix for both the neutrino and antineutrino cases based on the results of B.1 and B.2 in Table B1. Using the neutrino oscillation probability functions in Eq. (9) and the elements in Table B1, every neutrino oscillation probability is available. D 13 (e iΘ 13 , 1, e −iΘ 13 , 1) R 12 R 12 D 13 (e iΘ 13 , 1, e −iΘ 13 , 1)

Fig. B6.
(color online) Differences between and . The solid and dashed lines correspond to NH and IH, respectively. The input oscillation parameters are listed in Table 1.