Two-Zero Symmetric Neutrino Mass Matrices in Minimal Supersymmetric SO(10) S. Dev and Sanjeev Kumar ∗ † 8 0 0 Department of Physics, Himachal Pradesh University, Shimla 171005, 2 INDIA n a J 4 1 Abstract ] h The phenomenological neutrino mass matrix for two-zero symmetric texture has p - been obtained and used to rule out all possible two-zero symmetric neutrino mass p matrices obtained from Yukawa couplings to 10 and 126 Higgs representations within e h the framework of minimal supersymmetric SO(10). [ 1 v 1 Introduction 8 1 0 2 The origin of fermion masses and mixings alongwith the related problem of CP violation . 1 constitute a formidable challenge for elementary particle physics. Leaving apart extremely 0 small neutrino masses, even the charged fermion mass hierarchy ranges over at least five 8 0 orders of magnitudes. Since the fermion masses and the mixing angles are derived from v: the Yukawa couplings, which are free parameters within the Standard Model (SM), these i Yukawa couplings must span several orders of magnitude to accommodate the strongly hi- X erarchical pattern of fermion masses and mixings. However, the currently available data on r a fermion masses and mixing are insufficient for an unambiguous reconstruction of fermion mass matrices. To make matters worse, radiative corrections can obscure the underlying structure. Thus, the existing data cannot, without some additional assumptions, determine all the elements of the Yukawa coupling matrices for quarks and leptons. Some of these assumptions, invoked to restrict the form of fermion mass matrices include the presence of texture zeros [1], requirement of zero determinant [2] and zero trace condition [3] to name just a few. The main motivation for invoking different mass matrix ansatze is to relate fermion masses and mixing angles in a testable manner which reduces the number of free parameters in the Yukawa sector of SM. The recent evidence for non-zero neutrino masses and mixings leads to a further proliferation of free parameters in the Yukawa sector. In the absence of a significant breakthrough in the theoretical understanding of fermion flavors, the phenomenological approaches are bound to play a crucial role in interpreting new experi- mental data on quark and lepton mixing. These approaches are expected to provide useful ∗[email protected] †[email protected] 1 hints toward unaraveling the dynamics of fermion mass generation, CP violation and iden- tification of possible underlying symmetries of the fermion flavors from which viable models of fermion mass generation and flavor mixing could, hopefully, be constructed. The strong fermion mass hierarchy should be apparent in the fermion mass matrices them- selves with the contribution of smaller elements to physical masses and mixing angles ex- pectedtobenegligiblysmall. Thus, theseelements can, effectively, beneglectedandreplaced by zeros: the so-called texture zeros. However, the current neutrino oscillation data are con- sistent only with a limited number of texture schemes [1]. Specifically, the available neutrino oscillation data disallow all neutrino mass matrices with three or more texture zeros [1] in the flavor basis. The texture zeros in different positions in the neutrino mass matrix, in particular, and fermion mass matrices, in general, could be the consequence of some under- lying symmetry. Such universal textures of fermion mass matrices can be realized within the framework of Grand Unified Theories (GUTs). Though Grand Unification on its own does not shed any light on the flavor problem, the GUTs provide the optimal framework in which possible solutions to the flavor problem could be embedded. As mentioned earlier, texture zeros in different positions of the fermion mass matrices could result because of an underlying symmetry. Grand Unified models attempt to explain the masses and mixings in both the quark and lepton sectors simultaneously. The textures for the mass matrices obtained in these models can either be assumed at the very outset or can be derived from the observed mixing matrix in the flavor basis. Alternatively, the textures of the mass matrix can be obtained by embedding some family symmetry within the chosen Grand Unification group. One particularly interesting class of models is that based upon the SO(10) grand uni- fication group. There are two kinds of minimal models in this class: those based upon Higgs representation with dimension 10, 126, 126 and possibly also 120 and/or 210 and those based upon 10, 16, 16, and 45 representations. The former choice, generally, has symmetric and/or antisymmetric texture mass matrices while the latter type generally imply lopsided mass matrices forthe down type quarks andcharged leptons. In thepresent work, we restrict ourselves to the former choice i.e. we adopt the so called symmetric two-zero texture for the up-type quark mass matrices within the SUSY SO(10) GUT framework. Within this framework, we not only have a relation between down type quark mass matrices (M ) and d charged lepton mass matrices (M ) but also a relation between the up-type quark mass ma- l trices (M ) and Dirac neutrino mass matrices (M ). Thus, once we fix the representation u νD of the Higgs field corresponding to each element, M and M are uniquely determined from l νD M and M respectively. d u Neutrino mass matrices with two texture zeros in the charged lepton basis have only one degree of freedom [4, 5]. The presence of two texture zeros in the neutrino mass matrix imply four conditions on the nine free parameters of the model. The four parameters out of the remaining five parameters are determined by the neutrino data for the values of two squared-mass differences ∆m2 and ∆m2 and two mixing angles θ and θ . So, we are left 12 23 12 23 with only one free parameter in the neutrino mass matrix [4, 5]. However, we have one more experimental measurement by the CHOOZ experiment which establishes an upper bound on the third mixing angle θ . A lower bound on θ is inherent in the nature of two texture zero 13 13 neutrino mass matrices [4, 6]. So, even the remaining one degree of freedom is constrained. In fact, the neutrino mass matrix with two texture zeros can be completely determined [7] with the help of the presently available neutrino data for the seven allowed texture zero schemes of Frampton, GlashowandMarfatia [1]. Therefore, the present neutrino data, which 2 is sufficient to determine the neutrino mass matrix with two texture zeros, will be able to overrule various neutrino mass models with two texture zeros if the model parameters are known and to constrain the model parameters if they are unknown. In the class of GUT models under consideration, the mass matrices M for the Dirac νD neutrinos which gives rise to neutrino masses via see-saw mechanism can be taken identical toM except fortheaccompanying CGcoefficients which depend upontherepresentations of u thecouplingHiggsfield[8,9]. Therefore, thesemassmodelspredict theneutrinomassmatrix at GUT scale. This theoretical mass matrix should be consistent with the phenomenological mass matrix at weak scale calculated for the texture scheme consistent with the model. However, the neutrino mass matrix calculated at GUT scale has to be run down to the weak scale before comparing it with the neutrino mass matrix calculated phenomenologically from its texture scheme from the neutrino data at the weak scale. The effect of renormalization group (RG) running will be small for the neutrino mass matrices with normal hierarchy. However, it can be significantly large for the mass matrices with inverted or quasi-degenerate hierarchy. In the present work, we first derive the phenomenological neutrino mass matrix for a partic- ular texture scheme proposed by Frampton, Glashow and Marfatia [1] and confront it with a class of two-zero symmetric texture GUT models based upon minimal supersymmetric SO(10) [8, 9]. 2 Phenomenological neutrino mass matrix with a two zero symmetric texture We take the neutrino mass matrix to be the real symmetric matrix with two texture zeros: 0 A 0 M = A D B (1) ν 0 B C where C = m m +m D, A2 = m m m /C, 1 − 2 3 − 1 2 3 (2) B2 = (m +m D)(m m D)(m m +D)/C 3 1 3 2 2 1 − − − − since we are interested in obtaining only the relative magnitudes of the neutrino mass matrix elements without phases or signs. This is the texture of Frampton, Glashow and Marfatia 2 A [1]. The eigenvalues of M are m , m and m . The real orthogonal matrix O which ν 1 2 3 − diagonalizes the neutrino mass matrix M according to the relation OTMO = diag(m , m ,m ) (3) 1 2 3 − 3 is given by O11 = (m1 m2+mm2m33(Dm)3(m−m1+2m−D2))(m3 m1) r − − − OO1123 == rrOO22((12mm11==−−mmrr22++mm((mmmmmm3112mm3312++−−((12mmmm((DDmm3323))−+))32((((mm+−mmmmmm122131++−−−+11mm−+DDmmDD23))12))))))((mm33+−mm21)) (4) O = m3(m2 m1+D) where we have given only tOOOhe333231m===agrrrn2i(((t3mmmu111dmmm−−−emmm123rs(((222mmmo+++(223mfmmm−−−3t333mmm−h−−−112meDDD++−1)))eDDD−)(((l(mmm)))em(((m232mmm2++++333emmm+−+mn123mmm3t))))s(((121mmm−−−o313DDDf−+−)))Ommm121.)))This diagonalization scheme is taken from [10] where the references to earlier works can, also, be found. Many interesting relations can be extracted from the above expressions for the elements of the neutrino mass matrix and its eigenvalues. The two mass ratios m1 and m1 can be written m2 m3 as m O O 1 12 21 = (5) m O O 2 11 22 and m O O 1 13 21 = . (6) m O O 3 11 23 The elements A, B and C are given by O O 11 12 A = m m , (7) 1 2 s O O 21 22 (m +m )(m +m )(m m ) 1 2 2 3 3 1 B = − O O O (8) 12 21 23 m m 1 3 and O O 21 22 C = m . (9) 3 O O 11 12 Therefore, the ratios A and B become C C A O O O 11 12 13 = (10) C O O O 21 22 23 and B (m +m )(m +m )(m m ) 1 2 2 3 3 1 = O O O − . (11) 11 12 13 C m m m 1 2 3 Another important relation is √m m O O 1 2 11 12 O = O . (12) 13 23 m sO O 3 21 22 4 Eqs. (5)-(12) are exact and have been derived from Eqs. (4). These relations imply very interesting consequences when written approximately as Taylor series in the powers of s 13 coming from the orthogonal mixing matrix O which can be written as c c s c s 12 13 12 13 13 O = s c c s s c c s s s s c (13) 12 23 12 23 13 12 23 12 23 13 23 13 − − − s s c c s c s s c s c c 12 23 12 23 13 12 23 12 23 13 23 13 − − − where the symbols have their usual meaning. We substitute these elements in equations (5)-(12) and expand the results in the powers of s . The mass ratios in Eqs. (5) and (6) 13 are given by m s2 1 = 12 + (s ) (14) m c2 O 13 2 12 and m s c s 1 = 12 23 13 + (s2 ) (15) m c s O 13 3 12 23 which are consistent with our earlier results [4, 5]. The mass ratio m1 is much smaller than m2 one. Therefore, the neutrino mass matrix considered here has a normal hierarchy. This fact is important since the neutrino oscillation parameters change very little in the RG evolution from the weak scale to GUT scale for the normal hierarchy. The matrix elements A, B and C are given by √m m 1 2 A = + (s ), (16) 13 c O 23 B = m c s + (s ) (17) 3 23 23 13 O and C = m c2 + (s ). (18) 3 23 O 13 Moreover, the element D can be expressed as D = m s2 + (s ). (19) 3 23 O 13 The ratios A and B can be written as C C A s = 13 + (s2 ) (20) C c2 s O 13 23 23 and B = 1+ (s ). (21) 13 C O Eq. (12) for O gives 13 √m m s = 1 2 + (s2 ). (22) 13 m O 13 3 This equation is consistent with the analytical results obtained by Desai et al. [7]. Eq. (14) can be used to obtain the relation ∆m2 m = s2 12 (23) 1 12scos2θ 12 5 as shown in our earlier work [11]. Using Eq. (23) and neglecting the higher order terms in s , Eq. (22) can be written as 13 1 ∆m2 2 s c s 12 12 12 . (24) 13 ≈ ∆m213! c2 s2 12 − 12 q For the tri-bi-maximal value [12] s2 = 1, the trigonometric factor in Eq. (24) is given by 12 3 s c 2 12 12 = . (25) c2 s2 s3 12 − 12 q Substituting the present best fit values of the oscillation parameters in Eq. (24) (∆m2 = 12 7.9 10 5, ∆m2 = 2.4 10 3, and s = 0.3 [13]), we find that θ = 7.3o which is consistent × − 13 × − 12 13 with the earlier estimates [4, 5] for this mass matrix texture. Using the above leading order estimates of the elements A, B, C and D, the neutrino mass matrix given in Eq. (1) can be written as 0 ǫ 0 m 3 M = ǫ 1+ (ǫ) 1+ (ǫ) (26) 2 O O 0 1+ (ǫ) 1+ (ǫ) O O where the parameter ǫ is given by ǫ = 2√2s (27) 13 for maximal atmospheric mixing. For the present best fit values of the neutrino oscillation parameters, ǫ 0.36. So, the hierarchical structure of neutrino mass matrix should be as ∼ follows: 0 0.36 0 1 M ∆m2 0.36 1 1 (28) ∼ 2 13 q 0 1 1 3 A minimal supersymmetric SO(10) realization of the two-zero symmetric texture The phenomenological neutrino mass matrix obtained here has to be confronted with the neutrino mass matrices given by various neutrino mass models. As an example, we confront a class of minimal supersymmetric SO(10) GUT models with two-zero symmetric texture [8, 9] with the phenomenological mass matrix obtained above. In this class of models, each element of M is dominated by contribution either from 10 or 126 Higgs representation and u the Yukawa couplings of Dirac neutrinos ν are that of corresponding up quarks multiplied D by a factor of 1 or -3 respectively. Therefore, the Dirac neutrino mass matrix for this class of models can be written as 0 a 0 M = a b c m (29) νD t 0 c d 6 which is related to the up-quark mass matrix M u 0 a 0 u M = a b c m (30) u u u u t 0 c 1 u through the relations a = x a , b = x b , c = x c , d = x (31) 12 u 22 u 23 u 33 where √m m m m u c c u a = , b = , c = (32) u u u m m m t t s t and m , m and m are the masses of the up quarks u, c and t at the GUT scale. This mass u c t matrix is consistent with the quark masses and mixing data. The coefficients x , x , x 12 22 23 and x are the CG factors which connect the respective elements of M and M and are 33 u νD equal to 1 or -3 depending upon whether the corresponding Higgs representation is 10 or 126. The Higgs representation for M which has the best agreement with the data is given u by [8] 0 126 0 M : 126 10 10 . (33) u 0 10 126 For this M , the CG coefficients are given by u x = 3, x = 1, x = 1, x = 3. (34) 12 22 23 33 − − The link with the neutrino mass matrix is established through the see-saw mechanism using the right handed neutrino mass matrix matrix M : R 0 rm 0 R M = rm 0 0 . (35) R R 0 0 m R This particular form of M is consistent with the 126 dimensional representation for Higgs R field which couples with up quarks. It gives two degenerate right handed neutrino states of mass rm and a third neutrino state of mass m . Thus, r is the ratio of two right handed R R neutrino mass scales. The neutrino mass matrix M can be obtained from the right-handed ν Majorana mass matrix (M ) and the Dirac mass matrix (M ) via the see-saw relation R νD M = MT M 1M (36) ν νD R− νD and is given by 0 a2 0 r m2 M = a2 2ab +c2 c(a +1) t . (37) ν r r r m 0 c(a +1) d2 R r The above neutrino mass matrix, to a very good approximation, can be written as [8, 9] 0 a2 0 r m2 M = a2 2ab ac t . (38) ν r r r m 0 ac d2 R r 7 It will be convenient to parameterize M in the following way [8, 9] ν 0 β 0 d2m2 M = β α h t (39) ν m 0 h 1 R where 2ab x x 2a b 12 22 u u α = = , (40) rd2 x233 r ! a2 x2 a2 β = = 12 u (41) rd2 x233 r ! and ac x x a c 12 23 u u h = = . (42) rd2 x2 r 33 (cid:18) (cid:19) Since, α, β and h are all proportional to 1, we can eliminate r and work in term of the ratios r β and α which are given by α h β x 1 m 12 u = (43) α x23 2smc! and α x m 22 c = 2 . (44) h x23 √mumt! The ratio α is found to be of the order of unity and is in agreement with the value obtained h from the phenomenological neutrino mass matrix [Eq. (26) and Eq. (28)]. However, this ratio depends on the top quark mass m which has rather large uncertainty. Moreover, the t quark masses have to be run to the GUT scale and the RG running effects are significantly large for this ratio. Therefore, this ratio cannot serve as a good criterion to test this model. On the other hand, the ratio β depends on the relatively better known quark masses m and α u m and the RG running effects in calculating this ratio at the GUT scale are quite small. c In fact, this ratio is quite stable against the various input SUSY parameters (like tan2β and threshold corrections) which affect the RG evolution [14]. Substituting the value of the quark mass ratio mu at GUT scale [14]: mc m u = 0.0027 0.0006, (45) m ± c we find that the model prediction for the ratio β is about 0.078 0.009. α ± The ratio β can, also, be obtained from the phenomenological neutrino mass matrix since α β A A = = + (s ) ǫ (46) 13 α D C O ∼ and hence β 0.4. This value has to be evolved to the GUT scale by the RG running. α ∼ However, it can be seen that this ratio is quite stable against RG evolution which can lower this ratio by a factor of about 9 at the most when going from the weak scale to the GUT 10 scale. This can be seen from the relations [15] α β h α , β , h (47) → (1 ǫ )2 → (1 ǫ )(1 ǫ ) → (1 ǫ ) µ e µ µ − − − − 8 where the sign ‘ ’ denotes the RG evolution from the weak scale to the GUT scale and ǫ e → and ǫ are the RG parameters which can be at the most 0.1. Hence, the smallest possible µ value of the ratio β at the GUT scale is approximately 0.3 which is still larger thanthe model α value for this ratio by a factor of about 4. This disagreement between the phenomenological values of the ratio β calculated from the neutrino data and the values predicted by the GUT α model can not be reconciled even if the errors of the neutrino oscillation data are taken into account. We have done the exact numerical calculation and found that β = 0.36 0.05. This α ± value differs from the model value by more than 5 σ C.L.. Therefore, the two-zero symmetric texture of the neutrino mass matrix based upon the minimal supersymmetric SO(10) GUT is ruled out by more than 5 σ C.L.. The neutrino mass matrix M given in Eq. (39) corresponds to the up-quark mass matrix ν M with the Higgs representation given in Eq. (33) which is only one of the possible Higgs u representations. M , a symmetric matrix with two texture zeros, contains four independent u elements which can couple with either 10 or 126 dimensional Higgs fields. There are 16 such representations listed in Table 1 where we have adopted the classification scheme of [9] for comparison. Only eight textures out of the total sixteen were found to be allowed in earlier studies [9]. Just like the basic texture based upon the Higgs representation given in Eq. (33) for M , the other fifteen possible textures with different possible Higgs representation for M u u are also ruled out by the phenomenological values for the ratios β and α. The representation α h chosen for M in Eq. (33) is and so we can call the M obtained in Eq. (39) to be of u 1 ν S type . Similarly, one can calculate the neutrino mass matrix for all the sixteen Higgs field 1 S representations of M which will differ only in the values of the CG factors x , x , x and u 12 22 23 x . These factors for all the sixteen type of neutrino mass matrices have been tabulated 33 in Table 2 where we also calculate the parameters α, β and h in the units of 2aubu , a2u r r and aucu respectively and their ratios β and α in the units of 1 mu an(cid:16)d 2(cid:17) m(cid:16)c (cid:17) r α h 2 mc √mumt respe(cid:16)ctivel(cid:17)y. In these units, the ratio α should be about 1 and the(cid:18)ratqio β(cid:19)should(cid:16)be abou(cid:17)t h α 15 to explain the present neutrino data. It can be seen that the ratio α is 1 in agreement h with its phenomenological value only for the categories , and . These categories were S A B favored by the experimental data in some earlier studies [8, 9]. However, the ratio β is α approximately smaller by the factors of 5, 15 and 45 from its phenomenological value for the classes , and , respectively, ruling out these categories. The categories and are S A B C F ruled out in a straightforward manner since they do not reproduce either of the two ratios correctly. Hence, the entire class of two-zero symmetric texture of neutrino mass matrices is ruled out. 4 Conclusions In conclusion, we have obtained the phenomenological neutrino mass matrix for the two-zero symmetric texture and found that it is not compatible with all possible two-zero symmetric texture neutrino mass matrices obtained from the Yukawa couplings with 10 and 126 Higgs representations within the framework of minimal supersymmetric SO(10). Therefore, the theoreticalmotivationfortwo-zerosymmetric textureoftheneutrino massmatrix, nolonger, exists. 9 References [1] Paul H. Frampton, Sheldon L. Glashow and Danny Marfatia, Phys. Lett. B 536, 79 (2002). [2] G. C. Branco, R. Gonzalez Felipe, F. R. Joaquim, T. Yanagida, Phys. Lett. B 562, 265 (2003); B. C. Chauhan, J. Pulido, M. Picariello, Phys. Rev. D 73, 053003 (2006). [3] X.-G. He, A. Zee, Phys. Rev. D 68, 037302 (2003). [4] S. Dev, Sanjeev Kumar, Surender Verma and Shivani Gupta, Nucl. Phys. B 784, 103 (2007). [5] S. Dev, Sanjeev Kumar, Surender Verma and Shivani Gupta, Phys. Rev. D 76, 013002, (2007). [6] S. Dev, Sanjeev Kumar, Surender Verma and Shivani Gupta, Phys. Lett. B 656, 79 (2007). [7] Bipin R. Desai, D. P. Roy and Alexander R. Vaucher, Mod. Phys. Lett A 18, 1355 (2003). [8] M. Bando and M. Obara, hep-ph/0212242. [9] M. Bando and M. Obara, Prog. Theor. Phys., 109, 995 (2003); Masako Bando, Satoru Kaneko, Midori Obara and Morimitu Tanimoto, Phys. Lett. B 580 ,229 (2004). [10] Gulsheen Ahuja, Sanjeev Kumar, Monika Randhawa, Manmohan Gupta and S. Dev, Phys. Rev. D 76, 013006(2007). [11] S. Dev and Sanjeev Kumar, Mod. Phys. Lett. A 22, 1401 (2007). [12] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002); Phys. Lett. B 458, 79 (1999). [13] G. L. Fogli, E. Lisi, A. Marrone and A. Palazzo, Prog. Part. Nucl. Phys. 57, 742 (2006). [14] Graham Ross and Mario Serna, arXiv:0704.1248v1 [hep-ph]. [15] M. Bando and M. Obara, hep-ph/0307360; H. Haba and N.Okamura, hep-ph/9906481. 10