Squark and gluino production cross sections in pp collisions at $\sqrt{s}$ = 13, 14, 33 and 100 TeV

TL;DR

This work delivers state-of-the-art predictions for squark and gluino production cross sections at $pp$ colliders with $\sqrt{s}=13$, 14, 33, and 100 TeV by combining next-to-leading order SUSY-QCD corrections with next-to-leading-logarithmic soft-gluon resummation (NLO+NLL) and matching to preserve non-soft contributions. The results, computed with the $\text{NLL-fast}$ framework and cross-validated against $\text{CTEQ6.6}$ and $\text{MSTW2008}$ PDFs, provide central values and robust uncertainty bands (scale and PDF/αs) following PDF4LHC prescriptions for a wide range of SUSY masses and production channels, including gluino-gluino, squark-antisquark, squark-gluino, and direct stop/sbottom production. The paper also discusses special cases such as decoupled squarks or gluinos and simplified models, and makes the full set of results and interpolation tools publicly available to experimental collaborations for interpreting current and future searches. The findings underscore how cross sections grow with collider energy and how PDFs at large $x$ influence uncertainties, delivering essential benchmarks for high-mass SUSY reach at the LHC Run II/III and future hadron colliders. The work thus provides a practical, theoretically solid reference for SUSY cross sections and their uncertainties across a broad energy landscape.

Abstract

We present state-of-the-art cross section predictions for the production of supersymmetric squarks and gluinos at the upcoming LHC run with a centre-of-mass energy of $\sqrt{s} = 13$ and $14$ TeV, and at potential future $pp$ colliders operating at $\sqrt{s} = 33$ and $100$ TeV. The results are based on calculations which include the resummation of soft-gluon emission at next-to-leading logarithmic accuracy, matched to next-to-leading order supersymmetric QCD corrections. Furthermore, we provide an estimate of the theoretical uncertainty due to the variation of the renormalisation and factorisation scales and the parton distribution functions.

Figures (16)

• Figure 1: NLO+NLL production cross sections for the case of equal degenerate squark and gluino masses as a function of mass at $\sqrt{s}={13}$ TeV.
• Figure 2: NLO+NLL production cross sections for the case of equal degenerate squark and gluino masses as a function of mass at $\sqrt{s}={14}$ TeV.
• Figure 3: NLO+NLL production cross sections for the case of equal degenerate squark and gluino masses as a function of mass at $\sqrt{s}={33}$ TeV.
• Figure 4: NLO+NLL production cross sections for the case of equal degenerate squark and gluino masses as a function of mass at $\sqrt{s}={100}$ TeV.
• Figure 5: NLO+NLL gluino pair production cross section with squarks decoupled as a function of mass at $\sqrt{s}={13}$ TeV in the wider (upper plot) and narrower (lower plot) mass range. The different styled black (red) lines correspond to the cross section and scale uncertainties predicted using the CTEQ6.6 ( MSTW2008) PDF set. The yellow (dashed black) band corresponds to the total CTEQ6.6 ( MSTW2008) uncertainty, as described in the text. The green lines show the final cross section and its total uncertainty.