arxiv:hep-ph/ v4 1 Sep 2002

Transcription

1 DPNU-02-4 Exact Formulas and Simple CP dependence of Neutrino Oscillation Probabilities in Matter with Constant Density Keiichi Kimura,, Akira Takamura,,2, and Hidekazu Yokomakura, Department of Physics, Nagoya University, Nagoya, , Japan 2 Department of Mathematics, Toyota National Collage of Technology Eisei-cho 2-, Toyota-shi, , Japan arxiv:hep-ph/ v4 Sep 2002 We investigate neutrino oscillations in constant matter within the context of the standard three neutrino scenario. We derive an exact and simple formula for the oscillation probability applicable to all channels. In the standard parametrization, the probability for ν e ν µ transition can be written in the form P(ν e ν µ) = A eµcosδ+b eµsinδ+c eµ without any approximation using CP phase δ. For ν µ ν τ transition, the linear term of cos2δ is added and the probability can be written in the form P(ν µ ν τ) = A µτ cosδ+b µτ sinδ+c µτ +D µτ cos2δ. We give the CP dependences of the probability for other channels. We show that the probability for each channel in matter has the same form with respect to δ as in vacuum. It means that matter effects just modify the coefficients A, B, C and D. We also give the exact expression of the coefficients for each channel. Furthermore, we show that our results with respect to CP dependences are reproduced from the effective mixing angles and the effective CP phase calculated by Zaglauer and Schwarzer. Through the calculation, a new identity is obtained by dividing the Naumov-Harrison-Scott identity by the Toshev identity. PACS numbers: 4.60.Pq I. INTRODUCTION Strong evidences for neutrino oscillations have been accumulated and they indicate the finite masses and mixings of neutrinos. The first evidence is the observation of zenith-angle dependence of the atmospheric neutrino deficit [], which is consistent with ν µ ν τ transition with the mass difference and the mixing as ev 2, sin 2 2θ 23, () where ij = m 2 i m2 j and m i is neutrino mass. The second evidence is the solar neutrino deficit [2], which is consistentwithν e ν µ /ν τ transition. Recentreportsfor the measurementsofcc and NC interactionsin the SNO experiment [3] are consistent with the Standard Solar Model and strongly suggest the LMA MSW solution ev 2, sin 2 2θ (2) to the solar neutrino problem. Thus, the information on neutrino masses and mixings (Maki-Nakagawa-Sakata matrix) [4] has been increased considerably. The main physics goals in future experiments are the determination of the unknown parameters θ 3 (only upper bound sin 2 2θ 3 < 0. is obtained [5]) and the CP phase δ in addition to the high precision measurement of the parameters in () and (2). In particular, the observation of δ is quite interesting from the point of view that δ is related to the origin of the matter in the universe. In order to get at the CP phase δ, a considerable number of studies have been made on CP violation [6, 7, 8, 9, 0]. address: kimukei@eken.phys.nagoya-u.ac.jp address: takamura@eken.phys.nagoya-u.ac.jp address: yoko@eken.phys.nagoya-u.ac.jp Moreover several experiments using neutrino beam are planned [, 2]. In such experiments with high energy neutrino beam, matter effects cannot be neglected and prevent the precise measurement of pure CP effects. In this paper we derive an exact and simple formula of oscillation probabilities for all channels. Before giving our results, let us give the present status of exact formulas associated with oscillation probability in constant matter. The mass eigenvalues have been given by Barger et al. [3]. Zaglauer and Schwarzer [4] have calculated effective mixing angles and effective CP phase in the standard parametrization. Furthermore, the elements of the MNS matrix have been calculated by Xing [5] in parametrization independent way. However, these exact formulas obtained by diagonalizing the Hamiltonian, are complicated and it is not easy to obtain useful information on matter effects and the CP phase. Ohlsson and Snellman [6] have directly calculated the amplitude without diagonalizing the Hamiltonian although the CP phase is not considered (which is similar to our result as shown in Sec.III). Recently, by comparing the Hamiltonian in matter and that in vacuum, Naumov[7] and Harrison and Scott [8] have presented simple and elegant identity J = J, (3) where J is the Jarlskog factor [9] and tilde represents the quantities in matter. We have investigated the matter enhancement of J in Ref. [20] by using this identity. Furthermore, the enhancement of T-violating quantity P T = P(ν e ν µ ) P(ν µ ν e ) has been also investigated by Parke and Weiler [2] and Akhmedov et al. [22]. In our previous paper [23], we have found other identities on the product of the MNS matrix elements by comparing the Hamiltonian in matter and that in vacuum.

2 2 From these identities, we have presented the exact formula for P(ν e ν µ ) and we have definitely separated pure CP effects from matter effects. In the standard parametrization, we have also shown that the probability for ν e ν µ transition can be written in the form [33] [23] P(ν e ν µ ) = A eµ cosδ +B eµ sinδ +C eµ. (4) It has also been found that CP trajectory proposed by Minakata and Nunokawa [24] is exactly elliptic in biprobability space when δ changes from 0 to 2π. In this paper, extending our previous results [23], we derive an exact and simple formula of the oscillation probabilities for all channels. We also find that the probability for each channel in matter has the same form with respect to the CP phase δ as in vacuum. For example, the probability for ν µ ν τ transition can be written as P(ν µ ν τ ) = A µτ cosδ +B µτ sinδ +C µτ +D µτ cos2δ, (5) both in vacuum and in matter. Namely, matter effects just modify A µτ, B µτ, C µτ and D µτ. We also give the exact expression of these coefficients. Furthermore, we show that our results with respect to CP dependences are reproduced from the effective mixing angles and effective CP phase calculated by Zaglauer and Schwarzer [4]. Although the effective CP phase are complicated, the CP dependence of probability is simplified for all channels. This is guaranteed by the Toshev identity s 23 c 23 sin δ = s 23 c 23 sin δ. (6) vacuum H = is diagonalized as H ee H eµ H eτ H µe H µµ H µτ H τe H τµ H ττ (8) U HU = 0 2, (9) 2E 3 where U is the MNS matrix which connects flavor eigenstates with mass eigenstates [34]. Then, the amplitude of ν α to ν β transition after transporting the distance L is A(ν α ν β ) = 3 i= U αi e i i 2E L U βi, (0) where α and β denote flavor indices e, µ or τ. And the oscillation probability is given by P(ν α ν β ) = A(ν α ν β ) 2 = δ αβ 4 ReJ ij αβ sin2 ij ±2 J sin2 ij,() Finally, we obtain a new identity s 2 c 2 s 3 c 2 3 = s 2 c 2 s 3 c 2 3, (7) related to the -2 mixing and the -3 mixing by dividing the Naumov-Harrison-Scott identity by the Toshev identity. where J ij αβ U αiu βi(u αj U βj), J ImJ 2 eµ, (2) ij ijl 4E (m2 i m2 j )L. (3) 4E II. CP DEPENDENCES OF OSCILLATION PROBABILITIES IN VACUUM In this section, we review the CP dependences of the oscillation probabilities in vacuum. The Hamiltonian in The ± sign of the third term takes (+) in the case that (α,β) is given by (anti) permutation of (e,µ). The sum is over = (2),(23),(3). Next let us calculate ReJeµ ij and J to study the CP dependences of probabilities in the standard parametrization [25], c 2 c 3 s 2 c 3 s 3 e iδ U αi = s 2 c 23 c 2 s 23 s 3 e iδ c 2 c 23 s 2 s 23 s 3 e iδ s 23 c 3 s 2 s 23 c 2 c 23 s 3 e iδ c 2 s 23 s 2 c 23 s 3 e iδ c 23 c 3, (4) where s ij sinθ ij, c ij cosθ ij. ReJ ij eµ are given by ReJ 2 eµ = (c 2 2 s 2 2)J r cosδ +s 2 2c 2 2c 2 3(s 2 23s 2 3 c 2 23), (5) ReJeµ 23 rcosδ s 2 2 s2 23 s2 3 c2 3, (6) ReJeµ 3 rcosδ c 2 2 s2 23 s2 3 c2 3, (7)

3 and ReJeτ ij ReJeτ 2 = (c 2 2 s2 2 )J rcosδ +s 2 2 c2 2 c2 3 (c2 23 s2 3 s2 23 ), (8) ReJeτ 23 = J r cosδ s 2 2c 2 23s 2 3c 2 3, (9) ReJeτ 3 = J r cosδ c 2 2 c2 23 s2 3 c2 3, (20) and ReJ ij µτ are given by ReJ 2 µτ = s 2 2 c2 2 s2 23 c2 23 (+s2 3 +s4 3 ) (s2 2 c2 2 +s2 23 c2 23 )s2 3 (c 2 2 s2 2 )(c2 23 s2 23 )s 2c 2 s 23 c 23 s 3 (+s 2 3 )cosδ +2s2 2 c2 2 s2 23 c2 23 s2 3cos2δ, (2) ReJµτ 23 = (c 2 23 s2 23 )J rcosδ +s 2 23 c2 23 c2 3 (s2 2 s2 3 c2 2 ), (22) ReJµτ 3 = (c 2 23 s2 23 )J rcosδ +s 2 23 c2 23 c2 3 (c2 2 s2 3 s2 2 ), (23) 3 where J r s 2 c 2 s 23 c 23 s 3 c 2 3. We also obtain the Jarlskog factor, J = J r sinδ. (24) It should be noted that the replacements s 23 c 23, c 23 s 23 in ReJeµ ij leads to ReJij eτ. In summary, the CP dependences of transition probabilities are obtained by substituting (5)-(24) into () as P(ν e ν µ ) = A eµ cosδ +Bsinδ +C eµ, (25) P(ν e ν τ ) = A eµ cosδ Bsinδ +C eτ, (26) P(ν µ ν τ ) = A µτ cosδ +Bsinδ +C µτ +Dcos2δ.(27) P(ν e ν µ ) and P(ν e ν τ ) have only linear terms of cosδ and sinδ. On the other hand, P(ν µ ν τ ) has also a term proportional to cos2δ in addition to the terms of cosδ and sinδ. Next, let us consider survival probabilities. One can calculate the CP dependences from the unitarity and the transition probabilities. Those are given by P(ν e ν e ) = C ee, (28) P(ν µ ν µ ) = A µµ cosδ +C µµ Dcos2δ, (29) P(ν τ ν τ ) = A ττ cosδ +C ττ Dcos2δ, (30) where the coefficients can be calculated from (5)-(24). Note that the term proportional to sin δ disappears in survival probabilities. In particular, not only sin δ term but also cosδ term does not exist in P(ν e ν e ) [28]. In P(ν µ ν µ ) and P(ν τ ν τ ), the linear term of cos2δ has opposite sign to that in P(ν µ ν τ ). The probabilities for CP conjugate and T conjugate channels are obtained from (25)-(30) by the exchange δ δ. III. DERIVATION OF EXACT FORMULA FOR P(ν α ν β) In this section, we present three identities on the product of the MNS matrix elements. Then, we derive an exact formula for the oscillation probability. Furthermore, we confirm that the Naumov-Harrison-Scott identity can be also derived from our identities. The Hamiltonian in matter H is given by H = H + a 0, (3) 2E 0 where H is the Hamiltonian in vacuum and is diagonalized as (9), a 2 2G F N e E is the matter potential, G F is the Fermi constant and N e is the electron density in matter. Note that the second term in (3) is energyindependent because the matter potential a is proportional to the energy E. The Hamiltonian in matter is diagonalized by the effective MNS matrix as Ũ HŨ = λ λ2, (32) 2E λ3 where λ i is mass eigenvalue in matter. The amplitude and probability for ν α ν β transition in matter are obtained by the replacements m 2 i λ i and U Ũ in (0)-(3). As one of the approaches to derive the probability, one can calculate single Ũαi by directly diagonalizing H and ij J αβ calculate, which is the product of four Ũs. However, the expressions for the probabilities obtained in this approach are complicated and it is not easy to extract the information on matter effects and δ. Here, we derive the probability using another approach. Note that ŨαiŨ βi appears in the amplitude in (0) and the probability in (). So, we only have to calculate the product ŨαiŨ βi and the expression for single Ũαi is not necessary. We use the relation (3) in order to calculate ŨαiŨ βi. We introduce p αβ and q αβ as p αβ = 2E H αβ, (33) q αβ = (2E) 2 Hαβ = (2E) 2 ( H γβ Hαγ H αβ Hγγ ),(34) where (αβγ) = (eµτ),(µτe),(τeµ). Then, we obtain three kinds of identities on ŨαiŨ βi. First, from unitarity

4 4 we obtain Ũ α Ũ β +Ũα2Ũ β2 +Ũα3Ũ β3 = δ αβ. (35) Second, 2E H αβ = p αβ in (33), which is also seen in Ref. [29], is rewritten as λ Ũ α Ũ β +λ 2 Ũ α2 Ũ β2 +λ 3 Ũ α3 Ũ β3 = p αβ. (36) Third, (2E) 2 Hαβ = q αβ in (34), is rewritten as λ 2 λ 3 Ũ α Ũ β +λ 3 λ Ũ α2 Ũ β2 +λ λ 2 Ũ α3 Ũ β3 = q αβ, (37) where we use the relation (2E) 2 H = (2E) 2 H (det ( H) = Ũdiag,, λ λ 2 λ 3 ) Ũ λ λ 2 λ 3. (38) Eqs. (35)-(37) on ŨαiŨ βi can be simply solved as Ũ αi Ũ βi = p αβλ i + q αβ δ αβ λ i (λ j +λ k ) ji ki, (39) where (ijk) takes (23),(23),(32). Note that p αβ and q αβ become constants in the following cases, p αβ = p αβ (α e or β e), (40) q αβ = q αβ (α = e or β = e) (4) from (3). In other cases, Hee is included in p αβ and q αβ, and they depend on a explicitly. Similar expression for ŨαiŨ βi in (39) has been obtained by Ohlsson and Snellman in Ref. [6] and also by Harrison and Scott [30] although the method of derivation is different. We have given an exact and simple formula of P(ν e ν µ ) for the first time in Ref. [23]. Here, we give the general formula of oscillation probability P(ν α ν β ) applicable to all channels by calculating Re J ij αβ and J ij from J αβ = Ũ αi Ũβi (ŨαjŨ βj ). In the case of α β, the expression is given by P(ν α ν β ) = 4 Re J ij αβ sin2 ij ±2 J sin2 ij, where Re J ij αβ = p αβ 2 λ i λ j + q αβ 2 +Re( p αβ q αβ )(λ i +λ j ) ij , (42) (43) J = Im( p eµ q eµ) , ij ij L 4E. (44) The concrete expression for λ i [35] is calculated by the MNS matrix elements and the masses in vacuum as λ = 3 s [ s2 3t u+ 3 ] 3( u 2 ), (45) λ 2 = 3 s [ s2 3t u 3 ] 3( u 2 ), (46) λ 3 = 3 s+ 2 3 u s 2 3t, (47) where s,t,u are given by s = a, (48) t = 2 3 +a[ 2 ( s 2 2c 2 3)+ 3 ( s 2 3)], (49) [ ( 2s 3 9st+27a 2 3 c 2 )] u = cos 3 cos 2 c2 3,(50) 2(s 2 3t) 3/2 in the standard parametrization [3, 4, 5]. One can see that the effective masses in matter do not depend on the CP phase δ. Finally, we reproduce the Naumov-Harrison-Scott identity from (44) as J = (2E) 3Im[ H eµ ( H eτ Hτµ H eµ Hττ ) ] = (2E) 3Im[ H eµ Hµτ Hτe ] = (2E) 3Im[H eµh µτ H τe ] = J. (5) Thus, the Naumov-Harrison-Scott identity is easily derived from the identities which we used. IV. CP DEPENDENCES OF OSCILLATION PROBABILITIES IN MATTER In this section, we calculate the oscillation probabilities in all channels by using(42)-(44) and investigate their CP dependences. P(ν e ν µ) Let us consider the CP dependence of P(ν e ν µ ) which is the best channel for the observation of CP violation. First, we calculate p eµ and q eµ from (36) and (37). As p eµ and q eµ do not include H ee, they are equivalent to the constants p eµ and q eµ, which are represented by the quantities in vacuum. Therefore, we obtain p eµ = 2 U e2 Uµ2 + 3 U e3 Uµ3, (52) q eµ = 3 2 U e Uµ (53) from (36) and (37). Next, we divide p eµ and q eµ into the terms including and not including δ as p eµ = p a eµ e iδ +p b eµ, q eµ = q a eµ e iδ +q b eµ, (54) where p a eµ, pb eµ, qa eµ and qb eµ are real numbers given by p a eµ = ( 3 2 s 2 2)s 23 s 3 c 3, (55) p b eµ = 2 s 2 c 2 c 23 c 3, (56) q a eµ = 3 2 c 2 2 s 23s 3 c 3, (57) q b eµ = 3 2 s 2 c 2 c 23 c 3, (58)

5 5 in the standard parametrization. The oscillation probability can be written in the form P(ν e ν µ ) = Ãeµcosδ + Bsinδ + C eµ (59) from (42)-(44). Note that P(ν e ν µ ) is exactly linear in sinδ and cosδ. Next, we give the exact expressionfor Ãeµ, B and C eµ. We rewrite them in the product of the oscillation part which depends on L and (Ãr) ij, Br and ( C r ) ij as à eµ = B = C eµ = (Ãr) ij sin 2 ij, (60) B r sin2 ij, (6) ( C r ) ij sin 2 ij. (62) B is expressed in the form of the sum as (6). Under the condition x+y +z = 0, the relation sin2x+sin2y+sin2z = 4sinxsinysinz (63) holds and B from (6) is rewritten in the form of product as B = B r sin2 ij = 4 B r sin 2sin 23sin 3.(64) à eµ is rewritten in the same way. Under the same condition as in deriving B, the relation sin 2 x = (sinxsinycosz +sinxcosysinz) (65) holds and Ãeµ is rewritten as à eµ = (Ãr) ij sin 2 ij = [(Ãr) jk +(Ãr) ki ]cos ij sin jk sin ki.(66) (ijk) Substituting (55)-(58) into p eµ and q eµ in (42)-(44), à eµ, B and C eµ are rewritten with the masses and mixings as à eµ = (ijk) 8[J r 2 3 λ k (λ k 3 )+(Ãeµ) k ] 2 jk 2 ki cos ij sin jksin ki, (67) B = 8J r sin 2sin 23sin 3, (68) C eµ = 4[ 2 3 s2 3 s2 23 c2 3 λ iλ j +( C eµ ) ij ] ij sin 2 ij. (69) See Appendix for the expression of (Ãeµ) k and ( C eµ ) ij. Note that these expressions (67)-(69) are still exact. In the limit of small 2, (Ãeµ) k and ( C eµ ) ij are higher order in 2 and can be ignored. Finally, we obtain the well known approximate formula by neglecting the smallest effective mass. Considering the energy of neutrino (E 0GeV) in neutrino factory experiment and the earth matter density (ρ 2.8g/cm 3 ), matter potential is given by a = 2 2G F N e E = ev 2 ρ E gcm 3 GeV ev 2. (70) It means that the smallest effective mass λ is almost equivalent to 2. Other effective masses λ 2 and λ 3, correspond to a or 3. Accordingly, the coefficients are approximated by à eµ 8J r 2 3 a( 3 a) cos 3sina sin( 3 a), (7) B 8J r 2 3 a( 3 a) sin 3sina sin( 3 a), (72) C eµ s2 23 s2 3 c2 3 ( 3 a) 2 sin 2 ( 3 a). (73) Although the approximate formula derived here is in agreementwiththeonesseeninref. [7, 8], thederivation is rather simple. The coefficients Ãeµ and B responsible for δ are O( 2 s 3 ). On the other hand, Ceµ is O(s 2 3 ). Therefore, the terms related to δ become relatively large in the case of small s 3.

6 6 P(ν e ν τ) In this subsection, we calculate P(ν e ν τ ) in the same way as P(ν e ν µ ). We need to learn ŨeiŨ τi for this purpose. At first, we consider p eτ and q eτ. Since p eτ and q eτ do not include H ee, they are equivalent to the constants p eτ and q eτ. Therefore, we obtain p eτ = 2 U e2 U τ2 + 3U e3 U τ3, (74) q eτ = 3 2 U e U τ (75) from (36) and (37). Next, we divide p eτ and q eτ into the terms including and not including δ as p eτ = p a eτ e iδ +p b eτ, q eτ = q a eτ e iδ +q b eτ, (76) where p a eτ, p b eτ, q a eτ and q b eτ are real numbers given by p a eτ = ( 3 2 s 2 2 )c 23s 3 c 3, (77) p b eτ = 2 s 2 c 2 s 23 c 3, (78) qeτ a = 3 2 c 2 2 c 23s 3 c 3, (79) qeτ b = 3 2 s 2 c 2 s 23 c 3, (80) in the standard parametrization. Note that p eτ and q eτ are obtained by the replacements s 23 c 23,c 23 s 23 in p eµ and q eµ. The coefficients ofcosδ and sinδ havethe opposite sign by the replacements, from (67) and (68). Therefore, P(ν e ν τ ) can be also written in the form where C eτ = P(ν e ν τ ) = Ãeµcosδ Bsinδ + C eτ, (8) 4[ 2 3 s2 3 c2 23 c2 3 λ iλ j +( C eτ ) ij ] ij sin 2 ij,(82) and the concrete expression for ( C eτ ) ij is given in Appendix. As for Ãeµ and B, see (67) and (68). Note that this expression is also exact. In the limit of small 2, ( C eτ ) ij is high order in 2 and can be neglected. As in the case of P(ν e ν µ ), we obtain the approximate formula as C eτ 4 2 3c 2 23s 2 3c 2 3 ( 3 a) 2 sin 2 ( 3 a), (83) by neglecting the smallest effective mass λ 2. As the result, the leading term C eτ which is not related to δ is O(s 2 3 ) as in the case of P(ν e ν µ ). It means that CP effect becomes relatively large also in this channel. P(ν µ ν τ) In this subsection, we calculate P(ν µ ν τ ). The expression for the probability of this channel is useful to analyze the appearance experiments of MINOS [3] and CNGS [32]. It is also needed for the analysis of neutrino factory experiments in which ν τ is produced by the high energy ν µ beam. In order to derive P(ν µ ν τ ), we calculate ŨµiŨ τi. Let us first consider p µτ and q µτ in Ũ µi Ũ τi. Since q µτ includes H ee, it is not constant and depends on the matter potential a explicitly as (2E) 2 q µτ = H { ( µτ = H eτ H µe H ee + a } )H µτ 2E = (2E) 2(q µτ ap µτ ). (84) Then, ŨµiŨ τi is given by Ũ µi Ũ τi = p µτλ i + q µτ ji ki = p µτ(λ i a)+q µτ ji ki. (85) Here, the different point between P(ν e ν µ ) and P(ν e ν τ ) is that the coefficient of p µτ is not λ i but λ i a. Fortunately, the matter potential a does not depend on δ, so we have only to investigate p µτ and q µτ, in order to learn the CP dependence of the probability. In the standard parametrization, we obtain p µτ = p a µτ e iδ +p b µτ +pc µτ eiδ, q µτ = q a µτ e iδ +q b µτ +qc µτ eiδ, (86) where p a µτ = 2 s 2 c 2 c 2 23s 3, p b µτ = [ 3 c (c 2 2 s 2 2s 2 3)]s 23 c 23, p c µτ = 2 s 2 c 2 s 2 23s 3, (87) qµτ a = 3 2 s 2 c 2 c 2 23 s 3, qµτ b = 3 2 ( s 2 2 +c2 2 s2 3 )s 23c 23, qµτ c = 3 2 s 2 c 2 s 2 23 s 3. (88) And then, the transition probability can be written in the form P µτ = õτ cosδ + Bsinδ + C µτ + Dcos2δ, (89)

7 7 where we obtain the expression for the coefficients after some straightforward calculations as à µτ = C µτ = D = (ijk) (ijk) 8[J r 2 3 (λ k a)(λ k a 3 )(c 2 23 s2 23 )+(õτ) k ] 2 jk 2 ki cos ij sin jksin ki, (90) 4[ 2 3 s2 23 c2 23 c4 3 (λ i a)(λ j a)+( C µτ ) ij ] ij sin 2 ij, (9) (λ k a 3 ) 2 s 2 2 c2 2 s2 23 c2 23 s2 3 2 jk 2 ki cos ij sin jksin ki, (92) and the concrete expressions for (õτ) k and ( C µτ ) ij are given in Appendix. As for B, see (68). Note that these expressions (90)-(92) are also exact. In the limit of small 2, (õτ) k and ( C eτ ) ij are high order in 2 and can be neglected. Here, let us comment on the magnitude of each term. D is O( 2 2 s2 3 ) and hence it is difficult to observe the term proportional to cos2δ. à µτ is suppressed by 2 s 3 as in the previous cases. In addition, à µτ is also proportional to c 2 23 s2 23. Since the 2-3 mixing is almost maximal from the atmospheric neutrino experiments, à µτ is considered to be rather small. On the other hand, Cµτ has not any suppression factor. Accordingly, the leading contribution comes from C µτ and the next to leading comes from B. Other contributions from õτ and D are extremely small and we give the approximate formula for C µτ. As in the previous cases, we obtain C µτ 4s 2 23 c2 23 c4 3 sin2 3, (93) by neglecting the smallest effective mass λ 2. We conclude that it is considerably difficult to observe the CP effects in the standard scenario. P(ν e ν e), P(ν µ ν µ) and P(ν τ ν τ) In this subsection, we calculate survival probabilities using unitarity and the transition probabilities derived in the previous subsections. Let us first consider P(ν e ν e ). The expression for the probability of this channel is needed, for example in neutrino factory experiments with the high energy ν e ( ν e ) beam. In the standard scenario, the unitarity relation P(ν e ν e ) = P(ν e ν µ ) P(ν e ν τ ) (94) holds. From (67) and (68) the coefficients of cosδ and sinδ are proportional to J r which includes s 23 c 23. As P(ν e ν τ ) is obtained from P(ν e ν µ ) by the exchange s 23 c 23,c 23 s 23 (i.e. J r J r ), one can understand that the terms proportional to cosδ and sinδ disappear in (94). Therefore, P(ν e ν e ) is completely independent of δ and can be written in the form P(ν e ν e ) = C ee, (95) as first pointed out in Ref. [28]. Here, C ee = C eµ C eτ (96) = + 4[ 2 3 s2 3 c2 3 λ iλ j +( C ee ) ij ] ij sin 2 ij.(97) As for the concrete expression of ( C ee ) ij, see Appendix. Note that this expression is exact. In the limit of small 2, ( C ee ) ij are higher order in 2 and can be ignored. Furthermore, we obtain an approximate formula C ee 4 2 3s 2 3c 2 3 ( 3 a) 2 sin2 ( 3 a), (98) by neglecting the smallest effective mass λ 2. P(ν e ν e ) mainly depends on s 3. This channel may be interesting because s 3 is determined regardless of the ambiguity of δ. Next, we calculate P(ν µ ν µ ). From unitarity, the relation P(ν µ ν µ ) = P(ν µ ν e ) P(ν µ ν τ ) holds. Then, the term proportional to sin δ disappears and the probability can be written in the form P(ν µ ν µ ) = õµcosδ + C µµ Dcos2δ, (99) where the coefficients are given by à µµ = Ãeµ õτ Ãeµ, (00) C µµ = C eµ C µτ C µτ. (0) In the approximation (00) we use the fact that õτ is proportional to c 2 23 s2 23 from (90) and is very small in the case of maximal 2-3 mixing. In the approximation (0), we use the fact that C eµ is suppressed by s 2 3. As for D, it is equivalent to that in P(ν µ ν τ ). Finally, it may be useful to study the CP dependence of P(ν τ ν τ ) although this channel is very difficult to observe. From unitarity, the relation P(ν τ ν τ ) =

8 8 P(ν e ν τ ) P(ν µ ν τ ) also holds. Then, the probability can be written in the form P(ν τ ν τ ) = Ãττ cosδ + C ττ Dcos2δ. (02) The coefficients are given by à ττ = Ãeτ õτ Ãeτ, (03) C ττ = C eτ C µτ C µτ. (04) Here, D is also the same as that in P(νµ ν τ ). P(ν µ ν µ ) = õµcosδ + C µµ Dcos2δ, (09) P(ν τ ν τ ) = Ãττ cosδ + C ττ Dcos2δ. (0) Comparing these results with (25)-(30), we find that the probability for each channel in matter has the same form with respect to the CP phase δ as in vacuum. It means that the matter effects just modify the coefficients. Finally, we comment that the probabilities for T conjugate channels are obtained from (05)-(0) by the replacement δ δ and those for CP conjugate channels are obtained by δ δ and a a. Summary of CP dependences in Matter Here, we summarize the CP dependences of the oscillation probabilities for all channels in matter. They are given by P(ν e ν µ ) = Ãeµcosδ + Bsinδ + C eµ, (05) P(ν e ν τ ) = Ãeµcosδ Bsinδ + C eτ, (06) P(ν µ ν τ ) = õτ cosδ + Bsinδ + C µτ + Dcos2δ,(07) P(ν e ν e ) = C ee, (08) V. EFFECTIVE MIXING ANGLES AND EFFECTIVE CP PHASE In this section, we show that our results (05)-(0) with respect to the CP dependences are reproduced from the effective mixing angles and the effective CP phase calculated by Zaglauer and Schwarzer [4]. They have derived the effective mixing angles and CP phase from Ũαβ 2 as sin 2 θ3 = λ2 3 αλ 3 +β 3 23, () sin 2 θ2 = (λ 2 2 αλ 2 +β) 3 (λ 2 αλ +β) 32 (λ 2 2 αλ 2 +β) 3, (2) sin 2 θ23 = G2 s F 2 c GFs 23 c 23 cosδ G 2 +F 2, (3) e i δ (G 2 e iδ F 2 e iδ )s 23 c 23 +GF(c 2 23 = 23 ) (G2 s F2 c GFs 23c 23 cosδ)(g 2 c F2 s GFs 23c 23 cosδ), (4) where α, β, G and F are given by α = p µµ +p ττ, (5) β = q ee, (6) Gs 23 = p a eµ λ 3 +qeµ a, (7) Fc 23 = p b eµ λ 3 +q b eµ. (8) Let us here describe the derivation of the effective mixing angles [36]. We can easily calculate the -3 mixing from Ũe3 2 = s 2 3. We can also calculate the -2 mixing by eliminating c 3 from Ũe 2 = c 2 2 c2 3 and Ũe2 2 = s 2 2 c2 3. Furthermore, the 2-3 mixing can be calculated by eliminating c 3 from the relations, Ũµ3 2 = Ũe3Ũ µ3 2 Ũe3 2 = p eµ λ 3 +q eµ (λ 2 3 αλ 3 +β) = s2 23 c2 3, (9) Ũτ3 2 = Ũe3Ũ τ3 2 Ũe3 2 = p eτ λ 3 +q eτ (λ 2 3 αλ 3 +β) = c2 23 c 2 3. (20) Next, we consider the CP dependence. One may think that the CP dependence of (4) is very complicated due to cosδ in the denominator at a glance. However, our results (05)-(0) are very simple and the probabilities are linear in sinδ, cosδ and cos2δ. It suggests that the complicated CP dependence in e i δ cancels with that in other terms. From ()-(4), δ is only included in sin θ 23 except for e i δ. Furthermore, the numerator of sin θ 23 is partially in agreement with the denominator of e i δ. The product of s 23, c 23 and e i δ is calculated as s 23 c 23 e i δ = (G2 e iδ F 2 e iδ )s 23 c 23 +GF(c 2 23 s2 23 ) G 2 +F 2. (2)

9 9 Note that the CP phase in the denominator of e i δ completely disappears. The imaginary part of (2) is known as the Toshev identity [27] s 23 c 23 sin δ = s 23 c 23 sin δ, (22) which is independent of matter effects. Finally, we demonstrate how the CP dependences of the oscillation probabilities become simple taking the channel P(ν e ν µ ) as an example. In general, the oscillation probabilities in matter are also obtained from those in vacuum by the exchange θ ij θ ij and δ δ. Namely, from (5)-(7) and (24), P(ν e ν µ ) in matter is constituted by Re J 2 eµ = ( c2 2 s2 2 ) J r cos δ + s 2 2 c2 2 c2 3 ( s2 23 s2 3 c2 23 ), (23) Re J eµ 23 r cos δ s 2 2 s2 23 s2 3 c2 3, (24) Re J eµ 3 = J r cos δ c 2 2 s 2 23 s 2 3 c 2 3, (25) J = J r sin δ. (26) Wecanseethatsin δ (cos δ)alwaysappearstogetherwith J r in Re J ij eµ or in J. On the other hand, the CP dependence for the product s 23 c 23 e i δ become simple as shown in (2). Therefore, the CP dependence of the probability is simplified although that of the effective CP phase is complicated. Then, we obtain the same results P(ν e ν µ ) = Ãeµcosδ + Bsinδ + C eµ (27) as(59). As forother channels, weobtain the same results as those calculated in our approach. This is guaranteed by the Toshev identity. VI. NEW MATTER INVARIANT IDENTITY In the previous section, we rederived two matter invariant identities. One is the Naumov-Harisson-Scott identity (5) and second is the Toshev identity (22). We obtain a new identity s 2 c 2 s 3 c 2 3 = s 2 c 2 s 3 c 2 3 (28) related to the -2 mixing and the -3 mixing by dividing the Naumov-Harrison-Scott identity by the Toshev identity. Here, we give another proof of this new identity (28). We can independently derive the identity (28) without using other identities in this proof. At first, taking k,k 2 as the eigenvalues of sub-matrix ( ) Hµµ H µτ, (29) 2E H τµ H ττ we rewrite Ũei 2 obtained from (39) as Then, we obtain Ũei 2 = λ2 i (p µµ +p ττ )λ i +q ee ji ki = (λ i k )(λ i k 2 ) ji ki. (30) ( ) 2 Ũe 2 Ũe2 2 Ũe3 2 = (λ k i )(λ 2 k i )(λ 3 k i ) i=,2 = (m 2 k i )(m 2 2 k i )(m 2 3 k i ) = const., (3) i=,2 where we use the relations obtained from the characteristic equation of Hamiltonian λ +λ 2 +λ 3 = m 2 +m2 2 +m2 3 +a, (32) λ λ 2 +λ 2 λ 3 +λ 3 λ = m 2 m2 2 +m2 2 m2 3 +m2 3 m2 +a(k +k 2 ), (33) λ λ 2 λ 3 = m 2 m 2 2m 2 3 +ak k 2, (34) from the first line to the second line in (3). Note that k i is the eigenvalue of the sub-matrix of Hamiltonian in vacuum and then does not depend on matter effects. Therefore, (3) is a matter invariant quantity and it is rewritten as (28). VII. SUMMARY We summarize the results obtained in this paper. We consider the probabilities and the CP dependences for all channels within the context of standard three neutrino oscillations in constant matter.

10 0 (i) We have derived an exact and simple formula for the oscillation probability applicable to all channels. We have also investigated the CP dependences of various probabilities using this formula. As the results, P(ν e ν µ ) can be written in the form A eµ cosδ+bsinδ+c eµ. P(ν e ν τ ) has the same form as P(ν e ν µ ). On the other hand, the probabilityfor ν µ to ν τ transition can be written as P(ν µ ν τ ) = A µτ cosδ+bsinδ+c µτ +Dcos2δ. In survival probabilities P(ν e ν e ), P(ν µ ν µ ) and P(ν τ ν τ ), there is no sinδ term. As for P(ν e ν e ), cosδ term also disappears and the probability is completely independent of δ. We have found that the probability for each channel in matter has the same form with respect to the CP phaseas in vacuum. That is, the matter effects just modify the coefficients A, B, C and D. We have also given the exact expression for the coefficients in constant matter. (ii) We have shown that our results with respect to the CP dependences are reproduced from the effective mixing angles and CP phase derived by Zaglauer and Schwarzer. Although the effective CP phase are complicated, the CP dependence of probability is simplified for all channels. This is guaranteed by Toshev identity. Finally we have obtained a new identity related to the -2 mixing and the -3 mixing by dividing the Naumov-Harrison-Scott identity by the Toshev identity. Acknowledgments The authors would like to thank Prof. A. I. Sanda for making a number of helpful suggestions. We would like to thank Prof. H. Minakata and Prof. O. Yasuda for discussions and valuable comments. [] Super-Kamiokande Collaboration, Y. Fukuda et al., Phys. Rev. Lett. 8, 562 (998); Phys. Lett. B 433, 9 (998); Phys. Lett. B 436, 33 (998); Phys. Rev. Lett. 82, 2644 (999). [2] GALLEX Collaboration, W. Hampel et al., Phys. Lett. B 447, 27 (999); SAGE Collaboration, J. N. Abdurashitov et al., Phys. Rev. C (999); Homestake Collaboration, B. T. Cleveland et al., Astrophys. J. 496, 505 (998); Super-Kamiokande Collaboration, Y. Fukuda et al., Phys. Rev. Lett. 86, 565 (200); Super-Kamiokande Collaboration, S. Fukuda et al., hepex/ [3] SNOCollaboration, Q. R. Ahmad et al., Phys. Rev. Lett. 87, 0730 (200); nucl-ex/ ; nucl-ex/ [4] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 (962). [5] CHOOZ Collaboration, M. Apollonio et al., Phys. Lett. B 420, 397 (998); Phys. Lett. B 466, 45 (999). [6] For example J. Arafune and J. Sato, Phys. Rev. D 55, 653 (997); H. Minakata and H. Nunokawa, Phys. Lett. B 43, 369 (997); S. M. Bilenky, C. Giunti and W. Grimus, Phys. Rev. D 58, (998); A. De Rujula, M. B. Gavela andp.hernandez, Nucl.Phys.B547, 2(999); K.Dick, M. Freund, M. Lindner and A. Romanino, Nucl. Phys. B 562, 29 (999); M. Tanimoto, Phys. Lett. B 462, 5 (999); A. Donini, M. B. Gavela, P. Hernandez and S. Rigolin, Nucl. Phys. B 574, 23 (2000); A. Romanino, Nucl. Phys. B 574, 675 (2000); O. Yasuda, Acta Phys. Polon. B 30, 3089 (999); J. Sato, Nucl. Instrum. Meth. A 45, 36 (2000); M. Koike and J. Sato, Phys. Rev. D 6, (2000); Erratum ibid., D 62, (2000); H. Minakata and H. Nunokawa, Phys. Lett. B 495, 369 (2000); M. Koike, T. Ota and J. Sato, Phys. Rev. D 65, (2002); T. Miura, E. Takasugi, Y. Kuno and M. Yoshimura, Phys. Rev. D 64, (200); I. Mocioiu andr.shrock, JHEP 0, 050(200); J. PinneyandO. Yasuda, Phys. Rev. D 64, (200); P. Lipari, Phys. Rev. D 64, (200); M. Freund, P. Huber and M. Lindner, Nucl. Phys. B 65, 33 (200); J. Burguet- Castell et al., Nucl. Phys. B 608, 30 (200); V. Barger, D. Marfatia and K. Whisnant, Phys. Rev. D 65, (2002); T. Kajita, H. Minakata and H. Nunokawa, Phys. Lett. B 528, 245 (2002). [7] A. Cervera et al., Nucl. Phys. B 579, 7 (2000) [Erratumibid. B 593, 73 (2000)]. [8] M. Freund, Phys. Rev. D 64, (200). [9] J. Arafune, M. Koike and J. Sato, Phys. Rev. D 56, 3093 (997); Erratum ibid., D 60, 9905 (999). [0] H. Minakata and H. Nunokawa, Phys. Rev. D 57, 4403 (998). [] Y. Itow et al., hep-ex/00609; M. Aoki et al., hepph/004220; hep-ph/ [2] S. Geer, Phys. Rev. D 57, 6989 (998). [3] V. Barger, K.Whisnant, S.PakvasaandR. J.N. Phillips, Phys. Rev. D 22, 278 (980). [4] H. W. Zaglauer and K. H. Schwarzer, Z. Phys. C 40, 273 (988). [5] Z. Xing, Phys. Lett. B 487, 327 (2000). [6] T. Ohlsson and H. Snellman, J. Math. Phys. 4, 2768 (2000); Phys. Lett. B 474, 53 (2000). [7] V. A. Naumov, Int. J. Mod. Phys. D, 379 (992). [8] P. F. Harrison and W. G. Scott, Phys. Lett. B 476, 349 (2000). [9] C. Jarlskog, Phys. Rev. Lett. 55, 039 (985). [20] H. Yokomakura, K. Kimura and A. Takamura, Phys. Lett. B 496, 75 (2000). [2] S. J. Parke and T. J. Weiler, Phys. Lett. B 50, 06 (200). [22] E. Kh. Akhmedov, P. Huber, M. Lindner and T. Ohlsson, Nucl. Phys. B 608, 394 (200). [23] K. Kimura, A. Takamura and H. Yokomakura, Phys. Lett. B 537, 86 (2002). [24] H. Minakata and H. Nunokawa, JHEP 00, 00 (200); H. Minakata, H. Nunokawa and S. Parke, Phys. Lett. B 537, 249 (2002). [25] Particle Data Group, D. E. Groom et al., Eur. Phys. J. C 5, (2000).

11 [26] P. I. Krastev and S. T. Petcov, Phys. Lett. B 205, 84 (988). [27] S. Toshev, Mod. Phys. Lett. A 6, 455 (99). [28] T. K. Kuo and J. Pantaleone, Phys. Rev. D 35, 3432 (987). [29] Z. Z. Xing, Phys. Rev. D 64, (200). [30] P. F. Harrison and W. G. Scott, Phys. Lett. B (2002). [3] Numi MINOS project at Fermi National Accelerator Laboratory, [32] CERN Neutrinos to Gran Sasso, [33] It has been first pointed out by Krastev and Petcov [26] that CP odd term is proportional to sinδ. It is easy to understand from the Naumov-Harrison-Scott identity [7, 8]. [34] Overall phase does not affect the probability from (). By using this, the eigenvalues of the Hamiltonian are written with only mass differences subtracting m 2. [35] When the calculation from amplitude to probability in matter, overall phase does not affect the probability. We remove the matter effects from NC interaction in (3) and m 2 in (32). Therefore, the expression for λ i does not depend on m 2 but mass differences. [36] Zaglauer and Schwarzer derive effective mixing angles by using Ũαβ 2. In this paper, the expression for ŨαiŨ βi is used to avoid the complexity of calculation. APPENDIX A: SUBLEADING TERMS IN THE à αβ, C αβ Here, we describe the terms which include higher order of 2 explicitly in Ãαβ and C αβ. First, the terms in P(ν e ν µ ) are (Ãeµ) k = 2 2J r [ 3 λ k (c 2 2 s 2 2)+λ 2 ks c 2 2], (A) ( C eµ ) ij = 2 s 2 3 [ 3 { λ i (λ j s c 2 2 ) λ j(λ i s c 2 ] 2 )}s2 23 c [ (λ i 3 )(λ j 3 )s 2 ] 2 c2 2 c2 23 c s2 3 [ (λ i s c 2 2 )(λ js c 2 2 )s2 23 3] c2, (A2) where the second line of ( C eµ ) ij is relatively large in the case of 2 / 3 < s 3 so that there is no suppression factor s 2 3. Second, the terms related to P(ν e ν τ ) are ( C eτ ) ij = 2 s 2 3 [ 3 { λ i (λ j s c 2 2) λ j (λ i s c 2 2)}c 2 23c 2 ] [ (λ i 3 )(λ j 3 )s 2 2c 2 2s 2 23c 2 ] s 2 3 [ (λ i s c 2 2)(λ j s c 2 2)c 2 23c 2 3], (A3) where there is also no suppression factor s 2 3 in the second line. Third, the terms related to P(ν µ ν τ ) are (õτ) k = 2 2 s 3(c 2 23 s2 23 ) [ s 2 c 2 s 23 c 23 (λ k a 3 ){ 3 (s 2 2 s2 3 c2 2 )+(λ k a 3 )c 2 2 }] (A4) ( C µτ ) ij = 2 [ 3 { 3 (λ i +λ j 2a)(c 2 3 +c 2 2) 2(λ i a)(λ j a)c 2 2}s 2 23c 2 23c 2 ] s 2 3 [ 3 { 3 (λ i +λ j 2a)+2(λ i a)(λ j a)}s 2 2s 2 23c 2 23c 2 ] [ { 3 c 2 3 +(λ i a 3 )c 2 2}{ 3 c 2 3 +(λ j a 3 )c 2 2}s 2 23c 2 ] s 2 3 [ 3 (λ i +λ j 2a 2 3 )s 2 2s 2 23c 2 23c 2 3 +(λ i a 3 )(λ j a 3 ){s 2 2c 2 2(c 2 23 s 2 23)+s 4 2s 2 23c 2 23s 2 3} ]. (A5) Finally, the terms related to P(ν e ν e ) are ( C ee ) ij = 2 s 2 3 [ 3 { λ i (λ j s c 2 2) λ j (λ i s c 2 2)}c 2 ] [ (λ i 3 )(λ j 3 )s 2 2c 2 2c 2 ] s 2 3 [ (λ i s c 2 2)(λ j s c 2 2)c 2 3], (A6) where there is no suppression factor s 2 3 in the second line as in P(ν e ν µ ) or P(ν e ν τ ). Therefore, the

12 contribution is relatively large in the case of 2 / 3 < s 3. 2