Monojet Signatures from Heavy Colored Particles: Future Collider Sensitivities and Theoretical Uncertainties

In models with colored particle $\mathcal{Q}$ that can decay into a dark matter candidate $X$, the relevant collider process $pp\to \mathcal{Q}\bar{\mathcal{Q}}\rightarrow X\bar{X}+$jets gives rise to events with significant transverse momentum imbalance. When the masses of $\mathcal{Q}$ 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 $\sqrt{s}=14\,\mathrm{TeV}$ with $\mathcal{L}=3\,\mathrm{ab}^{-1}$ of data and the proposed High-Energy LHC at $\sqrt{s}=27\,\mathrm{TeV}$ with $\mathcal{L}=15\,\mathrm{ab}^{-1}$ of data. We estimate the reach for various Lorentz and QCD color representations of $\mathcal{Q}$. Identifying the nature of $\mathcal{Q}$ is very important to understanding the physics behind the monojet signature. Therefore, we also study the dependence of the observables built from the $pp\to\mathcal{Q}\bar{\mathcal{Q}} + j $ process on $\mathcal{Q}$ 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 $\mathcal{Q}$ 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 most 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 Rparity conservation, then 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 g/q g/q g/q Q Q Figure 1: Diagrammatic depiction of QQ pair production with an extra hard QCD radiation in pp collisions.
strongly constrained by LHC data if Q has a mass around 1 TeV and below [7][8][9][10]. Such constraints on Q, however, 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 * .
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 Tp = mg, the pair production cross sections for these particles exhibit the hierarchy σ(tt * ) σ(T pTp ) σ(gg).
(1.1) 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 − 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 the heavy colored particle (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 Sec. 2 we provide in-dept detail of our computational setup. In Sec. 3, we discuss observed and expected sensitivity of monojet searches at present and hypothetical future facilities, and address various theoretical uncertainties in Sec. 4. A brief outlook on the impact of this work is discussed in Sec. 5, and conclude in Sec. 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. The section continues as follows: In Sec. 2.1, we briefly summarize the Q we consider from representative BSM scenarios. In Sec. 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 Secs. 2.3 and 2.4, and our detector simulation in Sec. 2.5. In Sec. 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 GeV, 600 GeV, 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].
We treat the Q decay in the narrow-width approximation, thereby decoupling its production from its decay into an invisible particle X. 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. Fort, X is a neutral fermion andt are decayed via an off-shell top quark to the four-body final state,t + t * → X + bqq . Forg, X is a neutral fermion andg are decayed to the three-body 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 (relatively) soft. Hence, the SM decay products of Q fail the criteria needed to be identified as the leading jet.
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 high-p 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 [38][39][40] 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-order resummation accuracy of the parton shower below Q cut . In this work we use the MLM scheme [39] 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 partons. While genuine O(α s ) (and higher) corrections are included via this procedure, the calculation remains formally LO accurate (LO+LL after parton showering) due to missing virtual corrections. Some of the earlier studies on monojet spectra with LO multijet merging technique using then available Madgraph/MadEvent v4 [41], we refer [42][43][44][45][46].
• QQ Production at Next-to-Leading Order with Parton Shower Matching: The accuracy of fixed order (FO) matrix element calculations at NLO can be combined with resummed parton showers at LL by means of NLO+PS matching techniques [47][48][49], such as the MC@NLO formalism [47]. In this approach, the LO matrix element for an extra hard, wide-angle QCD emission in the final state is naturally included as part of the O(α s ) FO correction. Extra soft and/or collinear emissions enter through the parton shower. Potential double counting of O(α s ) soft/collinear contributions is avoided by the use of additional counter terms. It is worth noting that leading-jet observables in this calculation are at most LO+LL accurate, but are nonetheless well-defined at all p T .
• QQ + j Production at Next-to-Leading Order with Parton Shower Matching: 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/resumed Sudakov form factors for each NLO-accurate matrix element in consideration.
We use an extension of the Catani-Krauss-Kuhn-Webber (CKKW) [50][51][52] merging formalism as implemented in Sherpa by Refs. [53,54]. In the following, we employ MEPS@NLO multijet merging with up to one or two jets, meaning that samples will contain up two or three real radiations, respectively, beyond the lowest order process before parton showering. As in the LO(+LL) case, while O(α 2 s ) corrections are present in this calculation, the final result remains formally NLO(+LL) accurate.

Event Generation in MadGraph5 aMC@NLO + Pythia8
Cross section calculations and event generation with accuracy up to NLO in QCD are handled using MG5 aMC@NLO. For signal processes, we use the NLO-accurate UFO libraries described in Sec. 2.1 and listed in Table. 1 as inputs to MG5 aMC@NLO. Within the MG5 aMC@NLO framework, one-loop virtual corrections are evaluated numerically via MadLoop [55] and matched with real emissions using the Frixione-Kunszt-Signer (FKS) subtraction formalism [56], as implement by Ref. [57]. Decays of Q are then handled at LO accuracy with MadSpin [58]. The central value Q 0 for the renormalization scale Q R and factorization scale Q F is set to where H T is the scalar sum of the transverse energy of the final state partons and Q. In Sec. 3, where we discuss experimental searches and sensitivity to Q in monojet searches, we impose an analysis-level 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 j1 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 j1 T,cut , we apply higher p j1 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 j1 T,cut and apply the procedure iteratively. For the samples with the highest p j1 T,cut , we apply exponential biasing on p j1 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 [59,60] (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 threepoint 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 [61]. 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. [61]. Tree-level matrix elements of the calculation are provided by Sherpa's in-house matrix element generators AMEGIC [62] and COMIX [63]. One-loop amplitudes are treated by interfacing with OpenLoops [64]. Parton shower-matching is performed according to the MC@NLO formalism [47,65,66], using Sherpa's Catani-Seymour subtraction-based shower procedure [67,68].
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 high-p 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 + 2j. (2.5) We merge these samples with the fully inclusive pp → T pTp sample following the MEPS@NLO prescription [53,54]. 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 [68].
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.

Detector simulation and object reconstruction
For fast detector simulation, we use Delphes 3.3.3 [69] with the default ATLAS card. Jets are constructed from calorimeter tower elements using Fastjet 3.2.1 [70], according to the anti-k T jet clustering algorithm [71] 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. (2.7) Events are then accepted or rejected based the monojet selection criteria discussed in Sec. 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.

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 [74]. The model-independent limits derived in Sec. 3.1 are then applied in Sec. 3.2 to the heavy colored particles Q described in Sec. 2.1. The ability to determine and distinguish principle properties of Q is then discussed in Sec. 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,75]. 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 T thresholds spanning E miss T = 250 GeV to 1 TeV. The EM signal regions are defined in terms of E miss T binning. For example: signal region EM1 (EM5) selects for events with E miss T satisfying 250 GeV < E miss T < 300 GeV (500 GeV < E miss T < 600 GeV). The IM signal regions are defined in terms of minimum E miss T selections. For example: signal region IM1 (IM5) selects for events with E miss T satisfying E miss T > 250 (500) GeV. Additionally, events are required to satisfy the following selection criteria: • At least one analysis-level jet with p T > 250 GeV and |η| < 2.4.
• A maximum of four analysis-level jets satisfying p T > 30 GeV and |η| < 2.8.
• An azimuthal separation of ∆φ(jet, E miss T ) > 0.4 for each analysis-level jet and the E miss T vector.  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 .
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]. Non-statistical 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% to 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 modelindependent 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 † , 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 values 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 † Explicitly, the visible cross section is defined as the product of total cross section, acceptance, and efficiency.
[GeV]       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. 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 Sec. 2.) We overlay FO scale uncertainty, computed according to Eq. (2.3). 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 Sec. 4. We observe that in the high-p T bins both the systematic and statistical experimental uncertainties play a crucial role. As argued in the previous section, one expects 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 Tp 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 [76]. In light of this prospect, we briefly investigate the impact on the production of Q and the SM Z + jets background, i.e., the dominant background of the monojet signature, and estimate the experimental reach of such a collider. For representative masses, m Q , we show in Fig. 3 the cross sections for the pp → QQ + j and 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. We scale the Z + j cross section (double-dot-dashed) by 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 Fig. 2. In the lower panel, we show the 27 TeV-to-14 TeV cross section ratios.
We briefly note that while the signal processes are obtained at NLO in QCD accuracy, the Z + j background is considered only at LO ‡ . This is due to the fact that jet observables in Z + j events ‡ State-of-the-art high-precision calculations for Z + jets process are already available in literature, see [77].  can possess very large QCD corrections due to the opening of new kinematic configurations [78]. In the present case, the pp → Z + j process with TeV-scale p gen T indeed exhibits such behavior. The increase corresponds to the real emission of low-p T Z bosons in high-p T dijet events. However, these configurations have relatively small E miss T , and thus not relevant to our monojet study. We have checked that the QCD K-factor lies within 1.2 for events after requiring large missing momentum at the NLO+PS level. As a result, LO predictions for the Z + j process represent a more appropriate description of the background than Z + j at NLO, particularly after the aforementioned normalization procedure. 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.
Returning to the lower panel of Fig. 3, 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 [79][80][81] 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 Tp = 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 [82][83][84].
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 model-independent 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 . (3.1) In the above, σ pp→Zj ( √ 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 Sec. 3.1.
Arguably, the background modeling of a LO process appears naive at first glance. However, the dominant SM backgrounds for inclusive monojet searches are indeed electroweak boson production [18]. Such processes possess an initial-state parton composition and color structure comparable to Z +j, and also exhibit a similar dependence on collider energy. Hence, 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 compressed spectra scenarios featuring fermionic top partners with masses m Tp 1100 GeV, gluinos with masses mg 1800 GeV, and stops with masses mt 600 GeV.

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 Tp = 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 j1 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 gen T in a universal manner. Keeping to Fig. 3, one sees that while σ(QQ + j; p gen T = 600 GeV) are the same at √ s = 14 TeV for the (Q, m Q ) configurations under consideration, the relative size of σ(QQ+j) changes with p gen T . In other words, while σ(QQ+j; p gen T ) follows an anticipated power-law of σ(p gen T ) ∼ (p gen T ) −β , with β > 0, the precise value of the exponent is dependent on the color and spin structure of Q. In a particular extreme, the gluino rate is the smallest (largest) of the configurations for p gen T smaller (larger) than p gen T = 600 GeV, suggesting a smaller β than for other Q. Information on Q can be extracted from σ(p gen T ) by considering its ratio with respect to a benchmark σ(QQ + j; p gen T ). For example: for the benchmark processtt * + j with mt = 400 GeV, the T observed for gluinos, fermionic top partners, and stops arise from the fact that heavier particles give rise to harder, i.e. less steeply falling, p j1 T distributions. The benchmark cases having significantly different input masses, with mg = 900 GeV, m Tp = 600 and 800 GeV, and mt = 400 GeV, and makes a significant impact on the p j1 T dependence. To isolate this behavior, in Fig. 5 we plot the cross section double ratio as a function of p j1 T,cut , .

(3.2)
In the top (bottom) single ratio, both cross sections are with respect to the mass m Q (m * ) but different p j1 T,cut . This has the effect of canceling overall color and kinematic factors while isolating logarithmic terms of the form log(m Q /p j1 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 j1 T,cut ) for m Tp = 400 and 800 GeV. Quantitatively, one sees that the double ratio increases (decreases) by about 50% at p j1 T,cut = 800 GeV for m Tp = 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 j1 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 j1 T,cut , i.e., the differential probability is proportional to dP = (1/σ)dσ ∝ α s (m Q ) log(m 2 Q /p j2 T ). Hence, a fixed probability ∆P implies a fixed (m Q /p j1 T ) 2 ratio, and indicates that increasing m Q results in p j T increasing commensurately. Qualitatively, the emission of higher-p T QCD partons becomes easier for heavier Q because high-p T emissions become relatively soft as m Q increases. This results in a rightward shift of the so-called Sudakov shoulder [85,86]. For TeV-scale particles, the rightward shift of what constitutes "soft" is known to be large; see for example Refs. [87,88].
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  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 j1 T ) with respect to p j1 T does not significantly depend on either the color structure or spin of Q. To show this, we display in Fig. 6 a second cross section double ratio, R Q (p gen T , m Q ) = σ(QQ + j; m Q , p gen T )/σ(QQ + j; m Q , p gen T = p * T ) σ(tt * + j; mt = m Q , p gen T )/σ(tt * + j; mt = m Q , p gen The structure of R Q (p gen T , m Q ) here is analogous to P Q (p j1 T,cut , m Q ) in Eq. (3.2), 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 [89][90][91], and overlay sgluon behavior on the same figure § . 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.  Figure 7: Normalized scale uncertainty bands of the pp → T pTp + j cross section as a function of minimum p j1 T (p j1 T,cut ), for (a) the inclusive pp → T pTp process at NLO+PS (dash) and at LO+PS as calculated with up to two additional matrix element-level jets via multijet merging (dot); as well as (b) the inclusive pp → T pTp process at NLO+PS (dot) and the pp → T pTp + j process (no multijet merging) at NLO+PS (dash).

Theoretical Uncertainties of the Monojet Process
As investigated in Sec. 3.3, were one to discovery 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 j1 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 j1 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 within QCD; see Sec. 2 for details. As the relevant signal topology in this work is characterized by the presence of hard QCD radiation recoiling against the T pTp system, our primary benchmark observable is the pp → T pTp + j cross section as a function of the transverse momentum of the process' leading, i.e., highest p T , jet (p j1 T ). We begin with Fig. 7a, where we show the pp → T pTp +j cross section as a function of p j1 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 Sec. 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 Tp = 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 a dependence on Q F by matching low-p 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 a smaller scale dependence at low p j1 T than in the merged LO+PS calculation. The small scale uncertainty observed for the lowest p j1 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. [92,93], lowest order-accurate observables, e.g., the p j1 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 uncertainty.
We assess this uncertainty 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 j1 T,cut . We assume multiplicative variations of the default parton shower scale, given in Eq. (2.4), 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 increasing (decreasing) for larger (smaller) values of Q S over the range of p j1 T considered. The shower scale variation amounts to absolute deviations from the fNLO result up to about +60% and −35%. For vanishing p j T , one should take caution in 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 TpTp → 0. 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 j1 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 j1 T m Tp ; for p j1 T m Tp , however, the T pTp + j uncertainty is about 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 [85]. We safeguard against the need for k T -resummation on the leading jet by using the (somewhat conservative) CSS consistency requirement from Ref. [88]. For m Tp = 600 GeV, and hence m TpTp > 1.2 TeV, one needs p j1 T,cut 200 − 250 GeV to force collinear logarithms from t-channel exchanges, which scale as δσ ∼ α s log(m TpTp /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 j1 T . The pp → T pTp process at NLO+PS, on the other hand, extends the support for p j1 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 j1 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. The curves are normalized to the Sherpa rate. For p j1 T,cut > 200 GeV, one sees that the Sherpa result, much like the MG5 aMC@NLO+PY8 result, possesses a stable uncertainty envelop spanning roughly ±18%. Moreover, the central value of Sherpa remains only 3 − 11% below MG5 aMC@NLO+PY8, and is consistent with the difference choice of α s (Q s ) employed by Sherpa. For p j1 T < 200 GeV, the MEPS@NLO calculation reveals that the Q F and Q R dependence remains smooth, and tapers down to ±16% as p j1 T,cut is relaxed. In comparison to the inclusive pp → T pTp at NLO+PS [GeV] Figure 9: Comparison of the pp → T pTp + j cross section estimated using MG5 aMC@NLO+PY8 and Sherpa with scale uncertainty band.
in Fig. 7a, the uncertainty band here is O(5%) larger in both directions. The slightly larger uncertainty here is consistent with the presence of additional soft O(α 2 s ) radiation terms whose scale dependence would otherwise be stabilized by O(α 2 s ) virtual terms. As described in Sec. 2, to construct the MEPS@NLO sample, various pp → T pTp + j NLO+PS computations at increasing p gen T are necessary to populate the p j T tail for p j T p gen T . To further demonstrate the comparability of the MG5 aMC@NLO+PY8 and Sherpa curves, we show in Fig. 9b the same quantities plotted in Fig. 9a for p j1 T,cut > 200 GeV but overlaid with the individual, unmerged pp → T pTp + j samples. For the range of p j1 T,cut considered, we find that the difference between the MG5 aMC@NLO+PY8 NLO+PS prediction and the Sherpa NLO+PS predictions for p j1 T,cut > p gen T differ about by only 0 − 15% again within the Q F , Q R uncertainty.
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 sub-leading 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 j1 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 deviations up to ±25%, this represents significant reduction in Q S dependence. For the T pTp +2j cross section, we see a larger scale dependence in both the 1j and 2j 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 moderateto-small factorization, 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. [72]. For the range of p j1 T,cut investigated, we find that the band is stable, symmetric, and spans roughly ±6% at low p j1 T,cut to ±10% at high p j1 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 j1 T,cut , we observe variations below ±5% for p j1 T,cut < 400 GeV and below 10% for p j1 T,cut < 1 TeV, and is comparable to our MC uncertainty.

Outlook
In Sec. 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 j1 T,cut ≈ 600 GeV at √ s = 14 TeV. The degeneracy can be resolved, however, through additional cross section measurements with larger p j1 T,cut and/or at a higher center-of-mass 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 break 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 j1 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 Sec. 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. 13(b)  In ratios and double ratios, such as those discussed in previous sections, the normalization component of uncertainties cancel. As shown in figure 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. (5.1). To verify this, we calculate the scale uncertainty for the ratio of the cross section to the nominal scale choice: .

(5.2)
For the normalizing cross section, we consider the m Tp = 600 GeV and vary the factorization and renormalization scales in 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. (5.2), 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 j1 T,cut dependence, one should consider different functional forms of Q R and Q F as well. In Fig. 14   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 Tp is about 15% higher than Q 0 = H T /2. The alternative scale choice Q R = Q F = m Tp is, however, also quite extreme in the sense that it is completely static [94]. 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 j1 T,cut = 1 TeV implies √ s > 2 TeV. Therefore, the ratio √ s/Q R is factor of 7, which leads to large non-asymmetric effects as seen in 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. 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 + 2j 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 MG5 aMC@NLO scale scheme of Q F = Q R = H T /2.) The 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  Figure 15: Cross section for the pp → T pTp process with additional light QCD partons at the ME level at leading, next-to-leading and next-to-next-to-leading orders in α s provided by the available MC tools.
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 [95][96][97][98].

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 (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-to-leading 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-to-no 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 dependencies 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.