Neutron EDM constrains direct dark matter detection prospects

A non-relativistic effective field theory (NREFT) offers a bottom-up framework to classify Dark Matter (DM)-nucleon interactions relevant for scattering at direct detection experiments by organizing the interactions in powers of the momentum transfer $\vec{q}$ and DM velocity $\vec{v}$. This approach generates a number of operators including P-odd and T-odd operators; these can only be generated from a relativistic theory with CP violating interactions. We consider the leading order P-odd, T- odd operators viz. $\mathcal{O}_{10}$, $\mathcal{O}_{11}$ and $\mathcal{O}_{12}$ and compare the constraints on these operators from leading direct detection searches and from the bound on the neutron EDM (nEDM). We perform our analysis using simplified models with charged mediators and compute the loop diagrams contributing to the nEDM. We find that constraints on the DM scattering cross section from the bound on the nEDM are several orders of magnitude stronger than the limits from direct searches, and even well below the neutrino floor for such NREFT operators, for the entire sub-GeV to TeV DM mass range. This indicates that these operators need not be considered when analyzing data from present or future direct dark matter detection experiments.


I. INTRODUCTION
The particle identity of Dark Matter (DM) has not yet been understood despite overwhelming evidence for its gravitational interactions. Dedicated direct searches exploit the possible non-gravitational interactions of DM by looking for the recoil of target nuclei due to elastic scattering and have already probed a large fraction of DM mass and cross-section parameter space [1][2][3][4][5][6][7][8][9][10].
Much of the direct detection search strategy for roughly three decades has been focused on a category of DM candidates known as Weakly Interacting Massive Particles (WIMPs). WIMPs with weak-scale masses, produced in the early universe via freeze-out from the thermal plasma [11], can easily match the measured DM relic abundance via roughly electroweak-strength couplings to Standard Model (SM) particles. Most analyses assume DM-nucleon interactions to be dominated by only two operators, which describe Spin-Independent (SI) or Spin-Dependent (SD) interactions in the limit of vanishing WIMP velocity. The leading bounds on SI DM-nucleon cross section and SD DM-neutron cross section are currently set by XENON1T [1,3]; for a 30 GeV WIMP, the bounds are 4.1 × 10 −47 cm 2 and 6.3 × 10 −42 cm 2 , respectively. The best direct detection bound on the SD DM-proton cross section comes from the PICO-60 experiment [7]; for 25 GeV DM mass cross sections above 2.5 × 10 −41 cm 2 are excluded.
The standard SI and SD operators are actually the leading (zeroth order) terms of an EFT organized in powers of small expansion parameters such as DM velocity v and three-momentum transfer normalized to the nucleon mass | q | /m N . Since the DM velocity is relatively small in the solar neighborhood, v/c ∼ O(10 −3 ), the momentum exchange is limited to | q | < ∼ min[vm χ , vm A ] < ∼ O(100 MeV); here m χ is the WIMP mass and m A is the mass of the target nucleus. This is sizable on nuclear physics scales but certainly well below the electroweak scale as well as the range of WIMP masses to which these experiments are sensitive. Keeping only the leading order terms therefore a priori seems justified.
Nevertheless a non-relativistic effective field theory (NREFT) was developed which retains NLO and NNLO terms by modeling the nucleus after taking into account the finite spatial extent of its charge and spin densities [12][13][14][15]. Working up to second order then generates a total of 14 operators, some of which lead to novel nuclear responses motivating the inclusion of angular-momentum dependent (LD) as well as spin and angular-momentum dependent (LSD) interactions, along with SI and SD. Moreover, these operators can interfere with each other and can lead to appreciable differences from the usual exponentially falling recoil energy spectra. It has recently been pointed out that this affects the optimization of the experimental search strategy, in particular by allowing events with larger recoil energy [16]. Considerable effort has also been devoted to analyzing published data in this framework [17][18][19][20].
The NREFT contains a total of 28 free parameters if operators involving neutrons and protons are treated separately. This large parameter space can be difficult to probe in its entirety. There have been a number of global analyses of the full multi-dimensional parameter space defining the NREFT using constraints from ongoing and planned direct detection experiments [21][22][23][24][25][26][27][28]. One conclusion is that direct searches constrain momentum or velocity suppressed operators that are odd under P (parity) and T (time reversal) transformations about as strongly as zeroth order SD interactions, but the former have been less explored by both theorists and experimenters. This motivates analyses of particle physics aspects of models that can give rise to these operators in the non-relativistic limit.
From a particle physics perspective, the CPT theorem implies that all P-and T-odd NREFT operators must necessarily arise from a CP violating theory, where C refers to charge conjugation. However, CP violation (CPV) is highly constrained experimentally [29] and therefore bounds on CPV observables can be used to constrain such NREFT operators. The neutron EDM (nEDM) is one such observable which probes flavor diagonal CP violation and provides one of the most stringent bounds on CPV in extensions of the SM. Ongoing nEDM experiments are sensitive to a possible signal which is many orders of magnitude larger than the SM prediction [29,30]. The current constraint is |d n | < 2.9 × 10 −26 e · cm (90 % C.L.) . (1.1) Various works have linked well-motivated extensions of the SM to low energy NREFT operators [31][32][33][34][35]. In particular, Ref. [36] lists a set of simplified models for scalar, spinor and vector DM candidates and derives the full set of NREFT operators in terms of the parameters for each simplified model. In this letter, we investigate the P-and T-odd operators in the NREFT with the help of simplified models that respect SU (3) × U (1) em invariance and are renormalizable, as in [36]. The models in this work extend the SM by a WIMP candidate and a mediator particle. We show that in these models the Wilson coefficients of the CPV NREFT operators and of the neutron EDM scale with the same combination of Yukawa couplings. As a result, we find that the constraints from the nEDM on the scattering cross section are many orders of magnitude stronger for both CPV NREFT operators than the current XENON1T bound and even lie below the irreducible background from coherent neutrino scattering ("neutrino floor") [37][38][39][40]. This casts doubt on the importance of these operators for the analysis of future experiments [41] and for the formalism of NREFT itself.
The rest of this letter is organized as follows. In Section II we introduce the NREFT formalism, the matching procedure of a relativistic theory to NREFT operators and review the computation of the differential event rate for elastic WIMP-nucleus scattering. In Section III we consider two simplified models yielding these P-and T-odd CPV NREFT operators and compute the nEDM at 1-loop order for each case. Section IV contains our main numerical results, comparing the nEDM and the leading direct detection bound. We conclude in Section V.

II. NREFT FORMALISM
A non-relativistic effective field theory (NREFT) of elastic scattering between DM and nuclei which goes beyond the strict velocity v → 0 limit [12][13][14][15] offers a parametrically richer description than the traditional relativistic four-field effective operator approach. Since an incoming DM particle striking the target nuclei is fairly slow in the target rest frame, v/c ∼ O(10 −3 ), an NREFT is well suited to describe the scattering process. We briefly review this formalism for completeness.
We begin by discussing the basis of NREFT operators that describe the interactions of a WIMP χ with a nucleon N at low relative velocities. Galilean symmetry restricts the operator basis to consist of quantities which are invariant under Galilean boosts. The three-momentum transfer q and relative initial velocity v = v χ,in − v N,in satisfy this criterion, as does any quantity that is computed only from these three-vectors. Since any two relative velocities are sufficient to describe the scattering, it is helpful to consider orthogonal combinations such that the scalar product of the two invariants vanishes. The transverse velocity v ⊥ ≡ v + q 2µ N is by construction orthogonal to the momentum transfer q, where µ N = m χ m N /(m χ + m N ) is the reduced mass of the WIMP-nucleon system; v ⊥ · q = 0 follows from energy and momentum conservation.
Along with the two three-vectors q and v ⊥ , the scattering can be fully described by the spins of the target nucleon S N and of the DM particle S χ (when non-zero), which are also invariant under a Galilean transformation. The operator basis constructed so far is not so convenient for the construction of an interaction Hamiltonian since it is not entirely Hermitean; for instance the momentum transfer q is anti-Hermitean. Therefore, the most general non-relativistic interactions for elastic DM-nucleus scattering can be written down as functions of the following four Hermitean, Galilean invariant quantities: (2.1) A set of linearly independent operators O i can be constructed using the building blocks in Eq. (2.1). Table I lists the set of operators which are obtained when terms up to second order in momentum exchange q are retained. The operators O 1 and O 4 correspond to the traditional SI and SD interaction, respectively.  [14] by defining the operators normalized by the nucleon mass m N in order to have a dimensionless basis. We omit the invariant O 2 = v 2 ⊥ because it is a second order correction to the SI operator O 1 , as well as In deriving the NREFT operators, no requirement to obey discrete symmetries such as invariance under spatial parity P and time reversal T transformations has been imposed. 1 Therefore, the operators in Table I can be classified according to their P and T quantum numbers. Since velocities and three momenta are odd under both P and T whereas spin is odd under T but even under P, the following operators are both P-odd and T-odd: P-odd and T-odd ; (2.2) note that the anti-linear transformation T maps a complex number to its complex conjugate.
In order to match some UV-complete or simplified model for DM onto the NREFT, a two-step procedure is used. In the first step one integrates out the heavy mediator field(s), leading to a still relativistic but non-renormalizable description in terms of four-field DM-quark operators. If WIMPs are spin-1/2 fermions the nucleon and DM bilinears can be formed using a basis of 16 linear hermitean matrices Γ χ,N ∈ {I 4 , iγ 5 , γ µ , γ µ γ 5 , σ µν }. For scalar DM the Lorentz structure appearing in the scalar bilinear can be 2 {I, i ← → ∂ µ /m S } and the nucleon bilinear can be {I 4 , iγ 5 , γ µ , γ µ γ 5 }. In either case only combinations are allowed where all Lorentz indices are contracted between the two bilinears. By convention the scalar DM bilinears are multiplied with a factor of the scalar mass m S i.e. m S S † Γ S S such that the Wilson coefficients obtained have the same mass dimension as those of the fermionic WIMPs.
In the non-relativistic limit these dimension-6 four field operators composed of a nucleon bilinear and a DM bilinear reduce to O N and Oχ respectively; their products can be matched on to NREFT operators O i,NR of Table I with Wilson coefficients c N i (in the isospin basis) of mass dimension GeV −2 , with O i,NR ≡ Oχ · O N . Any relativistic theory with a DM candidate can thus be matched on to a particular linear combination of the NREFT operators. The rate, differential in the recoil energy E R , of elastic WIMP-nucleus scattering events per unit time and unit detector mass, dR/dE R , is then given by Here N T is the number of target nuclei in the detector, ρχ is the local DM density, mχ is the DM mass, m A is the mass of the target nuclei 3 , f (v) is the DM halo velocity distribution in the laboratory frame, v min is the minimum velocity required to cause a nuclear recoil E R , v max is related to the galactic DM escape velocity, and jχ and j A are the total spin of the DM and nucleus, respectively. Finally, the matrix element M sc is given by [14]: The result (in the isospin basis) has been factorized into WIMP response functions R k and nuclear response functions W k . Both response functions depend on the momentum transfer. The former also depend on the relative WIMPnucleon velocity and encode the particle physics of the scattering process, and thus explicitly depend on the product of Wilson coefficients c N i c N j and can be computed perturbatively; see Ref. [14] for a list of explicit expressions. The latter encode the target-dependent nuclear physics; they are functions of y ≡ (qb/2) 2 , where b is the nuclear size. The sum over k in Eq. (2.4) involves the enlarged set of five independent and two interfering responses: the standard SI M , the transverse SD Σ , the longitudinal SD Σ , the angular momentum dependent (LD) ∆ and the angular momentum and spin dependent (LSD) Φ whereas interference of the angular momentum and spin dependent response with the standard SI leads to Φ M and the transverse SD and the angular momentum response interfere to give rise to ∆Σ . Five of the seven responses are accompanied by a factor of q 2 /m 2 N , a parameter related to the relative velocities of nucleons bound in the nucleus, further reflecting that the NREFT operators associated with the responses are all suppressed by at least q 2 /m 2 N . If the Wilson coefficients in the final non-relativistic effective action are all of the same order of magnitude, the NREFT operators listed in Table I are not equally relevant for scattering. The traditional SI operator O 1 then dominates over the rest since it is enhanced by a factor of A 2 , where A is the nucleon number, and has no velocity suppression. Naively, one would conclude that the other remaining momentum/velocity independent operator O 4 would be the dominant contribution to scattering in the absence of O 1 ; however, that is not always the case. Recall that the momentum transfer involved in elastic DM-nucleon scattering can be of the order a hundred MeV, | q| ∼ 100 MeV. Once the target nucleus can be resolved into its constituent nucleons, there will be contributions which are only suppressed by the normalized momentum transfer q 2 /m 2 N , which can be O(10 −2 ), rather than by the squared velocity v 2 T ∼ O(10 −6 ). Therefore the contribution from momentum transfer suppressed NREFT operators frequently dominates over those from velocity suppressed ones, at least for relatively large recoil energies. As a result, in the absence of O 1 the SI P-and T-odd operator O 11 often dominates scattering instead of O 4 , if the corresponding Wilson coefficients are of similar size. This explains why numerical scans [21,22] found fairly stringent constraints on the Wilson coefficient of these operators.

III. DM SIMPLIFIED MODELS
Simplified models of DM provide a minimal framework to explore its phenomenology across direct, indirect and LHC searches. In these models the SM is augmented with a WIMP candidate and a mediator via which quarks and/or leptons can interact with WIMPs. Although renormalizable and invariant under SU (3) C × U (1) em gauge transformations, they usually ignore the SU (2) part of the SM gauge group; they are therefore not UV complete theories. However, since simplified models contain a relatively small number of parameters they can be powerful tools to interpret experimental results and explore the complementarity of different searches.
Ref. [36] has constructed a set of simplified models by exhaustively listing different WIMP and mediator spins and matched the resulting relativistic Lagrangians to the set of non-relativistic operators listed in Table I. In these models the WIMP is an SU (3) C × U (1) em singlet; the interactions are further constrained by the requirement that the WIMP cannot decay.
Here we focus on two simplified models with color-triplet mediators which in the NR limit give rise to the leading order P-and T-odd operators O 11 , O 10 and O 12 ; O 11 is spin-independent while O 10 and O 12 are spin-dependent. In the following subsections we describe these models and use the experimental upper bound on the neutron EDM to constrain their direct detection prospects.

A. Model I
Model I contains a complex spin-0 WIMP S and one or more heavy quark-like mediator(s) Q, both of which are odd under a discrete symmetry Z 2 to forbid dark matter decay. The SM quarks q can have two new Yukawa-type interactions with the WIMP and the mediator(s), via scalar couplings y q 1 and pseudoscalar couplings y q 2 , where the superscript denotes the quark flavor index q. The Lagrangian thus reads: U (1) em invariance implies that at least two mediators are required if the WIMP S is to couple to both up-and down-type quarks. Moreover, a mediator Q k with nonvanishing couplings to quarks of different generations will give rise to flavor changing neutral current (FCNC) processes at one-loop order. For example, a Q k coupling to both d and s quarks will contribute to K 0 − K 0 mixing via box diagrams with Q k and S in the loop. We will therefore assume that each mediator Q q couples only to a single SM quark q in the mass basis, so that the Yukawa coupling matrices appearing in eq.(3.1) reduce to diagonal matrices. The matrix element for s-channel scattering S(p 1 ) + q(p 2 ) −→ S(k 1 ) + q(k 2 ) is then given by: Here we have neglected the mass of the incoming quark as well as terms of order | q| 2 in the denominator of the Q q propagator, but we have kept these terms in the numerator; note that at tree-level only scattering on u, d and s quarks contribute to WIMP-nucleon scattering. Moreover, we have used the Dirac equation and 4-momentum conservation to write the matrix element in a form that is symmetric in WIMP momenta. This facilitates the construction of the corresponding relativistic effective Lagrangian, which results when the mediator Q q is integrated out. Defining the Hermitean derivative on the complex scalars as The subscripts on the Wilson coefficients denote the NREFT operator that the corresponding relativistic effective operator reduces to and the superscripts denote the mass dimension of the corresponding relativistic effective operator.
TABLE II: Non-relativistic reduction of effective operators in Model I. f N T q , N N q , ∆ N q and∆ N are coefficients arising due to promoting quark bilinear to nucleon bilinears [42,43]. We use the values as given in the Appendix of Ref. [36]. Table II   Qq m S ) contributions, the contribution ∝ (S † S) (qq) scales as the difference of the absolute value squared of the Yukawa couplings |y q 1 | 2 − |y q 2 | 2 whereas that ∝ i S † ← → ∂ µ S qγ µ γ 5 q scales as the sum of the absolute value squared of the Yukawa couplings |y q 1 | 2 + |y q 2 | 2 ; the latter contribution is suppressed by a relative factor m S /m Qq . The scalar-pseudoscalar operator (S † S) (qiγ 5 q) reduces to the q suppressed P-and T-odd SD operator O 10 and i S † ← → ∂ µ S qγ µ γ 5 q reduces to the v ⊥ suppressed P-even, T-odd SD operator O 7 .
O 1 is the leading order operator, because it does not suffer from any velocity suppression. It therefore contributes dominantly when compared to O 7 and O 10 , unless its Wilson coefficient is suppressed by a factor of 10 −3 or less. It should be emphasized that this suppression should occur for scattering on both neutrons and protons. 4 Since the relevant hadronic matrix elements are different for neutrons and protons, one cannot arrange both cancellations with only a single mediator. In other words, the P-and T-odd operators can only be significant in this model if one introduces at least two mediators. Suppressing all contributions from O 1 would then require two relations between the twelve couplings y q 1 , y q 2 and six masses m Qq to hold to better than one part in 10 3 . On the other hand, if one introduces mediators with charge 1/3 and with charge 2/3, one can require the Lagrangian (3.1) to respect strong isospin invariance. In this case the cancellations for protons and neutrons are in fact (almost) the same. In the most symmetric set-up where each SM quark q is assigned a mediator Q q , the couplings y q 1 ≡ y 1 and y q 2 ≡ y 2 are flavor-universal, and all mediators have the common mass m Q the contributions from O 1 vanish if Here f N T ≡ q N |qq |N is the contribution of light and heavy quarks to the nucleon mass (scalar nucleon bilinear) and N N ≡ q N |qγ µ q |N is the number of valence quarks in the nucleon (vector nucleon bilinear). For our assumption of flavor-universal couplings, f N T is given by where the heavy quark contribution is due to the trace anomaly of the energy momentum tensor [44]. 5 The vector nucleon bilinear coefficient N N is 3 for both neutrons and protons reflecting the number of valence quarks. Using the values for f N T q given in Ref. [45], the ratio of the nucleon bilinears appearing in Eq. (3.4) is For a detector containing a single isotope one would only need to suppress the coefficient c 1 for the specific linear combination of neutrons and protons determined by the target. However, the Xenon detectors, which currently give the best bounds, contain several isotopes, and unsuppressed O 1 constraints from Germanium are still significantly stronger than O 10 constraints from Xenon. 5 The heavy quark contributions to f N T in Eq.(3.5) assumes m Q mq so that the mediators can be integrated out consistently, leaving only SM quarks behind. This may not be a good approximation for the top quark and its mediator. For flavor-universal couplings the contributions from heavy SM quarks is in any case negligible. However, one might also entertain the possibility that the |y q i | 2 scale ∝ mq, so that the Q − S two-point function corrections to the SM quark masses scale like mq; in this case the contributions from heavy SM quarks would be comparable to that in models where the WIMP interacts with nucleons via Higgs exchange.
Note that the stability of the WIMP S requires m S ≤ m Q , hence eq.(3.4) relates the two Yukawa couplings via an O(1) factor. From the simplified model point of view there is a priori no reason why Eq.(3.4) (or a suitable modification thereof) should hold, i.e. in almost all of the parameter space the contribution from the traditional operator O 1 will in fact dominate. 6 Imposing Eq.(3.4), O 10 remains as dominant operator since its contribution to the scattering matrix element is only suppressed by q mn , as compared to a suppression v ⊥ in O 7 , as discussed in Section II. The quark EDM, d q , can be calculated as the coefficient of a dimension-5 P-and T-odd interaction term (−i/2)qσ µν γ 5 q F µν at vanishing momentum transfer. Since the mediator Q is charged under SU (3) C , it can couple to gluons, hence non-vanishing chromo-quark EDMs might also be generated. These are calculated similar to quark EDMs, by finding the coefficient of (−i/2)qσ µν t a γ 5 q G µν a at vanishing momentum transfer. From the left diagram in Fig. 1 we calculate the quark EDM to be where Q q is the electric charge of the quark q, and hence of the fermionic mediator (in units of the proton charge), and the loop function F (m 2 q , m 2 S , m 2 Q ) is given by The chromo-EDM,d q can be obtained by replacing the external photon with a gluon as in Fig. 1a. The chromo-quark EDM is thus also given by Eq. (3.7) after replacing eQ q with the strong coupling g s . In the limit m q → 0 the loop function (3.8) varies between 1/(2m 2 Q ) for m 2 S m 2 Q and 1/(3m 2 Q ) for m 2 S = m 2 Q .

B. Model II
Model II introduces a fermionic WIMP χ and one or more complex spin-0 mediator(s) Φ, both odd under a Z 2 symmetry which stabilizes χ. This model is similar to Model I, but the spins of the WIMP and the mediator have been exchanged, i.e. Model II resembles a supersymmetric model where χ is the lightest neutralino and the mediators are squarks. The renormalizable SU (3) C × U (1) em invariant Lagrangian is given by We have again assumed the two new types of Yukawa couplings, l 1 and l 2 , to be flavor-diagonal in order to avoid new one-loop contributions to FCNC processes due to mediator-WIMP loops. Of course, the electric charge of the mediator Φ q must again be the same as that of the corresponding quark q. Moreover, complex phases in the new Yukawa couplings will again lead to CP violation. The matrix element for WIMP-quark scattering via Φ q exchange in the s-channel is now given by Applying Fierz rearrangement identities from Appendix B to the spinor quadrilinears appearing in eq.(3.10), a number of relativistic effective operators are generated. They are listed in Table III, together with their matching to the NREFT operators in terms of the parameters of this model. 6 The contribution to the Wilson coefficient of O 1 which is suppressed by m S /m Q has not been included in [36]. If this contribution is neglected, the Wilson coefficient vanishes for |y 1 | = |y 2 |, which could be motivated via chiral symmetry. However, an otherwise unsuppressed contribution from (m S /m Q )O 1 would still dominate over contributions from O 7 and O 10 unless m Q > ∼ 10 3 m S in which case the scattering cross section is anyway very small. The left column gives the relativistic 4-fermion operators that appear in the matrix element for WIMP-nucleon scattering, and the right column gives the corresponding NREFT operators including the coefficient appearing in this scattering matrix element. f N T q , ∆q N , N N q , ∆ N q and δ N q are coefficients arising due to promoting quark bilinear to nucleon bilinears [42,43]. We use the values as given in the Appendix of Ref. [36].
We see that Model II generates the following NREFT operators: are all spin-dependent but are suppressed by at least one factor of v ⊥ ∼ 10 −3 or q/m N . They can therefore be neglected relative to O 4 , which did not contribute in Model I. Moreover, O 8 is suppressed by v ⊥ ∼ 10 −3 and can thus also be neglected. As in Model I the dominant term will come from O 1 unless its coefficient is tuned to be tiny. The contributions from O 4 and O 11 might be roughly comparable: the latter is suppressed by a single power of q/m N , but enhanced by the coherence factor A since it does not depend on the spin of the nucleus.
As in Model I, there are two contributions that reduce to O 1 : the scalar-scalar bilinearχχqq contributes proportional to the difference of the absolute squares of the Yukawa couplings |l q 2 | 2 − |l q 1 | 2 whereas the vector-vector bilinear The total coefficient of O 1 can therefore again be made to vanish by an explicit cancellation. Assuming universal couplings and masses for all mediators, the cancellation condition is This relation depends on the same ratio of nuclear matrix elements as the analogous relation (3.4) for Model I. We note that these universality assumptions forbid to also tune the coefficient of O 4 to zero. This could be arranged for general couplings; however, if the new interactions violate strong isospin, the coefficients of O p 1 and O n 1 would have be tuned to zero separately.
We see from Table III that the coefficient of the leading new operator O 11 is proportional to the relative phase between the two new Yukawa couplings; this is true also for the coefficients of the other two P-and T-odd operators, O 10 and O 12 . Since this phase signals CP violation, we expect it to show up in the the quark EDM and CEDMs. The relevant diagrams in this setup are depicted in Fig. 1b. The resulting quark EDM is given by Here Q q is again the electric charge of q, which is equal to the charge of Φ q , and the loop function G(m 2 q , m 2 Φ , m 2 χ ) is given by Again, the quark CEDM can be obtained by replacing eQ q by the strong coupling g s in Eq. (3.12). The EDMs in Eq. (3.7) and (3.12) are proportional to the mass of the fermion running in the loop, since the dipole operators violate chirality. Moreover, we see the same CP-violating combination of couplings as in the coefficient of the P-and T-odd NREFT operators. The loop function G varies from 1/(2m 2 Φ ) for m 2 χ m 2 φ to 1/(6m 2 Φ ) for m χ = m Φ . We note that our loop calculations in (3.7) and (3.12) are in complete agreement with earlier results in the literature, for instance the model independent calculation of EDMs of ref. [47].

IV. RESULTS AND DISCUSSION
We had seen in the previous Section that the same combination of Yukawa couplings appears in the expression of the nEDM and in the coefficients of the P-and T-odd operators contributing to WIMP-nucleon interactions. The upper bound on the nEDM [30] therefore leads to upper bounds on the WIMP-nucleon scattering cross sections due to these operators. In this section we discuss these constraints quantitatively, focusing on O 10 for Model I and O 11 for Model II which are the leading new operators (beyond O 1 and O 4 ) in these models. We compare the resulting constraints with the most recent XENON1T results [1] as well as with the irreducible "neutrino floor" background in Xenon experiments.
We compute the nEDM induced from quark EDM and CEDM dimension-5 contributions which were given in the previous Section. We assume that these are the only contributions to the (color) EDMs of the quarks. In order to calculate the value of the nEDM from quark EDM d q and CEDMd q , we use the following formula: Here the tensor charges g u T = −0.233 (28), g d T = 0.774(66) and g s T = 0.009 (8) have been calculated using lattice QCD [49,50] (see also Refs. [51][52][53]) at a renormalization scale of 2 GeV. We are not aware of a reliable lattice computation of the contribution of the chromo-EDMs to d n ; we therefore employ a computation using QCD sum rules evaluated at a renormalization scale of 2 GeV [54], although there is an O(50%) uncertainty in these results [55,56]. We will again assume that the new Yukawa couplings are universal. In this case the uncertainty in the coefficients in Eq.(4.1) might shift the boundary of the excluded region slightly, without affecting our results qualitatively. We also stress that these results are independent of the cancellation relation (3.4) or (3.11) that one has to impose if the contribution of the new operators is to be significant.
Using the computed expressions for quark EDM and CEDM in Eq.(3.7) and (3.12), we find limits on |Im(y 1 y † 2 )| and |Im(l 1 l † 2 )| respectively according to the stringent experimental bound [30] given in Eq.(1.1). We then convert these limits on the imaginary part of the product of Yukawa couplings to limits on the cross section, using the relations given in Tables II and III as well as the following expression for the "WIMP-nucleon scattering cross section" 7 [57]: We note that due to the chosen normalization of the operators, these expressions look very similar to those for the traditional operators O 1 and O 4 . In order to compare the nEDM-derived 90% c.l. limit with the XENON1T sensitivity, we assume a standard isothermal DM halo with ρ χ = 0.3 GeVcm −3 , v 0 = 220 km/s, v e = 232 km/s and v esc = 544 km/s. We use the Mathematica code dmformfactor [14] for computing the Xenon nuclear response functions W N N k but use our own routines to calculate the exclusion limits using the procedure outlined in Appendix A. This leads to the current exclusion limits on the WIMP-nucleon cross section for O 10 and O 11 at 90 % C.L. from the most recent XENON1T results [1]. We also show the irreducible background level from coherent neutrino-nucleus scattering ("neutrino floor") as estimated in Ref. [48] Figure 2 depicts the σχ n − m S exclusion involving the operator O 10 for Model I with four values of the common mediator mass m Q . We see that the current XENON1T exclusion contour does not go below 10 −2 pb. Since O 10 is a momentum-suppressed SD operator, the limits plotted are roughly q 2 /(m 2 N A 2 ) ∼ O(10 −8 − 10 −9 ) suppressed in comparison to the SI results quoted by the XENON1T collaboration.
In order to derive the constraint on the scattering cross section that results from the bound on the nEDM, one can replace the combination Im(y 1 y † 2 )m Q by a constant (proportional to the nEDM) divided by the loop function F of eq.  Fig. 2 only differ in the end point of the x axis, which is given by m Q . For m S m Q the constraint becomes somewhat weaker, partly because the loop function F becomes smaller, as we saw above, but mostly due to the 1/(m 2 Q − m 2 S ) 2 factor from the squared Q propagator in the scattering cross section. We see that this indirect constraint is more than twelve orders of magnitude below the present XENON1T sensitivity, and at least eight orders of magnitude below the neutrino floor. In other words, in this scenario the nEDM constraint implies that the contribution from O 10 to WIMP-nucleon scattering is totally negligible. Figure 3 shows analogous results for the O 11 operator in Model II. O 11 is independent of the nuclear spin but momentum suppressed; in fact, for small q = | q| the suppression must be of order q 2 /m 2 A rather than q 2 /m 2 N since a pointlike (unresolved) nucleus does not "know" about its nucleonic constituents. The XENON1T limit on the cross section is therefore weaker by three to four orders of magnitude compared to the case of the traditional SI operator O 1 .
Regarding the indirect constraint on the scattering cross section from the nEDM bound, the situation is slightly more complicated here than in Model I because according to Table III there Table III, i.e. the nEDM constraint on |c 11 | 2 scales like 1/m 2 χ . Since now µ χN m N the nEDM constraint on the cross section again scales ∝ 1/m 2 χ . Since the two contributions to c 11 have opposite sign, there is a cancellation for m χ ∼ m N , leading to a very strong upper bound on the cross section.
For m χ m Φ the nEDM constraint becomes a bit weaker again, chiefly due to the 1/(m 2 Φ − m 2 χ ) 2 factor from the squared mediator propagator. Altogether we nevertheless see that the nEDM constraint is again several orders of magnitude below the neutrino floor, i.e. the contribution from O 11 to the WIMP-nucleon cross section can be neglected.
Although we have shown above that for Model II, O 10 and O 12 are expected to be sub-leading due to the presence of O 4 and O 11 , in Fig. 4 we repeat the exercise of   Table III shows that in Model II the coefficient of O 10 also receives two contributions, one of which scales like m N /m χ ; hence the dependence of the nEDM constraint on m χ is similar to that in Fig. 3. However, since the contribution that is not proportional to m N /m χ has a much smaller coefficient, the cancellation now happens at significantly larger WIMP mass than in Fig. 3. On the other hand, O 12 only receives one contribution, which is independent of m χ as long as m 2 χ m 2 Φ . The shape of the nEDM constraint on the scattering cross section therefore resembles that for Model I shown in Fig. 2. Evidently in all cases the nEDM constraint once again implies that the contributions from these P-and T-odd operators to WIMP-nucleon scattering can safely be neglected.

V. CONCLUSIONS
We studied the detection prospects of P-and T-odd operators arising in the NREFT formalism for WIMP-nucleus elastic scattering. These operators occur only if one includes terms that are suppressed by powers of relative velocity or three-momentum transfer. They can therefore be expected to be significant only if the coefficient of the leading spin-independent operator O 1 is very small or vanishes. Moreover, T-odd NREFT operators can only arise in the low-energy limit of a relativistic theory if the latter violates CP invariance. This can give rise to very stringent constraints on the theory, in particular from electric dipole moments.
We analyzed these issues in the framework of two simplified models, taken from ref. [36]. They introduce one or more s−channel mediator(s) carrying electric and color charge, with spin 1/2 (Model I) or 0 (Model II); the dark matter particle then has spin 0 or 1/2, respectively. These models will cause FCNC at one-loop unless each mediator only couples to one quark. Moreover, we pointed out that generically the coefficient of O 1 is not suppressed in these models. Since by now strong constraints on spin-independent WIMP-nucleon scattering exist for several isotopes, the new operators will be significant only if the coefficient of O 1 is suppressed for both neutrons and protons. This generically requires two independent cancellations, which should hold to one part in 10 3 or better. We also note that one has to introduce several mediators for these cancellations to be possible at all, i.e. in minimal versions of these models, with a single mediator, the contribution from O 1 will always dominate.
We therefore assumed that the new couplings and mediator masses are universal for all generations. Since in this case the new terms in the Lagrangian respect strong isospin invariance the two cancellation conditions become almost the same. However, we showed that these scenarios give rise to one-loop contributions to the electric dipole moment of the neutron (nEDM) which are proportional to exactly the same combination of new Yukawa couplings which appears in the coefficients of the new P-and T-odd NREFT operators. The resulting indirect constraint lies well below the irreducible neutrino background. This implies that the contributions from the new T-odd operators on the WIMPnucleon scattering can safely be neglected in the analyses of both current and future experiments. This conclusion can only be avoided if one relaxes the universality assumption on the new couplings and mediator masses. However, one will then have to impose three independent conditions, in order to cancel the coefficients of O 1 for protons and neutrons as well as the neutron EDM. Since these cancellations are not enforced by any symmetry, we conclude that it is very unlikely indeed that the new T-odd operators can contribute significantly. We emphasize that, for a given size of the Wilson coefficients, these are the most important of the total eleven new operators. This indicates that at least in models with charged mediator coupling to light quarks, WIMP-nucleon scattering can safely be analyzed using the traditional operators O 1 and O 4 only.
Charged mediators coupling exclusively to heavy (c, b, t) quarks will contribute to WIMP-nucleon scattering only at one-loop level. On the other hand, the neutron EDM constrains the EDMs of heavy quarks only weakly; the nEDM may then receive dominant contributions at the two-loop level, e.g. via Weinberg's three-gluon operator [58]. In this case the contribution to the nEDM would still only occur at one order higher in perturbation theory than WIMP-nucleon scattering, as is the case in the scenarios we analyzed here, so the conclusions are likely to be similar. Finally, the P-and T-odd operators can also be generated from simplified models not containing charged mediators [36]. We leave the analysis of these scenarios to future work.