arxiv:hep-ph/ v3 8 Jul 2005

Transcription

1 KEK-TH-9 SO() Group Theory for the Unified Model Building arxiv:hep-ph/4v 8 Jul Takeshi Fukuyama Amon Ilakovac Tatsuru Kikuchi Stjepan Meljanac 4 and Nobuchika Okada Department of Physics Ritsumeikan University Kusatsu Shiga -877 Japan Department of Physics University of Zagreb P.O. Box Bijenička cesta HR- Zagreb Croatia Institut Rudjer Bošković Bijenička cesta 4 P.O. Box 8 HR- Zagreb Croatia Theory Group KEK Oho - Tsukuba Ibaraki -8 Japan February 8 Abstract The complete tables of Clebsch-Gordan (CG) coefficients for a wide class of SO() SUSY grand unified theories (GUTs) are given. Explicit expressions of states of all corresponding multiplets under standard model gauge group G = SU() C SU() L U() Y necessary for evaluation of the CG coefficients are presented. The SUSY SO() GUT model considered here includes most of the Higgs irreducible representations usually used in the literature: 4 4 and. Mass matrices of all G multiplets are found for the most general superpotential. These results are indispensable for the precision calculations of the gauge couplings unification and proton decay etc. fukuyama@se.ritsumei.ac.jp ailakov@rosalind.phy.hr rp9979@se.ritsumei.ac.jp 4 meljanac@irb.hr nobuchika.okada@kek.jp

2 Introduction A particularly attractive idea for the physics beyond the standard model (SM) is the possible appearance of grand unified theories (GUTs) []. The idea of GUTs bears several profound features. Perhaps the most obvious one is that GUTs have the potential to unify the diverse set of particle representations and parameters found in the SM into a single comprehensive and hopefully predictive framework. For example through the GUT symmetry one might hope to explain the quantum numbers of the fermion spectrum or even the origins of fermion mass. Moreover by unifying all U() generators within a non- Abelian theory GUTs would also provide an explanation for the quantization of electric charge. By combining GUTs with supersymmetry (SUSY) we hope to unify the attractive features of GUTs simultaneously with those of SUSY into a single theory SUSY GUTs []. The apparent gauge couplings unification of the minimal supersymmetric standard model (MSSM) is strong circumstantial evidence in favor of the emergence of a SUSY GUT near M G GeV []. While there are a priori many choices for such possible groups the list can be narrowed down by requiring groups of rank 4 that have complex representations. The smallest groups satisfying these requirements are SU() SU() SO() and E. Amongst these choices SO() is particularly attractive [4] because SO() is the smallest simple Lie group for which a single anomaly-free irreducible representation (irrep) (namely the spinor representation) can accommodate the entire SM fermion content of each generation. Once we fix SO() as the gauge group we have also many choices of the Higgs fields though they are limited by the gauge symmetry. The Higgs fields play an essential role in the spontaneous symmetry breaking of the SO() gauge group and as a source of the observed fermion masses. The SO() gauge group must be broken down to the standard model gauge group G = SU() C SU() L U() Y gauge group and each SO() irrep has to be decomposed into G multiplets. In this paper we make an explicit construction of the states of these G multiplets. Using these states we calculate the CG coefficients appearing in the mass matrices for the states belonging to G irreps and corresponding mass matrices in a wide class of the SO() models. The purpose of the present paper is to give detailed structures of the SO() GUTs based on a general model as far as has been possible and to serve a wide range of unified model builders. The paper is organized as follows. In Sec. we give a class of SUSY SO() GUTs and give an explicit form of the most general superpotential. In such a superpotential we postulate a renormalizability in order to keep the predictability [ 7]. However the result developed here is also applicable to the non-renormalizable models [8 9 4]. The symmetry breakings are considered in Sec.. Sec. 4 is devoted to present explicit forms of the states in the G multiplets for all SO() irreps. This is the central part of the present paper. Using these tables we give in Sec. the mass matrices with the CG coefficients for a class of SUSY SO() GUTs together with suitable tests and consistency checks for them. In Sec. we consider the quark and lepton mass matrices in

3 general SUSY SO() models. Sec. 7 is devoted to conclusion. We list the decompositions of each SO() irreps under G subgroup in Appendix A. In Appendix B we present the complete list of the CG coefficients for the G multiplets for all SO() irreps. A Class of SUSY SO() GUTs In this section we consider a class of renormalizable SUSY SO() models. They include three families of matter fields Ψ i (i = ) transforming as dimensional fundamental spinor representation gauge fields contained in the adjoint representation and set of SO() multiplets of Higgs fields enabling most general Yukawa couplings. The most general Yukawa couplings follow from decomposition of = + + i.e. they include the Higgs fields in H = D = = irreps respectively. Furthermore to consider as general case of the symmetry breaking of SO() to the standard model gauge group G as possible we add several Higgs fields containing G singlets. They are A = 4 = Φ = and E = 4 irreps. Of course that is not the most general case. However this set of Higgs fields is quite rich and gives rise to several realistic SUSY SO() models. Our aim is to give a systematic method for treatment of models with complicated Higgs sectors. This is a generalization of the method proposed in Refs. [] [] and [7]. We shall assume that the SUSY is preserved so that we consider the breaking of SUSY SO() to the MSSM. Then the Yukawa couplings are W Y = Y ij Ψ i HΨ j + Y ij Ψ i DΨ j + Y ij Ψ i Ψ j () where i j = denote the generation indices. Note that H is a fundamental SO() irrep and A D Φ and + are antisymmetric tensors of rank 4 and respectively. E is a symmetric traceless tensor of rank. The most general Higgs superpotential is given by: W = m Φ + m + m H + m 4A + m E + m D + λ Φ + λ Φ + ( λ + λ 4 ) HΦ + λ A Φ iλ A + λ 7 εaφ + E ( λ 8 E + λ 9 A + λ Φ + λ + λ + λ H ) + D (λ 4 E + λ Φ) + D { λ HA + λ 7 HΦ + ( λ 8 + λ 9 ) A + ( λ + λ ) Φ } ()

4 where SO() invariants are defined in the fundamental SO() basis 9 and in the Y diagonal basis 9 (which we will define in the next section) as follows: Φ Φ a b c d Φ a b c d = Φ abcdφ abcd a b c d e a b c d e = abcde abcde H H a H a = H a H a A A a b A a b = A aba ab E E a b E a b = E abe ab D D a b c D a b c = D abcd abc Φ Φ a b c d Φ a b e f Φ c d e f = Φ abcd Φ abefφ cdef Φ Φ a b c d a b e f g c d e f g = Φ abcd abefg cdefg HΦ a b c d e H a Φ b c d e = abcde H aφ bcde HΦ a b c d e H a Φ b c d e = abcde H aφ bcde A Φ A a b A c d Φ a b c d = A ab A cdφ abcd ia ia a b a c d e f b c d e f = ia ab acdef bcdef εaφ ε a A a a a 4 a a a 7 a 8 a 9 a a Φ a a Φ a 4 a a a 7 a 8 a 9 a = ε a a a a 4 a a a 7 a 8 a 9 a A a a Φ a a 4 a a Φ a7 a 8 a 9 a E E a b E a c E b c = E ab E ace bc EA E a b A a c A b c = E ab A aca bc EΦ E a b Φ a c d e Φ b c d e = E ab Φ acdeφ bcde E E a b a c d e f b c d e f = E ab acdef bcdef E E a b a c d e f b c d e f = E ab acdef bcdef EH E a b H a H b = E ab H ah b ED E a b D a c d D b c d = E ab D acdd bcd D Φ D a b c D a d e Φ b c d e = D abc D adeφ bcde DHA D a b c H a A b c = D abc H aa bc DHΦ D a b c H d Φ a b c d = D abc H dφ abcd D A D a b c a b c d e A d e = D abc abcdea de D A D a b c a b c d e A d e = D abc abcdea de D Φ D a b c a b d e f Φ c d e f = D abc abdefφ cdef D Φ D a b c a b d e f Φ c d e f = D abc abdefφ cdef. () Here a b c (a b c ) run over all the SO() vector (Y diagonal) indices and ε is 4

5 a totally antisymmetric SO() invariant tensor with ε = iε 4789 =. (4) Symmetry Breaking Here we first introduce Y diagonal basis (see also Ref. [8]): = + i = i = +4 i 4 = 4 i = + i = i 7 = 7 +8 i 8 = 7 8 i 9 = 9 + i = 9 i up to the normalization factor. It is more convenient since ( 7 9) transforms as -plet and ( 4 8 ) transforms as -plet under SU() U() X (for that reason Y diagonal basis could also be called SU() basis). Consequently ( ) and ( 4) are SU() L doublets with definite hypercharges Y = and Y = respectively. Similarly ( 7 9) and ( 8 ) transform under SU() C as and with definite hypercharges Y = and Y = respectively. Note that under the complex conjugation (c.c.) = = 4 = 7 = 8 9 = and vice versa. The SO() invariants are build in such a way that an index a is contracted (summed) with the corresponding c.c. index a for example T a T a. The basis in A = 4 D = Φ = and + = + dimensional spaces are defined by totally antisymmetric (unit) tensors (a b ) (a b c ) (a b c d ) and (a b c d e ) respectively and similarly in a b c d e indices in Y diagonal basis. The states of the and have additional properties iεā ā ā ā 4 ā ā ā 7 ā 8 ā 9 ā a a 7 a 8 a 9 a = ā ā ā ā 4 ā iεā ā ā ā 4 ā ā ā 7 ā 8 ā 9 ā a a 7 a 8 a 9 a = ā ā ā ā 4 ā () that allow one to project out the and states respectively from the antisymmetric states (abcde). The explicit expressions for antisymmetric tensors are for example etc. Important relations are (ab) = ab ba (abc) = abc + cab + bca bac acb cba () () = i( ) (4) = i( 4 ) () = i( ) (78) = i(7 8 ) (9) = i(9 ) (7) Symmetric E = 4 dimensional space is spanned by traceless symmetric states {a b } a b + b a (a b = 9 ) and a c a {a a } with a c a. Also important

6 relations are {} = + {4} = {} = + {78} = {9} = (8) Now the Higgs fields A E and Φ contain 8 directions of singlets under the G subgroup (see Appendix A). The corresponding vacuum expectation values (VEVs) are defined by: A = A i Âi (9) i= E = E Ê () = v R v R () = v R v R () Φ = Φ Φ i i () where unit directions Âi (i = ) Ê v R v R and Φ i (i = ) in the Y diagonal basis are: i= Â = Â() () = i ( + 4) (4) Â = Â() () = i ( ) () Ê = v R = v R = Φ = Φ = Ê() () = { [ + 4] [ ]} () () () = () () = (48) (7) (79) (8) Φ () () = 4 (4) (9) Φ () () = 7 ( ) () () Φ = Φ () = ([ + 4][ ]). () Here and hereafter the upper and the lower indices indicate the SU() C SU() L U() Y SU(4) SU() L SU() R quantum numbers respectively in the case of double indices. A

7 word about notation: the square brackets are used for grouping of indices. This grouping of indices is used to emphasize the SU() L and SU() C structures within the state vectors. The square brackets satisfy usual distributive law with respect to summation of indices and tensor product of indices e.g. ([ + 4][ ]) = ( ) = () + (78) + (9) + (4) + (478) + (49) ([ ][[78 + 9] 7[ + 9] 9[ + 78]]) = ( ) = ( ) + ( ). () Further the numerical factors which could be misinterpreted as additional SO() indices are written in italics. The unit directions appearing in VEVs satisfy the following orthonormality relations  i Âj = δ ij (i j = ) Ê = v R v R = v R v R = v R v R = Φ i Φ j = δ ij (i j = ). () Due to the D-flatness condition the absolute values of the VEVs v R and v R are equal v R = v R. (4) The superpotential of Eq. () calculated at the VEVs in Eqs. (4) () is: W = m Φ + m + m 4 A + m E + λ Φ + λ Φ + λ A Φ iλ A + λ 7 ε A Φ + E [ λ 8 E + λ 9 A + λ Φ ]. () Inserting the VEVs from Eqs. (4) () one obtains W = m [ Φ + Φ + ] Φ + m v R v R 7

8 + m ( ) 4 A + A + m E + λ [Φ 9 + Φ Φ [ + λ Φ + Φ + λ [A Φ + A Φ [ ( + λ A ) ( + A ] + Φ Φ 9 ] + Φ v R v R ] + A A Φ )] v R v R [ + λ 7 A Φ Φ + A Φ ] + A Φ Φ [ ( + λ 8 E + λ 9E A + A )] [ ( + λ E Φ + Φ ) ] + Φ 4. () The VEVs are determined by the following equation: { Φ Φ Φ v R v R A A = m Φ + λ Φ + λ v R v R + λ A + λ 7 A Φ } W =. (7) E From Eq. (7) we obtain seven equations for Φ Φ Φ A A E and v R v R. They are: λ Φ E = m Φ + λ Φ + λ Φ + λ v R v R + + λ 7 A Φ = m Φ + λ Φ Φ λ 7 A Φ [ = v R v R λ Φ E λ Φ Φ + λ Φ E λ A + + λ v R v R + + λ7 Φ A λ A A + λ 7 A Φ m + λ Φ + λ Φ + λ Φ λ ] A λ A λ A Φ λ A = m 4 A + Φ + λ v R v R + λ 7 Φ Φ 8 + λ9 A E

9 λ A Φ = m 4 A + + λ A Φ λ v R v R + λ 7 Φ + λ 7 Φ Φ λ 9A E λ8 E = m E + + λ9 A λ 9A λ Φ + λ Φ + λ Φ 4. (8) If we assume v R v R we obtain five quadratic equations and one linear equation for Φ i A i and E. For that case there are solutions. Two of them correspond to SU() symmetry and remaining solutions to G standard gauge group symmetry solutions. If we set v R = we find six quadratic equations for Φ Φ Φ A A and E with 4 solutions with symmetry groups having rank. They are isomorphic to G SU() C SU() L U() R U() B L. However there are solutions with higher symmetries. They are G G 4 G 4 and G. For general coupling constants λ λ m m 8 the solutions with higher symmetries are specified by following relations. Solutions with higher symmetries are characterized by:. SU() U() X and (SU() U()) flipped symmetry solutions { E = vr = Φ = ε Φ Φ = ε Φ A = ε A (9) where ε = and ε = correspond to the SU() U() X symmetric vacua and (SU() U()) flipped symmetric vacua respectively.. SU() symmetry solutions { E = Φ = Φ Φ = Φ A = A v R. (). G 4 SU(4) SU() L SU() R symmetry solutions { Φ = Φ = A = A = v R = Φ E. () 4. G SU() C SU() L SU() R U() B L symmetry solutions { Φ = A = v R = Φ Φ A E. () For reader s convenience we list the decompositions of each representation in Appendix A. 9

10 . G 4 SU(4) SU() L U() symmetry solutions { Φ = Φ = A = v R = Φ A E. (). G SU() C SU() L U() R U() B L symmetry solutions { vr = Φ i (i = ) A i (i = ) E. (4) The higher symmetry solutions given in Eqs. (9)-(4) lead to the crucial consistency checks for all results in this paper. 4 The states and Clebsch-Gordan coefficients 4. The states in the G multiplets In order to obtain and study the mass matrices it is convenient to decompose the Higgs representations under the G gauge group. The explicit decompositions of 4 4 and representations in Y diagonal basis are presented according to the G multiplets with the same quantum numbers which generally mix among themselves. The 8 singlets ( ) are already given in Eqs. (4) (). There are 4 = would-be NG modes. They are in the following multiplets: ( ) [( ) + c.c.] [( )+c.c.] [( )+c.c.] and [( )+c.c.]. The corresponding orthonormal states are listed in Table in Y -diagonal basis. The physically important modes so called Higgs doublets [( ) + c.c.] and color triplets [( ) + c.c.] states are listed in Tables and respectively. The remaining G multiplets are listed in Tables 4 and. There are altogether 9 states accommodated in G multiplets (see Table ). There are mass matrices containing NG modes containing doublets containing color triplets 9 containing the other modes. There are multiplets ( states) with zero mass and 9 G multiplets with masses different from zero (containing 9 = 8 states). Hence the mass spectrum contains 7 different mass eigenvalues. For the SU() solutions there are only different masses and G multiplets are grouped into multiplets transforming under the SU() group. For the G 4 solutions there are 7 different mass eigenvalues and G multiplets are grouped into multiplets transforming under the G 4 group (see Appendix A). The higher symmetries serve as a strong consistency check of our mass matrices and CG coefficients. We point out that the main basic blocks in all 9 states are SU() L irreps and

11 and SU() irreps and 8: SU() L [ + 4] () (4) (4) [ [ ] [ ] ][ { 4] [ 4 ] 4 4 } { 44 } 4 SU() C [ ] (79) (8) ( ) [ 7 9] (8 8) ( ) ([ 8 ] (79 9 7) ( ) ) [9 78] 8[ 9] [78 ] { } ( [9 78] 8 7[ 9] { 88 {8} {} {8} } [ [87 ] ) ] ( ) () () All states can be constructed combining and antisymmetrizing or symmetrizing the basic blocks. The basic blocks () and () which appear only in antisymmetric tensors are embraced by parenthesis the basic blocks which appear only in the symmetric tensors are embraced by curly brackets ({aa}/ = aa) while the basic blocks that appear both in symmetric and antisymmetric tensors are embraced by square brackets.

12 Table : States in the would-be NG modes ( ) + c.c. Â() () D () () () () Φ () () ( ) + c.c. Â( ) () D ( ) () ( ) () Φ ( ) () Φ ( ) () ( ) + c.c. Â( ) () Ê ( ) () Φ ( ) () Φ ( ) () ( ) + c.c. Â( ) () Ê ( ) () D ( ) () ( ) () ( ) () Φ ( ) () Φ ( ) () i () + c.c. (8) + c.c. 4 ([ + 4] 8) + c.c. 7 ( [ ]) + c.c. i (8 8) + c.c. ( [ 7 9]) + c.c. 4 ( [ [78 + 9]7[ + 9]9[ + 78]]) + c.c. 4 ( ) + c.c. 48 ([ + 4] [8 8]) + c.c. i ([ 4][ 7 9]) + c.c. {[ 4][ 7 9]} + c.c. 4 ([4 4][ 7 9]) + c.c. 48 ([ 4][ [78 + 9]7[ + 9]9[ + 78]]) + c.c. i ([ ][ 7 9]) + c.c. {[ ][ 7 9]} + c.c. ([ 4 ][8 8]) + c.c. 4 [ ([ 4 4][8 8]) + ([ 4 ][ ]) ] + c.c. [ ] 4 ([ 4 4][8 8]) ([ 4 ] [ ]) + c.c. 4 ([4 ][ 7 9]) + c.c. 48 ([ ][ [78 + 9]7[ + 9]9[ + 78]]) + c.c.

13 Table : States in the doublets [( ) + c.c.] ( ) + c.c. Ĥ( ) () [ ] + c.c. D ( ) () D ( ) () ( ) () ( ) () Φ ( ) () (4 ) + c.c. 8 ([ ][ ]) + c.c. [ ] 7 ([4 ][ ]) ([ ][ ]) + c.c. 7 [ ([4 ][ ]) + ([ ][ ]) ] + c.c. 4 ([ 4 ]8) + c.c. Table : States in the color triplets [( ) + c.c.] ( ) + c.c. Ĥ( ) () [ 7 9] + c.c. D ( ) () D ( ) () ( ) () ( ) () ( ) () Φ ( ) () ([ + 4] [ 7 9]) + c.c. ( [78 + 9] 7 [ + 9] 9 [ + 78]) + c.c. [ ] 4 (4 [ 7 9]) ( ) + c.c. [ ] 4 (4 [ 7 9]) + ( ) + c.c. ([ + 4][ [78 + 9]7[ + 9]9[ + 78]]) + c.c. 4 (4 [8 8]) + c.c.

14 Table 4: States in the other G multiplets including or ( ) + c.c. () () ( ) + c.c. Ê() () () () ( 4) + c.c. D ( 4 ) () ( 4 ) () ( 7 ) + c.c. D ( 7 ) () ( 7 ) () ( 7 ) () ( ) + c.c. D ( ) () ( ) () ( ) + c.c. Ê( ) ( ) ( ) () ( ) + c.c. D ( ) () ( ) () ( 4) + c.c. ( 4 ) () ( ) + c.c. ( ) () (8 ) + c.c. D (8 ) () (8 ) () (8 ) () (8) + c.c. [ {} ] + c.c ([ 4 4 ] 8 ) + c.c. (4 [ 7 9]) + c.c. 4 (4 [ [78 + 9]7[ + 9] 9 [ + 78]]) + c.c. ([ ][8 8]) + c.c. 4 ([4 ][8 8] [ ][ ]) +c.c 4 ([4 ][8 8] + [ ][ ]) +c.c ([ ] ) 4 4 [ 7 9] + c.c. ([ ] [ [78 + 9]7[ + 9]9[ + 78]]) + c.c. [ ] {79} {9} {7} + c.c. ( [ ]) [9 78] 8[ 9] [78 ] + c.c. ( [9 78] 8[ 9] [78 ] ) + c.c. ( [ ]) 4 [ + 4] [9 78] 8[ 9] [78 ] + c.c. ( [ ]) [9 78] + c.c. 8[ 9] ] ([ 4 4 [78 ] [ ]) [9 78] 8[ 9] [78 ] + c.c. ( [ ]) [ ] [9] + c.c. ( [ 4 [4 ] [9] ] + [ ][ ]) [78] c.c. ( [ ] 4 [4 ] [9] [ ][ ]) c.c [78]

15 Table : States in the other G multiplets including Φ ( ) + c.c. Φ ( ) () ( ) Â() () Ê () () Φ () () ( ) + c.c. Φ ( ) () ( ) + c.c. Φ ( ) () ( ) + c.c. Φ ( ) () ( ) + c.c. Φ ( ) () (8 ) Â(8) () Ê (8) ( ) Φ (8) () Φ (8) () 4 ([ ]8) + c.c. ( ) i 4 4 { } 4 4 ([ ] ) [ ] 4 ( [8 8]) + c.c. 4 ([ 4 4 ]) [8 8] + c.c. ( [ 4 [ 4] (9 78) ( [ 4 [ ] (9 78) 8( 9) ( i [ 78] (78 ) ]) + c.c. ]) 8( 9) (78 ) + c.c. [ + 78 [9]] ) { [ 78] [ + 78 [9]] } ( [9 789] [[78] 9 789]] ) ( [ 48 [ + 4] [ 78] [ + 78 [9]] ]) (8 ) + c.c. (8 ) Φ (8) () Φ (8) () ( [ [ 78] [ + 78 [9]] ]) + c.c. 4 ([ 4 4 ] [ [ 78] [ + 78 [9]] ])

16 Table : Number of G multiplets 8 ( ) 4 [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] ( ) [( ) + c.c.] [( 4 ) + c.c.] 7 [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] 4 [( ) + c.c.] 7 [( ) + c.c.] [( 7) + c.c.] [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] [( 4) + c.c.] [( ) + c.c.] [( ) + c.c.] [( ) + c.c.] 4 (8 ) [(8 ) + c.c.] [(8 ) + c.c.] (8 ) sum 74 multiplets NG multiplets 9 massive multiplets

17 4. H operators and Clebsch-Gordan coefficients Let us denote by R the sum of all representations 4 4 and R = I R I dim R = 9. (7) There are cubic invariants (see Eq. ()) I(R I R J R K ) I IJK R K R K (8) where R I R J = K R K. Let us denote H-operators (see Refs. [] and []) transforming as R K and H K (R I R J ) = H K (R J R I ) R K (9) H K (R I R J ) = N The normalization factor N is chosen so that For example in the Y -diagonal basis I IJK R K. (4) H K (R I R J ) R K = I IJK. (4) [H Φ (Φ Φ )] abcd = [ (Φ ) abef (Φ ) cdē f (Φ ) acef (Φ ) bdē f +(Φ ) adef (Φ ) bcē f + (Φ ) cdef (Φ ) abē f ] (Φ ) bdef (Φ ) acē f + (Φ ) bcef (Φ ) adē f (4) For invariants of the type I III there is only one H operator H I for invariants of the type I IIK there are two H operators H I and H K and for invariants of the type I IJK there are three H operators H I H J and H K. The H operators are symmetric in R I and R J H(R I R J ) = H(R J R I ). More generally I IJK = I JIK = I IKJ = I JKI = I KIJ = I KJI. (4) Specially we are interested in CG coefficients when at least one of the states transforms as a singlet under the G group. That requires the decomposition of each SO() irrep R I into G irreps R I = Ri I (44) i where indices i just enumerate the G irreps with fixed Y contained in a specific R I irrep. There are eight singlets denoted here shortly by S I i defined in Eqs. (4)-(). For 7

18 this choice of G -singlet states the action of H operators are reduced to the invariant subspaces of states with fixed G quantum numbers which are listed in Tables H K (ŜI i RJ j ) = k C IJK ijk RK k (4) where R J j and RK k transform as identical G irreps (have the same G quantum numbers). For example for the evaluation of the first column in the () mass matrix (see Eq. (74)) the following H operators have to be evaluated Note that H E (Ê Ê() () ) = H (v R Ê() () ) = Ê () () () (). (4) C IJK ijk = I(ŜI i ˆR J j ˆR K k ) ŜI i ˆR J j ˆR K k (47) (the second part of Eq. (47) defines shortand notation) where ˆR I i ˆR I i (48) is the complex conjugated irrep of ˆRI i. The CG coefficients are listed in tables in Appendix B. Note that the CG coefficients depend only on the indices i j that represent G multiplets and not on the specific states within the G multiplets. That can be used as a consistency check of the states belonging to the specific G multiplets. We choose all SO() invariants in Eq. () to be real except for invariants containing the or the field separately (for example HΦ or E ). In this last case the sum of the two invariants one containing field and the other having field in place of are real. Generally CG coefficients are complex (in our case they are either real or pure imaginary). Starting with cubic invariants in Eq. () our choice of phases of states in Tables is such that it leads to the minimal number of imaginary terms in the mass matrices. Symmetry relations imply C IJK ijk CG coefficients also satisfy hermiticity relations = C JIK jik. (49) C IJK ijk = C IKJ ikj. () Here the i represents the label assigned to the irrep complex conjugated to the irrep designated by i. 8

19 Furthermore following relations are valid j i i C IJJ ijj dim RJ j = j dim Rj J = dim R J () j C IJK ijk CIJK ijl dim R I i dim RJ j = C(I J K)δ kl dim R I i = dim R I where C(I J K) are constants depending on irreps R I R J and R K. dim Rj J = dim R J () j 9

20 Mass Matrices. Mass Matrices For the sum of all representations R (see Eq. (7)) we define fluctuations of the Higgs field around the VEVs as R = R + ri I. () Ii Then the matrix element of the mass matrix corresponding to G multiplets i j (R I i and R J j transform identically under G ) is: M IJ ij = M IJ ij W r I i r J j R= R (4) = m I δ IJ δ ij + λ IJK Si I CIJK ijk () Kk where λ III = λ p λ IIK = λ p and λ IJK = λ p for invariants containing three identical representations two identical representations (different from the third one) and three different representations respectively. Specifically {λ IJK } = {λ λ 8 λ λ 7 λ 9 λ λ λ λ λ 4 λ λ λ λ 4 λ λ λ 7 λ 8 λ 9 λ λ }. According to Eq. () the trace of the total mass matrix over all dim R states in R is Tr M = I m I dim R I. () The mass matrices are generally non-hermitian. So the squares of the physical masses are equal to the absolute values of eigenvalues of matrices M M and MM. (One can obtain hermitian matrices M M i.e. MM with the same spectra.) For a real superpotential that is for λ = λ 4 λ = λ λ 8 = λ 9 and λ = λ and all coupling constants and VEVs real the matrices are hermitian due to the hermiticity relation (). Now we are ready to present the explicit forms of the mass matrices calculated from the superpotential of Eq. (). Every matrix is designated with the corresponding G multiplet and appears dim Ri I times in the total mass matrix M. A mass matrix associated with a G multiplet and the mass matrix associated with the corresponding complex conjugated G multiplet are equal up to transposition and therefore for multiplets with Y we list only one of the two mass matrices. Of course when enumerating the total degrees of freedom one has to include all mass eigenvalues (9 in total). The basis designating the columns (c:) of the mass matrices listed below is given in the same way as shown in Tables 4 and while the rows are designated by the corresponding complex conjugated G multiplets (r:) :

21 ( ) c: Â() () Â() () Ê() () () () () () () () () Φ () Φ () Φ () r: Â() () Â() () Ê() () () () () () Φ () () () () Φ () Φ () where m () λ Φ λ 9A λ v R λ v R λ A λ 7 Φ m () 8 λ Φ m () λ 9A λ v R λ v R λ 7 Φ m () 7 m () 8 λ 9A λ 9A m () λ Φ λ Φ λ Φ λ v R λ v R m () 44 λ v R λ v R λ v R λ v R λ v R m () 44 λ v R λ v R λ v R λ A λ 7 Φ λ Φ λ v R λ v R m () λ 7 A λ Φ λ 7 Φ m () 7 λ Φ λ v R λ v R λ 7 A m () 77 m () 87 m () 8 m () 8 λ Φ λ v R λ v R λ Φ m () 87 m () 88 (7) m () m 4 + λ Φ + λ 9E m () m 4 + λ Φ m () m + λ 8E λ 9E m () 44 m + λ Φ + λ Φ + λ Φ m () m + λ E λ A λ A m () 7 λ A m () 77 m + m () 8 m () 8 + λ 7 Φ λ Φ λ E λ A + λ7 Φ λ A + 4 λ 7Φ m () 87 λ Φ + λ 7 A

22 [( ) + c.c.] m () 88 m + λ Φ λ Φ λ 7A + λ E. (8) c: Â() () () () D () () r: Â( ) () where m 4 + λ Φ ( ) ( ) D () () iλ 8 v R λ v R + λ 9 E λ A λ 7 Φ iλ 9v R m λ 4E + λ Φ iλ 9 A λ Φ λ v R m () m + λ Φ + λ Φ m () 44 m + λ Φ + λ Φ + 4 λ v R iλ 8A λ Φ m () λ v R λ A λ A λ 7 Φ λ v R λ v R m () 44 (9) λ 7A + λ E. () [ ( ) + c.c.] c: Â( ) () D ( ) ( () ) () Φ ( ) () Φ ( ) () r: Â( ) () D ( ) () ( ) Φ ( ) () () Φ ( ) () where m ( ) iλ 9 v R λ v R λ A λ 7 Φ λ A λ 7Φ iλ 8v R m ( ) m ( ) λ v R λ v R λ v R m ( ) m ( ) λ v R λ v R m ( ) λ Φ m 4 + m ( ) m + 4λ 4E + λ 9E m ( ) iλ 9A + λ Φ + λ Φ λ A λ 7 Φ λ v R λ v R m ( ) 44 m ( ) 4 λ Φ + λ Φ 9 λ A λ 7Φ λ v R λ v R m ( ) 4 m ( ) ()

23 m ( ) iλ 8A + λ Φ + λ Φ m ( ) m + λ Φ + λ Φ + λ Φ λ A λ A m ( ) 44 m + λ Φ λ E m ( ) 4 λ Φ + λ7 A m ( ) m + λ Φ + λ Φ + λ 7A + λ E. () [ ( ) + c.c.] c: Â( ) () Ê( ) () Φ ( ) () r: Â( ) () Ê( ) () Φ ( ) () Φ ( ) () Φ ( ) () where m 4 λ Φ + λ 9E λ 9A λ 9A λ A + λ A + λ 7Φ λ 7Φ λ 9A λ 9A m + λ 8E λ Φ + λ Φ 4 λ Φ + λ Φ λ A + λ 7Φ λ Φ + λ Φ 4 m λ Φ + 7λ E 4 λ Φ 4λ 7A λ 7Φ λ A + λ Φ + λ Φ λ Φ 4λ 7A m ( ) 44 () m ( ) 44 m + λ Φ λ Φ λ 7A 4 λ E. (4)

24 [ ( ) + c.c.] c: Â( ) () Ê( ) () D ( ) () ( ) ( () ) () Φ ( ) () Φ ( ) () r: Â( ) () Ê( ) where m ( ) λ 9 A λ 9A iλ 8 v R λ v R m ( ) m ( ) 7 () D ( ) () λ 9 A λ 9A m + λ 8E λ v R λ Φ 4 λ Φ λ Φ λ Φ ( ) () ( ) () Φ ( ) () Φ ( ) () iλ 9v R m ( ) m ( ) 4 m ( ) λ v R λ v R m ( ) m 4 + λ Φ m ( ) m λ 4E + λ v R m ( ) 4 m ( ) 44 λ E λ v R λ v R + λ 9E λ Φ 9 λ v R m ( ) λ E m ( ) + λ Φ 9 m ( ) λ Φ 4 λ Φ λ v R λ v R m ( ) λ Φ + 4λ 7A m ( ) 4 iλ 8A iλ 8A λ Φ 4 λ Φ λ Φ 4 m ( ) iλ 9A + iλ 9A λ Φ 4 + λ Φ λ Φ m ( ) 4 iλ 9A + iλ 9A λ Φ 4 λ Φ λ Φ 4 m ( ) 44 m + λ Φ + λ Φ λ A λ A m ( ) iλ 8A iλ 8A λ Φ 4 + λ Φ λ Φ m ( ) m + λ Φ + λ Φ + λ A m ( ) λ A m ( ) m + λ Φ m ( ) 7 λ A λ 7Φ + 7λ E 4 λ 7Φ 4 + λ A m ( ) 7 λ Φ λ Φ λ v R λ v R λ Φ + 4λ 7A m ( ) 77 ()

25 [ ( ) + c.c.] m ( ) 77 m + λ Φ + λ Φ + λ 7A 4 λ E () c: Ĥ( ) () D ( ) () D ( ) () ( ) ( () ) () Φ ( ) () r: Ĥ( ) () D ( ) () D ( ) () ( ) () ( ) () Φ ( ) () where m + λ E iλ A λ 7Φ iλ A λ 7Φ λ 4 Φ λ 4Φ λ Φ λ Φ λ v R iλ A λ 7Φ m + λ 4E λ Φ iλ 9 A iλ 8 A + λ Φ 4 + λ Φ 4 λ v R iλ A λ 7Φ λ Φ m ( ) m ( ) 4 m ( ) λ v R m ( ) m λ 4E + λ Φ 9 λ Φ λ Φ iλ 8A m ( ) 4 m ( ) 44 λ E + λ Φ 4 λ 4Φ λ 4Φ iλ 9A + λ Φ 4 m ( ) λ E m ( ) λ v R (7) λ 4v R λ v R λ v R λ v R m ( ) m ( ) 4 iλ 8A + iλ 8A + λ Φ 4 λ Φ m ( ) iλ 9A iλ 9A + λ Φ 4 + λ Φ m ( ) 4 iλ 9A iλ 9A + λ Φ 4 λ Φ m ( ) 44 m + λ Φ λ Φ + λ A m ( ) iλ 8A + iλ 8A + λ Φ 4 + λ Φ m ( ) m + λ Φ + λ Φ λ A m ( ) m + λ Φ + λ Φ + λ 7A 4 λ E. (8)

26 [ ( ) + c.c.] c: Ĥ( ) () D ( ) () D ( ) () ( ) () ( ) () ( ) () Φ ( ) () r: Ĥ( ) () D ( ) () D ( ) ( ) () () ( ) () ( ) Φ ( ) () () where m λ E m ( ) m ( ) m ( ) 4 m ( ) λ Φ λ v R m ( ) m ( ) λ Φ 9 m ( ) 4 m ( ) m ( ) λ v R m ( ) λ Φ 9 m ( ) m ( ) 4 m ( ) m ( ) λ v R m ( ) iλ A λ 7Φ m ( ) 4 m ( ) 4 m ( ) 4 m + λ A λ E m ( ) i λ A λ 7Φ m ( ) 4 λ Φ λ Φ m ( ) λ 4Φ λ 4Φ m ( ) iλ A λ 7Φ m ( ) m + 4λ 4E + m ( ) 4 iλ 8A + λ Φ m ( ) iλ 9A + λ Φ m ( ) iλ 9A + λ Φ m ( ) i λ A m ( ) m λ 4E + λ Φ λ 7Φ λ Φ 9 m ( ) m ( ) m ( ) λ E m λ A λ Φ λ v R λ 4Φ m ( ) m ( ) λ Φ m ( ) λ v R λ 4 v R λ v R λ v R λ v R λ v R m ( ) 77 (9)

27 [( ) + c.c.] m ( ) 4 iλ 8A λ Φ m ( ) iλ 9A + λ Φ m ( ) iλ 9A + λ Φ m ( ) 4 λ 4Φ λ 4Φ m ( ) 4 iλ 9A + λ Φ m ( ) 4 iλ 9A λ Φ m ( ) λ Φ λ Φ m ( ) iλ 8A + λ Φ m ( ) iλ 8A + λ Φ m ( ) iλ 8A + λ Φ m ( ) iλ 8A + λ Φ m ( ) m + λ Φ + λ Φ λ A m ( ) 77 m + λ Φ + λ Φ + λ Φ + λ 7A + λ E. (7) c: r: () () ( ) () m + λ Φ + λ Φ λ Φ + λ A λ A. (7) 7

28 [ ( ) + c.c.] c: Φ ( ) () r: Φ ( ) () ( ) c: Â() () Ê() () r: Â() () Ê() () Φ () () Φ () () m 4 λ Φ + λ 9E λ 9 A λ A + [( ) + c.c.] () c: Ê() () () r: Ê( ) () ( ) () [ ( 4 ) + c.c.] c: D ( 4 ) () 4 ( ) () r: D ( 4 ) () ( 4 ) () m + 4λ 4E m + λ Φ λ Φ λ 7A 4 λ E. (7) λ 9 A m + λ 8E λ7 Φ m + λ 8E λ v R + iλ 8 A + λ Φ λ Φ λ Φ λ Φ λ Φ 9 λ A + λ7 Φ λ Φ m λ Φ + λ Φ 4 λ v R m λ Φ + λ Φ + λ A λ 7A + λ E. (7). (74) iλ 9A + λ Φ λ Φ m + λ Φ + λ Φ λ Φ + λ A λ A. (7) 8

29 [ ( ) + c.c.] c: Φ ( ) () r: Φ ( ) () [ ( 7 ) + c.c.] c: D ( 7 ) () ( 7 ) 7 ( () ) () r: D ( 7 ) () where m + λ Φ + λ Φ λ Φ + λ 7A + λ E. (7) 7 ( ) () ( 7 ) () m ( 7 ) m ( 7 ) m ( 7 ) m ( 7 ) m ( 7 ) λ E m ( 7 ) m λ 4E + λ Φ 9 m ( 7 ) λ E m ( 7 ) λ Φ 9 m ( 7 ) iλ 8A + iλ 8A + λ Φ 4 + λ Φ λ Φ 4 (77) m ( 7 ) iλ 9A iλ 9A + λ Φ 4 λ Φ λ Φ m ( 7 ) iλ 9A iλ 9A + λ Φ 4 + λ Φ λ Φ 4 m ( 7 ) m + λ Φ λ Φ + λ A λ A m ( 7 ) iλ 8A + iλ 8A + λ Φ 4 λ Φ λ Φ m ( 7 ) m + λ Φ λ Φ λ A + λ A. (78) 9

30 [ ( ) + c.c.] c: D ( ) () ( ) () r: D ( ) () ( ) () [ ( ) + c.c.] c: Φ ( ) () r: Φ ( ) () λ Φ iλ 9 A + λ Φ m + 4λ 4E iλ 8A + λ Φ m λ Φ + λ Φ + λ A. (79) [ ( ) + c.c.] c: Ê( ) ( ) ( ) () r: Ê( ) ( ) ( ) () [ ( ) + c.c.] c: D ( ) ( () ) () r: D ( ) ( ) () () ( m λ 8E λ v R m λ Φ + λ Φ λ 7A + λ E. (8) m λ 4E λ Φ 9 iλ 8 A λ Φ λ v R m + λ Φ λ Φ + λ Φ + λ A + λ A iλ 9A λ Φ m + λ Φ λ Φ + λ A ). (8). (8)

31 [ ( 4 ) + c.c.] c: ( 4 ) () r: ( 4 ) () [ ( ) + c.c.] c: Φ ( ) () r: Φ ( ) () [ ( ) + c.c.] c: Φ ( ) () r: Φ ( ) () [ ( ) + c.c.] c: ( ) () r: ( ) () m + λ Φ λ Φ λ Φ λ A + λ A. (8) m λ Φ + λ Φ λ 7A 4 λ E. (84) m λ Φ λ Φ + λ 7A 4 λ E. (8) m λ Φ λ Φ λ A. (8)

32 (8 ) c: Â(8) () Ê(8) (8) (8) ( ) Φ () Φ () r: Â(8) () Ê(8) (8) (8) ( ) Φ () Φ () where m 4 λ Φ λ A λ 9E λ 9A + λ7 Φ λ A + λ 7Φ [(8 ) + c.c.] c: r: (8) Φ () (8 ) Φ () [ (8 ) + c.c.] c: D (8 ) () (8 ) (8 () ) () r: D (8 ) () λ 9A m λ 8E λ Φ λ Φ λ A + λ 7 Φ λ Φ m λ Φ λ E λ Φ λ7 A λ A + λ 7Φ λ Φ λ Φ λ7 A m (8) 44 (87) m (8) 44 m + λ Φ λ Φ λ 7A + λ E. (88) (8 ) () (8 ) () m + λ Φ λ Φ λ 7A + λ E. (89) where m (8 ) m (8 ) m (8 ) m (8 ) m λ 4E λ Φ 9 m (8 ) m (8 ) λ E m (8 ) λ E m (8 ) (9) m (8 ) iλ 8A iλ 8A + λ Φ 4 + λ Φ

33 m (8 ) iλ 9A + iλ 9A + λ Φ 4 λ Φ m (8 ) iλ 9A + iλ 9A + λ Φ 4 + λ Φ m (8 ) m λ Φ + λ Φ + λ A m (8 ) iλ 8A iλ 8A + λ Φ 4 λ Φ m (8 ) m λ Φ λ Φ λ A. (9) (8 ) c: r: (8) Φ () (8) Φ () m λ Φ λ Φ + λ 7A + λ E. (9). Tests and consistency checks of the total mass matrix The following consistency checks of the total mass matrix have been performed (see also Ref. ([7])):. The CG coefficients appearing in the total mass matrix have been found to satisfy the hermiticity relation ().. The trace of the total mass matrix has been evaluated and it has been found that it satisfies Eq. ().. For the G symmetric vacuum the number of the would-be NG modes have been found to be equal to the number of the broken generators i.e. massive gauge bosons 4 =. Here 4 represents the number of gauge bosons in the adjoint irrep of the SO() group and represents the number of the gauge bosons in the standard model (G group). The check has been performed numerically for several tens of randomly chosen sets of the parameters of the superpotential λ λ and m m constrained by VEV Eqs. (8). 4. For the SU() symmetric vacuum it has been found that the number of the different mass eigenvalues is. The number of independent SU() irreps contained in the total representation of the model R is (see Appendix A) and the would-be NG bosons are in two multiplets (()+( + )). Therefore there are in total mass

34 eigenvalues: different from zero and one equal to zero. All other states have to be accommodated in the SU() irreps. Further for the SU() symmetric vacuum all eigenvalues of the doublet Higgs matrix ( ) are contained in the spectrum of the triplet Higgs matrix ( ). The only mass eigenvalue of the triplet Higgs matrix spectrum not contained in the doublet mass spectrum belongs to the SU() multiplet leading to the following relation between determinants of these two matrices det(m ( ) λ ) = (m ( ) λ) det(m ( ) λ ). (9) The above relation gives a very strong test for these two matrices. All above checks have been performed numerically in the same way as explained in the previous item. Our mass matrices passed all checks. For the SU() symmetric vacuum all VEVs except E are different from zero and they are non trivially correlated through the VEV Eqs. (8) and SU() symmetry conditions (). Therefore these tests serve as a strong check of the total mass matrix up to terms which depend on E.. For the G 4 symmetric vacuum it has been found that the number of different mass eigenvalues is 7 what is just the number of the independent G 4 irreps (see Appendix A). The number of NG states contained in G 4 irrep () is 4 = 4. These checks have been performed numerically as explained above. All other states have been found to be accommodated in G 4 irreps. For the G 4 symmetric vacuum VEVs different from zero are Φ and E and therefore this serves as a check of the E-dependent parts of the total mass matrix.. Similar tests are satisfied for other higher symmetries G G 4 G and G. We stress that all mass matrices in Ref. [7] derived for minimal SO() model (R = ) satisfy all above consistency checks and are just a special case of the mass matrices in this paper. Further all results including CG coefficients and mass matrices were obtained analytically including checks. In the recent calculations in Refs. [9 ] that appeared after Ref. [7] the necessary condition (9) between doublets and triplets is not satisfied. The hermiticity condition () and the total trace relation () are also not satisfied in these references. Further none of the higher symmetry tests is satisfied. That is a consequence of different phase conventions in Ref. [7] and phase conventions in Refs. [9] and []. In this paper all results for the CG coefficients and mass matrices have been also obtained analytically. Checks have been performed numerically and it has been found that the mass matrices satisfy ALL consistency checks. 4

35 Mass Matrices of Quarks and Leptons After the symmetry breaking down to the G subgroup the electroweak symmetry breaking SU() L U() Y U() em can be achieved by the VEVs of doublets included in the fields H D and Φ. These fields are given by (see Table ) H u/d = H (± ) () D u/d = D (± ) () D u/d = D (± ) () u/d = (± ) () u/d = (± ) () Φ u/d = Φ (± ). () The Yukawa couplings of Eq. () including these doublets can be written as follows: ( W Y = Ui c Y ij H u + Y ij D u + Y ij D u + Y ij ) u Qj ( + Di c Y ij H d + Y ij D d + Y ij D d + Y ij ) d Qj ( + Ni c Y ij H u + Y ij D u Y ij D u Y ij ) u Lj ( + Ei c Y ij H d + Y ij D d Y ij D d Y ij ) d Lj. (94) where U c D c N c and E c are the right-handed SU() L singlet quark and lepton superfields Q and L are the left-handed SU() L doublet quark and lepton superfields respectively. This is a generalization of the renormalizable minimal SO() model [ 7] including. Note that the successful gauge-couplings unification is realized with only the MSSM particle contents. This means that only one pair of Higgs doublets remains light and others should be heavy ( M G ). Here we accept the simple picture that the low-energy superpotential is described by only one pair of light Higgs doublets (H u and H d ) in the MSSM. But in general these Higgs fields are admixtures of all Higgs doublets having the same quantum numbers in the original model such as: H u = α u H u + α u D u + α u D u + α 4 u u + α u u + α u Φ u H d = α d H d + α d D d + α d D d + α 4 d d + α d d + α d Φ d (9) where α ud i (i = ) denote elements of the unitary matrix which rotate the flavor basis in the original model into the (SUSY) mass eigenstates. As mentioned above the low-energy superpotential is described only by the light Higgs doublets H u and H d ( W Y = Ui c α u Y ij + αu Y ij + αu Y ij + αu 4 Y ) ij Hu Q j ( + Di c α d Y ij + αdy ij + αdy ij + αdy) 4 ij Hd Q j + N c i + E c i ( α u Y ij + αu Y ij αu Y ij αu 4 Y ) ij Hu L j ( α d Y ij + αd Y ij αd Y ij αd 4 Y ) ij Hd L j (9) where the formulae of the inverse unitary transformation of Eq. (9) H u = αu H u D u = αu H u D u = αu H u u = αu 4 H u H d = αd H d D d = αd H d D d = αd H d d = αd 4 H d (97)

36 have been used. Providing the Higgs VEVs H u = v sin β and H d = v cosβ with v 74. [GeV] the quark and lepton mass matrices can be read off as M u = c M + c M + c M + c M M d = M + M + M + M M D = c M + c M c M c M M e = M + M M M M R = c R M (98) where M u M d M D M e and M R denote the up-type quark down-type quark neutrino Dirac charged-lepton and right-handed neutrino Majorana mass matrices respectively. Note that the mass matrices at the right-hand side of Eq. (98) are defined as M = Y α d v cosβ M = Y α d v cos β M = Y α d v cosβ and M = Y α 4 d v cosβ respectively and the coefficients are defined as c = (α u /α d ) tanβ c = (α u /α d ) tanβ c = (α u/α d) tan β c = (α 4 u/α 4 d) tanβ and c R = v R /(α 4 dv cosβ). These mass matrices are directly connected with low-energy observations and are crucial to model judgement. 7 Conclusion We have presented a simple method for the calculation of CG coefficients. We have constructed all states for all antisymmetric and symmetric SO() tensor irreps. We list all tables for the CG coefficients for the SO() irreps 4 4 and for all possible cubic invariants. We have constructed all mass matrices for the corresponding Higgs-Higgsino sector is SUSY GUT SO() models. We have found a set of consistency checks for the CG coefficients and mass matrices which proved the correctness of all our results. The results obtained here are useful for a wide class of GUT models based on the SO() group. Acknowledgements The work of A.I. and S.M. is supported by the Ministry of the Science and Technology of the Republic of Croatia under contracts No. 9 and No. 98 respectively. The work of T.F. T.K. and N.O. is supported by the Grant in Aid for Scientific Research from the Ministry of Education Science and Culture. The work of T.K. is also supported by the Research Fellowship of the Japan Society for the Promotion of Science for Young Scientists. T.F would like to thank G. Senjanović and A.Y. Smirnov for their hospitality at ICTP.

37 Appendix A Decomposition of Representations under G Here we list the decompositions of 4 4 and representations under the chain of subgroups G 4 G G G in Tables and where U() groups are related to SU(4) SU() C U() B L SU() L SU() R SU() L U() R. 7 Note also that we may consider another chain of subgroups SO() SU() U() X and SU() G. The relations between the generators of U() XY and U() B LR are: where T R denotes the U() R generator. ( X + 4 Y ) = B L Y = B L + TR (99) 7 We use the same notation as Slansky [] but with the proper U() normalizations. 7

38 Table 7: Decomposition of the representation (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) ( ) ( ) ( ) () ; ; ; ( ) ( ) ( ) ( ) ( ) ; ; ; ( () (; ) ; ( ) ) ; ( ) ( ; ( ) ( ) ) ; Table 8: Decomposition of the representation (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) ( ) ( ) ( ) (4) ; ; ; ( ) ( ) ( ) ( ) ( ) ; ; ; ( ) ( ) ( 4 ; ; ) ( ) ( ) ; ( ; ) ( ) ; ( ) ( ) ( ; ; ) (; ) ( ) ( ) ; (; ) ( ) 8

39 Table 9: Decomposition of the representation 4 (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) () (; ) (; ) (; ) ( 4) (; ) (; ) ( ( ) ) (; ) (; ) 4 () (; ) (; ) (; ) (4 ) ( ) ( () ; ; ) ( ) ; ( 4) ( ; ) ( ) ; (4 ) ( ) ( ; ; ) ( ) ; (4 ) ( ) ( ) ( ) ; ; 4 () ( (; ) ) ( (; ) ) ( (; ) ) ( (4 ) ) ; ; ; 4 ( ) ( ) ( ) ; ; ; ( 4) (8; ) (8; ) (8; ) (4 ) Table : Decomposition of the representation 4 (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) () (; ) (; ) (; ) (4 ) () (; ) (; ) (; ) ( 4) (; ) (; ) ((4 ) ) (; ) (; ) 4 ( ) ( ) ( ) ( ) ( ) ; ; ; 4 ( ) ( ) ( ) ; ; ; ( 4) (8; ) (8; ) (8; ) (4 ) ( ) ( () ; ; ) ( ) ; ( 4) ( ; ) ( ) ; (4 ) ( ) ( ; ; ) ( ) ; (4 ) ( ) ( ) ( ) ; ; 4 9

40 Table : Decomposition of the representation (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) ( () (; ) ; ( ) ) ; ( ) ( ; ( ) ( ) ) ; ( ) () (; ) (; ) (; ) ( ) ( ) ( ) ; ; ; ( ) ( ) ( ) ( ) ( ) ; ; ; 4 ( ) (; ) (; ) (; ) ( ) ( ) ( ) ( ) ( ) ; ; ; ( ) ( ) ( ) ; ; ; (4 ) ( ) ( ) ( ) () ; ; ; (4 ) ( ) ( ) ( ) ( ) ; ; ; 4 ( ) ( ) ( ) ( ) () ; ; ; ( ) ( ) ; ; (4 ) ( ) ( ) ( ) ; ; 4 4 ( ) ( ) ( ) ; ; ; 4 (4 ) ( ) ( ) ( ) ; ; 4 ( ) ( ) ; ; ( ) ( () (; ) ; ( ) ) ; (4 ) ( ; ( ) ( ) ) ; 4 ( ) ( ; ; ) ( ) ( ) ; 7 4 ( ; ) ( ) ; ( ) ( ) ( ; ; ) ( ) ; 7 (4 ) ( ; ) ( ) ( ) ; ( (8; ) 8; ( ) ) 8; (4 ) ( 8; ( ) ( ) ) 8; 4 4

41 Table : Decomposition of the representation (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) ( ) ( ) ( ) () ; ; ; ( ) ( ) ( ) ( ) ( ) ; ; ; 4 ( ) (; ) (; ) (; ) ( ) ( ) ( ) ( ) ( ) ; ; ; 4 ( ) ( ) ( ) ; ; ; ( ) () (; ) (; ) (; ) ( ( ) ) (; ) (; ) ( ) (; ( ) ) (; ( ) ) ( ( ) ) ; ; ; ( ) ( ) ; ; ( ) ( ) ( ) ( ) ; ; 4 4 ( ) ( ) ( ) ; ; ; 4 ( ) ( ) ( ) ( ) ; ; 4 ( ) ( ) ; ; ( ) ( () (; ) ; ( ) ) ; ( ) ( ; ( ) ( ) ) ; 4 ( ) ( ; ; ) ( ) ( ) ; 7 4 ( ; ) ( ) ; ( ) ( ) ( ; ; ) ( ) ; 7 ( ) ( ; ) ( ) ( ) ; ( (8; ) 8; ( ) ) 8; ( ) ( 8; ( ) ( ) ) 8; 4 4

42 Table : Decomposition of the representation (4 L R ) ( C L R B L ) ( C L R B L ) ( C L Y ) ( X ) () (; ) (; ) (; ) ( ) () ((; ) ) ((; ) ) ((; ) ) ((4 ) ) ; ; ; 4 ( ) ( ) ( ) ; ; ; ( 4) (8; ) (8; ) (8; ) (4 ) ( ) ( () ; ; ) ( ) ; ( 4) ( ; ) ( ) ; (4 ) ( ) ( ; ; ) ( ) ; (4 ) ( ) ( ) ( ) ; ; 4 () ((; ) ) ((; ) ) ((; ) ) (4 ) ; ; ; (4 4) ( ) ( ) ( ) ( ) ; ; ; 4 4 (8; ) (8; ) (8; ) (7 ) () (; ) (; ) (; ) ( 4) (; ) (; ) ( (7 ) ) (; ) (; ) 4 ( ) ( ) ( ) ; ; ; (7 ) ( ) ( ) ; ; (4 4) ( ) ( ) ; ; ( 8) ( ) ( ) ( ) ( ) ; ; ; 8 ( ) ( ) ( ) ; ; 4 4 ( ) ( ) ; ; (7 ) ( ) (8; ) (8; ) (8; ) 4 4 (8; ) (8; ) (7 ) (8; ) (8; ) (4 4) 4

43 ( () (; ) ; ( ) ( ) ) ; 8 ( ; ( ) ( ) ) ; 4 4 ( ) ( ; ; ) ( ) ( ) ; 4 4 ( ; ) ( ) ; (7 ) ( ) ( ; ; ) ( ) ; (7 ) ( ; ) ( ) ; (4 4) ( ) ( (; ) ; ( ) ) ; (4 4) ( ; ( ) ) ; ( 8) ( ) ( ; ; ) ( ) ; (7 ) ( ; ) ( ) ; (4 4) ( ) ( ; ; ) ( ) ( ) ; 4 4 ( ) ( ) ; ; (7 ) B CG coefficients The CG coefficients for HH AA EE DD ΦΦ EA AΦ E E EΦ H H HΦ Φ Φ A A DH DA D D and DΦ combinations are listed in Tables and 8 respectively. Table 4: CG coefficients for the HH combination of fields H H ÊHH Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () 4

44 Table : CG coefficients for the AA combination of fields A A ÊAA AAΦ AAΦ AAΦ Â () () Â () () Â () () Â () () Â () () Â ( ) () Â ( ) () Â ( ) () Â (8) () Â () () Â () () Â () () Â ( ) () Â () () Â ( ) () Â ( ) () Â ( ) () 4 4 Â (8) () Table : CG coefficients for the EE combination of fields E E ÊEE Ê () () Ê () () Ê () () Ê ( ) () Ê ( ) () Ê ( ) ( ) Ê (8) ( ) Ê () () Ê () () Ê ( ) () Ê ( ) () Ê ( ) () 4 4 Ê ( ) ( ) Ê (8) ( ) 44

45 Table 7: CG coefficients for the DD combination of fields D D ÊDD DD Φ DD Φ DD Φ D ( ) () D () () D ( ) () D ( ) () D ( ) () D ( ) () D ( 4 ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( 7 ) () D ( ) () D ( ) () D (8 ) () D ( ) () 9 D ( ) () D ( ) () D ( ) () D ( 4 ) () D ( ) () D ( ) () 9 9 D ( ) () 9 D ( ) () D ( ) () D ( ) () D ( 7 ) () D ( ) () 9 9 D ( ) () 9 D (8 ) ()

46 Table 8: CG coefficients for the combination of fields Ê () () ( ) () ( ) () ( ) () ( ) () ( 7 ) () (8 ) () () ()) ( ) () ( ) () ( ) () ( ) () ( 7 ) () (8 ) () Table 9: CG coefficients for the combination of fields Ê () () ( ) () ( ) () ( ) () ( ) () ( 7 ) () (8 ) () () ()) ( ) () ( ) () ( ) () ( ) () ( 7 ) () (8 ) () 4

47 Table : CG coefficients for the combination of fields iâ iâ Φ Φ Φ () () ( ) () ( ) () ( ) () ( ) () () () ( 4 ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( 7 ) () ( 7 ) () ( ) () ( ) () ( ) () ( 4 ) () ( ) () (8 ) () (8 ) () () () () () () () ( ) () ( ) () ( ) () ( 4 ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( 7 ) () ( 7 ) () ( ) () ( ) () ( ) () ( 4 ) () ( ) () (8 ) () (8 ) () 47

48 Table : CG coefficients for the ΦΦ combination of fields Φ Φ εâφφ Φ () () Φ () () Φ () () Φ () () Φ () () Φ () () Φ () () Φ ( ) () Φ ( ) () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () εâφφ ÊΦΦ Φ ΦΦ Φ ΦΦ Φ ΦΦ () Φ () Φ () () Φ () () Φ () () Φ () () 9 Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () 4 4 Φ ( ) () 8 Φ ( ) () Φ ( ) () Φ ( ) () 7 Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () 7 Φ ( ) ()

49 Φ Φ εâφφ Φ ( ) () Φ ( ) () Φ (8) () Φ (8) () Φ (8) () Φ (8) () Φ (8) () Φ ( ) () 8 Φ ( ) () εâφφ ÊΦΦ Φ ΦΦ Φ ΦΦ Φ ΦΦ Φ (8) () Φ (8) () Φ (8) () Φ (8 ) () Φ (8) () Table : CG coefficients for the EA combination of fields E A EAÂ EAÂ Ê () () Ê () () Ê () () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê (8) ( ) Â () () Â () () Â () () Â ( ) () 4 Â ( ) () 4 Â ( ) () Â ( ) () 4 4 Â (8) () 49

50 Table : CG coefficients for the AΦ combination of fields A Φ ÂAΦ ÂAΦ εaφφ εaφφ Â () () Â () () Â () () Â () () Â () () Â () () Â () () Â ( ) () Â () () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â ( ) () Â (8) () Â (8) () Φ () () Φ () εaφφ () Φ () () Φ () () Φ () () Φ () () Φ ( ) () Φ () () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ (8) () Φ (8) ()

51 Table 4: CG coefficients for the E combination of fields E E Ê () () Ê () () Ê ( ) () Ê ( ) () Ê ( ) ( ) () () ( ) () ( ) () ( ) () ( ) () Table : CG coefficients for the E combination of fields E E Ê () () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) ( ) () () () () ( ) () ( ) () ( ) ()

52 Table : CG coefficients for the EΦ combination of fields E Φ EΦΦ EΦΦ EΦΦ Ê () () Ê () () Ê () () Ê () () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê ( ) () Ê (8) ( ) Ê (8) ( ) Φ () () Φ () () Φ () () 4 Φ () () 4 Φ ( ) () 4 Φ ( ) () Φ ( ) () 4 Φ ( ) () 4 Φ ( ) () Φ ( ) () 4 Φ ( ) () 4 Φ ( ) () Φ (8) () Φ (8) () Table 7: CG coefficients for the H combination of fields H H Φ H Φ H Φ ( ) () ( ) () ( ) () ( ) () ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) ()

53 Table 8: CG coefficients for the H combination of fields H H Φ H Φ H Φ ( ) () ( ) () ( ) () ( ) () ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Table 9: CG coefficients for the HΦ combination of fields H Φ Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () Ĥ ( ) () HΦ HΦ Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) ()

54 Table : CG coefficients for the Φ combination of fields Φ Φ Φ () () Φ () () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () () () () () () () () () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () 4

55 Table : CG coefficients for the Φ combination of fields Φ Φ Φ () () Φ () () Φ () () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () () () () () () () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) () ( ) ()

56 Table : CG coefficients for the A combination of fields A ia  () ()  () ()  ( ) ()  ( ) ()  ( ) ()  ( ) () () () () () () () ( ) () ( ) () ( ) () Table : CG coefficients for the A combination of fields A ia  () ()  () ()  () ()  ( ) ()  ( ) ()  ( ) () () () () () ( ) () ( ) () ( ) () ( ) ()

57 Table 4: CG coefficients for the DH combination of fields D H DH DH DH Φ DH Φ DH Φ D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () Ĥ ( ) () i Ĥ ( ) () i Ĥ ( ) () i Ĥ ( )) () i Ĥ ( ) () i Ĥ ( ) () i Ĥ ( ) () i Ĥ ( ) () i Table : CG coefficients for the DA combination of fields D A D A D A D () () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) ()  ( ) () i  () ()  ( ) () i  ( ) () i i  ( ) () i  ( ) () i 7

58 Table : CG coefficients for the D combination of fields D D Â D Â D Φ D Φ D Φ () () i D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( 4 ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( 7 ) () D ( 7 ) () D ( ) () D ( ) () D (8 ) () D (8 ) () ( ) () i ( ) () i ( ) () i ( ) () i i 4 i ( 4 ) () i ( ) () i ( ) () ( ) i ()) i ( ) () i ( ) () i ( ) () i ( ) () i ( ) () i ( ) () i ( 7 ) () i ( 7 ) () i i i i i ( ) () i ( ) () (8 ) () i (8 ) () 4 4 i i i i 4 4 8

59 Table 7: CG coefficients for the D combination of fields D D Â D Â D Φ D Φ D Φ ( ) () i D () () D ( ) () D ( ) () D ( ) () D ( ) () D ( 4 ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( 7 ) () D ( 7 ) () D ( ) () D ( ) () D (8 ) () D (8 ) () ( ) () i ( ) () ( ) () i ( ) () i i ( 4 ) () i ( ) () i ( ) () i i 4 i ( ) () i ( ) () i ( ) () i ( ) () i ( ) () i ( ) () i ( ) () i ( 7 ) () i ( 7 ) () i i i i i ( ) i () 4 4 ( ) () i (8 ) () i (8 ) () i i i 4 4 9

60 Table 8: CG coefficients for the DΦ combination of fields D Φ D Φ D Φ Φ ( ) () D () () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () D ( ) () Φ () () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) () Φ ( ) ()