The hybrid energy spectrum of Telescope Array’s Middle Drum Detector and surface array

The Telescope Array experiment studies ultra high energy cosmic rays using a hybrid detector. Fluorescence telescopes measure the longitudinal development of the extensive air shower generated when a primary cosmic ray particle interacts with the atmosphere. Meanwhile, scintillator detectors measure the lateral distribution of secondary shower particles that hit the ground. The Middle Drum (MD) fluorescence telescope station consists of 14 telescopes from the High Resolution Fly’s Eye (HiRes) experiment, providing a direct link back to the HiRes measurements. Using the scintillator detector data in conjunction with the telescope data improves the geometrical reconstruction of the showers significantly, and hence, provides a more accurate reconstruction of the energy of the primary particle. The Middle Drum hybrid spectrum is presented and compared to that measured by the Middle Drum station in monocular mode. Further, the hybrid data establishes a link between the Middle Drum data and the surface array. A comparison between the Middle Drum hybrid energy spectrum and scintillator Surface Detector (SD) spectrum is also shown.     

Prospects for radio detection of ultra-high energy cosmic rays and neutrinos

The origin and nature of the highest energy cosmic ray events is currently the subject of intense investigation by giant air shower arrays and fluorescent detectors. These particles reach energies well beyond what can be achieved in ground-based particle accelerators and hence they are fundamental probes for particle physics as well as astrophysics. One of the main topics today focuses on the high energy end of the spectrum and the potential for the production of high-energy neutrinos. Above about 10²⁰ eV cosmic rays from extragalactic sources are expected to be severely attenuated by pion photoproduction interactions with photons of the cosmic microwave background. Investigating the shape of the cosmic ray spectrum near this predicted cut-off will be very important. In addition, a significant high-energy neutrino background is naturally expected as part of the pion decay chain which also contains much information.Because of the scarcity of these high-energy particles, larger and larger ground-based detectors have been built. The new generation of digital radio telescopes may play an important role in this, if properly designed. Radio detection of cosmic ray showers has a long history but was abandoned in the 1970s. Recent experimental developments together with sophisticated air shower simulations incorporating radio emission give a clearer understanding of the relationship between the air shower parameters and the radio signal, and have led to resurgence in its use. Observations of air showers by the SKA could, because of its large collecting area, contribute significantly to measuring the cosmic ray spectrum at the highest energies. Because of the large surface area of the moon, and the expected excellent angular resolution of the SKA, using the SKA to detect radio Cherenkov emission from neutrino-induced cascades in lunar regolith will be potentially the most important technique for investigating cosmic ray origin at energies above the photoproduction cut-off.     

A design study of atmospheric Cerenkov radiation telescope for very high energy gamma-ray astronomy

Detection of cosmic sources of very high energy gamma rays based on the atmospheric Cerenkov technique is discussed. Very high energy gamma-rays initiate, on entering the terrestrial atmosphere, electron-photon cascade showers with in turn produce Cerenkov photons in the air. Parabolic reflectors are used to focus these photons onto fast photomultipliers. Two methods of deployment of parabolic reflectors are in vogue: one in which all the reflectors are located close to each other in a compact array and the other in which the reflectors are spread out farther apart forming a distributed array. In the latter mode, the arrival direction of individual showers can be determined accurately by using the measured relative arrival times between different detectors. Detailed studies with the distributed array helped us to understand the various parameters in the two designs and evaluate their relative merits in reaching the ultimate goals of lowering the energy threshold and improving the signal to background ratio for the detection of gamma-ray sources. It is found that the relative superiority among the two types of arrays is a function of the exponent assumed for the differential power law energy spectrum for the gamma ray source. It is also seen that with the type of reflectors commonly used in atmospheric Cerenkov work, lower energy thresholds can be achieved with use of larger aperture.     

Evidence for the charge-excess contribution in air shower radio emission observed by the CODALEMA experiment

CODALEMA is one of the pioneer experiments dedicated to the radio detection of ultra high energy cosmic rays (UHECR), located at the radio observatory of Nançay (France). The CODALEMA experiment uses both a particle detector array and a radio antenna array. Data from both detection systems have been used to determine the ground coordinates of the core of extensive air showers (EAS). We discuss the observed systematic shift of the core positions determined with these two detection techniques. We show that this shift is due to the charge-excess contribution to the total radio emission of air showers, using the simulation code SELFAS. The dependences of the radio core shift to the primary cosmic ray characteristics are studied in details. The observation of this systematic shift can be considered as an experimental signature of the charge excess contribution.     

Instruments of RT-2 experiment onboard CORONAS-PHOTON and their test and evaluation IV: background simulations using GEANT-4 toolkit

Hard X-ray detectors in space are prone to background signals due to the ubiquitous cosmic rays and cosmic diffuse background radiation that continuously bombards the satellites which carry the detectors. In general, the background intensity depends on the space environment as well as the material surrounding the detectors. Understanding the behavior of the background noise in the detector is very important to extract the precise source information from the detector data. In this paper, we carry out Monte Carlo simulations using the GEANT-4 toolkit to estimate the prompt background noise measured with the detectors of the RT-2 Experiment onboard the CORONAS-PHOTON satellite.     

Simulations of a giant hybrid air shower detector

We have used Monte Carlo simulations to investigate the capabilities of a giant air shower observatory designed to detect showers initiated by cosmic rays with energies exceeding 10¹⁹ eV. The observatory is to consist of an array of detectors that will characterise the air shower at ground level, and optical detectors to measure the fluorescence light emitted by the shower in the atmosphere. Using these detectors together in a ‘hybrid’ configuration, we find that precise geometrical reconstruction of the shower axis is possible, leading to excellent resolution in energy. The technique is also shown to provide very good reconstruction below 10¹⁹ eV, at energies where the ground array is not fully efficient.     

Comparison of measured and simulated lateral distributions for electrons and muons with KASCADE

Lateral distributions for electrons and muons in extensive air showers measured with the array of the KASCADE experiment are compared to results of simulations based on the high-energy hadronic interaction models QGSJet and SIBYLL. It is shown, that the muon distributions are well described by both models. Deviations are found for the electromagnetic component, where both models predict a steeper lateral shape than observed in the data. For both models the observed lateral shapes of the electron component indicate a transition from a light to a more heavy composition of the cosmic ray spectrum above the knee.     

Detection of high energy showers by the Astroneu extensive air shower array

The Astroneu Extensive Air Shower (EAS) array comprises autonomous detection stations; each station consisting of 3 large scintillator detectors and one or more Radio Frequency (RF) antennas. The scintillator detectors of a station are able to detect showers with a low energy threshold of $20-30$ TeV, at a rate of $10-20$ showers per hour depending on the station geometrical layout. The RF antennas are used to detect very high energy EAS (E > 10¹⁷eV) through their radio wave emission. This work focuses in reconstructing and studying showers that have been detected synchronously by the scintillator detectors of two distant Astroneu stations. The performance of the array to detect and reconstruct the direction of such high energy (E > 5 · 10¹⁵eV) showers is evaluated by comparing the experimental measurements to the predictions of a detail Monte Carlo (MC) simulation.     

Capabilities of a giant hybrid air shower detector

Recent reports of superhigh energy cosmic rays beyond the expected spectral cutoff have intensified interest in the unknown origin of the highest energy cosmic rays. There is a need for a much larger data base of more precisely measured air showers. This requires new sensitive detectors of enormous aperture. Combining a ground array of particle counters with an optical detector of atmospheric fluorescence yields a detector of outstanding capability. Such a hybrid detector provides far more accurate measurements of energies, arrival directions, and primary particle atomic masses than can be achieved by either type of detector separately.     

Cosmic ray composition and energy spectrum from 1–30 PeV using the 40-string configuration of IceTop and IceCube

The mass composition of high energy cosmic rays depends on their production, acceleration, and propagation. The study of cosmic ray composition can therefore reveal hints of the origin of these particles. At the South Pole, the IceCube Neutrino Observatory is capable of measuring two components of cosmic ray air showers in coincidence: the electromagnetic component at high altitude (2835 m) using the IceTop surface array, and the muonic component above ∼1 TeV using the IceCube array. This unique detector arrangement provides an opportunity for precision measurements of the cosmic ray energy spectrum and composition in the region of the knee and beyond. We present the results of a neural network analysis technique to study the cosmic ray composition and the energy spectrum from 1 PeV to 30 PeV using data recorded using the 40-string/40-station configuration of the IceCube Neutrino Observatory.     

Composition of Cosmic Rays with <Emphasis Type="Italic">E</Emphasis><Superscript>0</Superscript> ⩾ 10<Superscript>17</Superscript> eV Based on Data from the Ground Scintillation Detectors of the Yakutsk EAS Array

The lateral distribution of cascade particles in extensive air showers from cosmic rays with energy E⁰ ? 10¹⁷ eV has been studied at the Yakutsk array by the ground scintillation detectors over the period of continuous observations 1977–2017. The experimental data are compared with those computed with various models for the development of extensive air showers from the CORSIKA software package. The best agreement between the theory and experiment is observed for the QGSjet-01-d model. In the energy range (1?20)× 10¹⁷ eV there is a change in the cosmic-raymass composition from 〈lnA〉 ≈ 2.5 to the proton one.     

Changes of the cosmic-ray mass composition in the 10–10 eV energy range

Data taken with ten Cosmic Ray Tracking (CRT) detectors and the HEGRA air-shower array on La Palma, Canary Islands, have been analysed to investigate changes of the cosmic ay mass composition at the ‘knee’ of the cosmic-ray flux spectrum near 10¹⁵ eV energy. The analysis is based on the angular distributions of particles in air showers. HEGRA data provided the shower size, direction, and core position and CRT data the particle track information. It is shown that the angular distribution of muons in air showers is sensitive to the composition over a wide range of shower sizes and, thus, primary cosmic-ray energies with little systematic uncertainties. Results can be easily expressed in terms of ln A of primary cosmic rays. In the lower part of the energy range covered, we have considerable overlap with direct composition measurements by the JACEE collaboration and find compatible results in the observed rise of ln A. Above about 10¹⁵ eV energy we find no or at most a slow further rise of ln A. Simple cosmic-ray composition models are presented which are fully consistent with our results as well as the JACEE flux and composition measurements and the flux measurements of the Tibet ASγ collaboration. Minimal three-parameter composition models defined by the same power-law slope of all elements below the knee and a common change in slope at a fixed rigidity are inconsistent with these data.     

Measurement of the flux of ultra high energy cosmic rays by the stereo technique

The High Resolution Fly’s Eye (HiRes) experiment has measured the flux of ultrahigh energy cosmic rays using the stereoscopic air fluorescence technique. The HiRes experiment consists of two detectors that observe cosmic ray showers via the fluorescence light they emit. HiRes data can be analyzed in monocular mode, where each detector is treated separately, or in stereoscopic mode where they are considered together. Using the monocular mode the HiRes collaboration measured the cosmic ray spectrum and made the first observation of the Greisen–Zatsepin–Kuzmin cutoff. In this paper we present the cosmic ray spectrum measured by the stereoscopic technique. Good agreement is found with the monocular spectrum in all details.     

Measurement of the cosmic-ray energy spectrum above 1016 eV with the LOFAR Radboud Air Shower Array

The energy reconstruction of extensive air showers measured with the LOFAR Radboud Air Shower Array (LORA) is presented in detail. LORA is a particle detector array located in the center of the LOFAR radio telescope in the Netherlands. The aim of this work is to provide an accurate and independent energy measurement for the air showers measured through their radio signal with the LOFAR antennas. The energy reconstruction is performed using a parameterized relation between the measured shower size and the cosmic-ray energy obtained from air shower simulations. In order to illustrate the capabilities of LORA, the all-particle cosmic-ray energy spectrum has been reconstructed, assuming that cosmic rays are composed only of protons or iron nuclei in the energy range between ∼2 × 10¹⁶ and 2 × 10¹⁸ eV. The results are compatible with literature values and a changing mass composition in the transition region from a Galactic to an extragalactic origin of cosmic rays.     

A MC approach to simulate up- and down-going neutrino showers including local topographic conditions

The extragalactic flux of protons is predicted to be suppressed above the famous Greisen–Zatsepin–Kuzmin cut-off at about E_GZK  50 EeV due to the resonant photo-pion production with the cosmic microwave background. Current cosmic ray data do not give a conclusive confirmation of the GZK cut-off and the quest about the origin and the chemical composition of the highest energy cosmic rays is still open. Amongst other particles neutrinos are expected to add to the composition of the cosmic radiation at highest energies. We present an approach to simulate neutrino induced air showers by a full Monte Carlo simulation chain. Starting with neutrinos at the top of the atmosphere, the performed simulations take into account the details of the neutrino propagation inside the Earth and atmosphere as well as inside the surrounding mountains. The products of the interactions are input for air shower simulations. The mountains are modelled by means of a digital elevation map. To exemplify the potential and features of the developed tools we study the possibility to detect neutrino induced extensive air showers with the fluorescence detector set-up of the Pierre Auger Observatory. Both, down-going neutrinos and up-going neutrinos are simulated and their rates are determined. To evaluate the sensitivity, as a function of the incoming direction, the aperture, the acceptance and the total observable event rates are calculated for the Waxman–Bahcall (WB) bound.     

Gamma-Hadron Separation using Čerenkov Photon Density Fluctuations

In the atmospheric Čerenkov technique γ-rays are detected against the abundant background produced by hadronic showers. In order to improve the signal to noise ratio of theexperiment, it is necessary to reject a significant fraction of hadronic showers. Traditional background rejection methods based on image shape parameters have been extensively used for the data from imaging telescopes. However, non-imaging Čerenkov telescopes have to develop very different means of statistically identifying and removing cosmic ray events. Some of the parameters, which could be potentially important for non-imaging arrays, are the temporal and spectral differences, the lateral distributions and density fluctuations of Čerenkov photons generated by γ-ray and hadron primaries. Here we study the differences in fluctuations of Čerenkov photon density in the light pool at the observation level from showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the PACT array has been used to evaluate the efficiency of the new technique. Various types of density fluctuations like the short range and medium range fluctuations as well as flatness parameter are studied. The estimated quality factors reflect the efficiencies with which the hadrons can be rejected from the data. Since some of these parameters are independent, the cuts may be applied in tandem and we demonstrate that the proton rejection efficiency of ∼90% can be achieved. Use of density fluctuations is particularly suited for wavefront sampling observations and it seems to be a good technique to improve the signal to noise ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

Characterization of the atmospheric muon flux in IceCube

Muons produced in atmospheric cosmic ray showers account for the by far dominant part of the event yield in large-volume underground particle detectors. The IceCube detector, with an instrumented volume of about a cubic kilometer, has the potential to conduct unique investigations on atmospheric muons by exploiting the large collection area and the possibility to track particles over a long distance. Through detailed reconstruction of energy deposition along the tracks, the characteristics of muon bundles can be quantified, and individual particles of exceptionally high energy identified. The data can then be used to constrain the cosmic ray primary flux and the contribution to atmospheric lepton fluxes from prompt decays of short-lived hadrons.In this paper, techniques for the extraction of physical measurements from atmospheric muon events are described and first results are presented. The multiplicity spectrum of TeV muons in cosmic ray air showers for primaries in the energy range from the knee to the ankle is derived and found to be consistent with recent results from surface detectors. The single muon energy spectrum is determined up to PeV energies and shows a clear indication for the emergence of a distinct spectral component from prompt decays of short-lived hadrons. The magnitude of the prompt flux, which should include a substantial contribution from light vector meson di-muon decays, is consistent with current theoretical predictions.The variety of measurements and high event statistics can also be exploited for the evaluation of systematic effects. In the course of this study, internal inconsistencies in the zenith angle distribution of events were found which indicate the presence of an unexplained effect outside the currently applied range of detector systematics. The underlying cause could be related to the hadronic interaction models used to describe muon production in air showers.     

Study of dispersion of mass distribution of ultra-high energy cosmic rays using a surface array of muon and electromagnetic detectors

We consider a hypothetical observatory of ultra-high energy cosmic rays consisting of two surface detector arrays that measure independently electromagnetic and muon signals induced by air showers. Using the constant intensity cut method, sets of events ordered according to each of both signal sizes are compared giving the number of matched events. Based on its dependence on the zenith angle, a parameter sensitive to the dispersion of the distribution of the logarithmic mass of cosmic rays is introduced. The results obtained using two post-LHC models of hadronic interactions are very similar and indicate a weak dependence on details of these interactions.     

A composition dependent energy scale and the determination of the cosmic ray primary mass in the ankle region

At present there are still several open questions about the origin of the ultra high energy cosmic rays. However, great progress in this area has been made in recent years due to the data collected by the present generation of ground based detectors like the Pierre Auger Observatory and Telescope Array. In particular, it is believed that the study of the composition of the cosmic rays as a function of energy can play a fundamental role for the understanding of the origin of the cosmic rays.The observatories belonging to this generation are composed of arrays of surface detectors and fluorescence telescopes. The duty cycle of the fluorescence telescopes is ∼10% in contrast with the ∼100% of the surface detectors. Therefore, the energy calibration of the events observed by the surface detectors is performed by using a calibration curve obtained from a set of high quality events observed in coincidence by both types of detectors. The advantage of this method is that the reconstructed energy of the events observed by the surface detectors becomes almost independent of simulations of the showers because just a small part of the reconstructed energy (the missing energy), obtained from the fluorescence telescopes, comes from simulations. However, the calibration curve obtained in this way depends on the composition of the cosmic rays, which can introduce biases in composition analyses when parameters with a strong dependence on primary energy are considered. In this work we develop an analytical method to study these effects. We consider AMIGA (Auger Muons and Infill for the Ground Array), the low energy extension of the Pierre Auger Observatory corresponding to the surface detectors, to illustrate the use of the method. In particular, we study the biases introduced by an energy calibration dependent on composition on the determination of the mean value of the number of muons, at a given distance to the showers axis, which is one of the parameters most sensitive to primary mass and has an almost linear dependence with primary energy.     

A new method to reconstruct the energy and determine the composition of cosmic rays from the measurement of Cherenkov light and particle densities in extensive air showers

A Monte Carlo study to reconstruct energy and mass of cosmic rays with energies above 300 TeV using ground based measurements of the electromagnetic part of showers initiated in the atmosphere is presented. The shower properties determined with two detector arrays measuring the air Cherenkov light and the particle densities as realized at the HEGRA experiment are processed to determine the energy of the primary particle without the need of any hypothesis concerning its mass. The mass of the primary particle is reconstructed coarsely from the same observables in parallel to the energy determination.