Astroparticle Physics at the Forward Physics Facility

Dennis Soldin Department of Physics and Astronomy, University of Utah,
Salt Lake City, UT 84112, USA
Abstract

High-energy collisions at the High-Luminosity Large Hadron Collider (HL-LHC) will generate a substantial flux of particles along the beam collision axis that current LHC experiments cannot access. Multi-particle production in the far-forward region is particularly relevant to astroparticle physics. High-energy cosmic rays initiate extensive air showers (EASs) in the atmosphere, driven by hadron-ion collisions at low momentum transfer in the non-perturbative regime of QCD. Therefore, understanding high-energy hadronic interactions in the forward region is essential for interpreting EAS data and estimating backgrounds for astrophysical neutrino searches, among other applications. The Forward Physics Facility (FPF) is a proposal to construct a new underground cavern at the HL-LHC dedicated to hosting a range of experiments aimed at detecting particles beyond the current capabilities of LHC experiments. This article will provide an update on the current planning status of the FPF and emphasize its synergies with astroparticle physics.

1 Introduction

Cosmic rays with energies above 1011GeVsuperscript1011GeV10^{11}\leavevmode\nobreak\ \rm{GeV}10 start_POSTSUPERSCRIPT 11 end_POSTSUPERSCRIPT roman_GeV enter Earth’s atmosphere where they interact with air molecules, generating extensive air showers (EASs) that can be detected with large ground-based detector arrays. Determining the properties of initial cosmic rays, such as their energy and mass, relies on indirect measurements of particles detected at ground level by these arrays, which requires a simulation-based interpretation of the EAS development. The main challenge in the description of EASs lies in modeling hadronic interactions in the forward region over a wide range of energies [1, 2], which current collider facilities cannot directly probe. Figure 1 depicts simulated particle densities from proton-proton collisions (solid lines), overlaid with the pseudorapidity (η𝜂\etaitalic_η) ranges relevant to current LHC experiments [3]. Dashed lines show estimated muon densities, Nμsubscript𝑁𝜇N_{\mu}italic_N start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT, produced by these particles, assuming NμElab0.93proportional-tosubscript𝑁𝜇superscriptsubscript𝐸lab0.93N_{\mu}\propto E_{\rm{lab}}^{0.93}italic_N start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∝ italic_E start_POSTSUBSCRIPT roman_lab end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 0.93 end_POSTSUPERSCRIPT, where Elabsubscript𝐸labE_{\rm{lab}}italic_E start_POSTSUBSCRIPT roman_lab end_POSTSUBSCRIPT is the energy of secondary EAS particles in the laboratory frame. While mid-rapidity ranges minimally affect EAS particle production, the forward region (η>4𝜂4\eta>4italic_η > 4) significantly influences the EAS development, which is driven by relativistic hadron-ion collisions in the atmosphere at low momentum transfer within the non-perturbative regime of quantum chromodynamics (QCD). Given the limitations of describing hadron production from first principles, and the lack of data from existing collider experiments, simulations highly rely on various phenomenological models of hadronic interactions, yielding substantial uncertainties. Hence, precise measurements of hadronic interactions under controlled experimental conditions, particularly in the far-forward region, are crucial for validating and refining existing EAS models.

Refer to caption
Figure 1: Simulated densities of particles [3] in arbritrary units (solid lines) in proton-proton collisions using EPOS-LHC. Dashed lines show the estimated number of muons produced by these particles, assuming an equivalent energy for the fixed target collisions in the laboratory frame, Elabsubscript𝐸labE_{\rm{lab}}italic_E start_POSTSUBSCRIPT roman_lab end_POSTSUBSCRIPT, and NμElab0.93proportional-tosubscript𝑁𝜇superscriptsubscript𝐸lab0.93N_{\mu}\propto E_{\rm{lab}}^{0.93}italic_N start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ∝ italic_E start_POSTSUBSCRIPT roman_lab end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 0.93 end_POSTSUPERSCRIPT.

2 The Forward Physics Facility

The Forward Physics Facility [4, 5] (FPF) is a proposal to build a new underground cavern at CERN, designed to accommodate a suite of experiments focused on the far-forward region during the high-luminosity LHC (HL-LHC) era. Current LHC detectors have uninstrumented areas along the beam line, potentially missing physics opportunities arising from the substantial flux of particles produced in the forward direction. Without the FPF, the HL-LHC would lack sensitivity to neutrinos and many particles proposed by various models beyond the Standard Model, including new force carriers, sterile neutrinos, and axion-like particles, for example. However, recent pathfinder experiments operating in the forward region at the LHC have successfully observed collider neutrinos directly for the first time and have demonstrated the potential for discovering new physics [6, 7]. A diverse suite of experiments will be available at the FPF to exploit these physics opportunities by detecting neutrino interactions at the highest accelerator energies, advancing our understanding of particle interactions in the far-forward region.

The proposed site at CERN is situated at a depth of 88m88m88\,\rm{m}88 roman_m along the line-of-sight (LOS) of the ATLAS collision axis and is located 627m627m627\,\rm{m}627 roman_m west of the interaction point, shielded by over 200m200m200\,\rm{m}200 roman_m of rock. In an updated baseline layout, the facility will measure approximately 75m75m75\,\rm{m}75 roman_m in length and 11.8m11.8m11.8\,\rm{m}11.8 roman_m in internal width, providing the infrastructure needed to host a diverse set of experiments aimed at exploring various physics phenomena at pseudorapidities above η7similar-to𝜂7\eta\sim 7italic_η ∼ 7.

2.1 Experiments

Currently, plans are underway to host four experiments at the FPF, each utilizing distinct detector technologies tailored to specific physics objectives [5]. The proposed experiments are shown in the baseline layout in Fig. 2:

  • FASER2: A magnetic tracking spectrometer, designed to search for light and weakly-interacting states, including new force carriers, sterile neutrinos, axion-like particles, among others, and to distinguish ν𝜈\nuitalic_ν and ν¯¯𝜈\bar{\nu}over¯ start_ARG italic_ν end_ARG charged current scattering in the upstream detectors.

  • FASERν𝜈\nuitalic_ν2: An on-axis emulsion detector, with pseudorapidity range η>8.4𝜂8.4\eta>8.4italic_η > 8.4, that will detect TeV neutrinos with unparalleled spatial resolution, including tau neutrinos.

  • FLArE: A noble liquid fine-grained time projection chamber to detect neutrinos and search for light dark matter with high kinematic resolution and wide dynamic range.

  • FORMOSA: A detector composed of scintillating bars, with world-leading sensitivity to millicharged particles across a large range of masses.

These experiments will explore a wide range of physics phenomena [4, 5], such as new particles, neutrinos, dark matter, dark sectors, and QCD. In this article, we will emphasize the unique interdisciplinary opportunities at the FPF for studies at the intersection of high-energy particle physics and astroparticle physics.

Refer to caption
Figure 2: The proposed baseline layout of the Forward Physics Facility, located on the collision axis line-of-sight of the ATLAS experiment. The FPF will be 75m75m75\,\rm{m}75 roman_m long and 11.8m11.8m11.8\,\rm{m}11.8 roman_m wide and will house a diverse set of experiments to explore the many physics opportunities in the far-forward region.

3 Astroparticle Physics at the FPF

Measurements conducted at the FPF will probe high-energy hadronic interactions in the far-forward region. These measurements aim to enhance the modeling of such interactions in the atmosphere, thereby reducing uncertainties in air shower observations and contributing to a better understanding of the properties of the highest-energy cosmic rays [1, 3]. Moreover, atmospheric neutrinos generated in extensive air showers in the far-forward region constitute a significant background for searches for high-energy astrophysical neutrinos using large-scale neutrino telescopes [24, 25, 26]. Therefore, measurements at the FPF will play a crucial role in refining our understanding of the atmospheric neutrino flux and minimizing associated uncertainties in searches for high-energy astrophysical neutrinos. The direct connections between measurements at the FPF and astroparticle physics will be further explored in the following.

3.1 Light Hadron Production

Muons serve as tracers of hadronic interactions, making their measurement in extensive air showers crucial for testing hadronic interaction models. Over the past two decades, various EAS experiments have reported discrepancies between model predictions and experimental data, referred to as the muon puzzle in EASs [3]. In particular, evidence for a deficit in muon numbers compared to simulations has been observed in analyses from the Pierre Auger Observatory [8, 9] and a systematic meta-analysis of data from nine air shower experiments revealed an energy-dependent trend of these discrepancies with high statistical significance [10, 11, 12]. Recent studies suggest that these observed discrepancies may have a very complex nature [13] and indicate severe deficits in our understanding of particle physics which are currently not understood [3].

Muons detected in air shower experiments range typically from a few to tens of GeV. They are produced at the end of the cascade of hadronic interactions spanning up to approximately ten generations, predominantly through soft hadron production processes. Figure 1 shows that hadron production at forward pseudorapidities has the largest influence on muon generation in EASs. The sensitivity to these hadrons is substantial, where even small deviations of around 5%percent55\%5 % in their multiplicity or type can have sizable impact the resulting muon flux [14].

Refer to caption
Figure 3: Neutrino energy spectra [5, 15] for electron neutrinos (left) and muon neutrinos (right) passing through FASERν𝜈\nuitalic_ν2 for an integrated luminosity of 3ab13superscriptab13\,\rm{ab}^{-1}3 roman_ab start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. The different production modes are shown separately, i.e., pion decays (red), kaon decays (orange), hyperon decays (magenta), and charm decays (blue). The predictions are obtained from SIBYLL-2.3d (solid), DPMJET-III.2017 (short dashed), EPOS-LHC (long dashed), QGSJET-II.04 (dotted), and Pythia 8.2 using soft-QCD processes (dot-dashed) and with hard-QCD processes for charm production (double-dot-dashed).

Dedicated measurements at the FPF will significantly enhance our understanding of light hadron production in the far-forward region. The ratio of electron neutrino (νesubscript𝜈𝑒\nu_{e}italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT) to muon neutrino (νμsubscript𝜈𝜇\nu_{\mu}italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT) fluxes measured in FPF experiments serves as an indirect measure of the ratio of charged kaons (K𝐾Kitalic_K) to pions (π𝜋\piitalic_π). Electron neutrino fluxes are predominantly produced from kaons, while muon neutrinos originate from both pion and kaon decays. However, νesubscript𝜈𝑒\nu_{e}italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and νμsubscript𝜈𝜇\nu_{\mu}italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT have distinct energy spectra, enabling their differentiation. In addition, neutrinos from pion decay are more concentrated around the LOS than those from kaons, given that mπ<mKsubscript𝑚𝜋subscript𝑚𝐾m_{\pi}<m_{K}italic_m start_POSTSUBSCRIPT italic_π end_POSTSUBSCRIPT < italic_m start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT, and thus neutrinos from pions obtain less additional transverse momentum than those from kaon decays. Thereby, the closeness of the neutrinos to the LOS, or equivalently their rapidity distribution, can be used to disentangle different neutrino origins to get an estimate of the pion-to-kaon ratio. Figure 3 shows predictions of the neutrino energy spectra for νesubscript𝜈𝑒\nu_{e}italic_ν start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and νμsubscript𝜈𝜇\nu_{\mu}italic_ν start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT interacting in FASERν𝜈\nuitalic_ν2, assuming an integrated luminosity of 3ab13superscriptab13\,\rm{ab}^{-1}3 roman_ab start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. These predictions are derived from various models, including SIBYLL-2.3d [16], DPMJET-III [17], EPOS-LHC [18], QGSJET-II.04 [19], and Pythia 8 [20]. The differences in the resulting fluxes exceed a factor of two, which is much larger than the expected statistical uncertainties at the FPF [15]. Since the muon puzzle is assumed to be of soft-QCD origin, there is also a strong connection to the QCD program of the FPF and dedicated QCD measurements [5] will further help to understand particle production in EAS.

Refer to caption
Figure 4: Neutrino energy spectra [22, 23] for electron neutrinos (left) and muon neutrinos (right) passing through the FLArE detector (10-ton target). The vertical axis shows the number of neutrinos that pass the detector’s cross-sectional area of 1m21superscriptm21\,\mathrm{m}^{2}1 roman_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for an integrated luminosity of 3ab13superscriptab13\,{\rm ab}^{-1}3 roman_ab start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT: pion decays (red), kaon decays (orange), hyperon decays (magenta), and charm decays (blue). The different line styles correspond to predictions obtained from SIBYLL-2.3d by varying fssubscript𝑓𝑠f_{s}italic_f start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT.

A potential key to understanding the muon puzzle is the universal strangeness and baryon enhancement observed by ALICE in collisions at mid-rapidity [21]. This enhancement depends solely on the particle multiplicity and not on the details of the collision system, allowing predictions of hadron composition in EAS in a phase space beyond current collider capabilities. If the enhancement increases in the forward region, it will impact muon production in EASs and could be traced by the ratio of charged kaons to pions measured at the FPF. Figure 4 shows predictions from a simple toy model [22, 23] based on SIBYLL-2.3d, where strangeness enhancement is introduced by allowing to swap πK𝜋𝐾\pi\rightarrow Kitalic_π → italic_K with a probability fssubscript𝑓𝑠f_{s}italic_f start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT at large pseudorapidities (η>4𝜂4\eta>4italic_η > 4). For fs=0.1subscript𝑓𝑠0.1f_{s}=0.1italic_f start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.1 (fs=0.2subscript𝑓𝑠0.2f_{s}=0.2italic_f start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.2), the predicted electron neutrino flux is a factor of 1.6 (2.2) higher at its maximum than the baseline prediction. These differences are significantly larger than the anticipated uncertainties at the FPF. It has been shown that for 0.4<fs<0.60.4subscript𝑓𝑠0.60.4<f_{s}<0.60.4 < italic_f start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT < 0.6 this simple model can partially accommodate data from the Pierre Auger Observatory [22], thus providing a potential explanation for the muon puzzle. This example demonstrates how the FPF will be able to uniquely test and constrain hadronic interaction models, fundamentally improving our understanding of multi-particle production in EASs.

In addition, it might also be possible to use forward-going muons to constrain the production of pions and kaons. The muon flux at FASER2, for example, is estimated to be approximately 1kHz1kHz1\,\rm{kHz}1 roman_kHz per cm2superscriptcm2\rm{cm}^{2}roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Studies of these muons could provide complementary information to determine the ratio of charged pions to kaons. However, these measurements may be challenging because the origin of the muon flux at the FPF is not well understood and further investigations based on detailed simulations are currently ongoing.

3.2 Charm Hadron Production

High-energy neutrinos of astrophysical origin are routinely observed by large-scale neutrino telescopes. The upcoming generation of telescopes, such as IceCube-Gen2 [25] and KM3NeT [26], are expected to detect significantly more astrophysical neutrinos. However, atmospheric neutrinos produced in EASs in the atmosphere present an irreducible background for these searches. Completely eliminating atmospheric background events is experimentally not feasible, instead, these backgrounds must be estimated and subtracted. The acceptance for these events is usually determined through simulations, highlighting the need for precise experimental data to test and refine these simulations, thereby reducing associated uncertainties. Accurate understanding of cosmic sources thus requires a comprehensive understanding of neutrino production in high-energy hadron interactions, particularly from heavy hadron decays.

Atmospheric neutrinos are generated by the semileptonic decays of hadrons in EAS, typically from π𝜋\piitalic_π and K𝐾Kitalic_K decays (see also Section 3.1). This is known as the conventional neutrino flux, which decreases with increasing energy. At sufficiently high energies, however, atmospheric neutrinos are also generated from the semileptonic decay of heavy flavor hadrons, such as D𝐷Ditalic_D mesons, B𝐵Bitalic_B mesons, and ΛcsubscriptΛ𝑐\Lambda_{c}roman_Λ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT baryons, referred to as the prompt neutrino flux. Due to their very short decay lengths, these hadrons decay immediately to neutrinos upon production, causing the prompt neutrino flux to decrease more slowly with energy than the conventional flux. At high energies (Eν105106GeVsimilar-tosubscript𝐸𝜈superscript105superscript106GeVE_{\nu}\sim 10^{5}-10^{6}\,\rm{GeV}italic_E start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT roman_GeV), prompt atmospheric neutrinos become the main background for astrophysical neutrino searches. For neutrino energies of 106GeVsuperscript106GeV10^{6}\,\rm{GeV}10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT roman_GeV, the corresponding center-of-mass energy is around 8TeV8TeV8\,\rm{TeV}8 roman_TeV, a range that can be explored at the HL-LHC.

While Figs. 3 and 4 provide predictions of the charm contribution to the neutrino flux at the FPF based on various hadronic interaction models, Fig. 5 shows several predictions for the prompt atmospheric muon neutrino flux based on pQCD calculations [27], assuming a simple broken power law (BPL) cosmic-ray spectrum. The left panel of Fig. 5 shows that theoretical predictions of the prompt atmospheric neutrino flux have substantial uncertainties. These uncertainties arise from various input parameters in the calculation of prompt fluxes, such as cross-sections for heavy flavor production, parton distribution functions, and fragmentation functions. The right panel shows that an additional significant uncertainty originates from the cosmic-ray flux assumption. Moreover, Fig. 5 (right) illustrates the individual contributions from charm hadrons produced at different rapidity ranges. At energies where the prompt component dominates, the prompt atmospheric neutrino flux originates from charm hadrons produced at rapidities of y>4.5𝑦4.5y>4.5italic_y > 4.5, accessible by the proposed FPF experiments at the HL-LHC [28].

Refer to captionRefer to caption

Figure 5: Comparison of the prompt atmospheric muon neutrino flux from various recent flux calculations [27] (left). The incident cosmic-ray flux is approximated with a broken power law in all predictions. If available, the accociated uncertainties are shown as error band. Also shown is the prompt atmospheric muon neutrino flux produced in different collider rapidity ranges [29] (right) and a prediction of the conventional atmospheric neutrino flux [30]. The predictions are shown assuming a simple broken power law (BPL) and the Gaisser H3a flux model [31] for the initial cosmic-ray flux.

The FPF will provide data with unprecedented statistics in the far-forward region, at energies that are relevant to astrophysical neutrino searches. This will provide crucial information to reduce uncertainties and strongly constrain models of charm hadron production. Thus, the FPF will be pivotal in improving predictions of the prompt atmospheric neutrino flux, which will significantly impact searches for astrophysical neutrinos by reducing the associated uncertainties.

4 Conclusions

The Forward Physics Facility is a proposal to construct a new underground cavern at the high-luminosity LHC, which will host a variety of experiments focused on particle production in the far-forward region. These experiments will enable unique neutrino measurements with unprecedented statistics, probing hadron production in a phase space that is inaccessible to any existing collider experiment but crucial for modeling high-energy hadronic interactions in the atmosphere. Consequently, the measurements at the FPF are closely linked to open questions in modern astroparticle physics. They will provide unique tests to constrain hadronic interaction models, fundamentally improving our understanding of multi-particle production in extensive air showers, and thereby notably contribute to solving the muon puzzle. This will significantly reduce the uncertainties in ground-based cosmic-ray observations which typically rely on interpretations based on these models. Additionally, measurements at the FPF will provide crucial information about charm hadron production in the forward region which will improve predictions of the prompt atmospheric neutrino flux and reduce associated uncertainties in searches for astrophysical neutrinos. Thus, the experiments at the FPF will contribute to enhance our understanding of the origin and nature of the highest energy cosmic rays and astrophysical neutrinos in the context of future multi-messenger observations [32, 33].

References

References

  • [1] K.-H. Kampert and M. Unger, Astropart. Phys. 35, 660-678 (2012).
  • [2] R. Engel, D. Heck, and T. Pierog, Ann. Rev. Nucl. Part. Sci. 61, 467-489 (2011).
  • [3] J. Albrecht et al., Astrophys. Space Sci. 367, 27 (2022).
  • [4] L. A. Anchordoqui et al., Phys. Rept. 968, 1-50 (2022).
  • [5] J. L. Feng et al., J. Phys. G 50, 030501 (2023).
  • [6] H. Abreu et al., Phys. Rev. Lett. 131, 031801 (2023).
  • [7] R. Albanese et al., Phys. Rev. Lett. 131, 031802 (2023).
  • [8] A. Aab et al. (Pierre Auger Collaboration), Phys. Rev. D 91, 032003 (2015).
  • [9] A. Aab et al. (Pierre Auger Collaboration), Phys. Rev. Lett. 117, 192001 (2016).
  • [10] H. P. Dembinski et al. (EAS-MSU, IceCube, KASCADE Grande, NEVOD-DECOR, Pierre Auger, SUGAR, Telescope Array, Yakutsk EAS Array Collaborations), EPJ Web Conf. 210, 02004 (2019).
  • [11] L. Cazon (EAS-MSU, IceCube, KASCADE Grande, NEVOD-DECOR, Pierre Auger, SUGAR, Telescope Array, Yakutsk EAS Array Collaborations), PoS ICRC2019, 214 (2019).
  • [12] D. Soldin (EAS-MSU, IceCube, KASCADE Grande, NEVOD-DECOR, Pierre Auger, SUGAR, Telescope Array, Yakutsk EAS Array Collaborations), PoS ICRC2021, 349 (2021).
  • [13] J. C. Arteaga Velazquez, PoS ICRC2023, 466 (2023).
  • [14] S. Baur et al., Phys. Rev. D 107, 094031 (2023).
  • [15] F. Kling and L. J. Nevay Phys. Rev. D 104, 113008 (2021).
  • [16] F. Riehn et al., Phys. Rev. D 102, 063002 (2020).
  • [17] J. Ranft, R. Engel, and S. Roesler et al., Nucl. Phys. B Proc. Suppl. 122, 392-395 (2003).
  • [18] T. Pierog et al., Phys. Rev. C 92, 034906 (2015).
  • [19] S. Ostapchenko, EPJ Web Conf. 208, 11001 (2019).
  • [20] T. Sjöstrand et al., Comput. Phys. Commun. 191, 159-177 (2015).
  • [21] J. Adam et al. (ALICE Collaboration), Nature Phys. 13, 353-539 (2017).
  • [22] L. A. Anchordoqui et al., JHEAP 34, 19-32 (2022).
  • [23] S. J. Sciutto et al., PoS ICRC2023, 388 (2023).
  • [24] M. G. Aartsen et al. (IceCube Collaboration), JINST 12, P03012 (2017).
  • [25] M. G. Aartsen et al. (IceCube-Gen2 Collaboration), J. Phys. G 48, 060501 (2021).
  • [26] R. Coniglione et al. (KM3NeT Collaboration), JPCS 632, 012002 (2015).
  • [27] O. Zenaiev et al., JHEP 04, 118 (2020).
  • [28] Y. S. Jeong et al., PoS ICRC2023, 968 (2023).
  • [29] Y. S. Jeong et al., PoS ICRC2021, 1218 (2021).
  • [30] M. Honda et al., Phys. Rev. D 75, 043006 (2007).
  • [31] T. K. Gaisser et al., Astropart. Phys. 35, 801-806 (2012).
  • [32] F. G. Schröder et al., BAAS 51, 131 (2019).
  • [33] A. Coleman et al., Astropart. Phys. 149, 102819 (2023).