Development Of Air Showers Research Articles

Search for photons above 1018 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1017 eV. It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum (Xmax) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV (1 EeV=1018 eV) has been developed by combining the fluorescence detector-based measurement of Xmax with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km−2 sr−1 yr−1 above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society 2024

Lateral distribution cosmic ray muon coincidences up to 36 m

Primary cosmic ray particles comprise about 85 % protons, 12 % helium, 3 % iron, and heavier elements. These particles interact with the Earth's atmosphere, generating the Extensive Air Showers (EAS). Among the particles produced are pions and kaons, which decay into cosmic ray muons. In this research, the lateral distribution of cosmic ray muons was measured using two-fold coincidences. Four NaI (Tl) detectors and the associated electronics were used in the measurements of cosmic ray muons. The detectors were positioned from 0 to 36 m at regular intervals. The muon count rate was observed to decrease as the distance between the detectors increased. The measurements were fitted to the Nishimura–Kamata–Greisen (NKG) function to analyze the lateral distribution. Monte Carlo (MC) simulations of EAS were performed using the Cosmic Ray Simulations for the KAscade Grande (CORSIKA) program. The simulations made use of EPOS and GHEISHA models for high and lower energies respectively.•The measurements for the two-fold coincidence are consistent with the NKG function.•The simulated and measured data were found to be in agreement.•The knowledge gained from the lateral distribution of cosmic ray muons is essential for the understanding of the development of extensive air showers.

Results from recent analysis of KASCADE-Grande data

KASCADE and its extension array of KASCADE-Grande were devoted to measure individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The experiment has substantially contributed to investigate the energy spectrum and mass composition of cosmic rays in the transition region from galactic to extragalactic origin of cosmic rays as well as to quantify the characteristics of hadronic interaction models in the air shower development through validity tests using the multi-detector information from KASCADE-Grande. Although the data accumulation was completed in 2013, data analysis is still continuing. Recently, we investigated the reliability of the new hadronic interactions models of the SIBYLL version 2.3d only with the energy spectra from the KASCADE-Grande data. The evolution of the muon content of high energy air showers in the atmosphere is studied as well, using EPOS-LHC, SIBYLL 2.3, QGSJET-II-04 and SIBYLL 2.3c. In this talk, recent results from KASCADE-Grande and the update of the KASCADE Cosmic Ray Data Centre (KCDC) will be discussed.

Ultra-peripheral collisions of charged hadrons in extensive air showers

We discuss the electromagnetic collisions of high energy protons, pions and kaons with atmospheric nuclei. In particular, we use the equivalent photon approximation to estimate (i) the diffractive collisions where the projectile scatters inelastically off a nucleus, and (ii) the usual radiative processes (bremsstrahlung, pair production and photonuclear interactions) of these charged hadrons in the air. We then include the processes in the simulator AIRES and study how they affect the longitudinal development of extensive air showers. For 109-11 GeV proton primaries we find that they introduce a very small reduction (below 1%) in the average value of both X max and ΔX max. At a given shower age (relative slant depth from X max), these electromagnetic processes slightly increase the number of charged particles at the shower maximum and reduce the number of muons when it is old, decreasing by 1% the muon-to-(γ + e) near the ground level.

The Influence of the Strong Atmospheric Electric Fields on the Extensive Air Shower Development

The Influence of the Strong Atmospheric Electric Fields on the Extensive Air Shower Development

Study of longitudinal development of air showers in the knee energy range

Ground-based cosmic ray experiments detect cosmic ray mainly by measuring the longitudinal and lateral distribution of secondary particles produced in the extensive air shower (EAS). The EAS of cosmic ray in the knee energy region is simulated via CORSIKA software. Several simulation samples with different energy, composition and zenith angles were carried out to understand the longitudinal development of electron, muon and Cherenkov light in EAS. All the results presented were obtained assuming an observation plane at an altitude of 4400 m a.s.l. The differences of longitudinal development between electron and Cherenkov light were studied, and the reconstruction uncertainty of shower maximum for electron from Cherenkov light was estimated to be 10–15 g/cm2 for nuclei above 1 PeV. The performances of energy measurement and the composition discrimination ability based on longitudinal development were studied and compared with that from lateral distribution. It was found that number of electron per depth at its shower maximum has the smallest shower-to-shower fluctuations, but the shower-to-shower fluctuations of electron density measured at observation level was very close to it when the appropriate zenith angle was employed. The shower-to-shower fluctuations of shower maximum for electron is 50–55 g/cm2 for proton, and 20–25 g/cm2 for iron, but the composition discrimination ability between nuclei from muon density measured at observation level is much better than the shower maximum variable from longitudinal development. The hadronic model dependencies of the longitudinal development and lateral distribution were also discussed.

Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 1017 eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 1018 eV (the so-called “ankle” feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources.To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).

Testing the Compatibility of the Depth of the Shower Maximum Measurements performed at Telescope Array and the Pierre Auger Observatory

The Telescope Array and the Pierre Auger Observatory estimate the composition of ultra-high-energy cosmic rays by observing the distribution of depths of air-shower maxima, X max. Both experiments directly observe the longitudinal development of air showers using fluorescence telescopes with surface particle detectors used in conjunction to provide precision in determining air-shower geometry. The two experiments differ in the details of the analysis of events, so a direct comparison of X max distributions is not possible. The Auger – Telescope Array Composition Working Group presents their results from a technique to compare X max measurements from Auger with those of Telescope Array. In particular, the compatibility of the first two moments of the X max distributions of Auger with the data from the Black Rock Mesa and Long Ridge detectors of the Telescope Array is tested for energies above 1018.2 eV. Quantitative comparisons are obtained using air-shower simulations of four representative species made using the Sibyll 2.3d high-energy interaction model. These are weighted to fit the fractional composition seen in Auger data and reconstructed using the Telescope Array detector response and analysis methods.

Performance evaluation of LHCf-ATLAS ZDC joint measurement using proton beam

Measurements of forward neutrons in pp collisions will allow us to investigate π-p cross-section via one-pion exchange process, which are important for air shower development. However, the precision of these measurements is limited by the energy resolution of the LHCf detectors. To improve it, a joint measurement with the ATLAS ZDC was planned. In 2021, a beam test was conducted to evaluate the performance of the joint measurement of the LHCf-Arm1 and ZDC detectors using proton beams of 350 GeV at SPS. Combining the LHCf data with the ZDC data, we confirmed that the energy resolution improved from about 40% to 21.4%.

Air shower radio signal electric field orientation as measured with the Astroneu Cosmic Ray telescope.

The Astroneu cosmic ray telescope is a small scale hybrid array consisting of both scintillator counters and electromagnetic radiation detectors in the radio-wave frequency (RF) band (antennas). Astroneu was installed and operated in the area of the Hellenic Open University (HOU) campus near the city of Patras in Greece. In the present development phase, the Astroneu telescope includes two stations consisting of 3 scintillation detectors modules (SDM) and one RF antenna while a third station includes 3 particle detectors and 4 RF antennas (3SDM-4RF). In this context we present the resent results from both the 3SDM-4RF autonomous station and simulations related to the estimation of the direction of the electric field emitted during the air shower development. The electric field measured in the ground level is the superposition of the two dominant emission mechanism, the time depended transverse current induced by the geomagnetic field and the net negative charge variation at the shower front. Since the electric field emitted by the two contributions is polarized in different directions the measured electric field in the ground encloses information about the charge-excess to geomagnetic ratio (CGR). Furthermore the orientation of the electric field emitted by charge-excess mechanism is strongly depends from the distance to the shower core (the intersection of the shower axis with the ground level). In this study we use the core information as reconstructed using the radio data and simulations. The estimated charge-excess to geomagnetic ratio is in agreement with previous studies which reveals that the shower core reconstruction method is efficient. Finally we report on that CGR measurements can be used for an efficient noise rejection in a future self trigger mode

Simulating and Modeling the Extensive Air Showers Development through the Estimating the Energy of some Created Particles

Using the AIRES code and three hadronic interaction models, a simulation study of the energy of created particles due to the development of an extensive air shower was carried out (EPOS-LHC, QGSJET-II-04, and Sibyll2.3c). The main focus is put on the energy of the create particles as a function of depth for (12C,56Fe, p, and 28Si) primary particles with high primary energies (1017,1018,1019, and 1020) eV with two zenith angles 0° and 30°. The New parameters were obtained by fitting the energy of created particles curves using Gaussian and linear functions for created particles, initiated by primary particles at the energy 1020 eV. Good agreement was obtained by comparing the present results with results simulated by CORSIKA simulation for primary proton at the energy 1019 eV

Sensitivity of the large muon detector with the Tibet air-shower array to measure the primary proton spectrum between 40 and 630 TeV

Abstract The Tibet ASγ group has been continuously observing cosmic rays and cosmic gamma rays above several TeV using the muon detector array (MD) and high-density Tibet air-shower array (Tibet-III) installed on the Tibet plateau at an altitude of 4300 m. The MD is a water Cherenkov pool array with a large effective area of 3400 m2 and has an excellent capability of primary selection using the number of muons in the shower. We report the sensitivity of the proton spectrum measurements for energies 40–630 TeV obtained via Monte Carlo simulations for an air-shower experiment. It was found that protons could be separated with a purity of 90%, and the survival ratio of protons including model dependence was 14.2%–19.1% and 3.7%–7.4% at about 35 TeV and about 450 TeV, respectively. The maximum total systematic error of the proton flux depending on interaction models in air-shower development and composition models was ±37%. With a large effective area and high proton separation capability, the Tibet ASγ experiment can measure the proton spectrum in the energy range from tens to hundreds of TeV with high statistical accuracy.

Feasibility Studies of Charge Exchange Measurements in pp Collisions at the LHC

(1) Pions produced in the development of extended atmospheric cosmic ray air showers subsequently decay to muons. The measured yield of those muons is generally underestimated by current phenomenological models and event generators optimized for cosmic ray physics. The importance of those disagreements motivates the feasibility studies for testing these models at the Large Hadron Collider (LHC) energies, at the highest center-of-mass energies achievable in a laboratory. The interaction of a nucleus and a virtual pion created in a charge exchange reaction at the LHC is a similar process to those contributing to the development of air showers in case of cosmic rays. The crucial problem of such an analysis is the selection of charge exchange events with the highest possible efficiency and high purity from proton–proton collisions at the LHC. (2) For this we consider distributions of various measurable quantities given by event generators commonly used in cosmic ray physics. (3) We examine the expected distributions of energy deposited in different calorimeters of an LHC experiment. We consider the geometrical acceptance and energy resolution of the detectors at the Compact Muon Solenoid (CMS) experiment, as an example. We determine a working point cut from the various options for event selection, and compare signal and background predictions using different models for a representative simple observable, such as average transverse momentum or charge particle yield. (4) A set of event selection cuts along these considerations is proposed, with the aim of achieving optimal efficiency and purity.

Simulation of in-ice cosmic ray air shower induced particle cascades

We present simulations of the development of high-energy cosmic-ray air showers penetrating high-altitude ice layers that can be found at the polar regions. We use a combination of the corsika Monte Carlo code and the geant4 simulation toolkit, and focus on the particle cascade that develops in the ice to describe its most prominent features. We discuss the impact of the ice layer on the total number of particles in function of depth of the air shower, and we give a general parametrization of the charge distribution in the cascade front in function of ${X}_{\mathrm{max}}$ of the cosmic ray air shower, which can be used for analytical and semianalytical calculations of the expected Askaryan radio emission of the in-ice particle cascade. We show that the core of the cosmic ray air shower dominates during the propagation in ice, therefore creating an in-ice particle cascade strongly resembling a neutrino-induced particle cascade. Finally, we present the results of simulations of the Askaryan radio emission of the in-ice particle cascade, showing that the emission is dominated by the shower core, and discuss the feasibility of detecting the plasma created by the particle cascade in the ice using RADAR echo techniques.

Probing hadronic interactions with measurements from the Pierre Auger Observatory

The Pierre Auger Observatory is the largest facility in the world to study ultra-high-energy cosmic rays. It has a hybrid detection technique that combines the observation of the longitudinal development of extensive air showers and the measurement of their particles at the ground. This capability has opened the possibility to probe hadronic interactions taking place at energies well beyond those accessible by human-made accelerators. In this report, we present a selection of the latest results on hadronic interactions with measurements from the Pierre Auger Observatory. These data span over three decades in energy, showing the tension between data from the muon component of air showers and predictions based on the most updated hadronic interaction models.

Symmetries in the Superposition Model of Extensive Air Shower Development

According to the superposition principle, an extensive air shower initiated by a nucleus with energy E and mass number A can be approximated as the superposition of A proton-initiated showers each with energy E/A. The superposition principle for interactions of atomic nuclei proposes to describe nucleus-initiated extensive air showers using simulations performed for proton showers. Single detectors and systems working in tight coincidence mainly register events initiated by particles with very low energies, which are affected by major statistical fluctuations, such as those used in high schools for education and outreach purposes. Verifying whether the superposition principle is still a good approximation in the low-energy region is important for the validity of the interpretation of such measurements. We present results of the comparison of results of the superposition model with detailed simulations of showers with the CORSIKA program from the energy of 10 GeV. While the energy dependence of the mean shower parameters satisfies the superposition principle, the higher moments do not. A modification of the superposition model based on the wounded nucleon model, reducing these discrepancies, is proposed. The semi-analytical description of showers in the modified superposition model can give the density spectrum of cosmic ray particles, which is consistent with the measurements. In this paper, we present results both consistent with the superposition model and indicating the need for its modification. This modification is proposed and tested.

Testing the superposition model in small CORSIKA shower simulations

The idea of superposition in the high-energy interactions of cosmic ray nuclei and the development of extensive air showers initiated by them has been known for more than half a century. It has been thoroughly and successfully tested in a number of simulations for primary energies around 1015 and above. In this work, we will investigate its applicability to lower energies. At the lowest energies, when the shower contains on average about one charged particle (or even less), deviations from the superposition model can be seen in the simulation results. Fluctuations of higher moments of the main shower parameters are systematically broader than expected. Further studies to confirm superposition in particular in the shower longitudinal profile are in progress. A correct description of the longitudinal development of the small shower, a precise description of its fluctuations on the observational level with a correct implementation of the superposition principle will enable to construct a simple and fast phenomenological algorithm generating small showers indispensable for the interpretation of measurements made using a small local shower array and determination of the flux of single, incoherent, secondary cosmic ray particles.

Direct measurement of upward-going ultrahigh energy dark matter at the Pierre Auger Observatory

Abstract It is assumed that two types of dark matter particles exist: superheavy dark matter particles (SHDM), the mass of which ∼ inflaton mass, and light fermion dark matter (DM) particles, which are the ultrahigh energy (UHE) products of the decay of SHDM. The Earth will be taken as a detector to search for the UHE DM particles directly. These upward-going particles, which pass through the Earth and air and interact with nuclei, can be detected by the fluorescence detectors (FD) of the Pierre Auger Observatory (Auger), via fluorescent photons due to the development of an extensive air shower. The numbers and fluxes of expected UHE DM particles are evaluated in the incoming energy range between 1 EeV and 1 ZeV with the different lifetimes of decay of SHDM and mass of Z′. According to the Auger data from 2008 to 2019, the upper limit for UHE DM fluxes is also estimated at 90% confidence limit with the FD of Auger. Finally, it is reasonable to make a conclusion that UHE DM particles could be directly detected in the energy range between O(1 EeV) and O(10 EeV) with the FD of Auger. This might prove whether SHDM particles exist in the Universe.

Muon Tomography sites for Colombian volcanoes

By using a very detailed simulation scheme, we have calculated the cosmic ray background flux at 13 active Colombian volcanoes and developed a methodology to identify the most convenient places for a muon telescope to study their inner structure. Our simulation scheme considers three critical factors with different spatial and time scales: the geo-magnetic effects, the development of extensive air showers in the atmosphere, and the detector response at ground level. The muon energy dissipation along the path crossing the geological structure is modeled considering the losses due to ionization, and also contributions from radiative Bremßtrahlung, nuclear interactions, and pair production. By examining each particular volcano topography and assuming reasonable statistics for different instrument acceptances, we obtained the muon flux crossing each structure and estimated the exposure time for our hybrid muon telescope at several points around each geological edifice. After a detailed study from the topography, we have identified the best volcano to be studied, spotted the best points to place a muon telescope and estimated its time exposures for a significant statistics of muon flux. We have devised a mix of technical and logistic rules –the “rule of thumb” criteria– and found that only Cerro Machín, located at the Cordillera Central (4°29'N 75°22'W), can be feasibly studied today through muography. Cerro Negro and Chiles could be good candidates shortly.

Performance and science reach of the Probe of Extreme Multimessenger Astrophysics for ultrahigh-energy particles

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a potential NASA Astrophysics Probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. POEMMA will monitor colossal volumes of the Earth's atmosphere to detect extensive air showers (EASs) produced by extremely energetic cosmic messengers: UHECRs above 20 EeV over the full sky and cosmic neutrinos above 20 PeV. We focus most of this study on the impact of POEMMA for UHECR science by simulating the detector response and mission performance for EAS from UHECRs. We show that POEMMA will provide a significant increase in the statistics of observed UHECRs at the highest energies over the entire sky. POEMMA will be the first UHECR fluorescence detector deployed in space that will provide high-quality stereoscopic observations of the longitudinal development of air showers. Therefore, it will be able to provide event-by-event estimates of the calorimetric energy and nuclear mass of UHECRs. The particle physics in the interactions limits the interpretation of the shower maximum on an event by event basis. In contrast, the calorimetric energy measurement is significantly less sensitive to the different possible final states in the early interactions. We study the prospects to discover the origin and nature of UHECRs using expectations for measurements of the energy spectrum, the distribution of arrival direction, and the atmospheric column depth at which the EAS longitudinal development reaches maximum. We also explore supplementary science capabilities of POEMMA through its sensitivity to particle interactions at extreme energies and its ability to detect ultra-high energy neutrinos and photons produced by top-down models including cosmic strings and super-heavy dark matter particle decay in the halo of the Milky Way.