INT-PUB-17-031 Improved limits for Higgs-portal dark matter from LHC searches Martin Hoferichter,1 Philipp Klos,2, 3 Javier Men´endez,4 and Achim Schwenk2, 3, 5 1Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA 2Institut f¨ur Kernphysik, Technische Universit¨at Darmstadt, 64289 Darmstadt, Germany 3ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum f¨ur Schwerionenforschung GmbH, 64291 Darmstadt, Germany 4Center for Nuclear Study, The University of Tokyo, 113-0033 Tokyo, Japan 5Max-Planck-Institut f¨ur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany Searches for invisible Higgs decays at the Large Hadron Collider constrain dark matter Higgs- portal models, where dark matter interacts with the Standard Model fields via the Higgs boson. While these searches complement dark matter direct-detection experiments, a comparison of the two limits depends on the coupling of the Higgs boson to the nucleons forming the direct-detection nuclear target, typically parameterized in a single quantity fN . We evaluate fN using recent phe- nomenological and lattice-QCD calculations, and include for the first time the coupling of the Higgs boson to two nucleons via pion-exchange currents. We observe a partial cancellation for Higgs-portal models that makes the two-nucleon contribution anomalously small. Our results, summarized as fN = 0.308(18), show that the uncertainty of the Higgs–nucleon coupling has been vastly overesti- mated in the past. The improved limits highlight that state-of-the-art nuclear physics input is key to fully exploiting experimental searches. PACS numbers: Dark matter, WIMPs, chiral Lagrangians Keywords: 95.35.+d, 14.80.Ly, 12.39.Fe INTRODUCTION matrix elements for the scalar current N mN fq = hN|mqqq¯ |Ni (1) After the initial discovery of a new particle by the AT- LAS and CMS collaborations at the Large Hadron Col- for the light quarks q = u,d,s (with quark masses mq lider (LHC) in 2012 [1, 2], its mass has been further con- and nucleon mass mN ), as well as through heavy-quark strained to sub-GeV accuracy, mh = 125.09(24) GeV [3], loops coupling to gluon fields, see Fig. 1. Integrating out and so far neither studies of its spin and parity [4, 5] nor the heavy quarks at leading order in the strong coupling its branching fractions [6] have revealed significant devi- constant αs leads to an effective coupling [19] ations from the Standard Model (SM) prediction for the 2 7 Higgs boson h. An interesting consequence is that events f = f N = + f N (2) N q 9 9 q with large missing transverse momentum would be ex- q=u,d,s,c,b,tX q=Xu,d,s pected if the Higgs boson decayed to long-lived non-SM states that do not leave a signature in the detector, so- to describe the Higgs–nucleon interaction. Both re- called invisible Higgs decays. The absence of such ob- cent LHC analyses constraining Higgs-portal dark mat- servations sets limits on the invisible decay channels [7– ter from invisible Higgs decays [12, 13] use a central 13], which provide stringent constraints on beyond-SM value fN = 0.326 [20] with variations covering fN = physics. 0.260 ... 0.629 (motivated by [21], see [15]; see also [22]). Given recent progress in lattice-QCD calculations [23– One particular example concerns Higgs-portal models arXiv:1708.02245v2 [hep-ph] 20 Oct 2017 26] and pion–nucleon (πN) phenomenology [27–30], this for dark matter, see, e.g., [14–17], in which the dark- large range in fN no longer reflects the current knowl- matter candidate χ, a weakly interacting massive particle edge of the scalar couplings of the nucleon. Consequently, (WIMP), interacts with the SM fields via exchange of the the limits for σχN derived from the searches for invisible Higgs boson. With the nature of dark matter still a major Higgs decays can be significantly improved. In the first puzzle, the quest to understand its composition is being part of this paper, we provide a more detailed assessment vigorously pursued, besides collider signatures, in direct- of the current situation, including corrections to Eq. (2) detection experiments looking for the WIMP scattering regarding isospin violation and higher orders in αs for off atomic nuclei, as well as indirect searches [18]. As the heavy-quark loops. long as the WIMP mass fulfills mχ < mh/2, limits on In addition, the standard decomposition (2) does not the invisible decay width of the Higgs boson are in a one- reflect that atomic nuclei are strongly interacting many- to-one correspondence with direct-detection limits on the nucleon systems, so that not only the one-body (single- WIMP–nucleon cross section σχN , provided the nuclear nucleon) matrix elements enter. As first pointed out physics input is sufficiently under control. in [31], corrections to this picture from two-body cur- Such comparison relies on the coupling of the Higgs rents, where the Higgs couples to a pion exchanged be- boson to a single nucleon, which proceeds via the nucleon tween two nucleons, can become important. Such cor- 2 χ q χ g Modifications of the leading result (2) arise from two Q sources: isospin-breaking and perturbative corrections. For u- and d-quarks the dominant isospin-breaking effects h h N N cancel in the sum fu + fd , with remaining effects of the χ q χ g size [45] str FIG. 1: Higgs-mediated interaction of the WIMP χ with the (mp − mn) fp − fn = , (6) light quarks q = u,d,s (left) and the gluon field g by closing mN a heavy-quark loop Q = c,b,t (right). str where (mp − mn) denotes the portion of the proton– neutron mass difference proportional to md − mu. Ina rections are conveniently addressed within chiral effec- nucleus, this would lead to a relative correction tive field theory (EFT) [32, 33], an expansion around IV ∆fN Z − N fp − fn < the chiral limit of QCD in terms of momenta and quark = ∼ 0.1%, (7) masses, as successfully demonstrated in previous applica- fN A 2fN tions to WIMP–nucleus scattering [34–41] (related work where Z/N refer to proton/neutron number, A = Z +N, using chiral EFT for WIMP–nucleon interactions [39] is and we have inserted values for a typical Xe isotope. restricted to the coupling to a single nucleon [42, 43] or Thus, isospin-violating effects in fN are small [46] and focuses on the calculation of the nuclear states based on will be neglected in the following. The perturbative cor- chiral EFT forces [44]). In this work, we extend previ- rections have been worked out in detail in [47]. Here, we ous calculations limited to xenon [40] to a broad range of use the full result for the c quark target nuclei and show that the corrections from the cou- N N 4 −1 pling to two nucleons are largely coherent and can thus fc =0.083 − 0.103 fq + O αs,mc , (8) be absorbed into a redefinition of the single-nucleon fN . q=Xu,d,s The combination of the one- and two-body contributions as well as the O(α ) correction for q = b,t. to fN then allows us to compare the limits on Higgs- s portal models from collider and direct-detection experi- In this way, we are left with the determination of the ments using state-of-the-art nuclear physics input. scalar couplings of the light quarks in the isospin limit. In general, these couplings cannot be measured in ex- periment due to the absence of a scalar source, leav- SCALAR COUPLINGS OF THE NUCLEON ing lattice QCD as the only viable option. However, an important exception concerns the scalar couplings of u- The evaluation of the quark-level operators from Fig. 1 and d-quarks, which the Cheng–Dashen low-energy the- between nucleon states requires knowledge of the scalar orem [48, 49] relates to the πN scattering amplitude, of- couplings defined in Eq. (1). To arrive at Eq. (2), the cou- fering a unique opportunity to test the results of lattice- plings fQ of the heavy quarks Q = c,b,t are eliminated QCD calculations with experiment. More precisely, the in favor of the light quarks’ according to the following low-energy theorem allows one to extract the πN σ-term N N procedure [19]: at O(αs), the QCD trace anomaly for Nf σπN = mN (fu + fd ) from the isoscalar amplitude eval- active degrees of freedom reads uated at a particular kinematic point in the subthreshold region. The required analytic continuation can be per- Nf formed in a stable manner using constraints from ana- µ 2Nf αs a µν θ = mqi q¯iqi − 11 − G G , (3) µ 3 8π µν a lyticity and unitarity in the form of Roy–Steiner equa- Xi=1 tions [28, 50–55], with the result that ultimately the a value of σπN is fully determined by the amplitude at with gluon field-strength tensor Gµν . Integrating out the heavy quarks then yields for each flavor threshold, i.e. the πN scattering lengths. These, in turn, are known to high accuracy from measurements in pio- N αs a µν nic atoms [56–60]. The result σπN = 59.1(3.5) MeV [28] mN f = hN|mQQQ¯ |Ni = − hN|G G |Ni Q 12π µν a agrees with calculations based on chiral perturbation the- 2 µ ory, once the low-energy constants are extracted from a = hN|θ |Ni− hN|mqqq¯ |Ni , (4) 27 µ phase-shift solution consistent with the pionic-atom scat- q=Xu,d,s tering lengths [27].