Extensive Air Showers Research Articles

Energy spectra of elemental groups of cosmic rays with the KASCADE experiment data and machine learning

We report the reconstruction of the mass component spectra of cosmic rays (protons, helium, carbon, silicon and iron) and their mean mass composition, at energies from 1.4 to 100 PeV. The results are derived from the archival data of the extensive air shower experiment KASCADE. We use a novel machine learning technique developed specifically for this reconstruction, and post-LHC hadronic interaction models: QGSJet-II.04, EPOS-LHC and Sibyll 2.3c. We have found an excess of the proton component and a deficit of intermediate and heavy nuclei components compared to the original KASCADE results. The spectra of protons and helium show a knee-like behavior at ∼ 4.4 PeV and ∼ 11 PeV, with significances 5.2σ and 3.9σ, respectively. The spectrum of the iron component has a hint (2.4σ) of a hardening at ∼ 4.5 PeV, which can be interpreted as a counterpart of a hardening in the proton spectrum at 166 TeV, recently reported by the GRAPES-3 experiment. The systematic uncertainties of our analysis were found to be smaller than those of the original KASCADE, as well as those of IceTop and TALE experiments, over the most part of the energy range studied. We also estimated separately the uncertainty related to the difference between the three mentioned hadronic interaction models. We also compute a mean logarithm mass of CR flux as a function of energy. It is in agreement with the results of IceTop, TALE and LHAASO within the uncertainties.

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.

Muon measurements at the Pierre Auger Observatory

The Pierre Auger Observatory is the world’s largest detector for observation of ultra-high-energy cosmic rays (UHECRs) (above the energy of 10^{17}1017 eV). It consists of a Fluorescence Detector (FD) and an array of particle detectors known as the Surface Detector (SD). Observations of extensive air showers by the Observatory can be used to probe hadronic interactions at high energy, in a kinematic and energy region inaccessible to experiments at man-made accelerators and to measure the muon component of the shower. Air showers induced by different primaries have different muon contents. With increasing mass of the primary cosmic ray particle, it is expected that the muon content in the corresponding air showers should also increase. Therefore, the determination of the muon component in the air shower is crucial to infer the mass of the primary particle. This is a key ingredient in the searches conducted to pinpoint the sources of UHECRs. Recent results obtained from the Pierre Auger Observatory and other experiments indicate that all the shower simulations underestimate the number of muons in the showers compared to the data. This is the so-called muon deficit. In this paper we briefly review the muon measurements, and present in more detail recent results on fluctuations in the muon number. These results provide new insights into the origin of the muon deficit in air shower simulations and constrain the models of hadronic interactions at ultrahigh energies. With the current design of the surface detectors it is also difficult to reliably separate the contributions of muons to the SD signal from the contributions of photons, electrons, and positrons. Therefore, we also present a new method to extract the muon component of the signal time traces recorded by each SD station using recurrent neural networks. The combination of such algorithms, with the future data collected by the upgraded Pierre Auger Observatory, will be a major step forward, as we are likely to achieve an unprecedented resolution in mass estimation on an event-by-event basis.

Air showers and hadronic interactions with CORSIKA 8

The CORSIKA 8 project is a collaborative effort aiming to develop a versatile C++ framework for the simulation of extensive air showers, intended to eventually succeed the long-standing FORTRAN version. I present an overview of its current capabilities, focusing on aspects concerning the hadronic and muonic shower components. In particular, I demonstrate the “cascade lineage” feature and its application to quantify the importance of certain phase-space regions in hadronic interactions for muon production. Additionally, I show first results using Pythia 8.3, which as of late is usable as interaction model in cosmic-ray applications and is currently being integrated into CORSIKA 8.

Sibyll[formula omitted]

In the last decade, an increasing number of datasets have revealed a consistent discrepancy between the number of muons measured in ultra-high-energy extensive air showers (EAS) and the numbers predicted by simulations. This gap persists despite incorporating Large Hadron Collider (LHC) data into the tuning of current hadronic interaction models, leading to the phenomenon often termed the “muon puzzle”. To gain a deeper understanding of the potential origins of this muon puzzle, we have developed Sibyll★, a series of phenomenologically modified versions of Sibyll 2.3d. In these models, we have increased muon production by altering ρ0, baryon–antibaryon pair, or kaon production in hadronic multiparticle production processes. These variants remain within bounds from provided by accelerator measurements, including those from the LHC and fixed-target experiments, notably NA49 and NA61, showing a level of consistency comparable to Sibyll 2.3d. Our findings show that these modifications can increase the muon count in EAS by up to 35%, while minimally affecting the depth of shower maximum (Xmax) and other shower variables. Additionally, we assess the impact of these modifications on various observables, including inclusive muon and neutrino fluxes and the multiplicities of muon bundles in deep underground and water/ice Cherenkov detectors. We aim for at least one of these model variants to offer a more accurate representation of EAS data at the highest energies, thereby enhancing the quality of Monte Carlo predictions used in training neural networks. This improvement is crucial for achieving more reliable data analyses and interpretations.

Probing atmospheric effects using GRAPES-3 plastic scintillator detectors

The GRAPES-3 extensive air shower (EAS) array has been designed to study cosmic rays from 1013\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{13}$$\\end{document}–1016\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{16}$$\\end{document} eV. It employs 400 scintillator detectors spread across 25,000 m2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{2}$$\\end{document}, mainly of cone-type and fiber-type, each covering a 1 m2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{2}$$\\end{document} area. These detectors record EAS particle densities and arrival times, which are crucial for determining primary particle energy and direction. A decade (2013–2022) of EAS data is analyzed to investigate the dependence of particle densities on ambient temperature and atmospheric pressure. Notably, ambient temperature exhibits a delayed response, with a more pronounced delay in fiber-type detectors, while cone-type detectors exhibit a higher observed temperature coefficient. In contrast, atmospheric pressure instantly and uniformly affects both detector types, with Monte Carlo simulations backing the observed pressure coefficient. These findings established a reliable pressure coefficient for EAS within this distinctive energy range and contributed to the refinement of correction algorithms, ultimately improving particle density precision for more accurate shower parameter estimates.

Extensive air showers and atmospheric electric fields. Synergy of space and atmospheric particle accelerators

Various particle accelerators operate in the space plasmas, filling the Galaxy with high-energy particles (primary cosmic rays). Reaching the Earth's atmosphere, these particles originate extensive air showers (EASs) consisting of millions of elementary particles (secondary cosmic rays), covering several km2 on the ground. During thunderstorms, strong electric fields modulate the energy spectra of EAS secondary particles, changing the shower size (number of EAS electrons) and altering the primary particle's estimated energy and frequency of the surface array triggers. Impulse amplifications of particle fluxes (the so-called thunderstorm ground enhancements, TGEs) manifest themselves as large peaks in the time series of count rates of particle detectors located on the Earth’s surface. Free electrons are abundant at any altitude in the atmosphere, from small to large EASs. These electrons serve as seeds for electron accelerators, which operate in the thunderous atmosphere and send particle avalanches in the direction of Earth's surface and into space (terrestrial gamma flashes, TGFs). EAS cores randomly hitting arrays of particle detectors also generate short bursts of relativistic particles. For years, particle detectors, electric field sensors, and lightning locators have gathered information about the complex interactions of secondary particle fluxes, electric fields, and lightning flashes. This information is crucial for establishing a field of high-energy physics in the atmosphere. Plain language summaryCorrelated measurements of particle fluxes modulated by strong atmospheric electric fields, registration of broadband radio and optical emission from atmospheric discharges, and registration of electric fields and various meteorological parameters lead to a better understanding of the complex processes of particle-field interactions in the terrestrial atmosphere. The cooperation of cosmic rays and atmospheric physics has led to the development of models of the origin of particle bursts recorded on the Earth's surface, vertical and horizontal profiles of electric fields, initiation of lightning flashes, etc. Interdisciplinary atmospheric science primarily requires monitoring particle fluxes around the clock by synchronized networks of identical sensors that record and store multidimensional data in databases with open, fast, and reliable access. The advances in multidimensional measurements over the past decade significantly intensified the development of new integrated models of atmospheric electricity and electron acceleration, giving more insight into understanding the modulation effects posed on EAS particles in the strong atmospheric electric fields.

Science with the ASTRI Mini-Array: From Experiment to Open Observatory

Although celestial sources emitting in the few tens of GeV up to a few TeV are being investigated by imaging atmospheric Čerenkov telescope arrays such as H.E.S.S., MAGIC, and VERITAS, at higher energies, up to PeV, more suitable instrumentation is required to detect ultra-high-energy photons, such as extensive air shower arrays, as HAWC, LHAASO, Tibet AS-γ. The Italian National Institute for Astrophysics has recently become the leader of an international project, the ASTRI Mini-Array, with the aim of installing and operating an array of nine dual-mirror Čerenkov telescopes at the Observatorio del Teide in Spain starting in 2025. The ASTRI Mini-Array is expected to span a wide range of energies (1–200 TeV), with a large field of view (about 10 degrees) and an angular and energy resolution of ∼3 arcmin and ∼10 %, respectively. The first four years of operations will be dedicated to the exploitation of Core Science, with a small and selected number of pointings with the goal of addressing some of the fundamental questions on the origin of cosmic rays, cosmology, and fundamental physics, the time-domain astrophysics and non γ-ray studies (e.g., stellar intensity interferometry and direct measurements of cosmic rays). Subsequently, four more years will be dedicated to Observatory Science, open to the scientific community through the submission of observational proposals selected on a competitive basis. In this paper, I will review the Core Science topics and provide examples of possible Observatory Science cases, taking into account the synergies with current and upcoming observational facilities.

Astrophysical approach to search for heavy neutrino decay

AbstractCosmic rays are a very valuable tool of multi‐messenger astrophysics, as they provide a very different picture of the sky. During the past decades, a large number of active astrophysical objects in our Galaxy and beyond have been discovered through the detection of gamma‐rays with Cherenkov telescopes. Cosmic rays, neutrinos have been successfully supplementing the astronomical view. Also, cosmic rays may offer to investigation of the elementary particle properties. Neutrino telescope detects the Cherenkov radiation generated in water or ice by the passage of relativistic charged particles produced by neutrino collisions with nucleons in the detector volume. Some alternative approaches have been proposed. One of them is using earth matter or mountains as a target volume for the conversion of neutrinos to leptons which then initiate extensive air showers (EAS) in the atmosphere, then showers can be detected by the Cherenkov telescope. Investigations with SHALON Cherenkov telescope have included observations of EAS from the sub‐horizontal direction . Five EAS of ∼10 TeV energies were detected with SHALON from the sub‐horizontal directions in the conditions with the zero expected number of showers. These events may be caused by the decay of a long‐lived penetrating particle entering the atmosphere from the ground and decaying in front of the telescope. As a possible explanation, two scenarios with an unstable neutrino of mass GeV and m is discussed. Remarkably, one of these models has been proposed to explain an excess of electron‐like neutrino events at MiniBooNE.

SWGO: a wide-field of view gamma-ray observatory in the southern hemisphere

Abstract The recent LHAASO and HAWC results opened the way to the search of gamma ray sources emitting at energies above 100 TeV. Both detectors are in the northern hemisphere; the need for such an observatory in the southern hemisphere is therefore clear. The goal of the SWGO collaboration is the construction of a wide field of view, high duty cycle observatory to explore the Southern hemisphere sky searching for gamma ray sources at energies above 100 GeV. Such an array will detect extensive air showers particles and must be able to select the photon originated showers from the background of the hadronic ones. The experiment must be located in a site at latitude between 10° and 30° degrees south and at an altitude above 4400 m a.s.l. The baseline detection technique chosen by the collaboration is Water Cherenkov Detectors. The array will have a central region with high fill factor (>60%) and a large (about 1 km2) outer region with a much lower fill factor (around 4–5%).

Mini-EUSO on Board the International Space Station: Mission Status and Results

The telescope Mini-EUSO has been observing, since 2019, the Earth in the ultraviolet band (290–430 nm) through a nadir-facing UV-transparent window in the Russian Zvezda module of the International Space Station. The instrument has a square field of view of 44°, a spatial resolution on the Earth surface of 6.3 km and a temporal sampling rate of 2.5 microseconds. The optics is composed of two 25 cm diameter Fresnel lenses and a focal surface consisting of 36 multi-anode photomultiplier tubes, 64 pixels each, for a total of 2304 channels. In addition to the main camera, Mini-EUSO also contains two cameras in the near infrared and visible ranges, a series of silicon photomultiplier sensors and UV sensors to manage night-day transitions. Its triggering and on-board processing allow the telescope to detect UV emissions of cosmic, atmospheric and terrestrial origin on different time scales, from a few microseconds up to tens of milliseconds. This makes it possible to investigate a wide variety of events: the study of atmospheric phenomena (lightning, transient luminous events (TLEs) such as ELVES and sprites), meteors and meteoroids; the search for nuclearites and strange quark matter; and the observation of artificial satellites and space debris. Mini-EUSO is also potentially capable of observing extensive air showers generated by ultra-high-energy cosmic rays with an energy above 1021 eV and can detect artificial flashing events and showers generated with lasers from the ground. The instrument was integrated and qualified in 2019 in Rome, with additional tests in Moscow and final, pre-launch tests in Baikonur. Operations involve periodic installation in the Zvezda module of the station with observations during the crew night time, with periodic downlink of data samples, and the full dataset being sent to the ground via pouches containing the data disks. In this work, the mission status and the main scientific results obtained so far are presented, in light of future observations with similar instruments.

Particle correlation in the cosmic ray showers around the spectral knee

In this paper, we answer the question to what extent the particles in showers initiated by particles with energies of 1014 – 1017 eV observed at the edges of extensive air showers, where average particle densities are low (< a few per m2), are correlated. The particle number multiplicity distributions in the surface detectors of the shower arrays are affected by the existence of such a correlation. They are well described by a negative binomial distribution rather than Poisson. Positive correlations cause the distributions to become broader and, in general, the variance of the observed particle number is larger than is usually assumed in the analysis of the shower registration by the Extensive Air Shower experiments. We show how this increase can alter the results of shower axis localisation algorithms.

Demonstrating Agreement between Radio and Fluorescence Measurements of the Depth of Maximum of Extensive Air Showers at the Pierre Auger Observatory.

We show, for the first time, radio measurements of the depth of shower maximum (X_{max}) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence dataset, and between a subset of air showers observed simultaneously with both radio and fluorescence techniques, a measurement setup unique to the Pierre Auger Observatory. Furthermore, we show radio X_{max} resolution as a function of energy and demonstrate the ability to make competitive high-resolution X_{max} measurements with even a sparse radio array. With this, we show that the radio technique is capable of cosmic-ray mass composition studies, both at Auger and at other experiments.

Radio measurements of the depth of air-shower maximum at the Pierre Auger Observatory

The (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of 17 km2 with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the 30–80 MHz band. Here, we report the AERA measurements of the (Xmax), a probe for mass composition, at cosmic-ray energies between 1017.5 and 1018.8 eV, which show agreement with earlier measurements with the fluorescence technique at the Pierre Auger Observatory. We show advancements in the method for radio Xmax reconstruction by comparison to dedicated sets of / air-shower simulations, including steps of reconstruction-bias identification and correction, which is of particular importance for irregular or sparse radio arrays. Using the largest set of radio air-shower measurements to date, we show the radio Xmax resolution as a function of energy, reaching a resolution better than 15 g cm−2 at the highest energies, demonstrating that radio Xmax measurements are competitive with the established high-precision fluorescence technique. In addition, we developed a procedure for performing an extensive data-driven study of systematic uncertainties, including the effects of acceptance bias, reconstruction bias, and the investigation of possible residual biases. These results have been cross-checked with air showers measured independently with both the radio and fluorescence techniques, a setup unique to the Pierre Auger Observatory. Published by the American Physical Society 2024

Study of Atmospheric Depth Profiles at LHAASO Using the MSISE-90 Model

High Altitude Cosmic Ray Observatory (LHAASO) is located at Haizi Mountain in Daocheng County, Sichuan Province. Its Wide Field of view Cherenkov Telescope Array (WFCTA) primarily studies cosmic rays by observing the Cherenkov light signals produced during extensive air showers. Calibration, simulation, and reconstruction of WFCTA are all related to atmospheric depth. The currently used atmospheric depth model is the US Standard Atmosphere Depth Profile Model. In this study, the US Standard Atmosphere Depth Profile Model is compared with the atmospheric depth profile recorded by the infrared radiometer SABER carried by the satellite TIMED at LHAASO from 14 km to 50 km, as well as with the atmospheric depth recorded by the ground meteorological station at LHAASO. The atmospheric depth obtained from the US Standard Atmosphere Model is consistently smaller. The MSISE-90 atmospheric model describes the neutral temperature and density from the Earth's surface to the thermosphere. Further research shows good consistency between the MSISE-90 atmospheric model and the atmospheric depth recorded by TIMED/SABER and the ground standard meteorological station at LHAASO. According to the MSISE-90 atmospheric model, the mean atmospheric depth profile at LHAASO is lowest in January, followed by February, March, April, November, and December, which are also the optimal observation months for WFCTA operation. The atmospheric boundary layer is highest in April, and there is about 2% diurnal variation in atmospheric depth. Using the functional form of the US Standard Atmosphere Model, the atmospheric depth profiles from 4.4 to 100 kilometers were fitted for each month, obtaining monthly atmospheric depth profile models at the LHAASO site. And compared the lateral distribution of the Cherenkov photons produced by 100 TeV cosmic-ray protons incident at a zenith angle of 30° in the MSISE-90 atmospheric model and the US Standard Atmosphere Model,with the maximum difference reaching approximately 20%.

Main parameters of extensive air showers at mountain altitudes

In high-energy physics one of the main directions is extensive air showers (EAS). This phenomenon currently remains poorly studied: scientists do not know the nature of their occurrence, mechanisms of their formation in the higher layers of the atmosphere, etc. In this paper, the main characteristics of air showers and their main properties were described. The basis of EAS research was formed by numerous experiments conducted at the Hadron-55 complex facility. The showers were recorded and processed by the ionization calorimeter method, the essence of which is to measure energy of EAS by the total ionization of particles. This calorimeter with iron and lead absorbers provides high accuracy of measurements, the experimental error of which is only 10%. Experimental data were processed by a program developed in C++. The experiments showed that the most common particles in the showers are electron-photon and hadronic components. Their share is 95-98% of the particles in the shower stem. The program processed 8717 showers for a time equal to 485 hours. The analysis of the energy range of the registered showers confirms the theory of their galactic origin. Also, the calculation of the main parameters makes it possible to carry out further studies of EAS in a more extended form.

Effects of thunderstorms electric field on secondary photons of cosmic ray at large high altitude air shower observatory

Large high altitude air shower observatory (LHAASO) is a complex of extensive air shower (EAS) detector arrays, located on the Mt. Haizi (29°21' N, 100°08' E) at an altitude of 4410 m a. s. l., Daocheng, Sichuan Province, China. The information about primary cosmic rays can be obtained by using data from secondary particles measured at LHAASO, with photons make up the majority among these secondary particles. During thunderstorms, the atmospheric electric field can affect secondary charged particles (mainly positrons and electrons), thus changing the information of photons on the ground. In this work, Monte Carlo simulations are performed to investigate the effects of near-ground thunderstorm electric fields on cosmic ray secondary photons at LHAASO. A simple model with a vertical and uniform atmospheric electric field in a layer of atmosphere is used in our simulations. During thunderstorms, the number and energy of photons are found to significantly change and strongly depend on the electric field strength. In a field of –1000 V/cm (below the threshold of the relativistic runaway electron avalanche (RREA) process), the number of photons is increased by 23%. Also, the spectrum of photons softens, and the increased number of photons with energy less than 2 MeV exceeds 29%. In an electric field of –1700 V/cm (above the threshold of the RREA process), the number of photons experiences exponential growth, with an increase of 279%. The spectrum of photons becomes softer than that at –1000 V/cm, and the increased number with energy less than 2 MeV is more than 361%. It is consistent with the theory of RREA. For these phenomena of photons at LHAASO, the main factor is that the number of positrons and electrons are increased due to the acceleration of negative electric field on electrons, with increase of 65% in –1000 V/cm and 992% in –1700 V/cm, and the spectrum of positrons and electrons soften. Newborn free positrons/electrons may undergo bremsstrahlung and deposit part of their energy into photons, causing the change of number and energy of photons to follow roughly the same pattern as positrons and electrons. The simulation results can provide the information for understanding the variations of the data detected by LHAASO during thunderstorms and the acceleration mechanisms of secondary charged particles caused by an atmospheric electric field.

Methods of machine learning for the analysis of cosmic rays mass composition with the KASCADE experiment data

We study the problem of reconstruction of high-energy cosmic rays mass composition from the experimental data of extensive air showers. We develop several machine learning methods for the reconstruction of energy spectra of separate primary nuclei at energies 1–100 PeV, using the public data and Monte-Carlo simulations of the KASCADE experiment from the KCDC platform. We estimate the uncertainties of our methods, including the unfolding procedure, and show that the overall accuracy exceeds that of the method used in the original studies of the KASCADE experiment.

Proof of principle for template synthesis approach for the radio emission from vertical extensive air showers

The radio detection technique of cosmic ray air showers has gained renewed interest in the last two decades. While the radio experiments are very cost-effective to deploy, the Monte-Carlo simulations required to analyse the data are computationally expensive. Here we present a proof of concept for a novel way to synthesise the radio emission from extensive air showers in simulations. It is a hybrid approach which uses a single microscopic Monte-Carlo simulation, called the origin shower, to generate the radio emission from a target shower with a different longitudinal evolution, primary particle type and energy. The method employs semi-analytical relations which only depend on the shower parameters to transform the radio signals in the simulated antennas. We apply this method to vertical air showers with energies ranging from 1017eV to 1019eV and compare the results with CoREAS simulations in two frequency bands, namely the broad [20, 500] MHz band and a more narrow one at [30, 80] MHz. We gauge the synthesis quality using the maximal amplitude and energy fluence contained in the signal. We observe that the quality depends primarily on the difference in Xmax between the origin and target shower. After applying a linear bias correction, we find that for a shift in Xmax of less than 150 g/cm2, template synthesis has a bias of less than 2% and a scatter up to 6%, both in amplitude, on the broad frequency range. On the restricted [30, 80] MHz range the bias is similar, but the spread on amplitude drops down to 3%. These fluctuations are on the same level as the intrinsic scatter we observe in Monte-Carlo ensembles. We therefore surmise the observed scatter in amplitude to originate from intrinsic shower fluctuations we do not explicitly account for in template synthesis.

EAS Cherenkov LDF Analysis: CORSIKA Simulations for Tunka-133 and Chacaltaya Arrays

This study leverages CORSIKA simulations to analyze the Lateral Distribution Function (LDF) of Cherenkov photons in Extensive Air Showers (EAS) within the knee region of the cosmic ray energy spectrum (10^15 - 10^16 eV). Focusing on primary particles like helium, proton, oxygen, and iron nuclei at varying zenith angles (0˚, 20˚, and 45˚), we aimed to reconstruct Cherenkov photons' LDF using an Exponential Function model, tailored as a function of primary energy. Our approach involved a comparative analysis of simulated LDFs with experimental data from the Tunka-133 and Chacaltaya arrays. The results exhibit a high degree of concordance between simulated and observed data, affirming the validity of our method. We developed a set of approximation functions for different primary particles and zenith angles, enhancing our ability to identify the particle type in EAS events and accurately determine its energy. The primary contribution of our work lies in its potential to rapidly compile a comprehensive LDF pattern library, instrumental for analyzing real EAS array events and reconstructing the mass composition and primary cosmic ray energy spectrum. This advancement in CORSIKA-based simulation methods marks a significant stride in cosmic ray research, offering a robust tool for detailed EAS analysis.Highlights : Validation of Simulation Accuracy: Demonstrated high concordance between CORSIKA simulated LDFs and experimental data from Tunka-133 and Chacaltaya arrays. Enhanced Particle Identification: Development of approximation functions for various primary particles and zenith angles, improving accuracy in identifying particle types in EAS events. Advancement in Cosmic Ray Research: Potential to create a comprehensive LDF pattern library, significantly aiding in the reconstruction of mass composition and primary cosmic ray energy spectrum. Keywords : CORSIKA Simulations, Cherenkov LDF, EAS Analysis, Cosmic Ray Spectrum, Particle Identification