Bled Workshops in Physics Vol. 19, No. 2 A Proceedings to the 21 st Workshop What Comes Beyond ... (p. 247) Bled, Slovenia, June 23-July 1, 2018 11 Electric Dipole Moment and Dark Matter in a CP Violating Minimal Supersymmetric SM T. Shindou * Division of Liberal-Arts, Kogakuin University Nakanomachi 2665-1, Hachioji, 192-0015, Tokyo, Japan Abstract. We consider a dark matter scenario in the minimal supersymmetric standard model with CP violation where the Bino-like neutralino is a dark matter and its annihilation cross section is enhanced enough to reproduce the observed relic abundance of the dark matter through heavy Higgs bosons exchange. In this benchmark scenario, we examine the electric dipole moments of the electron, the mercury, and the neutron. We also consider the spin-independent cross section for the dark matter scattering with nuclei. We show that the electric dipole moments will be very powerful tool to explore the parameter space in this model, even when most of the new particles are very heavy. Povzetek. Avtor obravnava model za temno snov v okviru minimalnega supersimetricnega standardnega modela s krsitvijo CP, v katerem temno snov tvori vrsta nevtralina z dovolj velikim sipalnim presekom za anihilacijo z izmenjavo teZkih Higsovih bozonov, da da njegova gostota ustreze izmerjeni pogostosti temne snovi. V tem modelu oceni električne dipolne momente elektrona, jedra zivega srebra in nevtrona. Obravnava od spina neodvisne sipalne preseke za sipanje te temne snovi na jedrih. Ugotovi, da je elektricni dipolni moment elektrona koristno orodje za raziskavo prostora parametrov tega modela tudi v primeru, ce je vecina delcev v tem supersimetricnem modelu zelo masivnih. Keywords: dark matter, neutralino, EDM, MSSSM 11.1 Introduction Though there is no evidence of supersymmetry (SUSY) at the LHC experiments, SUSY is still an attractive candidate of physics beyond the Standard Model (SM). There are several motivations to consider the minimal SUSY Standard Model (MSSM) than it in the SM. For example, (i) the gauge coupling unification is improved in the MSSM, (ii) quadratic divergence in the scalar sector is cancelled, (iii) spin-0 scalar fields are naturally introduced, (iv) MSSM provides a well-defined ultraviolet picture of type-II two Higgs doublet model, (v) If R-parity is unbroken, the lightest SUSY particle (LSP) can be a dark matter (DM) candidate, and so on. Among such attractive motivations, we focus on the point (v). In the SM, there are several unsolved problems and one of the most serious problems is absence of * E-mail: shindou@cc.kogakuin.ac.jp 248 T. Shindou the DM candidate. In the MSSM, all the SM particles are R-parity even and all the SUSY partner particles are R-parity odd, so that the lightest R-parity odd particle cannot decay. Therefore, unbroken R-parity guarantees the stability of the LSP which can be a DM. Several different candidates can be considered in the MSSM such as the neutralino, the gravitino, the axino, the saxion and the sneutrino. In this talk, we briefly review the analysis studied in Ref. [1] where a neutralino DM scenario is considered. In the neutralino DM scenario, the relic abundance of the LSP tends to be much more than the observed value. In order to realise the observed relic abundance of the DM, a mechanism to enhance the annihilation cross section of LSP is necessary. For example, following scenarios are sometimes considered: (i) neutralinos annihilate significantly through SU(2) gauge interaction, or (ii) annihilation cross section of Bino-like neutralino is enhanced with a particular mass spectrum of other associated particles. In the former class, one possible case is the Higgsino-like neutralino DM scenario with the mass of about 1 TeV. In this scenario, phenomenology such as the direct detection of DM, contribution to the EDMs, and collider signals have been studied in Ref. [2]. There is another possibility that a neutralino DM whose main component is Bino annihilates through heavy Higgs boson resonance [3-7]. We, here, focus on the second case. In this scenario, masses of the heavy Higgs boson are about twice of the mass of the neutralino DM. This Bino-like neutralino also contains small Higgsino component so that the neutralino can directly be searched through Higgs bosons exchange by the spin-independent scattering off nucleus [8]. We consider the MSSM with CP violating phases. In this case, the CP violating phases can significantly affect the electric dipole moments (EDM). Therefore the EDMs are powerful tools to explore the CP violating phases in the model. In this talk, we examine the electron EDM, the nucleon EDM, and the mercury EDM. CP phases can also contribute to the DM-nucleon spin-independent scattering cross section. Since the pseudo scalar exchange process is strongly suppressed in the non-relativistic limit, the spin-independent cross section is suppressed with a significant size of CP phase. 11.2 The benchmark of our analysis The superpotential and the soft SUSY breaking terms in the MSSM are given by[9] W = eab[(ye)ijHaLtbE j + (yd)ijHaQtbD j + (yu)ijHa QbUj (11.1) - ^HaHb], 11 Electric Dipole Moment and Dark Matter... 249 and Lsoft =- M B B - W aw a - M23 g Ag A - m2H, H|aHa + m2H2 H?aH2a - q?La(M2 )tJq?L - ^a(M?)ij??L -iJ-iR(MUu)ijHiR - diR(M?)ijd?R - eiR(Mj2)ijgiR - £ab [(Te)ijHa£tbLejR + (TdlijHfqtbLdjR +(Tu)ijHaqbL-UjR + m2HaHb + h.c.] , (11.2) respectively. In the following, we ignore the Yukawa couplings except for the third generation quarks and leptons. Then yt, yb, and yT denote the Yukawa couplings of top, bottom, and tau, respectively. We also neglecting the flavor mixing in the soft SUSY breaking terms, we take flavor diagonal soft scalar masses as M?. = (M? )ii, M?. = (M?)ii, M2. = (M2)ii, M2. = (M2 )ii, and M^. = (M?)ii. For the trilinear couplings, A parameters defined by (Tu)33 = ATyt, (Td)33 = ATyb, and (Te )33 = ATyT are used. Since we consider the CP violating case, the each parameter in the above superpotential and the soft SUSY breaking Lagrangian can be a complex number. The mass of the SM-like Higgs boson in the MSSM is calculated by the input parameters in the superpotential and the SUSY breaking Lagrangian. In order to reproduce the observed mass value mh = 125 GeV, we take tan |3 := (H2)/(H } = 30 and we fix the stop mass parameters as M?3 = 7 TeV, Mt := Mu3 = 7 TeV and At = 10 TeV. The other SUSY particles are irrelevant to the mass of the SM-like Higgs boson as well as the DM relic density. Therefore we can take their masses much heavier than stop. In such a case, they are decoupled from low energy observables. Here we take masses of the other sfermions as 100 TeV and M2 = M3 = 10 TeV. In our analysis, we focus on the Bino-like DM with the Higgs funnel scenario so that the heavy Higgs boson mass is close to twice the mass of the DM. In the scenario, the Bino-like neutralino rapidly annihilate through the heavy Higgs bosons resonance and the appropriate cosmic abundance for DM is reproduced. In addition, the masses of heavier neutral Higgs bosons, mH and mA, are close to the charged Higgs boson mass mH± in the MSSM. Thus we fix mH± to be twice of Bino mass parameter Mi. Note that the X-X-Higgs boson coupling depends on non-vanishing Higgsino component in the neutralino. We choose to reproduce the correct amount of DM relic density as nDMh2 = 0.1198 ±0.0015 [10]. As a consequence of these fact, both the Bino mass |M11 and the Higgsino mass should be of the order of TeV. We consider M1 as a free parameter and solve from the measured dark matter energy density. In the following, we summarise our benchmark parameter set: |M21 = |M3| = 10 TeV, M?,,2 = Mai,2 = Maii2,3 = M^,2,3 = Meli2,3 = 100 TeV, Mq3 = Mt = 7 TeV, At = 10 TeV, mH± = 2M1, tan | = 30. (11 3) (11 4) (11 5) (11 6) (11 7) (11 8) 250 T. Shindou The other A-terms are zero. With this parameter set, CP phases in the five parameters, M1, M2, M3,At), may be relevant to our analysis of EDMs and the spin-independent cross section. The CP phases of these parameters are described as $Ml, $M2 , $M3 , $At), respectively, where each phases of a quantity X are defined by X = |X|eT*X. Note that some of those CP phases are unphysical. It is known that there is a rephasing degree of freedom in the MSSM. Actually, all the physical quantities are described by the following combinations of the parameters, By using the rephasing degree of freedom, without loss of generality, we can take the basis of CP phases as $M3 = 0. We also take $At = 0 for simplicity. In general, the CP phase $At also significantly contributes to the predictions of the EDMs. However, in our benchmark parameter set given in Eqs. (11.3) - (11.8), the contribution from $At is strongly suppressed because the mass splitting between two stops is small. Therefore we scan the following four parameters, 11.3 Numerical analysis In calculations of dark matter thermal relic density and the Higgs mass, we use micrOMEGAs 4.3.5 [11] with CPsuperH2.3 [12]. The Higgs mass is almost fixed to be 125 GeV in our benchmark point. When we scattered the parameters, we pick up the parameter sets which reproduce the correct DM relic abundance and the correct Higgs mass. Then we calculate the electron EDM, the neutron EDM, and the mercury EDM. We also discuss the scattering cross section for the direct detection experiments. Since the sfermions are too heavy to contribute to the EDMs via one-loop diagrams, the two-loop Barr-Zee diagrams provide dominant contributions unless Wino, stop, and sbottom masses are heavy enough to be decoupled. In Fig. 11.1, we show our numerical results. We can see the M1 dependence by comparing the left panels and the right panels where M1 = 1 TeV and 2 TeV, respectively. It is easily seen that larger Mi weaken the constraint from EDM experiments. The Bino mass M1 is approximately identified to be the mass of the dark matter neutralino. Then for larger M1, heavy Higgs bosons and Higgsinos become heavier, and the contributions to the EDMs become smaller. We also discuss the ^ and $M2 dependence of the EDMs. The left panels in Fig. 11.1 shows the electron EDM, the mercury EDM, and the neutron EDM with $Ml = 0. The shaded regions are already excluded by the current upper bound on the EDMs. We find the combination of the electron EDM and the mercury EDM exclude the arg(MiMi) , arg(MtAi) , arg(^Mt) , arg(^At) , (i,j = 1,2,3) . (11.9) (|Mi|, $ (11.10) 11 Electric Dipole Moment and Dark Matter... 251 large region of the parameter space. Both ^ and §m.2 cannot be large. We also find that the electron EDM strongly depends on . On the other hand, dependence of the mercury EDM and the neutron EDM are milder. Fig. 11.2 displays the dependence. Taking into account the constraint from the mercury EDM, we find that the mercury EDM and the neutron EDM are almost independent of . On the other hand, the dependence of the electron EDM on is mild but visible. From these figures, one can see that the neutron and the mercury EDMs are sensitive to and also weakly depend on . On the other hand, the electron EDM is sensitive to + and weakly depend on . Most of the parameter space in Figs. 11.1 and 11.2 are within the future prospects of the electron EDM and the neutron EDM. In Summer of 2018, the constraint on the electron EDM is updated to be |de/e| <1.1 x 10-29e cm by ACME collaboration[13]. With this new constraint, the allowed regions in Figs. 11.1 and 11.2 become very thin stripes. Thus the correlation among the EDMs in future experiments provide a strong hint to explore the CP phases in the SUSY breaking sector. Let us discuss DM-nucleon scattering cross section. Since we consider the Higgs funnel scenario, the DM neutralino couples to neutral scalar bosons. Through these couplings, the DM neutralino and nucleon interact with each other. Though the couplings are rather small in the scenario, the couplings lead to a significant size of the spin-independent cross section and it will be within future prospects of the DM direct detection experiments. In Figure 11.3, the and ^ dependence of ctsi is shown. In this figure, the parameter choice is the same as in Fig. 11.1. Figure 11.4 displays the and ^ dependence of o$i with the same parameter choice as in Fig. 11.2. The spin-independent cross section is found to be smaller than the current upper bound [14-16] in all the region of the parameter space. However is is within the future prospects of the DARWIN [17], the DarkSide-20k [18], and the LZ [19]. Note that the scattering cross section depends on + and the dependence is not significant. 252 T. Shindou Fig. 11.1. The EDMs for tan ( = 30 and = 0°. The left (right) panels are for Mi = 1 TeV (Mi = 2 TeV). The contours in the top, the center, and the bottom panels are those of the electron EDM, the mercury EDM, and the neutron EDM, respectively. The dashed lines show the negative values. The red and blue shaded regions are excluded by the electron EDM and the mercury EDM, respectively. The figures are taken from Ref. [1]. 11 Electric Dipole Moment and Dark Matter... 253 Fig. 11.2. The EDMs for Mi = 1 TeV and tan ( = 30. In the left (right) panels, <£m2 = 0° (30°). The shadings and contours are the same as in Fig. 11.1. The figures are taken from Ref. [1]. 254 T. Shindou Fig. 11.3. The DM-nucleon scattering cross sections for tan p = 30 and = 0°. The left (right) panel is for |Mi | = 1 TeV (2 TeV). The shadings are the same as in Fig. 11.1. The figures are taken from Ref. [1]. Fig. 11.4. The DM-nucleon scattering cross sections for Mi = 1 TeV and tan p = 30. The left (right) panel is for |^m2 | = 0° (30°). The shadings are the same as in Fig. 11.1. The figures are taken from Ref. [1]. 11.4 Summary In this talk, we have considered the MSSM with CP phases, and we have focused on a DM scenario where the Bino-like neutralino is a DM whose annihilation cross section is enhanced enough through heavy Higgs bosons exchange so that the observed relic abundance of the DM can be explained. In this benchmark scenario, we have examined several EDMs and the spin-independent cross section for DM scattering with nuclei. We have shown that the EDMs are very powerful tool to explore the parameter space in the MSSM with CP phases even when most of the SUSY particles are very heavy. 11 Electric Dipole Moment and Dark Matter... 255 Acknowledgements This work was supported by JSPS KAKENHI Grant Number 17H05408 and Ko- gakuin University Grant for the project research. References 1. T. Abe, N. Omoto, O. Seto and T. Shindou, Phys. Rev. D 98 (2018) no.7, 075029. 2. N. Nagata and S. Shirai, JHEP 1501, 029 (2015) [arXiv:1410.4549 [hep-ph]]. 3. M. Drees and M. M. Nojiri, Phys. Rev. D 47, 376 (1993) [hep-ph/9207234]. 4. H. Baer and M. Brhlik, Phys. Rev. D 53, 597 (1996) [hep-ph/9508321]. 5. H. Baer and M. Brhlik, Phys. Rev. D 57, 567 (1998) [hep-ph/9706509]. 6. V. D. Barger and C. Kao, Phys. Rev. D 57, 3131 (1998) [hep-ph/9704403]. 7. J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. B 510, 236 (2001) [hep-ph/0102098]. 8. G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) [hep-ph/9506380]. 9. B. C. Allanach et al., Comput. Phys. Commun. 180, 8 (2009) [arXiv:0801.0045 [hep-ph]]. 10. P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 594 (2016) A13 doi:10.1051/0004-6361/201525830 [arXiv:1502.01589 [astro-ph.CO]]. 11. D. Barducci, G. Belanger, J. Bernon, F. Boudjema, J. Da Silva, S. Kraml, U. Laa and A. Pukhov, Comput. Phys. Commun. 222, 327 (2018) [arXiv:1606.03834 [hep-ph]]. 12. J. S. Lee, M. Carena, J. Ellis, A. Pilaftsis and C. E. M. Wagner, Comput. Phys. Commun. 184,1220 (2013) [arXiv:1208.2212 [hep-ph]]. 13. V. Andreev et al. [ACME Collaboration], Nature 562 (2018) no.7727, 355. doi:10.1038/s41586-018-0599-8 14. D. S. Akerib et al. [LUX Collaboration], Phys. Rev. Lett. 118 (2017) no.2, 021303 doi:10.1103/PhysRevLett.118.021303 [arXiv:1608.07648 [astro-ph.CO]]. 15. E. Aprile et al. [XENON Collaboration], Phys. Rev. Lett. 119 (2017) no.18, 181301 doi:10.1103/PhysRevLett.119.181301 [arXiv:1705.06655 [astro-ph.CO]]. 16. X. Cui et al. [PandaX-II Collaboration], Phys. Rev. Lett. 119 (2017) no.18, 181302 doi:10.1103/PhysRevLett.119.181302 [arXiv:1708.06917 [astro-ph.CO]]. 17. J. Aalbers et al. [DARWIN Collaboration], JCAP 1611, 017 (2016) [arXiv:1606.07001 [astro-ph.IM]]. 18. C. E. Aalseth et al, Eur. Phys. J. Plus 133, no. 3,131 (2018) [arXiv:1707.08145 [physics.ins-det]]. 19. D. S. Akerib et al. [LUX-ZEPLIN Collaboration], arXiv:1802.06039 [astro-ph.IM].