Search for dark matter in events with a leptoquark and missing transverse momentum in proton-proton collisions at 13 TeV

A search is presented for dark matter in proton-proton collisions at a center-of-mass energy of $\sqrt{s}=$ 13 TeV using events with at least one high transverse momentum ($p_\mathrm{T}$) muon, at least one high-$p_\mathrm{T}$ jet, and large missing transverse momentum. The data were collected with the CMS detector at the CERN LHC in 2016 and 2017, and correspond to an integrated luminosity of 77.4 fb$^{-1}$. In the examined scenario, a pair of scalar leptoquarks is assumed to be produced. One leptoquark decays to a muon and a jet while the other decays to dark matter and low-$p_\mathrm{T}$ standard model particles. The signature for signal events would be significant missing transverse momentum from the dark matter in conjunction with a peak at the leptoquark mass in the invariant mass distribution of the highest $p_\mathrm{T}$ muon and jet. The data are observed to be consistent with the background predicted by the standard model. For a first benchmark scenario, dark matter masses up to 430 GeV are excluded for leptoquark masses $M_\mathrm{LQ}\approx$ 1000 GeV, and up to 300 GeV for $M_\mathrm{LQ}\approx$ 1170 GeV. For a second benchmark scenario, dark matter masses up to around 480 GeV are excluded for $M_\mathrm{LQ}\approx$ 1100 GeV.


Introduction
Dark matter (DM) has been a subject of intense interest for decades. Extensive astrophysical evidence for DM exists [1][2][3], such as from observations of the dynamics of galaxy clusters and measurements of anisotropies in the cosmic microwave background. Nonetheless, the nature of DM remains unknown and it has not been observed outside the astrophysical context. Its relic density is determined to be Ω DM = (0.1186 ± 0.0020)/h 2 [4,5], where h is the Hubble constant.
Dark matter could potentially be created in high-energy proton-proton (pp) collisions, such as at the CERN LHC. Because of its presumed weakly interacting nature, a DM particle produced at the LHC would escape unobserved, manifesting itself as missing transverse momentum p miss T in the reconstructed events. The most generic signal for DM at the LHC thus consists of an excess, relative to the standard model (SM) expectation, of events with sizable p miss T recoiling against a visible SM object such as a jet, a photon, or an electroweak boson. Such searches have been conducted at the LHC by the ATLAS and CMS Collaborations [6-11] but with no evidence, to date, for DM [12]. The absence of a signal in these generic searches suggests that alternative strategies should be pursued.
In this Letter, we present a search for DM at the LHC using a new approach, based on the coannihilation paradigm introduced in Ref. [13]. The data, corresponding to an integrated luminosity of 77.4 fb −1 of pp collisions, were collected by the CMS Collaboration in 2016 (35.9 fb −1 ) and 2017 (41.5 fb −1 ) at a center-of-mass energy of √ s = 13 TeV. The coannihilation process considered arises within a general class of simplified models in which a DM particle is either annihilated or produced in conjunction with a so-called coannihilation partner, denoted "X". At the LHC, DM could thus be produced through a reaction like M → DM + X, where M is a mediator representing either an SM or a beyond-the-SM particle. To be consistent with the observed value of Ω DM , the fractional mass difference ∆ X,DM ≡ (m X − m DM )/m DM between the X and DM particles should be less than ≈0.2 [13].
The considered coannihilation paradigm introduces many DM signatures that are not covered by current searches. Here, we consider the principal case-study scenario of Ref. [13], in which the mediator M is a scalar leptoquark (LQ) and the particle X is a new Dirac fermion. An LQ is a hypothetical color-triplet, fractionally charged boson that carries both lepton and baryon quantum numbers. Leptoquarks appear in many extensions of the SM, such as grand unification theories [14-17] and models with composite quarks and leptons [18]. To be consistent with experimental constraints on flavor changing neutral currents, we assume the LQ to couple to a single SM flavor generation only [19], taken to be the second generation for this study. We choose the second generation because muons provide a clear experimental signature. We further assume pair production of LQs, as is predominantly expected in pp collisions. In the examined scenario, one LQ decays to a muon and a charm (c) quark while the other decays, through the coannihilation paradigm, to DM and X. We assume that LQs can decay only through these two channels. The X particle subsequently decays through a crossed coannihilation process to a DM particle and an off-shell LQ, where the decay products of this latter particle, also a muon and a c quark, have low transverse momentum (p T ) because of the smallness of ∆ X,DM and are potentially undetected.
An example Feynman diagram is shown in Fig. 1. The final state consists of a high-p T muon and a high-p T jet from the decay of the on-shell LQ, p miss T from the DM particles, and low-p T SM objects from the decay of the off-shell LQ. Note that, in the analysis, we do not employ c quark tagging criteria but rather-to improve the signal event selection efficiency-utilize generic untagged jets as the c quark jet candidates. The existence of the signal process is inferred by a Figure 1: An example Feynman diagram for the signal process considered in this study, where g is a gluon, LQ a leptoquark, DM a dark matter particle, and X a new Dirac fermion. The superscript "*" indicates an off-shell particle.
peak at the LQ mass m LQ , in the invariant mass m µj distribution of the high-p T muon and jet, in conjunction with significant p miss T from the DM. This peak at the LQ mass provides a unique signature in the search for signal processes containing DM. In contrast, generic searches for DM, in which there are no new particles other than DM and intermediate mediator states, mostly rely on a mere enhancement in the tail of the p miss T distribution.
The most restrictive lower limit on the mass of a pair-produced second-generation LQ, assuming an LQ branching fraction B(LQ → cµ) = 100%, is currently 1530 GeV [20]. However, when the decay to DM and X is allowed, the LQ → cµ branching fraction is reduced and the limit on m LQ becomes weaker. The branching fraction then also depends on m DM , ∆ X,DM , and B 0 , where B 0 is defined as the B(LQ → cµ) in the limit of massless X and DM particles and is related to the LQ coupling strengths y q (to the lepton and quark) and y D (to the DM and X particles) with the following formula: Following Ref.
[13], we set y D = 0.1 and ∆ X,DM = 0.1. We consider two values for B 0 : 0.5 and 0.1. Recasting the results of Ref.
[20] in terms of the scenario of Fig. 1, and taking m DM = 300 GeV as a representative value, the lower limit on a second-generation LQ is reduced to 1340 GeV for B 0 = 0.5 and to 960 GeV for B 0 = 0.1.
Again following Ref.
[13], the DM particle is assumed to be a Majorana fermion with the gauge group structure (1, 1, 0), where the numbers in parentheses indicate the color SU(3) C , weak isospin SU(2) L , and weak hypercharge U(1) Y multiplet dimensions, respectively. We use the convention Q = T 3 + Y/2 for the electric charge Q of the particle, with T 3 the third component of weak isospin and Y the weak hypercharge. The corresponding assignments for both the X and LQ particles are (3, 2, 7/3). For the LQ, T 3 = 1/2 and Q = 5/3.
The principal SM backgrounds in this search arise from events with a W boson and jets (W+jets) or with a top quark-antiquark (tt) pair: in both cases, the leptonic decay of a W boson can yield a high-p T muon and neutrino, where the neutrino can lead to significant p miss T . Events with single top quark or diboson (WW, WZ, and ZZ) production similarly can enter the background, al-though at a lower level. Other smaller sources of SM background arise from quantum chromodynamics (QCD) events, namely events with a multijet final state produced exclusively through the strong interaction, and from events with a Z boson and jets (Z+jets). A QCD event can enter the background if a muon and a neutrino are produced through the semileptonic decay of a quark, or if a jet is erroneously identified as a muon in conjunction with spurious p miss T arising from the mismeasurement of jet p T . Events with Z+jets production can enter the background if one of the leptons in Z → µ + µ − decays is not reconstructed or lies outside the acceptance of the detector, leading to p miss T , or if p miss T arises because of misreconstructed jet p T .

The CMS detector and trigger
The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip inner tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Extensive forward calorimetry complements the coverage provided by the barrel and endcap detectors. Muons are detected in gas-ionization chambers embedded in the steel flux-return yoke outside the solenoid. A detailed description of the CMS detector, together with a definition of the coordinate system and the relevant kinematic variables, is given in Ref. [21].
Events of interest are selected using a two-tiered trigger system [22]. The first level, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz. The second level, based on an array of microprocessors running a version of the full event reconstruction software optimized for fast processing, reduces the event rate to around 1 kHz before data storage. The set of triggers used for this analysis requires events to contain a muon with p T > 50 GeV.

Event reconstruction and selection
Individual particles are reconstructed with the CMS particle-flow (PF) algorithm [23], which identifies them as charged hadrons, neutral hadrons, muons, electrons, or photons. Muon reconstruction is performed by matching a track segment reconstructed in the inner tracker with a track segment reconstructed in the muon detector and performing a global fit of the hits from the two track segments. The candidate muons are required to satisfy the tight selection criteria of Ref.
[24], to pass within 2 mm of the primary event vertex in the direction along the beam axis and within 0.45 mm in the plane perpendicular to that axis, and to have p T > 60 GeV and |η| < 2.4, where η is the pseudorapidity. Electrons are reconstructed by matching energy deposits in the electromagnetic calorimeter with track segments in the inner tracker [25] and are required to have p T > 15 GeV and |η| < 2.5. Hadronically decaying τ leptons (τ h ) are reconstructed using the hadrons-plus-strips algorithm described in Ref. [26] and are required to have p T > 20 GeV and |η| < 2.3. The electron and τ h candidates are used to veto events, as described below, and are selected using loose [25, 26] identification criteria.
The primary event vertex is defined to be the reconstructed interaction vertex with the largest value of summed physics-object p 2 T . The physics objects considered for this purpose are the jets found by clustering the charged-particle tracks assigned to the vertex, using the anti-k T jet finding algorithm [27,28] with a distance parameter 0.4, and the associated missing transverse momentum, taken as the negative of the vector p T sum of those jets. The transverse momentum imbalance p miss T in an event is calculated as the negative of the vector p T sum of all PF candidates. Its magnitude is denoted p miss T . Events are required to have p miss T > 100 GeV.
To suppress the contributions of muons that arise from hadron decays, muon candidates are subjected to an isolation requirement. The scalar p T sum of charged hadron, neutral hadron, and photon PF candidates within a cone of radius ∆R ≡ √ (∆η) 2 + (∆φ) 2 = 0.4 around the muon direction, where φ is the azimuthal angle, is calculated. The expected contributions of neutral particles from additional pp interactions in the same or neighboring bunch crossings (pileup) are subtracted [29]. An isolation variable I is defined by dividing this sum by the muon p T . The isolation requirement is I < 0.15. At least one isolated muon candidate is required to be present in the event.
The reconstruction of jets is performed by clustering PF candidates using the anti-k T algorithm with a distance parameter of 0.4, excluding charged-particle tracks not associated with the primary vertex. The jet energies are corrected for the combined response function of the calorimeters [30] and to account for the expected contributions of neutral particles from pileup [29,31]. Jets are required to appear within |η| < 2.4. Bottom (b) quark jets are identified (b tagged) from this sample using the combined secondary vertex (CSVv2) algorithm at the tight working point [32], which yields a b quark jet identification efficiency of approximately 60% and a misidentification probability for gluon and light-flavored quark jets of about 0.5%. Jets tagged as b jets are required to have p T > 30 GeV.
The leading (highest p T ) jet in an event is required to have p T > 100 GeV and to be separated by ∆R > 0.5 from the leading muon candidate. The leading muon and leading jet are then combined to form the LQ candidate. Studies with simulated signal events establish that, for the values of model parameters used in the present study, this matching identifies the correct combination over 98% of the time.
To suppress background from tt production, events are rejected if they contain a b-tagged jet, an electron candidate, or a τ h candidate. The veto on events with a b-tagged jet reduces the background from events with a top quark by more than a factor of 2 while reducing the signal efficiency by only around 10%. The vetoes on electron and τ h candidates also suppress the W+jets and Z+jets background. The W+jets background is further suppressed by requiring the transverse mass m T [33] formed from the p T vector of the leading muon and p miss T to exceed 100 GeV. The Z+jets background is further reduced by eliminating events with a loosely identified and isolated (I < 0.25) muon candidate if that muon candidate has p T > 10 GeV, an opposite charge to the leading muon, and forms a dimuon mass with the leading muon within 10 GeV of the Z boson mass.
Background from QCD events mostly arises when the p T of one of the highest p T jets is underestimated or when a hadron in a jet undergoes a semileptonic decay, introducing p miss T that is aligned with that jet. To suppress this background, the angular difference ∆φ between the leading jet and p miss T , and between the leading muon and p miss T , is required to exceed 0.5.
The above requirements are referred to as the "preselection" criteria, and form the basis for the definition of several control regions used to evaluate background, as described in Section 5. The final selection criteria, corresponding to the signal region, are the same as the preselection criteria except for a more stringent requirement on m T , m T > 500 GeV. This condition is determined from optimization studies utilizing simulated signal and SM event samples. For m LQ > 800 GeV, the signal efficiency for events satisfying the preselection criteria is around 73%, essentially independent of m LQ . For the signal region criteria, the signal efficiency varies from 47 to 63% as m LQ increases from 800 to 1500 GeV.

Event simulation
Monte Carlo (MC) simulation of signal and background processes is used to validate the analysis procedures, evaluate background, and determine the signal efficiency. Simulation of tt and single top quark events is performed at next-to-leading-order (NLO) accuracy with the POWHEG 2.0 [34][35][36][37][38][39] event generator. To describe W+jets and Z+jets production, the MAD-GRAPH5 aMC@NLO 2.2.2 [40,41] program at leading order (LO) is used. A K factor, calculated as described in Ref. [42], is applied as a function of boson p T to account for next-to-NLO (NNLO) corrections. The statistical precision of our available MADGRAPH5 aMC@NLO W+jets sample is low for off-shell W boson masses above around 100 GeV. Therefore, to describe W+jets production for W boson masses above 100 GeV, we use the LO PYTHIA 8.212 [43] program with the CUETP8M1 tune [44], rather than MADGRAPH5 aMC@NLO, with a K factor to account for NNLO corrections applied as a function of the W boson mass. This K factor is determined using the FEWZ  For simulated samples at LO (NLO), the NNPDF3.0LO (NNPDF3.0NLO) [50] parton distribution functions (PDFs) are used. All samples are interfaced to PYTHIA 8.212 to describe parton showering and hadronization. Pileup interactions are modeled using simulated minimum bias event samples generated with PYTHIA, with the distribution of pp interactions per bunch crossing adjusted to reproduce the observed spectrum.
The response of the CMS detector is modeled, for both the signal and background samples, using the GEANT4 [51] suite of programs. Small differences between the data and simulation in the trigger, particle identification, and muon isolation efficiencies, and in the jet p T resolution and p miss T , are accounted for through the application of scale factor corrections.

Background estimation
Background from tt production is evaluated from simulation, with the events reweighted to reproduce the observed distribution of the top quark p T [52, 53]. The normalization is determined using a data-to-simulation scale factor derived from a tt-enhanced control sample. The control sample is defined using the preselection criteria of Section 3 except that events are required to contain at least one b-tagged jet. The purity of tt events in the control sample is estimated to be 85%. The remaining 15% of the sample is composed primarily of events with W+jets and single top quark production. The scale factor is defined as the ratio of the number of events in the control sample from data, after subtraction of the non-tt components, estimated from simulation, to the number in the corresponding simulated tt control sample scaled to the same integrated luminosity. The scale factor is found to be 0.95 ± 0.01 (stat) for the 2016 data and 1.16 ± 0.01 (stat) for the 2017 data and is essentially independent of the m T selection requirement. The difference in the scale factors between 2016 and 2017 arises from changes in the running conditions, the reconstruction procedures, and the tuning of the simulation programs. Figure 2 (left) shows the distribution of m µj in the combined 2016+2017 tt control sample for data and simulation.
The systematic uncertainty in the scale factor is derived using an orthogonal tt-enhanced control sample, selected with the preselection criteria except that events must contain at least one The background from W+jets production is similarly estimated from simulation, with a correction to the normalization obtained from a W+jets-dominated control sample. The control sample is defined using the preselection criteria except with the additional requirement 50 < m T < 150 GeV. The purity of W+jets events in this sample is about 80%. The main remaining contribution is from tt production. A data-to-simulation normalization scale factor is obtained by subtracting the non-W+jets contribution, estimated from simulation, from the control sample in data, and dividing the resulting number of events by the number of events in the simulated control region normalized to the same integrated luminosity. The scale factor is found to be 1.02 ± 0.01 (stat) and 1.11 ± 0.01 (stat) for the 2016 and 2017 data, respectively. The distribution of m µj in the combined 2016+2017 W+jets control sample, for data and simulation, is shown in Fig. 2 (right).
To evaluate a systematic uncertainty in the W+jets background prediction, we examine the level of agreement between the data and simulation for the normalization of the m µj distribution in a Z(→ µ + µ − )+jets control sample for which one of the two muons in each event is removed to emulate a sample of W+jets events. This control sample is selected by requiring two oppositely charged isolated muons with p T > 30 GeV, with a dimuon invariant mass between 80 and 100 GeV, but with otherwise similar selection criteria to those of the preselection. The Z boson candidate is boosted to its rest frame, the dimuon mass is scaled to the mass of the W boson, and the system is then boosted back to the original laboratory frame. One of the two muons is randomly removed to simulate a neutrino from W boson decay. The resulting missing momentum is added to the p miss T of the event and the value of m T recalculated before applying the signal region selection criteria. On the basis of the level of agreement in the resulting distribution of m µj between data and simulation, a 20% uncertainty is assigned to the W+jets background prediction. Additional checks are performed through comparison of the simulated and observed m µj distributions as the requirement on m T is tightened. These checks validate the assignment of a 20% uncertainty, also for large values of the off-shell W boson mass.
The background from QCD processes is expected to be small. However, since the QCD background primarily arises as a consequence of jet mismeasurement, it is not well modeled by simulation. Thus, we evaluate the QCD background using a method based primarily on data. A control sample is defined using the signal region criteria of Section 3 except that muon candidates are required to fail the isolation condition. To estimate the QCD background in the signal region, the events in the control sample are weighted as a function of muon p T by a muon misidentification probability, called the jet-to-muon misidentification rate, determined in a QCD-enriched event sample denoted the "low-∆φ" sample. The low-∆φ sample is defined in the same manner as the signal region except the angle ∆φ between p miss T and the leading jet is required to be ∆φ < 0.5 rather than ∆φ > 0.5 and we require m T > 100 GeV rather than m T > 500 GeV. The jet-to-muon misidentification rate is defined, from this sample, as the ratio of the number of events that satisfy the muon relative isolation criterion I < 0.15 to the number of events with no requirement on the isolation, after subtracting the non-QCD components, evaluated with simulation, from both the numerator and denominator. The purity of QCD events in the numerator is about 25%, while that in the denominator is approximately 70%. The jet-to-muon misidentification rate is parameterized in terms of the muon p T using an analytical function and varies from 5% for muon p T = 60 GeV to 50% for p T > 300 GeV.
A 50% uncertainty is assigned to the QCD background prediction to account for the uncertainty in the jet-to-muon misidentification rate. This uncertainty primarily arises from the uncertainties in the normalization of the non-QCD processes subtracted from the numerator and denominator in calculating the misidentification rate. This uncertainty is sufficiently large to account for uncertainties from non-QCD background processes. Other sources of uncertainty, such as the choice of the analytic function used to parameterize the jet-to-muon misidentification rate, or the uncertainties in the values of the fit parameters, are negligible in comparison.
The backgrounds from events with diboson, single top quark, and Z+jets production are estimated from simulation. An uncertainty of 15% is assigned to the diboson background prediction based on the level of agreement between data and simulation in a diboson-enhanced control sample selected by requiring events to contain three leptons (e or µ), two of which are consistent with arising from Z boson decay. Uncertainties of 15 and 10% are assigned to the single top quark and Z+jets backgrounds, respectively, based on the results of Refs.

Systematic uncertainties
We evaluate systematic uncertainties that affect the normalization or shape of the m µj spectrum, either in the signal or background predictions. Uncertainties specific to individual background components were presented in Section 5.
The uncertainty in the trigger efficiency is 5%. Those in the muon reconstruction and isolation efficiency [56], the electron reconstruction efficiency [56], and the τ h reconstruction efficiency [57, 58] are 5, 2, and 2%, respectively. The uncertainty related to the b-tagging misiden-tification is 1% [32]. The uncertainty in the pileup description in simulation is assessed by varying the total inelastic cross section by 4.6% [59], and is found to be 1%. Statistical uncertainties related to the limited number of events in the data control samples are accounted for as described in Ref.
[60], while those related to the limited number of events in simulation are accounted for by allowing the content in each bin of the simulated distributions to vary within its statistical uncertainty. The uncertainty in the integrated luminosity is 2.5% [61] for the 2016 data and 2.3% [61] for the 2017 data. Uncertainties associated with the jet energy scale, the jet energy resolution, and p miss T are also evaluated [62,63]. These uncertainties affect both the shape and normalization of the simulated signal and background distributions. Uncertainties related to the top quark p T reweighting in simulated tt events are evaluated by varying the reweighting parameters between zero and twice their nominal values [52,53].
Uncertainties related to the PDFs, evaluated for the signal acceptance, are determined following the PDF4LHC prescription [64] and are found to be 3%. Those related to the renormalization and factorization scales, evaluated for the signal yields and for the tt and W+jets backgrounds, are estimated by varying each scale independently, and also coherently, by a factor of 2.0 and 0.5. The largest upward and downward variations in the results are used to define an uncertainty envelope. The uncertainty associated with the method chosen to determine the K factor for the W+jets background evaluation is estimated to be 5%.

Results
The observed distribution of m µj is presented in Fig. 3. The results are shown in comparison to the predictions for the SM background. For purposes of illustration, the predictions of two signal models with m LQ = 1000 GeV and m DM = 400 GeV are also shown: one with B 0 = 0.5 and the other with B 0 = 0.1. Numerical values are given in Table 2.
The data are found to be consistent with the SM predictions within the uncertainties. There is a small excess of events above the SM prediction in the m µj region between 1600 and 1900 GeV, consistent with a statistical fluctuation. This excess corresponds to a statistical significance of around 1.5 standard deviations. Thus, we do not obtain evidence for DM or LQ production.
A binned maximum likelihood fit is applied to the m µj distribution in the signal region. The fitted parameters are the yields of the individual background components listed in Table 2, the signal yield, and various nuisance parameters. The nuisance parameters are introduced to treat systematic uncertainties. Log-normal probability distributions are used for nuisance parame- Obs./Exp. ters that affect the normalizations of the signal and background yields. Gaussian probability distributions are used for nuisance parameters that affect the shape of the m µj distribution. All uncertainties are assumed to be fully correlated across bins except those that are statistical in origin, which are assumed to be uncorrelated. Upper limits at 95% confidence level (CL) are determined on the product of the signal production cross section and branching fraction. These limits are calculated using a modified frequentist approach with the CL s criterion [65,66] and an asymptotic approximation for the test statistic [67,68]. The limits are determined as a function of m LQ and m DM . [20] are seen to be even more restrictive, excluding most of the m DM -m LQ combinations in Fig. 4 (left).
The results of the present analysis are of particular interest in the region of model parameter space corresponding to smaller values of B 0 . The results for B 0 = 0.1 are shown in Fig. 4 (right). In this case, the present analysis is more restrictive than the recast of Ref.
[20] for DM masses in the range from 300 and 430 GeV. In this interval, the signal scenario of Fig. 1 is excluded for m LQ values up to 1140 GeV. More generally, the maximum value of dark matter mass that can be excluded is 480 GeV, corresponding to an LQ mass of ≈1100 GeV.

Summary
A search has been performed for dark matter in events containing a muon, a jet, and significant missing transverse momentum. The study is conducted using proton-proton collision data at √ s = 13 TeV recorded with the CMS detector, corresponding to an integrated luminosity of 77.4 fb −1 . It is assumed that dark matter is produced through the production of a leptoquark pair, with one leptoquark decaying to a muon and a jet, and the other to dark matter and low-p T standard model particles. The analysis is performed by searching for a peak in the leptoquark candidate invariant mass m µj distribution formed from the highest p T muon and jet in an event, with the requirement of significant missing transverse momentum, as is expected from the presence of dark matter. The observation of such a peak in this novel search would provide strong evidence for the existence of both dark matter particles and leptoquarks. The data are observed to agree with the standard model background predictions within the uncertainties. Upper limits on the product of the cross section and branching fraction are obtained at 95% confidence level as a function of the leptoquark and dark matter particle masses. For the first benchmark scenario, dark matter masses up to 430 GeV are excluded for leptoquark masses m LQ ≈ 1000 GeV, and up to 300 GeV for m LQ ≈ 1170 GeV. For the second benchmark scenario, dark matter masses up to 480 GeV are excluded for m LQ ≈ 1100 GeV.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. [ [20] CMS Collaboration, "Search for pair production of second-generation leptoquarks at √ s = 13 TeV", (2018). arXiv:1808.05082. Submitted to Phys. Rev. D.