Arxiv:1610.07314V2 [Hep-Ph]

Home , Chargino

arXiv:1610.07314v2 [hep-ph] 23 Dec 2016 ewrs Pvoaig lcrcdpl oet oa ag symmetr gauge local moment, dipole electric violating, CP Keywords: lar reach can numbers: EDMs PACS t and c the for assu results some numerical with The analyzed space. the are account results into take numerical We The t. and diagrams. c quarks the Consi of (EDM). EDMs moment the analyze dipole electric s quark violating the CP to new tributions are there symmetries(BLMSSM), lepton and ∗ † h n opcnrbtost ()eeti ioemmn in moment dipole electric c(t) to contributions loop one The [email protected] [email protected] nteC iltn uesmercetnino h standar the of extension supersymmetric violating CP the In 1 eateto hsc,HbiUiest,Boig010,C 071002, Baoding University, Hebei Physics, of Department 13.r 12.60.Jv,14.80.Cp 11.30.Er, ee nvriyo esine,Siizun 501 Chi 050031, Shijiazhuang Geosciences, of University Hebei h-i Zhao Shu-Min a-i Zhang Hai-Bin 2 colo ahmtc n Science, and Mathematics of School PvoaigBLMSSM violating CP 1 Dtd ac ,2018) 8, March (Dated: ∗ a-uFeng Tai-Fu , 1 igXn Dong Xing-Xing , Abstract 1 1 † hn-u Yang Zhong-Jun , 1 Tao-Guo , eigteC iltn hss we phases, violating CP the dering pin nterlvn parameter relevant the on mptions oe ihlclgue baryon gauged local with model d y evalues. ge ucswihcngv e con- new give can which ources otiuin rmteoeloop one the from contributions 2 1 , na hina the I. INTRODUCTION

The CP violation found in the K- and B system [1] can be well explained in the standard model. It is well known that, the electric dipole moment(EDM) of elementary particle is a clear sinal of CP violation[2]. The Cabbibo-Kobayashi-Maskawa(CKM) phase is the only source of CP violation in the SM, which has ignorable effect on the EDM of elementary particle. In the SM, even to two loop order, the EDM of a fermion does not appear, and there are partial cancelation between the three loop contributions[3]. If EDM of an elementary fermion is detected, one can confirm there are new CP phases and physics beyond the SM. Though SM has obtained large successes with the detection of the lightest CP-even Higgs h0[4], it is unable to explain some phenomena. Physicists consider the SM should be a low energy effective theory of a large model. The minimal supersymmetric extension of the standard model(MSSM) is very favorite and people have been interested in it for a long time[5]. There are also many new models beyond the SM, such as µνSSM[6]. Generally speaking, the new models introduce new CP-violating phases that can affect the EDMs of fermions, B0 B¯0 mixing et al. The EDMs of electron and neutron are strict constraints on − the CP-violating phases[7]. In the models beyond SM, there are new CP-violating phases which can give large contributions to electron and neutron EDMs[8]. To make the MSSM predictions of electron and neutron EDMs under the experiment upper bounds, there are three possibilities[9]: 1 the CP-violating phases are very small, 2 varies contributions cancel with each other in some special parameter spaces, 3 the supersymmetry particles are very heavy at several TeV order. Taking into account the local gauged B and L, people obtain the minimal supersymmetric extension of the SM, which is the so called BLMSSM[11]. The authors in Ref.[10] first proposed BLMSSM, where they studied some phenomena. At TeV scale, the local gauge symmetries of BLMSSM breaks spontaneously. Therefore, in BLMSSM R-parity is violated and the asymmetry of matter-antimatter in the universe can be explained. We have studied the lightest CP-even Higgs mass and the decays h0 VV, V = (γ,Z,W )[12] in → the BLMSSM, where some other processes[13] are also researched. Taking the CP-violating phases with nonzero values, the neutron EDM, lepton EDM and B0 B¯0 mixing are re- − searched in this model[14]. From neutron experimental data, the bounding of top EDM is analyzed[15]. Taking into

2 account the precise measurements of the electron and neutron EDMs, the upper limits of heavy quark EDMs are also discussed[16]. The upper limits on the EDMs of heavy quarks + are researched from e e− annihilation[17]. In the CP-violating MSSM, the authors study c quark EDM including two loop gluino contributions[18]. There are also other works on the c quark EDM[19]. Considering the pre-existing works, the upper bounds of EDMs for 17 15 c and t are about d < 5.0 10− e.cm and d < 3.06 10− e.cm. In this work, we c × t × calculate the EDMs of charm quark and top quark in the framework of the CP-violating BLMSSM. At low energy scale, the quark chromoelectric dipole moment(CEDM) can give important contributions to the quark EDM. So, we also study the quark CEDM with the renormalization group equations. In Section 2, we brieﬂy introduce the BLMSSM and show the needed mass matrices and couplings, after this introduction. The EDMs(CEDMs) of c and t are researched in Section 3. In Section 4, we give out the input parameters and calculate the numerical results. The last Section is used to discuss the results and the allowed parameter space.

II. THE BLMSSM

Considering the local gauge symmetries of B(L) and enlarging the local gauge group of the SM to SU(3) SU(2) U(1) U(1) U(1) one can obtain the BLMSSM C ⊗ L ⊗ Y ⊗ B ⊗ L model[11]. In the BLMSSM, there are the exotic superfields including the new quarks Qˆ (3, 2, 1/6, B , 0), Uˆ c (3¯, 1, 2/3, B , 0), Dˆ c (3¯, 1, 1/3, B , 0), Qˆc 4 ∼ 4 4 ∼ − − 4 4 ∼ − 4 5 ∼ (3¯, 2, 1/6, (1+ B ), 0), Uˆ (3, 1, 2/3, 1+ B , 0), Dˆ (3, 1, 1/3, 1+ B , 0), and − − 4 5 ∼ 4 5 ∼ − 4 the new leptons Lˆ (1, 2, 1/2, 0, L ), Eˆc (1, 1, 1, 0, L ), Nˆ c (1, 1, 0, 0, L ), 4 ∼ − 4 4 ∼ − 4 4 ∼ − 4 Lˆc (1, 2, 1/2, 0, (3 + L )), Eˆ (1, 1, 1, 0, 3+ L ), Nˆ (1, 1, 0, 0, 3+ L ) to 5 ∼ − 4 5 ∼ − 4 5 ∼ 4 cancel the B and L anomalies. With the detection of the lightest CP even Higgs h0 at LHC[4], Higgs mechanism is very convincing for particle physics, and BLMSSM is based on the Higgs mechanism. The introduced Higgs superfields Φˆ (1, 1, 0, 0, 2), ϕˆ (1, 1, 0, 0, 2) and L − L Φˆ (1, 1, 0, 1, 0), ϕˆ (1, 1, 0, 1, 0) break lepton number and baryon number spon- B B − taneously. These Higgs superfields acquire nonzero vacuum expectation values (VEVs) and provide masses to the exotic leptons and exotic quarks. To make the heavy exotic quarks unstable the superfields Xˆ(1, 1, 0, 2/3+B , 0), Xˆ ′(1, 1, 0, (2/3+B ), 0) are introduced 4 − 4

3 in the BLMSSM.

The SU(2)L doublets Hu, Hd obtain nonzero VEVs υu, υd,

H+ 1 υ + H0 + iP 0  u   √2 d d d  Hu = , Hd = . (1) 1 0 0  υu + Hu + iPu   Hd−   √2    The SU(2)L singlets ΦB, ϕB and ΦL, ϕL obtain nonzero VEVs υB, υB and υL, υL respectively,

1 0 0 1 0 0 ΦB = υB + Φ + iP , ϕB = υB + ϕ + iP . √2 B B √2 B B 1 0 0 1 0 0 ΦL = υL + Φ + iP , ϕL = υL + ϕ + iP . (2) √2 L L √2 L L Therefore, the local gauge symmetry SU(2) U(1) U(1) U(1) breaks down to the L ⊗ Y ⊗ B ⊗ L electromagnetic symmetry U(1)e. We show the superpotential of BLMSSM [12]

= + + + , WBLMSSM WMSSM WB WL WX = λ Qˆ Qˆc Φˆ + λ Uˆ cUˆ ϕˆ + λ Dˆ cDˆ ϕˆ + µ Φˆ ϕˆ WB Q 4 5 B U 4 5 B D 4 5 B B B B ˆ ˆ ˆ c ˆ ˆ ˆ c ˆc ˆ ˆ ˆc ˆ ˆ +Yu4 Q4HuU4 + Yd4 Q4HdD4 + Yu5 Q5HdU5 + Yd5 Q5HuD5 , = Y Lˆ Hˆ Eˆc + Y Lˆ Hˆ Nˆ c + Y LˆcHˆ Eˆ + Y Lˆc Hˆ Nˆ WL e4 4 d 4 ν4 4 u 4 e5 5 u 5 ν5 5 d 5 c c c +YνLˆHˆuNˆ + λN c Nˆ Nˆ ϕˆL + µLΦˆ LϕˆL ,

c c c = λ QˆQˆ Xˆ + λ Uˆ Uˆ Xˆ ′ + λ Dˆ Dˆ Xˆ ′ + µ XˆXˆ ′, (3) WX 1 5 2 5 3 5 X with representing the superpotential of the MSSM. The soft breaking terms WMSSM Lsoft of the BLMSSM are collected here[11, 12].

MSSM 2 ˜ c ˜ c 2 ˜ ˜ 2 ˜ c ˜ c 2 ˜ c ˜ c soft = (m c )IJ N ∗N m Q† Q4 m U ∗U m D ∗D L Lsoft − N˜ I J − Q˜4 4 − U˜4 4 4 − D˜ 4 4 4 2 c c 2 2 2 2 c c m Q˜ †Q˜ m U˜ ∗U˜5 m D˜ ∗D˜ 5 m L˜†L˜4 m N˜ ∗N˜ − Q˜5 5 5 − U˜5 5 − D˜5 5 − L˜4 4 − ν˜4 4 4 2 c c 2 c c 2 2 2 m E˜ ∗E˜ m L˜ †L˜ m N˜ ∗N˜5 m E˜∗E˜5 m Φ∗ ΦB − e˜4 4 4 − L˜5 5 5 − ν˜5 5 − e˜5 5 − ΦB B 2 2 2 m ϕ∗ ϕB m Φ∗ ΦL m ϕ∗ ϕL mBλBλB + mLλLλL + h.c. − ϕB B − ΦL L − ϕL L − ˜ ˜ c ˜ ˜ c ˜c ˜ ˜c ˜ + Au4 Yu4 Q4HuU + Ad4 Yd4 Q4HdD + Au5 Yu5 Q HdU5 + Ad5 Yd5 Q HuD5 n 4 4 5 5 c c c +ABQλQQ˜4Q˜ ΦB + ABU λU U˜ U˜5ϕB + ABDλDD˜ D˜ 5ϕB + BBµBΦBϕB + h.c. 5 4 4 o ˜ ˜c ˜ ˜ c ˜c ˜ ˜c ˜ + Ae4 Ye4 L4HdE + Aν4 Yν4 L4HuN + Ae5 Ye5 L HuE5 + Aν5 Yν5 L HdN5 n 4 4 5 5 c c c +AN YνLH˜ uN˜ + AN c λN c N˜ N˜ ϕL + BLµLΦLϕL + h.c. o c c c + A1λ1Q˜Q˜ X + A2λ2U˜ U˜5X′ + A3λ3D˜ D˜ 5X′ + BX µX XX′ + h.c. . (4) n 5 o

4 MSSM are the soft breaking terms of MSSM. Lsoft

A. mass matrix

From the soft breaking terms and the scalar potential, we deduce the mass squared matrix for superﬁelds X.

2 µ + S µ∗ B∗ X  X X X X    X =(X∗ X′) | | − , −L 2  µX BX µX SX   X′∗   − | | −    g2 2 S = B ( + B )(v2 v¯2 ). (5) X 2 3 4 B − B

We diagonalize the mass squared matrix for the superﬁelds X through the unitary transfor- mation,

X X µ 2 + S µ B m2 0  1     X X X∗ X∗   X1  = ZX† , ZX† | | − ZX = . (6) 2 2  X2   X′∗   µX BX µX SX   0 mX       − | | −   2 

ψX and ψX′ are the superpartners of the scalar superﬁelds X and X′. ψX and ψX′ can composite four-component Dirac spinors, whose mass term are given out[14]

ψ mass ¯˜ ˜ ˜  X  X˜ = µX XX, X = , (7) −L ¯ ′  ψX   

with µX denoting the mass of X˜.

In the BLMSSM, there are the new baryon boson, the SU(2)L singlets ΦB and ϕB. Their

superpartners are respectively λB, ψΦB and ψϕB , and they mix together producing 3 baryon neutralinos. In the base (iλB, ψΦB , ψϕB ), the mass mixing matrix MBN is obtained and diagonalized by the rotation matrix ZNB [20].

2mB vBgB v¯BgB  −  k0 M = , χ0 =  Bi  , BN  vBgB 0 µB  Bi  − −  k¯0    Bi   v¯ g µ 0     B B − B  1i 0 2i 0 3i 0 iλB = ZNB kBi , ψΦB = ZNB kBi , ψϕB = ZNB kBi . (8)

0 χBi (i =1, 2, 3) represent the mass eigenstates of baryon neutralinos.

5 The exotic quarks with charged 2/3 is in four-component Dirac spinors, whose mass matrix reads as[12]

1 λ v 1 Y v t mass  √2 Q B √2 u5 d   4′ L  t′ =(t¯′ R t¯′ R) − , (9) −L 4 5 1 1  Yu4 vu λuv¯B   t′ L   − √2 √2   5  Using the unitary transformations, the two mass eigenstates of exotic quarks with charged

2/3 are obtained by the rotation matrices Ut and Wt,

 t4L   t4′ L   t4R   t4′ R  = Ut† , = Wt† ,  t5L   t′ L   t5R   t′ R     5     5  1 λ v 1 Y v  √2 Q B √2 u5 d  Wt† − Ut = diag(mt4 , mt5 ). (10) 1 1  Yu4 vu λuv¯B   − √2 √2  2 The mass squared matrix for charged 2/3 exotic squarks ′ is obtained in our previous Mt˜ 2 work[12]. For saving space in the work, we do not show it here. ′ is diagonalized by Z˜′ Mt˜ t 2 2 2 2 2 through the formula Z†′ ′ Z˜′ = diag(m , m , m , m ). t˜ t˜ t ˜1 ˜2 ˜3 ˜4 M U U U U

B. needed couplings

To study quark EDMs, the couplings between photon (gluon) and exotic quarks(exotic squarks) are necessary. We derive the couplings between photon (gluon) and exotic quarks.

2 2 2e µ a µ a ′ ′ ¯ ¯ γ(g)q q = ti+3 γ ti+3 Fµ g3 ti+3 T γ ti+3 Gµ, (11) L − 3 Xi=1 − Xi=1 a a with Fµ and Gµ representing electromagnetic ﬁeld and gluon ﬁeld respectively. T (a = 1, , 8) are the strong SU(3) gauge group generators. Similarly, the couplings between ··· photon (gluon) and exotic squarks are also deduced

4 4 2 µ a a µ ′ ′ ˜ ˜ ˜ ˜ ˜ ˜ γ(g)˜q q˜ = e δjβFµ j∗i∂ β g3T δjβGµ j∗i∂ β. (12) L −3 j,βX=1 U U − j,βX=1 U U

From the superpotential WX , one can ﬁnd there are interactions at tree level for quark, exotic quark and X. The needed Yukawa interactions can be deduced from the superpotential WX . The couplings of quark-exotic quark-X are shown in the mass basis,

2 L I R I ′ = ( ) X t¯ P u +( ) X t¯ P u + h.c. Xt u t′ ij j i+3 L t′ ij j i+3 R L i,jX=1 N N L R ( ′ ) = λ1 (Wt†) 1 (Z )1 , ( ′ ) = λ∗(U †) 2 (Z )2 . (13) Nt ij − i X j Nt ij − 2 t i X j 6 From the superpotential WX , in the same way we can also obtain another type Yukawa couplings (quark-exotic squark-X˜)[14].

4 ˜ ˜ u¯X˜ ˜ = u¯ λ1(Zt˜′ )3∗iPL + λ2(Zt˜′ )4iPR X i. (14) L U − Xi=1 U Beyond the MSSM, there are couplings for baryon neutralino, quarks and squarks. They are deduced in our previous work[20], and can give new contributions to the quark EDMs.

3 6 0 √2 1i Ij 1i (I+3)j I (χ qq˜)= gBχ¯ 0 (Z Z ∗PL Z ∗ Z ∗PR)u U˜ ∗ + H.c. (15) B Bi NB U˜ NB U˜ j L I,iX=1 jX=1 3 −

III. FORMULATION

Using the eﬀective Lagrangian[21] method, one obtains the fermion EDM df from i = d fσµν γ fF , (16) LEDM −2 f 5 µν

with Fµν representing the electromagnetic field strength, f denoting a fermion field. It is obviously that this effective Lagrangian is CP-violating. In the fundamental interactions, this CP-violating Lagrangian can not be obtained at tree level. Considering the CP-violating electroweak theory, one can get this effective Lagrangian from the loop diagrams. The a µν a chromoelectric dipole moment (CEDM) fT σ γ5fGµν of quark can also give contribution a to the quark EDM. Gµν denotes the gluon field strength. To describe the CP-violating operators obtained from loop diagrams, the effective method is convenient. The coefficients of the quark EDM and CEDM at the matching scale µ should be evolved down to the quark mass scale with the renormalization group equations. At matching scale, we can obtain the effective Lagrangian with the CP-violating operators. The effective Lagrangian containing operators relating with the quark EDM and CEDM are

eff = Ci (Λ) i (Λ) , L Xi O = qσµνP qF , = qσµν P qF , O1 L µν O2 R µν = qT aσµν P qGa , = qT aσµν P qGa . (17) O3 L µν O4 R µν

with Λ representing the energy scale, where the Wilson coeﬃcients Ci(Λ) are evaluated. In our previous work[14], we have studied the neutron EDM in the CP-violating BLMSSM, where the contributions from baryon neutralino-squark and X˜-exotic squark are

7 neglected, because they are all small in the used parameter space. Here we take into account all the contributions at one loop level to study the c and t EDMs. In the CP-violating BLMSSM, the one-loop corrections to the quark EDMs and CEDMs can be divided into six types according to the quark self-energy diagrams. We divide the quark self-energy diagrams according to the virtual particles as: 1 gluino-squark, 2 neutralino-squark, 3 chargino-squark, 4 X-exotic quark, 5 baryon neutralino-squark, 6 X˜-exotic squark. From the quark selfenergy diagrams, one obtains the needed triangle diagrams by attaching a photon or gluon on the internal lines in all possible ways. After the calculation, we obtain the eﬀective Lagrangian contributing to the quark EDMs and CEDMs. The BLMSSM is larger than MSSM and includes the MSSM contributions. In Fig.(6), we plot all the one loop self energy diagrams of the up-type quark.

FIG. 1: In the BLMSSM, one loop self energy diagrams are collected here, and the corresponding triangle diagrams are obtained from them by attaching a photon or a gluon in all possible ways.

In this section, we show the one loop corrections to the quark EDMs (CEDMs). The one loop chargino-squark contributions are

6 2 m ± γ I eα D D χk ± † † dχ (u )= 2 VUDVDU Im (AC )k,i(BC )i,k 2 k 4πs X X m W i k D˜i 2 2 1 mχ± mχ± k + k , ×h − 3B m2 A m2 i D˜i D˜ i 2 6 2 m ± m ± g I g3α D D χk χk ± † † dχ (u )= 2 VUDVDU Im (AC )k,i(BC )i,k 2 2 , k 4πs X X m B m W i k D˜i D˜ i D muI Ji 2k (AC )k,i = (ZD˜ ) (Z+) , √2mW sβ D mdI (J+3)i 2k Ji 1k (BC )k,i = (ZD˜ ) (Z ) (ZD˜ ) (Z ) . (18) √2mW cβ − − −

2 Here α = e /(4π),sW = sin θW , cW = cos θW , θW is the Weinberg angle, V is the CKM matrix. We deﬁne the one loop functions (r) and (r) as[9] A B

2 1 (r) = [2(1 r) ]− [3 r + 2 ln r/(1 r)], A − − − 2 1 (r) = [2(r 1) ]− [1 + r +2r ln r/(1 r)]. (19) B − −

mD˜ i (i = 1 ... 6) are the squarks masses and mχ0 (k = 1, 2, 3, 4) denote the eigenvalues k of neutralino mass matrix. ZD˜ i is the rotation matrix to diagonalize the mass squared matrix for the down type squark. Z and Z+ are the rotation matrices to obtain the mass − eigenstates of charginos.

9 We show the gluino-squark corrections to the quark EDMs and CEDMs 6 2 γ I 4 (I+3)i Ii iθ3 mg˜ mg˜ dg˜(u )= eαs Im ( U˜ ) ( U˜ ) ∗e− | 2 | | 2| , −9π X Z Z m B m i=1 U˜i U˜i 6 2 g I g3αs (I+3)i Ii iθ3 mg˜ mg˜ dg˜(u )= Im ( U˜ ) ( U˜ ) ∗e− | 2 | | 2| , (20) 4π X Z Z m C m i=1 U˜i U˜i 2 2 with α = g /(4π). Z ˜ is the matrix for the up type squarks, with the deﬁnition † m ˜ = s 3 U ZU˜ U˜ ZU diag(m2 ,...,m2 ). The concrete form of the loop function (r) is[9] q˜1 q˜6 C 2 1 (r) = [6(r 1) ]− [10r 26 (2r 18) ln r/(r 1)]. (21) C − − − − − To check the functions (r), (r) and (r) in the Ref.[9], we calculate the one loop A B C triangle diagrams using the eﬀective Lagrangian method. In the calculation, we use the approximation 1 2k p + p2 4(k p)2 =1 · + · , (22) (k + p)2 m2 − k2 m2 (k2 m2)2 − − − with k representing the loop integral momentum and p representing the external momentum. It is reasonable because the internal particles are at the order of TeV, and the external quark 2 mt is lighter than TeV, even for t quark. The ratio 2 0.03 is small enough to use the 1T eV ∼ approximation formula. For the diagram that the photon is just attached on the internal charged Fermions, our result corresponding to (r) is the function a[F,S] A Λ2 k2 a[F,S]= NP dk4 iπ2 Z (k2 m2 )3(k2 m2 ) − F − S F 2 +2S2 log(F ) 4FS +3S2 2S2 log(S) = − − , − 2(F S)3 − (23)

2 2 mF mS with the deﬁnition F = 2 and S = 2 . λNP represents the the energy scale of the new ΛNP ΛNP 2 2 2 mF physics. In order to compare with the function (r), we use λNP = mS,S 1, F m2 = r A → → S and obtain r2 4r + 2 log(r)+3 a[r, 1] = − = (r). (24) − 2(r 1)3 A − When the photon is just emitted from the internal charged scalars, our result for (r) is B Λ2 m2 b[F,S]= NP dk4 S iπ2 Z (k2 m2 )(k2 m2 )3 − F − S F 2 2FS log(F )+2FS log(S) S2 = − − , (25) 2(F S)3 − 10 With the same approach as that of a[F,S], b[F,S] turns into the form r2 2r log(r) 1 b[r, 1] = − − = (r). (26) 2(r 1)3 B − (r) is obtained from the diagrams that the photon is attached on both the internal C charged Fermions and charged scalars. Therefore, (r) is the linear combination of (r) C A and (r), B 1 (r)= (r) 3 (r). (27) C 3B − A From the above discussion, the results in Ref.[9] are the same with our results. In our calculation, we do not ignore the mass of the external fermion. So, it is clear that the analytical expressions for the quark EDM in this work are practicable for both c and t. Similarly, the contributions from the one loop neutralino-squark diagrams are also obtained 2 6 4 m 0 m 0 γ I eα χk χk † dχ0 (u )= 2 2 Im (AN )k,i(BN )i,k 2 2 , k 12πs c X X m B m W W i=1 k=1 U˜i U˜i 2 6 4 m 0 m 0 g I g3α χk χk † dχ0 (u )= 2 2 Im (AN )k,i(BN )i,k 2 2 , k 8πs c X X m B m W W i=1 k=1 U˜i U˜i 4 (I+3)i 1k muI cW Ii 4k (AN )k,i = sW (ZU˜ ) (ZN ) + (ZU˜ ) (ZN ) , −3 mW sβ Ii sW 1k 2k muI cW (I+3)i 4k (BN )k,i =(ZU˜ ) ( (ZN ) ∗ + cW (ZN ) ∗)+ (ZU˜ ) (ZN ) ∗. (28) 3 mW sβ

ZN is the mixing matrix to get the eigenvalues m 0 (k = 1, 2, 3, 4) of neutralino mass χk matrix. In the MSSM, there are also the front three type contributions Eqs.(18)(20)(28). At one loop level, there are three new type corrections beyond MSSM. The corrections from the virtual X and exotic up-type quark has been deduced in the work[14]

2 2 γ I eλ1λ2 mti+3 mti+3 d (u )= Im (Wt)1i(ZX )∗ (Ut)∗ (ZX )2j , Xj 2 2 1j 2i 2 24π i,jX=1 mXj A mXj 2 2 g I g3λ1λ2 mti+3 mti+3 dX (u )= Im (Wt)1i(ZX )∗ (Ut)∗ (ZX )2j , (29) j 2 2 1j 2i 2 16π i,jX=1 mXj A mXj mti+3 and mXi (i=1,2) are mass eigenvalues of the exotic up type quarks and X superﬁelds.

Wt, Ut and ZX are the mixing matrices deﬁned in the Eqs.(6)(10). The one loop baryon neutralino and up-type squark contributions read as

2 6 2 2 m i m i γ I egB 1i 2 (I+3)j Ij χB χB ∗ dχB (u )= Im (ZNB ) Z ˜ Z ˜ 2 2 , −108π2 X X U U m B m i=1 j=1 U˜j U˜j

11 2 6 2 2 m i m i g I g3gB 1i 2 (I+3)j Ij χB χB ∗ dχB (u )= Im (ZNB ) Z ˜ Z ˜ 2 2 . (30) −72π2 X X U U m B m i=1 j=1 U˜j U˜j m i (i=1,2,3) are the eigenvalues of baryon neutralino masses. χB The exotic up-type squark and X˜ can also contribute to the c(t) EDM and CEDM 2 2 γ eλ1λ2 3i 4i mX˜ mX˜ d = Im (Zt˜′ ) ∗(Zt˜′ ) ∗ , X˜ 24π2 m2 Bm2 Xi=1 t˜i+3 t˜i+3 2 2 g g3λ1λ2 3i 4i mX˜ mX˜ d = Im (Zt˜′ ) ∗(Zt˜′ ) ∗ , (31) X˜ 16π2 m2 Bm2 Xi=1 t˜i+3 t˜i+3 ˜ with mX˜ and mt˜i+3 denoting the masses of X and exotic up-type squark respectively. Using the renormalization group equations[22], we evolve the coeﬃcients of the quark EDM and CEDM at matching scale µ down to the quark(c, t) mass scale[23]

γ γ g g dq (Λχ )=1.53dq (Λ), dq (Λχ )=3.4dq (Λ), (32)

The quark CEDMs can contribute to the quark EDMs at low energy scale. Therefore, they must be taken into account in the numerical calculation and the formula is[24] e d = dγ + dg. (33) c c 4π c

IV. THE NUMERICAL RESULTS

Here, the results are studied numerically. We take into account not only the experiment constraints from Higgs and neutrino, but also our previous works in this model. From ATLAS collaboration, m 1460 GeV is the updated bound on the gluino mass[25]. The g˜ ≥ parameters are supposed as

2 2 m1 = m2 = ABQ = ABU = 1 TeV, BX = 500 GeV, mD = δij TeV , (i, j =1, 2, 3),

Au = Ad = Au′ = Ad′ = 500 GeV,Yd4 = Yd5 =0.7Yb, λQ = λu =0.5,

2 2 2 2 2 3 m = m = m = m = 1 TeV , B4 = , Au4 = Au5 = 500 GeV. (34) Q˜4 Q˜5 U˜4 U˜5 2

A. c quark EDM

For the c quark EDM, we use the following parameters as

tan β = 10, µ = 800 GeV, mg˜ = 1600 GeV,

tan βB =2, vBt = 3 TeV,Yu4 = Yu5 =0.7Yt. (35)

12 The baryon neutralino and squarks can give contributions to c quark EDM, which is

relevant to the parameters gB and mB. mB representing baryon gaugino masses, can have

nonzero CP-violating phase θmB . Both mB and gB inﬂuence the baryon neutralino masses.

Furthermore, gB is the coupling constant for the quark-squark-baryon neutralino. So, with θ = 0.5π, λ = λ = 0.1, µ = µ = 3 TeV, m2 = m2 = δ TeV2 for (i, j = mB − 1 2 X B Q U ij 1, 2, 3), we study c quark EDM versus gB. If we do not mention the other CP-violating phases, it indicates the other CP-violating phases are zero. In Fig.2, the numerical results iθ corresponding to m = (1, 2, 3) e mB TeV are plotted by the solid line, dotted line and B × dashed line respectively. In the whole, the three lines are all increasing functions of gB. The dotted line and the dashed line vary slightly. The solid line largen quickly with the increasing g . These three lines also imply the results are suppressed by large m . The B | B|

obtained numerical results from the nonzero CP-violating phase θmB are at the order of 21 10− e.cm, which are four orders smaller than the upper bound of c quark EDM.

L 4 e.cm H

21 3 10 ´ c d 2

0 0.0 0.2 0.4 0.6 0.8 gB

FIG. 2: The one loop corrections to c EDM versus gB with θm = 0.5π, the solid line, dotted B − iθm line and dashed line corresponding to mB = (1, 2, 3) e B TeV respectively. ×

Here, we discuss the eﬀects from the new phase θµB of µB, it can inﬂuence the baryon iθ neutralino masses. Based on the supposition µ = 2 e µB TeV (θ = 0.5π), λ = λ = B × µB 1 2 2 2 2 2 0.5, µX = 3 TeV, mQ = 4δij TeV , mU = 2δij TeV (i, j = 1, 2, 3), the results versus 1 mB are shown as the solid line for gB = 3 . The solid line decreases quickly in the region

1000GeV < mB < 1300GeV. When mB > 1300GeV, the change extent of the solid line is 1 small. The dotted line and the dashed line respectively represent the results for gB = 5 and

13 1 gB = 10 , and they are both the slowly decreasing functions of mB. Generally speaking, the 21 results are around 10− e.cm that are at the same order of Fig.2.

L 5 e.cm

H 4 21 10

´ 3 c d 2

0 1500 2000 2500 3000 3500 mBGeV

FIG. 3: The one loop corrections to c EDM versus mB with θµB = 0.5π, the solid line, dotted line

and dashed line corresponding to gB = (1/3, 1/5, 1/10) respectively.

λ1 and λ2 are important parameters for the couplings: quark-exotic quark-X and quark- exotic squark-X˜. Therefore, the numerical results maybe be inﬂuenced obviously by the 1 varying λ1 and λ2. For simplicity, we suppose λ1 = λ2 = Lam, gB = 3 , mB = 1TeV, µB = 2 2 2 3TeV, mQ = mU = δijTeV for (i, j = 1, 2, 3). With the nonzero CP-violating phase

iθX θX = (0.5π, 0.3π, 0.1π) and µX = e TeV, the results versus Lam are denoted respectively by the solid line, dotted line and dashed line. The three lines are the increasing functions 17 of Lam and the results are around 10− e.cm. As Lam > 0.8, the solid line even exceeds its 17 the upper bound 5.0 10− e.cm. From the Figs.(2, 3, 4), one can ﬁnd the eﬀects from θ × X

are much larger than the eﬀects from θµB and θmB .

B. t quark EDM

To calculate the t quark EDM numerically, we use the parameters as

1 g = , m = 1 TeV, µ = 3 TeV. (36) B 3 B B

The quark-gluino-squark coupling corrections to t quark EDM are shown in Eq.(20). tan β is important, because it inﬂuences the masses of chargino, neutralino, squark and so on.

14 8

6 L e.cm H 4 17 10 ´ c d 2

0 0.0 0.2 0.4 0.6 0.8 1.0 Lam

FIG. 4: The one loop corrections to c EDM versus Lam, while the solid line, dotted line and dashed

i(0.5π,0.3π,0.1π) line denoting the results with µX = e TeV respectively.

The mixing matrices of squarks and exotic squarks have relation with tan β. The absolute value of gluino mass also inﬂuence the results obviously from Eq.(20). Using the parameters θ = 0.5π, µ = 800 GeV, tan β = 2,Y = Y = 0.7Y , λ = λ = 0.5, µ = 3 − B u4 u5 t 1 2 X 2 2 2 2 2 1 TeV, VBt = 3300 GeV, mQ = δij1500 GeV , mU = δijTeV with (i, j = 1, 2, 3), for t quark EDM we plot the results versus mg˜. In Fig.(5), the solid line (tan β = 5), dotted line

(tan β = 10) and dashed line (tan β = 15) are all decreasing functions of mg˜. During the m region (1500 2000) GeV, the results of the three lines shrink quickly. Near the point g˜ ∼ 16 mg˜ = 1480 GeV, the theoretical predictions are at the order of 10− e.cm and even reach 15 10− e.cm. On the other hand, the extent of the influence from tan β is not large. The effects to the t quark EDM from the µ parameter are also of interest. µ is included in 2 2 the mass matrices of chargino and neutralino. On the other hand mQ and mu can affect the masses and mixings of the up-type squark. For simplification of the numerical discussion, 2 2 2 we adopt the relation mQ = mU = Mus and use the parameters tan β = 10, mg˜ =

1600 GeV, tan βB =1.5, VBt = 3600 GeV,Yu4 = Yu5 =0.7Yt, λ1 = λ2 =0.1, µX = 1 TeV. As θ = 0.5π and µ = (1500, 2500, 3500)eiθµ GeV, the numerical results for t quark µ − EDM are plotted respectively by the dotted line, solid line and dashed line. The results are 18 all decreasing functions of Mus and at the order of 10− e.cm as Mus < 3500 GeV. The absolute value of µ also inﬂuences the results, and the extent is small.

The exotic squarks and exotic quarks are in connection with Yu4 and Yu5 , and the related

15 10

1 L e.cm H 17

10 0.1 ´ t d

0.01

1600 1800 2000 2200 2400 mgGeV

FIG. 5: The one loop corrections to t EDM versus mg with θ = 0.5π, the solid line, dotted line ˜ 3 − and dashed line corresponding to tan β = (5, 10, 15).

L 5 e.cm

H 4 18 10 ´

t 3 d 2

0 2000 3000 4000 5000 MusGeV

FIG. 6: The one loop corrections to t EDM versus Mus with θµ = 0.5π, the dotted line, solid − line and dashed line corresponding to the results with µ = (1500, 2500, 3500)eiθµ GeV.

contributions are shown in Eqs.(29,31). As discussed in the front subsection, nonzero θX can give large contributions. In Fig.7 we plot the solid line, dotted line and dashed line versus Yu45 with Y = Y = Yu45 Y t, θ = (0.5π, 0.2π, 0.05π) and µ = 1 eiθX TeV. The u4 u5 ∗ X X other used parameters are tan β = 10, µ = 800GeV, mg˜ = 1600 GeV, tan βB = 2, VBt = 2 2 2 3 TeV, λ1 = λ2 = 0.5, mQ = mU = δijTeV with (i, j = 1, 2, 3). They are all increasing functions of Yu45 and as Yu45 > 0.8 the numerical results largen quickly. At the point 17 Yu45 = 0.7, the solid line and dotted line are larger than 1 10− e.cm. The three lines ×

16 18 are at the order of 10− e.cm, when Yu45 is small. For c and t quark EDMs, in our used parameter space the CP-violating phases θ3 and θX are important and can provide large contributions.

8 L 6 e.cm H 17 10

´ 4 t d

0 0.0 0.2 0.4 0.6 0.8 1.0 Yu45

iθ FIG. 7: The one loop corrections to t EDM versus Yu45 with µX = e X TeV, the solid line, dotted line and dashed line corresponding to the results with θX = (0.5π, 0.2π, 0.05π).

V. DISCUSSION

In the CP-violating BLMSSM, there are new CP-violating phases θX , θµB , θmB beyond MSSM. For the c quark EDM, we consider the conditions θ = 0, θ = 0 and θ = 0 X 6 µB 6 mB 6 respectively. The results show θX can give large contributions, even reach the experiment 17 upper bound(5 10− e.cm) for c EDM. The eﬀects produced from θ and θ are at the × mB µB 21 order of 10− e.cm, which are much smaller than those from θX . For the t quark EDM the CP-violating phases θ , θ and θ are studied. In both θ = 0 and θ = 0 conditions, we 3 µ X X 6 3 6 17 ﬁnd dt at the order of 10− e.cm. Especially for nonzero θ3 with mg˜ near its lower bound, t 16 quark EDM can reach 10− e.cm and even larger. They are both larger than the results for θ = 0. µ 6 In BLMSSM, at one loop level there are three type contributions (1. the virtual X and exotic up-type quark, 2. baryon neutralino and up-type squark, 3. exotic up-type squark and X˜) to quark EDM beyond MSSM. For the contributions beyond MSSM, to obtain large dc and dt the CP-violating phase θX should be nonzero in our used parameter space. With

17 17 only θ = 0, the nonzero contributions come from Eqs.(23)(25). In the Fig.4 d 10− e.cm X 6 c ∼ 17 and in the Fig.7 d 10− e.cm, the EDMs d and d are of the same order of magnitude. t ∼ c t In the other works, the EDMs dc and dt should be of diﬀerent order. What is the reason in this work? The reason can be found from the couplings in Eqs.(13,14). The coupling constants λ1 and λ2 are important parameters. In our numerical calculation, we adopt that

the values of λ1(λ2) for c quark are same with the values of λ1(λ2) for t quark. That is to

say, for the up type quark generation 2 and generation 3, the adopted values of λ1(λ2) are the same.

From these numerical results and our previous work on neutron EDM, we ﬁnd θ3 and 2 2 θX are important CP-violating phases. tan β, Yu4 ,Yu5 , mg˜, mQ, mU , λ1, λ2, VBt are also important. In the whole, our numerical results are large to be detected in the future. The work can conﬁne the parameter space in this model and possess meaning to the relevant experiments for c and t quark EDMs.

Acknowledgments

The work has been supported by the National Natural Science Foundation of China (NNSFC) with Grant No. 11275036, No. 11535002, No. 11605037, No. 11647120, the Fund of Natural Science Foundation of Hebei Province(A2011201118), Natural Science Fund of Hebei University with Grant No. 2011JQ05, No. 2012-242. Hebei Key Lab of Optic- Electronic Information and Materials, the midwest universities comprehensive strength pro- motion project.

[1] J.M. Christensen, J.W. Cronin, V.L. Fitch, R. Turlay, Phys. Rev. Lett. 13 (1964) 138; K. Abe, et al., [Belle Collaboration], Phys. Rev. Lett. 87 (2001) 091802; B. Aubert, et al., [BaBar Collaboration], Phys. Rev. Lett. 89 (2002) 201802. [2] J. Ellis, S. Ferrara, and D.V. Nanopoulos, Phys. Lett. B. 114 (1982) 231; J. Polchinski and M. B. Wise, ibid. 125 (1983) 393; P. Nath, Phys. Rev. Lett. 66 (1991) 2565; Y. Kizukuri and

18 N. Oshimo, Phys. Rev. D. 46 (1992) 3025; 45 (1992) 1806. [3] E. Shabalin, Sov. J. Nucl. Phys. 28 (1978) 75; F. Sala, JHEP 1403 (2014) 061; I. Khriplovich, Phys. Lett. B 173 (1986) 193-196; A. Czarnecki and B. Krause, Phys. Rev. Lett. 78 (1997) 4339-4342. [4] CMS Collaboration, Phys. Lett. B 716 (2012) 30; ATLAS Collaboration, Phys. Lett. B 716 (2012) 1. [5] J. Rosiek, Phys. Rev. D 41 (1990) 3464; arXiv:hep-ph/9511250; H.P. Nilles, Phys. Rept. 110 (1984) 1; H.E. Haber and G.L. Kane, Phys. Rept. 117 (1985) 75-263. [6] N. Escudero, D.E. Lopez-Fogliani, C. Munoz and R.R. de Austri, JHEP 12 (2008) 099; Hai- Bin Zhang, Tai-Fu Feng, Fei Sun, et al., Phys. Rev. D 89 (2014) 115007. [7] A. Pilaftsis, Phys. Rev. D 58 (1998) 096010; M. Carena, J. Ellis, A. Pilaftsis and C.E.M. Wagner, Nucl. Phys. B 586 (2000) 92; Tai-Fu Feng, Lin Sun, Xiu-Yi Yang, Nucl. Phys. B 800 (2008) 221; Tai-Fu Feng, Lin Sun, Xiu-Yi Yang, Phys. Rev. D 77 (2008) 116008. [8] I. Altarev et al., Phys. Lett. B 276 (1992) 242; K. Smith et al., Phys. Lett. B 234 (1990) 191; C. Patrignani et al. (Particle Data Group), Chin. Phys. C 40 (2016) 100001. [9] T. Ibrahim and P. Nath, Phys. Rev. D 57 (1998) 478; Erratum Phys. Rev. D 58 (1998) 019901; Phys.Lett. B 418 (1998) 98; Erratum: Phys.Lett. B 460 (1999) 498. [10] P.F. Perez, Phys. Lett. B 711 (2012) 353; J.M. Arnold, P.F. Perez, B. Fornal, and S. Spinner, Phys. Rev. D 85 (2012) 115024; P.F. Perez, M.B. Wise. Phys. Rev. D 84 (2011) 055015. [11] P.F. Perez and M.B. Wise, JHEP 1108 (2011) 068; Phys. Rev. D 82 (2010) 011901. [12] Tai-Fu Feng, Shu-Min Zhao, Hai-Bin Zhang, Yin-Jie Zhang, Yu-Li Yan, Nucl.Phys. B 871 (2013) 223. [13] S.M. Zhao, T.F. Feng, H.B. Zhang, et al., JHEP 11 (2014) 119. S.M. Zhao, T.F. Feng, H.B. Zhang, et al., Phys. Rev. D 92 (2015) 115016. Shu-Min Zhao, Tai-Fu Feng, Xing-Xing Dong, et al., Nucl. Phys. B 910 (2016) 225C239. [14] S.M. Zhao, T.F. Feng, B. Yan, et al., JHEP 10 (2013) 020; F. Sun, T.F. Feng, S.M. Zhao, et al., Nucl. Phys. B 888 (2014) 30; S.M. Zhao, T.F. Feng, X.J. Zhan, et al., JHEP 07 (2015) 124. [15] A. Cordero-Cid, J.M. Hernandez, G. Tavares-Velasco, J.J. Toscano, J.Phys.G 35 (2008) 025004. [16] A.G. Grozin, I.B. Khriplovich, A.S. Rudenko, Nucl. Phys. B 821 (2009) 285-290.

19 [17] A.E. Blinov, A.S. Rudenko, Nucl. Phys. Proc. Suppl. 189 (2009) 257-259. [18] Xiao-Jun Bi, Tai-Fu Feng, Xue-Qian Li, et al., hep-ph/0412360. [19] Filippo Sala, JHEP 1403 (2014) 061; Z.Z. Aydin, U. Erkarslan, Phys. Rev. D 67 (2003) 036006. [20] Xing-Xing Dong, Shu-Min Zhao, Hai-Bin Zhang, et al., Chin. Phys. C 40 (2016) 093103. [21] Tai-Fu Feng, Xue-Qian Li, Jukka Maalampi, Xin-min Zhang, Phys. Rev. D 71 (2005) 056005. [22] R. Arnowitt, J. Lopez, and D. V. Nanopoulos, Phys. Rev. D 42 (1990) 2423; R. Arnowitt, M. Duﬀ and K. Stelle, ibid. 43 (1991) 3085. [23] R. Arnowitt, J. Lopez, and D. V. Nanopoulos, Phys. Rev. D. 42 (1990) 2423; R. Arnowitt, M. Duﬀ and K. Stelle, ibid. 43, 3085(1991). [24] A. Manohar, and H. Georgi, Nucl. Phys. B. 234 (1984) 189. [25] M. Aaboud et.al,ATLAS Collaboration, Phys. Rev. D 94 (2016) 052009.