Measurements of proton and nuclear collisions at the Large Hadron Collider at nucleon-nucleon c.m. energies up to $\sqrt{s_{NN}}=$ 13 TeV, have improved our understanding of hadronic interactions at the highest energies reached in collisions of cosmic rays with nuclei in the earth atmosphere, up to $\sqrt{s_{NN}}\approx 450$ TeV. The Monte Carlo event generators (EPOS, QGSJET, and SIBYLL) commonly used to describe the air showers generated by ultrahigh-energy cosmic rays (UHECR, with $E_{CR}\approx 10^{17}$--$10^{20}$ eV) feature now, after parameter retuning based on LHC Run-I data, more consistent predictions on the nature of the cosmic rays at the tail of the measured spectrum. However, anomalies persist in the data that cannot be accommodated by the models. Among others, the total number of muons (as well as their maximum production depth) remains significantly underestimated (overestimated) by all models. Comparisons of EPOS, QGSJET, and SIBYLL predictions to the latest LHC data, and to collider MC generators such as PYTHIA, indicate that improved description of hard multiple minijet production and nuclear effects may help reduce part of the data--model discrepancies, shed light on the UHECR composition approaching the observed $E_{CR}\approx 10^{20}$ eV cutoff, and uncover any potential new physics responsible of the observed anomalies.

The PYTHIA 6 Monte Carlo (MC) event generator, commonly used in collider physics, is interfaced for the first time with a fast transport simulation of a hydrogen atmosphere, with the same density as air, in order to study the properties of extended atmospheric showers (EAS) produced by cosmic ray protons with energies E$_{CR}\approx 10^{14}$--$10^{20}$ eV. At variance with the hadronic MC generators (EPOS-LHC, QGSJET, and SIBYLL) commonly used in cosmic-rays physics, PYTHIA includes the generation of harder hadronic jets and heavy (charm and bottom) quarks, thereby producing higher transverse momentum final particles, that could explain several anomalies observed in the data. The electromagnetic, hadronic, and muonic properties of EAS generated with various settings of PYTHIA 6, tuned to proton-proton data measured at the LHC, are compared to those from EPOS-LHC, QGSJET 01, QGSJET II, and SIBYLL 2.1. Despite their different underlying parton dynamics, the characteristics of the EAS generated with PYTHIA 6 are in between those predicted by the rest of MC generators. The only exceptions are the muonic components at large transverse distances from the shower axis, where PYTHIA predicts more activity than the rest of the models. Heavy-quark production, as implemented in this study for a hydrogen atmosphere, does not seem to play a key role in the EAS muon properties, pointing to nuclear effects as responsible of the muon anomalies observed in the air-shower data.

The framework to compute the cross sections for the production of particles with high mass and/or large transverse momentum in double- (DPS), triple- (TPS), and in general $n$-parton scatterings, from the corresponding single-parton ($\sigma_{\rm SPS}$) values in high-energy proton-proton, proton-nucleus, and nucleus-nucleus is reviewed. The basic parameter of the factorized $n$-parton scattering ansatz is an effective cross section $\sigma_{\rm eff}$ encoding all unknowns about the underlying generalized $n$-parton distribution in the proton (nucleon). In its simplest and most economical form, the $\sigma_{\rm eff}$ parameter can be derived from the transverse parton profile of the colliding protons and/or nucleus, using a Glauber approach. Numerical examples for the cross sections and yields expected for the concurrent DPS or TPS production of heavy-quarks, quarkonia, and/or gauge bosons in proton and nuclear collisions at LHC and Future Circular Collider (FCC) energies are provided. The obtained cross sections are based on perturbative QCD predictions for $\sigma_{\rm SPS}$ at next-to-leading-order (NLO) or next-to-NLO (NNLO) accuracy including, when needed, nuclear modifications of the corresponding parton densities.

A generic expression to compute triple parton scattering (TPS) cross sections in high-energy proton-nucleus (pA) collisions is derived as a function of the corresponding single-parton cross sections and an effective parameter encoding the transverse parton profile of the proton. The TPS cross sections are enhanced by a factor of $9\,A\approx 2000$ in pPb compared to those in proton-nucleon collisions at the same center-of-mass energy. Estimates for triple charm ($c\overline{c}$) and bottom ($b\overline{b}$) production in pPb collisions at LHC and FCC energies are presented based on next-to-next-to-leading order (NNLO) calculations for $c\overline{c}, b\overline{b}$ single-parton cross sections. At $\sqrt{s_{NN}} = 8.8$ TeV, about 10% of the pPb events have three $c\overline{c}$ pairs produced in separate partonic interactions. At $\sqrt{s_{NN}} = 63$ TeV, the pPb cross sections for triple-J$/\psi$ and triple-$b\overline{b}$ are ${\cal O}$(1--10 mb). In the most energetic cosmic-ray collisions observed on earth, TPS $c\overline{c}$-pair cross sections equal the total p-Air inelastic cross section.

The global properties of the final states produced in hadronic interactions of protons at centre-of-mass energies of future hadron colliders (such as FCC-hh at CERN, and SppC in China), are studied. The predictions of various Monte Carlo (MC) event generators used in collider physics (PYTHIA 6, PYTHIA 8, and PHOJET) and in ultrahigh-energy cosmic-rays studies (EPOS, and QGSJET) are compared. Despite their different underlying modeling of hadronic interactions, their predictions for proton-proton (p-p) collisions at $\sqrt{s}$ = 100 TeV are quite similar. The average of all MC predictions (except PHOJET) for the different observables are: (i) p-p inelastic cross sections $\sigma_{\rm inel}$ = 105 $\pm$ 2 mb; (ii) total charged multiplicity $\rm N_{_{\rm ch}}$ = 150 $\pm$ 20; (iii) charged particle pseudorapidity density at midrapidity $\rm dN_{ch}/d\eta|_{\eta=0} = 9.6 \pm 0.2$; (iv) energy density at midrapidity $\rm dE/d\eta|_{\eta=0} = 13.6 \pm 1.5$ GeV, and $\rm dE/d\eta|_{\eta=5} = 670 \pm 70$ GeV at the edge of the central region; and (v) average transverse momenta at midrapidities $\rm < p_{T}> = 0.76 \pm 0.07$ GeV/c. At midrapidity, EPOS and QGSJET-II predict larger per-event multiplicity probabilities at very low ($\rm N_{\rm ch}<3$) and very high ($\rm N_{\rm ch}>100$) particle multiplicities, whereas PYTHIA 6 and 8 feature higher yields in the intermediate region $\rm N_{\rm ch}\approx$ 30--80. These results provide useful information for the estimation of the detector occupancies and energy deposits from pileup collisions at the expected large FCC-hh/SppC luminosities.

An overview of the impact of the first three years of LHC operation on two of the most important open questions in astroparticle physics is presented. Measurements in proton-proton collisions at the energy frontier that provide valuable information on the identity of the highest-energy particles in the cosmos as well as new constraints on the nature of dark matter, are summarized.

The workshop on "Hadron-Hadron and Cosmic-Ray Interactions at multi-TeV Energies" held at the ECT* centre (Trento) in Nov.-Dec. 2010 gathered together both theorists and experimentalists to discuss issues of the physics of high-energy hadronic interactions of common interest for the particle, nuclear and cosmic-ray communities. QCD results from collider experiments -- mostly from the LHC but also from the Tevatron, RHIC and HERA -- were discussed and compared to various hadronic Monte Carlo generators, aiming at an improvement of our theoretical understanding of soft, semi-hard and hard parton dynamics. The latest cosmic-ray results from various ground-based observatories were also presented with an emphasis on the phenomenological modeling of the first hadronic interactions of the extended air-showers generated in the Earth atmosphere. These mini-proceedings consist of an introduction and short summaries of the talks presented at the meeting.

First data on inclusive particle production measured in proton-proton collisions at the Large Hadron Collider (LHC) are compared to predictions of various hadron-interaction Monte Carlos (QGSJET, EPOS and SIBYLL) used commonly in high-energy cosmic-ray physics. While reasonable overall agreement is found for some of the models, none of them reproduces consistently the sqrt(s) evolution of all the measured observables. We discuss the implications of the new LHC data for the modeling of the non-perturbative and semihard parton dynamics in hadron-hadron and cosmic-rays interactions at the highest energies studied today.

The determination of the primary energy and mass of ultra-high-energy cosmic-rays (UHECR) generating extensive air-showers in the Earth's atmosphere, relies on the detailed modeling of hadronic multiparticle production at center-of-mass (c.m.) collision energies up to two orders of magnitude higher than those studied at particle colliders. The first Large Hadron Collider (LHC) data have extended by more than a factor of three the c.m. energies in which we have direct proton-proton measurements available to compare to hadronic models. In this work we compare LHC results on inclusive particle production at energies sqrt(s) = 0.9, 2.36, and 7 TeV to predictions of various hadronic Monte Carlo (MC) models used commonly in cosmic-ray (CR) physics (QGSJET, EPOS and SIBYLL). As a benchmark with a standard collider physics model we also show PYTHIA (and PHOJET) predictions with various parameter settings. While reasonable overall agreement is found for some of the MC, none of them reproduces consistently the sqrt(s) evolution of all the observables. We discuss implications of the new LHC data for the description of cosmic-ray interactions at the highest energies.

Proceedings of the 13th International Conference on Elastic and Diffractive Scattering (Blois Workshop) - Moving Forward into the LHC Era

We present and compare the predictions of various cosmic-ray Monte Carlo models for the energy (dE/deta) and particle (dN/deta) flows in p-p, p-Pb and Pb-Pb collisions at sqrt(s) = 14, 8.8, and 5.5 TeV respectively, in the range covered by forward LHC detectors like CASTOR or TOTEM (5.2<|eta|<6.6) and ZDC or LHCf (|eta|>8.1 for neutrals).