Simplified dark matter models with loop effects in direct detection and the constraints from indirect detection and collider search

  • We reexamine the simplified dark matter (DM) models with fermionic DM particle and spin-0 mediator. The DM-nucleon scattering cross sections of these models are low-momentum suppressed at tree-level, but receive sizable loop-induced spin-independent contribution. We perform one-loop calculation for scalar-type and twist-2 DM-quark operators and complete two-loop calculation for scalar-type DM-gluon operator. By analyzing the loop-level contribution from new operators, we find that future direct detection experiments can be sensitive to a fraction of parameter space. The indirect detection and collider search also provide complementary constraints on these models.
  • 加载中
  • [1] C. Boehm, M. J. Dolan, C. McCabe et al, JCAP, 1405: 009 (2014), arXiv:1401.6458[hep-ph
    [2] S. Ipek, D. McKeen, and A. E. Nelson, Phys. Rev D, 90: 055021 (2014), arXiv:1404.3716[hep-ph
    [3] M. Drees and M. Nojiri, Phys. Rev. D, 48: 3483 (1993), arXiv:hep-ph/9307208[hep-ph
    [4] M. Freytsis and Z. Ligeti, Phys. Rev. D, 83: 115009 (2011), arXiv:1012.5317[hep-ph
    [5] U. Haisch and F. Kahlhoefer, JCAP, 1304: 050 (2013), arXiv:1302.4454[hep-ph
    [6] G. Arcadi, M. Lindner, F. S. Queiroz et al, JCAP, 1803: 042 (2018), arXiv:1711.02110[hep-ph
    [7] N. F. Bell, G. Busoni, and I. W. Sanderson, JCAP, 1808: 017 (2018), [Erratum: JCAP1901, no.01, E01(2019)], arXiv: 1803.01574[hep-ph]
    [8] T. Li, Phys. Lett. B, 782: 497 (2018), arXiv:1804.02120[hep-ph
    [9] T. Han, H. Liu, S. Mukhopadhyay et al, JHEP, 03: 080 (2019), arXiv:1810.04679[hep-ph
    [10] T. Abe et al (LHC Dark Matter Working Group), LHC Dark Matter Working Group, (2018), arXiv:1810.09420[hep-ex
    [11] K. Ghorbani and P. H. Ghorbani, (2018), arXiv:1812.04092[hep-ph
    [12] W. Chao, G.-J. Ding, X.-G. He et al, (2018), arXiv:1812.07829[hep-ph
    [13] K. A. Mohan, D. Sengupta, T. M. P. Tait et al, (2019), arXiv:1903.05650[hep-ph
    [14] T. Abe, M. Fujiwara, and J. Hisano, JHEP, 02: 028 (2019), arXiv:1810.01039[hep-ph
    [15] J. Abdallah et al, Phys. Dark Univ., 9-10: 8 (2015), arXiv:1506.03116[hep-ph doi: 10.1016/j.dark.2015.08.001
    [16] G. Busoni et al, (2016), arXiv:1603.04156[hep-ex
    [17] F. Ertas and F. Kahlhoefer, (2019), arXiv:1902.11070[hep-ph
    [18] L. Lopez-Honorez, T. Schwetz, and J. Zupan, Phys. Lett. B, 716: 179 (2012), arXiv:1203.2064[hepph
    [19] A. Beniwal, F. Rajec, C. Savage et al, Phys. Rev. D, 93: 115016 (2016), arXiv:1512.06458[hep-ph
    [20] S. Baek, P. Ko, and J. Li, Phys. Rev. D, 95: 075011 (2017), arXiv:1701.04131[hep-ph
    [21] P. Athron et al, Eur. Phys. J. C, 79: 38 (2019), arXiv:1808.10465[hep-ph
    [22] T. Abe and R. Sato, Phys. Rev. D, 99: 035012 (2019), arXiv:1901.02278[hep-ph
    [23] G. Belanger, F. Boudjema, A. Pukhov et al, Comput. Phys. Commun., 185: 960 (2014), arXiv:1305.0237[hep-ph doi: 10.1016/j.cpc.2013.10.016
    [24] A. L. Fitzpatrick, W. Haxton, E. Katz et al, JCAP, 1302: 004 (2013), arXiv:1203.3542[hep-ph
    [25] H. H. Patel, Comput. Phys. Commun., 197: 276 (2015), arXiv:1503.01469[hep-ph doi: 10.1016/j.cpc.2015.08.017
    [26] J. Pumplin, D. R. Stump, J. Huston et al, JHEP, 07: 012 (2002), arXiv:hep-ph/0201195[hep-ph
    [27] J. Hisano, K. Ishiwata, and N. Nagata, Phys. Rev. D, 82: 115007 (2010), arXiv:1007.2601[hep-ph
    [28] J. Hisano, R. Nagai, and N. Nagata, JHEP, 05: 037 (2015), arXiv:1502.02244[hep-ph
    [29] J. Hisano, K. Ishiwata, and N. Nagata, JHEP, 06: 097 (2015), arXiv:1504.00915[hep-ph
    [30] R. J. Hill and M. P. Solon, Phys. Rev. D, 91: 043504 (2015), arXiv:1401.3339[hep-ph
    [31] R. J. Hill and M. P. Solon, Phys. Rev. D, 91: 043505 (2015), arXiv:1409.8290[hep-ph
    [32] F. Bishara, J. Brod, B. Grinstein et al, (2017), arXiv:1708.02678[hep-ph
    [33] F. D’Eramo, B. J. Kavanagh, and P. Panci, JHEP, 08: 111 (2016), arXiv:1605.04917[hep-ph
    [34] J. Aebischer, J. Kumar, and D. M. Straub, Eur. Phys. J. C, 78: 1026 (2018), arXiv:1804.05033[hep-ph
    [35] C. Arina, E. Del Nobile, and P. Panci, Phys. Rev. Lett., 114: 011301 (2015), arXiv:1406.5542[hep-ph doi: 10.1103/PhysRevLett.114.011301
    [36] C. Balázs, J. Conrad, B. Farmer et al, Phys. Rev. D, 96: 083002 (2017), arXiv:1706.01505[astro-ph.HE
    [37] T. Li, JHEP, 01: 151 (2018), arXiv:1708.04534[hep-ph
    [38] G. Bélanger, F. Boudjema, A. Goudelis et al, Comput. Phys. Commun., 231: 173 (2018), arXiv:1801.03509[hep-ph doi: 10.1016/j.cpc.2018.04.027
    [39] K. M. Zurek, Phys. Rev. D, 79: 115002 (2009), arXiv:0811.4429[hep-ph
    [40] S. Profumo, K. Sigurdson, and L. Ubaldi, JCAP, 0912: 016 (2009), arXiv:0907.4374[hep-ph
    [41] P. A. R. Ade et al (Planck), Astron. Astrophys, 594: A13 (2016), arXiv:1502.01589[astro-ph.CO doi: 10.1051/0004-6361/201525830
    [42] M. Ackermann et al (Fermi-LAT), Phys. Rev. Lett., 115: 231301 (2015), arXiv:1503.02641[astroph.HE doi: 10.1103/PhysRevLett.115.231301
    [43] A. Albert et al (Fermi-LAT, DES), Astrophys. J., 834: 110 (2017), arXiv:1611.03184[astro-ph.HE doi: 10.3847/1538-4357/834/2/110
    [44] M. Cirelli, G. Corcella, A. Hektor et al, JCAP, 1103: 051 (2011), [Erratum: JCAP1210, E01(2012)], arXiv: 1012.4515[hep-ph]
    [45] F. James and M. Roos, Comput. Phys. Commun., 10: 343 (1975) doi: 10.1016/0010-4655(75)90039-9
    [46] M. Aaboud et al (ATLAS), JHEP, 01: 126 (2018), arXiv:1711.03301[hep-ex
    [47] A. M. Sirunyan et al (CMS), Phys. Rev. D, 97: 092005 (2018), arXiv:1712.02345[hep-ex
    [48] M. Aaboud et al (ATLAS), Eur. Phys. J. C, 78: 18 (2018), arXiv:1710.11412[hep-ex
    [49] M. Aaboud et al (ATLAS), JHEP, 06: 108 (2018), arXiv:1711.11520[hep-ex
    [50] A. M. Sirunyan et al (CMS), Eur. Phys. J. C, 77: 845 (2017), arXiv:1706.02581[hep-ex
    [51] M. Aaboud et al (ATLAS), ATLAS, (2019), arXiv:1903.01400[hep-ex
    [52] N. Anand, A. L. Fitzpatrick, and W. C. Haxton, Phys. Rev. C, 89: 065501 (2014), arXiv:1308.6288[hep-ph
    [53] E. Aprile et al (XENON), Phys. Rev. Lett., 119: 181301 (2017), arXiv:1705.06655[astro-ph.CO doi: 10.1103/PhysRevLett.119.181301
    [54] E. Aprile et al (XENON), Phys. Rev. Lett., 121: 111302 (2018), arXiv:1805.12562[astro-ph.CO doi: 10.1103/PhysRevLett.121.111302
  • 加载中


Get Citation
. doi: 10.1088/1674-1137/43/11/113102
.  doi: 10.1088/1674-1137/43/11/113102 shu
Received: 2019-07-08
Article Metric

Article Views(8)
PDF Downloads(0)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌,
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article


Simplified dark matter models with loop effects in direct detection and the constraints from indirect detection and collider search

    Corresponding author: Tong Li,
    Corresponding author: Peiwen Wu,
  • 1. School of Physics, Nankai University, Tianjin 300071, China
  • 2. School of Physics, KIAS, 85 Hoegiro, Seoul 02455, Republic of Korea

Abstract: We reexamine the simplified dark matter (DM) models with fermionic DM particle and spin-0 mediator. The DM-nucleon scattering cross sections of these models are low-momentum suppressed at tree-level, but receive sizable loop-induced spin-independent contribution. We perform one-loop calculation for scalar-type and twist-2 DM-quark operators and complete two-loop calculation for scalar-type DM-gluon operator. By analyzing the loop-level contribution from new operators, we find that future direct detection experiments can be sensitive to a fraction of parameter space. The indirect detection and collider search also provide complementary constraints on these models.


    1.   Introduction
    • Although the existence of Dark Matter (DM) has been established by substantial cosmological and astronomical observations, the microscopic nature of DM particles is still unknown. An appealing candidate of DM is the Weakly Interacting Massive Particle (WIMP) arising from various extensions of the Standard Model (SM). The experimental searches for WIMP consist of four main categories, i.e. the direct detection (DD) of possible scattering between DM and SM target materials, the indirect detection (ID) looking for signals of DM annihilation/decay products from the sky, the collider searches for signals from DM productions at high energy accelerators, and the gravitational and/or cosmological effects originating from the DM in the early and/or the current Universe.

      Among the aforementioned four categories, the DD experiments have achieved significantly improved sensitivity in the past two decades, but yield null results up to now and very stringent bounds on the WIMP-nucleon scattering cross section. A natural explanation of the absence of a confirmed DM signal is that the scattering rate is highly suppressed by the small typical value of transfer momentum of the process and/or the relative velocity between DM and nucleon. A simple but compelling scenario resulting in the suppressed rate at tree-level is that fermionic DM particles $ \chi $ scatter off the target nucleon N through a pseudo-scalar mediator in the t-channel scattering process [1, 2]. The corresponding tree-level DM-nucleon contact interaction can reduce down to a non-relativistic contact operator

      $ \bar{\chi}i\gamma_5\chi N i\gamma_5 N \to ({{s}}_\chi\cdot {{q}})({{s}}_N\cdot {{q}}), $


      in the non-relativistic limit. Here, $ {{s}}_\chi $ ($ {{s}}_N $) is the DM (target nucleon) spin, and the scattering exchange momentum $ {{q}} $ is only of order 10 MeV. As a result, this scenario leads to a momentum-suppressed spin-dependent (SD) scattering cross section and thus an undetectable signal rate.

      The suppression of tree-level scattering rate makes it appealing to further scrutinize the high-order effect from one-loop induced processes (see early discussions in e.g. [3, 4]) and the possibly detectable signal at the upgrades of DD experiments. Integrating out the one-loop diagrams can induce distinct scalar-type operators giving non-momentum-suppressed spin-independent (SI) scattering cross section. The enhancement of loop-level SI cross section by the squared total nucleon number in a nucleus competes with the loop suppression and may dominate the WIMP-nuclei cross section over the suppressed tree-level scattering. This one-loop effect in direct DM detection has been investigated in both simplified frameworks and UV complete models [2, 5-13]. Recent progresses on the pseudo-scalar mediator scenario go beyond the one-loop processes for scalar-type DM-quark operator $ m_q\bar{\chi}\chi \bar{q}q $ and include the dedicated contributions from two-loop scattering diagrams for scalar-type DM-gluon operator $ {\alpha_s\over \pi}\bar{\chi}\chi GG $ after integrating out both heavy quarks and the mediator. It has been shown in Ref. [14], in the framework of Two Higgs Doublet Model, that the full two-loop calculation deviates sizably from the result obtained by the conventional relation between the scalar-type current of heavy quarks and that of the gluon. This discrepancy is caused by the failure of the quark momentum expansion for heavy quarks and an ignored two-loop diagram for gluon emission when one utilized the relation for DM-gluon scattering.

      In this work we revisit the loop effect in the DD of simplified DM models including either a pseudo-scalar mediator or maximal CP violation. We consider the hypotheses with one spin-0 mediator only coupled to the SM quarks and fermionic DM particles [15, 16], giving momentum-suppressed WIMP-nuclei scattering cross sections at tree-level. The latest approaches of dedicated loop calculations are utilized in the high-order contributions to the cross sections, together with the estimate of the Renormalization Group Equations (RGE) running effects. We find that the differential event rates with respect to recoil energy can provide particular spectral shapes at different energies induced by loop correction. Besides the loop effect on the direct DM detection for pseudo-scalar mediated models also discussed in the recent works [14, 17], we in addition take into account the constraints from other DM detection categories, e.g. the DM relic abundance, the ID constraint in terms of gamma-ray emission, as well as the current status of collider search. These synergistic studies are regarded as a more complete improvement of the Refs. [14, 17].

      The paper is outlined as follows. In Sec. 2 we describe the simplified dark matter models. Then we give the effective DM-nucleon interactions at tree-level and the corresponding DM-nucleus scattering cross sections. In Sec. 3, we present the effective operators and Wilson coefficients at loop-level for the DM-nucleon cross section. The numerical results are given in Sec. 5. Our conclusions are drawn in Sec. 6. Some technical details for loop calculations are collected in the Appendix.

    2.   Simplified dark matter hypothesis
    • In this work, we focus on the simplified DM frameworks which consist of Majorana fermion DM $ \chi $ and a spin-0 mediator a coupled to $ \chi $ and the SM quarks with strength $ g_\chi $ and $ g_q $ respectively. We consider each of the following three scenarios at a time

      $ {\cal{L}}_{\rm D2} = -{g_{\chi}\over 2} a \bar{\chi} i\gamma_5 \chi -g_q {m_q\over v_0} a \bar{q} q, $


      $ {\cal{L}}_{\rm D3} = -{g_{\chi}\over 2} a \bar{\chi} \chi -g_q {m_q\over v_0} a \bar{q} i\gamma_5 q, $


      $ {\cal{L}}_{\rm D4} = -{g_{\chi}\over 2} a \bar{\chi} i\gamma_5 \chi -g_q {m_q\over v_0} a \bar{q} i\gamma_5 q. $


      Here the $ a\bar{q}q $ coupling is also scaled by the SM-like Yukawa coupling with $ v_0 = 246 $ GeV being the SM Higgs vacuum expectation value. The model D4 with a being a pure pseudo-scalar is designated as the pseudo-scalar mediator DM model in most of literatures. Models D2 and D3 are induced by more specific UV complete models with CP violation [17-22] and correspond to the cases with maximal CP violation.

      Based on the DM interactions with quarks and gluons at tree-level, the DM-nucleon contact interactions are described by the effective Lagrangians as follows

      $\begin{split}& {\cal{L}}^{\rm eff}_{\rm D2} = {C_N^{\rm tree}({\rm D2})\over 2m_{a}^2} \bar{\chi} i\gamma_5 \chi \bar{N} N, \\ &{\cal{L}}^{\rm eff}_{\rm D3} = {C_N^{\rm tree}({\rm D_3})\over 2m_a^2} \bar{\chi} \chi \bar{N} i\gamma_5 N, \\ &{\cal{L}}^{\rm eff}_{\rm D4} = {C_N^{\rm tree}({\rm D_4})\over 2m_a^2} \bar{\chi} i\gamma_5 \chi \bar{N} i\gamma_5 N, \end{split} $


      where the tree-level coefficients are defined as

      $ C_N^{\rm tree}({\rm D2}) = \sum\limits_{q = u,d,s}{m_N\over m_q} C_q f_q^{N} + \sum\limits_{q = c,b,t}{m_N\over m_q} C_q {2\over 27} f_G^{N}, $


      $ \begin{split} C_N^{\rm tree}({\rm D4}) =& \sum\limits_{q = u,d,s}{m_N\over m_q} \left(C_q-C\right) \Delta_q^{N}, \\ C =& \bar{m}\sum\limits_{q = u,\cdots,t}{C_q\over m_q}, \\ \bar{m}^{-1} =& \sum\limits_{q = u,d,s}m_q^{-1}. \end{split}$


      The coefficient $ C_N^{\rm tree}({\rm D_3}) $ is equal to $ C_N^{\rm tree}({\rm D_4}) $ as models D3 and D4 share the same quark bilinear form $ \bar{q}\gamma_5 q $. Here $ \Delta_q^{N} $, $ f_q^{N} $ and $ f_G^{N} $ are quark/gluon-nucleon form factors as numerically used in micrOMEGAs [23]. The quark-level constant is defined as $ C_q = g_\chi g_q{m_q\over v_0} $. Consequently, their differential DM-nucleus scattering cross sections read as

      $\begin{split} {{\rm d}\sigma_{\rm SI}({\rm D2})\over {\rm d}E_R} =& {1\over 32\pi} {m_T\over m_\chi^2 m_N^2 v^2} {4m_N^2q^2\over m_a^4}\\&\times \sum\limits_{N,N' = p,n}C_N^{\rm tree}({\rm D_2}) C_{N'}^{\rm tree}({\rm D_2}) F_{M}^{(N,N')}(q^2) , \end{split}$


      $\begin{split} {{\rm d}\sigma_{\rm SD}({\rm D3})\over {\rm d}E_R} =& {1\over 32\pi} {m_T\over m_\chi^2 m_N^2 v^2} {4m_\chi^2 q^2\over m_a^4} \\&\times\sum\limits_{N,N' = p,n}C_N^{\rm tree}({\rm D_3}) C_{N'}^{\rm tree}({\rm D_3}) F_{\Sigma''}^{(N,N')}(q^2) , \end{split}$


      $\begin{split} {{\rm d}\sigma_{\rm SD}({\rm D4})\over {\rm d}E_R} =& {1\over 32\pi} {m_T\over m_\chi^2 m_N^2 v^2} {q^4\over m_a^4}\\&\times \sum\limits_{N,N' = p,n}C_N^{\rm tree}({\rm D_4}) C_{N'}^{\rm tree}({\rm D_4}) F_{\Sigma''}^{(N,N')}(q^2) , \end{split}$


      where $ m_T $ is the nucleus mass, v is the DM speed in Earth's frame, $ E_R $ is the nuclear recoil energy and $ F_{\Sigma''}^{(N,N')}(q^2), F_{M}^{(N,N')}(q^2) $ are the form factors defined in Ref. [24]. The tree-level WIMP-nucleus scattering cross sections of the above simplified models are all dependent on the transfer momentum $ q = \sqrt{2E_R m_T} $. As seen above, they are suppressed by $ m_N^2q^2/ m_a^4 $, $ m_\chi^2 q^2/ m_a^4 $ and $ q^4/m_a^4 $ for models D2, D3 and D4, respectively.

    3.   Loop effect in direct detection
    • In this section, we derive the loop-level effect in direct DM detection of the above simplified models, followed by the estimate of the scale effects in terms of renormalization group evolutions.

    • 3.1.   Loop effect from scalar-type quark/gluon operators

    • The general Lagrangian for the non-momentum-suppressed DM-nucleon SI cross section is given by

      $\begin{split} {\cal{L}}_{\rm eff} =& {1\over 2}\sum_{q = u,d,s}C'_q m_q \bar{\chi}\chi \bar{q}q + {1\over 2}C_G\left(-{9\alpha_s\over 8\pi}\bar{\chi}\chi G^a_{\mu\nu}G^{a\mu\nu}\right) \\& + {1\over 2}\sum_{q = u,d,s,c,b}\left[C_q^{(1)}\bar{\chi}i\partial^\mu \gamma^\nu \chi {\cal{O}}^q_{\mu\nu}+C_q^{(2)}\bar{\chi}i\partial^\mu i\partial^\nu \chi {\cal{O}}^q_{\mu\nu}\right], \end{split}$


      where $ {\cal{O}}^q_{\mu\nu} = {i\over 2}\bar{q}\left(\partial_\mu\gamma_\nu+\partial_\nu\gamma_\mu-{1\over 2}g_{\mu\nu}{\not\!\! \partial }\right)q $ is the twist-2 operator. For the models we consider, the Wilson coefficients in Eq. (11) are all zero at tree-level but can be generated at loop-level, denoted by $ C'_q = C_q^{\rm box} $, $ C_G = C_G^{\rm box} $, $ C_q^{(1)} = C_q^{\rm (1)box} $, $ C_q^{(2)} = C_q^{\rm (2)box} $. The coefficients for scalar-type DM-quark operator and the twist-2 operator, i.e. $ C_q^{\rm box} $, $ C_q^{\rm (1)box} $ and $ C_q^{\rm (2)box} $, are generated by the box diagrams in the top panels of Fig. 1. The two-loop diagrams in Fig. 1 with only heavy quark Q in the quark loop contribute to the scalar-type DM-gluon operator and the coefficient $ C_G^{\rm box} $.

      Figure 1.  Loop diagrams for the DM-quark currents (top) and DM-gluon current (bottom).

      Following the non-relativistic limit used in Ref. [14], we expand the small momentum of valence quarks in the amplitude of the DM-quark scattering box diagrams. The coefficients $ C_q^{\rm box} $, $ C_q^{\rm (1)box} $ and $ C_q^{\rm (2)box} $ are then obtained by reading out the DM-quark effective operators. For the DM-gluon coefficient $ C_G^{\rm box} $, one needs to calculate the amplitude of two-loop diagrams and find the effective operator $ \bar{\chi}\chi G^a_{\mu\nu}G^{a\mu\nu} $. The complete two-loop calculations ensure the validity of the obtained $ C_G^{\rm box} $ for any values of mediator mass $ m_a $. For model D4, the above Wilson coefficients are equivalent to those for pseudo-scalar mediator model as derived in Ref. [14]

      $ \begin{split} C_q^{\rm box}({\rm D4}) =& {-m_\chi\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^2} [6X_{001}(m_\chi^2,m_\chi^2,0,m_a^2)\\&+m_\chi^2X_{111}(m_\chi^2,m_\chi^2,0,m_a^2) -6X_{001}(m_\chi^2,m_\chi^2,m_a^2,0)\\&-m_\chi^2X_{111}(m_\chi^2,m_\chi^2,m_a^2,0)], \end{split} $


      $\begin{split} C_q^{\rm (1)box}({\rm D4}) =& {-8\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^2} [X_{001}(m_\chi^2,m_\chi^2,0,m_a^2)\\&-X_{001}(m_\chi^2,m_\chi^2,m_a^2,0)], \end{split}$


      $ \begin{split} C_q^{\rm (2)box}({\rm D4}) =& {-4m_\chi\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^2} [X_{111}(m_\chi^2,m_\chi^2,0,m_a^2)\\&-X_{111}(m_\chi^2,m_\chi^2,m_a^2,0)], \end{split}$


      $ C_G^{\rm box}({\rm D4}) = \sum_{Q = c,b,t}{-m_\chi\over 432\pi^2}\left({m_Q\over v_0}\right)^2 g_\chi^2 g_q^2 {\partial F(m_a^2)\over \partial m_a^2}, $


      where the loop functions $ X_{001}, X_{111}, F $ are given in Ref. [14] and the references therein. Following the same procedure, for models D2 and D3, we obtain the corresponding Wilson coefficients which are related to those in model D4

      $\begin{split} C_q^{\rm box}({\rm D2}) =& C_q^{\rm box}({\rm D4}) +{4m_\chi\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^2}\bigg[C_2(m_\chi^2,m_a^2,m_\chi^2)\\&+{1\over m_a^2}B_1(m_\chi^2,0,m_\chi^2)-{1\over m_a^2}B_1(m_\chi^2,m_a^2,m_\chi^2)\bigg], \end{split} $


      $ \begin{split} C_q^{\rm box}({\rm D3}) =& C_q^{\rm box}({\rm D4})+{-8m_\chi\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^4}\bigg[X_{00}(m_\chi^2,m_a^2,m_\chi^2)\\&+{m_\chi^2\over 4}X_{11}(m_\chi^2,m_a^2,m_\chi^2)\bigg], \end{split} $


      $ \begin{split} C_q^{\rm (1)box}({\rm D2}) = &C_q^{\rm (1)box}({\rm D3}) = C_q^{\rm (1)box}({\rm D4}), \\ C_q^{\rm (2)box}({\rm D2}) =& C_q^{\rm (2)box}({\rm D4}), \end{split}$


      $ C_q^{\rm (2)box}({\rm D3}) \!=\! C_q^{\rm (2)box}({\rm D4})\!+\!{-8m_\chi\over (4\pi)^2}\left({m_q\over v_0}\right)^2 {g_\chi^2 g_q^2\over m_a^4}X_{11}(m_\chi^2,m_a^2,m_\chi^2), $


      $ C_G^{\rm box}({\rm D2}) = \sum_{Q = c,b,t}{-m_\chi\over 432\pi^2}\left({m_Q\over v_0}\right)^2 g_\chi^2 g_q^2 {\partial F'(m_a^2)\over \partial m_a^2}, $


      $ C_G^{\rm box}({\rm D3}) = \sum_{Q = c,b,t}{-m_\chi\over 432\pi^2}\left({m_Q\over v_0}\right)^2 g_\chi^2 g_q^2 {\partial F''(m_a^2)\over \partial m_a^2}. $


      The new loop functions here are collected in Appendix. In our numerical calculation, we use Package-X [25] to compute the above loop functions.

      Based on the above effective operators for SI DM-nucleon scattering and the corresponding Wilson coefficients, we define the DM-nucleon constant at loop-level

      $ \begin{split} C_N^{\rm loop} =& m_N\bigg[\sum_{q = u,d,s}C_q' f_{q}^N + C_G f_{G}^N \\&+ {3\over 4}\sum_{q = u,d,s,c,b}\left(m_\chi C_q^{(1)}+m_\chi^2C_q^{(2)}\right)\left(q^N(2)+\bar{q}^N(2)\right)\bigg], \end{split} $


      where the second moments of the parton distribution functions for quarks $ q^N(2) $ and anti-quarks $ \bar{q}^N(2) $ are taken from the CTEQ PDFs [26]. The SI cross section of the DM interaction with nucleon is thus given by

      $ \sigma_{\rm SI} = {1\over \pi}\left({m_\chi m_N\over m_\chi + m_N}\right)^2 |C_N^{\rm loop}|^2. $


      In terms of the form factor function $ F_{M}^{(N,N')}(q^2) $, the differential SI cross section of the DM interaction with nucleus with mass $ m_T $ is

      $ {{\rm d}\sigma_{\rm SI}\over {\rm d}E_R} = {1\over 2\pi} {m_T\over v^2} \sum_{N,N' = p,n}C_N^{\rm loop} C_{N'}^{\rm loop} F_{M}^{(N,N')}(q^2). $

    • 3.2.   Loop effect from RGE running

    • Another manifestation of loop effect is the mixing of operators according to the RGE. This RGE effect can be important if one considers DM phenomenology at vastly different energy scales. For instance, the DM annihilation typically happens at the electroweak scale and the DM particles are possibly produced near TeV scale at colliders. The energy scale for the DM-nucleon scattering in DD experiments is of the order of the hadron scale $ \mu_{\rm had} $. One usually starts with a gauge-invariant renormalizable DM model defined near or above the electroweak (EW) scale $ \mu_{\rm EW}\simeq m_Z $, but studies the non-relativistic DM-nucleon scattering rate with characteristic scale $ \mu_{\rm had}\simeq 1 $ GeV. A series of effective field theories (EFTs) should be properly constructed by integrating out particles heavier than the current EFT scale $ \mu_{\rm EFT} $ and reasonably matched when passing thresholds of particles lighter than $ \mu_{\rm EFT} $ where they are integrated out in a similar way. Between the thresholds, the evolutions and mixings of the EFT operators should be performed according to RGE. The scale of the first EFT constructed in the whole analysis determines the procedures of RGE and threshold matching. The above procedures have been well elaborated in e.g. [27-31] and implemented in packages such as DirectDM [32], runDM [33], Wilson [34] for specific or generic models.

      An important difference should be emphasized between the scalar-type and twist-2 operators in Eq. (11). The scalar-type form factors $ f_q^{N} $ for light quarks $ q = u,d,s $ are attributed to the non-perturbative QCD effects with energy scale around or below 1 GeV, and are obtained from the lattice QCD simulations. Thus, the scalar-type Wilson coefficients of light quarks $ u,d,s $ and gluon must take values around 1 GeV. Depending on the scale of the first EFT constructed in the whole analysis, e.g. at $ \mu = m_Z $, this implies the procedures of RGE and threshold matching when calculating the scalar-type operator contributions. At the scale of about 1 GeV, the heavy quarks $ Q = c,b,t $ have been integrated out into the scalar-type gluon operator using the full two-loop calculations, as emphasized in Refs. [14, 17]. In turn, the twist-2 form factors $ q^N(2), \bar{q}^N(2) $ in Eq. (11) can be calculated perturbatively using parton PDFs at various scales [3, 27-31], e.g. 1 GeV or $ m_Z $. One can choose a convenient scale to calculate the twist-2 contributions, with the proper active field contents (e.g. 5 flavor quarks $ u,d,s,c,b $ at $ \mu = m_Z $) and the Wilson coefficients and form factors evaluated at that scale. Note that we ignored the negligible contributions from twist-2 gluon operator, since its Wilson coefficient is suppressed by an additional $ \alpha_s/\pi $ due to the operator definition [28, 29].

      We note that in Refs. [14, 17], the values of $ q^N(2), \bar{q}^N(2) $ are evaluated at $ \mu = m_Z $ in the calculations. Since the same (similar) box diagrams and model parameters are used to obtain the scalar-type and twist-2 Wilson coefficients for quarks (gluons), a more consistent implementation should involve the RGE running effects for the scale-type operators from $ m_Z $ to 1 GeV as discussed above. Since the coupling between the mediator and the SM quarks is chosen to mimic the SM Yukawa structure, the DM-nucleon constant from scalar-type DM-gluon interaction $ C_G f_{G}^N $ in Eq. (22) is dominated by the top quark loop. The constant $ C_G f_{G}^N $ also dominates over the scalar-type and twist-2 DM-quark interactions. To have a conservative estimate of the scale effects on the scalar type DM-gluon operator $ -{9\alpha_s\over 8\pi}\bar{\chi}\chi G^a_{\mu\nu}G^{a\mu\nu} $, we utilize the package DirectDM to perform its RGE running from $ m_Z $ to 1 GeV. We find that the scale effects give a negative correction of 1%~2%, and thus do not affect our main conclusions in this work.

    4.   Other dark matter constraints
    • In this section, we consider DM constraints on the above simplified models from other categories mentioned in the Introduction, including the relic abundance, indirect detection and collider search.

      Assuming that DM particles have frozen out in the early Universe, as standard thermal relics, they acquire their present abundance through annihilation processes. The pairs of DM particle $ \chi $ in the simplified models can either annihilate into SM quark or gluon pairs via s-channel processes $ \chi \chi\to a\to q\bar{q}, gg $ or annihilate into two mediators $ \chi \chi\to a a $ when kinematically allowed [15, 16, 35-37]. The amplitudes of the two annihilation channels are governed by $ g_\chi g_q $ and $ g_\chi^2 $, respectively. We assume all kinematically accessible final states of the DM annihilation and use micrOMEGAs 5.0 [38] to calculate the relic abundance. Note that the WIMP candidate here may account only a fraction of the total DM of the Universe, referred as multi-component DM scenario [39, 40]. In this scenario, the DM energy density measured by PLANCK [41] is imposed as an upper limit on the WIMP relic abundance.

      Dwarf galaxies are the search targets for DM annihilation into gamma rays. The Fermi Large Area Telescope (LAT) has detected no excess of gamma ray emission from the dwarf spheroidal satellite galaxies (dSphs) of the Milky Way. Thus, an upper limit on the DM annihilation cross section can be placed from a combined analysis of multiple Milky Way dSphs [42, 43]. For individual dwarf galaxy target, Fermi-LAT provided tabulated values of delta-log-likelihood as a function of the energy flux bin-by-bin. The gamma ray energy flux from DM annihilation for the jth energy bin and the kth dwarf is given by

      $ \Phi^E_{j,k}(m_{\chi},\langle \sigma v\rangle,J_k) = \frac{\langle \sigma v\rangle}{16\pi m_{\chi}^2}J_k\int^{E^{\rm max}_j}_{E^{\rm min}_j}E\frac{{\rm d}N_\gamma}{{\rm d}E}{\rm d}E, $


      where $ J_k $ is the J factor for the kth dwarf. The energy flux only depends on $ m_{\chi} $, $ \langle \sigma v\rangle $ and $ J_k $, and is thus calculable for the DM annihilation process from the above simplified models. We use the PPPC4DMID package [44] to obtain the spectrum of photons $ {\rm d}N_\gamma/{\rm d}E $. The likelihood for kth dwarf is

      $ {\cal{L}}_k(m_{\chi},\langle \sigma v\rangle,J_k) = {\cal{L}}_J(J_k|\bar{J}_k,\sigma_k)\prod_j {\cal{L}}_{j,k}(\Phi^E_{j,k}(m_{\chi},\langle \sigma v\rangle,J_k)), $


      where $ {\cal{L}}_{j,k} $ is the tabulated likelihood provided by Fermi-LAT for each dwarf and energy flux. The uncertainty of the J factors is taken into account by profiling over $ J_k $ in the likelihood below [42]

      $ {\cal{L}}_J(J_k|\bar{J}_k,\sigma_k) = {1\over \ln(10)J_k\sqrt{2\pi}\sigma_k}\times e^{-(\log_{10}(J_k)-\log_{10}(\bar{J}_k))^2/2\sigma_k^2}, $


      with the measured J factor $ \bar{J}_k $ and error $ \sigma_k $. A joint likelihood for all dwarfs can then be performed as

      $ {\cal{L}}(m_{\chi},\langle \sigma v\rangle,\mathbb{J}) = \prod_k {\cal{L}}_k(m_{\chi},\langle \sigma v\rangle,J_k), $


      where $ \mathbb{J} $ is the set of J factors $ J_k $. In our numerical implementation, we adopt the corresponding values of $ {\cal{L}}_{j,k} $ and $ \bar{J}_k, \sigma_k $ for 19 dwarf galaxies considered in Ref. [43].

      According to the maximum likelihood analysis adopted by Fermi-LAT, the delta-log-likelihood is given by

      $ -2\Delta \ln {\cal{L}}(m_{\chi},\langle \sigma v\rangle) = -2\ln\left({{\cal{L}}(m_{\chi},\langle \sigma v\rangle,\widehat{\widehat{\mathbb{J}}})\over {\cal{L}}(m_{\chi},\widehat{\langle \sigma v\rangle},\widehat{\mathbb{J}})}\right), $


      where $ \widehat{\langle \sigma v\rangle} $ and $ \widehat{\mathbb{J}} $ maximize the likelihood at any given $ m_{\chi} $, and $ \widehat{\widehat{\mathbb{J}}} $ maximizes the likelihood for given $ m_{\chi} $ and $ \langle \sigma v\rangle $. The 95% C.L. upper limit on the annihilation cross section for a given $ m_{\chi} $ is determined by demanding $ -2\Delta \ln{\cal{L}}(m_{\chi},\langle \sigma v\rangle)\leqslant 2.71 $. We perform the likelihood analysis and obtain the upper limit using Minuit [45]. Once the annihilation cross section calculated by a certain set of model parameters is larger than the limit, we claim the corresponding parameter values are excluded by Fermi-LAT dSphs.

      The Large Hadron Collider (LHC) performed the search for DM in association with energetic jet [46, 47] or the third generation quarks [48-50] for simplified DM models with spin-0 mediator at $ \sqrt{s} = 13 $ TeV collisions. The most severe limits are from final states with $ t\bar{t} $ and missing transverse momentum [51]. For model D4 with pseudo-scalar mediator, assuming unitary couplings $ g_\chi = g_q = 1 $, the range of mediator mass between 15 and 25 GeV is excluded [49]. This limit is valid for all DM masses as long as $ m_a>2m_\chi $ and closely related models D2 and D3 should have very similar collider constraint.

    5.   Results
    • By combining the theoretical calculations of DM-nucleus scattering cross sections with a certain velocity distribution for DM particles, we can calculate WIMP signal rates for DD experiments. The differential event rate with respect to the recoil energy is given by

      $ {{\rm d}N\over {\rm d}E_R} = {\rho_\chi\over m_\chi}\int {\rm d}^3v vf(\vec{v}){{\rm d}\sigma(v,E_R)\over {\rm d}E_R}, $


      where $ \rho_\chi $ is the local DM density which is fixed to be $ 0.3 \ {\rm GeV}/{\rm cm}^3 $. We take a cut-off Maxwell-Boltzmann velocity distribution for $ f(\vec{v}) $ with the escape velocity as $ v_{\rm esc} = 544 \ {\rm km}/{\rm s} $. Together with the differential scattering cross sections obtained above, the differential event rate can be evaluated and in practice we employ DMFormFactor [24, 52] for the numerical calculation on xenon nucleus $ ^{129}{\rm Xe} $. The recoil energy spectra of models D2, D3 and D4 are shown in Figs. 2, 3 and 4, respectively. Model D2 leads to SI scattering cross section with a strong enhancement for large nuclei. Thus, although there is momentum suppression at tree-level, the spectrum at tree-level for model D2 still dominates over the loop-level contribution. As being suppressed by $ q^4 $ at tree-level, in turn, the loop-level spectrum is much greater than that at tree-level for model D4. Although models D3 and D4 both lead to SD cross section at tree-level, this discrepancy is smaller in model D3 as its tree-level scattering cross section is suppressed by $ q^2 $ only. The loop-level spectrum of model D3 dominates over the tree-level one only in the range of low recoil energy.

      Figure 2.  (color online) Recoil energy spectra of model D2 with $ m_a = 10 $ GeV, $ g_\chi g_q = 1 $ and $ m_\chi = 40 $ GeV (left) or $ m_\chi = 400 $ GeV (right). The tree-level and loop-level spectra are denoted by red and black curves, respectively

      Figure 3.  (color online) Recoil energy spectra of model D3, as labeled in Fig. 2.

      Figure 4.  (color online) Recoil energy spectra of model D4, as labeled in Fig. 2.

      As the recoil spectrum of the SI scattering induced by loop diagrams is dominant in model D4, the prediction of SI DM-nucleon cross section in Eq. (23) can be compared directly to the limits set by DD experiments to yield a bound on $ g_\chi g_q $. As shown in Fig. 5, for $ g_\chi g_q \leqslant 1 $ in model D4, the SI scattering cross sections in the red region are below the XENON1T exclusion limit [53, 54] but above the neutrino floor. The green region gives cross sections below the neutrino floor. Future DD experiments can thus be more sensitive to the region of $ g_\chi g_q > 0.4 $ and $ m_\chi < 200 $ GeV in the case of $ m_a = 10 $ GeV.

      Figure 5.  (color online) The region of $ g_\chi g_q $ vs. $ m_\chi $ below the Xenon1T exclusion limit and above the neutrino floor (red) and the region below the neutrino floor (green) for model D4, with $ m_a = 10 $ GeV (left) and $ m_a = 30 $ GeV (right). Assuming $ g_\chi = 1 $, the light blue values are excluded by Fermi-LAT dSphs and the DM relic abundance satisfies $ \Omega h^2\leqslant 0.12 $ above the black curves.

      The cross section of DM annihilation into SM quarks for model D4 is proportional to $ m_q^2/v_0^2 $, thus the $ t\bar{t} $ channel dominates if kinematically allowed. The annihilation to mediator pairs is governed by $ m_a/m_\chi $ and plays a crucial role in small $ m_\chi $ range. The Fermi-LAT dSphs exclude a majority of parameter region for $ m_{\chi}\lesssim 100 $ GeV and $ m_{\chi}\gtrsim m_t $, as shown in light blue region in Fig. 5. The region around $ m_{\chi}\simeq 100 $ GeV evades the ID constraint due to the fact that the $ \chi\chi\to a\to t\bar{t} $ channel is not kinematically allowed and the annihilation to mediator pairs is suppressed [36, 37]. Above the black contours in Fig. 5, the DM relic abundance satisfies $ \Omega h^2\leqslant 0.12 $. In the multi-component DM scenario, other DM particles are required to explain the measured value of the thermal DM energy density in the region above the black curves. Then the dwarf spheroidal galaxies contain not only $ \chi $ but also others. The number density of $ \chi $ is thus smaller than the value estimated from the J-factors. We expect the ID constraint on $ \chi $ in the multi-component scenario is less stringent than the exclusion shown in Fig. 5. The detailed study of the multi-component DM scenario is beyond the scope of this work.

    6.   Conclusion
    • We reexamined the loop-level correction to the WIMP-nucleon scattering cross section in the framework of simplified DM models. A spin-0 mediator is assumed to couple with fermionic DM particle and the SM quarks in each model. The simplified models include either pseudo-scalar mediator or maximal CP violation. The cross sections of these models are low-momentum suppressed at tree-level, but receive sizable loop-induced SI contribution. Following the recent progress on the loop-level correction, we perform one-loop calculation for scalar-type and twist-2 DM-quark operators and complete two-loop calculation for scalar-type DM-gluon operator. The differential event rates with respect to recoil energy can provide particular spectral shapes induced by loop correction for models with pseudo-scalar mediator or CP violation. By including the loop-level SI cross section, we find that future DD experiments can be sensitive to a fraction of parameter space with $ g_\chi g_q > 0.4 \ (0.8) $ and $ m_\chi < 200 \ (80) $ GeV for the pseudo-scalar mass being 10 (30) GeV which gives no detectable signal with only tree-level contribution. The sensitivity of DD experiments to these models is also complementary to the constraints from ID and collider search.

      T. L. would like to thank Tomohiro Abe, Motoko Fujiwara and Junji Hisano for helpful discussion. P. W. would like to thank Chengcheng Han for the beneficial argument.

    Appendix A: Loop functions
    • The used loop functions for models D2 and ${\rm{D3}} $ are

      $\tag{A1} \begin{split} \int {{\rm d}^D \ell \over (2\pi)^D} {\ell_\mu \ell_\nu \over [(\ell+p)^2-m_\chi^2](\ell^2-m_a^2)^2\ell^4} = &{i\over (4\pi)^2 m_a^4} [ g_{\mu\nu} X_{00}(p^2,m_a^2,m_\chi^2)\\&+p_\mu p_\nu X_{11}(p^2,m_a^2,m_\chi^2)] \end{split} $



      $ \tag{A2}\begin{split} F'(m_a^2) =& \int^1_0 {\rm d}x [3Y_1(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)-m_Q^2{-x^2-3x\over x^2(1-x)^2}Y_2(m_\chi^2,m_\chi^2,m_a^2,m_Q^2) \\& +4m_Q^4{x^2(1-2x)\over x^3(1-x)^3}Y_3(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)], \end{split} $


      $\begin{split} F''(m_a^2) =& \int^1_0 {\rm d}x [3Y_1(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)-m_Q^2{9x-5x^2\over x^2(1-x)^2}Y_2(m_\chi^2,m_\chi^2,m_a^2,m_Q^2) \\& -2m_Q^4{2x^2\over x^3(1-x)^3}Y_3(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)] \\& +2\int^1_0 {\rm d}x [3Z_1(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)-m_Q^2{9x-5x^2\over x^2(1-x)^2}Z_2(m_\chi^2,m_\chi^2,m_a^2,m_Q^2) \end{split}$

      $ \tag{A3}\begin{split} -2m_Q^4{2x^2\over x^3(1-x)^3}Z_3(m_\chi^2,m_\chi^2,m_a^2,m_Q^2)], \end{split}$



      $\tag{A4} \int {{\rm d}^D\ell\over (2\pi)^D}{1\over [(\ell+p)^2-m_\chi^2](\ell^2-m_a^2)\left[\ell^2-\displaystyle{m_Q^2\over x(1-x)}\right]} = {i\over (4\pi)^2}Z_1(p^2,m_\chi^2,m_a^2,m_Q^2), $


      $\tag{A5} \int {{\rm d}^D\ell\over (2\pi)^D}{1\over [(\ell+p)^2-m_\chi^2](\ell^2-m_a^2)\left[\ell^2-\displaystyle{m_Q^2\over x(1-x)}\right]^2} = {i\over (4\pi)^2}Z_2(p^2,m_\chi^2,m_a^2,m_Q^2), $


      $\tag{A6} \int {{\rm d}^D\ell\over (2\pi)^D}{1\over [(\ell+p)^2-m_\chi^2](\ell^2-m_a^2)\left[\ell^2-\displaystyle{m_Q^2\over x(1-x)}\right]^3} = {i\over (4\pi)^2}Z_3(p^2,m_\chi^2,m_a^2,m_Q^2). $


      Here the $ Y_i (i = 1,2,3) $ functions are defined in Ref. [14] and the references therein.

Reference (54)



DownLoad:  Full-Size Img  PowerPoint