Monojet signatures from heavy colored particles: future collider sensitivities and theoretical uncertainties

In models with colored particle Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} that can decay into a dark matter candidate X, the relevant collider process pp→QQ¯→XX¯+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp\rightarrow \mathcal {Q}\bar{\mathcal {Q}}\rightarrow X\bar{X}\,+\,$$\end{document}jets gives rise to events with significant transverse momentum imbalance. When the masses of Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} and X are very close, the relevant signature becomes monojet-like, and Large Hadron Collider (LHC) search limits become much less constraining. In this paper, we study the current and anticipated experimental sensitivity to such particles at the High-Luminosity LHC at s=14\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=14$$\end{document} TeV with L=3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {L}=3$$\end{document} ab-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document} of data and the proposed High-Energy LHC at s=27\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=27$$\end{document} TeV with L=15\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {L}=15$$\end{document} ab-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document} of data. We estimate the reach for various Lorentz and QCD color representations of Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document}. Identifying the nature of Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} is very important to understanding the physics behind the monojet signature. Therefore, we also study the dependence of the observables built from the pp→QQ¯+j\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp\rightarrow \mathcal {Q}\bar{\mathcal {Q}} + j $$\end{document} process on Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} itself. Using the state-of-the-art Monte Carlo suites MadGraph5_aMC@NLO+Pythia8 and Sherpa, we find that when these observables are calculated at NLO in QCD with parton shower matching and multijet merging, the residual theoretical uncertainties are comparable to differences observed when varying the quantum numbers of Q\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}$$\end{document} itself. We find, however, that the precision achievable with NNLO calculations, where available, can resolve this dilemma.


Introduction
The nature of dark matter (DM) remains one of the outstanding mysteries in the particle physics today. Among the many possible particle candidates known [1], weakly interacting massive particles (WIMP) are arguably the most theoretically motivated and well-studied scenarios; for a recent review, see [2]. WIMPs with masses in the range of 10 GeV-20 TeV [3] can be stable on the age of the universe, and once they are in thermal equilibrium in the early universe remain so, even after decoupling occurs. Moreover, such stable particles are naturally present in many beyond the Standard Model (BSM) frameworks. For example: In weak-scale supersymmetric (SUSY) theories [4,5], if one assumes R-parity conservation, then lightest SUSY particle is stable, and hence is a potential component of DM. In the universal extra dimension (UED) models [6], the lightest Kaluza-Klein excitations of neutral electroweak bosons can be viable DM candidates.
Notably, a number of these new physics scenarios also involve additional heavy colored particles, Q, that couple to DM candidate(s). Hence, if Q and DM are kinematically accessible at the Large Hadron Collider (LHC), or its potential successors, such as the High Energy (HE)-LHC, then it may be possible to study DM in a laboratory setting. In particular, once produced and if allowed, Q can decay into DM and SM particles leading to a plethora of interesting signatures at the LHC.
At hadron colliders, search strategies for these hypothetical colored particles usually involve investigating jets and leptons produced in association with final-state DM candidates manifesting as large missing transverse energy (E miss T ). In the context of simplified SUSY models, such signatures are now strongly constrained by LHC data if Q has a mass around 1 TeV and below [7][8][9][10]. Such constraints on Q, how-g/q g/q g/q Q Q Fig. 1 Diagrammatic depiction of QQ pair production with an extra hard QCD radiation in pp collisions ever, can be circumvented. One of the most celebrated examples of this is the compressed spectrum scenario [11][12][13]. In this situation, the DM and new colored particles have a small mass splitting. Consequently, the visible decay products in the Q →DM+SM process do not have sufficient momenta to be readily distinguished from SM backgrounds. In other words, the compression of the mass spectrum constrains the visible decay products of Q to possess such low momenta that they fail experimental selection criteria. This leads to significantly smaller selection and acceptance efficiencies, and hence significantly weaker bounds on heavy particle masses.
Even though a compressed mass spectrum represents a special corner of a typical BSM parameter space, the most attractive feature of this situation is that it allows relatively light, colored DM partners in light of present-day LHC data. This is particularly true for stops (t) and gluinos (g) in SUSY [14][15][16][17]. If a compressed scenario is realized in nature, then one can experimentally resolve the soft, i.e., low p T , visible decays of Q by recoiling against a relatively hard, i.e., high p T , electroweak or QCD radiation that, in its own right, is sufficiently energetic to satisfy trigger criteria. One such process, shown diagrammatically in Fig. 1 and the focus of this study, is the inclusive monojet plus E miss T collider signature. 1 Were evidence for the new particle Q established at the LHC, or a successor experiment, it would be crucial to determine Q's properties, especially its mass, spin, and color representation. Generically, such a program would involve investigating various collider observables that can discriminate against possible candidates for Q. For example: cross sections are highly sensitive to the aforementioned spacetime and internal quantum numbers. Consider the cases of a scalar topt (color triplet), a spin-1/2 top partner T p (color triplet), and spin-1/2 gluinog (color octet). For a fixed mass, i.e., mt = m T p = mg, the pair production cross sections for these particles exhibit the hierarchy Conversely, for a fixed cross section, i.e., σ (tt * ) = σ (T pTp ) = σ (gg), one finds that mg > m T p > mt . (2) This implies, however, that were a monojet cross section measured, the result could be replicated by different Q scenarios by a simple tuning of mass m Q . In another way, one cannot constrain the mass of Q from cross section measurements alone without first asserting its color representation and spin. Quantitatively, this is more nuanced due to fact that leading order (LO) calculations are poor approximations for QCD processes, even when using sensible scale choices. Using Ref. [19], one can easily verify that, like the top quark [20,21], QCD corrections at next-to-leading order (NLO) increase the production cross section σ (T pTp ) by O(50)% for TeV-scale T p . This is the case at both √ s = 14 and 27 TeV, and despite scale uncertainties at LO and NLO spanning O (20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30)% and O(10)%, respectively, for σ (T pTp ). Moreover, it is well-known that next-to-next-to-leading order (NNLO) corrections are non-negligible for SM top production [22]. It is also known that such large theoretical uncertainties can greatly limit the interpretation of the experimental results, particularly in searches for so-called top-philic dark matter [23].
The situation, however, is more hopeful following the advent of general-purpose, precision Monte Carlo (MC) event generators. With software suites such as Herwig [24], MadGraph5_aMC@NLO+Pythia8 (MG5_aMC@NLO+PY8) [19,25], and Sherpa [26], automated event generation at NLO in QCD with parton shower (PS) matching is now possible for both SM and BSM [27] processes. Not only can one now readily include potentially important O(α s ) corrections to cross sections normalizations, but parton showers augment fixed order predictions with resummed corrections to at least the leading logarithmic (LL) level. As a consequence, such observables like the associated jet multiplicity in the monojet process, an exclusive observable that is critical to search strategies, is automatically modeled at LO+LL accuracy. This is the lowest order at which the quantity is qualitatively correct. In light of the availability of such sophisticated technology, one is now in position to systematically investigate the impact of QCD corrections on the inclusive monojet process.
In this report, we perform such a dedicated precision study on the inclusive monojet signature in the context of a compressed mass spectrum. As mentioned, observables associated with this process are highly sensitive to the mass, spin, and color representation of the mediating states. Hence, we consider benchmark models with representative mass, spin and color configurations for Q, with Q ∈ {t, T p ,g}. Our study is aimed at the HL-LHC, assuming L = 3 ab −1 at √ s = 14 TeV, and the HE-LHC, assuming L = 15 ab −1 at √ s = 27 TeV. We show how the monojet search strategy can be used to identify the nature of Q from the p T dependence of the leading jet in QQ pair production. We also estimate the precision required to distinguish these new physics scenarios. We quantitatively discuss various sources of theoretical uncertainty, including event generator dependence. Although we do not discuss the second jet distribution intensively in this paper, azimuthal angle correlation of the first and second jet contains information of the spin of Q, as discussed in [28,29]. The remainder of this study continues in the following manner: In Sect. 2 we provide in-dept detail of our computational setup. In Sect. 3, we discuss observed and expected sensitivity of monojet searches at present and hypothetical future facilities, and address various theoretical uncertainties in Sect. 4. A brief outlook on the impact of this work is discussed in Sect. 5, and we conclude in Sect. 6.

Computational and theoretical setup
Systematic studies of QCD radiation in the production of hypothetical, TeV-scale colored particles are now possible due to the availability of precision, general-purpose MC event generators. In practice, this nontrivial task is handled by using several individually published formalisms and software packages that have largely been integrated into a single framework or well-specified tool chain sequence. In this section, we describe our computational and theoretical setups for modeling QQ production with various associated jet multiplicities at LO+PS and NLO+PS in MG5_aMC@NLO+PY8 and Sherpa. For numerical results, readers can go directly to Sect. 3.
The section continues as follows: In Sect. 2.1, we briefly summarize the Q we consider from representative BSM scenarios. In Sect. 2.2, we enumerate the several methods for incorporating additional QCD radiation into QQ production that we employ, and briefly note their main features and formal accuracies. We describe our setup for MG5_aMC@NLO+PY8 and Sherpa, respectively, in Sects. 2.3 and 2.4, and our detector simulation in Sect. 2.5. In Sect. 2.6, we summarize the SM inputs.

Framework for new heavy colored particles
In this analysis, we consider three benchmark BSM candidates for Q: a stop squarkt, a gluinog, and a fermionic top partner T p . To model these states in pp collisions at our desired accuracy, we use the NLO in QCD-accurate Universal FeynRules Object (UFO) [30,31] model libraries available from the FeynRules model database [32]. The O(α s ) counterterms required for NLO computations and contained in these libraries are generated with FeynRules [30], using NLOCT [27] and FeynArts [33]. For illustrative purpose, we choose three mass values for Q, namely m Q = 400, 600, and 800 GeV. We note that as the spacetime and SU(3) c quantum numbers for T p are identical to those of the SM top quark, several publicly available calculations can be adapted in straightforward ways for T p . This includes total cross section predictions for inclusive pp → T pTp production at NNLO in QCD, which we obtain using the HATHOR package [34]. In Sect. 3.3, where we study the dependence of the selected observables on the color and spin representations of Q, we also briefly consider the well-studied [35][36][37] case of a scalar gluon (σ ) at m σ = 600 GeV.
Furthermore, in order to focus on the ISR from QQ production and also to use traditional analysis techniques, we assume that the decay of Q is prompt, with its characteristic lifetime (τ Q ) satisfying d 0 = βcτ Q 100 μm, where β denotes the relative velocity of Q. As a result, Q's total width must respect the boundary Γ Q =h/τ Q 2 meV. In realizations of the Q we consider, e.g., Refs. [38][39][40], this stipulation on Γ Q is readily satisfied by large regions of the models' phenomenologically relevant parameter spaces. We also assume that the Q → X + SM decay is Q's dominant decay mode. By virtue of DM being weakly coupled, this implies that Q's total width scales as Γ Q ∼ g 2 Q m Q , with an effective Q-DM coupling g Q 1. This ensure that the width-to-mass adheres to the inequality Subsequently, we decouple the production of Q from its decay into an invisible particle X using the narrow width approximation (NWA), as implemented in MadSpin [41]. For technical simplicity, we neglect spin correlation between the production and decay of Q. Fort and σ decays, this has no impact due to their scalar nature. For T p andg decays, this can impact the p T and η distributions of X and subleading jets associated with Q's decay, and hence the true E miss T distribution. However, we have checked that this approximation has a relatively small impact (up to 5%) on our specific study since the E miss T observable we use is built from visible objects, and is dominated by contributions from the hard radiation.
More specifically, for each Q, the particle nature of the DM particle X is chosen in accordance with its underlying theory: For T p , X is a neutral scalar and T p are decayed to the two-body final state X +q, where q is a light QCD quark. For t, X is a neutral fermion andt are decayed via an off-shell top quark to the four-body final state,t → t * X → X +bqq . Forg, X is a neutral fermion andg are decayed to the threebody final state X + qq . We enforce a compressed mass spectrum by fixing a small mass gap between Q and X to be Δm = m Q − m X = 20 GeV. Since the mass gap is (relatively) small, the SM decay products of Q are forced to be  [32,40] (relatively) soft. Hence, the SM decay products of Q fail the criteria needed to be identified as the leading jet, further justifying our neglect of spin-correlation in the decays of T p and g. Qualitatively and quantitatively, our results are expected to be robust against varying Δm factors of 2 as this does not change substantially the likelihood of Q decay promptly, nor the likelihood of its decay products being tagged as subleading jets in the event. For pathologically smaller Δm Q , the lifetime of Q can be extended beyond 100 μm, but does not change the prompt monojet collider signature, nor our conclusions. For further discussion on displaced (monojet) signatures, see [42,43]. For hierarchically larger Δm Q , the SM decay products of Q become more energetic and the process transitions to a multi-jet+E miss T collider signature. The signature associated with the latter scenario is outside the scope of our study and hence will not be discussed further.
The choices of Q, their relevant quantum numbers and decay path, and the corresponding UFO library references are summarized in Table 1.

Multi-leg matching and merging prescriptions
The collider signature considered in this work is characterized by the presence of a highp T jet recoiling against the QQ system. It is thus necessary to include at least one QCD radiation at the matrix element level beyond what is modeled in Born-level QQ production. With presently available MC technology, this can be achieved in different ways and at various formal accuracies. We now briefly describe the several prescription used in this study.
-Leading Order Multijet Merging: The LO multijet/multileg merging techniques [44][45][46] outline how parton shower emissions can be augmented with full matrix elements. The emissions are classified according to their hardness, i.e., p T , and in terms of a dimensionful variable Q cut . Emissions above a hardness threshold Q cut are described at LO accuracy using the appropriate matrix elements while preserving the all-orders resummation accuracy of the parton shower below Q cut . In this work we use the MLM scheme [45] as implemented in MadGraph5_aMC@NLO 2.6.0 (MG5_aMC@NLO) [19]. We take into account matrix element corrections to QQ pair production in association with up to two QCD In order to achieve NLO+LL accuracy for leading-jet observables, the above NLO+PS matching technique must be applied to the QQ + j process. This requires explicitly introducing a regularizing p T selection on the leading jet at the matrix element level. Since no Sudakov form factor is present for this jet, the p T selection must also be well above the Sudakov shoulder of the inclusive QQ-system to ensure that the FO is perturbatively valid. -Next-to-Leading Order Multijet Merging: The LO multijet merging technique described above can be extended to describe jet observables at NLO precision for jets above a hardness scale Q cut . Analogous to the LO case, matrix element merging with parton shower matching at next-to-leading order, known colloquially as MEPS@NLO, is achievable by introducing additional all-orders/resummed Sudakov form factors for each NLO-accurate matrix element in consideration.
We use an extension of the Catani-Krauss-Kuhn-Webber (CKKW) [56][57][58] merging formalism as implemented in Sherpa by Refs. [59,60]. In the following, we employ MEPS@NLO multijet merging with up to one or two jets, meaning that samples will contain up to two or three real radiations, respectively, beyond the lowest order process before parton showering.
where H T is the scalar sum of the transverse energy of the final state partons and Q. In Sect. 3, where we discuss experimental searches and sensitivity to Q in monojet searches, we impose an analysislevel selection on the leading jet ( j 1 ) in an event. More specifically, we set the ptj variable in MG5_aMC@NLO's run_card.dat file to ptj= p j 1 T,cut − 100 GeV (and jet radius R = 0.6). Then, to enhance yields at relatively high p T , we generate events by binning the phase space in p T . After preparing event samples for a particular p j 1 T,cut , we apply higher p j 1 T,cut until the statistical uncertainty of the MC samples is no longer negligible. At this point, we prepare another event sample based on p j 1 T,cut and apply the procedure iteratively. For the samples with the highest p j 1 T,cut , we apply exponential biasing on p j 1 T at the event-generation level to enhance the tail of jet p T distributions.
After the QQ pair have been decayed, events are passed to Pythia 8.2.26 [25] for parton showering and hadronization. We choose a shower starting scale Q S small enough such that light, colored final-state partons in the matrix element remain the hardest emissions in the full process, if the parton exists. Namely, the parton shower is restricted to operate below the scale [19] Here, d * is the minimum d i, j distance measure as calculated with the k T algorithm [64,65] (for R = 1) over all momentum recombinations of light, colored partons during the clustering phase in fixed order event generation. The events are essentially categorized by whether or not a hard O(α s ) emission is present.
To quantify and estimate the size of missing, higher order QCD corrections, we compute the three-point scale-variation envelope. This is obtained in the usual fashion, i.e., by varying discretely and jointly the factorization and renormalization scales Q R and Q F over the range, Where necessary, we also consider the uncertainty associated with the parton shower starting scale. We quantify this by discretely and independently computing the scale variation over the range,

Event generation in Sherpa
For cross-validation of the fermionic top quark partner T p , we employ Sherpa 2.2.4 [26]. At LO in QCD, arbitrary BSM models can be simulated through Sherpa's generic UFO model [31] interface [66]. At NLO, processes involving T p can be simulated by slightly modifying the default SM model file, and setting the top quark mass to the mass of T p . For the decay into a scalar dark matter particle, an additional decay vertex is added using the methodology of Ref. [66]. Tree-level matrix elements of the calculation are provided by Sherpa's in-house matrix element generators AMEGIC [67] and COMIX [68]. One-loop amplitudes are treated by interfacing with OpenLoops [69]. Parton shower-matching is performed according to the MC@NLO formalism [53,70,71], using Sherpa's Catani-Seymour subtraction-based shower procedure [72,73].
To study potential improvements to modeling fermionic top quark partners, we also employ multijet merging at NLO in QCD with Sherpa. To account for additional highp T QCD emissions at the matrix-element level, beyond what is already present in inclusive T pTp production at NLO, we include the NLO-accurate matrix elements for the processes pp → T pTp + j and pp → T pTp + 2 j.
We merge these samples with the fully inclusive pp → T pTp sample following the MEPS@NLO prescription [59,60]. In this scheme, the nominal values for the factorization, renormalization, and parton shower scale are determined through a backward clustering procedure that maps higher multiplicity configurations to a 2 → 2 configuration. We set all scales to the invariant mass of the top partner pair. As a nominal value for the merging scale Q cut , we set Q cut = 120 GeV. For parton showering we employ one of Sherpa's dipole showers, which is published in [73].
In addition to MEPS@NLO merged samples we also simulate pp → T pTp + j at NLO+PS with Sherpa. We do not add matrix element corrections to these samples beyond what is already present at NLO+PS. For the generation of the pp → T pTp + j process, we use the scale schemes, where m(T pTp ) and p T (T pTp ) denote the (T pTp )-system's invariant mass and transverse momentum respectively.

Detector simulation and object reconstruction
For fast detector simulation, we use Delphes 3.3.3 [74] with the default ATLAS card. Jets are constructed from calorimeter tower elements using Fastjet 3.2.1 [75], according to the anti-k T jet clustering algorithm [76] with jet radius R = 0.4. Analysis-quality jets are required to satisfy the fiducial and kinematic criteria, p j T > 20 GeV and |η j | < 4.5.
Events are then accepted or rejected based the monojet selection criteria discussed in Sect. 3.

Standard model inputs
We assume n f = 5 active/massless quark flavors and a diagonal Cabbibo-Kobayashi-Maskawa (CKM) quark mixing matrix with unit entries. Relevant SM inputs used in our study include, The evolution of parton distribution functions (PDFs) and the strong coupling constant α s (μ R ) are extracted using the LHAPDF 6.1.6 [77] libraries. As discussed in Sect. 2.2, LO (NLO) multijet merging with two (one) additional partons accounts for new kinematic channels and configurations that first arise at NNLO. (In principle, all one would need to achieve NNLO accuracy are the missing two-loop virtual corrections.) However, such O(α 2 s ) contributions are already accounted for in the normalizations of NLO PDFs. Hence, to minimize potential double counting of initialstate contributions, we use the NNPDF 3.0 NNLO PDF set (lhaid=261000) [78] for all signal process calculations. For the LO SM pp → Z + j calculation in Sect. 3, we use the NNPDF 3.0 LO PDF.

Monojet searches at the HL-and HE-LHC
At hadron colliders, the term "monojets" represents a broad class of sensitive collider signatures and search strategies that assume varying degrees of particle multiplicity and inclusiveness. In this section, we consider specifically the inclusive monojet signature, as implemented by ATLAS during Run II of the LHC's operations after collecting L =36.1 fb −1 of integrated luminosity at √ s = 13 TeV [18]. After discussing various sources of experimental and theoretical uncertainties, we report the observed and expected sensitivity of the channel at current and proposed pp colliders using the CL s modified frequentist approach [79]. The model-independent limits derived in Sect. 3.1 are then applied in Sect. 3.2 to the heavy colored particles Q described in Sect. 2.1. The ability to determine and distinguish principle properties of Q is then discussed in Sect. 3.3.

Monojet searches, LHC data, and model-independent limits
The ATLAS and CMS collaborations have both reported on their search of early Run II data for anomalous events with significant transverse momentum imbalance and at least one energetic jet [18,80]. For the case under consideration [18], the ATLAS collaboration has investigated two overlapping signal regions, categorized as exclusive modes (EM) and inclusive modes (IM), based on various E miss  For the remainder of this study, we focus on the IM monojet signal regions.
In Table 2, we display the expected number of SM (background) events passing all selection criteria with uncertainties (statistical and total) in each of the inclusive mode signal regions (IM1-IM10), as reported by ATLAS [18]. Nonstatistical errors include both experimental and theoretical uncertainties. Sources of systematic uncertainty estimated in Ref. [18] include: dependencies on parton shower and PDF modeling, which span ± 0.7 to ± 0.8%; uncertainties in jet energy and E miss T scales, which range ± 0.5% (IM1) to ± 5.3% (IM10); jet quality and pileup descriptions additionally provide uncertainties ranging ± 0.8 to ± 1.8%; for more detailed discussions, see Ref. [18]. We note that the total errors for IM1-5 are nearly flat, with 2.2-2.6%, indicating that these signal regions' uncertainties are systematics dominated. Statistical and systematic uncertainties are much larger in the higher E miss T regions. However, with the HL phase of LHC, one expects to collect significantly more data that will correspondingly reduce statistical errors for the high-E miss T regions. Additionally, the analysis' control sample will also increase during the HL run, therefore also reducing systematic uncertainties. Thus, one anticipates that total uncertainties will shrink for future inclusive monojet searches at LHC.
The non-observation of data deviating significantly from SM predictions enables one to set model-independent upper limits on the production cross section of new particles. In the second and third columns of Table 3, we tabulate the expected and observed limits at the 95% confidence level (CL) on the visible cross section, 2 respectively, for each IM signal region, as reported in [18] with L =36.1 fb −1 of data. To address the prospect of the HL-LHC with L = 3 ab −1 at √ s = 13 TeV, we calculate the expected upper limits by scaling the number of events at the 13 TeV run of LHC. We choose two val-2 Explicitly, the visible cross section is defined as the product of total cross section, acceptance, and efficiency. ues of the total systematic uncertainty, namely δ Sys. =2.5% and 1%. The former is a pessimistic assertion that systematic uncertainties will not be reduced beyond present levels (see Table 2), even after 15-20 years of LHC operations. The latter is an optimistic, but benchmark, assumption. The likelihoods of background only and signal-plus-background hypotheses are set as Gaussian, with a standard deviation set to the total uncertainty. We have checked that our likelihoods are in agreement with those reported by Ref. [18] for L = 36.1 fb −1 . We report the scaled limits in the last two columns of Table 3. We now discuss the impact of these limits on the production of QQ pairs in pp collisions.

HL-and HE-LHC sensitivity to heavy colored particles Q
We now compare and apply the model-independent upper limits on the cross sections derived in the previous section to the production of QQ pairs at the LHC. In Fig. 2, we show the model-independent 95% CL upper limits along with NLO+PS-accurate cross section for Q produced in association with a hard jet at the matrix-element level, viz. QQ + j at NLO. (For details of our computational setup, see Sect. 2.) We overlay FO scale uncertainty, computed according to Eq. (6). For all cases of Q, we find that the scale uncertainty at NLO is a dominant source of uncertainty and hence take it as a representative measure of the total uncertainty. A dedicated and in-depth discussion of this uncertainty is given in Sect. 4.
From inclusive monojet searches at √ s = 13 TeV with L = 36.1 fb −1 of data, we find that the lower limits on Q masses stand at around m T p = 400 GeV for the fermionic top partner and mg = 600 GeV for the gluino, while no constraint on stop masses is found within the range under consideration. We observe that in the highp T bins both the systematic and statistical experimental uncertainties play a crucial role. As argued in the previous section, one expects Table 3 Model-independent 95% CL upper limit on the visible cross section for each inclusive mode (IM) signal region, after L =36.1 fb −1 of data at √ s = 13 TeV, as reported by Ref. [18], and the estimated CL s limit assuming L =3 ab −1

IM
Observed limit sensitivity to improve at the HL-LHC due to a much larger dataset, leading to better control on both uncertainties. From the scaled limits, we find that fermionic top partners with masses m T p 800 GeV, gluinos with mg 1000 GeV, and stops with masses mt 600 GeV, in a compressed spectrum scenario, can be excluded at 13 TeV with L = 3 ab −1 , using the inclusive monojet signature. Along with more data, a possibility that can greatly push the sensitivity to heavy colored particles is increasing the beam energy of the LHC itself. Presently, community discussions are underway on upgrading the LHC's magnet system to handle a center-of-mass energy up to √ s = 27 TeV [81]. In light of this prospect, we briefly investigate the impact of a higher beam energy on the production of Q and the SM Z + jets background, i.e., the dominant background of the monojet signature [18], and estimate the experimental reach of such a collider.
Before this, however, we shortly digress to describe our modeling of the SM Z + jets background. In particular, we note that while the signal processes are consistently determined at NLO in QCD with parton shower matching, we consider the Z +jets background only at LO with an experimentally determined normalization factor. For our purposes, we believe this provides a sufficiently reliable description of the SM backgrounds after selection cuts.
The motivation comes precisely from the fact that proper modeling of the SM EW boson+jets background for monojet searches is highly nontrivial [82], particularly in comparison to simulating the inclusive W/Z + jets process. The technical difficulty is due, in part, to strong phase space restrictions (cuts) on the final-state jets (see above Sect. 3.1 for the list of cuts). Present implementations of the MC@NLO formalism into general purpose event generators require that one integrates over the entire phase space of additional real radiation at NLO to ensure infrared pole cancellation. In the present case, this renders event generation at NLO inefficient. For example: the pp → Z + j process with p j T 250 − 1000 GeV is known [83,84] to have giant QCD corrections stemming from the opening of new kinematic configurations. In this instance, large corrections are driven by the highp T dijet process with a relatively soft Z emission off a finalstate quark, a configuration that would otherwise fail the Δφ( jet, E miss T ) and minimum E miss T selection criteria. Moreover, after careful consideration, we find that the EW boson + jets background, that survives the selection analysis, is dominated by Born-like configurations. For example: after cuts, partonic channels such as qg → q Z contribute much more to the pp → Z +jets background than channels like qq → qq Z.
Once selection criteria have been applied, the difference then between our background modeling and a much more precise determination, e.g., Ref. [82], is largely an overall normalization. This does not necessarily hold true when taking into account EW corrections at NLO and beyond. However, such corrections are beyond the claimed accuracy of our work and we refer readers to Refs. [83,84] Fig. 2 The pp → QQ + j cross section as a function of minimum E miss T after the experimental selection criteria at 13 TeV, for Q = T p ,g, andt, with current 95% CL s limits after L = 36.1 fb −1 of data at the 13 TeV LHC. Also shown is the estimated sensitivity with L = 3 ab −1 , assuming δ Syst. =2.5 and 1% systematical errors of the pp → Z + j rates at NLO+PS and LO+PS after selection cuts are applied is roughly a constant K NLO ∼ 1.2. This is consistent with the size of finite virtual corrections at NLO in QCD to the DY process [85], and supports our arguments that the residual EW+jets background exhibits Born-like kinematics. Now, to achieve a reliable normalization of the EW+jets background, we scale the generator-level pp → Z + j cross section at LO by (approximately) a factor 1/20 so that the ratio of the curve to the signal cross section is normalized with respect to the post-event selection limit in  Table 4, we list the cross sections for the Z + j process at √ s = 13, 14, and 27 TeV for representative p gen T selections. We stress that this procedure is only an estimation of the SM background; we do not advocate that this is a suitable replacement for full NLO+PS (or more accurate) computations in experimental searches. For representative masses, m Q , we show in Fig. 3 the cross sections for the pp → QQ + j and (normalized) Z + j processes as a function of the leading jet p T . More specifically, the cross sections are calculated as a function of a generator-level p T threshold ( p gen T ) on the light jet. In the lower panel, we show the 27 TeV-to-14 TeV cross section ratios. There, one sees that the production cross section of Q increases faster than the SM background with increasing center-of-mass energy. The enhancement follows from the well-documented [86][87][88] growth in PDF luminosities for fixed partonic mass scales but increasing collider beam energy. Quantitatively, for p gen T = 1 TeV, the Z + j cross section increases by a factor of 10 with respect to the change of √ s, while the fermionic top partner cross section for m T p = 600 (800) GeV increases by approximately 24 × (28×), the gluino rate for mg = 900 GeV by 40×, and the stop rate for mt = 400 GeV by 20×. Although other sources of SM backgrounds for the monojet signature have been presently neglected, the signal over background ratio (S/N ) still increases significantly at higher collider energies due to the larger luminosity enhancement for the signal process than dominant SM backgrounds. Subsequently, the HE-LHC enables ones to investigate parameter regions that are not accessible at the LHC. We do emphasize, however, that S/N ratios can change drastically if additional information is provided to enhance the separation of the signal events from the backgrounds. For example: proposals exist on how to utilize soft leptons, jets, and displaced vertices associated with decays of Q that can further reduce SM background rates [89][90][91].
As shown in Fig. 4, once the scaled limits for the SM backgrounds are determined, we can compare the predicted cross sections for pp → QQ + j NLO process and estimate the expected reach at the 27 TeV HE-LHC. We find that with L = 3 − 15 ab −1 , one is sensitive to compressed spectra scenarios featuring fermionic top partners with masses m T p 1100 GeV, gluinos with masses mg 1800 GeV, and stops with masses mt 600 GeV.
We end this discussion by providing an estimate of the anticipated sensitivity of new colored particles Q at the 27 TeV HE-LHC. Here, we assume that the SM background is dominated by the Z + j process and simply scale the modelindependent 95% CL upper limit at √ s = 13 TeV by the 27-to-13 TeV production cross section ratio. In other words, the SM background cross section at a given √ s and E miss T,cut , denoted as σ ( √ s; E miss T,cut ), is obtained from the relation In the above, σ pp→Z j ( √ s; p gen T ) is the LO Z + j cross section with p gen T at a collider energy √ s. We further assume that the detector acceptance and efficiencies are the same at 13 and 27 TeV. This assumption is not as strong as one may anticipate in more general circumstances. The HE-LHC project proposes to refit, replace, and/or upgrade the current LHC magnet system and detector experiments. As the detector experiment caverns themselves cannot physically grow, one is forced to adopt a detector fiducial volume at 27 TeV that is largely unchanged from 13 TeV. Similarly, we also assume systematic uncertainties of 2.5% and 1%, the same considered in Sect. 3

Properties determination of heavy colored particles
We now turn to the possibility of extracting properties of the heavy colored particle Q from jet behavior within the monojet signature. As briefly discussed in the introduction, asserting color representation and spin of Q is required to infer information on its mass from cross section measurements (or limits). Consequently, a single cross section measurement of a particular monojet signal region does not help much in determining the nature of Q. For example: in Fig. 3, one sees that the production cross section for the process pp → QQ + j with p gen T = 600 GeV, for a top squarkt with mt = 400 GeV, a gluinog with mg = 900 GeV, and a fermionic top partner T p with m T p = 600 GeV are roughly within Δσ ∼ 5-10 fb of one another. However, despite this ambiguity, it is still possible to extract information from the pp → QQ + j cross section as a function of the leading jet p T , which can be measured directly, since it obeys a distinguishing pattern for each Q hypothesis. That the nature of Q is, in part, encoded in this observable reflects a nontrivial interplay between Q's mass, m Q , its color representation and spin, and the dimensionless ratio (m Q / p j 1 T ). This interplay is what we now discuss.
The first discerning observation is that the pp → QQ + j cross sections do not depend on p In the top (bottom) single ratio, both cross sections are with respect to the mass m Q (m * ) but different p j 1 T,cut . This has the effect of canceling overall color and kinematic factors while isolating logarithmic terms of the form log(m Q / p j 1 T,cut ) 2 . The double ratio, then, is a measure of this logarithmic dependence with respect to a baseline mass m * and minimum transverse momentum p * T . For Q = T p , we choose m * = 600 GeV and p * T = 200 GeV, and plot P Q ( p j 1 T,cut ) for m T p = 400 and 800 GeV. Quantitatively, one sees that the double ratio increases (decreases) by about 50% at p j 1 T,cut = 800 GeV for m T p = 800 (400) GeV. This feature is universal for particles in the same color representation and follows from the nature of massless gauge boson emission in scattering processes.
To better understand this behavior, consider the pp → QQ + j process. The t-channel propagators gives rise, after phase space integration, to the aforementioned logarithms log(m Q / p j 1 T,cut ) 2 . In the context of parton shower resummation, this dependence can be interpreted as the likelihood of emitting an additional QCD parton with transverse momentum p j 1 T,cut , i.e., the differential probability is proportional to Hence, a fixed probability ΔP implies a fixed (m Q / p j 1 T ) 2 ratio, and indicates that increasing m Q results in p j T increasing commensurately. Qualitatively, the emission of higherp T QCD partons becomes easier for heavier Q because highp T emissions become relatively soft as m Q increases. This results in a rightward shift of the so-called Sudakov shoulder [92,93]. For TeV-scale particles, the rightward shift of what constitutes "soft" is known to be large; see for example Refs. [94,95].
The discrimination power of cross section ratios extends if one considers the additional dependency on a collider's beam energy. In particular, we find that the cross section ratios shift from unity for p gen T = 600 GeV at √ s = 14 TeV, to σ (g)/σ (t) = 2.1 and σ (T p )/σ (t) = 1.35 at 27 TeV. Hence, a measurement of the ratio of the cross section at the √ s = 14 and 27 TeV with 30% accuracy can resolve our benchmark gluino, fermionic top parter, and stop scenario. Moreover, increasing the beam energy can also significantly improve the signal-over-background ratio, thereby enabling measurements of the QQ+ j cross section for different minimum p j T over a wide range of p T . Measuring the signal cross section for several p j 1 T,cut at √ s = 27 TeV with 10% accuracy allows one to distinguish the benchmark scenario by narrowing down the mass of Q without √ s = 14 TeV information. In light of this, it is necessary to emphasize that theoretical predictions on total and differential cross sections, as well as their ratios, must have the requisite accuracy to make these measurements.
Lastly, we note that the slope of σ ( p . The structure of R Q ( p gen T , m Q ) here is analogous to P Q ( p j 1 T,cut , m Q ) in Eq. (13), but differs in that the normalizing process is fixed to pp →tt * + j with p gen T = p * T = 300 GeV for √ s = 14 TeV, and we vary pp → QQ + j in the upper ratio. Subsequently, while overall color and kinematic multiplicative factors cancel, the relative dependence of q q g, and gg scattering within an individual process does not cancel and is inherently dependent on the color representation.
From Fig. 6, one can observe that fermionic top partners and stops with same mass have the almost same slope over a wide range of p T . On the other hand, because of color and matrix element effects from gluinos, the ratio Rg( p gen T ) tend to decrease with p T but again independent to the mass of the heavy particle. For validation, we consider the case for scalar gluons, known as sgluons [35][36][37], and overlay sgluon behavior on the same figure. 3 As the plot suggests, the distribution of the hard radiation associated with the colored particle production follows universally for a given color representation of Q.

Theoretical uncertainties of the monojet process
As investigated in Sect. 3.3, were one to discover new colored particles at the LHC, or a potential successor experiment, a measurement of the pp → QQ+ j cross section as a function of the leading jet's minimum transverse momentum ( p j 1 T,cut ) can help establish the quantum numbers of Q. Ascertaining such information, however, requires accurate BSM signal predictions. For the benchmark scenarios listed in Table 1, we find one needs theoretical uncertainties no larger than O(30%) and O(5%), respectively, on the total inclusive cross section normalizations and on the change of the cross section for σ ( pp → QQ+ j) as a function of p j 1 T,cut per 100 GeV. In this section, we discuss and quantify theoretical uncertainties associated with the monojet signal process. We particularly investigate (potential) sources of uncertainties when employing the state-of-the-art MC suites MG5_aMC@NLO+PY8 and Sherpa.
To investigate theoretical uncertainties associated with the QQ + j process, we focus on fermionic top partners T p . We choose this benchmark because of wide implementability across different event generators as well as its comparability to tt + j production in the SM. We compare several simulation setups and techniques at different levels of precision We begin with Fig. 7a, where we show the pp → T pTp + j cross section as a function of p j 1 T as derived from the fully inclusive pp → T pTp calculation at NLO+PS (dash), as obtained from the MC@NLO formalism. We also show the calculation at LO+PS (dot), as obtained with up to two additional matrix element-level jets via multijet merging. The curves are shown with their factorization (Q F ) and renormalization (Q R ) variation envelopes and are normalized to the respective nominal prediction. As discussed in Sect. 2, the leading jet p T distribution for both calculations is only LO accurate for p j T comparable to the scale of the hard process. Hence, the scale uncertainties for p j T Q F , Q R are large, asymmetric, equal for the two calculations, and span approximately + 45 to − 30%. For smaller p j T , namely p j T m T p = 600 GeV, one observes that the two uncertainty envelopes begin to differ. Whereas the uncertainty for the multijet calculation reduce only slightly for decreasing p j T , the MC@NLO uncertainty reduces to roughly the ± 10% level. The difference originates from virtual corrections present in the NLO+PS calculation, which soften dependencies on Q F , Q R , but are obviously absent in the LO+PS calculation. The Sudakov-like factor in the multijet merging prescription only partially reduces the dependence on Q F by matching lowp T QCD emissions in the hard matrix element with those in the PDF. Hence, even for observables that are formally of the same precision, the presence of all O(α s ) terms in the NLO+PS calculation leads to The small scale uncertainty observed for the lowest p j 1 T suggests high theoretical precision is achievable with NLO+PS computations obtained via the MC@NLO prescription. However, unlike pure FO calculations, calculations matched to parton showers possess the additional dependence on the parton shower starting scale Q S . Within the MC@NLO formalism, Q S controls whether the leading O(α s ) emission beyond the Born process is included in the FO matrix element or the all-orders parton shower. Loosely speaking, QCD radiations with p j T above (below) Q S originate from the hard matrix element (parton shower). As pointed out in Refs. [96,97], lowest order-accurate observables, e.g., the p j 1 T distribution when the pp → T pTp process is evaluated at NLO+PS, and processes that possess large virtual corrections suffer from ambiguities when choosing Q S . This manifests as a strong dependence on Q S , and hence a large deviation.
We assess this behavior in Fig. 8 by plotting the pp → T pTp + j cross section, derived from the inclusive pp → T pTp calculation at NLO+PS, as a function of p j 1 T,cut . We assume multiplicative variations of the default parton shower scale, given in Eq. (7), but fix the factorization and renormalization scales to their central values. Rates are normalized to the fixed order NLO (fNLO) prediction. The Q S variation of the NLO+PS result with respect to the default Q S choice spans roughly ± 25%, with the dependence increas- interpreting the vanishing shower scale uncertainty. In this limit, the FO calculation is unphysical. The FO calculation possess an integrable singularity at p j T = 0 that leads to arbitrarily large cross sections. As a result, the ratio of the NLO+PS-level cross section, a finite and physical quantity, to the fNLO cross section, the unphysical quantity, vanishes as p j T /m T pTp → 0. For the curves with shower scale multiplier 0.5 and 1, we observe the curves are flat beyond ∼ 300 GeV and 600 GeV, respectively, as expected. However, the curves with shower scale multiplier larger than 1, we observe a crossing at p T ∼ 700 GeV. The large shower starting scale introduces high p T initial state radiations, which is neither soft nor collinear to the incoming parton. If the high p T jets generated from the parton shower of Born-type S events immersed in this region, then the cross section is not reliable. Therefore, the curves (and the crossing) are not trustworthy beyond the top partner mass, or equivalently √ŝ /2. A better way to estimate the pp → T pTp + j cross section (and associated shower scale uncertainty) would be to consider the process pp → T pTp + j at NLO+PS.
In light of the large theoretical uncertainties in the pp → QQ + j cross section stemming from Q F , Q R , and Q S , it is clear that the precision achieved with the aforementioned methods is insufficient for distinguishing Q candidates. To explore if such precision is yet still possible with presently available general-purpose Monte Carlo technology, we consider the cross section obtained from the pp → T pTp + j process itself at NLO+PS. In Fig. 7b we plot the normalized pp → T pTp + j cross section with its Q F , Q R envelope again as a function of p j 1 T,cut , but as derived from the pp → T pTp (dot) and pp → T pTp + j (dash) calculations at NLO+PS. We observe that the uncertainty reduces to a largely uniform band that is just shy of ± 20% for the T pTp + j calculation. Both uncertainties are comparable for p j 1 T m T p ; for p j 1 T m T p , however, the T pTp + j uncertainty is about 40-60% smaller.
It is worth reiterating that the FO T pTp + j calculation requires a minimum p j T selection to regulate matrix element poles and, since no Sudakov suppression is applied to low p T QCD emissions in this case, to ensure perturbative stability, in the Collins-Soper-Sterman (CSS) sense [92]. We safeguard against the need for k T -resummation on the leading jet by using the (somewhat conservative) CSS consistency requirement from Ref. [95]. For m T p = 600 GeV, and hence m T pTp > 1.2 TeV, one needs p j 1 T,cut 200 − 250 GeV to force collinear logarithms from t-channel exchanges, which scale as δσ ∼ α s log(m T pTp / p j T ) 2 , to be much smaller than 1.
There is now an apparent conflict between theoretical needs and computational capabilities: While pp → T pTp + j at NLO+PS provides improved control and stability over Q F and Q R , the calculation is only meaningful for sufficiently large p j 1 T . The pp → T pTp process at NLO+PS, on the other hand, extends the support for p j 1 T to low p T but at the cost of a larger uncertainties. Such demands, however, are precisely resolved with multijet merging at NLO+PS, i.e., MEPS@NLO. As the CKKW MEPS@NLO prescription is natively available in the event generator Sherpa, for the remainder of this section we use Sherpa to further quantify scale uncertainties in T pTp production. However, before proceeding, we non-trivially demonstrate that at least up to the QCD order presently being investigated, both MG5_aMC@NLO+PY8 and Sherpa predictions are in agreement with one another.
In Fig. 9a, we plot the pp → T pTp + j cross section and its scale uncertainty band as a function of p j 1 T,cut as derived from the pp → T pTp + j process at NLO+PS using MG5_aMC@NLO+PY8, as well as from the pp → T pTp process at NLO+PS merged with up to one additional NLO-accurate, matrix element-level jet using Sherpa. A theoretical difficulty of the inclusive monojet signature is ensuring perturbative control over predictions, despite complicated, multi-scale requirements on jet momenta and particle multiplicity. Hence, it is important to ensure that subleading jets in the monojet processes can be modeled as well as the leading jet. To investigate whether such uncertainties for subleading jets are attainable, we compare MEPS@NLO uncertainties for the pp → T pTp process when merged with different jet multiplicities. As a control, we show in Fig. 10a the (normalized) pp → T pTp + j cross section and uncertainty band as a function of p j 1 T,cut for the pp → T pTp process at NLO+PS merged with up to one (dark curve) or two (light curve) additional NLO-accurate, matrix element-level jets. No significant difference between the two curves is observed. One does not expect to see such deviations for this observable as the two calculations are identical at this order of α s . The difference appear at one order higher: At the matrix element level, the two calculations possess up to two and three hard QCD emissions, respectively. In the 1 j MEPS@NLO calculation, the second jet is LO+LL (alternatively, LO+PS) accurate. In the 2 j MEPS@NLO calculation, the second jet is NLO+LL accurate and the third jet is LO+LL accurate. Therefore, in Fig. 10b, we plot the pp → T pTp + 2 j cross section as a function of minimum p j 2 T for the two calculations. In analogy to the comparison of T pTp and T pTp + j production at NLO+PS in Fig. 7b, we observe here that the 1 j MEPS@NLO prediction suffers from an uncertainty spanning roughly + 45 to − 30% whereas the 2 j MEPS@NLO rate exhibits a largely uniform uncertainty of about ± 22%. For vanishing p j 2 T,cut , the uncertainties become equal since the observable is no longer sensitive to such a high order of α s . With successive MEPS@NLO matching to higher jet multiplicities, one expects a comparable reduction of dependence on Q F , Q R in the corresponding cross section. However, one should not expect the uncertainty to drop below the observed O(10 − 15%) without first accounting for missing two-loop terms at O(α 2 s ). We now turn to the issue of parton shower scale uncertainties for processes evaluated with MEPS@NLO. Following the same variation procedure as before, we plot in Fig. 11a, b, respectively, the Q S dependence for the T pTp + j and T pTp + 2 j cross sections as a function of minimum p T for the leading ( p j 1 T ) and subleading ( p j 2 T ) jet. The upper (lower) panel corresponds to the 1 j (2 j) MEPS@NLO prediction. All curves are normalized to the appropriate cross section evaluated at the default scale choice. We find for both calculations that the uncertainty in the T pTp + j cross section is very small, reaching maximally to ± 4% over the range of p j 1 T considered, and comparable to our MC, i.e., statistical, uncertainty. Relative to the T pTp at NLO+PS uncertainty in Fig. 8, which exhibits deviations up to ± 25%, this represents significant reduction in Q S dependence. For the T pTp + 2 j cross section, we see a larger scale dependence in both the 1 j and 2 j calculations. The uncertainties for the two rates are comparable, reaching about + 10 to − 5%, and again comparable to our MC uncertainty.
Having established that multijet-merged simulations with NLO matrix elements feature moderate-to-small factoriza- tion, renormalization, and shower scale dependencies, we lastly investigate the uncertainty associated with our PDF input and the MEPS@NLO merging scale. In Fig. 12a, we show again the normalized pp → T pTp + j cross section, as derived from the T pTp process with MEPS@NLO merging to one additional jet, and the associated 1σ PDF variation band. The uncertainty is derived from PDF replicas, following the procedure of Ref. [77]. For the range of p j 1 T,cut investigated, we find that the band is stable, symmetric, and spans roughly ± 6% at low p j 1 T,cut to ± 10% at high p j 1 T,cut . At 27 TeV, we largely investigate a comparable range of Bjorken-x as we do at 13 TeV, and hence expect the PDF uncertainty to remain the same due to the stability of DGLAP evolution. In Fig. 12b, we show the uncertainty associated with the MEPS@NLO jet merging scale, Q cut . Over a considerable range of p j 1 T,cut , we observe variations below ± 5% for p j 1 T,cut < 400 GeV and below 10% for p j 1 T,cut < 1 TeV, and is comparable to our MC uncertainty.

Outlook
In Sect. 3.3, we discussed the dependence of the pp → QQ + j cross section on the mass, spin, and color representation of Q. A single signal cross section measurement does not, of course, constrain the nature of Q uniquely because different spin and color hypotheses can lead to identical cross sections if the mass is chosen/tuned accordingly. For example: in the representative case above, a stop of mass 400 GeV, a fermionic top partner of mass 600 GeV, and a gluino of mass 800 GeV have practically the same pp → QQ+ j cross section for p j 1 T,cut ≈ 600 GeV at √ s = 14 TeV. The degeneracy can be resolved, however, through additional cross section measurements with larger p j 1 T,cut and/or at a higher center-ofmass energy. Due to the different masses of the representative cases, the cross sections scale differently with increasing jet transverse momentum and the energy. This leads to a break in the degeneracy, thereby enabling one to discriminate against the various Q hypotheses.
In the left panel of Fig. 13, we show that all three cases could be distinguished through a secondary measurement of the pp → QQ + j cross section at p j 1 T,cut = 1 TeV were our predictions were precise at the 10% level. The nominal difference between the gluino and fermionic top partner case is roughly that magnitude, while for the case of scalar top, the difference could reach 50%. In light of the uncertainties estimated in Sect. 4, it is unrealistic to claim that one can distinguish a gluino from top partner in this manner, with renormalization and factorization scale uncertainties alone reaching about 20%. Additionally one can also make use of the higher energy collisions, viz. the potential 27 TeV upgrade of LHC, to distinguish these three particles. As shown in Fig. 13b  Alternatively, instead of considering cross sections themselves, one could also consider ratios of the form As discussed in Sect. 3.3, ratios of cross sections measured at √ s = 14 TeV and √ s = 27 TeV can also effectively discriminate against various Q hypotheses if theoretical uncertainties on these ratios are smaller than about 30%.
In ratios and double ratios, such as those discussed in previous sections, the normalization component of uncertainties cancel. As shown in Fig. 9a, for example, variations of the factorization and renormalization scales affect mainly the overall normalization, not the shape/ p j T dependence. Therefore, uncertainties estimated in this way are expected to drop out in Eq. (15). To verify this, we calculate the scale uncertainty for the ratio of the cross section to the nominal scale choice: . (16) For the normalizing cross section, we consider the m T p = 600 GeV and vary the factorization and renormalization scales in the same manner and choose p = 300 GeV. We observe that the effect of varying μ R,F between 1/2 and 2 with Q 0 = H T /2 is marginal (only 2%) at 1 TeV. The PDF uncertainty, however, impacts the cross section more, though the effect remains no more than 10% at 1 TeV.
In Eq. (16), we assume the scale choices are correlated for both the numerators and denominators. Fairly, one might argue that in order to estimate uncertainties on cross section ratios through the shape of the p j 1 T,cut dependence, one should consider different functional forms of Q R and Q F as well. In Fig. 14 we show that the effect of choosing a nondynamical nominal scale can in fact induce significant slope variations. Unsurprisingly, deviations are observed for both MG5_aMC@NLO+PY8 and Sherpa events. The ratio of the cross section at 1 TeV and 300 GeV for Q 0 = m T p is about 15% higher than Q 0 = H T /2. The alternative scale choice Q R = Q F = m T p is, however, also quite extreme in the sense that it is completely static [98]. For example: the upper envelope for the scale variation (i.e., μ R = μ F = 0.5) is calculated with Q R = Q F = 300 GeV, but the event with p j 1 T,cut = 1 TeV implies √ s > 2 TeV. Therefore, the ratio √ s/Q R is a factor of 7, which leads to large non-asymmetric effects as seen in the right panel of Fig. 14 We conclude that cross section ratios, in general, will provide a more robust way for discrimination but also that the assessment of theoretical uncertainties on such ratios is a highly nontrivial issue. It is natural to expect a majority of the scale uncertainty will cancel in the cross section ratio. However, in general, an NNLO calculation can induce correction terms that scale with jet p T , and result in an effect may not be small. We have seen in Fig. 9a that different reference scale choices affect the cross section ratio at the level of 10%. The assessment of theoretical uncertainty on such ratio can only be estimated reliably by incorporating higher order calculations. A more quantitative study of this problem is, however, beyond the scope of this paper.
As a final comment, we note that more precise predictions for cross sections can in any case help to improve the discriminating power of the measurements proposed here. Presently, NNLO in QCD predictions for the relevant processes (which include an additional jet in the final state) are not available. For the inclusive process pp → T pTp , however, NNLO predictions can be obtained from HATHOR [34] and top++ [22]. In Fig. 15 we show that the uncertainties on the total cross section are approximately halved when going from NLO to NNLO in α s . Note that the uncertainty for the pp → T pTp + 2 j LO matched cross-section is around 25-30%, while the uncertainty reduces to 10% at NLO. Moreover, the uncertainty in NLO cross-section estimated by three independent sources also show excellent agreement. (Note that for Sherpa, we use the scale scheme Q F = Q R = m(T pTp )/2 to approximate the  Cross section for the pp → T pTp process with additional light QCD partons at the ME level at leading, next-to-leading and nextto-next-to-leading orders in α s provided by the available MC tools NNLO uncertainty on the total cross section is only around 5%. Uncertainties on monojet signal cross section are expected to be larger than the uncertainty on the inclusive pp → T pTp process. In Fig. 7, we show that renormalization and factorization scale uncertainties are around 40% at LO and 20% at NLO, which is twice as high as for the inclusive pp → T pTp cross section. This is to be expected because the lowest order term of the total cross section is proportional to α 2 s in the inclusive case, while the monojet signal cross section is proportional to α 3 s . Assuming that a comparable relative gain can be achieved for pp → T pTp + j when going to NNLO, it might be possible to get to a theoretical precision of about 10% when NNLO calculations become available. Considering the large momentum transfers involved, EW corrections might have to be added as well in order to capture the effect of (real and virtual) EW Sudakov logarithms [84,[99][100][101].

Summary
We have studied the monojet collider signature arising from the process pp → QQ + j, where Q is a heavy colored particle that decays into an invisible particle X with mass m X close to the mass of Q (m Q ). In the limit where (m Q − m X )/m Q is small, the full process signature is monojet-like, and the discovery prospects for Q and X become much more challenging, even at the HL-LHC or proposed HE-LHC, due to large SM backgrounds.
In this context, we have investigated the feasibility of extracting the properties of Q were a monojet signature discovered, focusing on the observation that the monojet signal cross section σ ( pp → QQ + j) is sensitive to the mass, spin, color representation of the particle. Due to an interplay of these quantities, one cannot readily and uniquely determine the nature of the particle Q from a single cross section measurement alone without assuming the mass of the particle.
We have studied several processes calculated at next-toleading order in QCD with parton shower matching and multijet merging using the state-of-the-art Monte Carlo suites MG5_aMC@NLO+PY8 and Sherpa. It has been pointed out that the dependence on p j T in the pp → QQ + j cross section, up to overall normalization, depends considerably on the mass of Q but exhibits a much milder dependence on its color representation. Moreover, we find that there is little-tono spin dependence in the shape of the cross section curves. The dependencies of the overall normalization and p j T on the underlying parameters can, in principle, allow efficient discrimination of the nature of Q, and hence the underlying physics behind the monojet signature.
Were the process discovered, there are many obstacles to extracting the cross section of the monojet signature and its p T dependence. First, the signal is overwhelmed by the background. Therefore, large statistics and a precise understanding of the background is required to observe the signature at the HL-LHC for the parameter region that has not been excluded by current LHC data. Increasing the collider energy improves sensitivity significantly, and represents a promising possibility. In addition, theoretical predictions of the overall cross section normalization need to be sufficiently accurate in order to differentiate against various candidate scenarios. For the relevant model parameter regions investigated, we observe that the overall cross sections need to be predicted within an uncertainty of Δσ/σ < 30%, while p j T depen-dencies should be known within δσ/δp j T < 10%/100 GeV. Finally, we find that the achievable precision for such quantities is comparable to the required precision.