Ultrahigh-energy cosmic rays: Anomalies, QCD, and LHC data

Measurements of proton and nuclear collisions at the Large Hadron Collider at nucleon-nucleon c.m. energies up to 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 ≈ 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 ECR ≈ 1017-1020 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 ECR ≈ 1020 eV cutoff, and uncover any potential new physics responsible for the observed anomalies.