OSU-HEP-05-10 July 005 Model of Geometric Neutrino Mixing arxiv:hep-ph/050717v 0 Jul 005 K.S. Babu 1 and Xiao-Gang He, 1 Oklahoma Center for High Energy Physics Department of Physics, Oklahoma State University Stillwater, OK 74078, USA Department of Physics, Nankai University, Tianjin, China NCTS/PTE, Department of Physics, National Taiwan University, Taipei Abstract Current neutrino oscillation data from solar, atmospheric, and reactor experiments are consistent with the neutrino mixing matrix elements taking values sin θ 1 = 1/, sin θ = 1/, and sin θ 1 = 0. We present a class of renormalizable gauge models which realize such a geometric mixing pattern naturally. These models, which are based on the non Abelian discrete symmetry A 4, place significant restrictions on the neutrino mass spectrum, which we analyze. It is shown that baryogenesis via leptogenesis occurs quite naturally, with a single phase (determined from neutrino oscillation data) appearing in leptonic asymmetry and in neutrinoless double beta decay. Such predicted correlations would provide further tests of this class of models.
Introduction Our understanding of the fundamental properties of neutrinos has improved dramatically over the last few years. Atmospheric and solar neutrino experiments have by now firmly established occurrences of neutrino flavor oscillations [1]. In order for neutrinos to oscillate, they must have non degenerate masses. In addition, different neutrino flavor states must mix with one another. When positive evidence for oscillations from solar and atmospheric neutrinos are combined with results from reactor data [, ], one obtains the following neutrino mass and mixing pattern (with σ error bars) []: m = m m 1 = 7.9 10 5 (1 ± 0.09) ev, (1) m atm = m m = ±.4 10 (1 +0.1 0.61 ) ev, () sin θ 1 = 0.14(1 +0.18 0.15 ), sin θ = 0.44(1 +0.41 0. ), sin θ 1 = 0.9 +0. 0.9 10. () Here m i are the (positive) neutrino mass eigenvalues, and θ ij are the neutrino mixing angles. m m 1 > 0 in Eq. (1) is necessary for MSW resonance to occur inside the Sun. The sign of m atm, which is physical, is currently unknown. A remarkable feature of the oscillation data is that they are all consistent with a geometric neutrino mixing pattern defined by the parameters sin θ 1 = 1/, sin θ = 1/, and sin θ 1 = 0. In fact, these geometric mixing angles are very close to the central values of Eq. (). We observe that unlike the quark mixing angles, which are related to the quark mass ratios in many models (eg: θ C m d /m s ), the neutrino mixing angles seem to be unrelated to the neutrino mass ratios. The purpose of this Letter is to provide a derivation of such a geometric neutrino mixing based on renormalizable gauge theories. The neutrino mixing matrix (the MNS matrix) that we will derive has the form [4, 5, 6, 7] U MNS = 1 0 1 6 1 1 1 6 1 1 P. (4) Here P is a diagonal phase matrix which is irrelevant for neutrino oscillations, but relevant for neutrinoless double beta decay. Eq. (4) yields the desired values, U e = 1/, U µ = 1/, U e = 0. 1
Our derivation of Eq. (4) will be based on the non Abelian discrete symmetry A 4, the symmetry group of a regular tetrahedron. This symmetry group has found application in obtaining maximal atmospheric neutrino mixing [8] and in realizing quasi degenerate neutrino mass spectrum [9]. No successful derivation of Eq. (4) has been achieved to our knowledge (based on A 4 or other symmetries) in a renormalizable gauge theory context. For attempts along this line see Ref. [4, 10, 11]. In Refs. [4, 6], Eq. (4) was suggested as a phenomenological ansatz. In Ref. [10], a higher dimensional set up is used to motivate Eq. (4). Ref. [11] analyzes special cases of an A 4 derived neutrino mass matrix towards obtaining the structure of Eq. (4). A large number of models in the literature have derived maximal atmospheric mixing based on non Abelian symmetries [1], but in most models the solar mixing angle is either maximal (now excluded by data) or is a free parameter. We will see that the derivation of Eq. (4) places strong restrictions on the neutrino mass pattern. We find that the out of equilibrium decay of the lightest right handed neutrino generates lepton asymmetry at the right level to explain the observed baryon asymmetry of the universe. A single phase appears in leptonic asymmetry as well as in neutrinoless double beta decay, thus providing some hope for testing high scale phenomena via low energy experiments. The Model We work in the context of low energy supersymmmetry, which is motivated by a solution to the gauge hierarchy problem as well as by the observed unification of gauge couplings. The gauge group of our model is that of the Standard Model, SU() C SU() L U(1) Y. We augment this symmetry with a non Abelian discrete symmetry A 4. This order 1 group is the symmetry group of a regular tetrahedron. A 4 has a unique feature in describing the lepton sector: It has one triplet and three inequivalent singlet representations, thus allowing for assigning the left handed lepton fields to the triplet and the right handed charged lepton fields to the three inequivalent singlets. Denoting the three singlets of A 4 as (1, 1, 1 ), with the 1 being the identity representation, we have 1 1 = 1, 1 1 = 1, and 1 1 = 1. Furthermore, = 1 + 1 + 1 + s + a. Specifically, for the product of two triplets we have (a 1, a, a ) (b 1, b, b ) = (a 1 b 1 +a b +a b ) (1); (a 1 b 1 +ω a b +ωa b ) (1 ); (a 1 b 1 +ωa b +ω a b ) (1 );
(a b + a b, a b 1 + a 1 b, a 1 b + a b 1 ) ( s ). Here ω = e iπ/. In addition to the A 4 symmetry, we assume a Z 4 Z discrete symmetry. The Z 4 is an R symmetry under which the superpotential carries units of charge. The Z 4 and Z symmetries are broken softly in the superpotential via the lowest dimensional operators. The lepton and Higgs fields transform under A 4 Z 4 Z as follows. L : (, 1, 0), e c : (1 + 1 + 1,, 0), ν c : (, 0, 1), E : (, 1, 0), E c : (, 1, 0), H u : (1, 1, ), H d : (1, 0, 0), χ : (,, 0), χ : (,, 1), S 1, : (1,, 1). (5) Here in the fermion sector we have introduced new vector like iso singlet fields E and E c transforming under the SM gauge group as (1,1,-1) and (1,1,1), respectively, which will acquire large masses and decouple. H u and H d are the usual Higgs fields of MSSM, while χ, χ, S 1, are all SM singlet fields needed for achieving symmetry breaking. The quark fields (Q, u c, d c ) are all singlets of A 4 with Z 4 Z charges of Q(1, 1); u c (0, 0) and d c (1, ), so that the usual quark Yukawa couplings Qd c H d + Qu c H u are allowed in the superpotential. 1 The superpotential terms relevant for lepton masses consistent with the symmetries is W Yuk = M E E i E c i + f e L i E c i H d + h e ijke i e c jχ k + 1 f S l ν c i ν c i S l + 1 f ijkν c i ν c jχ k + f ν L i ν c i H u.(6) Here the flavor structure of the three independent h e ijk couplings and the one independent f ijk coupling can be easily obtained from the A 4 multiplication rules given earlier. The Higgs superpotential of the model is given by W Higs = λ χ χ 1 χ χ + λ χ s(χ 1 + χ + χ )S 1 + λ χ χ 1 χ χ + λ s11 S 1 + λ s1 S 1S + λ s1 S 1 S + λ s S + µ 1 S 1 + µ S + µ χ (χ 1 + χ + χ ). (7) Here χ = (χ 1, χ, χ ), and χ = (χ 1, χ, χ ). The last three terms in the last line of Eq. (7) break the Z 4 and Z symmetries softly. The µ 1, terms are the lowest dimensional terms that break the Z symmetry softly, while leaving Z 4 unbroken. The µ χ term is the lowest dimensional term that breaks the Z 4 symmetry softly. Such soft breaking can be 1 A bare mass term µh u H d is not allowed in the superpotential by the symmetries, but the Kahler potential, which is assumed to not respect these symmetries, allows a Planck mass suppressed term, L Hu H d Z d 4 θ/m Pl, generating the required µ term.
understood as spontaneous breaking occurring at a higher scale. We have chosen without loss of generality the combination of S 1 and S that couples to χ as simply S 1 in Eq. (7). Minimizing the potential derived from Eq. (7) in the supersymmetric limit, we obtain the following vacuum structure: S = v s, S 1 = 0; χ 1 = χ = χ = v χ ; χ = v χ, χ 1 = 0; χ = 0. (8) with v χ = µ χ /λ χ, v s = µ /(λ s ), and v χ = (λ s 1 µ λ s µ 1 )/(λ s λ χ s). Electroweak symmetry breaking is achieved in the usual way by H u = v u, H d = v d. We emphasize that vanishing of certain VEVs is a stable result, owing to the discrete symmetries present in the model. This is important for deriving the MNS matrix of Eq. (4). We observe that there are no pseudo-goldstone modes, as can be seen by directly computing the masses of the Higsinos from Eq. (7). The mass matrices M ee for the charged leptons and M νν c for the neutral leptons resulting from Eqs. (6) and (8) are given by (in the notation L = (e, E) M ee (e c, E c ) T ) 0 0 0 f e v d 0 0 0 0 0 0 f e v d 0 0 0 0 0 0 f M ee = e v d h e 1v χ h e v χ h e, v χ M E 0 0 h e 1 v χ h e ωv χ h e ω v χ 0 M E 0 h e 1v χ h e ω v χ h e ωv χ 0 0 M E 0 0 0 f ν v u 0 0 0 0 0 0 f ν v u 0 0 0 0 0 0 f M νν c = ν v u. (9) f ν v u 0 0 f s v s 0 f χ v χ 0 f ν v u 0 0 f s v s 0 0 0 f ν v u f χ v χ 0 f s v s Since the E and the E c fields acquire large masses, of order the GUT scale, they can be readily integrated out. The reduced mass matrices for the light charged leptons is 4
given by M e = U L m e 0 0 0 m µ 0 0 0 m τ, U L = 1 1 1 1 1 ω ω 1 ω ω, (10) where m i = (f e v d v χ /M E )h e i(1 + (h i v χ ) )/M E) 1/. The light neutrino mass matrix is found to be 1 0 x Mν light = m 0 0 1 x 0, (11) x 0 1 where m 0 = f νv uf s v s /(f s v s f χ v χ ), and x = f χ v χ /(f s v s ). We define x = x e iψ. M light ν can be diagonalized by the transformation Mν light = Uν D νu ν with 1 0 1 U ν = 1 0 1 + x 0 P ; D ν = m 0 1 x 1 0 1 1 x P is a diagonal phase matrix given by. (1) P = diag{e iφ 1/, e i(φ 1+φ )/, e i(φ +π)/ }, φ 1 = arg(1 + x), φ = arg(1 x). (1) These Majorana phases will not be relevant for neutrino oscillations, but they will appear in neutrinoless double beta deacy and in leptogenesis. The MNS matricx is given by U MNS = U T L U ν which has the form given in Eq. (4). (In making this identification, we make a field redefinition of the left handed charged lepton fields, e L = diag.{1, ω, ω }e L.) This is the main result of this paper. Constraints on neutrino masses From Eq. (1), the expressions for the mass eigenvalues can be inverted to obtain the following relations for the parameters m 0, x and ψ = arg(x): m 0 = m 1m x = m 1 m1 m [ m 1 m + m m m 1 m ] 1/, cosψ = (m m 1 )m m 1 m [m 1m + m m m 1m 1/. (14) ] 5
Here m 1 = m 0 (1 + x), m = m 0 (1 x ), m = m 0 (1 x). There are restrictions arising from the conditions that x be real and cosψ 1, which we analyze now. Because of the observed hierarchy m atm m, and the requirement of MSW resonance for solar neutrinos, two possible neutrino mass ordering are allowed. (i) m 1 < m < m (normal mass ordering) and (ii) m < m 1 < m (inverted mass ordering). If the neutrino masses are strongly hierarchical, m 1 m m, then from Eq. (14) one sees that cosψ 1 cannot be satisfied, since cosψ m /( m 1 ) 1 in this case. Similarly m m 1 m is also not allowed. We find that only two possibilities can arise, depending on (a) m m 1 m + m 1 (normal ordering), and (b) m m 1 m + m 1 (inverted ordering). We consider these cases in turn. These conclusions can also be arrived at by analyzing the neutrino mass matrix in the flavor basis, i.e., in a basis where the charged lepton mass matrix is diagonal: Mν flavor = m + x x x x x x 0 x x x x x x. (15) x x x x x x When x = 1 + q with q 1, we see that entries in the first row and column of Eq. (15) become small. This will be the case of normal ordering of masses. (a) Normal ordering This case is realized when m m 1 m + m 1. The condition cosψ 1 can be satisfied only if m 1 /m 1/. Making expansions in small m 1 /m and small (m 1 /m 1), we find ( ) x 1 + m 1 m1 1, m m [ { m ( ) }] cos ψ 1 1 m1 1. (16) m We see a further restriction that m 1 /m m 1 /m 1. m Neutrinoless double beta decay is sensitive to the effective mass m ββ = U ei m i = m 0 x 1 + x. (17) i In the case under study this takes the value m ββ 4 m 1 4 ( ) m 1/ 0.0068 ev. (18) 6
The effective neutrino mass measurable in tritium beta decay m νe is given by which in this case takes the value m νe [ ] 1/ m νe = U ei m i (19) i [ m ] 1/ 0.007 ev. (0) Here we made use of the fact that m 1 m. Finally, the sum of neutrino masses, sensitive to cosmological measurements, is given by m i m 1 + m i ( m ) 1/ + m atm 1/ 0.064 ev. (1) Cos Ψ mass ev -0.97-0.975-0.98-0.985-0.99-0.995 0.1 0.05 0.0 0.01 0.005 a 0.8 0.9 1 1.1 1. c 0.8 0.9 1 1.1 Ue Ε mi m j 6 5 4 1 0 b 0.8 0.9 1 1.1 1. d 0.07 0.065 0.06 0.055 0.05 0.045 0.04 0.05 0.8 0.9 1 1.1 1. Figure 1: Various quantities as functions of x for normal mass hierarchy case. (a) cos ψ vs. x ; (b) m 1 /m (dashed) and m /m (solid) vs. x ; (c) m ee (solid (green)), m νe (dotted (red)), m (dashed (black)) and m i (dot-dashed (blue)) (in ev unit) vs. x ; (d) RG running correction to U e /ǫ vs. x. 7
We have plotted in Fig. 1a - 1c various masses and mass ratios as functions of x using the exact expressions with central values for m and m atm. These results confirm our analytical solutions. (b) Inverted mass ordering In this case, m m 1 m + m 1. The small parameter expansion is different from (a), since the denominator of Eq. (16) in cosψ becomes small. Here m 1 and m are nearly equal, and m will turn out to be much smaller than m 1. In order to satisfy m m atm it is necessary that cosψ x /, which then requires x. Writing we find cosψ = x (1 + q), q 1, () m = m m 1 m 0 x (1 + x )q m atm = m m m 0 x. () It becomes clear that m < m (and thus m < m 1 ), that is, this case corresponds to an inverted hierarchy. With negative q 1, the hierarchy in the two oscillation parameters can be accommodated. In this case the effective mass for double beta decay is given by m ββ m [ ] 1/ [ ] 1/ 9 + m atm 1 m atm (4) m m with m atm being negative. Here we used the fact that m 0 m. The value of m is not determined by oscillation data. If m m atm, we have three fold degeneracy of masses and m ββ m. As m (1/9) m atm, the double beta decay amplitude vanishes. Furthermore we have for this case m νe m 1 m atm + m, i m i m 1 + m m atm + m + m. (5) The three fold degenerate case is obtained from this case by setting m much larger than m atm (or equivalently, x 1). This case also coincides with the leading results of Ref. [9]. The exact results for various masses and mass ratios are plotted in Fig. a - c as functions of x. These results confirm the analytical approximations presented here. 8
Cos Ψ mass ev 1 0.8 0.6 0.4 0. a 0 0.5 1 1.5 c 0.5 1 0. 0.1 0.05 0.0 0.01 0 0.5 1 1.5 mi m j Ue Ε 1 0.8 0.6 0.4 1 0.8 0.6 0.4 0. b 0 0.5 1 1.5 d 0 0.5 1 1.5 Figure : Various quantities as functions of x for inverted mass hierarchy case. (a) cosψ vs. x ; (b) m 1 /m (dashed) and m /m (solid) vs. x ; (c) m ee (solid (blue)), m νe m (dashed (red)) and m i (dot-dashed (black)) (in ev unit) vs. x ; (d) RG running correction to U e /ǫ vs. x. Stability of U e A distinctive feature of the geometric mixing pattern is that U e = 0 at the scale of A 4 symmetry breaking, which we have taken to be of order the GUT scale. When running from this high scale to low energy scale (M EW ), the mixing matrix may change, in particular U e may not be zero any more. One should ensure that the pattern of Eq. (4) is not destabilized, which can happen if the induced U e is too large. We demonstrate this stability now. The leading flavor-dependent effect of the running from high scale to low scale is given by the one loop RGE [1] dmν e d lnt = 1 π [Me ν Y e Y e + (Y e Y e) T Mν e ] +... (6) This leads to correction, to the leading order, to the entries M 1, (1 ǫ) and M (1 ǫ) 9
with ǫ Y τ ln(m GUT/M EW )/π. One obtains to order ǫ, U e ǫx x + cosψ + i sin ψ cosψ( x + cosψ). (7) One obtains the induced U e for the normal and inverted hierarchies by inserting the corresponding expressions for cosψ and x given earlier. The results are shown in Figs 1d (normal mass ordering) and in d (inverted ordering) where we plot U e /ǫ as a function of x. We see that the induced U e is small, too small to be measured by near future experiments for the normal mass hierarchy case in the whole allowed x range. For the inverted mass hierarchy case for x larger than about 0., U e remains small. For smaller values of x, with ǫ of order one (corresponding to Y τ 1), U e can be as large as 0.1 which may be measured in the future. In this case, all three neutrinos are nearly degenerate and the cosmological mass limit on neutrinos will be nearly saturated. We conclude that the structure of the mixing matrix derived is not upset by radiative corrections. Leptogenesis Leptogenesis occurs in a simple way in this model via the decay of the right-handed neutrinos [14]. The heavy Majorana mass matrix of ν c is given in the model as (see Eq. (9)) M ν c = M R 1 0 x 0 1 0 x 0 1. (8) The Dirac neutrino Yukawa coupling matrix is proportional to an identity matrix at the scale of A 4 symmetry breaking which we take to be near the GUT scale. The ν c fields will remain light below that scale, down to the scale M R. Renormalization group effects in the momentum range M R < µ < M GUT will induce non-universal corrections to the Dirac neutrino Yukawa coupling matrix. Without such non-universality no lepton asymmetry will be induced in the decay of right handed neutrinos. The effective theory in this momentum range is the MSSM with the ν c fields. From the renormalization group equation dy ν dt = 1 16π Y ν(y l Y l) +... (9) 10
where W = ey l LH d + ν c Y ν LH u +.., we obtain at the scale M R, Y ν = Yν 0 diag(1, 1, 1 δ) with δ (Yτ /16π )ln(m GUT /M R ). Yν 0 at the GUT scale. We diagonalize M ν c by the rotation ν c = OQN, where is the value of the universal Dirac Yukawa coupling 1 0 1 O = 1 0 0, Q = diag{e iφ 1/, 1, e iφ/ } (0) 1 0 1 so that the N fields are the mass eigenstates with real and positive mass eigenvalues: M N = M R diag( 1 + x, 1, 1 x ). In the basis where the heavy ν c fields have been diagonalized, the Dirac neutrino Yukawa coupling matrix takes the form Ŷν = QO T Y ν, so that where ν Ŷ ν Ŷ ν = Y 0 1 + (1 δ) 0 e i(φ φ 1 )/ {(1 δ) 1} 0 0 e i(φ 1 φ )/ {(1 δ) 1} 0 1 + (1 δ) The CP asymmetry arising from the decay of the field N i is given by As y 1, f(y) / y. ǫ i = 1 8π 1 [ŶνŶ ν ] ii f(y) = y j ( ) M Im{[ŶνŶ ν ] ij }f j Mi ( y 1 + log1 + y ) y. (1) (). () In the normal hierarchy case, m 1 /m = M 1 /M, so the lightest N field is N 1. In this case we have ǫ 1 = Y 0 8π ± Y 0 8π ( ν δ m1 m ν δ ( m1 m ) sin(φ φ 1 ) ) [ 1 (m 1/m 1) ] 1/. (4) (m 1 /m ) To see the numerical value of ǫ 1, we note that Y 0 ν can be of order one, δ (0.1Y τ ), and m 1 /m [ m solar / m atm ]1/ 0.1. For very large value of tanβ, Y τ 1, and we find ǫ 1 10 4. Even for moderate values of tanβ 0, we find that ǫ 1 10 6 is possible. The negative sign will also ensure the correct sign of baryon asymmetry. The induced lepton asymmetry is converted to baryon asymmetry through electroweak sphaleron 11
processes. The baryon asymmetry is given by Y B Y L /, where Y L = κǫ 1 /g, where g 00 is the effective number of degrees of freedom in equilibrium during leptogenesis, and κ is the efficiency factor obtained by solving the Boltzman s equations. A simple approximate formula for κ is [15] where [ ] 0.01 1.1 κ 10 (5) m 1 ev m 1 = v u M 1 [ŶνŶ ν ] 11 (6) For δ 0.1 and M 1 10 14 GeV, we obtain Y B 7 10 11, in good agreement with observations. For the case of inverted mass hierarchy, N is lighter than N 1, so we focus on ǫ. It is given by ν ǫ Y 0 4π Again we see that reasonable lepton asymmetry is generated. ( ) δ m x 1 x /4. (7) m 1 1 + x In summary, we have presented a class of renormalizable gauge models based on the non Abelian discrete symmetry A 4 which realize the geometric neutrino mixing pattern of Eq. (4) naturally. The resulting constraints on the neutrino masses have been outlined. We have also highlighted an intriguing connection between high scale leptogenesis and low energy neutrino experiments. Acknowledgments The work of KSB is supported in part by the US Department of Energy grant #DE-FG0-04ER46140 and #DE-FG0-04ER4106. The work of X-G.H is is supported in part by a grant from NSC. KSB would like to thank NCTS/TPE at the National Taiwan University for hospitality where this work was initiated. References [1] Q.R. Ahmad et al., (SNO Collaboration), Phys. Rev. Lett. 89, 01101 (00); Phys. Rev. Lett. 89, 0110 (00); S. Fukuda et al., (Super-Kamiokande Collaboration), 1
Phys. Lett. B59, 179 (00); B.T. Cleveland et al., Astrophys. J. 496, 505 (1998); R. Davis, Prog. Part. Nucl. Nucl. Phys., 1(1994); D. N. Abdurashitov et al., (SAGE Collaboration), Phys. Rev. D60, 055801 (1999); W. Hampel et al., (GALLEX Collaboration), Phys. Let. B447, 17 (1999); C. Cattadori, (GNO Collaboration), Nucl. Phys. B111 (Proc. Suppl.), 11 (00). [] M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 1 (004); S. Goswami and A. Y. Smirnov, arxiv:hep-ph/041159; S. Goswami, A. Bandyopadhyay and S. Choubey, Nucl. Phys. Proc. Suppl. 14, 11 (005); M. C. Gonzalez-Garcia, arxiv:hep-ph/041000; H. Back et al., arxiv:hep-ex/041016. [] G. Fogli et al., hep-ph/050608. [4] P.F. Harrison, D. H. Perkins and W.G. Scott, Phys. Lett. B458, 79 (1999); Phys. Lett. B50, 167 (00). [5] Z.-Z. Xing, Phys. Lett. B5, 85 (00). [6] X. G. He and A. Zee, Phys. Lett. B 560, 87 (00); Phys. Rev. D 68, 070 (00). [7] For a related but different ansatz, see L. Wolfenstein, Phys. Rev. D18, 958 (1978). [8] E. Ma and G. Rajasekaran, Phys. Rev. D64, 1101 (001). [9] K. S. Babu, E. Ma and J. W. Valle, Phys. Lett. B55, 07 (00). [10] G. Altarelli and F. Feruglio, hep-ph/0504165. [11] E. Ma, hep-ph/050509. [1] W. Grimus, A. S. Joshipura, S. Kaneko, L. Lavoura and M. Tanimoto, JHEP 0407, 078 (004); W. Grimus and L. Lavoura, arxiv:hep-ph/050415; K. S. Babu and J. Kubo, Phys. Rev. D 71, 056006 (005); G. Seidl, arxiv:hep-ph/001044; R. N. Mohapatra, JHEP 0410, 07 (004); K.S. Babu and S.M. Barr, Phys. Lett. B55, 89(00). [1] K.S. Babu, C.N. Leung and J.T. Pantaleone, Phys. Lett. B19, 191 (199); P.H. Chankowski and Z. Pluciennik, Phys. Lett. B16, 1 (199). 1
[14] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). [15] W. Buchmuller, P. Di Bari and M. Plumacher, Nucl. Phys. B 64, 67 (00); G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl. Phys. B 685, 89 (004). 14