Recently the atmospheric muon spectra at high energies were reconstructed for two ranges of zenith angles, basing on the events collected with the IceCube detector. These measurements reach high energies at which the contribution to atmospheric muon fluxes from decays of short-lived hadrons is expected. Latest IceCube measurements of the high-energy atmospheric muon spectrum indicate the presence of prompt muon component at energies above 500 TeV. In this work, the atmospheric conventional muon flux in the energy range 10 GeV - 10 PeV is calculated using a set of hadronic models in combination with known parameterizations of the cosmic ray spectrum by Zatsepin $\&$ Sokolskaya and by Hillas $\&$ Gaisser. The calculation of the prompt muons with use of the quark-gluon string model (QGSM) reproduces the muon data of the IceCube experiment. Nevertheless, an additional contribution to the prompt muon component is required to describe the IceCube muon spectra in case if a charm production model predicts the appreciably lower prompt lepton flux as compared with QGSM. This addition, apparently originating from rare decay modes of the short-lived unflavored mesons $\eta, \eta^\prime, \rho, \omega, \phi$, might ensure the competing contribution to the high-energy atmospheric muon flux.

We calculate the atmospheric neutrino fluxes in the energy range $100$ GeV -- $10$ PeV with the use of several known hadronic models and a few parametrizations of the cosmic ray spectra which take into account the knee. The calculations are compared with the atmospheric neutrino measurements by Frejus, AMANDA, IceCube and ANTARES. An analytic description is presented for the conventional ($\nu_\mu+\bar\nu_\mu$) and ($\nu_e+\bar\nu_e$) energy spectra, averaged over zenith angles, which can be used to obtain test data of the neutrino event reconstruction in neutrino telescopes. The sum of the calculated atmospheric $\nu_\mu$ flux and the IceCube best-fit astrophysical flux gives the evidently higher flux as compared to the IceCube59 data, giving rise the question concerning the hypothesis of the equal flavor composition of the high-energy astrophysical neutrino flux. Calculations show that the transition from the atmospheric electron neutrino flux to the predominance of the astrophysical neutrinos occurs at $30-100$ TeV if the prompt neutrino component is taken into consideration. The neutrino flavor ratio, extracted from the IceCube data, does not reveal the trend to increase with the energy as is expected for the conventional neutrino flux in the energy range $100$ GeV - $30$ TeV. A depression of the ratio $R_{\nu_\mu/\nu_e}$ possibly indicates that the atmospheric electron neutrino flux obtained in the IceCube experiment contains an admixture of the astrophysical neutrinos in the range $10-50$ TeV.

High-energy neutrinos from decays of mesons, produced in collisions of cosmic ray particles with air nuclei, form unavoidable background for detection of astrophysical neutrinos. More precise calculations of the high-energy neutrino spectrum are required since measurements in the IceCube experiment reach the intriguing energy region where a contribution of the prompt neutrinos and/or astrophysical ones should be discovered. Basing on the referent hadronic models QGSJET II-03, SIBYLL 2.1, we calculate high-energy spectra, both of the muon and electron atmospheric neutrinos, averaged over zenith-angles. The computation is made using three parameterizations of cosmic ray spectra which include the knee region. All calculations are compared with the atmospheric neutrino measurements by Frejus and IceCube. The prompt neutrino flux predictions obtained with thequark-gluon string model (QGSM) for the charm production by Kaidalov & Piskunova do not contradict to the IceCube measurements and upper limit on the astrophysical muon neutrino flux. Neutrino flavor ratio, $\phi_{\nu_ mu}/\phi_{\nu_e}$, extracted from IceCube data decreases in the energy range $0.1 - 5$ TeV energy contrary to that one might expect from the conventional neutrino flux. Presumable reasons of such behavior are: i) early arising contribution from decays of charmed particle, differing from predictions of present models, ii) revealed diffuse flux of astrophysical electron neutrinos. The likely diffuse flux of astrophysical neutrinos related to the PeV neutrino events, detected in the IceCube experiment, leads to a decrease of the flavor ratio at the energy below 10 TeV, that is in qualitative agreement with a rough approximation for theflavor ratio obtained from the IceCube data.

High-energy neutrinos, arising from decays of mesons produced through the collisions of cosmic ray particles with air nuclei, form the background in the astrophysical neutrino detection problem. An ambiguity in high-energy behavior of pion and especially kaon production cross sections for nucleon-nucleus collisions may affect essentially the calculated neutrino flux. We present results of the calculation of the energy spectrum and zenith-angle distribution of the muon and electron atmospheric neutrinos in the energy range 10 GeV to 10 PeV. The calculation was performed with usage of known hadronic models (QGSJET-II-03, SIBYLL 2.1, Kimel & Mokhov) for two of the primary spectrum parametrizations, by Gaisser & Honda and by Zatsepin & Sokolskaya. The comparison of the calculated muon neutrino spectrum with the IceCube40 experiment data make it clear that even at energies above 100 TeV the prompt neutrino contribution is not so apparent because of tangled uncertainties of the strange (kaons) and charm (D-mesons) particle production cross sections. An analytic description of calculated neutrino fluxes is presented.

High-energy neutrinos, arising from decays of mesons that were produced through the cosmic rays collisions with air nuclei, form unavoidable background noise in the astrophysical neutrino detection problem. The atmospheric neutrino flux above 1 PeV should be supposedly dominated by the contribution of charmed particle decays. These (prompt) neutrinos originated from decays of massive and shortlived particles, $D^\pm$, $D^0$, $\bar{D}{}^0$, $D_s^\pm$, $\Lambda^+_c$, form the most uncertain fraction of the high-energy atmospheric neutrino flux because of poor explored processes of the charm production. Besides, an ambiguity in high-energy behavior of pion and especially kaon production cross sections for nucleon-nucleus collisions may affect essentially the calculated neutrino flux. There is the energy region where above flux uncertainties superimpose. A new calculation presented here reveals sizable differences, up to the factor of 1.8 above 1 TeV, in muon neutrino flux predictions obtained with usage of known hadronic models, SIBYLL 2.1 and QGSJET-II. The atmospheric neutrino flux in the energy range $10-10^7$ GeV was computed within the 1D approach to solve nuclear cascade equations in the atmosphere, which takes into account non-scaling behavior of the inclusive cross-sections for the particle production, the rise of total inelastic hadron-nucleus cross-sections and nonpower-law character of the primary cosmic ray spectrum. This approach was recently tested in the atmospheric muon flux calculations [1]. The results of the neutrino flux calculations are compared with the Frejus, AMANDA-II and IceCube measurement data.

In the near future the energy region above few hundreds of TeV may really be accessible for measurements of the atmospheric muon spectrum by the IceCube array. Therefore one expects that muon flux uncertainties above 50 TeV, related to a poor knowledge of charm production cross sections and insufficiently examined primary spectra and composition, will be diminished. We give predictions for the very high-energy muon spectrum at sea level, obtained with the three hadronic interaction models, taking into account also the muon contribution due to decays of the charmed hadrons.

We study the influence of hadron interaction features on the high-energy atmospheric neutrino spectrum. The 1D calculation is performed with use of the known high-energy hadronic models, SIBYLL 2.1, QGSJET-II, Kimel and Mokhov, for the parameterizations of primary cosmic ray spectra issued from the measurement data. The results are compared with the Frejus data and AMANDA-II measurements as well as with other calculations. Sizable difference of the neutrino fluxes (up to the factor of 1.8 at 1 TeV) that are obtained with the SIBYLL and QGSJET-II appears to be rather unexpected keeping in mind the hadron and muon flux calculations in the same energy region.

A new calculation of the atmospheric fluxes of cosmic-ray hadrons and muons in the energy range 10-10^5 GeV has been performed for the set of hadron production models, EPOS 1.6, QGSJET II-03, SIBYLL 2.1, and others that are of interest to cosmic ray physicists. The fluxes of secondary cosmic rays at several levels in the atmosphere are computed using directly data of the ATIC-2, GAMMA experiments, and the model proposed recently by Zatsepin and Sokolskaya as well as the parameterization of the primary cosmic ray spectrum by Gaisser and Honda. The calculated energy spectra of the hadrons and muon flux as a function of zenith angle are compared with measurements as well as other calculations. The effect of uncertainties both in the primary cosmic ray flux and hadronic model predictions on the spectra of atmospheric hadrons and muons is considered.

We calculate the cosmic ray muon flux at ground level using directly the primary cosmic ray spectrum and composition measured in the ATIC-2 balloon experiment. In order to extend the calculations to high energies, up to 100 TeV, we use the data of the GAMMA experiment as well as Zatsepin and Sokolskaya model. This model, supported by the ATIC-2 data, comprises contributions to the cosmic ray flux of three classes of astrophysical sources - the shocks from exploding stars, nova and supernova of different types. The muon flux computation is based on the method for solution of atmospheric hadron cascade equations in which rising total inelastic cross-sections of hadron-nuclear interactions as well as non-power-law character of the primary cosmic ray spectrum are taken into account. The calculated muon spectrum agrees well with measurements of L3+Cosmic and BESS-TeV, CAPRICE, Frejus, MACRO, LVD as well as other experiments.

We discuss an effective numerical method for solving the transport equations for cosmic ray nucleons in the atmosphere. It is demonstrated that the nucleon attenuation lengths are strongly energy and depth dependent due to the non-power-law behavior of the primary spectrum, growth of the total inelastic cross sections with energy, and scaling violation in the nucleon-nucleus interactions. The numerical results are compared with the available experimental data.

The vertical sea-level muon spectrum at energies above 1 GeV and the underground/underwater muon intensities at depths up to 18 km w.e. are calculated. The results are particularly collated with a great body of the ground-level, underground, and underwater muon data. In the hadron-cascade calculations, the growth with energy of inelastic cross sections and pion, kaon, and nucleon generation in pion-nucleus collisions are taken into account. For evaluating the prompt muon contribution to the muon flux, we apply two phenomenological approaches to the charm production problem: the recombination quark-parton model and the quark-gluon string model. To solve the muon transport equation at large depths of homogeneous medium, a semi-analytical method is used. The simple fitting formulas describing our numerical results are given. Our analysis shows that, at depths up to 6-7 km w. e., essentially all underground data on the muon intensity correlate with each other and with predicted depth-intensity relation for conventional muons to within 10%. However, the high-energy sea-level data as well as the data at large depths are contradictory and cannot be quantitatively decribed by a single nuclear-cascade model.

The dependence of the semileptonic kaon decay form factors on the 4-momentum transfer to the lepton pair causes certain changes in the differential decay rates and, hence, it should affect the conventional atmospheric neutrino flux and neutrino flavour ratio. As a case in point, we included these rates into the calculation of the atmospheric neutrino flux at energies 1 to 100 TeV. The calculated neutrino spectra are between the earlier predictions while the neutrino flavor ratio, R, is somewhat affected by the kaon form factors. The R proves to be very sensitive to the contribution of neutrinos from decay of charmed particles, providing an additional method to test the charm production models in future experiments with large-volume neutrino telescopes.