We provide expressions for the nonperturbative matching of the effective field theory describing dark matter interactions with quarks and gluons to the effective theory of nonrelativistic dark matter interacting with nonrelativistic nucleons. We give the leading and subleading order expressions in chiral counting. In general, a single partonic operator already matches onto several nonrelativistic operators at leading order in chiral counting. Thus, keeping only one operator at the time in the nonrelativistic effective theory does not properly describe the scattering in direct detection. Moreover, the matching of the axial--axial partonic level operator, as well as the matching of the operators coupling DM to the QCD anomaly term, naively include momentum suppressed terms. However, these are still of leading chiral order due to pion poles and can be numerically important. We illustrate the impact of these effects with several examples.

We present the effective field theory for dark matter interactions with the visible sector that is valid at scales of O(1 GeV). Starting with an effective theory describing the interactions of fermionic and scalar dark matter with quarks, gluons and photons via higher dimension operators that would arise from dimension-five and dimension-six operators above electroweak scale, we perform a nonperturbative matching onto a heavy baryon chiral perturbation theory that describes dark matter interactions with light mesons and nucleons. This is then used to obtain the coefficients of the nuclear response functions using a chiral effective theory description of nuclear forces. Our results consistently keep the leading contributions in chiral counting for each of the initial Wilson coefficients.

We provide a Mathematica package, DirectDM, that takes as input the Wilson coefficients of the relativistic effective theory describing the interactions of dark matter with quarks, gluons and photons, and matches it onto an effective theory describing the interactions of dark matter with neutrons and protons. The nonperturbative matching is performed at leading order in a chiral expansion. The one-loop QCD and QED renormalization-group evolution from the electroweak scale down to the hadronic scale, as well as finite corrections at the heavy quark thresholds are taken into account. We also provide an interface with the package DMFormFactor so that, starting from the relativistic effective theory, one can directly obtain the event rates for direct detection experiments.