Abstract
We find exact analytical expressions for mixing angles in matter in the context of three generation neutrino oscillations in matter to discuss the role of resonances in this phenomenon. We show that some knowledge from conventional two neutrino MSW effect, which has been extended to approximated solutions to three neutrino oscillations, has to be abandoned in this exact approach. We observe that maximal values for the mixing angles in matter are found in nonresonant regions and stationary phases do not coincide anymore with resonances in this simple extension of the MSW effect. We present a general way to identify a resonance and discuss what we can physically expect in these regions
On Resonances and Mixing Angles in Three Neutrino Oscillations in Matter
J. Bellandi and M.M. Guzzo
Instituto de Física 'Gleb Wataghin'
Universidade Estadual de Campinas, Unicamp
13083-970 Campinas, São Paulo, Brasil
V.M. Aquino
Departamento de Física, Universidade Estadual de Londrina, UEL
86000-970 Londrina, Paraná, Brasil
Received 20 May, 1996
We find exact analytical expressions for mixing angles in matter in the context of three generation neutrino oscillations in matter to discuss the role of resonances in this phenomenon. We show that some knowledge from conventional two neutrino MSW effect, which has been extended to approximated solutions to three neutrino oscillations, has to be abandoned in this exact approach. We observe that maximal values for the mixing angles in matter are found in nonresonant regions and stationary phases do not coincide anymore with resonances in this simple extension of the MSW effect. We present a general way to identify a resonance and discuss what we can physically expect in these regions.
I. Introduction
Resonant regions are previleged zones for neutrino conversion. Concerning solar neutrinos, the importance of a resonance can be appreciated remenbering that the standard MSW solution to the solar neutrino problem requires values for the mixing angle in vacuum q and for the squared mass difference such that sin22q < 10- 3 and [1] which imply a resonance in the neutrino trajectory inside the sun when the approximately exponentially decreasing standard solar matter distribution is assumed [2]. This is the so-called nonadiabatic solution to the solar neutrino problem and the role of the resonance is evident in such situation once that it is well known that the adiabaticity parameter [3] presents its smallest values in a resonance region, which imply that neutrino transitions are less adiabatic in that region.
Resonances in two family MSW effect [3,5] are associated with maximum mixing between the two flavor eigenstates. This can be appreciated investigating the behavior of the matter mixing angle when the relevant matter density varies along the neutrino trajectory. The mixing angle in matter is introduced as the parameter that characterizes a rotation of the two-dimensional neutrino space from the basis of the current eigenstates (ne, nm) to the basis of the physical eigenstates (n1, n2) :
It can be calculated [4]:
where E is the neutrino energy and GFNe(t) is the consequence of electron neutrino coherent forward scattering from electrons in matter, the number density of which at the region reached by neutrinos at instant t is Ne(t) .
From Eq. (2) it is possible to see that is substantially modified by the neutrino coherent scattering from the medium. If and we recover vacuum expressions. When Ne(t) is extremely large, and n1 ® nm while n2 ® ne. An interesting intermediate case occurs when
and the brackets in the denominator of Eq. (2) vanishes. In this point the mixing of flavor eigenstates is maximal, i.e., from Eq. (1) we see that the probability of finding an electron or a muon neutrino in any of the mass eigenstates is 1/2. This feature has been used to characterize a resonance: the maximum of the bell-shaped graph indicates a resonance.
The resonance condition given by Eq. (3) coincides also with the position where the difference of the two squared matter eigenvalues of the corresponding time evolution matrix in matter is a minimum, suggesting that the resonance is the region where transitions between matter eigenstates are most likely to happen.
Finally, it was noticed in reference [6] that the resonance condition (3) coincides also with the condition of existence of a stationary phase [7] in the two neutrino time evolution equations. Such fact allows to investigate the evolution of this neutrino system around a resonance calculating, through the stationary phase method [7], the related Green function. Employing this method it was possible to evaluate [8] the level crossing probability, i.e., the probability of nonadiabatic transitions between matter eigenstates n1 and n2 as an alternative approach to Landau-Zener [9] or Petcov [10] methods.
In this paper we investigate how is the behavior of mixing angles in matter and how to identify a resonance in the context of a three neutrino system oscillating in matter. We assume standard electroweak interactions of neutrinos with matter as well as nonvanishing vacuum mixing angles and nondegenerated mass eigenstates (in vacuum). Therefore we are analysing the simplest extension of the conventional MSW effect [3,5] to the case where three families are present. We verify that the above mentioned three criteria usually used to define a resonance in two neutrino matter oscillations, namely, maximal mixing angles in matter, minimal eigenvalue difference and the presence of a stationary phase, do not lead anymore to the same region in the neutrino trajectory. Note also that these same criteria have been used in approximated solutions to three neutrino oscillations in matter [11,12]. Consequently some of them have to be abandoned. We present, therefore, based on exact analytical expressions for mixing angles in matter, how we can use our previous knowledge coming from two neutrino matter oscillations to arrive to a solid condition defining resonances in three neutrino oscillations and, therefore, an accurate analytical description of the physical consequences around such regions.
II. Analytical solution
A general time evolution equation describing a three level system can be written as an equation for a three-component spinor :
where the hamiltonian h(t) is a 3 x 3 matrix which elements are specified according to the dynamical situation from which a boundary condition F(to) is given. A general solution of Eq. (4) can be written in the the form
where the symbol Exp represents a sum of multiple time ordered integrals [13].
For a time-independent hamiltonian, the solution of Eq. (4) can be obtained by means of the Laplace transformation. Introducing the Laplace transformed Y(p) = L[F(t)], then
and
The solution F(t) depends on the elements of the h matrix and on the roots of the characteristic polynomial of the h matrix
In the particular case we are interested in, where a three neutrino system oscillates in matter, interacting with it through standard electroweak interactions, the h matrix is given by
where M 2 is a diagonal matrix given by
mi2 are the three neutrino squared mass eigenvalues in vacuum,
is the 3 x 3 mixing matrix where Di are the Gell-Mann matrices, , f and w are the mixing angles in vacuum and G is a matrix containing complex phases that we will ignore since we assume CP conservation ().
Since we consider here only standard neutrino interactions with ordinary matter, A matrix has its first element A11 given by
and all others are zero. Note that neutral current contributions to A are proportional to the unit matrix, giving only irrelevant overall phases to the final solution of Eq. (4). GF , E and Ne were previously introduced.
For neutrino propagating in vacuum, A = 0 , and the solution of Eq. (4) is trivial and simply given by
where m 2 is a diagonal matrix with elements
The solution of Eq. (4) in matter, with A being a time-dependent matrix, is given by Eq. (5) and it depends on the specific Ne function describing the electron density. However, when A can be considered a constant matrix, as it is supposed in the adiabatic approximation, Eq. (4) has an exact analytical solution, obtained by Laplace transformation. Furthermore, the A matrix is invariant under a rotation, then, introducing now
we observe that Y(t) satisfies the following differential equation
with boundary condition Y(to) = F(to) and
which can be explicitly written as
where .
On the Laplace space we have
To calculate Y(t) we have to obtain the roots of the characteristic polynomial of the matrix H , det[p1 + iH] = 0, which are given by [14]
where
and
Note now that in vacuum we have
and the Laplace anti-transformation of Eq. (19) reproduces the corresponding solution given by Eq. (13). The roots li of the characteristic polynomial are the squared mass eigenvalues in matter. Because of the arbitrariness in the choice of the order of the roots, we use the above vacuum limit to order the roots in terms of the squared mass eigenvalues in the matter. We define:
Finally, we can write the solution of Eq. (4) in terms of a T transition matrix such that
where the elements of the T matrix, given in terms of the li roots and of the elements of the H matrix, can be written as:
i) diagonal elements:
ii) non diagonal elements (Tij = Tji) :
where
with .
Note also that all well known results for a two neutrino system oscillating in matter can be straightforwardly obtained from the solution given by Eq. (28).
III. Mixing angles in matter
It is well known that the knowledge of the mixing angles in matter is important to study resonant transitions between flavor neutrino states [4]. In order to explicitly write an exact expression for these angles, we define as the mixing angles in the matter. We can write therefore the final solution of Eq. (4) in terms of mixing angles in matter in analogy with what we did in the vacuum case, Eq.(13), using now the final solution given by Eq. (28). This solution can be written in the following way
where
and
In order to get the matter mixing angles we simply compare Eq. (32) with Eq. (28), and after some algebra, we obtain
and
IV. Analysis of the results
Similar results for the mixing angles in matter were found in Ref. [15]. Note however that Eqs. (35)-(37) are general and exact. Furthermore their physical consequences have never been completely analysed. Based on these equations, we show in the following that some knowledge from conventional two neutrino MSW effect, which has been extended to approximated solutions to three neutrino oscillations, has to be abandoned in this exact approach. We observe that maximal values for the mixing angles in matter are found in nonresonant regions and stationary phases do not coincide anymore with resonances in this simple extension of the MSW effect. We present a general way to identify a resonance and discuss what we can physically expect in these regions.
In Fig. 1 it is presented a comparison of the behavior of the quadratic matter eigenvalues and the relevant matter mixing angles and as a function of the parameter A for specially chosen values of vacuum parameters (see the corresponding caption for details. A is given in units of m12). There are two resonances clearly indicated by the minimum difference between the shown quadratic masses. We observe that presents a maximum value in the lower resonance (where is minimum) while shows a maximum in the region of the higher resonance ( is a minimum). Interesting enough, differently from what is expected in the two flavor neutrino oscillations in matter, the conventional MSW effect, and also from what was found in previous approximated analyses of the three neutrino oscillations [11,12], a second peak for the mixing angle is found after the higher resonance [16]. In Fig. 2 we show the same graphs presented in Fig. 1 to evidenced this unexpected behavior of the mixing angle for larger values of A . It is clear from this figure that the criterion of defining a resonance by means of localizing the maximal mixing angle in matter, which can be safely used in two neutrino conventional MSW effect, leads to some ambiguity in the context of three neutrino oscillations and therefore has to be abandoned.
Instead, we can improve this criterion analysing the content of Figs. 3 and 4. Note that the admixture of flavor eigenstates in each of the matter eigenstates can be obtained through , where i = 1,2,3 , a = e, m, t and Uia is given by Eq. (11). Let us write now, as an example, the linear combination of flavor eigenstates in the first matter eigenstates:
Figure 1. Squared matter eingenvalues and the relevant matter mixing angles and as a function of the parameter A are presented. The values of the parameters f, w, and mi2 in vacuum: m32 = 5m22 = 25m12; sin2f = 5 x 10- 4; sin2wcos2f = 5 x 10- 2 were chosen in order to well demonstrate the behavior of these parameters as a function of A . The eingenvalues and the quantity A are given in units of m12.
In Figs. 3 and 4 we show therefore the coefficients of this admixture (values of the vacuum parameters are shown in the corresponding captions). From Fig. 3 we observe that in the lower resonance the mixing of electronic and muonic flavor eigenstates is maximal (when ), while, from Fig. 4, we see that the higher resonance coincides with the maximum admixture of nm and nt, when .
Figure 2. Mixing angles and as a function of the A for larger values of A . The parameters f, w, and mi2 are the same as that of the Fig. 1.
Therefore, although we detected maxima of the mixing angles in matter in regions far from resonances, it is still possible to identify a resonance region searching for maximal mixing between flavor eigenstates. Note also that such maximum are not anymore related with values of for flavor coefficients |Uia| in the way it happened in the conventional MSW phenomenon. This is because there could be nonnegligible contributions from the flavor eigenstate that does not participate in the resonant process. From the unitarity of the mixing matrix, we know that , for i = 1,2,3 . Therefore, in the case where i = 1 and is not vanishing, the maximal mixing is such that . A similar situation occurs for the higher resonance where we obtain . We can say that in three neutrino oscillation phenomenon the mixing between flavor eigenstate around a resonant region is as maximal as possible, although not in the same way as in two neutrino oscillations, where maximum mixing implies that each one of the neutrino flavor eigenstate participating in the resonant process contributes with 50% to the matter eigenstates.
Figure 3. The squared mass difference , and the quantities and quantities are presented as a function of the energy of neutrinos. The values of the parameters f, w and mi2 in vacuum are: m32 = 1.445 x 10- 4eV 2; m22 = 10- 8eV 2; m12 = 0; sin2f = 5 x 10- 4; sin2w = 0.050025. The squared masses are given in units of m22.
A final issue to be discussed is the criterion of identifying a resonance looking for a stationary phase in the neutrino evolution equations (4), in the same way it was proposed in reference [6] in the context of two neutrino MSW effect. A stationary phase is given by the smallest difference of any two diagonal elements of the relevant evolution matrix when one of these elements is time dependent. As an example, we quote solar neutrinos where the matter density considerably varies along the neutrino trajectory from the center of the sun, where neutrinos are created, to the solar surface. Although in two neutrino oscillations this criterion can be safely used, it does not work anymore in the presently analysed three neutrino MSW effect. Stationary phases do not coincide with the minimum squared mass differences or maximum flavor admixture.
Figure 4. The squared mass difference , the effective and quantities are presented as a function of the energy of neutrinos. The values of the parameters f, w and mi2 in vacuum are the same ones used to draw Fig. 3.
Note however that it is still possible to use the stationary phase method to calculate level crossing probabilities in the three neutrino oscillations. Making convenient SU(3) transformations on the evolution matrix (9) it is possible to conciliate resonances and stationary phases. This is because resonances are invariant under similarity transformations, while stationary phases do not. Therefore the matrix
presents a stationary phase for the minimum of coinciding with the minimum of the squared mass difference , and it can be used to calculate the level crossing probability [8] around the lower resonance. To obtain the correct stationary phase to analyse the higher resonance, we rotate the evolution matrix given in Eq. (4) in the following way:
Now the minimum of the difference H33 - H11 indicates a stationary phase which now coincides with the required resonance.
V. Conclusions
Resonances represent a crucial region in the time evolution of neutrinos oscillating in matter. They are closely related with the nonadiabatic character of the oscillation. We investigated a general criterion to define a resonant region when three neutrino are present in the oscillation phenomenon. We observed that two of the three commonly employed criteria to identify a resonance in two neutrino oscillations are not valid anymore in its simplest extension to three neutrino MSW effect. For instance, mixing angles can present maximal values far from resonant regions and therefore this criterion to define a resonance has to be abandoned. Furthermore, stationary phases do not necessarily coincides with resonant regions. The safest way to identify such resonance regions is to investigate the behavior of the squared matter eigenvalue differences, looking for their minimum values.
Acknowledgements
The authors would like to thank Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for several financial supports.
References
[1] V. Barger, R.J.N. Phillips and K. Whisnant, Phys. Rev. D43 (1991) 1110;
[2] J.N. Bahcall and R. Ulrich, Rev. Modern Physics 60, 297 (1988).
[3] S.P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 1913 (1985); Nuovo Cimento 9C, 17 (1986).
[4] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59, 671 (1987).
[5] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); 20, 2634 (1979).
[6] M.M. Guzzo and J. Bellandi, Phys. Lett. B 294, 243 (1992).
[7] A. Erdélyi, Asymptotic Expansions, Dover Pub. Inc. NY (1956).
[8] M.M. Guzzo, J. Bellandi, and V.M. Aquino, Phys. D 49, 1404 (1994).
[9] C. Zener, Proc. R. Soc. London A 137, 696 (1932); L. Landau, Phys. Z. Sow. 2, 46 (1932).
[10] S.T. Petcov, Phys. Lett. B 191, 299 (1987).
[11] T.K. Kuo and James Pantaleone, Phys. Rev. D 35, 3432 (1987).
[12] V. Barger, K. Whisnant, S. Pakvasa, and R.J.N. Phillips, Phys. Rev. D 22, 2718 (1980).
[13] R.P. Feynman, Phys. Rev. 84, 108 (1951).
[14] Harry Bateman, Tables of Integral Transforms, Mc Graw-Hill, Book Company Inc. (1954).
[15] H.W. Zaglauer and K.H. Schwarzer, Phys. Lett. B 198, 556 (1987); Z. Phys. C 40, 273 (1988).
[16] Some physical consequences of this fact have been analysed in M.M. Guzzo, O.L.G. Peres, V. Pleitez and R. Zukanovich Funchal, "Nonresonant MSW solution to the solar neutrino problem'', preprint IFT-P.003/95, IFUSP-P/1133.
- [1] V. Barger, R.J.N. Phillips and K. Whisnant, Phys. Rev. D43 (1991) 1110;
- [2] J.N. Bahcall and R. Ulrich, Rev. Modern Physics 60, 297 (1988).
- [3] S.P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 1913 (1985);
- [4] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59, 671 (1987).
- [5] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978);
- [6] M.M. Guzzo and J. Bellandi, Phys. Lett. B 294, 243 (1992).
- [7] A. Erdélyi, Asymptotic Expansions, Dover Pub. Inc. NY (1956).
- [8] M.M. Guzzo, J. Bellandi, and V.M. Aquino, Phys. D 49, 1404 (1994).
- [9] C. Zener, Proc. R. Soc. London A 137, 696 (1932);
- L. Landau, Phys. Z. Sow. 2, 46 (1932).
- [10] S.T. Petcov, Phys. Lett. B 191, 299 (1987).
- [11] T.K. Kuo and James Pantaleone, Phys. Rev. D 35, 3432 (1987).
- [12] V. Barger, K. Whisnant, S. Pakvasa, and R.J.N. Phillips, Phys. Rev. D 22, 2718 (1980).
- [13] R.P. Feynman, Phys. Rev. 84, 108 (1951).
- [14] Harry Bateman, Tables of Integral Transforms, Mc Graw-Hill, Book Company Inc. (1954).
- [15] H.W. Zaglauer and K.H. Schwarzer, Phys. Lett. B 198, 556 (1987);
Publication Dates
-
Publication in this collection
04 Feb 1999 -
Date of issue
Sept 1997
History
-
Received
20 May 1996