arxiv: v1 [hep-ph] 24 Jul 2009

Transcription

1 OSU-HEP July 4, 009 Flavor Violation in Supersymmetric Q 6 Model arxiv: v1 [hep-ph] 4 Jul 009 K.S. Babu 1 and Yanzhi Meng Department of Physics Oklahoma State University Stillwater, OK 74078, USA Abstract We investigate flavor violation mediated by Higgs bosons and supersymmetric particles in a predictive class of modelsbased on thenon Abelian flavor symmetryq 6. Thesemodels, which aim to reduce the number of parameters of the fermion sector and to solve the flavor changing problems of generic SUSY setup, assume three families of Higgs bosons and spontaneous/soft violation of CP symmetry. Tree level contributions to meson antimeson mixings mediated by Higgs bosons are shown to be within experimental limits for Higgs masses in the (1-5) TeV range. Calculable flavor violation induced by SUSY loops are analyzed for meson mixing and lepton decays and found to be consistent with data. Significant new SUSY contributions arise in B s B s mixing, but non-standard CP violation is suppressed. A simple solution to the SUSY CP problem is found, which requires light Higgsinos. 1 babu@okstate.edu yanzhi.meng@okstate.edu 1

2 1 Introduction The gauge interactions of the standard model (SM) fermions are invariant under separate U(3) L U(3) R transformations. Thisglobalsymmetry isbrokenexplicitly bythefermionyukawacouplings. For the light fermions violation of this symmetry is small, being proportional to their masses. This feature has played a crucial role in the success of the SM in the flavor sector. In extensions of the SM this property is generally lost, often leading to excessive flavor changing neutral current (FCNC) processes. A case in point is the supersymmetric standard model which is the subject of this paper. While the gauge interactions of the SUSY SM respect the U(3) L U(3) R global symmetry, there are new sources of violation of this symmetry, in the soft SUSY breaking sector. Indeed, generic soft SUSY breaking scenarios lead to excessive FCNC in processes such ad K 0 K 0 mixing, B 0 B 0 mixing, D 0 D 0 mixing, and flavor changing leptonic decays such as µ eγ [1]. This problem is most severe in the K 0 K 0 system. SUSY box diagrams involving gluino and squarks modify the successful SM prediction for M K and ǫ K, leading to the following constraints for the real and imaginary parts of the amplitude []: ( ) (Re, Im)(δ d LL ) 1 (δrr d ) 1/ 1 ( , m ). (1.1) 500 GeV Here (δ AB ) ij = (m AB ) ij/ m is a flavor violating squark mass insertion parameter, for (A,B) = (L,R), with m being the average mass of the relevant squarks ( d and s in this case). For this estimate the gluino mass was assumed to equal the average squark mass. Now, the natural magnitude of the mixing parameters (δ d LL ) 1 and (δ d RR ) 1, in the absence of additional symmetries, should be of order the Cabibbo angle, 0.. Since the parameters (δ AB ) ij split the masses of the squarks, one sees from Eq. (1.1) that a high degree of squark mass degeneracy is needed for consistency. Analogous limits from Bd 0 B0 d mixing are less severe, as given by [3]: ( ) (Re, Im)(δLL d ) 13(δRR d ) 1/ 13 (.1 10, m ). (1.) 500 GeV Note that the natural value of this mixing parameter, in the absence of other symmetries, is V ub The constraints from Eq. (1.) are well within limits. B s B s mixing provides even weaker constraints. It can be argued that a natural explanation for solving this problem is to enhance the symmetry of the SUSY SM by assuming a non Abelian symmetry G (a subgroup of the U(3) L U(3) R ) that pairs the first two families into a doublet, with the third family transforming trivially [5]. 3 3 Grouping all three families into an irreducible triplet representation of G is also possible. The large top quark mass however reduces the original U(3) L U(3) R symmetry to U() L U() R, so we find it is easier to work with ( + 1) assignment.

3 Invariance under G will then lead to degeneracy of squarks, as needed for phenomenology. A variety of such models have been proposed in the literature [5], [6], [7], [8], [9]. In Ref. [5], SU() family symmetry and its variants were proposed to solve the SUSY FCNC problem. If the symmetry is global, one has to deal with the Goldstone bosons associated with its spontaneous breaking. Global symmetries are susceptible to violations from quantum gravity. Local gauge symmetries are more natural, but in the SUSY context there would be new FCNC processes arising from the family SU() D terms [10]. Exceptions to this generic problem are known to exist [6]. A more natural solution to the problem is perhaps to choose G to be a non Abelian discrete symmetry group [7]. In this case there would be no D term problem, since there are no gauge bosons associated with G. Spontaneous breaking of such symmetries will not lead to Goldstone bosons. If the symmetry breaking occurs before the inflationary era, such models should also be safe from potential cosmological domain wall problems. Such non Abelian discrete symmetries have found application in understanding the various puzzles associated with the quark and lepton masses and mixing angles with or without supersymmetry [11], more recently for understanding the tri bimaximal neutrino mixing pattern [1]. It would be desirable to find a symmetry that sheds light on the fermion mass and mixing puzzle, and at the same time solves the SUSY FCNC problem. The supersymmetric standard model has another problem. In the flavor conserving sector CP violation is generically too large. Neutron and electron electric dipole moments (EDM) receive new contributions from SUSY loops. Unless the new phases in the SUSY breaking sector are small or conspire to be small, experimental limits on the EDM of the neutron (d n ), electron (d e ), and atoms will be violated by two to three orders of magnitude (depending on the squark and slepton masses) [13], [14]. The imaginary parts of the left right squark mixing parameters must satisfy the constraints (from the experimental constraints d n < e-cm, d e < e-cm) [15] ( ) ( ) Im[(δLR d ) 11] m, Im[(δLR e 500 GeV ) 11] m, (1.3) 100 GeV assuming that the gluino/bino has the same mass as the squark/slepton. Now, since these mixing parameters are expected to be suppressed by fermion helicity factors (but enhanced by the MSSM parameter tanβ) the natural values for these mixing parameters are of order (1 10 4, ) respectively, (for tanβ = 10 and assuming order one phases). This implies that the CP violating phases arising from the soft SUSY breaking sector must satisfy θ d 1/53, θ e 1/63 (for gluino (Bino) mass of 500 GeV (100 GeV). Why this is so, while the Kobayashi-Maskawa phase takes order one value, is the SUSY CP puzzle. It would be desirable to resolve this puzzle based on a symmetry principle in the same context where the SUSY FCNC problem is solved. ThepurposeofthispaperistostudyarecentlyproposedSUSYmodelbasedonthenon Abelian symmetry group Q 6 [8] which addresses these issues. Q 6 is a finite subgroup of SU() with twelve elements. Apart from providing a solution to the SUSY flavor problem, this class of models can 3

4 also constrain the quark masses and mixings. It was shown in Ref. [8] that with the assumption of spontaneous (or soft) CP violation, there is a non-trivial relation between quark masses and mixings in this model. This sum rule was found to be consistent with experimental data. A crucial aspect of the Q 6 model relevant for the quark mixing sum rule is that CP violation occurs either spontaneously or softly. This can help ameliorate the SUSY CP problem mentioned above. CP invariance requires that the gaugino masses, the µ terms and the trilinear A terms be all real. In the Q 6 model of Ref. [8] it was found that there is a phase alignment mechanism that makes the phases of the sfermion mixing terms arising from the A terms to align with the phases of the fermion masses. So SUSY CP violation is suppressed to a large extent. However, spontaneously induced complex VEVs do lead to non-zero contributions to EDM. Here we analyze these contributions. Since these complex VEVs are accompanied by the Higgsino µ terms, a simple solution to the problem is found by making the Higgsinos to be lighter than the squarks. Adequate suppression of EDM is obtained for µ 100 GeV, while squark masses are of order 500 GeV. This suggestion obviously has testable implications for physics that will be probed at the LHC. The fermion mass matrices that allow for a non-trivial prediction and the phase alignment is a generalization of well studied models [16]. The mass matrices for up and down quarks and the charged leptons take the following form: M = 0 C 0 ±C 0 B 0 B A. (1.4) The main feature of such mass matrices is that the phases can be factorized, i.e., M = P M 0 Q, with M 0 being real and P,Q being diagonal phase matrices. This feature, when combined with the Q 6 symmetry, has an the interesting consequence that CP violation induced by SUSY loops are suppressed. This will be discussed in more detail in Sec. 4. The form of Eq. (1.4) can be obtained in renormalizable theories based on Q 6 symmetry. This requires the introduction of three families of Higgs doublets, which fall into + 1 representations of the Q 6 group, very much like the quarks and leptons. With multiple Higgs fields coupling to fermions, invariably there will be tree-level FCNC mediated by the Higgs bosons. The flavor changing Higgs couplings are not arbitrary, but can be computed in terms of the fermion masses and mixings. We will show that these FCNC processes are within acceptable range, provided that the Higgs boson masses lie in the (1 5) TeV range (except of course for the standard model like Higgs boson, which has a mass in the ( ) GeV range). While Higgsinos are naturally light in this scenario, in the bosonic sector only the lightest SM like Higgs will be accessible to LHC experiments. One of our major results is that non standard CP violation is highly suppressed in this class of models. The phase factorizability of the fermion mass matrices implies that much of the SUSY induced CP violation is small. The structure of the Yukawa couplings in the model implies that the 4

5 amplitudes for tree level FCNC induced by neutral Higgs bosons are nearly real (see discussions in Sec. 5). While there can be significant new contributions to meson antimeson mixings, there is very little CP violation beyond the standard model. Our analysis is similar in spirit to that of Ref. [9]. Our approach is slightly different, with some differences in analytical results, fits, spectrum, and conclusions. In particular, we have presented complete analytical results for the Higgs boson spectrum, and we have a new proposal to solve the SUSY EDM problem, which requires light Higgsinos. We have also derived generalized constraints on SUSY FCNC parameters for the B d,s B d,s system appropriate for a (+1) mass spectrum. The plan of the paper is as follows. In Sec. we describe the SUSY Q 6 model, lay out the parameter choice, and summarize the prediction for the quark sector. In Sec. 3 we analyze the Higgs potential involving the three pairs of Higgs doublets. We provide analytic expressions for the mass spectrum of Higgs bosons as well as numerical fits. Consistency of symmetry breaking and spontaneous CP violation will be established here. In Sec. 4 we address tree level FCNC processes mediated by the heavy Higgs bosons. Sec. 5 is devoted to analysis of the SUSY flavor violation and EDM within the model. In Sec. 6 we conclude. Supersymmetric Q 6 Model Q 6 is the binary dihedral group, a subgroup of SU(), of order 1. It has the presentation {A,B;A 6 = E,B = A 3,B 1 AB = A 1 }. (.5) The 1 elements of Q 6 can be represented as {E,A,A,...,A 5,B,BA,BA,...,BA 5 }. (.6) In the two dimensional representation the generators are given in a certain basis by ( ) ( ) A = cos π sin π i B = sin π cos π 0 i 3 3. (.7) The irreducible representation of Q 6 fall into,,1,1,1,1, where the is complex valued but pseudoreal, while the is real valued. (Q 6 is the simplest group with two distinct doublet representations, which is very useful for model building.) The 1 and 1 are real representations, while 1 and 1 are complex conjugates to each other. The group multiplication rules are given as 1 1 = 1, 1 1 = 1, 1 1 = 1, 1 1 = 1, 1 1 = 1, 1 1 = 1 (.8) 1 =, 1 =, 1 =, 1 =, 1 =, 1 = (.9) = 1+1 +, = 1+1 +, = (.10) 5

6 The Clebsch Gordon coefficients for these multiplication can be found in Ref. [8]. The fermions of all sectors (up quark, down quark, charged leptons) are assigned to + 1 represtations of Q 6. The model assumes three families of Higgs bosons, which are also assigned to +1 under Q 6. Their transformation properties are given by ( ) ( ) ψ1 ψ ψ = =, ψ c c = 1 =, ψ 3 = 1 ψ3 c = 1, (.11) ψ H = ( H1 H ψ c ) =, H 3 = 1. (.1) Here ψ generically denotes the fermion fields, and H denotes the up type and the down type Higgs fields which are doublets of SU() L. Due to the constraints of supersymmetry, H u and H3 u couple only to up quarks, while H d and H3 d couple to down type quarks and leptons. The Yukawa couplings of the model in the down quark sector arise from the superpotential W = α d ψ 3 ψ c 3H 3 +β d ψ T τ 1 ψ c 3H β dψ 3 ψ ct iτ H +δ d ψ T τ 1 ψ c H 3 +h.c. (.13) with similar results for up type quarks and charged leptons. This leads to the mass matrix for the down quarks given by M d = 0 δ d v d3 β d v d δ d v d3 0 β d v d1 β d v d β d v d1 α d v d3. (.14) Here v d1, v d, v d3 are the vacuum expectation values of H d 1,,3 fields, which break the Q 6 symmetry. Now, the potential of the Q 6 model admits an unbroken S symmetry which interchanges H u,d 1 H u,d. This unbroken symmetry allows us to choose a VEV pattern v u1 = v u, v d1 = v d,. (.15) Consequently, a 45 0 rotation of the matrix in Eq. (.14) in the 1- plane can be done both in the up and the down quark sectors without inducing CKM mixing. This will bring the mass matrices to the desired form of Eq. (1.4). By using the unbroken S symmetry, we make a 45 0 rotation on the Higgs fields, Ĥ 1, = (H 1 ±H )/, so that Ĥ1 acquires a VEV, while Ĥ = 0. We shall drop the hat on these redefined fields, and simply denote the VEV of the redefined H 1 as v 1. We assume that CP is a good symmetry of the Lagrangian, and that it is broken spontaneously by the VEVs of scalar fields. If the full theory contains SM singlet Higgs fields, spontaneous CP violation in the singlet sector will show up as soft CP violation in the Higgs doublet sector. Explicit examples of this sort have been given in Ref. [8]. For now we simply assume that the Yukawa couplings in Eq. (.13) are real, and the CKM CP violation has a spontaneous origin, via complex VEVs of the Higgs doublet fields. We denote the phase of these (redefined) VEVs as θ u = arg(v u3 ) arg(v u1 ), θ d = arg(v d3 ) arg(v d1 ). (.16) 6

7 We make an overall 45 rotation on the Q 6 doublets, Q, D c and U c, and then a phase rotations on these fields: U P u U, U c P u cu c (.17) and similarly for D and D c fields, where P u, d = 0 exp(i θ u, d ) exp(i θ u, d ) exp( i θ u, d ) 0 0 P u c, d c = exp( i θ u, d ),. (.18) This will make the originally complex mass matrices of Eq. (1.4) real, which we parametrize as M u, d = m 0 t, b 0 q u, d /y u, d 0 q u, d /y u, d 0 b u, d 0 b u, d yu, d. (.19) These real mass matrices can be diagonalized by the following orthogonal transformations: m u, d 0 0 Ou, T d M u, dmu, T d O u, d = 0 m c, s 0, O T u c, d cmt u, dm u c, d co u c, d c = 0 0 m t, b m u, d m c, s 0. (.0) 0 0 m t, b The CKM matrix V CKM is then given by where with θ q = θ d θ u. V CKM = O T up q O d, (.1) P q = P up d = e iθq e iθq (.) Now it is clear how the Q 6 setup reduces the number of parameters in the quark sector. The total number of parameters in the quark sector is nine (four real parameters each in M u and M d, plus a single phase θ q ), which should fit ten observables. Spontaneous CP violation is crucial for this reduction of parameters. With explicit CP violation, there would have been one more phase 7

8 parameter. The single prediction of this model was numerically studied in Ref. [8], and shown to be fully consistent with data. Here we present a numerical fit to all the quark sector observables, which deviates somewhat from the fit given in Ref. [8]. The difference arises since here we have attempted to be consistent with the recent lattice determination of light quark masses. An excellent fit to the quark masses and mixings, including CKM CP violation, is obtained with the following choice of parameters at a momentum scale of µ = 1 TeV. m 0 t = GeV, m0 b =.5515 GeV, θ q = θ d θ u = 1.40, q u = , b u = , b u = , y u = , q d = , b d = , b d = , y d = , (.3) This choice yields at µ = 1 TeV, the following masses and mixings for the quarks: m u = 1.13 MeV, m c = GeV, m t = GeV, m d =.53 MeV, m s = 50.99MeV, m b =.43 GeV, V CKM = , η W = , (.4) where η W is the CP violation parameter in the Wolfenstein parametrization. These values, when extrapolated to lower energy scales, give extremely good agreement with data [18]. We have computed the orthogonal matrices that diagonalize M u and M d. These rotation matrices will be relevant for our discussion of Higgs induced flavor violation, as well as FCNC arising via SUSY loop diagrams. We find O d = O u = , O d c =, O u c = ,. (.5) In the case of charged leptons, there is some arbitrariness in the values of (A, B, B, C) l of Eq. (1.4), since we have three observables (charged lepton masses) and four parameters (without including the neutrino sector). We shall present a fit with a simplifying assumption B l = B l. At µ = 1 TeV, a consistent fit for all the lepton masses is found with the following input values: A l = GeV, B l = B l = GeV, C l = GeV. (.6) 8

9 These yield the following eigenvalues at µ = 1 TeV: m e = MeV, m µ = MeV, m τ = MeV. (.7) These values correspond to the central values of charged lepton masses when extrapolated down to their respective mass scales [18]. The orthogonal matrix that diagonalizes M e is given by O e = , (.8) with O e c obtained from the above by flipping the signs in the first row and column. 3 Symmetry breaking and the Higgs boson spectrum We now turn to the discussion of symmetry breaking and the Higgs boson spectrum in the model. We shall confine here to the case of having three pairs of Higgs doublets, and no Higgs singlets in the low energy theory. It is however, assumed that singlet fields are present in the full theory, so that spontaneous Q 6 breaking in the singlet sector appears as soft breaking in the doublet sector. As shown in Ref. [8], it is possible to realize such a scenario while preserving the 1 interchange symmetry for members (1, ) inside Q 6 doublets. We seek a consistent picture where CP violating phases are generated in the Higgs doublet VEVs. As it turns out, CP also has to be softly broken in the bilinear soft SUSY breaking terms, or else there would be no CP phases in the VEVs. The superpotential that we consider is the most general one consistent with softly broken Q 6 symmetry, but preserving the S interchange symmetry: W eff = µ 1 (H u 1 Hd 1 +Hu Hd )+µ 3H u 3 Hd 3 +µ 13(H u 1 +Hu )Hd 3 +µ 31 H u 3(H d 1 +H d )+µ 1 (H u 1H d +H d 1H u ). (3.9) As mentioned earlier, we make a 45 rotations in H1 d, Hd and Hu 1, Hu space, with Ĥu 1, = Hu 1 ±Hu and Ĥd 1, = Hd 1 ±Hd, so that the superpotential becomes W eff = (µ 1 +µ 1 )Ĥu 1Ĥd 1 +(µ 1 µ 1 )Ĥu Ĥd +µ 3 Ĥ u 3Ĥd 3 + µ 13 Ĥ u 1Ĥd 3 + µ 31 Ĥ u 3Ĥd 1. (3.30) The redefined fields have Ĥu = Ĥd = 0. We work in the hatted basis from now on, and drop the hat on the new fields. The soft SUSY breaking Lagrangian is given, in the rotated basis, as V soft = (b 1 +b 1 )H1ǫH u 1 d +(b 1 b 1 )HǫH u d +b 3 H3ǫH u 3 d + b 13 H1 u ǫhd 3 + b 31 H3 u ǫhd 1 +h.c. +m d1 ( Hd 1 + H d )+m d3 Hd 3 +m u1 ( Hu 1 + H u )+m u3 Hu 3, (3.31) 9

10 where ǫ = iσ. The full scalar potential including the soft terms, the F terms and the D terms has the form V = M d1 ( Hd 10 + H d 1 )+M d3 ( Hd 30 + H d 3 ) +M u1( H u H u 1 + )+M u3( H u H u 3 + ) +{M13 (Hd 10 H d0 3 +H d 1 H d 3 )+M 31 (H3 u0 H1 u0 +H3 u+ H1 u+ )+h.c.} +M d ( Hd 0 + H d )+M u ( Hu 0 + H u + ) +{b 1(H u 1 + H d 1 H u 1 0 H d 10 )+b3 (H u 3 + H d 3 H u 3 0 H d 30 ) + b 13 (H u 1 + H d 3 H u 1 0 H d 30 )+ b31 (H u 3 + H d 1 H u 3 0 H d 10 ) +b (H u + H d H u 0 H d 0 )+h.c.} (g 1 +g )( Hu H1 u0 + H3 u+ + H3 u0 + H u+ + H u0 H1 d H10 d H3 d H30 d H d H d 0 ) + 1 g Hu + 1 H10 d +H u 0 1 H1 d +H u + 3 H30 d +H u 0 3 H3 d Here we have redefined new effective parameters for convenience as +H u + H d 0 +H u 0 H d. (3.3) M d1 = µ 1 +µ 1 + µ 31 +m d1, M d3 = µ 3 + µ 13 +m d3, M u1 = µ 1 +µ 1 + µ 13 +m u1, M u3 = µ 3 + µ 31 +m u3, M d = µ 1 µ 1 +m d1, M u = µ 1 µ 1 +m u1, M 13 = (µ 1 +µ 1 ) µ 13 + µ 3 µ 31, M 31 = (µ 1 +µ 1 )µ 31 + µ 3 µ 13, b 1 = b 1 +b 1, b = b 1 b 1. (3.33) Before analyzing the spectrum, let us note that the potential should be bounded from below along all D flat directions. The following conditions should be satisfied: M d1 +M u1 b 1 > 0, M d1 +M u > 0, M d1 +M u3 b 31 > 0 M d +M u1 > 0, M d +M u b > 0, M d +M u3 > 0, M d3 +M u1 b 13 > 0, M d3 +M u > 0, M d3 +M u3 b 3 > 0. (3.34) In our numerical analysis, we shall verify that these conditions are indeed met. We parametrize the VEVs of the four neutral Higgs fields as v u1 = vsinβsinγ u e iθ u1, v u3 = vsinβcosγ u e iθ u3, v d1 = vcosβsinγ d e iθ d1, v d3 = vcosβcosγ d e iθ d3. (3.35) 10

11 Thus we have v u1 + v u3 + v d1 + v d3 = v = (174 GeV). γ u(d) reflect the orientation of the VEVs in the H u(d)1 H u(d)3 space, while tanβ is analogous to the up/down VEV ratio of MSSM. We can rewrite the potential of the H 1 H 3 sector of the neutral Higgs fields which acquire VEVs in a compact form: V (1 3) N = ( H1 u0 H3 u0 ( +[ H1 u0 H3 u0 ) ( M u1 M 31 M31 M u3 )( H1 u0 H3 u0 ) ( b 1 b (g 1 +g )[ ( H u0 1 H u0 3 b 31 b 3 ) ( H u0 1 H u0 3 ) ( + )( ) ( H d0 1 H d0 3 H1 d0 H3 d0 ) +h.c. H d0 1 H d0 3 ] ) ( M d1 M 13 ) ( H d0 1 H d0 3 M13 M d3 )( H d0 1 H d0 3 ) )]. (3.36) This suggests a unitary transformation that would diagonalize the first two matrices in Eq. (3.36), while leaving the D term unaffected. With such a rotation we have ( ) ( )( ) V (1 3) m N = h 1 h 1 0 h ( ) ( )( ) 1 m + h 0 m 3 h 3 0 h 3 4 h 0 m 4 h 4 [ ( ) ( )( ) ] m 13 m 14 h 3 + h 1 h +h.c. m 3 m 4 h ( ) 8 (g 1 +g)[ ( ) h ( ) ( )] h 1 h 1 h h 3 h 3 4. (3.37) The unitary transformations to go from Eq. (3.36) to Eq. (3.37) are defined as ( ) ( ) ( )( )( h 1 H1 u cosω u sinω u e iφu 0 = U U = Q h H3 u u sinω u cosω u 0 e i(φu+θ M 31 ) ( ) ( ) ( )( )( h 3 H1 d cosω d sinω d e iφ d 0 = U D = Q d h 4 sinω d cosω d 0 e i(φ d+θ M13 ) H d 3 h h 4 H u 1 H u 3 H d 1 H d 3 ) ),, (3.38) with θ M31 = arg(m31 ), θ M 13 = arg(m13 ) and ω u = 1 ( ) M 31 tan 1, ω Mu3 d = 1 ( ) M 13 M u1 tan 1. (3.39) Md3 M d1 The two phases φ u and φ d here are arbitrary. φ u φ d does not appear in the potential (being proportional to U(1) Y charges). φ u + φ d can be used to remove one phase of the bilinear terms in the potential. Q u,d are arbitrary diagonal phase matrices. If desired, one can take advantage of these phases to remove all but one phase from the parameters of the potential. Since we are interested in going back to the original basis from this rotated basis, we find it convenient to set Q u,d to be identity. 11

12 The other parameters of this transformation are m 1, = 1 ] [M (M u3 +Mu1 ± u3 Mu1 ) +4 M 31, m 3,4 = 1 [ ] Md3 +Md1 ± (M d3 Md1 ) +4 M 13. (3.40) and ( ) ( ) m 13 m 14 m 3 m 4 = U U b 1 b 13 b 31 b 3 U D. (3.41) If we choose [ φ u +φ d = π +arg b 1 sinω usinω d + b 31 cosω u sinω d e iθ M 31 + b 13 sinω u cosω d e iθ M 13 +b3 cosω u cosω d e i(θ M 31 +θ M13 ) ], (3.4) m 4 is real and positive (with Q u,d set to identity). We shall adopt this phase convention in our numerical study. However, we shall present analytical results that hold in an arbitrary phase convention. The task at hand is somewhat simplified, since Eq. (3.37) is relatively simple to analyze. The eight real neutral Higgs bosons in H u,d 1,3 can be conveniently parametrized as [ h 1 = e iδ 1 v (φ 1 +ieφ 5 +iaφ 7 +i v ] 1 v G), [ h = e iδ v + 1 (φ +ifφ 6 +ibφ 7 +i v ] v G), [ h 3 = e iδ 3 v (φ 3 +igφ 5 +icφ 7 i v ] 3 v G), h 4 = v (φ 4 +ihφ 6 +idφ 7 i v 4 G). (3.43) v Here v i (i = 1,, 3, 4) are the magnitudes of the VEVs of the redefined fields h i, and δ i are their phases. Without loss of generality we have taken v 4 to be real. G in Eq. (3.43) is the Goldstone field eaten up by the Z gauge boson. We shall work in the unitary gauge and set G = 0. We have checked explicitly that the G field does not mix with other scalar fields, and that its mass is exactly zero. The coefficients of various fields in Eq. (3.43) are functions of the v i s: a = v 1 v +v4 v, b = v v 1 +v3 v1 +v3 v, c = v 3 v +v4 v +v4 v, d = v 4 v 1 +v3 v1 +v3 v, v +v4 v 3 v e = 4 v, f = 1 v, g =, h =. (3.44) v 1 +v3 v +v4 v 1 +v3 v +v4 1

13 We shall allow for the soft SUSY breaking parameters (b i ) in the Higgs potential to be complex. Phase rotations cannot remove all phases from the potential, one phase is unremovable. Without this phase, the model cannot induce complex VEVs to the doublets, as shown in Ref. [19] by a geometric argument. For the case when all parameters in the Higgs potential are real, we have numerically verified that the CP violating extremum would generate two massless modes, signalling inconsistency with symmetry breaking [19]. We take the soft bilinear terms m 13, m 14, m 3, m 4 of Eq. (3.37) to be complex, and denote the phase of m ij as θ ij. The minimization conditions then read as m 1 v 1 + m 13 v 3cos(θ 13 +δ 1 +δ 3 )+ m 14 v 4cos(θ 14 +δ 1 )+ 1 4 (g 1 +g )v 1(v 1 +v v 3 v 4 ) = 0, m v + m 3 v 3cos(θ 3 +δ +δ 3 )+ m 4 v 4cos(δ +θ 4 )+ 1 4 (g 1 +g )v (v 1 +v v 3 v 4 ) = 0, m 3 v 3 + m 13 v 1cos(θ 13 +δ 1 +δ 3 )+ m 3 v cos(θ 3 +δ +δ 3 ) 1 4 (g 1 +g )v 3(v 1 +v v 3 v 4 ) = 0, m 4 v 4 + m 14 v 1cos(θ 14 +δ 1 )+ m 4 v cos(δ +θ 4 ) 1 4 (g 1 +g )v 4(v 1 +v v 3 v 4 ) = 0, m 13 (v 1 +v 3)sin(θ 13 +δ 1 +δ 3 )+ m 3 v 1 v sin(θ 3 +δ +δ 3 )+ m 14 v 3 v 4 sin(θ 14 +δ 1 ) = 0, m 4 (v +v 4 )sin(θ 4 +δ )+ m 14 v 1v sin(θ 14 +δ 1 )+ m 3 v 3v 4 sin(θ 3 +δ +δ 3 ) = 0, m 14 v 1 v 4 sin(θ 14 +δ 1 ) m 3 v v 3 sin(θ 3 +δ +δ 3 ) = 0. (3.45) Denoting the squared matrix for φ i, i = 1,,...7 from the H 1 H 3 sector as we obtain M 11 = λv 1 +κv v 4 [cot(θ v1 14 +δ 1 ) cot(θ 13 +δ 1 +δ 3 )], M = λv +κv 4 [cot(θ 3 +δ +δ 3 ) cot(θ 4 +δ )], v M 33 = λv 3 +κv v 4 [cot(θ v3 3 +δ +δ 3 ) cot(θ 13 +δ 1 +δ 3 )], M 44 = λv 4 +κv [cot(θ 14 +δ 1 ) cot(θ 4 +δ )], v 4 M 0,(1 3) = M ij, (3.46) M 55 = κ v v 4 [ v 3 cot(θ v1 +v 3 v1 14 +δ 1 )+ v 1 cot(θ v3 3 +δ +δ 3 ) (v 1 +v3) cot(θ v1 13 +δ 1 +δ 3 )], v 3 M 1 66 = κ v v 4 (v +v4) [v4 cot(θ 14 +δ 1 )+v4cot(θ 4 3 +δ +δ 3 ) (v +v4) cot(θ 4 +δ )], M 77 = κv v 4 (v1 +v +v 3 +v 4 ) [cot(θ (v1 +v3)(v +v4) 14 +δ 1 )+cot(θ 3 +δ +δ 3 )], M 1 = λv 1v, M 13 = λv 1v 3 +κ v v 4 v 1 v 3 cot(θ 13 +δ 1 +δ 3 ), 13

14 M 14 = λv 1v 4 κ v cot(θ 14 +δ 1 ), M 15 v = κ v v 4, M 1 v 3 v 1 +v3 16 = κ v, v 1 v +v4 M 17 = κ v v 4 v 1 +v +v3 +v4, M v 1 v 1 +v3 v +v4 3 = λv v 3 κ v 4 cot(θ 3 +δ +δ 3 ), v 3 M 4 = λv v 4 +κcot(θ 4 +δ ), M v 1 v 4 5 = κ, M v v 3 v 1 +v3 6 = κ, v +v4 M 7 = κ v 4 v 1 +v +v3 +v4, M v 1 +v3 v +v4 34 = λv 3 v 4, M v v 4 35 = κ, v 1 v 1 +v3 M 36 = κ v4, M v 3 v +v4 37 = κv v 4 v 1 +v +v 3 +v 4, M v 3 v 1 +v3 v +v4 45 = κ v v 3, v 1 v 1 +v3 M 46 = κ v 4, M v +v4 47 = κv v 1 +v +v 3 +v 4, v 1 +v3 v +v4 M 1 56 = κ [v v 1 v 3 v 1 +v3 v v 3cot(θ 14 +δ 1 )+v1v 4cot(θ 3 +δ +δ 3 )], +v4 M 57 = κ v v 4 v 1 +v +v3 +v4 v 1 v 3 (v1 +v 3 ) [v3cot(θ 14 +δ 1 ) v1cot(θ 3 +δ +δ 3 )], v +v 4 v M 67 = κ 1 +v +v 3 +v 4 v 1 +v3(v +v4) [v cot(θ 14 +δ 1 ) v4 cot(θ 3 +δ +δ 3 )]. (3.47) Here we have defined λ = (g 1 +g )/ = M Z /v and κ = m 4 sin(θ 4 +δ ). The potential of the H u H d fields which do not acquire VEVs is V () N = M u H u +M d H d {b H u H d +h.c.} + g 1 +g ( H u 8 H d + v u1 + v u3 v d1 v d3 ). (3.48) The corresponding squared mass matrix for the scalars in the basis (ReH u, ImH u, ReH d, ImH d ) is M 0() = Mu m Z cosβ 0 Reb Imb 0 Mu m Z cosβ Imb Reb Reb Imb Md + m Z cosβ 0 Imb Reb 0 Md + m Z cosβ. (3.49) This matrix has two pairs of degenerate eigenstates, owing to an unbroken U(1) symmetry. TheH u,d 1 H u,d 3 sectorchargedhiggsbosonmassmatrixis,inthebasis{h1 u+, H3 u+, H1 d, H d 3 }, with M ±(1 3) = (M ) ij, M 11 = M u1 1 m Z cosβ + 1 g v d1, M = M u3 1 m Z cosβ + 1 g v d3, 14

15 M 33 = Md1 + 1 m Z cosβ + 1 g v u1, M 44 = Md3 + 1 m Z cosβ + 1 g v u3, M 1 = M 1 = M g vd1v d3, M 13 = M 31 = b g vu1v d1, M 14 = M41 = b g v u3 v d1, M 3 = M3 = b g v u1 v d3, M 4 = M4 = b g v u3 v d3, M 34 = M43 = M g v u1vu3. (3.50) Finally, the H u Hd sector charged Higgs mass matrix is, in the basis {Hu +, H d }, ( ) M ±() = Mu 1 m Z cosβ b Md + 1 m Z cosβ b (3.51) Now we present two sets of numerical fits (cases (1) and ()) which show the consistency of symmetry breaking. We are interested in choosing the SUSY breaking parameters (including the µ terms) around the TeV scale, guided by arguments of naturalness. At the same time we wish the spectrum to be consistent with FCNC constraints arising from meson antimeson mixings. We have explored parameter space of the Higgs potential where both these constraints are met. For the FCNC constraint, we allow the new Higgs exchange contribution to M be not more than the experimentally measured values. Case (1) The parameters in the original Higgs potential of Eq. (3.36) are taken to have the following values. M d1 = TeV, M d3 = TeV, M u1 = 4.78 TeV, M u3 =.15 TeV, M 31 =.336 e i0.79 TeV, M 13 = e i1.05 TeV, b 1 = ei.963 TeV, b 3 = e i.064 TeV, b 31 = 4.05 e i.186 TeV, b 13 = e i3.109 TeV, M u = TeV, M d = TeV b = TeV, (3.5) In the representation of Eq. (3.37) this choice corresponds to m 1 = TeV, m = TeV, m 3 = 3.93 TeV,m 4 = TeV, m 13 = 1.851e i1.437 TeV, m 14 =.736e i0.73 TeV, m 3 =.44 e i1.347 TeV, m 4 =.104 TeV. (3.53) For completeness we also give values of other parameters, ω u = 0.70, ω d = 0.6, φ u +φ d = We obtain numerically the VEV parameters to be tanβ =.00, θ d = 0.03, θ u = 1.37, tanγ d =.50, tanγ u = (3.54) 15

16 The mass eigenvalues of the Higgs bosons in the H 1 H 3 sector are found to be M h0 = (99.4, 115.1) GeV, M 1 = 3.99 TeV, M M 1 = 0.6 GeV, M 3 = TeV, M 4 M 3 = GeV, M 5 = 5.14 TeV, M 6 M 5 = GeV. (3.55) Note the appearance of nearly degenerate states (M 1, M ) etc, with their mass splitting being proportional to m Z /4. The Higgs bosons from the H sector have degenerate masses given by M 7 = M 8 = TeV M 9 = M 10 = TeV. (3.56) The charged Higgs bosons are nearly degenerate with its neutral partner, so we list the mass splittings: M ±1 M 1 = 0.53 GeV, M ± M 3 = GeV, M ±3 M 5 = 0.03 GeV. (3.57) In the (H u Hd ) sector, the two charged Higgs bosons are degenerate with the neutral ones given in Eq. (3.56). The mass eigenstates H i are mixtures of h i, i = 1,,... 7 states in the (1-3) sector. The orthogonal transformation that diagonalizes the mass matrix of Eq. (3.46) is H k = h 1 h h 3 h 4 h 5 h 6 h 7, (3.58) with k = 0, 6. Since b = 0 in this case, the H0 mass matrix is diagonal, and thus the mass eigenstates are the original state. Case () Here we take the input parameters corresponding to Eq. (3.36) to be M d1 = TeV, M d3 = 5.41 TeV, M u1 =.765 TeV, M u3 = 3.69 TeV, M 31 =.85 e i0.781 TeV, M 13 = e i0.949 TeV, b 1 = ei1.495 TeV, b 3 = e i1.5 TeV, b 31 = 7.40 e i.48 TeV, b 13 = e i.77 TeV, M u = TeV, M d = TeV, b = 1.34e i1.56 TeV. (3.59) 16

17 This choice corresponds to parameters in Eq. (3.37) to be m 1 = TeV, m = TeV, m 3 = TeV,m 4 = TeV, m 13 = 3.81e i1.71 TeV, m 14 = 1.70e i0.974 TeV, m 3 = e i0.501 TeV, m 4 =.36 TeV, (3.60) with ω u = 0.501, ω d = 0.606, φ u +φ d = The Higgs VEV parameters are found for this input to be tanβ =.40, θ d = 0.06, θ u = 1.34, tanγ d = 1.80, tanγ u = (3.61) The mass spectrum of Higgs boson in the H 1 H 3 sector is M h0 = (104.1, 119.) GeV, M 1 =.869 TeV, M M 1 = 0.35 GeV M 3 = TeV, M 4 M 3 = 0.13 GeV, M 5 = TeV, M 6 M 5 = GeV, (3.6) while the mass eigenvalues of Eq. (3.49) are M 7 = M 8 = TeV M 9 = M 10 = TeV. (3.63) The charged Higgs boson masses are given by M ±1 M 1 = 0.5 GeV, M ± M 3 = 0.18 GeV, M ±3 M 5 = GeV. (3.64) with the remaining two charged Higgs bosons being degenerate with the neutral ones given in Eq. (3.63). The orthogonal matrix that diagonalizes Eq. (3.46) is h h h 3 H k = h 4, (3.65) h h h 7 The matrix diagonalizing Eq. (3.49) is H k = Re(H u) Im(H u) Re(H d) Im(H) d, (3.66) 17

18 with k = 7, 10. In these fits, M h0 is the light standard model like Higgs boson mass, for which radiative corrections are significant. In our computation we have included known two loop corrections. The two values listed for M h0 correspond to zero and maximal left right stop mixing (X t = 0 or 6). We have taken m t = 174 GeV, M SUSY = 1.5 TeV and α s (m t ) = for these evaluations and used the analytic approximation given in Ref. [0]. An interesting feature of these two fits is that the diagonal entries of the quadratic mass matrix of the potential of Eq. (3.36) are all positive. This of course does not preclude some soft squared masses turning negative as in the MSSM via large top quark Yukawa coupling (since the diagonal entries also receive µ term contributions), however, this is not necessary for symmetry breaking to be triggered. Yet, one of the eigenvalues of this matrix is negative, which facilitates symmetry breaking. For the two cases we find these eigenvalues to be Case (1) : {(5.13 TeV), (4.161 TeV), (3.300 TeV), (38.68 GeV) }, Case () : {(5.658 TeV), (5.115 TeV), (.869 TeV), (45.40 GeV) }. (3.67) The conditions for boundedness of the potential listed in Eq. (3.34) are found to be satisfied for both cases. 3.1 Neutralino and Chargino masses The symmetry breaking parameters do not fully determine the masses of the neutralinos and the charginos. Here we present analytical results for their mass matrices. is The mass matrix of H 1 H 3 sector neutralino in the basis of { B, W 0, Hu 1 0, Hu 3 0, Hd 1 0, Hd 3 0 } M χ 0 (13) = M B 0 g v u1 0 M W g 1v u1 g v u1 g v u3 g v d1 g v d3 g v u3 g 1v u3 g v d1 g 1 v d1 g v d3 g 1 v d3 g 1v u1 0 0 (µ 1 +µ 1 ) µ 13 g 1v u3 0 0 µ 31 µ 3 g 1 v d1 (µ 1 +µ 1 ) µ g 1 v d3 µ 13 µ (3.68) The mass matrix of the H sector in the basis { H u 0 0, Hd } is ( 0 (µ 1 µ 1 ) M χ 0 () (µ 1 µ 1 ) 0 ). (3.69) Themassmatrixofcharginosofthe H 1 H 3 sectorinthebasis{ W +, Hu 1 +, Hu 3 +, W, H d 1, Hd 3 } 18

19 has a block diagonal form: M χ ± (13) = ( 0 X T X 0 ) (3.70) with X = M W g 1 v u1 g 1 v u3 g 1 v d1 µ 1 +µ 1 µ31. (3.71) g 1 v d3 µ13 µ 3 The chargino mass matrix in the H u H d sector in the basis of { H u +, Hd } is ( ) 0 µ 1 µ 1 M χ ± () = µ 1 µ 1 0 (3.7) 4 Tree level Higgs induced FCNC processes In this section we discuss various FCNC processes mediated by tree level neutral Higgs boson exchange. 4.1 Neutral meson mixing via Higgs exchange Accurate measurements exist [1] for neutral meson antimeson mixings in the K 0 K 0, B 0 d B0 d, B 0 s B 0 s and in D 0 D 0 sectors. In the Q 6 model there are new contributions to these mixings arising through tree level Higgs exchange. These new contributions will modify the SM predictions, which are all in good agreement with data. Here we compute these new contributions, following the analysis of Ref. [], with updated QCD corrections and hadronic matrix elements. The Yukawa coupling α u, d, β u, d, β u, d, δ u, d of Eq. (.14) can be determined from the mass matrix Eq. (.19): α u, d = m0 t, b y u, d, β u, d = m0 t, b b u, d v u, d3 v u, d1 β u, d = m0 t, b b u, d v u, d1 Using the input values given in Eq. (.3) we get for the two cases Case (1) Case (), δ y, d = m0 t, b q u, d/y u, d, (4.73) v u, d3 α d = , β d = , β d = , δ d = , α u = , β u = , β u = 0.167, δ u =

20 are α d = 0.056, β d = , β d = , δ d = , α u = , β u = 0.037, β u = 0.074, δ u = After 45 rotation in the Q 6 doublet space, the Yukawa coupling matrices in the down sector β d Y d1 = Od T P d 0 0 β d P d co d c, Y d = Od T P d β d 0 β d δ d 0 Y d3 = Od T P d δ d α d P d co d c, P d co d c, (4.74) where P d, P d c are defined in Eq. (.18). The Yukawa couplings in the up quark sector and the charged lepton sector are similar. The new Higgs mediated contributions to F = Hamiltonian, responsible for the neutral meson antimeson mixings has the form [] H eff = 1 [ ( q i M k Y k ij 1+γ 5 +Y k ji ] ) 1 γ 5 q j. (4.75) Here q i,j are the relevant quark fields contained in the meson. Y k ij are the Yukawa couplings of q i, q j with Higgs mass eigenstate H k mediating FCNC interactions, k = 1,,...10 in our model, 6 from the (H 1 H 3 ) sector and 4 from the H sector. (The light standard model like Higgs boson has practically no FCNC couplings.) Y k ij can be obtained via inverse transformations, Eq. (3.38), (3.43) and (3.58) or (3.65). We obtain M φ 1 = φ H eff φ = f φ m [ φ M k 5 m ( ) B φ Y k 4(m qi +m qj ) ij +Y k ji η (µ) ( 1 +Yij k Y ji k ) m φ (m qi +m qj ) ] B 4 η 4 (µ). (4.76) Here φ is the neutral meson (K 0, B 0 d, B0 s, D 0 ). For our numerical study we use the modified vacuum saturation and factorization approximation results for the matrix elements [, 3] ( ) 1 φ f i (1±γ 5 )f j fi (1 γ 5 )f j φ = f φ m φ 6 + m φ B (m qi +m qj ) 4, φ f i (1±γ 5 )f j fi (1±γ 5 )f j φ = 5 6 f φ m φ m φ (m qi +m qj ) B. (4.77) B and B 4 are equal to one in the vacuum saturation approximation, but are found to be slightly different from one in lattice simulations. We use (B, B 4 ) = (0.66, 1.03) for the K 0 system, 0

21 (0.8, 1.16) for the Bd 0 and B0 s systems, and (0.8, 1.08) for the D 0 system []. In Eq. (4.76) η (µ), η 4 (µ) are QCD correction factors of the Wilson coefficients C and C 4 of the effective F = Hamiltonian in going from the SUSY scale M s to the hadronic scale µ. These factors are computed as follows. The F = effective Hamiltonian has the general form where H F= eff = 5 C i Q i + i=1 3 i=1 C i Q, (4.78) Q 1 = q i α L γ µ q j α L q i β L γν q j β L, Q = q i α R q j α L q i β R q j β R, Q 3 = q i α R q j β L q i β R q j α L, Q 4 = q i α R q j α L q i β L q j β R, with Q 1,,3 obtained from Q 1,,3 by the interchange L R. Q 5 = q i α R q j β L q i β L q j α R, (4.79) For computing η,4 we take the SUSY scale M s to be 1 TeV. All the supersymmetric particles and heavy Higgs bosons are integrated out at 1 TeV. The Wilson coefficients evolve from M s down to the hadron scale µ according to the equations C r (µ) = i s (b (r,s) i +ηc (r,s) i )η a i C s (M s ), (4.80) Here η is defined as η = α s (M s )/α s (m t ). The magic numbers a i, b (r,s) i and c (r,s) i can be found in Ref. [] for the K system, in Ref. [3] for the B d,s system and in Ref. [4] for the D system. With M s = 1 TeV and α s (m Z ) = 0.118, and m t (m t ) = GeV we find η = α s (1 TeV)/α s (m t ) = 0.088/0.108 = At the SUSY scale, the neutral Higgs bosons in our model generate only operators Q and Q 4. Consequently, at the hadron scale, for the K 0 system, we find C (µ) = C (M s ) (.54), C 4 (µ) = C 4 (M s ) (4.81), C 3 (µ) = C (M s ) ( ), C 5 (µ) = C 4 (M s ) (0.186), (4.81) leading to η (µ) =.54, η 4 (µ) = Although operator mixings induce non-zero C 3 and C 5 at the hadronic scale, their coefficients are found to be rather small. For the Bd,s 0 system, following the same procedure, we find C (µ) = C (M s ) (.00), C 4 (µ) = C 4 (M s ) (3.1), C 3 (µ) = C (M s ) ( ), C 5 (µ) = C 4 (M s ) (0.0874). (4.8) And for the D 0 system we have C (µ) = C (M s ) (.31), C 4 (µ) = C 4 (M s ) (3.99), C 3 (µ) = C (M s ) ( ), C 5 (µ) = C 4 (M s ) (0.144). (4.83) 1

22 In all cases we see that the induced operators C 3 and C 5 are negligible. K 0 K 0 mixing constraint: In the K 0 system, tree level neutral Higgs boson exchange contributes to K L K S mass difference, as well as to the indirect CP violation parameter, modifying the successful SM predictions. The mass difference is computed from m K = ReM1 K, while the CP violation parameter is ǫ K ImMK 1 mk. We seek consistency with the precisely measured experimental values m K /m K (7.1 ± 0.014) and ǫ K In our calculation, we choose m K = 498 MeV and f K = 160 MeV. For the two numerical fits we find the new contributions to be Case (1) : ( m K /m K ) new = , ǫ new K = , Case () : ( m K /m K ) new = , ǫ new K = (4.84) The contributions from H 0 1 H 0 3 sector and H 0 sector to Re(M K 1) are respectively ( , ) GeV for case (1) and ( , ) GeV for case (). We see that the new contributions to the mass difference is significant, but consistent with data. New contributions to CP violation is suppressed, which is a generic feature of Higgs exchange in this class of models. We elaborate on this issue later in this section. B 0 d B0 d mixing constraint: FortheB 0 d B 0 d system Weuseasinputm B d = 5.81 GeV, f Bd = 40 MeV andseekconsistency with the experimental value m Bd = GeV. We find for the Higgs induced contribution Case (1) : ( m Bd ) new = GeV, Case () : ( m Bd ) new = GeV. (4.85) The contributions from H 0 1 H0 3 sector and H0 sector to Mb d 1 are ( , ) GeV for case (1) and ( , ) GeV. Again, we see consistency with experimental values. CP violation parameter is found to be extremely tiny, 10 5, from Higgs boson exchange. B 0 s B 0 s mixing constraint: For the B 0 s B0 s system, we use m B s = 5.37 GeV, f Bs = 95 MeV and compare the new contributions with m Bs = GeV. Case (1) : ( m Bs ) new = GeV, Case () : ( m Bs ) new = GeV, (4.86) The H 0 1 H0 3 sector and the H0 sector contribute to MB0 s 1 given by ( , ) GeV for case (1) and ( , ) GeV for case (). These new contributions

23 are within experimentally allowed range. Higgs mediated CP violation is again found to be highly suppressed. D 0 D 0 mixing constraint: For the D 0 D 0 mixing we use m D = GeV, f D = 00 MeV and compare the new contribution with m D = GeV. Case (1) : ( m D ) new = GeV, Case () : ( m D ) new = GeV, (4.87) The H1 0 H0 3 sector contribution has different sign from that of the H0 sector. We find for MD0 1 these contributions to be ( , ) GeV for case (1) and ( , ) GeV for case (). Again these limits are within experimental range. We have found that new sources of CP violation through tree level Higgs is very small in meson antimeson mixings with typical values Im(M 1 ) 10 4 Re(M 1 ). This can be understood heuristically as follows. There are two types of contributions to the meson mixing as given in Eq. (4.76). The first term, proportional to B respects a global U(1) symmetry (strangeness in the K 0 system), which is only broken by the mass splittings in the neutral Higgs boson spectrum between a pair of particles. However, this splitting is very small, of order m Z in the squared mass, see Eqs. (3.55). The couplings of the nearly degenerate Higgs in each pair differ by a factor i, owing to the U(1) symmetry, and the two contributions cancel, in the limit of exact degeneracy. For both the real and imaginary parts of M 1 the contributions from the first term is suppressed by a factor m Z /(4M k ). Such a suppression is absent in the second term of Eq. (4.76), since the operator Q 4 explicitly breaks the U(1) symmetry. Thus, although the first term has CP violation, in relation to the CP conserving second term, it is suppressed by a factor m Z /(4M k ) Now, the second term, while it has no suppression factor, it is purely real. This can be seen from the following observation. In the mass basis of fermions in the original basis we have the relation (owing to the vanishing of off diagonal mass terms in the mass eigenbasis) (Y d3 ) ij H d 3 = (Yd1 ) ij H d 1 (4.88) for i j. The couplings of mass eigenstates of Higgs boosn to down type quarks are simply linear combinations of H1 d and Hd 3. Since we assume CP to be spontaneously broken, all components of (Y k ) ij with i j have the same phase. As a result the second term of Eq becomes real. The constraint imposed by SUSY, that Hu fields do not couple to down type quarks, and the fact that only two of the down type Higgs bosons acquire VEVs is very crucial for this result. 3

24 4. Neutron electric diploe moment from Higgs exchange Higgs boson exchange can generate non-zero electric dipole moments for the fermions. These diagrams are however suppressed by the light fermion Yukawa couplings. For the d quark EDM arising from neutral Higgs boson exchange at the one loop level we find [] d d = Q de 16π Im(Y dqy k qd) k m [ ( )] q 3 M Mk ln k ξ m d, (4.89) q where ξ d = (α s (M k )/α s (µ)) 16/3 0.1, and q is summed over d, s and b. The neutron EDM is determined using the quark model via We find D n = 4d d /3 d u /3. (4.90) Case (1) : D n = e cm, Case () : D n = e cm, (4.91) which are well within experimental limits. The EDM of the electron is similarly found to be extremely small from the Higgs boson exchange diagrams. 4.3 µ 3e and τ 3µ decays Tree level Higgs boson exchange can lead to flavor violating leptonic decays such as τ 3µ and µ 3e. The effective weak interaction mediating such decays can be parametrized as G eff = (Y e ) k 11 (Y e) k 1 1 k M k (4.9) The effective couplings are found for µ 3e for the two cases to be Case (1) : G eff = G F, Case () : G eff = G F. (4.93) And the couplings for τ 3µ decay are Case (1) : G eff = G F, Case () : G eff = G F. (4.94) Such small effective couplings will lead to negligible contributions to the decay branching ratios. For example, the branching ratio for τ 3µ is of order 10 15, well below the experimental sensitivity. We conclude that Higgs mediated FCNC in the lepton sector are all safe. 4

25 5 FCNC mediated by SUSY particles In this section we turn attention to the flavor changing processes mediated by the supersymmetric particles. The main motivation for the non Abelian Q 6 model was to bring such processes under control by a symmetry reason. Here we analyze meson antimeson mixings, flavor violating leptonic decays, and the EDM of the neutron and the electron. We present our proposal to suppress SUSY contributions to the EDM by making the Higgsinos of the model light, with masses of order 100 GeV. Owing to the Q 6 symmetry, the first two family squarks (and similarly sleptons) are degenerate in mass, while the third family, which is a Q 6 singlet has a different mass. In the fermion sector Q 6 symmetry is broken, which means that there will be SUSY loop induced flavor violation in the model. Constraints on such flavor violation has been listed in Ref. [ 4] assuming all three families of squarks are degenerate. While these results are applicable for the K 0 and D 0 system in our model, they do not work well for the B 0 d,s system. This is because the masses of the b and d,s masses are not the same. 5.1 Generalized constraints for B d system We have generalized the results of Ref. [3] by allowing for b mass to be different from the masses of d,s. We define new parameters y d A,B = ( m b ) A,B m d A,B (5.95) for A,B = L,R. We expect these y parameters to be of order one, but not very close to one. Taking account of y 1 we have generalized the constraints on the squark mixing parameters from B 0 d system as follows. The effective F = Hamiltonian for B d,s system can be written as H eff = 5 C i Q i + i=1 3 i=1 C i Q = α s {(δ 16m 13) d LL(4Q d 1 xf 6 (x,y)+66q 1 f6 (x,y))+(δ13) d RR(4 Q 1 xf 6 (x,y)+66 Q 1 f6 (x,y)) +(δ13 d ) LL(δ13 d ) RR(504Q 4 xf 6 (x,y) 7Q 4 f6 (x,y)+4q 5 xf 6 (x,y)+10q 5 f6 (x,y)) +(δ d 13) RL(04Q xf 6 (x,y) 36Q 3 xf 6 (x,y))+(δ d 13) LR(04 Q xf 6 (x,y) 36 Q 3 xf 6 (x,y)) +(δ d 13 ) LR(δ d 13 ) RL ( 13Q 4 f(x,y) 180Q 5 f6 (x,y))}, (5.96) The functions f 6 (x,y) and f 6 (x,y) are f 6 (x,y) = 1 [ lnx(x+y +xy 3x 3 )(y 1) 3 (x 1) 3 (y 1) 3 (z 1) 3 5