Direct detection of dark matter polarizability

We point out that the direct detection of dark matter via its electro-magnetic polarizability is described by two new nuclear form factors, which are controlled by the 2-nucleon nuclear density. The signature manifests a peculiar dependence on the atomic and mass numbers of the target nuclei, as well as on the momentum transfer, and can differ significantly from experiment to experiment. We also discuss UV completions of our scenario.


Motivations
Dark matter (DM) with non-vanishing couplings to ordinary matter may be probed in underground direct-detection experiments. Such couplings can arise from short-range interactions with protons and neutrons, or via weak interactions with photons.
The latter are particularly relevant whenever the DM field X directly couples to messengers that carry electro-weak charges, but couple only very weakly to gluons, the standard model (SM) fermions, and the Higgs boson. In these scenarios, complex DM with spin will generically acquire electro-magnetic dipole moments. These lead to very large direct detection (DD) signatures, and have already been studied by many authors [1] [2][3] [4][5] [6]. Here we are interested in the alternative scenarios with self-conjugate X (real scalar, Majorana fermion, real vector, etc.) in which DD is controlled by the DM electro-magnetic polarizability. We define the latter according to and a similar operator with F µν F µν replaced by µναβ F µν F αβ , with F µν the photon field strength. 1 For definiteness we assumed X is a Majorana fermion, but our results apply to self-conjugate DM of any spin. Self-conjugate DM also couples to photons via the anapole operator Xs µ X∂ ν F µν , with s µ the DM spin. This has lower dimensionality and generically dominates over O γ unless additional assumptions are made. Since the anapole violates separately C and P , while O γ does not, a natural way to suppress its effects is to assume that the dark sector approximately respects either C or P . Similarly, O γ = µναβ F µν F αβ iXγ 5 X can dominate over the anapole if the dark sector is approximately C/P invariant. We thus conclude that, under reasonable and generic conditions, the couplings of self-conjugate DM to photons are controlled by O γ , O γ .
Direct detection via DM polarizability was first studied in [2] in the limit in which the interaction is described by a DM wave propagating in the electro-magnetic field of an infinitely heavy target nucleus. Later, the authors of ref. [7] emphasized that DM scattering for arbitrary masses proceeds via a photon loop, and estimated the rate using an effective field theory for the nucleus. More recent work on the DD signatures of (1.1) can be found in [8] and [9].
In this paper we present a detailed analysis of the DD signature induced by O γ and O γ . After a qualitative discussion of the nucleon/target effective field theory (EFT) in section 2, our main results for O γ are presented in sec. 3. A numerical study in section 4 emphasizes the unique nature of the corresponding DD signature. A comparison between our results and the existing literature is given in Appendix A. The operator O γ is discussed in section 5. In section 6 we emphasize some important features that characterize UV complete models with unsuppressed O γ , O γ , and comment on the coupling XXH † H. A summary of our results is presented in section 7.

EFT at a qualitative level
We start with an analysis of O γ , whereas O γ will be discussed later on.
There are two types of direct detection signatures that (1.1) can lead to: an elastic scattering XT → XT (here T stands for the target nuclei) or an inelastic process XT → XT γ [7]. The first is numerically a loop effect. The latter process arises at tree-level, but its rate is suppressed at least by a factor v 2 4π/α ∼ 10 −3 (v is the incoming DM velocity) compared to the former, and is therefore completely negligible.

The nucleon Lagrangian
The rate for the elastic scattering XT → XT may be found exploiting the hierarchy of scales with q the momentum transfer, 1/Q 0 the radius of T , m N the nucleon mass, and m T the target mass. One first performs the RG evolution from the new physics scale ∼ Λ to the scale ∼ m c . Here one finds that O γ mixes with the quark mass operators O q = qHqX 2 at one loop, and the latter with the gluon operator O G = αs 4π G 2 µν X 2 via an additional QCD loop. In addition, one should take care of the top and bottom quark thresholds. Once this is done, the EFT at leading order in q/m c reads i=γ,u,d,s,G with α = e 2 /4π the fine structure constant. In section 6 we argue that the natural expectation in realistic models is The Wilson coefficients C q,G,γ (m c ) can be calculated using standard perturbation theory (see [8] [9] for a discussion of the case C q,G (Λ) = 0). Alternatively, one can derive the leading non-derivative interactions of X by simply observing that (1.1) renormalizes the QED gauge coupling. By a formal redefinition (A, e) → (A eff , e eff ), with we can remove X from the Lagrangian (up to momentum-suppressed terms). The EFT at the lower scale is now a function of e eff (X), whereas by gauge invariance eA = e eff A eff does not depend on the DM. This trick for example implies in agreement with an explicit loop analysis. Next one should match the quark EFT onto a theory for the nucleons N = n, p. The leading DM couplings now are: It is understood that all couplings and operators are renormalized at ∼ m N . The remainder O(q/m N ) also includes the chiral corrections discussed in [10] [11]. Importantly, the coefficients C p,n (m N ) receive, besides the familiar contributions from O u,d,s,G (m c ) (see [12] for a recent NLO analysis), also a correction induced by O γ (m c ) of order: This latter RG effect can be seen, for example, proceeding along the lines discussed around (2.3). (We emphasize that both proton and the neutron masses are corrected by QED at 1-loop, so that C n (m N ) is also affected despite the neutron has no net charge.) The crucial difference compared to the RG evolution at higher scales is that now the analog of eq.(2.4) is violated by non-negligible higher derivative operators of order m 2 N /Λ 2 QCD ∼ 1. In terms of a heavy baryon EFT these higher-derivative operators correspond to O(α) corrections to the nucleon masses, and more generally the two-nucleon Lagrangian: their main effect is a modification of the pion-nucleon coupling at the percent level. Unfortunately, with our current knowledge of QCD we cannot determine the Wilson coefficients of these operators, and thus C p,n (m N ), with an accuracy better than O(1), even under the assumption (unlikely, according to section 6) that C q,G (Λ) = 0.

EFT for the target nucleus
To determine the scattering rate for the process XT → XT one can proceed in two equivalent ways. The first is based on an EFT for the target nucleus defined at scales ∼ Q 0 , and will be qualitatively discussed in this subsection. The second, which is the one we will adopt in this paper, will be analyzed in section 3.
At scales µ Q 0 m N the target nucleus T is effectively a point-like particle of mass m T Q 0 and one should be allowed to use a heavy nucleus Lagrangian. Up to O(q/Q 0 ), the EFT at µ ∼ Q 0 includes O γ as well as the contact operator where we ignored numerical coefficients for simplicity. The contact operator mixes with O γ at one-loop under the RG, as seen from arguments completely analogous to those discussed above. The terms of order Z, A also receive corrections from C N (m N ) in δL m N . From (2.7) one immediately reads a short distance contribution to the amplitude for XT → XT . There is also a correction coming from a UV-sensitive one-loop diagram involving O γ , which scales as the O(Z 2 ) term in (2.7) [7]. The two contributions are individually scheme-dependent; only their sum is physical. For example, using a massindependent renormalization scheme the loop diagram vanishes at q = 0, and the O(Z 2 ) effect comes dominantly from the counterterm (2.7). 2

The 2-body process
The approach followed in section 2.2 is intuitive from a physical standpoint, but not very convenient. One reason is that it depends on several unknown parameters, even in the optimistic (and unrealistic) case in which only O γ is present at µ ∼ m N . More importantly, though, it obscures the accuracy of the perturbative expansion. For example, are we allowed to ignore QED vertex corrections to the one-loop diagram of [7]? These are naively of order αZ 2 /4π, and apparently not negligible for heavy targets. Feynman diagrams for the 1-proton and 2-proton processes. Both diagrams contribute to the contact DM-nuclei interaction, whereas the one on the right is also related in a schemedependent way to the one-loop diagram of [7] (see Appendix A for details).
In this section we will approach the problem from the point of view of the "fundamental" nucleon EFT. In practice we take (2.5) as our starting point, derive a multi-body effective theory for the nucleons, and finally take the appropriate nuclear matrix element. Using this formalism all the unknowns will be encoded in measurable nuclear form factors. Furthermore, within this formalism the perturbative expansion becomes manifest. For instance, an inspection of the O(αZ 2 /4π) "vertex corrections" mentioned above shows that these are secretly a renormalization of the nuclear wave-function, and hence already included in the nuclear potential.

The 2-proton form factor
We now want to calculate the amplitude for XT → XT from the nucleon Lagrangian (2.5). The discussion in section 2.2 shows that this process receives contributions from the operators O n,p in (2.5) as well as O γ (m N ). The former may be treated using standard methods. Our main focus here will be on O γ . In section 4 we will study in detail the interplay between all Wilson coefficients C γ,p,n .
Working at leading order in the DM operator we find that the dominant DM-nucleon diagrams contributing to XT → XT are those shown in fig. 1. The loop diagram on the left vanishes for q = 0 in any mass-independent renormalization scheme, which is the natural regulator in our EFT. Therefore, only the 2-body process in the right of figure 1 is relevant at leading q/m N order. The process is tree-level when using, say, MS.
The non-relativistic amplitude for the 2-proton process ppX → ppX is where q i,j are the three-momentum transferred to the nucleons, whereas s i,j,i ,j are spin indices for the nucleons (the DM spin indices are not shown because they cancel out in the cross section when summing and averaging over final and initial states). We did not add the crossed diagram because it will be automatically included when convoluting M 2 with the anti-symmetric nuclear wave function. Following [10][11], we used a nonrelativistic normalization for the 1-particle states. In practice, this corresponds to divide the relativistic amplitude by (2m N ) 2 (2m X ). With this convention, the formula (3.1) also applies to real scalars (with C γ a parameter with dimensions of a mass).
The DM-nucleon potential in "mixed coordinates" (spacial for the nucleons and Fourier for X) reads: The term linear in q arises because the transition is mediated by a massless particle. The contribution of the 2-body process to the amplitude for XT → XT , with the target remaining in the ground state, is given by [10][11] Hereρ (2) is the (diagonal) 2-proton nuclear density matrix, the sum extends over all proton pairs, and ⊗ indicates a contraction of the spin indices. In the second line we used the fact that protons are indistinguishable. The factor Z(Z − 1) signifies that this is truly a two-body effect, and as such it vanishes for Z = 1. The relevant spin-singlet quantity is the projection: Finally, employing eq. (3.3), and making a trivial coordinate transformation, we arrive at an expression for the dominant, spin-independent part of the amplitude: where we defined the 2-proton form factor 3 Eqs. (3.6) and (3.7) are one of the main result of the present paper. Before moving on it is important to establish that we have a consistent perturbative expansion. First, the relative size between two-body and one-body contributions induced by O p,n is ∼ ZQ 0 m N (see also eq.(2.7)), which is not negligible for heavy nuclei. On the other hand, consistently with what done in (2.5) we will neglect corrections O(q 2 /m 2 N ) to both the one-body and two-body terms because of order Q 2 0 /m 2 N = few %. 4 Furthermore, higher order corrections to the diagrams in figure 1 fall into two classes: corrections to the external nucleon legs and to the vertex. The former are encoded into the nuclear wave-function and are not relevant to us. The vertex corrections that leave the order n of the graph unaltered ("n" from "n-body") are down by q 2 /m 2 N or α/4π and are thus negligible. There are also diagrams that increase the order by attaching a pion or a photon to the internal proton line in the one-loop diagram (see left of fig. 1) and a nucleon pair on the other end. These are respectively of order q 2 m 2 N × A Z and α 4π compared to the dominant 2-body, and can also be neglected.

The role of proton-proton correlations
The 2-proton density F (2) appearing in (3.7) plays an important role in many nuclear reactions, such as electro-disintegration processes, (e, e N ) and (e, e N N ). In [10][11] [13], the relevance of 2-body densities in DM detection has been pointed out, although as a subleading effect in the chiral counting. In this paper we find another interesting application for DM direct detection, where the two-body term plays a dominant role.
The exact form of F (2) is not known, but some of its basic properties can be qualitatively understood. Without loss of generality we write is the (one-proton) charge form factor, and F corr a measure of the correlation between the two nucleons (protons in our case). From the normalization of the density distributions follows that F corr (0, 0) = 0. In practice this means that one may neglect the correlation when both |q|, |k| are much smaller than, say, the pion mass m π . As q 2 , k 2 ∼ m 2 π , nucleon-nucleon correlations become non-negligible and the approximation F corr = 0 is violated. At even larger momenta one anticipates a universal shape for F (2) , dominated by the repulsive pion exchange.
From this simple consideration follows that F corr cannot be ignored in general, because the integral in (3.7) probes a regime where the correlation is presumably non-negligible. In particular, F corr (k, 0) is likely to result in an overall O(1) correction in the nuclear matrix element. To see this more explicitly, we take the following phenomenological expression for the charge form factor and parametrize F (2) with where c i are order one numbers (in general functions of A, Z) mimicking the effect of a short distance correlation. The expression (3.7) can now be solve exactly: In the second line we used the expansion of f given in (3.3). As expected, we see that c 1 changes the overall rate, whereas both c 1,2 modify the momentum-dependence of the form factor (c 3 does not contribute for our choice of F (2),pheno ). We will present a numerical study in the next section. Interestingly, when c i → 0 the form factor F pheno pp reduces to that derived in [7] [8]. This correspondence is elucidated in Appendix A.
We conclude this section observing that the proton-proton correlation becomes parametrically small if one models the nuclear potential with a mean field approximation. Indeed, in that case the nuclear wave-function is described by a single Slater determinant, and it is a trivial exercise to show that F corr = O(1/Z) is entirely due to the Pauli exclusion principle. (We also numerically verified this expectation using a shell model.) In reality the mean field potential is no more than an intuitive picture of the nucleus, and F corr cannot be neglected.

Signatures in direct detection experiments
The study of the DD signature induced by O γ is complicated by two obvious hurdles. First, we do not have a reliable estimate of the 2-proton form factor (3.7). Second, large hadronic uncertainties make it impossible to precisely determine the actual relation among the Wilson coefficients in (2.5) and the fundamental parameters.
We thus take a bottom-up, phenomenological approach. For the form factors we use the phenomenological expressions given in (3.9) and (3.10). The spin-independent differential rate is then given by We simplified eq.(4.1) assuming that protons and neutrons have the same mass and 1particle densities. We use Q 0 = 0.5(0.3 + 0.9A 1/3 ) −1 GeV. Now, the discussion in section 6 (see also section 2.1) suggests that, in generic theories with unsuppressed C γ , the natural expectation is C N /C γ ∼ α/π, and hence f pp /f p , f n /f p = O(1). In the following we will therefore treat these latter as independent parameters of order unity. Because f pp already accounts for changes in the overall normalization of the two-body effect, it makes sense to work with form factors that are normalized to one at q = 0; for this reason we will set c 1 = 0 and vary only c 2 in our numerical analysis (c 3 has no effect on (3.10)). Comparing to [8], our phenomenological 2-body form factor for c 1 → 0 reduces to F pheno pp → F Ray (1 + c 2q 2 ) (see Appendix A).

Numerical analysis
The impact of the new form factor may be significant, both in the spectrum and the total rate. When both f pp /f p , f pp /f n have the same magnitude and sign, the new term f pp is expected to dominate for heavy nuclei (at least as long as the form factors are positive and non-vanishing). However, an opposite sign in either f pp /f p , f n /f p can generically result in destructive interference among the various contributions to (4.1), and hence lead to qualitatively new effects. 5 Because of the different momentum-dependence of the form factors, the suppression in dσ T /dE R will occur at a specific recoil energy. This energy is a strong function of the parameters f pp /f p , f n /f p and the form factors, as well as of the target, due to the peculiar dependence on A, Z.
To assess the effect of f pp /f p , f n /f p on the total rate, the relevant quantity to consider is the number of events within a certain signal region ∆: Here ρ X is the local DM density, Ex the detector exposure, N T the number of nuclear targets (summed over all isotopes), η eff (s)p(E R , s) the signal efficiency, and f lab the DM velocity distribution in the lab frame. While for f pp = 0 the surface N ∆ = const is defined by parallel lines f n ∝ f p , in the more general case it becomes an ellipses. To see this one can focus on a single isotope, in which case, denoting by the integrals in (4.3), it is clear that for x ∝ f pp /f p and y ∝ ( Z A−Z + f n /f p ). By the Cauchy-Schwarz inequality the last term is always positive, implying that the isocurves must be ellipses. We present our numerical results for three experiments: Xenon100, LUX and CDMS-Ge in Figures 2, 3. For the experimental details of the analysis we refer the reader to [13] and references therein. The color coding is Xenon100 (red), LUX (green), CDMS-Ge (blue).
The impact of f pp is most conveniently discussed in terms of the ratio is defined with f pp = 0, f n = f p . This is sometimes referred to as the "degradation plot" in the literature. In Figure 2 we present degradation plots for f pp /f p = −1 (short dashed lines), f pp /f p = 0 (solid lines) and f pp /f p = 1 (long dashed lines). We consider the case of light DM, m X = 10 GeV (left column), and m X = 100 GeV (right column). The top, center and bottom plots correspond to variations of the parameter c 2 (see our phenomenological twobody form factor) from −1, 0, to +1.
For the case f pp /f p = 0 we find the familiar LO degradation plot (see for instance [14] and references therein). By turning on non-zero values of f pp /f p , i.e. the two-body term, the entire plot is shifter and the minimum increases. This can be qualitatively understood looking at (4.4), from which we see that the minimum of the rate is always found at f pp = 0 (x = 0), and gets displaced when f pp = 0.
Importantly, for f pp = 0 the Xenon100 and LUX experiments (and more generally other Xe-based detectors) observe different rates, despite having the same target nuclei. The reason is that the rate depends on the integrals in (4.3), which are themselves functions of the energy threshold and the associated efficiency. If DM is detected, in principle this feature may be used to infer important information about the parameters f pp /f p , f n /f p . However, Fig. 2 suggests that the difference between the rates at different experiments significantly depends on the two-body form factor, so an accurate determination of F (2) would be required to draw any conclusion on f pp /f p , f n /f p .
To better appreciate this, in Figure 3 we show contour lines of N ∆ /N ref ∆ = 0.01 (similar contours are obtained for different values, as clear from Figure 2) as a function of f pp /f p , f n /f p . As anticipated, we see that the amount of overlap among the (green, red, and blue) ellipses is critically sensitive to the two-body form factor.
In Figure 3 we also observe a qualitative difference between light and heavy DM. For light DM the ellipses are very elongated, while for heavier masses they are more circular in shape. Again, this can be understood from (4.4). The astrophysical function in (4.3) is more steeply falling for light DM, and this basically forces E R to acquire values close the lowest bin. In the approximation that the integrand in (4.3) is a delta function in E R , the last term in (4.4) is very small and the domain N ∆ /N ref ∆ = const is determined by a line x ∝ y with an experiment-dependent slope. For heavier DM the recoil spectrum has effectively a larger range and the second term in (4.4) becomes important.
The nature of a possible cancellation in the event rate is analogous to that invoked in isospin-violating DM (see e.g [15][16]), as it arises from fine-tuning parameters that are naturally of the same order, but is qualitatively different for at least two reasons. First, the cancellation depends on E R . Second, it involves two parameters -i.e. f n /f p , f pp /f p -rather than one, so in principle it is possible to suppress dσ T /dE R in two experiments simultaneously.

The 2-nucleon form factor for F µν F µν
To study DM scattering induced by F µν F µν one can proceed in complete analogy with the O γ operator: determine the non-relativistic amplitude for N N X → N N X and derive the corresponding nucleon potential.
The crucial difference is that in the present case also neutrons contribute. The matrix element T f |F µν F µν |T i is dominated by the coherent scattering of one of the two photons on the proton charge, as in section 3.1, and the incoherent spin-dependent scattering of the second photon on the nucleon magnetic moment.
The multi-nucleon potential can again be written as in section 3.1: where now f ij depends non-trivially on the nucleon spin. Up to O(q 2 /m 2 N ) we find: where we defined − → µ ± ij = s i µ m i e j ± s j µ m j e i , with s i = σ/2 the nucleon spin operator, and µ m i , e i the nucleon magnetic moment and electric charge. To keep our discussion general, we did not specify the DM bilinear coupling to F µν F µν . This model-dependent contribution (generally momentum-dependent) is included in the coefficient V 0 . In other words, what we want to discuss here is the nuclear matrix element T f |F µν F µν |T i .
The nuclear form factor follows immediately from T f | i,jṼ ij |T i . One of the nucleon indices runs over the proton charge, and results in a single power of Z in the amplitude. The remaining (spin-dependent) sum involves both neutrons and protons and is similar to that found in ordinary spin-dependent interactions.
Similarly to O γ , the RG evolution of O γ will induce contact (spin-dependent) interactions that reduce in the non-relativistic limit to the familiar 1-nucleon potential ∝ q · s i . The ratio between long and short distance contributions to the amplitude for XT → XT scales as ZQ 0 /m N , as we found in section 3. An equivalent way to check this is to follow a logic similar to that of section 2.2, and add to the (spin-dependent) DM-target contact interaction a 1-loop diagram analogous to that evaluated in [7], but now with the nuclear magnetic moment in one of the two γ − T vertices.

Realistic UV completions
In this section we briefly comment on possible UV completions of O γ ( O γ ). We take X to be a Majorana fermion for simplicity, but keep in mind that our results generalize to self-conjugate DM of any spin (baring naturalness issues in the case of scalar DM).

New physics at the weak scale
The existing literature assumes that O γ is generated in isolation at some high scale m * by loops of some heavy mediator with electro-weak charges. However, generic field theories will also induce a lower dimensional coupling to the Higgs mass operator H † H (and the quark mass operator). On the basis of simple dimensional analysis, the Lagrangian will generically contain a term of order (from an additional Z 0 loop) along with C q ∼ α 4π C γ . We are here referring to a contribution at the matching scale m * , not to an RG effect below m * . Now, while C H is suppressed compared to C γ at the cutoff, the latter contributes in DD at one-loop order. A back of the envelop calculation tells us that the ratio between the DM-nucleon cross section induced by O γ and that of O H is parametrically suppressed by (m h /m * ) 4 . Numerically, the coupling O γ will be relevant (i.e. f pp /f p ∼ 1) only if the electro-weak charged mediators have masses m * below a few hundred GeV. This holds irrespective of the actual magnitude of the signal, and is a priori independent from the usual WIMP miracle and arguments based on the naturalness of the SM: generic models in which f pp /f p ∼ 1 will have new physics accessible at colliders. The strongest constraint on the messenger mass scale in these scenarios come from direct searches at LEP and the LHC, and for relatively short lived particles are compatible with new physics at the weak scale. 6 The situation is unchanged in UV-complete models for O γ , since XXH † H is allowed by all relevant symmetries.
It should be clear that, in analogy with the more familiar hierarchy problem, one can evade the natural expectation (6.1) by fine-tuning, or invoking Supersymmetry and/or Higgs compositeness at scales m * . In this sense we claim that a UV completion for (1.1) at the weak scale is generic, but not strictly necessary.

Large DD rates and suppressed indirect signatures
Much can be learnt about possible UV completions by making a rough estimate of the DD rate generated by O γ .
First of all, in a healthy theory (1.1) will be generated at one-loop level with g * a typical coupling (the power of m X is expected if m X m * ). Using the results of section 3 we estimate a nucleon-DM cross section (on Xe) of order with µ N the reduced DM-nucleon mass, and 1/Q 0 ∼ A 1/3 fm a measure of the radius of the target nuclei. Together with the collider bounds mentioned at the end of the previous section, this estimate strongly disfavors DM detection in forthcoming DD experiments. This is even more true for O γ , which has a much lower rate.
However, this estimate fails if the heavy messengers interact with DM primarily via a light scalar. Consider first O γ and postulate the dark sector couples to a light dilaton φ of mass m φ m X , m * . In this case the high energy theory contains a (loop-order) coupling F 2 µν φ and, provided the trilinear XXφ is unsuppressed (typically of order g * m X /m * ), O γ will be dominated by the tree-level exchange of φ. For φ heavier than ∼ 100 MeV the rate (6.3) receives an enhancement of ∼ (m 2 * /m 2 φ ) 2 . Now nucleon cross sections of order 10 −44 cm 2 become plausible for m X , m * ∼ 100 GeV provided m φ /g * is in the few GeV range. Still, natural considerations suggest m 2 φ e 2 m 2 * /16π 2 , from which follows that σ N s significantly larger than a few times 10 −44 cm 2 (100 GeV/m * ) 6 may be interpreted as indirect signature of an unnaturally light scalar.
It is hard to imagine a mechanism to enhance the rate of O γ up to similar values. Even if we add an axion a with mass m 2 a q 2 the cross section will be typically too small to be detected. The reason is that the power of 1/q 2 from the propagator is compensated by a factor q · s DM from the axion-DM coupling, and q · s T from the axion-target coupling. Yet, axion-mediated DM interactions can naturally suppress the "dangerous" DM-Higgs coupling of the previous subsection, and typically have f pp /f p ∼ 1.
An obvious question is whether indirect gamma-ray signatures of, say O γ , can be relevant. Importantly, for m φ m X the s-channel exchange of φ will lead to: Because of the much larger momentum flowing in the φ propagator, the result is a factor ∼ (m φ /m X ) 4 smaller than that obtained assuming that the same scale suppressing O γ in DD experiments also controls indirect signals. Indirect detection rates close to the current sensitivities are still possible (see for instance [21]), especially for scalar DM that does not suffer from the velocity suppression in (6.4).

Conclusions
We presented a detailed analysis of the DD signature induced by the DM polarizability operator O γ , see eq.(1.1), and derived the associated nuclear form factor as a function of the 2-nucleon density. The relevance of 2-body densities was previously pointed out in [10] [11], but in those cases the effect was subleading in the chiral expansion and for generic choices of the parameters. To the best of our knowledge, O γ provides the first known example of DM scattering dominated by a 2-body nuclear form factor. The long range force mediated by the photon is key to our result. The multi-body nature of the interaction implies that the scattering rate varies significantly from experiment to experiment, and that DD experiments with the same target nuclei are expected to measure different rates. The presence of destructive interference in a non-negligible portion of the parameter space makes these scenarios especially interesting in light of current anomalies.
DM scattering via F µν F µν was also discussed. The novel feature here is that the 2-body interaction actually describes an example of coherent and simultaneously spin-dependent DM scattering.
One expects F µν F µν and F µν F µν to be important to DD experiments in a large class of models with self-conjugate DM. We showed that realistic scenarios in which these operators are relevant to DM detection have new physics at scales accessible at colliders. Furthermore, large DD rates are indirect evidence of exotic light scalars, and are typically accompanied by suppressed gamma-ray signatures.
the DM couplings to nucleons was derived. However, our discussion in section 2.1 suggests that this last step is plagued by O(1) uncertainties arising from the RG evolution down to µ ∼ m N , which forced us to treat C p,n (m N ) as "free" parameters. As in [7], the authors of [8] [9] did not include proton-proton correlations. In particular, when c i = 0 the form factor F pheno pp defined in eq.(3.10) reduces to the quantity F Ray introduced in [8]. In general we find: