High-energy cosmic ray muons in the Earth’s atmosphere

We present the calculations of the atmospheric muon fluxes at energies 10–107 GeV based on a numerical-analytical method for solving the hadron-nucleus cascade equations. It allows the non-power-law behavior of the primary cosmic ray (PCR) spectrum, the violation of Feynman scaling, and the growth of the total inelastic cross sections for hadron-nucleus collisions with increasing energy to be taken into account. The calculations have been performed for a wide class of hadron-nucleus interaction models using directly the PCR measurements made in the ATIC-2 and GAMMA experiments and the parameterizations of the primary spectrum based on a set of experiments. We study the dependence of atmospheric muon flux characteristics on the hadronic interaction model and the influence of uncertainties in the PCR spectrum and composition on the muon flux at sea level. Comparison of the calculated muon energy spectra at sea level with the data from a large number of experiments shows that the cross sections for hadron-nucleus interactions introduce the greatest uncertainty in the energy region that does not include the knee in the primary spectrum.