PandaX-II Constraints on Spin-Dependent WIMP-Nucleon Effective Interactions

Jingkai Xiaa, Abdusalam Abdukerimb, Wei Chena, Xun Chena, Yunhua Chenc, Xiangyi Cuia, Deqing Fangd, Changbo Fua, Karl Gibonia, Franco Giuliania, Linhui Gua, Xuyuan Guoc, Zhifan Guoe, Ke Hana, Changda Hea, Shengming Hec, Di Huanga, Xingtao Huangf, Zhou Huanga, Peng Jig, Xiangdong Jia,h,1, Yonglin Jue, Shaoli Lia, Heng Lina, Huaxuan Liue, Jianglai Liua,h, Yugang Mad, Yajun Maoi, Kaixiang Nia, Jinhua Ningc, Xiangxiang Rena, Fang Shia, Andi Tanj, Anqing Wangf, Cheng Wange, Hongwei Wangd, Meng Wangf, Qiuhong Wangd, Siguang Wangi, Xiuli Wange, Xuming Wanga, Zhou Wange, Mengmeng Wug, Shiyong Wuc, Mengjiao Xiaoj,k, Pengwei Xieh, Binbin Yanf, Jijun Yanga, Yong Yanga, Chunxu Yug, Jumin Yuanf, Jianfeng Yuec, Dan Zhangj, Hongguang Zhanga, Tao Zhanga, Li Zhaoa, Qibin Zhengl, Jifang Zhouc, Ning Zhoua,∗, Xiaopeng Zhoui, Wick C. Haxtonm,n,∗∗

aSchool of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai Key Laboratory for Particle Physics and Cosmology, Shanghai 200240, China bSchool of Physics and Technology, Xinjiang University, Ur¨umqi830046,¨ China cYalong River Hydropower Development Company, Ltd., 288 Shuanglin Road, Chengdu 610051, China dShanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China eSchool of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China fSchool of Physics and Key Laboratory of Particle Physics and Particle Irradiation (MOE), , Jinan 250100, China gSchool of Physics, Nankai University, Tianjin 300071, China hTsung-Dao Lee Institute, Shanghai 200240, China iSchool of Physics, , Beijing 100871,China jDepartment of Physics, University of Maryland, College Park, Maryland 20742, USA kCenter of High Energy Physics, Peking University, Beijing 100871, China lSchool of Medical Instrument and Food Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China mDepartment of Physics, University of California, Berkeley, CA 94720, USA nLawrence Berkeley National Laboratory, Berkeley, CA 94720, USA

Abstract We present PandaX-II constraints on candidate WIMP-nucleon effective interactions involving the nucleon or WIMP spin, including, in addition to standard axial spin-dependent (SD) scattering, various couplings among vector and axial currents, magnetic and electric dipole moments, and tensor interactions. The data set corresponding to a total exposure of 54-ton-days is reanalyzed to determine constraints as a function of the WIMP mass and isospin coupling. We obtain WIMP-nucleon cross section bounds of 1.6 × 10−41cm2 and 9.0 × 10−42cm2 (90% c.l.) for neutron-only SD 2 and tensor coupling, respectively, for a mass MWIMP ∼ 40 GeV/c . The SD limits are the best currently available for 2 MWIMP > 40 GeV/c . We show that PandaX-II has reached a sensitivity sufficient to probe a variety of other candidate spin-dependent interactions at the weak scale. Keywords: PandaX-II experiment, WIMP , Spin-dependent effective interactions

Astrophysical and cosmological observations indicate periments. In the analysis of direct detection experiments, that a large amount of non-luminous dark matter (DM) frequently it is assumed that the scattering arises from the exists in our universe, constituting ∼ 27% of the closure light-quark-level (u, s, d) effective interaction arXiv:1807.01936v2 [hep-ex] 1 Apr 2019 density. However, the exact of DM remains a mys- GF X  VV µ AA µ  tery. One intriguing DM candidate, a weakly-interacting L ∼ √ χγ¯ µχ c qγ¯ q +χγ ¯ µγ5χ c qγ¯ γ5q (1) 2 q q massive particle (WIMP), arises naturally in many exten- q sions of the standard model [1, 2]. Many WIMP searches which can be reduced to a nucleon-level operator useful have been performed, including direct detection of their in analyzing the nuclear response to WIMP scattering [3]. scattering off target nuclei, indirect detection of their de- Limits on the vector spin-independent (SI) and axial spin- cay or annihilation, and their production in collider ex- dependent (SD) WIMP-proton and WIMP-neutron cross SI,SD sections σp,n can then be derived. The recent strin- ∗Corresponding author gent direct detection null results obtained successively by ∗∗Corresponding author LUX [4, 5], PandaX [6, 7] and XENON [8, 9, 10] have Email addresses: [email protected] (Xiangdong Ji), significantly tighten the bounds on σSI,SD [11]. [email protected] (Ning Zhou), [email protected] (Wick C. p,n Haxton) The interaction of Eq. (1) was motivated by supersym- 1Spokesperson of PandaX Collaboration metric DM candidates, like the , that can natu-

Preprint submitted to Physics Letters B April 2, 2019 7 µ ¯ 5 mN rally account for the DM relic density. The motivation to Lint ≡ χγ¯ χNγµγ N → −2O7 + 2 O9 focus exclusively on such candidates has weakened due in mχ particular to collider constraints [12]. An alternative ap- proach, effective field theory (EFT) [13, 14, 15], has gained 13 µ 5 ¯ Lint ≡ χγ¯ γ χNγµN → 2O8 + 2O9 favor because it allows one to do an analysis [16, 17, 18, 19] free of theory assumptions. One selects an EFT scale – 15 µ 5 ¯ 5 e.g., the light-quark or a nucleon scale – and constructs Lint ≡ χγ¯ γ χNγµγ N → −4O4 . (2) a complete basis of effective operators to a given order, ~ ~ taking into account all general symmetries limiting that with mN the nucleon mass and O1 = 1χ1N , O4 = Sχ · SN , ⊥ ⊥ ~q O7 = S~N ·~v , O8 = S~χ ·~v , and O9 = i(S~χ ×S~N )· the basis. The underlying UV theory of DM will reduce at mN 5 that scale to some definite combination of the basis opera- nonrelativistic operators of [14, 15]. While Lint generates tors, regardless of its nature. Experimentalists can explore the standard SI interaction, the other interactions involve the sensitivities of their detectors to the basis operators, an axial coupling and thus depend on spin. to make sure they are probing all possibilities. The EFT One can equally well start with a basis of light-quark approach has shown 1) relative experimental sensitivities effective operators, reducing these via chiral EFT to their depend on the operator choice, and 2) direct detection is nonrelativistic nucleon equivalents [20, 21]. The spin- potentially more powerful than might appear from SI/SD dependent nucleon-level operators arising from the axial analyses, as six (not two) independent constraints on DM part of Eq. (1) (the standard SD interaction) and from can be obtained, in principle [15]. the light-quark tensor interaction will be considered here. The PandaX-II detector, located in the China Jinping The dimension-five operators coupling the WIMP mag- Underground Laboratory (CJPL), is a dual-phase netic or electric dipole moments with the nucleon’s vector time-projection chamber with 580 kg of liquid xenon in current, and the dimension-six operator coupling WIMP the sensitive target volume. When the incoming WIMP and nucleon magnetic moments, are examples of other po- scatters off a xenon nucleus, both the prompt scintillation tential sources of spin-dependent scattering, photons (S1) in the liquid and the delayed proportional 9 µν qν ¯ scintillation photons (S2) in the gas are collected by 55 top Lint ≡ χiσ¯ χNγµN mM and 55 bottom Hamamatsu R11410-20 3-inch photomulti- 2 2 plier tubes. The experiment has an accumulated exposure ~q 2mN 2mN ~q → − O1 + O5 − ( 2 O4 − O6) of 54-ton-days [6]. The previously reported analysis for 2mχmM mM mM mN the standard (isoscalar) SI interaction yielded a 90% ex- clusion limit on σSI of 8.6 × 10−47cm2 for a WIMP mass 17 µν qν 5 ¯ 2mN of 40 GeV/c2. (These bounds were recently superseded Lint ≡ iχiσ¯ γ χNγµN → O11 mM mM by XENON1T results [10].) In this paper we present

PandaX-II constraints a variety of candidate interactions α 2 2 that depend on nucleon or WIMP spin. 10 µν qν ¯ q ~q mN Lint ≡ χiσ¯ χNiσµα N → 4( 2 O4 − 2 O6) Spin-dependent interactions other than SD can arise mM mM mM mM from WIMP magnetic and electric dipole moments, vector- (3) axial interference terms, tensor interactions, etc. While ~q ⊥ ~q ~q Here O5 = iS~χ ·( ×~v ), O6 = (S~χ · )(S~N · ), and the associated nucleon-level effective operators are conven- mN mN mN ~q tionally expressed in covariant form, they can be rewrit- O11 = iS~χ · . The dependence on ~q implies nuclear form mN ten in terms of the Galilean-invariant EFT basis used here, factors that peak at larger momentum transfers, influenc- convenient for nonrelativistic shell model (SM) treatments ing experimental strategies for optimally constraining such of the nuclear physics. This basis consists of fourteen interactions. We set mM ≡ mN , normalizing both WIMP operators generated at next-to-next-to-leading-order from and nucleon moments to the nucleon scale. i the nucleon-WIMP perpendicular relative velocity oper- Each operator Lint can have independent couplings to ator ~v⊥, the momentum transfer ~q, and the WIMP and protons and neutrons, or equivalently to isospin [15], nucleon spins, S~ and S~ [15]. χ N 0 1 p n We specialize to the scattering of a spin-1/2 WIMPs off di di di 1 + τ3 di 1 − τ3 2 + 2 τ3 = 2 + 2 (4) a natural xenon target, exploring four dimension-four and mV mV mV 2 mV 2 three higher dimension effective interactions, selected from where τ is the Pauli isospin operator. The couplings d Table 1 of Ref. [15]. The operator dimension is defined as 3 i are dimensionless, defined relative to the weak scale 4 + number of powers of mM in the denominator, where m is a scale that governs the strength of the WIMP − 1 M mV ≡ hvi = (2GF ) 2 = 246.2 GeV and nucleon moments being coupled. The dimension-four operators are the V-A interactions with hvi the Higgs vacuum expectation value. We will 1 0 consider isoscalar (di = 0) and isovector (di = 0) inter- 5 µ ¯ p p Lint ≡ χγ¯ χNγµN → O1 n n actions, for which di = di and −di = di , respectively; 2 1 0 as well as couplings only to protons (di = di ) or neu- evaluating Eq. (6) for a local WIMP mass density of ρχ = 1 0 2 3 trons (di = −di ). For our higher-dimension operators, 0.3 GeV/c /cm , assuming a Maxwellian WIMP velocity the di’s also encode information about the absolute size distribution peaked at v0 = 220 km/s and truncated at of the WIMP electric and dipole moments, given that we the galactic escape velocity vesc = 544 km/s. have normalized to nucleonic values (mM ≡ mN ). The nuclear response functions for DM elastic scatter- As the spin response of Xe is largely governed by the ing must be calculated before experimental limits can be 129,131 unpaired neutrons in Xe, one expects the 1−τ3 pro- converted to bounds on the operator coefficients of Eq. jection of nucleon-spin operators to dominate. (4). For each contributing Xe isotope and needed opera- tor, we performed full-basis shell-model calculations using the GCN5082 [22] interaction (so named because the SM 1011 5 1011 Lint isoscalar isovector 7 8 8 valence space resides between the shell closures at nucleon 10 Lint isoscalar isovector 10 13 Lint isoscalar isovector numbers 50 and 82). The calculations were done with- 4 15 4 10 Lint isoscalar isovector 10 (1/kg/day/keV)

R out truncation, using the SM code BIGSTICK to treat 1 1 bases that ranged to 9 billion Slater determinants [23]. dR/dE

10-4 10-4 Full GCN5082 Xe isotopes response functions are used in the Mathematica script of [24], an update of [15]. This -8 -8 10 10 script and the associated library of one-body nuclear den- mΧ=40GeV mΧ=400GeV 10-12 10-12 sity matrices are available on request from the authors 1011 10 20 30 40 50 60 70 80 901011 10 20 30 40 50 60 70 80 90 9 of [24]. The scripts of [24] and [15] were carefully cross- 8 Lint isoscalar isovector 8 10 10 10 checked against one another, to verify their consistency. Lint isoscalar isovector 4 17 4 10 Lint isoscalar isovector 10 Figure 1 shows computed recoil energy spectra. The up- per panels include the coherent isoscalar (N+Z) and isovec- 1 1 5 tor (N-Z) SI responses (operator Lint), which we show for 10-4 10-4 normalization. Though we used SM results, simple phe- 10-8 10-8 nomenological forms [25] also work well, below the diffrac-

mΧ=40GeV mΧ=400GeV 10-12 10-12 tion minimum, as the form factor is governed by its known 10 20 30 40 50 60 70 80 90 10 20 30 40 50 60 70 80 90 q = 0 value and the nuclear radius.

ER (keV) The remaining curves corresponding to interactions with 7 13 15 nucleon spin (Lint), WIMP spin (Lint), or both (Lint), Figure 1: Recoil energy spectra for the scattering of spin- 1 WIMPs 2 at two WIMP masses. All curves correspond to weak- on xenon nuclei, for WIMP masses of 40 GeV/c2 and 400 GeV/c2, p n p n 0 1 0 1 scale interactions: with di = di ≡ 1 (di = −di ≡ 1) for and unit isoscalar (di = 1, di = 0) and isovector (di = 0, di = 1) coupling. Top panel: dimension-four V-A interactions. Bottom isoscalar (isovector) coupling. The lower panels show the panel: higher dimension magnetic and electric dipole interactions. corresponding results for the magnetic and electric dipole 9 10 17 moment interactions, Lint, Lint, and Lint. The differential cross section for elastic scattering is The WIMP couplings di are constrained by PandaX-II rate limits. We use data from two low-background physics dσ(v, E ) dσ(v, ~q 2) R runs with a total exposure of 54-ton-days, Run 9 with 79.6 = 2mT 2 dER d~q live days in 2016 and Run 10 with 77.1 live days in 2017.   Calibration data from an AmBe source outside the cryo- 2mT 1 1 X 2 = 2  |M|  (5) genic vessel and tritium decays from CH3T injected into 4πv 2Jχ + 1 2J + 1 spins the Xe provided tests of the detector response to the nu- clear (NR) and electron (ER) recoil events, respectively [6]. where the square of the Galilean invariant amplitude M A NEST-based Monte Carlo (MC) simulation (Ref. [26] for is a product of WIMP and nuclear matrix elements and ER, Ref. [27] for NR) of the data in the (S1, S2) distri- is a function of initial WIMP velocity v (dimensionless, bution is optimized by adjusting the initial excitation-to- in units of c) and the three-momentum transfer ~q [15]. ionization ratio and the recombination fluctuation. A reli- Here J is the WIMP spin and J the nuclear ground state χ able tuned NEST model for PandaX-II is obtained for S1 angular momentum. In the long wavelength limit the nu- up to 50 photoelectrons (PEs). The S1 and S2 signal dis- clear response functions corresponding to matrix elements tributions for a given EFT interaction are simulated using of simple operators familiar from weak interactions, such the tuned NEST model and PandaX-II detector response as 1(i), ~σ(i), ~`(i), and ~σ(i) · ~`(i). The differential event parameters [6]. rate with respect to nuclear recoil energy is Our event selection criteria follow Ref. [6]: S1 from 3 Z to 45 PEs, S2 from 100 (raw, not corrected for electron dR ρχ dσ(v, ER) 3 = vf(~v)d v, (6) lifetime) to 10000 PEs, events lying within the 99.99% NR dER mχ dER acceptance, and the total fiducial target of 329±16 kg. The where the f(~v) is the normalized velocity distribution of backgrounds in Run 9 and Run 10, estimated in Ref. [6], the WIMP particles. We calculate WIMP signal rates by are dominated by 127Xe, tritium, other flat ER background 3 (85Kr, radon and detector gamma background), acciden- section, and compared with that from a large number of tal, and neutron contributions. toy MC calculations to derive the upper limits of the signal This search window can include nuclear recoil event en- yields at 90% confidence level (C.L.) [29, 30]. ergies up to 100 keVnr (NR energy) due to the smearing of S1 and S2, although the efficiency drops below 50% above 35 keV on average. For dimension-four operators, the 10-3 nr σ -40 2 n=10 cm , neutron coupled only majority of signal events have nuclear recoil energy below 10-4 35 keVnr. The overall signal selection efficiency is between 10-5 40% and 50%, very similar to that for the SI analysis in 10-6 (1/kg/day/keV)

Ref. [6]. For dimension-five or dimension-six operators, R -7 the analysis is more complicated due to the sharper form 10 -8 factor momentum dependence. For instance, the signal ef- dR/dE 10 9 -9 ficiency for the isovector Lint operator decreases from 65% 10 SD, m =40GeV to 6.5% as the WIMP mass increases from 40 to 400 GeV. -10 χ 10 SD, mχ=400GeV Thus efficiencies can be improved by adjusting search win- Tensor, m =40GeV 10-11 χ dows according to operator type and WIMP mass, as LUX Tensor, mχ=400GeV 10-12 has described [28]. Such strategies will be explored in fu- 0 20 40 60 80 100 ture PandaX analyses: detector calibration studies for S1 ER (keV) above 50 PE will be needed to implement such window adjustments. 10-5 σ -40 2 p=10 cm , proton coupled only 10-6 5 10-7 10 5 L isoscalar isovector int -8 L7 isoscalar isovector 10 int (1/kg/day/keV)

13 R 3 Lint isoscalar isovector -9 10 15 10 Lint isoscalar isovector -10 dR/dE 10 (dimensionless) i -11 d 10 10 SD, m =40GeV -12 χ 10 SD, m =400GeV -1 χ 10 Tensor, m =40GeV 10-13 χ Tensor, mχ=400GeV 10-14 10-3 0 20 40 60 80 100

ER (keV) 10-5 2 3 4 5 1 10 10 10 10 10 Figure 3: Recoil energy spectra for the scattering of spin- 2 WIMPs WIMP mass (GeV/c2) on xenon nuclei, for WIMP masses of 40 GeV/c2 and 400 GeV/c2, for SD and tensor interactions. The top (bottom) panel corresponds 105 to neutron-only (proton-only) coupling. 9 Lint isoscalar isovector L10 isoscalar isovector 3 int 17 The upper panel of Fig. 2 gives PandaX-II bounds on 10 L isoscalar isovector int 5 the coherent SI interaction Lint as well as the spin-spin 15 (dimensionless) i interaction Lint, related to the traditional SD interaction d 10 discussed below. They are well below the nominal weak

-1 scale d ∼ 1 over all or almost all WIMP mass range illus- 10 0 −5 trated. The isoscalar coupling limit |d5| < 5.6 · 10 at

-3 maximum sensitivity (∼ 40 GeV WIMP) corresponds to 10 SI −47 2 the PandaX-II cross section bound σp,n < 8.6 · 10 cm . But limits on the axial (WIMP) - vector (nucleus) cou- -5 10 pling d0,1 are also at the weak scale, illustrating the power 10 102 103 104 105 13 WIMP mass (GeV/c2) of current-generation experiments to probe momentum- or velocity-dependent interactions. The lower panel gives Figure 2: Exclusion limits on the coefficients of dimension-four (top the corresponding bounds for interaction involving WIMP panel) and higher dimension (bottom panel) operators. magnetic and electric dipole moments. The bounds are near and in one case (the WIMP electric dipole moment 17 We constrain the WIMP-nucleon EFT couplings di as a coupling to the nucleon vector current Lint) below the function of the WIMP mass, following the likelihood anal- nominal weak scale. Bounds given in the lower panel also ysis of Ref. [6]. A standard profile likelihood test statistic include an implicit dimensionless factor representing any was determined as a function of WIMP mass and cross needed rescaling of the WIMP moments to their physi- 4 cal values, relative to the nucleon scale we adopted via where µA (µp,n) is the reduced mass for scattering off the A 2 mM ∼ mN . nucleus (nucleon), and Sp,n(q ) is the nuclear spin struc- ture function that we take from [3], including exchange current corrections. Figure 3 shows recoil spectra for the -34 ) 2 10 SD and tensor cross sections. 10-35 The 90% C.L. cross section limits are shown in Fig. 4. PandaX-II (Tensor) 2 -36 The optimal bounds, obtained at MWIMP ∼ 40 GeV/c , 10 PandaX-II SD SD LUX limit σn (neutron-only) and σp (proton-only) to 1.6 × 10-37 −41 2 −40 2 XENON100 10 cm and 4.4 × 10 cm , respectively. The results 10-38 modestly improve existing LUX bounds [4] for MWIMP > 40 GeV/c2. 10-39 Models of asymmetric dark matter favor Dirac fermion 10-40 WIMPs, where candidate dimension-4 effective operators µν 10-41 include the SD and tensorχσ ¯ µν χ qσ¯ q operators. Con- SD WIMP-neutronSD cross section (cm straints from direct detection are particularly competitive 10-42 3 5 for the latter [38]. The tensor interaction generates, after 10 102 10 104 10 WIMP mass (GeV/c2) a leading-order chiral reduction [21], the nucleon-level op- erator 8O4. The dashed curve in the top panel of Fig. 4 -32 ) T 2 10 PandaX-II IceCube(ν+ν-) shows the neutron-only PandaX-II limits. Bounds on σn T −42 2 −38 2 LUX SuperK(τ+τ-) and σ of 9.0×10 cm and 2.2×10 cm , respectively, 10-34 p XENON100 CMS are found at MWIMP ∼ 40 GeV. PICO60 ATLAS In conclusion, we have presented new limits on a candi- 10-36 date spin-dependent WIMP interactions, using PandaX- 10-38 II Run 9 and Run 10 data with an exposure of 54-ton- days. In addition to the standard SD interaction, we con- 10-40 sidered vector-axial vector interferences, interactions gen- erated by WIMP magnetic and electric dipole moments, 10-42 and tensor interactions. We showed that PandaX-II has SD WIMP-protonSD cross section (cm achieved sufficient sensitivity to probe certain velocity- -44 10 and momentum-dependent interactions at the weak scale. 10 102 103 104 105 WIMP mass (GeV/c2) We obtained the most stringent upper limits to date on SD 2 σn for MWIMP above 40 GeV/c , with a lowest excluded Figure 4: PandaX-II exclusion limits on the WIMP-nucleon cross value of 1.6 × 10−41cm2 at 40 GeV/c2, 90% c.l. The section for the standard SD interaction assuming neutron-only (top corresponding proton and tensor interaction constraints panel) and proton-only (bottom panel) coupling. The 1σ sensi- are σSD < 4.4 × 10−40cm2, σT < 9.0 × 10−42cm2, and tivity bands are shown in green. Also shown are recent results p n T −38 2 from LUX [5], XENON100 [9], ATLAS [31], CMS [32], PICO- σp < 2.2 × 10 cm . 2L [33],PICO-60 [34, 35], IceCube [36] and Super-K [37]. The dashed line (top panel) gives PandaX-II limits on tensor WIMP-neutron cou- plings. Acknowledgments

An often-used measure of experimental sensitivity to This project has been supported by a 985-III grant spin-dependent WIMP scattering is provided by the SD from Shanghai Jiao Tong University, grants from National axial nucleon-level operator obtained from Eq. (1), Science Foundation of China (Nos. 11435008, 11455001, 11505112, 11525522, 11775141 and 11755001), a grant 2 2 gA(q ) mN gP (q ) from the Ministry of Science and Technology of China OSD ≡ O4 − O6 (7) gA(0) 2gA(0) (No. 2016YFA0400301). We thank the Office of Sci- ence and Technology, Shanghai Municipal Government with the specific combination of nucleon-level axial and (No. 11DZ2260700, No. 16DZ2260200, No. 18JC1410200) induced pseudoscalar operators dependent on isospin and the Key Laboratory for Particle Physics, Astrophysics AA through the Cq . The pseudoscalar coupling is enhanced and Cosmology, Ministry of Education, for important in the isovector channel by pion-pole dominance. Com- support. This work is supported in part by the Chi- parisons are typically made by selecting proton-only (p) nese Academy of Sciences Center for Excellence in Par- or neutron-only (n) couplings. The nuclear cross section ticle Physics (CCEPP) and Hongwen Foundation in Hong SD A SD σp,n is related to the nucleon cross section σp,n by Kong. WH is supported by the US Department of En- ergy (de-sc00046548 and DE-AC02-05CH11231) and the  2 SD A 2 4π µA A 2 SD National Science Foundation (PHY-1630782). We would σp,n (q ) = Sp,n(q )σp,n (8) 3(2J + 1) µp,n like to thank Dr. Yue-Lin Sming Tsai and Dr. Zuowei 5 Liu for the useful discussions on the EFT models in di- [19] S. Kang, S. Scopel, G. Tomar, J.-H. Yoon, Present and pro- rect detection experiment. Finally, we thank the CJPL jected sensitivities of dark matter direct detection experiments administration and the Yalong River Hydropower Devel- to effective wimp-nucleus couplings, arXiv:1805.06113. [20] M. Hoeferichter, P. Klos, A. Schwenk, Chiral power counting of opment Company Ltd. for indispensable logistical support one- and two-body currents in direct detection of dark matter, and other help. Phys. Lett. B 746 (2015) 410. [21] F. Bishara, J. Brod, B. Grinstein, J. Zupan, From quarks to nucleons in dark matter direct detection, JHEP 1711 (2017) 059. References [22] J. Menendez, A. Poves, E. Caurier, F. Nowacki, Disassembling the nuclear matrix elements of the neutrinoless ββ decay, Nucl. References Phys. A 818 (2009) 139, see also L. Vietze et al., Phys. Rev. D 91, 043520 (2015) for a comparison of GCN5082 Xe isotope spectroscopic predictions to experiment and for full-space cal- [1] J. L. Feng, Dark matter candidates from particle physics and culations of certain isotopes (132Xe, 134Xe, 136Xe) that provide methods of detection, Ann. Rev. Astro. Astrophys. 48 (2010) valuable numerical cross checks on present work; in other cases 495–545. doi:10.1146/annurev-astro-0082708-101659. Vietze et al. had to make truncations, which are avoided in the [2] G. Bertone, D. Hooper, J. Silk, Particle dark matter: evidence, present work. doi:10.1016/j.nuclphysa.2008.12.005. candidates and constraints, Phys. Rep. 405 (2005) 279. doi: [23] C. W. Johnson, Current status of very-large-basis hamilto- 10.1016/j.physrep.2004.08.031. nian diagonalizations for nuclear physics, arXiv:1809.07869To [3] P. Klos, J. Menendez, D. Gazit, A. Schwenk, Large-scale nuclear be published in the Proc. of CIPANP18; C. W. Johnson et al., structure calculations for spin-dependent wimp scattering with arXiv:1801.08432. chiral effective field theory currents, Phys. Rev. D 88 (2013) [24] A. L. Fitzpatrick, W. C. Haxton, C. Johnson, K. S. McElvain, 083516, see also J. Engel, S. Pittel, and P. Vogel, Int. J. Mod. in preparation for Ann. Rev. Nucl. Part. Sci. Phys. E 1, 1 (1992). doi:10.1103/PhysRevD.88.083516. [25] J. D. Lewin, P. F. Smith, Review of mathematics, numerical [4] D. Akerib, et al., Results from a search for dark matter in the factors, and corrections for dark matter experiments based on complete lux exposure, Phys. Rev. Lett. 118 (2017) 021303. elastic nuclear recoil, Astropart. Phys. 6 (1996) 87. doi:10. doi:10.1103/PhysRevLett.118.021303. 1016/S0927-6505(96)00047-3. [5] D. S. Akerib, et al., Limits on spin-dependent wimp-nucleon [26] M. Szydagis, A. Fyhrie, D. Thorngren, M. Tripathi, Enhance- cross section obtained from the complete lux exposure, Phys. ment of nest capabilities for simulating low-energy recoils in liq- Rev. Lett. 118 (2017) 251302. doi:10.1103/PhysRevLett.118. uid xenon, JINST 8 (2013) C10003. doi:10.1088/1748-0221/ 251302. 8/10/C10003. [6] X. Cui, et al., Dark matter results from 54-ton-day exposure [27] B. Lenardo, K. Kazkaz, A. Manalaysay, J. Mock, M. Szydagis, of pandax-ii experiment, Phys. Rev. Lett. 119 (2017) 181302. M. Tripathi, A global analysis of light and charge yields in liquid doi:10.1103/PhysRevLett.119.181302. xenon, IEEE Trans. Nucl. Sci. 62 (2015) 3387. doi:10.1109/ [7] C. Fu, et al., Spin-dependent weakly-interacting-massive- TNS.2015.2481322. particle–nucleon cross section limits from first data of pandax-ii [28] N. Larsen, the LUX Collaboration, Lux sensitivity to effective experiment, Phys. Rev. Lett. 118 (2017) 071301. doi:10.1103/ field theory interactions, in: Proceedings of CIPANP18 (to be PhysRevLett.118.071301. published by SLAC eConf), CIPANP18, 2018. [8] E. Aprile, et al., First dark matter search results from the [29] A. L. Read, Presentation of search results: the cls technique, J. xenon1t experiment, Phys. Rev. Lett. 119 (2017) 181301. doi: Phys. G 28 (2002) 2693. doi:10.1088/0954-3899/28/10/313. 10.1103/PhysRevLett.119.181301. [30] T. Junk, Confidence level computation for combining searches [9] E. Aprile, et al., Xenon100 dark matter results from a com- with small statistics, Nucl. Instrum. Meth. A 434 (1999) 435. bination of 477 live days, Phys. Rev. D 94 (2016) 122001. doi:10.1016/S0168-9002(99)00498-2. doi:10.1103/PhysRevD.94.122001. [31] M. Aabound, et al., Search for dark matter and other new phe- [10] E. Aprile, et al., Dark matter search results from a one nomena in events with an energetic jet and large missing trans- tonne×year exposure of xenon1t, Phys. Rev. Lett. 121 (2018) verse momentum using the atlas detector, JHEP 01 (2018) 126. 111302. doi:10.1103/PhysRevLett.121.111302. doi:10.1007/JHEP01(2018)126. [11] J. Liu, X. Chen, X. Ji, Current status of direct dark matter [32] A. M. Sirunyan, et al., Search for dark matter produced with√ an detection experiments, Nature Phys 13 (2017) 212. doi:10. energetic jet or a hadronically decaying w or z boson at s = 1038/NPHYS4039. 13 tev, JHEP 07 (2017) 014. doi:10.1007/JHEP07(2017)014. [12] M. E. Peskin, Supersymmetric dark matter in the harsh light [33] C. Amole, et al., Improved dark matter search results from of the large hadron collider, PNAS 112 (2015) 12256. -2l run 2, Phys. Rev. D 93 (2016) 061101. doi:10.1103/ [13] J. Fan, M. Reece, L.-T. Wang, Non-relativistic effective theory PhysRevD.93.061101. of dark matter direct detection, JCAP 1011 (2010) 042. doi: [34] C. Amole, et al., Dark matter search results from the pico-60 . 10.1088/1475-7516/2010/11/042 cf3i bubble chamber, Phys. Rev. D 93 (2016) 052014. doi: [14] A. L. Fitzpatrick, W. C. Haxton, E. Katz, N. Lubbers, Y. Xu, 10.1103/PhysRevD.93.052014. Th effective theory of dark matter detection, JCAP 02 (2013) [35] C. Amole, et al., Dark matter search results from the pico- 004. 60 c3f8 bubble chamber, Phys. Rev. Lett. 25 (2017) 251301. [15] N. Anand, A. L. Fitzpatrick, W. C. Haxton, Weakly interacting doi:10.1103/PhysRevLett.118.251301. massive particle-nucleus elastic scattering response, Phys. Rev. [36] M. G. Aartsen, et al., Improved limits on dark matter an- C 89 (2014) 065501. doi:10.1103/PhysRevC.89.065501. nihilation in the sun with the 79-string icecube detector and [16] K. Schneck, et al., Dark matter effective field theory scattering implications for supersymmetry, JCAP 04 (2016) 022. doi: in direct detection experiments, Phys. Rev. D 91 (2015) 092004. 10.1088/1475-7516/2016/04/022. doi:10.1103/PhysRevD.91.092004. [37] K. Choi, et al., Search for neutrinos from annihilation of cap- [17] E. Aprile, et al., Effective field theory search for high-energy tured low-mass dark matter particles in the sun by super- nuclear recoils using the xenon100 dark matter detector, Phys. kamiokande, Phys. Rev. Lett. 114 (2015) 141301. doi:10.1103/ Rev. D 96 (2017) 042004. doi:10.1103/PhysRevD.96.042004. PhysRevLett.114.141301. [18] Z. Liu, Y. Su, Y. S. Tsai, B. Yu, Q. Yuan, A combined analysis [38] M. A. Buckley, Asymmetric dark matter and effective operators, of pandax, lux and xenon1t experiments within the framework Phys. Rev. D 84 (2011) 043510. doi:10.1103/PhysRevD.84. of dark matter effective theory, JHEP 11 (2017) 024. doi:10. 043510. 1007/JHEP11(2017)024. 6