arxiv: v2 [hep-ph] 12 May 2007

Transcription

1 On the exact formula for neutrino oscillation probability by Kimura Takamura and Yokomakura Osamu Yasuda arxiv:74.5v [hep-ph] May 7 Department of Physics Tokyo Metropolitan University Minami-Osawa Hachioi Tokyo 9-97 Japan Abstract The exact formula for the neutrino oscillation probability in matter with constant density which was discovered by Kimura Takamura and Yokomakura has been applied mostly to the standard case with three flavor neutrino so far. In this paper applications of their formula to more general cases are discussed. It is shown that this formalism can be generalized to various cases where the matter potential have off-diagonal components and the two non-trivial examples are given: the case with magnetic moments and a magnetic field and the case with non-standard interactions. It is pointed out that their formalism can be applied also to the case in the long baseline limit with matter whose density varies adiabatically as in the case of solar neutrino. PACS numbers: 4.6.Pq 4.6.St Electronic address: yasuda at phys.metro-u.ac.p Typeset by REVTEX

2 I. INTRODUCTION Neutrino oscillations in matter (See e.g. Ref. [] for review.) have been discussed by many people in the past because the oscillation probability has non-trivial behaviors in matter and due to the matter effect it may exhibit non-trivial enhancement which could be physically important. Unfortunately it is not easy to get an analytical formula for the oscillation probability in the three flavor neutrino scheme in matter and investigation of its behaviors has been a difficult but important problem in the phenomenology of neutrino oscillations. In Kimura Takamura and Yokomakura derived a nice compact formula [ ] for the neutrino oscillation probability in matter with constant density. Basically what they showed is that the quantity Ũ αũβ which is a factor crucial to express the oscillation probability analytically can be expressed as a linear combination of U α U β where Ũα and U α stand for the matrix element of the MNS matrix in matter and in vacuum respectively. However their formula is only applicable to the standard three flavor case. In this paper we show that their result can be generalized to various cases. We also show that their formalism can be applied also to the case with slowly varying matter density in the limit of the long neutrino path. In Sect. II we review briefly some aspects of the oscillation probabilities including a simple derivation for the formula by Kimura Takamura and Yokomakura which was given in Ref. [4] because these are used in the following sections. Their formalism is generalized to the various cases where the matter potential has off-diagonal components and we will discuss the case with large magnetic moments and a magnetic field (Sect. III) and the case with non-standard interactions (Sect. IV). In Sect. V we summarize our conclusions. II. GENERALITIES ABOUT OSCILLATION PROBABILITIES A. The case of constant density It has been known [5] (See also earlier works [6 7 8].) that after eliminating the negative energy states by a Tani-Foldy-Wouthusen-type transformation the Dirac equation for neutrinos propagating in matter is reduced to the familiar form: where E diag(e E E ) i dψ dt = [ UEU +A(t) ] Ψ () A(t) G F diag(n e (t) N n (t)/ N n (t)/ N n (t)/) Ψ T (ν e ν µ ν τ ) is the flavor eigenstate U is the Maki-Nakagawa-Sakata (MNS) matrix E m + p ( = ) is the energy eigenvalue of each mass eigenstate and the matter effect A(t) at time (or position ) t is characterized by the density N e (t) of electrons and the one N n (t) of neutrons respectively. Throughout this paper we assume for simplicity that the density of matter is either constant or slowly varying so that its derivative is negligible. The matrix on the right hand side of Eq. () can be formally diagonalized as: UEU +A(t) = Ũ(t)Ẽ(t)Ũ (t) ()

3 where Ẽ(t) diag ( Ẽ (t)ẽ(t)ẽ(t) ) is a diagonal matrix with the energy eigenvalues Ẽ(t) in the presence of the matter effect. First of all let us assume that the matter density A(t) is constant. Then all the t dependence disappears and Eq. () can be easily solved resulting the flavor eigenstate at the distance L: Ψ(L) = Ũ exp( iẽl) Ũ Ψ(). () Thus the oscillation probability P(ν α ν β ) is given by P(ν α ν β ) = [Ũ exp( iel) Ũ ] βα = δ 4 Re ( ) ) k sin ( Ẽ k L <k + Im ( ) ( k sin Ẽ k L ) (4) <k where we have defined ŨαŨ β Ẽk Ẽ Ẽk and throughout this paper the indices αβ = (eµτ) and k = () stand for those of the flavor and mass eigenstates respectively. Once we know the eigenvalues Ẽ and the quantity the oscillation probability can be expressed analytically. B. The case of adiabatically varying density Secondly let us consider the case where the density of the matter varies adiabatically as in the case of the solar neutrino deficit phenomena. In this case instead of Eq. () we get ] Ψ(L) = Ũ(L)exp [ i L Ẽ(t) dt Ũ() Ψ() where Ũ() and Ũ(L) stand for the effective mixing matrices at the origin t = and at the end point t = L. The oscillation probability is given by [ { ] L P(ν α ν β ) = Ũ(L)exp i Ẽ(t) dt }Ũ() βα [ ] = Ũ(L) β Ũ(L) βkũ() αũ() αkexp k i L Ẽ(t) kdt. (5) InthestandardcasewiththreeflavorsofneutrinosinmattertheenergyeigenvaluesẼ canbeanalytically obtained by the root formula for a cubic equation [9]. So the only non-trivial problem in the standard case is to obtain the expression for and this was done by Kimura Takamura and Yokomakura [ ]. In general cases however the analytic expression for Ẽ is very difficult or impossible to obtain and we will discuss below only examples in which the analytic expression for Ẽ is known.

4 Eq. (5) requires in general the quantity like Ũ(t) βũ (t) βk which has the same flavor index β but different mass eigenstate indices k and it turns out that the analytical expression for Ũ(t) βũ (t) βk is hard to obtain. However if the length L of the neutrino path is very large and if L Ẽ(t) kdt is satisfied for k as in the case of the solar neutrino deficit phenomena after averaging over rapid oscillations Eq. (5) is reduced to P(ν α ν β ) = ββ (L) () where we have defined (t) Ũ(t) α. In the case of the solar neutrinos deficit process ν e ν e during the daylight ββ (L) at the end point t = L and () at the origin t = correspond to X ββ in vacuum and [ ] at the center of the Sun respectively where X [ ] U α U β [ Ũ α Ũ β are bilinear products of the elements of the mixing matrices in vacuum and at the center of the Sun respectively. Thus we obtain P(ν e ν e ) = X ee ] [ ] ee. Hence we see that evaluation of the quantity in the presence of the matter effect is important notonlyinthecase ofconstant matterdensity butalso inthecase ofadiabatically varying density. C. Another derivation of the formula by Kimura Takamura and Yokomakura In this subsection a systematic derivation of their formula is given because such a derivation will be crucial for the generalizations in the following sections. The arguments are based on the trivial identities. From the unitarity condition of the matrix Ũ we have δ = [ ŨŨ ] = Ũ α Ũ β =. (6) Next we take the (αβ) component of the both hand sides in Eq. (): [ UEU +A ] = [ ŨẼŨ ] = Ũ α Ẽ Ũ β = Ẽ (7) The argument here is the same as that in Ref. [4]. Since this derivation does not seem to be widely known it is reviewed here. 4

5 Furthermore we take the (αβ) component of the square of Eq. (): [ (UEU +A ) ] = [ ŨẼ Ũ ] = Ũ α Ẽ Ũ β = Ẽ (8) Putting Eqs. (6) (8) together we have Ẽ Ẽ Ẽ Ẽ Ẽ Ẽ = δ [UEU +A] [(UEU +A) ] which can be easily solved by inverting the Vandermonde matrix: = (ẼẼ (Ẽ +Ẽ) ) Ẽ Ẽ (ẼẼ (Ẽ +Ẽ) ) Ẽ Ẽ (ẼẼ (Ẽ +Ẽ) ) Ẽ Ẽ δ [UEU +A] [(UEU +A) ]. (9) [(UEU +A) ] ( = ) on the right hand side are given by the known quantities: [ UEU +A ] = [ (UEU +A ) ] = E X +Aδ αe δ βe E X +A ( E δαe X eβ ) +δ βe X αe +A δ αe δ βe. It can be shown that Eq. (9) coincides with the original results by Kimura Takamura and Yokomakura [ ]. A remark is in order on Eq. (9). Addition of a matrix c to Eq. () where c is a constant and is the identity matrix or in other words the shift E E +c ( = ) () should give the same result for ( = ) since Eq. () only affects the overall phase of the oscillation amplitude and the phase has to disappear in the probability. It is easy to show that the shift () indeed gives the same result as Eq. (9). The proof is given in Appendix A. In practical calculations below we will always put c = E i.e. we will consider the mass matrix U(E E )U +A instead of the original one UEU +A since all the diagonal elements (E E ) = E = m /E are expressed in terms of the relevant variables m and therefore calculations become simpler. To save space however we will use the matrix UEU +A in most of the following discussions. D. The case with arbitrary number of neutrinos It is straightforward to generalize the discussions in sect. II C to the case with arbitrary number of neutrinos where the matter potential is diagonal in the flavor eigenstate. The 5

6 scheme with number of sterile neutrinos is one of the example of these cases [4 ]. The time evolution of such a scheme with N neutrino flavors is described by i dψ N dt = ( U N E N U N +A N) ΨN where Ψ T N (ν α ν α ν αn ) is the flavor eigenstate E N diag(e E E N ) () is the energy matrix of the mass eigenstate A N diag(a A A N ) is the potential matrix for the flavor eigenstate and U N is the N N MNS matrix. As in the previous sect. by taking the αβ components we get Ẽ m = [( U N E N UN +A ) m ] N for m = N which leads to the simultaneous equation. Ẽ Ẽ Ẽ N Ẽ N.. Ẽ N ẼN N. N = δ [UN EU N +A N ] [. (UN ) ] EUN N +A N. () Eq. () can be solved by inverting the N N Vandermonde matrix V N :. N = V N δ [UN E N U N +A N ] [. (UN E N UN +A N ) N ]. () The determinant of V N is the Vandermonde determinant <k Ẽk and therefore V can be analytically obtained as long as we know the value of Ẽ. The factors [(U N E N UN + A N ) ] ontherighthandsideofeq. ()canbeexpressed asfunctionsoftheenergye the quantityx invacuumandthematterpotentiala γ sincethematrix(u N E N UN +A N) isa sum of products of the matrices [(U N E N UN ) l ] γδ = kex l γδ k ( l ) and [(A N ) m ] ǫη = (A ǫ ) m δ ǫη ( m ). From Eq. () it is clear that enhancement of the oscillation probability due to the matter effect occurs only when some of Ẽk becomes small. 6

7 III. THE CASE WITH LARGE MAGNETIC MOMENTS AND A MAGNETIC FIELD So far we have assumed that the potential term is diagonal in the flavor basis. We can generalize the present result to the cases where we have off-diagonal potential terms. One of such examples is the case where there are only three active neutrinos with magnetic moments and the magnetic field (See e.g. Ref. [] for review.). The hermitian matrix with ( ) UEU B M B U E(U ) B Bµ is the mass matrix for neutrinos and anti-neutrinos without the matter effect where neutrinos have the magnetic moments µ in the magnetic field B. Here we assume the magnetic interaction of Maorana type (4) µ ν α F λκ σ λκ νβ c +h.c. (5) and in this case the magnetic moments µ are real and anti-symmetric in flavor indices: µ = µ βα. If the magnetic field is constant then the oscillation probability can be written as P(ν A ν B ) = δ AB 4 + J<K J<K Im ( AB J Re ( AB J AB K ) ) K AB sin ( Ẽ JK L ) ( sin Ẽ JK L ) (6) where AB run eµτē µ τ and JK run 6 respectively and U AJ UBJ. Ẽ J (J = 6) are the eigenvalues of the 6 6 matrix M. On the other hand if the magnetic field varies very slowly and if the length L of the baseline is so long that ẼJKL is satisfied for J K then the oscillation probability is given by AB J P(ν A ν B ) = 6 J= J BB AA (L) J (). (7) Following the same arguments as before the quantity Vandermonde matrix V 6 : AB AB. AB 6 = V 6 δ AB [M] AB [. (M) 5 ] AB AB J is given by inverting the 6 6. (8) See [5] for derivation of Eq. (4) from the Dirac Eq. 7

8 As in the previous sections [(M) J ] AB (J = 5) on the right hand side of Eq. (8) can be expressed in terms of the known quantities XK AB and B CD and Eqs. (6) and (8) are useful only when we know the eigenvalues ẼJ. To demonstrate the usefulness of these formulae let us consider the case where the magnetic field is large at origin but is zero at the end point and the magnetic field varies adiabatically. For simplicity we assume that θ and all the CP phases vanish. 4 In this case the 6 6 matrix M in Eq. (4) becomes real and we obtain the following oscillation probabilities: P(ν α ν β ) = P( ν α ν β ) = (U β ) [ReŨ() α] = P(ν α ν β ) = P( ν α ν β ) = (U β ) [ImŨ() α] (9) = where Ũ() the unitary matrix which diagonalizes the matrix UEU +ib() at the origin: UEU +ib() = Ũ()Ẽ()Ũ (). In this example the energy eigenvalues are degenerate i.e. the 6 6 energy matrix becomes diag(ẽẽ) and the oscillation probability differs from Eq. (7) because the condition ẼJKL (J K) is not satisfied (e.g. ẼJK = not only for J = K = but also for J = K = 4). Each probability in Eqs. (9) itself is not expressed in terms of () but we find that the following relation holds: P(ν α ν β )+P( ν α ν β ) = (U β ) Ũ() α = = = X ββ (). () Eq. () is a new result and without the present formalism it would be hard to derive it. The details of derivation of Eq. (9) and explicit forms of () are given in Appendix B. Eq. () may be applicable to the case where high energy astrophysical neutrinos which are produced in a relatively large magnetic field are observed on the Earth on the assumption that the fluxes of neutrinos and anti-neutrinos are almost equal. IV. THE CASE WITH NON-STANDARD INTERACTIONS Another interesting application is the oscillation probability in the presence of new physics in propagation [ ]. In this case the mass matrix is given by UEU +A NP () 4 In the presence of the magnetic interaction(5) ofmaoranatype the twocp phases which are absorbed by redefinition of the charged lepton fields in the standard case cannot be absorbed and therefore become physical. Here however we will assume for simplicity that these CP phases vanish. 8

9 where A NP G F N e +ǫ ee ǫ eµ ǫ eτ ǫ eµ ǫ µµ ǫ µτ ǫ eτ ǫ µτ ǫ ττ. Thedimensionless quantitiesǫ standforpossibledeviationfromthestandardmattereffect. Also in this case the oscillation probability is given by Eqs. (4) and (9) where the standard potential matrix A has to be replaced by A NP. The extra complication compared to the standard case is calculations of the eigenvalues Ẽ and the elements [(UEU +A NP ) m ] (m = ). Again to demonstrate the usefulness of the formalism here we will discuss for simplicity the case in which the eigenvalues are the roots of a quadratic equation. It is known [] that the constraints on the three parameters ǫ ee ǫ eτ ǫ ττ from various experimental data are weak and they could be as large as O(). In Ref. [4] it was found that large values ( O()) of the parameters ǫ ee ǫ eτ ǫ ττ are consistent with all the experimental data including those of the atmospheric neutrino data provided that one of the eigenvalues of the matrix () at high energy limit i.e. A NP becomes zero. Simplifying even further here we will neglect the parameters ǫ eµ ǫ µµ ǫ µτ which are smaller than O( ) and we will consider the potential matrix +ǫ ee ǫ eτ A NP = A ǫ eτ ǫ ττ () where A G F N e the three parameters ǫ ee ǫ eτ ǫ ττ are constrained in such a way that two of the three eigenvalues become zero. We will assume that N e is constant and we will take the limit m. The oscillation probability P(ν µ ν e ) in this case can be analytically expressed and is given by P(ν µ ν e ) = 4Re ( µe 4Re ( µe µe µe ) ) sin ( Λ L 4Re ( µe ] ) sin [ (Λ+ Λ L)L µe ) ) sin ( Λ+ L + 8A( E ) Λ + Λ (Λ + Λ ) ǫ eτx eµ X µτ sin(arg(ǫ eµ )+δ) ( ) ( ) [ ] Λ L Λ+ L (Λ+ Λ )L sin sin sin. () Eq. () is another new result and it would be difficult to obtain it without using the present formalism. The details of derivation of Eq. () explanation of the notations and the explicit forms of all the variables in Eq. () are described in Appendix C. V. CONCLUSIONS The essence of the exact formula for the neutrino oscillation probability in constant matter which was discovered by Kimura Takamura and Yokomakura lies in the fact that the combination Ũα Ũ β of the mixing matrix elements in matter can be expressed 9

10 as polynomials in the same quantity X U α U β in vacuum. In this paper we have discussed applications of their formalism to more general cases. We have pointed out that their formalism can be useful for the cases in matter not only with constant density but also with density which varies adiabatically as in the case of the solar neutrino problem after taking the limit of the long neutrino path. We have shown that their formalism can be generalized to the cases where the matter potential has off-diagonal components. As concrete non-trivial examples we discussed the case with magnetic moments and a magnetic field and the case with non-standard interactions. The application of the present formalism to the case with unitarity violation has been discussed elsewhere [5]. The formalism by Kimura Takamura and Yokomakura is quite general and can be applicable to many problems in neutrino oscillation phenomenology. APPENDIX A: PROOF THAT EQ. () GIVES THE SAME (9) In this appendix we show that Eq. () gives the same result for value of ( = ) for Ũ ( Ẽ +c ) Ũ = UEU +A+c ( = ). The becomes at most quadratic 5 in c and all one has to do is to show that the coefficients of the terms linear and quadratic in c vanish. Let us introduce the notation Ẽ +c Ẽ +c Ẽ +c (Ẽ +c) (Ẽ +c) (Ẽ +c) δ [UEU +A+c] [(UEU +A+c) ] (V ) () +c(v ) () +c (V ) () B () +c B () +c B () where V (k) is the coefficient of the inverted Vandermonde matrix which is k-th order in c and B (k) is the coefficient of the vector (UEU +A+c) which is k-th order in c. Then the terms linear in c are given by = (V ) () B () +(V ) () B () (Ẽ +Ẽ ) Ẽ Ẽ ( (Ẽ +Ẽ) + ) Ẽ Ẽ (Ẽ +Ẽ ) Ẽ Ẽ δ [UEU +A] [(UEU +A) ] 5 Notice that all the factors Ẽk are invariant under the shift () and the only change by this shift comes either from the terms ẼẼk or from Ẽ +Ẽk in the inverse of the Vandermonde matrix (cf. Eq. (9)). Hence the difference by Eq. () is at most quadratic in c.

11 + Ẽ Ẽ Ẽ Ẽ Ẽ Ẽ and the terms quadratic in c are given by (+ẼẼ (Ẽ +Ẽ) +) ( ẼẼ +(Ẽ +Ẽ) ) δ = [UEU +A] (+ẼẼ (Ẽ +Ẽ) +) = (V ) () B () +(V ) () B () +(V ) () B () (+ ) Ẽ Ẽ δ ( ) [UEU [(UEU +A] Ẽ Ẽ +A) ] (+ ) Ẽ Ẽ (Ẽ +Ẽ ) Ẽ Ẽ + ( (Ẽ +Ẽ) + ) δ Ẽ Ẽ [UEU +A] (Ẽ +Ẽ ) Ẽ Ẽ (+ẼẼ (Ẽ +Ẽ) +) Ẽ Ẽ + ( ẼẼ +(Ẽ +Ẽ) ) =. Ẽ Ẽ δ (+ẼẼ (Ẽ +Ẽ) +) Ẽ Ẽ Thus ( = ) is independent of c as is claimed. APPENDIX B: DERIVATION OF EQ. (9) The matrix (4) can be rewritten as M = ( )( i UEU +ib i UEU ib )( ) i i so the problem of diagonalizing the 6 6 matrix (4) is reduced to diagonalizing the matrices UEU ± ib. Since we are assuming that θ and all the CP phases vanish all the matrix elements U α and B = B βα are real UEU ±ib can be diagonalized by a unitary matrix and its complex conugate: UEU +ib = ŨẼŨ UEU ib = Ũ Ẽ(Ũ ). Therefore we can diagonalize M by a 6 6 unitary matrix Ũ as (Ẽ ) M = Ũ Ẽ Ũ

12 where Ũ = ( )(Ũ ) i = ( i Ũ ) Ũ iũ iũ Ũ. We note in passing that the reason why diagonalization of the 6 6 matrix is reduced to that of the matrix is because the two matrices UEU and B are real. On the other hand without a magnetic field the 6 6 unitary matrix U is given by U = ( U U ) = ( U U where the CP phase δ has dropped out because θ =. From these we can integrate the equation of motion and we get the fields at the end point: ( ) ( ) ( ) Ψ(L) e Ψ c = (L) Ũ(L) iφ Ψ() e iφ Ũ() Ψ c () = ( )( ) Ue iφ Ũ +U e iφ (Ũ ) i(ue iφ Ũ U e iφ (Ũ ) ) Ψ() i(ue iφ Ũ U e iφ (Ũ ) ) Ue iφ Ũ +U e iφ (Ũ ) Ψ c () where Φ L Ẽ(t) dt and we have assumed that a large magnetic field exists at the origin whereas there is no magnetic field at the end point. Thus the oscillation probabilities for the adiabatic transition are give by: [ P(ν α ν β ) = P( ν α ν β ) = lim L Ue iφ Ũ +U e iφ (Ũ ) ] = U β [ Re(Ũα) ] = L P( ν α ν β ) = P(ν α ν β ) = lim Hence we obtain the following relation: = = ) [ Ue iφ Ũ U e iφ (Ũ ) ] U β [ Im(Ũα) ]. P(ν α ν β )+P( ν α ν β ) = P(ν α ν β )+P(ν α ν β ) = U β Ũα. = To get Ũα we need the explicit expression for the eigenvalues and the quantity in the presence of a magnetic field. In the following we will subtract E from the energy matrix E because it will only change the phase of the oscillation amplitude. For simplicity we will put θ = θ = π/4 and we will consider the limit m. Defining E k m k /E and B = Bµ p q p r q r

13 we have the eigenvalue equation = λ U(E E )U ib = λ E λ (p +q +r )λ+ E (p q). The three roots of the cubic equation (B) are given by (B) λ = Rcosϕ+ E where λ = Rcos(ϕ+ π)+ E R [( E /) +(p +q +r )/] / λ = Rcos(ϕ π)+ E ϕ (/)cos [ {( E /) + E (p +q +r )/6 E (p q) /4}/R ]. The quantity in the presence of a magnetic field is given by = λ λ (λ λ (λ +λ ) ) λ λ (λ λ (λ +λ ) ) λ λ (λ λ (λ +λ ) ) Y Y (B) where Y = [ U(E E )U +ib ] = E X { (α = e) = E / (α = µτ) [ {U(E Y = E )U +ib } ] = q +r (α = e) r +p +( E ) / (α = µ) p +q +( E ) / (α = τ) = ( E ) X (B ) (B). (B4) In evaluating Y we have used the facts θ = θ = π/4 E = B = B βα and that U(E E )U is a symmetric matrix. Using all these results it is straightforward to obtain the explicit form for P(ν α ν β )+P( ν α ν β ) by plugging the results of Eqs. (B) (B) (B4) into the following (although calculations are tedious): P(ν α ν e )+P( ν α ν e ) = c P(ν α ν β )+P( ν α ν β ) = c where s sinθ c cosθ. +s + s + (β = µτ)

14 APPENDIX C: DERIVATION OF EQ. () The oscillation probability () is obtained in two steps. First we will obtain the eigenvalues of thematrix() witheq. () andthenwe will plug the expressions for theeigenvalues into Eq. (9) with A replaced by A NP given in Eq. (). Let us introduce notations for hermitian matrices: λ λ i i λ 5 λ 9 i i λ 7 i i whereλ λ 5 andλ 7 arethestandardgell-mannmatriceswhereasλ andλ 9 arethenotations which are defined only in this paper. Simple calculations show that the matrix A NP in Eq. () can be rewritten as where A NP = Ae iγλ 9 e iβλ 5 λ +ǫ ee +ǫ ττ (+ǫee ) ǫ ττ +λ 9 + ǫ µτ e iβλ 5 e iγλ 9 (C) β tan ǫ eτ +ǫ ee ǫ ττ γ arg(ǫ eµ). From Eq. (C) we see that the two potentially non-zero eigenvalues λ e and λ τ of the matrix () are given by ( λe λ τ ) = A +ǫ ee +ǫ ττ (+ǫee ) ǫ ττ ± + ǫ µτ. In order for this scheme to be consistent with the atmospheric neutrino data particularly at high energy which are perfectly described by vacuum oscillations λ τ has to vanish [4]. In this case we have Thus we have tanβ = ǫ eτ +ǫ ee ǫ ττ = ǫ eτ +ǫ ee λ e = A(+ǫ ee ) [ + ǫ eτ ] (+ǫ ee ) = A(+ǫ ee) cos β. A NP = Ae iγλ 9 e iβλ 5 diag(λ e )e iβλ 5 e iγλ 9. (C) 4

15 If we did not have β and γ Eq. (C) would be the same as the standard three flavor scheme in matter which was analytically worked out in Ref. [6] in the limit of m. It turns out that by redefining the parametrization of the MNS matrix Eq. (C) can be also treated analytically in the limit of m as was done in Ref. [6]. The mass matrix can be written as [ UEU +A NP = e iγλ 9 e iβλ 5 e iβλ 5 e iγλ 9 UEU e iγλ 9 e iβλ 5 +diag(λ e ) ] e iβλ 5 e iγλ 9. Here we introduce the following two unitary matrices: U e iβλ 5 e iγλ 9 U diag(e iargu τ )U diag(e iargu e e iargu e ) where U is the MNS matrix in the standard parametrization [7] and U was defined in the second line in such a way that the elements U e U e U τ be real to be consistent with the standard parametrization in Ref. [7] 6. Then we have [ UEU +A NP = e iγλ 9 e iβλ 5 diag(e iargu τ ) U EU +diag(λ e ) ] diag(e iargu τ )e iβλ 5 e iγλ 9. Before proceeding further let us obtain the expression for the three mixing angles θ k and the Dirac phase δ in U. Since U = c β e iγ U e +s β e iγ U τ c β e iγ U e +s β e iγ U τ c β e iγ U e +s β e iγ U τ U µ U µ U µ c β e iγ U τ s β e iγ U e c β e iγ U τ s β e iγ U e c β e iγ U τ s β e iγ U e where c β cosβ s β sinβ we get θ = sin U e = sin c β e iγ U e +s β e iγ U τ θ = tan (U e/u e) = tan ( c β e iγ U e +s β e iγ U τ / c β e iγ U e +s β e iγ U τ ) θ = tan (U µ /U τ ) = tan ( U µ / c β e iγ U τ s β e iγ U e ) δ = argu e = arg(c βe iγ U e +s β e iγ U τ ). As was shown in Ref. [6] in the limit m the matrix on the right hand side of Eq. (C) can be diagonalized as follows: U EU +diag(λ e ) E = e iθ λ 7 Γ δ e iθ [ e iθ = e iθ λ 7 Γ δ λ 5 Γ δ eiθ λ diag( E )e iθ λ Γ δ e iθ λ 5 Γ δ e iθ λ 5 diag( E )+diag(λ e ) ] Γ δ e iθ λ 7 = e iθ λ 7 i θ Γ δ e λ 5 i θ diag(λ Λ + )e λ 5 Γ δ eiθ λ (C) λ 7 +diag(λ e ) where Γ δ diag(e iδ ) E m /E we have used the standard parametrization [7] U e iθ λ 7 Γ δ e iθ λ 5 Γ eiθ λ and the eigenvalues Λ ± are defined by Λ ± = ( E +λ e )± δ ( E cosθ λ e ) +( E sinθ ). 6 The element U τ has to be also real but it is already satisfied because U τ = U τ. 5

16 Having obtained the eigenvalues by plugging these into Eq. (9) with A A NP Ẽ Λ Ẽ Ẽ Λ + we obtain µe : µe µe µe = Λ (Λ + Λ ) ( Λ + ) Λ + Λ ( Λ + Λ (Λ + +Λ ) ) Λ + (Λ + Λ ) ( Λ ) Y µe Y µe = Y µe +Λ + Y µe Λ (Λ + Λ ) Y µe (Λ + +Λ )Y µe Λ + Λ Y µe Λ Y µe Λ + (Λ + Λ ) where Y µe are defined by [ (UEU ] Y µe +A NP) µe and are given by Y µe = E X µe Y µe = [( E ) +A(+ǫ ee ) E ]X µe +A E ǫ eτ Xµτ. Furthermore by introducing the notations ξ [( E ) +A(+ǫ ee ) E ]U µ U e η A E ǫ eτ U µ U τ ζ E U µ U e we can rewrite Y µe = ζe iδ and Y µe = ξe iδ + ηe iγ where δ is the Dirac CP phase of the MNS matrix U so we have µe = e iδ Λ (Λ + Λ ) [ξ +ηe i(γ+δ) Λ + ζ] µe = eiδ Λ + Λ [ξ +ηe i(γ+δ) (Λ + +Λ )ζ] µe = e iδ Λ + (Λ + Λ ) [ξ +ηe i(γ+δ) Λ ζ]. Notice that the phase factor e iδ in front of each drops out in the oscillation probability P(ν µ ν e ) because P(ν µ ν e ) is expressed in terms of µe µe k and the oscillation probability () depends only on the combination γ +δ = arg(ǫ eµ )+δ. In the present case the matrix Ũ is unitary and because of this three flavor unitarity all the T violating terms are proportional to one factor: Im ( µe <k = Im ( µe = 8Im ( µe µe ) ( k sin Ẽ k L ) µe µe ) ( [sin Ẽ L ) sin ( ẼL ) +sin ( ẼL ) ] ( ) ( ) ( ) ) Ẽ L Ẽ L Ẽ L sin µe 6 sin sin.

17 µe This modified Jarlskog factor Im( µe ) in matter can be rewritten as Im( µe µe ) = This completes derivation of Eq. (). µe Im(Y Y µe ) = Λ + Λ (Λ + Λ ) ηζsin(γ +δ) Λ + Λ (Λ + Λ ) = A( E ) Λ + Λ (Λ + Λ ) ǫ eτx eµ X µτ sin(arg(ǫ eµ )+δ). ACKNOWLEDGMENTS The author would like to thank Alexei Smirnov for bringing my attention to Refs. [ ]. He would also like to thank He Zhang for calling my attention to Refs. [4 ] which were missed in the first version of this paper. This research was supported in part by a Grant-in-Aid for Scientific Research of the Ministry of Education Science and Culture #946. [] J. N. Bahcall R. Davis P. Parker A. Smirnov and R. Ulrich Reading USA: Addison-Wesley (995) 44 p. (Frontiers in physics. 9) [] K. Kimura A. Takamura and H. Yokomakura Phys. Lett. B () [arxiv:hep-ph/99]. [] K. Kimura A. Takamura and H. Yokomakura Phys. Rev. D () [arxiv:hep-ph/595]. [4] Z. z. Xing and H. Zhang Phys. Lett. B 68 (5) [arxiv:hep-ph/58]. [5] W. Grimus and T. Scharnagl Mod. Phys. Lett. A 8 94 (99). [6] A. Halprin Phys. Rev. D 4 46 (986). [7] P. D. Mannheim Phys. Rev. D 7 95 (988). [8] R. F. Sawyer Phys. Rev. D 4 98 (99). [9] V. D. Barger K. Whisnant S. Pakvasa and R. J. N. Phillips Phys. Rev. D 78 (98). [] H. Zhang arxiv:hep-ph/664. [] M. M. Guzzo A. Masiero and S. T. Petcov Phys. Lett. B 6 54 (99); [] E. Roulet Phys. Rev. D (99). [] S. Davidson C. Pena-Garay N. Rius and A. Santamaria JHEP () [arxiv:hep-ph/9]. [4] A. Friedland and C. Lunardini Phys. Rev. D 7 59 (5) [arxiv:hep-ph/564]. [5] E. Fernandez-Martinez M. B. Gavela J. Lopez-Pavon and O. Yasuda arxiv:hep-ph/798. [6] O. Yasuda Proceedings of Symposium on New Era in Neutrino Physics (Universal Academy Press Inc. Tokyo eds. H. Minakata and O. Yasuda) p (999) [arxiv:hep-ph/9895]. [7] S. Eidelman et al. [Particle Data Group] Phys. Lett. B 59 (4). 7