Threshold Resummation for Dark-Matter Production at the LHC

We derive precision predictions for the production of dark-matter particles recoiling against a jet with large transverse momentum at the LHC. The dark-matter fermions are described within a simplified model and couple to the Standard Model via a vector mediator. Our predictions for the mono-jet signature include the resummation of the leading and next-to-leading threshold logarithms. The corresponding matching coefficient is evaluated at NLO. The resummed result is matched to the fixed-order NLO cross section obtained from the MadGraph framework. We discuss numerical results for several benchmark scenarios at the LHC.


Introduction
The nature of dark matter remains one of the big unanswered questions in modern physics. Weakly interacting massive particles (WIMPs) are well-motivated dark matter candidates and are predicted by many extensions of the Standard Model (SM) [1]. One of the main goals of the Large Hadron Collider (LHC) is to test the WIMP paradigm, see e.g. Ref. [2]. Many LHC searches for WIMP dark matter rely on the mono-X signature, where one would observe a SM particle recoiling against the missing transverse energy caused by WIMPs escaping the detector.
To interpret the LHC data, precise predictions for the dark-matter production cross section and for the corresponding SM backgrounds are essential. The main irreducible SM background is Z-boson production in association with jets, where the Z boson decays invisibly into neutrinos. A thorough discussion of the theoretical description of this dark-matter background can be found in Ref. [3].
Predictions for the dark-matter signal depend on the concrete model under consideration. We employ the simplified-model approach, Ref. [4,5], and derive precision predictions for the process pp → χχ + jet in a simplified dark-matter model with a fermionic dark-matter particle χ and an s-channel vector mediator Y which couples to SM quarks but not to leptons [6]. Hence, stringent bounds from resonance searches in the Drell-Yan process are avoided. The LHC has also placed bounds on the simplified model parameter space employing resonance searches in the di-jet final state [7][8][9].
Mono-X production of dark matter has already been investigated in simplified models at next-to-leading order (NLO) in QCD [10][11][12][13][14][15][16][17]. The simplified dark-matter model with a vector mediator has been analyzed at NLO including jet merging and partonshower effects [13]. The model has been implemented in the MadGraph5_aMC@NLOframework [18] using FeynRules [19,20] to obtain NLO predictions in a fully automated way. The NLO corrections result in a moderate enhancement of the integrated cross section and the missing transverse momentum distribution, depending in detail on the mediator and dark-matter masses. Moreover, a sizeable reduction of the scale dependence of the cross section is reported improving the predictions from LO to NLO, as expected. Given the current WIMP exclusion limits from the LHC [7,8] and from direct-detection experiments, see e.g. Ref. [21], future LHC searches will focus on large dark-matter and/or mediator masses. In this case, a significant part of the cross section can be attributed to the threshold region, where large logarithmic corrections from soft and collinear parton emission beyond NLO should be taken into account.
In this work, we refine the theoretical prediction of Ref. [13] for inclusive observables in the process pp → χχ+jet by resumming potentially large threshold logarithms up to nextto-leading logarithmic (NLL) accuracy. More specifically, we perform the resummation for the double-differential cross section with respect to the transverse momentum and the invariant mass of the dark-matter particles in the final state. To calculate the anomalous dimensions for the renormalization-group evolution of the soft and the jet functions as well as the soft and collinear NLO contributions to the matching coefficient, we use the strategy of regions [22]. The NLL expressions are identical to those presented in Ref. [23] after adapting the color-factors. The corresponding matching coefficient is calculated at NLO in the strong coupling α s . We call these predictions NLL' to distinguish them from NLL predictions with a LO matching coefficient. The hard function, which only enters the matching coefficient, is obtained numerically from MadLoop within the MadGraph5_aMC@NLO-framework [18]. Combining these numerical one-loop results with the analytic results obtained using the strategy of regions provides a convenient and straightforward way to calculate NLL' resummed cross sections. Since the invariant mass of the dark-matter particles is not accessible experimentally, we integrate the resummed distribution over the invariant mass. The resummed predictions are matched to the full fixed-order NLO prediction obtained from the MadGraph5_aMC@NLO-framework [18].
The resummation of threshold logarithms has been achieved for closely related SM processes. In particular, threshold resummation in the direct QCD approach was derived to NLL accuracy for direct photon hadroproduction in Refs. [24][25][26][27][28]. The resummed results in one-particle inclusive kinematics were later used to obtain approximate higherorder corrections for direct photon hadroproduction as well as for the production of a massive gauge boson at high transverse momentum [29][30][31][32]. A general discussion of threshold resummation for single inclusive cross sections at NLL and NNLL can be found in Refs. [26,33,34]. For the hadroproduction of Higgs bosons at large transverse momentum, threshold resummation is known to NLL' accuracy [23]. Recently, a similar analysis for the closely related Drell-Yan process at a fixed invariant mass and high transverse momentum has been performed in Ref. [35]. Within the formalism of Soft-Collinear-Effective-Field-Theory (SCET) threshold resummation to NNLL accuracy has been achieved in Ref. [36] for direct photon production, in Ref. [37] for the production of W, Z bosons, and in Ref. [38] for Higgs-boson production at high transverse momentum.
This work is organized as follows: In Section 2, we briefly introduce the simplified dark-matter model. The kinematics and the partonic channels for dark-matter production in association with a jet are discussed in Section 3. The resummation of the threshold logarithms is discussed in Section 4 where we derive analytical results for the soft and the jet function and calculate the hard function numerically. In Section 5, numerical results for dark-matter production at the LHC are presented for several benchmark scenarios and the effects of resummation are investigated. Finally, we conclude in Section 6 and present the analytic LO cross sections in Appendix A.

Dark-Matter Model
The simplified dark-matter model which we investigate in this work consists of a massive s-channel vector mediator Y and a Dirac dark-matter fermion χ. The interaction part of the Lagrangian in the simplified dark matter model reads The first term in the Lagrangian describes the interaction of the mediator with the dark-matter fermions. We allow for a vector-and an axial-vector coupling. The second term describes the interaction of the mediator with the quark fields q. Again we allow for a vector-and an axial-vector coupling. In addition, we assume that the couplings are flavour diagonal so that no flavour mixing is introduced. Consequently, our simplified model has six free parameters, i.e.
The benchmark points considered in the numerical analysis are given by the fixed couplings and we vary the masses of the mediator and the dark-matter particle. We stick to these values for the couplings which are often used as benchmark points in the literature although large regions of parameter space for the mediator and dark-matter masses are already excluded [7,8]. However, viable models can be obtained by choosing smaller couplings, e.g. g V SM = 0.1 or g A SM = 0.1. As discussed in Section 5, the exact choice for the couplings does not have a relevant impact on the relative size of the NLO or matched NLO+NLL' corrections.

Partonic Channels and Kinematic Considerations
We calculate predictions for the production of a pair of dark-matter fermions χ recoiling against a jet j. At LO, the corresponding partonic processes are where g denotes a gluon, q a quark, andq an anti-quark. In Section 4, we label the different contributions to the NLL' cross section by the final-state parton at LO, e.g. the cross section dσ g refers to the first process.
The central quantity for threshold resummation is the double differential cross section with respect to the transverse momentum p T of the pair of the dark-matter particles in the final state and their invariant mass M 2 , i.e.
where the squared partonic center-of-mass energy is given byŝ = x 1 x 2 s in terms of the squared collider energy s. The parton distribution functions (PDFs) are denoted by f i and the sums run over all contributing initial-state partons i, j. The threshold energy √ŝ T to produce a final state with a given p T and M 2 is given by and the corresponding partonic threshold variable is defined bŷ The LO cross section can be written as a function of the squared threshold variableŷ 2 T and the ratio r = p T /m T . The explicit results for the LO cross sections of the different partonic channels are given in Appendix A. Since the dark-matter particles are not accessible experimentally, we have simplified the otherwise fully differential cross section by integrating over the dark-matter phase space. Hence, we obtain an effective 2 → 2 process which, for a given mass M 2 , has the same kinematics as 2 → 2 SM processes for which threshold resummation has been obtained in the past [23,37]. We perform resummation in Mellin space. The Mellin moments are defined as We do not require an analytic expression for the full NLO partonic cross section in this work. However, it is instructive to inspect its generic structure written in terms of a threshold expansion Here, s H denotes the hadronic mass of the final state, and Q 2 = p T √ŝ T is a high-energy scale of this process which naturally appears in addition tô s T . For initial-state parton combinations i, j present at LO, the z-integration over the plus-distributions results in log 2 (1 −ŷ 2 T ) and log(1 −ŷ 2 T ) contributions at NLO. These threshold logarithms are resummed at NLL accuracy to all orders in this work. The function d ij denotes all terms which vanish in the threshold limit and are not needed at NLL' accuracy. The leading power coefficients of the threshold expansion c ij , and parts of c (1) ij are derived analytically, while the so-called hard contributions to c (1) ij are obtained numerically as discussed in Section 4. Note that This relation can be used in the soft and jet functions derived in Section 4 to relate the two high-energy scales. As long as r is of O(1), i.e. p T is of the same order as the mass of the dark-matter system,ŝ T and Q 2 are of the same order and standard threshold resummation holds without the need to take care of an additional large scale ratio.
As mentioned above, the distribution in M 2 is not accessible experimentally. The transverse momentum distribution or any fiducial cross section is simply obtained by integrating the resummed double-differential cross section, which is derived in Section 4.

Resummed Cross Section
In this section, we outline our strategy to calculate the resummed cross section at NLL' accuracy. For brevity, the double differential cross section introduced in Section 3 is denoted by dσ in the following. We use the method of regions to derive all the one-loop ingredients of the resummation formula. The resummation is performed in Mellin space.
Using the method of regions, there are the soft, the collinear, and the hard region which contribute to the cross section near threshold at leading power. To define the different regions, we decompose the momentum p of emitted radiation into components (n + · p, n − · p, p ⊥ ), where we have introduced light-like four-vectors n ± with spatial components along and opposite to the jet direction (n 2 ± = 0, n + · n − = 2). For the soft region p scales according to (λ, λ, λ), where the expansion parameter λ = 1 −ŷ 2 T is small near threshold. The collinear and the hard regions scale according to (1, λ, √ λ) and (1, 1, 1), respectively.
To obtain the resummed cross section we have to calculate the leading-power result in each region. For the soft and the collinear region, our calculation leads to simple analytical results. The contribution of the hard region is obtained numerically as explained below.
In the soft region, the leading-power expansion of the squared matrix elements for the real-emission contribution is identical to the Eikonal approximation for soft-gluon emission. Using the leading-power expansion also for the phase-space integration, one can write the soft contribution up to NLO for each partonic channel as a convolution with the LO cross section according to where the so-called soft function reads The color factors are given by for the partonic channels, i.e. each initial-state quark contributes C F , an initial-state gluon C A and for each final-state particle the respective color factor is subtracted. Moreover, we work in d = 4 − 2ε dimensions, and µ is the scale introduced in dimensional regularization to keep the running coupling dimensionless. In the calculation, we have used that the jet momentum is perpendicular to the beam direction at threshold. The term (1 − z) −(1+2 ) can be expanded in in terms of a δ-function and plus-distributions. In the following we consider the MS-renormalized quantity At NLO, the soft functions is, by construction, identical to the SCET soft function of [36], where a dimensionful convolution variable is used instead.
Since the Mellin transform converts convolutions like (10) into ordinary products, it is beneficial to go to Mellin space. The Mellin transform of the renormalized NLO soft function reads where we drop all terms suppressed by 1/N , since power-suppressed terms also contribute at O(1/N ), and useN = N exp(γ E ) with Euler's constant γ E . For the Mellin-space functions we drop the label R, since we always consider renormalized quantities. We observe that the characteristic scale of the soft function is given by µ 2 s =ŝ T /N 2 . The virtual contributions to the soft function consist of scaleless integrals and do not contribute.
In analogy, we can consider the collinear region which corresponds to real radiation collinear with the outgoing parton. The leading-power expansion of the squared matrixelement leads to the Born expression multiplied by the relevant splitting function for each partonic channel. Also expanding the phase space, in analogy to the soft case, one finds with the so-called jet function for the outgoing parton a, where for a quark or an anti-quark in the final state and for a gluon, which can also split into N f = 5 light quark flavours.
In analogy to the soft function, the renormalized NLO jet functions for outgoing partons in Mellin space read and We observe that the characteristic scale of the jet function is given by µ 2 c = Q 2 /N . The virtual contributions to the jet functions again consist of scaleless integrals and do not contribute. The calculated jet functions agree with the SCET jet functions in Ref. [36] at NLO.
Concerning the hard corrections, the real emission of a hard gluon is kinematically impossible near threshold. In contrast, the hard virtual contribution does not vanish. We do not attempt to derive or use analytical results for the hard contribution. Instead, we obtain this contribution numerically from the NLO MadGraph framework as detailed in the following. To access the virtual corrections at threshold, a MadLoop library for the process at hand is most suitable. To obtain the one-loop hard function for given p T and M 2 , at leading power we only have to consider phase-space points at threshold. As for the LO cross section, we integrate the MS-renormalized virtual part of the NLO result over the dark-matter phase-space and divide the result by the respective LO result. Alternatively, one can use the process with a mediator of mass M 2 in the final state. At threshold, the phase-space is completely fixed up to an azimuthal symmetry and no integration is needed. One calculates the interference contribution of the MS-renormalized one-loop and the tree-level amplitude and normalizes by the corresponding squared tree-level amplitude.
The hard function for given p T and M 2 does neither depend on the dark-matter mass nor on the mediator mass. It depends only on the variables Q 2 andŝ T representing the hard scale µ h , the scale µ and the renormalization scale µ R since the renormalization of the strong coupling is part of the virtual corrections. Hence, the hard function includes logarithms of the scale ratio µ 2 /µ 2 h . One can write the hard function as a scale independent piece given numerically and the scale-dependent terms which are available analytically. At the level of the complete matching coefficient, this is discussed at the end of the chapter, where we choose the hard scale µ 2 h = Q 2 and the ratio r = p T /m T as independent variables instead of Q 2 andŝ T . Since the hard function is a purely virtual correction, it is proportional to δ(1 − z) and, therefore, it has no N -dependence in Mellin space.
In Mellin space, the threshold limit corresponds to large N . Hence, in the threshold limit, large logarithms can only be avoided in one of the three contributions (soft, collinear or hard) by a corresponding scale choice but not in all of them. The NLO result therefore necessarily contains large logarithms, no matter how the scale µ is chosen. Since the factorization scale, where the PDFs are evaluated, is usually taken as equal to the hard scale, we choose to use µ ∼ µ h . To control the large logarithms in the soft and the jet function, we use renormalization-group techniques as introduced in the following.
It can be shown that the cross section at leading power factorizes to all orders in perturbation theory [24,36] and the Mellin-space cross section is given by an ordinary product of the LO cross section and the soft, the jet and the hard function. Using SCET convention [36] one writes dσ a (N ) = S a (N + 1) J a (N + 1) H a dσ LO a (N ).
Note that in this approach the physics of the incoming jets is completely contained in the soft function and only the hard collinear emissions of the outgoing jet give rise to the jet function. Hence, the assignment of the different physics contributions to the factorized functions differs compared to the direct QCD approach, cf. [23]. The factorized structure implies that the different pieces obey renormalization-group equations. For the soft and the jet function, one has where the anomalous dimensions γ = i (α S (µ 2 )/π) i γ (i) are given as a power series in the strong coupling. The O(α s ) coefficients can be directly read off the NLO expressions explicitly given before, i.e.
for the soft functions and for the jet functions. At NLL, the only necessary two-loop pieces are the O(α 2 s ) coefficients of γ S A a and γ J A q,g . They are related to the process-independent cusp-anomalous dimension according to where γ (1) cusp = 1 and γ (2) The soft function in Mellin space is perturbatively well-behaved for scales close to its natural scale, i.e. µ 2 s ∼ŝ T /N 2 , because no large logarithms appear. Given that Q 2 ∼ŝ T , we choose Q 2 /N 2 as the starting scale for the renormalization-group evolution. With this choice, our results for the exponentials derived below agree with Ref. [23]. The solution of the renormalization-group equation reads Separating and expanding the exponential containing onlyN -independent terms to the desired order in α s , one can cast the resummed soft function into the form where the exponential function for the resummed large logarithms reads and one has S a (µ 2 ) = S NLO a (µ 2 /N 2 ) at O(α s ). In analogy, the jet function in Mellin space is perturbatively well-behaved for µ 2 c ∼ Q 2 /N . We choose Q 2 /N as the starting scale for the renormalization-group evolution. The resummed jet function is then given by where where one needs γ S A a , γ J A a up to O(α 2 s ) and γ S B a , γ J B a up to O(α s ). The matching coefficient can be truncated at O(α s ). The O(α s ) contributions to C a (µ 2 ) due to S a (µ 2 ) and J a (µ 2 ) are known analytically as discussed before, while the hard function is obtained numerically with the help of MadLoop. The Mellin transform of the leading order cross section d 2σLO a (N )/dp 2 T dM 2 is calculated analytically and can be expressed in terms of hypergeometric functions, see Appendix A.
The resummed exponentials in Mellin space agree with the results in Ref. [23], where one can also find analytic expressions for the integrals in the exponential. As in Ref. [23], for our numerical results, we resum large logarithms of N instead ofN , i.e. N -independent exponentials containing γ E are expanded to O(α s ) and lead to additional terms proportional to γ E and γ 2 E in S a (µ 2 ) and J a (µ 2 ). The scale µ is physically relevant because we define the PDFs at µ F . Hence, we have to identify µ = µ F . We have checked that the µ F -dependence of the matching coefficient is given by the known logarithms introduced via the PDF renormalization. In addition, we have to define the numerical input for the strong coupling constant at some scale µ R . Within the matching coefficient, all terms of O(α s ) are simply evaluated using α S (µ R ). Hence the matching coefficient depends also on the choice of µ R . The renormalization of the strong coupling in the LO cross section is part of the hard function. Also R exp depends on the input value α S (µ R ) from which α S at another scale is derived using its renormalization-group running (see e.g. Ref. [23]). The matching coefficient can be written as where P g = 3 2 C F and P q = Pq = 3 4 C F + πb 0 . The constants P a are related to PDF renormalization, where the C F -pieces originate from the renormalization of the quark PDF and the piece proportional to b 0 from the renormalization of the gluon PDF. The terms proportional to γ E arise because we only keep log N instead of logN in the exponentials, as mentioned above. The functionsĈ a contain all scale-independent terms and include the hard contributions which are only available numerically. The scale independent partsĈ a of the matching coefficients only depend on the ratio r = p T /m T and are presented in Figure 1.
To obtain numerical results from the resummed expression, we have to perform the inverse Mellin transform of the Mellin-space cross section, i.e.
where y T = (p T + m T )/ √ s is the hadronic threshold variable and f i (N ) are the Mellin moments of the PDFs. Following the minimal prescription [39], the constant c is chosen to the left of the poles due to the running of the strong coupling α s in the exponents of dσ NLL a (N ), but to the right of all other poles of the integrand. To use the standard PDF implementations, which are not available in Mellin space, we rewrite the resummed hadronic cross section according to [40] dσ NLL Here, dσ NLL a (ŝ) is the inverse Mellin transform of dσ NLL a (N )/N 2 with respect to the partonic threshold variable. Moreover, thef i are the inverse Mellin transforms of N f i (N + 1), which are given by derivatives of the standard PDF functions according tõ , and the sums run over all initial-state partons i, j contributing to the LO process with parton a in the final state. Introducing the convergence factor 1/N 2 before performing the inverse Mellin transform of the resummed partonic cross section makes the numerical integration feasible.

Numerical Results
In this section, we discuss numerical results for benchmark scenarios within the simplified model discussed in Section 2. We compute the full NLO predictions using Mad-Graph5_aMC@NLO [18] using a FeynRules [19,20] implementation of the investigated dark-matter model which is publicly available [13,41]. The NLL' prediction which has been discussed in detail in Section 4 is matched to this NLO result. To avoid double counting, we add the NLO and the NLL' prediction and subtract the expansion of the resummed result to next-to-leading order in the strong coupling α s . Hence, only the contributions of the resummed cross section beyond NLO are taken into account.
We use the PDF set PDF4LHC15_nlo_100 as provided by LHAPDF [42] throughout our analysis, i.e. for LO, NLO and NLO+NLL' predictions. Our central scale choice is given by where the sum runs over all final state jets. This is a common scale choice in fixedorder calculations for processes involving jet production, cf. [13]. Moreover, H T /2 is a scale choice close to the two hard threshold variables √ŝ T and Q, since at threshold √ŝ T = m T +p T = H T and Q 2 = r 1+r H 2 T ≤ H 2 T /2. We use α s (µ R ) as provided by LHAPDF for the given PDF set. Moreover, we use N f = 5 as the number of light flavours. We present predictions for the LHC at the center-of-mass energy √ s = 13 TeV. We use a fixed width for the mediator which is calculated at LO as discussed in Appendix A.
The theoretical uncertainty of the NLO and the matched NLO+NLL' results due to missing higher-order corrections is estimated in the usual way by independently varying the factorization scale µ F and the renormalization scale µ R . To be precise, we use a standard 7-point variation, i.e. we make predictions for and show the minimal and maximal prediction as an error band around the nominal prediction (ξ F , ξ R ) = (1, 1). As discussed in Section 2, the vector and axial-vector couplings of the mediator to the dark-matter particles and the standard-model quarks are set to respectively. As can be seen in the analytic LO results presented in Appendix A, the couplings to the SM quarks g V,A SM enter the cross section explicitly as the sum of the squared vector and axial-vector couplings and in the mediator width. Hence, changing g V,A SM mainly rescales the cross section by an overall factor which depends on whether the mediator is produced on-shell or off-shell. To a good approximation, the cross section also only depends on the squared sum of g V,A χ if the mediator mass is much larger than the dark-matter mass and the production is dominated by mediators close to their mass shell. For m Y < m χ , the production involves off-shell mediators and it is velocity suppressed near M 2 = 4m 2 χ for axial-vector couplings. Hence, for axial-vector couplings the differential cross-section with respect to M 2 peaks at higher M 2 , the cross section will be smaller compared to a model with a vector coupling of the same size but no qualitative differences are expected for the radiative corrections.
The observable of interest is the missing transverse-momentum distribution. To obtain the p T,miss spectrum, the double-differential resummed result discussed in Section 4 is numerically integrated with respect to the invariant mass M 2 of the dark-matter system.
At NLO precision, the model has been already investigated in Ref. [13]. Using the setup of [13], we have reproduced the NLO cross section at the level of 0.5% or better. For the LO results and the corresponding K-factors we have observed substantial differences. It turned out that the LO cross sections in Ref. [13] have not been obtained with the PDF set stated in the paper.
We discuss three benchmark points in detail. We have chosen two parameter sets for the mediator mass m Y and the dark-matter mass m χ which lead to a large invariant mass of the dark-matter pair. In this case, the resummed cross section can be expected to reliably estimate the higher-order corrections beyond NLO accuracy. The first benchmark point is (m Y , m χ ) = (100, 500) GeV, i.e. a relatively light mediator with a heavier dark-matter particle leading to off-shell production. The second benchmark point is (m Y , m χ ) = (1000, 50) GeV, i.e. a heavy mediator with a relatively light dark-matter particle, which can be produced from an on-shell mediator. For these parameter points the resummed prediction expanded to NLO is slightly larger than the full NLO result and reproduces it within 5% throughout the considered range for p T,miss . In addition, we investigate the model with the mass combination (m Y , m χ ) = (100, 1) GeV, i.e. a relatively light mediator with almost massless dark-matter particles. This mass combination is kinematically very close to the main SM-background process, i.e. Z-boson production in association with a high-p T jet with subsequent invisible Z-boson decays into neutrinos. Since the final-state dark-matter system is rather light in this case, the difference between the NLO-expanded resummed result and the full NLO prediction is approximately -10% The p T,miss distributions are shown in Figure 2 and Figure 3 in the range between 150 GeV and 2 TeV. The NLO K-factor for the heavy dark-matter systems in Figure 2 does not vary much as a function of p T,miss and amounts to 30-40%. For the lighter dark-matter system in Figure 3 the NLO K-factor varies from 40 to 50%. For mediator and dark-matter masses below 100 GeV, the K-factor grows and the accuracy of the expanded NLL' result is getting worse. For example, one finds K ∼ 1.7 for (m Y , m χ ) = (10, 1) GeV for a p T,miss cut of 500 GeV. For smaller masses, di-jet final-states start to sizeably contribute, where the mediator is emitted from one of the jets. These contributions cannot be expected to be well approximated by the threshold behavior of the mono-jet final-state. Hence, in the following, we do not consider final states, for which the invariant mass of the dark-matter system is small. For the benchmark points with a massive dark-matter final state shown in Figure 2, the resummation adds roughly 5% of the LO predictions on top of the NLO prediction in the region where the LHC is most sensitive. In the very tails of the distributions, the corrections due to resummation increases up to 12% . Hence, as expected, the largest impact is observed in the tail of the p T,miss -distribution at 2 TeV, where the dark-matter production takes place closer to threshold. For (m Y , m χ ) = (100, 1) GeV, the correction is almost negligible at small p T,miss but also rises to almost 10% of the LO result in the tail of the distribution at 2 TeV. Note that it is essential to take into account the O(α s ) terms of the matching coefficient in (32) on top of the resummed exponentials. Using NLL instead of NLL' results, the full NLO result is not approximated well by the expanded resummed result and the correction beyond NLO is much smaller.
The matched NLO+NLL' prediction also considerably reduces the scale uncertainty at large transverse momenta. At large p T,miss above 1 TeV the NLO error band is reduced by roughly 50% for all three benchmark scenarios. For smaller values of p T,miss the reduction of the size of the error band decreases with decreasing p T,miss , as shown in Figure 3. As expected, these results reproduce known results for vector-boson plus jet production in the Standard Model [37]. Note that the reduction of the scale error is less pronounced if we define the scale-variation error by a nine-point variation, i.e. also including (ξ F , ξ R ) = (0.5, 2), (2, 0.5). The NLO error band is not affected by the different definition, while the width of the matched NLO+NLL' scale-variation band is increased at the upper end. Using the nine-point variation, there is no relevant improvement in the scale uncertainty for small transverse momenta. For (m Y , m χ ) = (100, 1) GeV, the width of the matched NLO+NLL' error band is even bigger than the NLO error band. However, the error band of the fixed-order prediction grows to more than 15% at transverse momenta close to 2 TeV, while the width of the error band for the matched predictions stays more or less constant at 10%, using the nine-point scale variation.

Conclusions
We have computed the resummed cross section for the hadroproduction of a dark-matter fermion pair in a simplified dark-matter model to NLO+NLL' accuracy. We have derived the corresponding soft and jet functions analytically using the strategy of regions. The hard function has been computed numerically using MadLoop. The resummation has been performed in Mellin-space, and we have matched the resummed results to the full fixed-order NLO predictions obtained with MadGraph5_aMC@NLO.
We have presented numerical results for dark-matter production at the LHC for various mediator and dark-matter masses. For dark-matter pairs produced with a large invariant mass in the TeV range, resummation increases the NLO prediction by approximately 5% to 10%. For smaller invariant masses of the dark-matter pair down to 100 GeV, the increase of the cross section beyond the NLO prediction is slightly reduced. For even smaller invariant masses, the threshold resummation as a means to estimate the leading higher-order corrections becomes questionable. Note that it is essential to include the NLO matching coefficient within the resummed result. We also observe a significant reduction in the scale uncertainty of the predictions, where the size of the NLO error band is reduced by 50% at large transverse momentum.
The combination of simple analytical formulas with numerical results from the Mad-Graph framework allows one to obtain the NLO+NLL' predictions for dark-matter production at the LHC in a straightforward manner.
where the theta-functions determine which decay channels are kinematically accessible for a given mediator mass, and Γ SM Y is obtained from Γ χ Y by replacing the couplings to the dark-matter particle by the SM couplings to the quarks. To be specific, we use m t =172 GeV, and the width of the mediator for the benchmark points is given by Γ Y = (0.056, 0.025, 0.051) m Y for (m Y , m χ ) = (1000, 50), (100, 500), (100, 1) GeV, respectively.
For the resummed cross section, we need the Mellin-space representation of the cross section according to d 2σLO a (N ) dp 2 it can be expressed in terms of hypergeometric functions, where z = (r−1) 2 (r+1) 2 .