Average Event Properties from LHC to FCC-hh

In the context of design studies for future $pp$ colliders, we present a set of predictions for average soft-QCD event properties for $pp$ collisions at $E_\mathrm{CM} = 14$, $27$, and $100$ TeV. The current default Monash 2013 tune of the PYTHIA 8.2 event generator is used as the baseline for the extrapolations, with uncertainties evaluated via variations of cross-section parametrisations, PDFs, MPI energy-scaling parameters, and colour-reconnection modelling, subject to current LHC constraints. The observables included in the study are total and inelastic cross sections, inelastic average energy and track densities per unit pseudorapidity (inside $|\eta|\le 6$), average track $p_\perp$, and jet cross sections for 50- and 100-GeV anti-$k_T$ jets with $\Delta R=0.4$, using aMC@NLO in conjunction with PYTHIA 8 for the latter.


Motivation
The long-term future of physics at the energy frontier is looking increasingly towards designs of pp colliders operating at CM energies of up to 100 TeV (FCC-hh or SppC) [1][2][3][4][5][6]. In the context of detector design efforts and basic phenomenology studies for such machines, as well as for the 14-TeV LHC and 27-TeV HE-LHC options, a set of predictions (or, at the very least, extrapolations) are required for elementary quantities characterising event properties such as overall rates (total and inelastic cross sections), tracker occupancies (charged-track densities) and total energy depositions (dE/dη). Going beyond predictions for single events, superpositions of minimumbias events are also used to estimate the dependence of the overall activity on the number of events per bunch crossing (pile-up), which is projected to reach O(1000) for peak luminosities of 3 × 10 35 cm −2 s −1 targeted at both HE-LHC and FCC-hh [6,7].
The above are all observables which are dominated by soft QCD effects, and which can therefore not be computed perturbatively. Instead, one is forced to rely either on parametric fits or on explicit physics models such as those implemented in (soft-inclusive) Monte Carlo event generators. In the context of the latter, we here consider the modelling offered by the Pythia 8 event generator [8], and the degree to which its extrapolations are affected by several known sources of theoretical uncertainties. For definiteness, we focus on predictions for an instrumented region corresponding to |η| ≤ 6, illustrated in fig. 1.
After verifying that the default Monash 2013 tune [9] of the Pythia 8.2 event generator still gives an acceptable description of soft-inclusive observables at the LHC including new measurements at 13 TeV, we consider its extrapolations to higher CM energies and perform several salient variations, e.g., of total cross sections [10], colour reconnections, PDF sets, and multi-parton interactions. This study updates and extends the Snowmass white paper in [11] which was based on Pythia 6.4 [12]. It is complementary to the that of [13] which considered extrapolations of the default predictions (i.e., without parameter variations) of several qualitatively different MC models of soft-inclusive QCD reactions, including Epos [14], Phojet [15], and Qgsjet [16]. We verify that we obtain consistent results for the common reference model considered in both studies (default Pythia 8.2), and note that most of the models considered in [13] exhibit a rather similar scaling behaviour of, e.g., the central charged-particle densities over the extrapolated region, predicting that it should grow by about a factor of 2 from 10 to 100 TeV. The exception is the significantly slower scaling exhibited by Phojet 1.12 which is however already in strong conflict with the LHC measurements, hence we do not consider it a realistic variation 1 .
We are thus reasonably confident that the Monash tune provides an acceptable starting point for extrapolations up to 100 TeV, and that useful uncertainty estimates can be obtained by variations of its salient parameters. As a desirable side effect, Pythia's universal modelling of both hard and soft processes then also allows the same models and variations to be (re)used in the context of studies of hard high-p ⊥ processes as well. In the context of such hard (perturbative) processes at 100 TeV, the reader may also be interested in the dedicated study performed in [18].
The report is organised as follows. In section 2, we briefly recapitulate the physics of soft QCD processes in the context of Pythia 8, with emphasis on inclusive soft-QCD cross sections and the modelling of MPI and colour reconnections. In section 3, we summarise the variations that we have performed for uncertainty estimates, along with the LHC constraints imposed to limit the variations. The main set of results for FCC-hh are presented in section 4, with a brief "executive summary" of our main conclusions provided in section 5.

Overview of soft QCD in Pythia 8
We shall start by providing a brief review of the modelling of soft QCD processes in Pythia 8 [8], placing particular emphasis on those areas relevant to the (re)tuning performed in section 3, in which variants of the default tune settings are determined. These variants will be used to provide uncertainty ranges on predictions extrapolated to FCC-hh energies in section 4.

Modelling of the total, elastic and diffractive cross sections
The total cross section is given by the sum of elastic (el) and inelastic cross sections; the latter is further broken down into the sum of non-diffractive (ND), single diffractive (SD) and double diffractive (DD) cross sections, so that the total cross section is expressed as a sum of four components, σ tot (s) = σ el (s) + σ ND (s) + σ SD (s) + σ DD (s) . (1) The distinction between the three components of the inelastic cross section is in principle slightly model dependent (related to exactly how one distinguishes high-mass diffractive excitations from non-diffractive events), but this is beyond the level of accuracy required for this study. We note also that Pythia allows for a small central-diffractive component to be added as well, typically at the level of O(1%) of the total cross section; however this option is switched off by default and has not so far been included in any tunes. Since it is primarily relevant for dedicated central-exclusive studies, beyond the level of detail aimed for in this study, it is neglected in the following discussion. In Pythia 8.2 by default the total cross section is modelled by the Donnachie-Landshoff (DL) parameterisation [19]: where ε = 0.0808, η = 0.4525, A and B denote the initial state particles and X AB , Y AB depend upon these states. The default elastic and diffractive cross sections are described by the Schuler-Sjöstrand (SaS) framework [20]. This is a diagrammatical approach based on the exchange of pomerons and reggeons and application of the optical theorem, similar to the Ingelman-Schlein approach [21]. The non-diffractive cross section, which is calculated from the total, elastic, and diffractive ones by solving for σ ND in eq. (1), not only determines the total rate of events that will be modelled as non-diffractive (typically starting from a partonic 2 → 2 QCD interaction process), but is also used for determining the expected number of multi-parton interactions, which we discuss in greater detail in section 2.2. The DL and SaS parameterisations worked well up to Tevatron energies; however the diffractive cross section grows somewhat too quickly with energy and already overshoots the data at LHC energies [22]. For this reason, some alternatives to these models have now been implemented in Pythia 8. First there is the Minimum Bias Rockefeller (MBR) implementation [23] of a model based on Regge theory [24,25]; this primarily addresses diffraction but also provides parameterisations of the total and elastic cross sections. Recently [10], two further options for the total and elastic cross section were implemented: the ABMST model [26] and a parameterisation by the COMPAS group found in the Review of Particle Physics 2016 (RPP) [27] 2 . The former is an extension of DL, allowing four types of single exchanges, and all combinations thereof for double exchanges; the RPP model is further complicated still, allowing six different single exchanges and certain combinations of double exchanges.
Additionally, the ABMST model also addresses single diffraction, although it was observed that this also predicts a diffractive cross section which grows faster than the total cross section with energy. Therefore some modifications to this model giving a more reasonable high energy behaviour were made in the context of the implementation in Pythia 8 [10].

Modelling of multi-parton interactions
In Pythia 8 the modelling of multi-parton interactions (MPI) [28][29][30] represents an effort to extend perturbative QCD down to low transverse momenta. Perturbatively, the hard scattering cross section above some cutoff p ⊥min is given by with the differential cross section obtained from the factorisation formula: dσ dp 2 The leading partonic differential cross sections are t-channel exchanges which behave as 1/p 4 ⊥ in the p ⊥ → 0 limit, so σ hard is divergent for p ⊥min → 0 and -at current LHC energies -exceeds the total cross section for p ⊥min ∼ 5 GeV (see e.g. [9, fig. 34]). This apparent contradiction is reconciled by recognising that σ hard represents the cross section for not one but any number of parton-parton interactions, of which there may be several in a single hadron-hadron collision.
In the limit that the pairwise interactions are independent, the number of parton-parton scatterings per hadron-hadron collision would be Poisson distributed, with an average given by the ratio of the hard to non-diffractive (ND) cross section 3 : Within Pythia, this expression is interpreted as the kernel of an ordered sequence of partonparton interactions [28], such that the interaction with the highest x ⊥ = p ⊥ / √ s is generated first. This allows to build in modifications to the PDFs for subsequent MPIs to account for correlations due to momentum and flavour conservation and to interleave the generation of MPI with the shower evolution in a single common sequence of ordered p ⊥ values [31]. All of these aspects generate deviations from the simple Poissonian ansatz, with momentum conservation being the most important.
While eq. (5) resolves the issue of why σ hard may exceed the total cross section, the former still exhibits a divergence as p ⊥min → 0. This is a problem because the average scattering energy does not decrease with p ⊥min as quickly as the average number of scatterings rises. The divergence in σ hard therefore naively implies that the total scattering energy becomes infinite. However, use of the partonic cross-section assumes free incoming partons, which is not valid for such low transverse momenta where partons become strongly bound. Instead one can employ a model of effective colour-screening, where colour charges cannot be resolved over distances d ∼ 1/p ⊥0 . The partonic cross section is thus regularised according to: and the effective screening scale p ⊥0 is given an energy dependence: The arbitrary reference energy E ref CM is typically set to 7 TeV in modern tunes, while the parameters p ref ⊥0 and E pow CM should be tuned. Essentially, p ref ⊥0 still sets the average number of MPI per (inelastic, non-diffractive) hadron-hadron collision at the given reference CM energy (for given σ ND and σ hard parameters), while E pow CM controls the scaling of that number to other CM energies, with lower values of E pow CM producing faster scalings of n . The number of scatterings should also reflect the fact that the incoming hadrons are extended objects, with more peripheral collisions expected to have fewer scatterings on average. This can be quantified through the impact parameter b; although the matter profile of the hadrons is not known a priori, one option 4 is to parameterise the time-integrated overlap function directly where the power p is a free parameter to be tuned, with e.g. p = 2 corresponding to a Gaussian density profile. Empirically, values between 1.5 and 2 are preferred. The average number of interactions (per event having at least one scattering) at impact parameter b is given by and k is a constant to be found by integrating over b and comparing to eq. (5). This may be used to find a probability distribution dependent on both x ⊥ and b. Thus, when the first (hardest) MPI is generated, an overall b is selected as well, and the corresponding overlap factor is used to enhance (for smaller-than-average b) or suppress (for larger-than-average b) the probability for all subsequent MPI as well.

Modelling of colour reconnections
The Lund model of hadronisation implemented in Pythia interprets the linear behaviour of the quark-antiquark potential at long distances as a string which can fragment into hadrons through quantum mechanical tunnelling; see e.g. [32]. In particular, for a given event the strings are identified as occurring between colour-connected partons, with quarks (and antiquarks) as string endpoints and gluons appearing as transverse "kinks" on these strings. The picture of colour connections is particularly simple in the leading colour approximation used in nearly all current parton showers 5 , where each successive emission generated is colour connected to its parent emitters. Models of colour reconnection allow for strings to form between partons beyond these simple leading colour topologies effectively allowing for either interference between different flows ('static' colour reconnection) or soft gluon exchanges ('dynamic' colour reconnection). Such effects were already seen to play a role at LEP in measurements of the W boson mass [39]; they become even more important at the LHC due to the existence of vastly more possibilities for coloured initial state partons, beam remnants and multi-parton interactions, all allowing for non-trivial (reconnections of) coloured topologies. The leading assumption is that colour reconnection models should give rise to a greater number of shorter strings, thereby affecting the number and energy density distributions of final state particles. There are by now many indications that such effects are essential for modelling soft QCD effects in general, and there are important open questions concerning their impact on precision observables such as the hadronic top quark mass [40,41].
In Pythia 8 several models for colour reconnections have been developed and implemented over the years, however it remains one of the least well understood phenomena in hadron collisions, and in particular we know little about how such effects scale with energy [42]. Therefore in assessing uncertainties on predictions for future colliders, it is important not just to include variations on the baseline model for colour reconnections, but also consider some qualitatively different models which may exhibit a different scaling.
In the oldest model (used by default), gluons from MPI with low transverse momentum scales may be successively inserted into the colour flow of a higher transverse momentum MPI, in such that the total string length is minimised by brute force. A more recent model has also been implemented [43], which combines the earlier string-length minimisation arguments with selection rules based on the colour algebra of SU(3). Another option known as the 'gluon-move' model [41] is also available, which we mention here for completeness; however in the following section we only consider tunes using the former two. In this model, any pair of gluons can be reconnected, irrespective of whether they were produced in the same or different MPI. A string-length-reduction measure is calculated, and pairs of gluons maximising this measure (or minimising change in string length) are sequentially reconnected until the incremental reduction in string-length is below some threshold.

Variants of Monash 201Tune and Validation
In the next section, we present predictions for FCC-hh of certain distributions sensitive to soft QCD, using the Monash 2013 tune of Pythia 8.2 as the baseline. Despite being a few years old by now, this tune remains a useful benchmark; as we show below, it remains in acceptable agreement with soft-inclusive LHC measurements up to 13 TeV, and as the default tune for Pythia 8.2 there is a large body of existing validations and further complementary studies can be done without requiring any special parameter settings. Asides from being a convenient choice, we also estimate that its predictions for future energies are likely to be a bit on the conservative side, since already at 13 TeV it appears to slightly overestimate the number density of charged particles (see section section 3.1 and fig. 4).
Since there are substantial uncertainties arising from ambiguities in the modelling of nonperturbative physics, however, we do not restrict ourselves to a single tune. Instead, in this section we explore variants of the baseline tune in which the parameters associated with the modelling aspects discussed in the previous section (modelling of PDFs, MPI, total, elastic and diffractive cross sections, and colour reconnections), are modified, one by one. We shall then use the baseline tune together with these variants to define an uncertainty envelope for the extrapolations to higher CM energies in section 4.

PDF variation
We recall from section 2.2 that the probability distribution for multiple scatterings is dependent on dσ/dp 2 ⊥ and hence on the PDFs through eq. (4). MPI are primarily driven by the low-x behaviour of the gluon PDF, which is not tightly constrained and can vary significantly between different sets. In addition to the baseline NNPDF 2.3 QCD+QED LO set (with α s (m Z ) = 0.13), the following three LO PDF sets (interfaced via LHAPDF 6 [44]) were chosen as spanning a reasonable range of qualitatively different behaviours: • NNPDF31 LO as 0130 [45] • HERAPDF15LO EIG [46] • MMHT2014lo68cl [47] In fig. 2, we show a comparison of the gluon distributions of these sets at a relatively low scale, Q = 2 GeV, representative of a "typical" scale probed by MPI. Note that the CT14lo [48] set is also shown 6 but was not found to add qualitatively to the range already spanned by the other sets.   . Plots generated using the APFEL web interface [49,50]).
To define Monash tune variants using the alternative PDF sets, we re-tune a few of the parameters which control the amount of MPI to reobtain a basic level of agreement with LHC measurements; see below. At the very least, we consider the p ⊥0 parameter discussed in section 2.2 to be intimately tied to the choice of PDF and hence obligatory to re-tune when the PDF is modified. For slightly more flexibility, we allow the following three parameters to be re-tuned when changing the PDF set:  7) and p in eq. (8) respectively. For the tuning, a variety of analyses measuring the underlying event and minimum bias at the LHC at 7 TeV, implemented in Rivet [55] were considered. The value of Multiparton-Interactions:pT0Ref was tuned such that the number of charged tracks in η, dN ch /dη was compatible with the level measured, for example, by CMS [51], as shown in fig. 3a. The energy flow at high η [52] is also shown, in fig. 3b, since that is one of the observables we will be considering in the extrapolations to 100 TeV. The value for MultipartonInteractions:expPow was tuned to the density of charged tracks as a function of the transverse momentum of the hardest track, for example, as measured by ATLAS [53], cf. fig. 3c. Although shape differences arise in this observable, these are no more significant than for the default Monash tune. Finally, we also include a comparison to the average track p ⊥ as function of the track multiplicity, in fig. 3d, again as a reference for one of the observables we shall consider in our extrapolations to higher energies.
We briefly note that for the setup in fig. 3d, namely minimum bias events having at least one charged track with transverse momentum exceeding 500 MeV, the tunes all undershoot the data by about 5%. It is known that the transverse momentum of charge tracks is imperfectly modelled: see for example figure 18 in [9], where the ratio to data of the differential multiplicity with respect to transverse momentum is not flat. This is likely due to an incorrect distribution in the transverse momentum kicks that result from string-breaking during hadronisation (currently modelled by a Gaussian). The result is that a slightly different selection at trigger level can result in a different normalisation relative to data in the average transverse momentum (see for example figure 20 in [9]). To improve upon this situation will require improvements in the modelling of hadronisation, which we do not discuss further here.
MultipartonInteractions:ecmPow controls the energy scaling of the p ⊥0 regularisation scale for MPI, as discussed in section 2.2. This was tuned by considering a particularly sensitive observable, namely the relative increase in the density of tracks between 0.9-2.36 TeV and 0.9-7 TeV measured by ALICE [56]. We here extend this plot to include recent 13 TeV data [57], though without the benefit of internal ALICE systematics studies, we cannot take correlations between the 900-GeV and 13-TeV uncertainties into account and merely combine them using simple quadratures, yielding a conservative overestimate which is fine for our purpose. While it is quite easy to tune MultipartonInteractions:ecmPow to capture the relative increase between from 0.9 to 2.36 and 7 TeV, such tunes typically exhibit some tension with the 13 TeV data. Therefore, for each of the changed PDF sets we determine two alternatives -a minimum and maximum value -for this parameter. The results are shown in fig. 4a, with two points shown at each energy for each of the HERAPDF, MMHT, and NNPDF 3.1 variants; the lower (higher) point corresponding to the minimum (maximum) value for MultipartonInteractions:ecmPow.
For completeness, the plot in fig. 4a also shows the predictions for PHOJET [58], as these were also originally included in [56]. These predictions clearly do not capture the rise in the  Figure 3: A selection of distributions measuring the underlying event level at the LHC at 7 TeV which were used in the tuning of variants of Monash 2013 with different PDF selection, described in the text. The differential number of charged tracks in pseudorapidity as measured by CMS [51] is shown in a. The differential energy deposit at large rapidities as measured by CMS [52] is shown in b. The density of charged tracks as a function of the transverse momentum of the hardest track as measured by [53] is shown in c. The average transverse momentum of charged tracks as a function of the number of charged tracks measured by ATLAS [54] is shown in d. Plots were produced using Rivet [55]. Central Charged Track Multiplicity Relative Increase 0.9-2.36 TeV ALICE data 0.9-7.0 TeV ALICE data 0.9-13.0 TeV ALICE data (b) Cross-Section and CR Variations Figure 4: The relative increase in the number density of charged tracks between 0.9-2.36 TeV and 0.9-7 TeV measured by ALICE [56] and updated to include 13 TeV data [57]. This observable was used to tune to MultipartonInteractions:ecmPow for variants of the Monash 2013 tune with a alternative choice for the PDF, and b different choices for the modelling of the total cross section and colour reconnections. For comparison, the predictions for PHOJET [58] originally included in [56] are also shown in a. In all cases, where two points are shown, these reflect two choices for the parameter MultipartonInteractions:ecmPow, with the lower point corresponding to the minimum value and the upper point corresponding to the maximum value.
charged track density, a trend which we expect to be exacerbated when going to 100 TeV. We therefore discount PHOJET as a realistic variation.

Variation of cross-section parameterisations
As discussed in section 2.1 the default parameterisations of the total, elastic and diffractive cross sections in Pythia 8 (namely the DL/SaS parameterisations) are not able to faithfully reproduce all aspects of the energy scalings seen in measurements, and a set of alternatives were recently implemented [10]. For the total cross section we consider the MBR, ABMST and RPP models (selected by modifying SigmaTotal:mode). In the case of the former two, the same models also offer alternative modelling of the diffractive cross section (selected by modifying  SigmaDiffractive:mode). Note that in the case of the choice of ABMST for modelling the diffractive cross section, the modifications proposed in [10] are selected (by using the default value for SigmaDiffractive:ABMSTmodeSD).
We only found it necessary to retune one parameter, the energy-scaling parameter Multipar-tonInteractions:ecmPow, for the above alternative parameterisations of the total/elastic/diffractive cross sections. As in the previous section it was found that the 13 TeV data highlight tensions in the predictions, and accordingly we selected two distinct values instead of a single "average" value, as shown in table 2. The effect of these variations upon the relative increase in the density of charged tracks is shown in fig. 4b.

Variation of modelling of colour reconnections
To gauge the uncertainty arising from the still poorly understood issue of colour reconnections, both the default model (based on string minimisation alone) and a more recent model (based on QCD selection rules) were considered, see the discussion in section 2.3. For the latter, the modifications to tuned parameters proposed in [43] were found to be sufficient. In addition, and as in previous sections, two values for MultipartonInteractions:ecmPow were determined, the effect of which is shown in fig. 4b. The full set of changed parameters for these variants is shown in

Predictions for FCC-hh
In this section we present predictions for FCC-hh obtained using the default Monash 2013 tune of Pythia 8 and the variants of this tune described in section 3.
In fig. 5, we show extrapolations of the total and inelastic cross sections predicted by the selection of parameterisations described in section 3.2, compared to measurements by TOTEM and ALICE [59][60][61][62] at LHC energies. (See appendix A for tabulated values.) The default parameterisations used in Monash 2013 are in good agreement with the measurements of the total inelastic cross section but are in conflict with the total cross section; this is due to the elastic cross section being too small (by about 10mb). The alternative models were selected for their good agreement with both the total and inelastic cross sections. We note also that the largest inelastic cross section at 100 TeV (corresponding to the highest inelastic event rate), among the models considered here, is predicted by the default parameterisation of the Monash 2013 tune, followed by the RPP/SaS, ABMST, and MBR parameterisations, respectively.
We now turn to the observables characterising the properties of average inelastic events. Below, in figs. 6 -9, we show predictions for the following quantities: • The total number of charged tracks inside |η| = 6.0, shown in fig. 6a.
• The average transverse momentum of charged tracks, shown in fig. 6b. Also shown, in fig. 7, is the differential average transverse momentum of charged tracks as a function of absolute pseudorapidity.
• The charged track density in slices of pseudorapidity, fig. 8. • The energy density of all final state particles in slices of pseudorapidity, fig. 9.
The events counted in distributions of charged tracks (unless specified otherwise) were required to have at least one charged track inside |η| = 6.5. Furthermore events were generated such that only those particles having proper lifetimes cτ 0 < 10 mm are decayed (by setting: Par-ticleDecays:limitTau0 = On and ParticleDecays:tau0Max = 10.0). In the distribution of energy density, all final state particles excluding neutrinos and muons were counted. As a cross check, the fully inclusive distributions were also calculated, but the inclusion of neutrinos and muons only give rise to per-mille level differences. We note that our benchmark predictions for the Monash tune of Pythia are consistent 7 with the earlier study of [13]. We emphasise that although [13] considered a similar set of observables to ours, they were mainly concerned with the differences between a range of qualitatively different generators, while our study primarily addresses how the modelling uncertainties for Pythia scale with energy.
So far, we have focused on the modelling of nonperturbative QCD effects as the dominant source of uncertainty on observables sensitive to the low-scale physics which dominates the bulk of inelastic events. By contrast, uncertainties for hard infrared-safe observables are mainly perturbative in origin, and are usually typified by performing variations of the renormalisation and factorisation scales. As an example, we here provide predictions for the inclusive jet cross sections at next-to-leading order in QCD for p ⊥j ≥ 50 GeV and p ⊥j ≥ 100 GeV using MG5 -aMC@NLO [63] matched to Pythia. Jets are defined using the anti-k T jet-finding algorithm with R = 0.4. Scale variations are performed by varying µ R and µ F independently by between factors as 2.0 and 0.5 8 . The central scale choice used is the default in MG5 aMC@NLO, namely half the sum of transverse masses over all particles m T 2 (equivalent to H T /2 for massless particles). In addition we provide predictions for two choices of PDF set, namely NNPDF 2.3 (which 7 In [13] there were some inconsistencies in the setup for the predictions of pT between that described in the main text and that shown on the plot; we find agreement if we use the setup shown on the plot. 8 Scale variations were performed using the conventional method of 7-point variation. is the default both for MG5 aMC@NLO and the Monash 2013 tune) and HERAPDF 1.5; these were selected as having the largest shape differences over a range of x values. Note that in MG5 aMC@NLO we use the NLO set for these PDFs, while in the shower and generation of MPI in Pythia we use the LO set so as not to break the underlying-event tuning 9 .
The results are shown in fig. 10, for fully hadronised final-state particles. (See appendix A for tabulated values.) For the modelling of non-perturbative effects we select the tune appropriate to the LO PDFs 10 , with the default settings for NNPDF 2.3, and the modifications shown in the fourth column of table 1 for HERAPDF 1.5. In the results for HERAPDF we considered both the upper and lower variation for the parameter MultipartonInteractions:ecmPow; the difference was found to be at most 1% for p T > 50 GeV at 100 TeV and otherwise considerably less. As this is notably less than the size of scale uncertainties, we therefore show only the results for the lower variation. The difference between the two PDF sets is somewhat larger, but nevertheless much smaller than scale uncertainties.
As an explicit check that the inclusive cross sections are indeed relatively insensitive to non-perturbative uncertainties, we note that the extreme variation obtained by switching off hadronisation and colour-reconnection in Pythia only modifies these cross sections by approximately 3-4% at √ s = 14 TeV and by 1-3% at √ s = 100 TeV. We note that the earlier study of [18] also made predictions for inclusive jet cross sections for a 100 TeV pp collider; however, this focused primarily on the scaling of jet rates and jet substructure (and did not include NLO predictions for the pure QCD single inclusive jet cross section). In addition, [18] considered several different hard processes and their impact for BSM searches.

Summary and Conclusions
We have presented predictions for the expectation values of some basic observables characterising inelastic event properties at future hadron colliders, focusing on CM energies of 14, 27, and 100 TeV. The predictions are based on extrapolations made with the Monash 2013 tune of the Pythia 8.2 event generator as the baseline and uncertainties were estimated by variations of its PDFs, cross-section parametrisations, colour-reconnection modelling, and MPI p ⊥0 regularisation scale, subject to constraints imposed by current LHC measurements at √ s = 0.9-13 TeV.
The total proton-proton cross section is expected to grow from about 108 mb at √ s = 14 TeV to about 145 mb at √ s = 100 TeV, though we note that the 1992 Donnachie-Landshoff fit used by the baseline Monash tune underestimates these numbers by about 5% mainly due to its elastic cross section being too small and already in conflict with LHC data. We therefore do not advise to use the Monash 2013 extrapolations for total and elastic cross sections. The inelastic proton-proton cross section is expected to rise from 77-80 mb at √ s = 14 TeV to 83-88 mb at √ s = 27 TeV, to 97-108 mb at √ s = 100 TeV, with the Monash 2013 predictions being the highest among the parametrisations considered in this study and hence relatively conservative within the context of this study. The total number of charged particles inside |η| < 6 (with cτ 0 ≥ 10 mm) per inelastic event grows by slightly less than a factor of 2 -from about 70 at √ s = 14 TeV to about 120 at √ s = 100 TeV, with estimated uncertainties of ∼ 10%. The average p ⊥ of those particles also increases, albeit only slightly, starting from p ⊥ ∼ 460 MeV at √ s = 14 TeV and increasing by 50 -100 MeV when extrapolating to √ s = 100 TeV.
In the central region of the detector, the amount of energy deposited per inelastic event grows by about a factor of 2 -from ∼ 7 GeV per unit η at √ s = 14 TeV to ∼ 15 GeV per unit η at √ s = 100 TeV. At high rapidities, much closer to the beam, the total amount of energy deposited is of course much larger, and it is also predicted to scale faster. At |η| = 6, we estimate about 600 GeV of total energy per unit rapidity at 14 TeV, while we predict ∼ 1700 GeV per unit rapidity at 100 TeV. For the highest rapidities, it is also worth remarking that our choice of reference model, the Monash 2013 tune, lies towards the upper limit of the range spanned by our uncertainty estimates, cf. fig. 9, hence its predictions of the maximum energy densities that would be faced by detectors in this region can be considered relatively conservative at least within the context of the variations we have studied.
The cross section for anti-k T jets with ∆R = 0.4 and p ⊥ ≥ 50 GeV is predicted to increase faster than the relative increase in CM energy, by a factor of ∼ 13 from √ s = 14 TeV to √ s = 100 TeV, while the cross section for a cut of p ⊥ ≥ 100 GeV is expected to increase by about a factor of 21 over the same range.
We hope this investigation, and the LHC-vetted Pythia parameter settings that we have developed in the course of it, will prove useful to the exploration of phenomenology and detector design concepts for future hadron colliders.       Only particles with proper lifetimes cτ 0 < 10 mm are considered to decay. Uncertainties are shown as percentages and are found as described in the text.