We examine the uncertainty of the calculation of the atmospheric neutrino flux and present a way to reduce it using accurately measured atmospheric muon flux. Considering the difference of the hadronic interaction model and the real one as a variation of hadronic interaction, we find a quantitative estimation method for the error of the atmospheric neutrino flux calculation from the residual of the reconstruction of the atmospheric muon flux observed in a precision experiment, by the study of atmospheric neutrino and muon fluxes response to the variation of hadronic interaction. However, the efficiencty of this method is largely dependent on the observation site of the atmospheric muon flux, as the relation of the error of the atmospheric neutrino flux calculation and the residual of the reconstruction of the atmospheric muon flux is also largely dependent on the muon observation site, especially for the low energy neutrinos. We calculate several observation sites, near Kamioka at sea level, same but 2770m a.s.l.., Hanle India (4500m a.s.l.), and at Balloon altitude ($\sim$ 32km). Then we estimate how stringently can the atmospheric muon reduce the error in the calculation of the atmospheric neutrino flux. We also discuss on the source of error which is difficult to reduce by only the observation of atmospheric muon.

In this paper, we extend the calculation of the atmospheric neutrino flux~\cite{hkkm2004,hkkms2006,hkkm2011} to the sites in polar and tropical regions. In our earliest full 3D-calculation~\cite{hkkm2004}, we used DPMJET-III~\cite{dpm} for the hadronic interaction model above 5~GeV, and NUCRIN~\cite{nucrin} below 5~GeV. We modified DPMJET-III as in Ref.~\cite{hkkms2006} to reproduce the experimental muon spectra better, mainly using the data observed by BESS group~\cite{BESSTeVpHemu}. In a recent work~\cite{hkkm2011}, we introduced JAM interaction model for the low energy hadronic interactions. JAM is a nuclear interaction model developed with PHITS (Particle and Heavy-Ion Transport code System)~\cite{phits}. In Ref.~\cite{hkkm2011}, we could reproduce the observed muon flux at the low energies at balloon altitude with DPMJET-III above 32 GeV and JAM below that better than the combination of DPMJET-III above 5~GeV and NUCRIN below that. Besides the interaction model, we have also improved the calculation scheme according to the increase of available computational power, such as the "virtual detector correction" introduced in Ref.~\cite{hkkms2006} and the optimization of it in Ref.~\cite{hkkm2011}. The statistics of the Monte Carlo simulation is also improved at every step of the work.

We present the calculation of the atmospheric neutrino fluxes for the neutrino experiments proposed at INO, South Pole and Pyh\"asalmi. Neutrino fluxes have been obtained using ATMNC, a simulation code for cosmic ray in the atmosphere. Even using the same primary flux model and the interaction model, the calculated atmospheric neutrino fluxes are different for the different sites due to the geomagnetic field. The prediction of these fluxes in the present paper would be quite useful in the experimental analysis.

We present the calculation of the atmospheric neutrino fluxes with an interaction model named JAM, which is used in PHITS (Particle and Heavy-Ion Transport code System). The JAM interaction model agrees with the HARP experiment a little better than DPMJET-III. After some modifications, it reproduces the muon flux below 1~GeV/c at balloon altitudes better than the modified-DPMJET-III which we used for the calculation of atmospheric neutrino flux in previous works. Some improvements in the calculation of atmospheric neutrino flux are also reported.

Using the ``modified DPMJET-III'' model explained in the previous paper, we calculate the atmospheric neutrino flux. The calculation scheme is almost the same as HKKM04 \cite{HKKM2004}, but the usage of the ``virtual detector'' is improved to reduce the error due to it. Then we study the uncertainty of the calculated atmospheric neutrino flux summarizing the uncertainties of individual components of the simulation. The uncertainty of $K$-production in the interaction model is estimated by modifying FLUKA'97 and Fritiof 7.02 so that they also reproduce the atmospheric muon flux data correctly, and the calculation of the atmospheric neutrino flux with those modified interaction models. The uncertainties of the flux ratio and zenith angle dependence of the atmospheric neutrino flux are also studied.

We have studied the hadronic interaction for the calculation of the atmospheric neutrino flux by summarizing the accurately measured atmospheric muon flux data and comparing with simulations. We find the atmospheric muon and neutrino fluxes respond to errors in the $\pi$-production of the hadronic interaction similarly, and compare the atmospheric muon flux calculated using the HKKM04 code with experimental measurements. The $\mu^++\mu^-$ data show good agreement in the 1$\sim$30 GeV/c range, but a large disagreement above 30 GeV/c. The the $\mu^+/\mu^-$ ratio shows sizable differences at lower and higher momenta for opposite directions. As the disagreements are considered to be due to assumptions in the hadronic interaction model, we try to improve it phenomenologically based on the quark parton model. The improved interaction model reproduces the observed muon flux data well. The calculation of the atmospheric neutrino flux will be reported in the following paper

We have revised the calculation of the flux of atmospheric neutrinos based on a 3-dimensional scheme with the realistic IGRF geomagnetic model. The primary flux model has been revised, based on the AMS and BESS observations, and the interaction model updated to DPMJET-III. With a fast simulation code and computer system, the statistical errors in the Monte Carlo study are negligible. We estimate the total uncertainty of the atmospheric neutrino flux prediction is reduced to $\lesssim$ 10 % below 10 GeV. The `3-dimensional effects' are found to be almost the same as the study with the dipole magnetic field, but the muon curvature effect remains up to a few 10 GeV for horizontal directions. The uncertainty of the absolute normalization of the atmospheric neutrino is still large above 10 GeV due to the uncertainty of the primary cosmic ray flux above 100 GeV. However, the zenith angle variation is not affected by these uncertainties.

Precise knowledge of the hadronic interaction between primary cosmic rays and atmospheric nuclei is very important and fundamental to study atmospheric neutrinos and their oscillations. We studied hadronic interaction models using the data of primary and secondary cosmic rays observed by BESS experiments. By comparing the observed spectra with the ones calculated by several interaction models, we found DPMJET-III is most favored among them.