Impact of the inelastic proton -- nucleus cross section on the prompt neutrino flux

The description of the inelastic proton -- nucleus cross section at very high energies is still an open question. The current theoretical uncertainty has direct impact on the predictions of the cosmic ray and neutrino physics observables. In this paper we consider different models for the treatment of $\sigma_{inel}^{pA}$, compare its predictions at ultrahigh cosmic ray energies and estimate the prompt neutrino flux at the neutrino energies that have been probed by the IceCube Observatory. We demonstrate that depending of the model used to describe $\sigma_{inel}^{pA}$, the predictions for the prompt neutrino flux can differ by a factor of order of three. Such result demonstrate the importance of a precise measurement of the inelastic proton -- nucleus cross section at high energies.

The description of the conventional and prompt atmospheric neutrino fluxes, produced by cosmic-ray interactions with nuclei in the atmosphere of the Earth, has been the theme of a series of studies in the last years [1][2][3][4]. Such analysis were strongly motivated by the detection of astrophysical neutrinos by the IceCube Observatory [6][7][8]. In order to determine the cosmic neutrino flux, it is fundamental to have a precise knowledge of the atmospheric neutrino flux (For a recent review see e.g. Ref. [5]). Currently, the description of the prompt contribution is a subject of intense debate, since it strongly depends on the modelling of the heavy quark production at large energies and very forward rapidities, beyond those probed at the LHC [4]. Although the LHC data on the prompt heavy quark cross sections (see e.g. Refs. [9,10]) helped us to improve the description of heavy meson production at forward rapidities and significantly reduced some of the theoretical uncertainties, the predictions obtained by different groups can still differ by a factor ≥ 2 depending on the treatment of heavy quark production at high energies and of the QCD dynamics at small values of the Bjorken -x variable [1][2][3][4]. These previous studies have mainly focused in the calculation of the prompt neutrino flux considering different approaches for the factorization of the heavy quark cross section, distinct parametrizations for the parton distribution functions as well different models for the primary incident nucleon flux. Another important ingredient in the calculation of the prompt neutrino flux at high energies is the inelastic proton -Air cross section (σ pAir inel ), which determines the magnitude of the Z -moments and, consequently, the evolution of the hadronic cascades in the atmosphere. The modelling of σ pA inel has been discussed by several authors during the last years (See e.g. Refs. [11][12][13][14]), with its predictions at high energies being largely distinct (See Fig. 1). In this paper we will estimate the impact of these different models for σ pAir inel on the prompt neutrino flux at high neutrino energies, as those probed by the IceCube Observatory and future neutrino telescopes. As we will demonstrate in what follows, the current uncertainty associated to the treatment of the proton -Air cross section at high energies is a factor of order of 3, independent of the model used to describe the primary nucleon flux. Such result demonstrate the importance of a precise measurement of σ pAir inel at Cosmic Ray energies. Initially let's present a brief review of the formalism used to estimate the prompt neutrino flux and refer the reader to Refs. [3,4] for more details. As in Ref. [4], we will calculate the prompt neutrino flux using the semi-analytical Zmoment approach, proposed many years ago in Ref. [15] and discussed in detail e.g. in Refs. [2,16]. In this approach, a set of coupled cascade equations for the nucleons, heavy mesons and leptons (and their antiparticles) fluxes is solved, with the equations being expressed in terms of the nucleon-to-hadron (Z N H ), nucleon-to-nucleon (Z N N ), hadron-tohadron (Z HH ) and hadron-to-neutrino (Z Hν ) Z-moments. These moments are inputs in the calculation of the prompt neutrino flux associated with the production of a heavy hadron H and its decay into a neutrino ν in the low-and high-energy regimes. As a example we present the definition of the heavy hadron Z-moment, which can be expressed as follows where E is the energy of the produced particle (heavy hadron), x F is the Feynman variable, φ N is the primary cosmic ray flux, σ pAir inel is the inelastic proton-Air cross section and dσ/dx F is the differential cross section for heavy hadron production. In what follows we will focus on vertical fluxes and will assume that the cosmic ray flux φ N can be described by a broken power-law (BPL) spectrum [17] or by the H3a spectrum proposed in Ref. [18], with the incident flux being represented by protons. As in Ref. [4], we will assume that the charmed hadron Z-moments can be expressed in terms of the charm Z-moment as follows: where f H is the fraction of charmed particle which emerges as a hadron H, which will be assumed to be: .080 and f Λc = 0.094 [16]. Moreover, we will disregard nuclear effects assuming that σ(pAir → cc) = 14.5 × σ(pp → cc). We will use the collinear factorization formalism to describe the heavy quark production and consider the CT14 parametrization [19] to describe the quark and gluon distributions in the proton. All other moments have been estimated as discussed in detail in Ref. [4].
We are interested in the calculation of the prompt neutrino flux at high energies. One have that the production of neutrinos at a given neutrino energy, E ν , is determined by collisions of cosmic rays with nuclei in the atmosphere at energies that are a factor of 100-1000 larger. Therefore, the magnitude of the prompt neutrino flux measured in the kinematical range that is probed by the IceCube Observatory and future neutrino telescopes is directly associated to the the modelling of the inelastic cross section [See Eq. (1)]. As briefly pointed out above, the treatment of this quantity at high energies is still an open question. The measurement of the inelastic proton-Air cross section performed by the Pierre Auger collaboration [20] at √ s = 57 TeV helped to improve our understanding of the hadronic interactions at high energies and constrain its description. In what follows we will consider some phenomenological models commonly used in the literature and the corresponding prompt neutrino flux will be estimated. In particular, we will consider four different hadronic interaction models that derive the inelastic proton-Air cross section using the Glauber formalism [21]. All of them are based on the eikonal representation but differ in the way the eikonal functions are constructed. In the QGSJET01 model [22] all eikonal functions are described via independent Pomeron exchanges, following the Gribov's reggeon framework [23,24]. An improved version of this model, denoted QGSJETII.03 model, was proposed in Ref. [25], which differs from its predecessor by including non-linear effects, which are described in terms of Pomeron-Pomeron interactions in the employed framework. Such effects are important in the description of the hadronic interactions at high energies and small impact parameters. On the other hand, the Sibyll 2.1 model [26] is based on a two-channel ("soft" + "hard") eikonal function: while the energy dependence of the "soft" channel is modelled as a sum of two power laws, one related to Pomeron exchange and another related to Reggeon exchange, as in Regge theory [27], the energy dependence of the "hard" component is modelled as a "minijet" model [30]. In this model, the steep increasing of the parton distributions at low values of Bjorken-x (i. e. high energies) in the colliding hadrons leads to the increase of the hadronic cross section. In addition, we also have considered the QCD-inspired model from ref. [28] (denoted EMM hereafter) which also parametrizes the "hard" channel as a minijet model (employing newer parton distribution functions) but handle the "soft" channel and the p min T different from Sibyll: while both quantities increase with the energy in the Sibyll 2.1, Ref. [28] assumes a constant p min T and a soft channel that is practically energy-independent, so that the hadronic cross section increases only due to hard (partonic) interactions. Finally, for completeness, we also consider the analytical parametrization for σ pAir inel proposed in Ref. [29], denoted by MFEK, which has been used in early estimates of prompt neutrino fluxes.
In Fig. 1 we present a comparison between the predictions of the phenomenological models with the experimental data for the inelastic proton -Air cross section. We have that its predictions are similar at low values of the proton energy in the laboratory frame (E lab ), with the EMM one overestimating the data in this kinematical regime and the QGSJETII.03 being a lower bound. On the other hand, at large energies they differ by ≈ 40%, with the MFEK (EMM) prediction being an upper (lower) bound. Although the experimental uncertainty is still large, one have that P. Auger (2012) data is not described by the MFEK and Sybill 2.1 models. It is important to emphasize that the more recent version of the Sybill Monte Carlo implies smaller values for σ pAir inel at large energies, with its predictions being similar to those of the QGSJET model. In what follows we will estimate the impact of these differences on the predictions for the prompt neutrino flux.
In Fig. 2 we present our predictions for the energy dependence of the prompt neutrino flux, normalized by a factor E 3 ν , which have been calculated assuming the BPL (left panel) and H3a (right panel) models for the primary cosmic ray flux. We have that the different models for the inelastic proton -Air cross section predict distint values for the neutrino flux in the kinematical range currently covered by the IceCube observatory (E ν 3 × 10 6 GeV), with the MFEK prediction being a lower bound and the QGSJETII.03 the upper one. Such result is expected from the analysis of the Eq. (1), which show that the Z pH moment is inversely proportional to σ pAir inel . A surprising aspect is the fact that the EMM prediction is the second smaller one. From the analysis of the Fig. 1 at high energies, we will expect that the resulting flux would be the larger one and similar to the QGSJETII.03 one. However, we should to take into account that in the cascade evolution, the contribution of small energies is also important. As the EMM prediction overestimate the data at low energies, it reduces the value of the moments and implies a smaller prediction of the neutrino flux. Such result indicates that the modelling of σ pAir inel should also be under control at low energies in order to obtain realistic predictions for the neutrino flux. At very large neutrino energies (beyond the IceCube energies), the predictions of all models, excluding the MFEK one, are similar for both primary fluxes considered. In order to estimate the magnitude of the difference between the predictions, in Fig. 3 we present our results for the ratio between the prompt neutrino flux predicted by the distinct models and that obtained considering the QGSJETII.03 one. The predictions are almost independent of the primary cosmic ray flux considered. We have that the two different versions of the QGSJET model differ by approximately 50 % at low energies. On the other hand, the EMM one predicts a flux that is smaller than the QGSJETII.03 by ≈ 40%. Finally, the QGSJET01 and MFEK predictions differ by ≈ 3 in the kinematical range covered by the IceCube Observatory. Such results indicate that the modelling of the inelastic proton -Air cross section is an important source of uncertainty in the predictions of the prompt neutrino flux. Certainly, a more precise measurement of σ pAir inel will be useful to constrain the description of the hadronic interactions in the cascade evolution equations.
Finally, let's summarize our main results and conclusions. One of the main shortcomings on the determination of the atmospheric neutrino flux, as many other observables related to cosmic ray physics, is associated to the fact the predictions for the observables are strongly dependent on the modelling of the strong interactions at high energies. In particular, the prompt neutrino flux at the IceCube Observatory and future neutrino telescopes depends on our knowledge about the hadronic interactions at the relevant energies, as well as several quantities, such as the primary cosmic ray spectrum, the longitudinal momentum distribution (x F ) of the incident particles and the inelastic proton-Air cross section. In this paper we have complemented previous studies and investigated the dependence on the modelling of σ pAir inel . We have considered some examples of phenomenological models that are largely used in the literature and demonstrated that the corresponding predictions for the prompt neutrino flux can differ by a factor of 3 in the kinematical range covered by the IceCube Observatory. Such result demonstrate that a future precise measurement of the inelastic proton -Air cross section is fundamental in order to derive realistic predictions of the prompt neutrino flux.