Cherenkov τ shower earth-skimming method for PeV–EeV ντ observation with Ashra

We describe a method of observation for PeV–EeV τ neutrinos using Cherenkov light from the air showers of decayed τs produced by τ neutrino interactions in the Earth. Aiming for the realization of neutrino astronomy utilizing the Earth-skimming τ neutrino detection technique, highly precise determination of arrival direction is key due to the following issues: (1) clear identification of neutrinos by identifying those vertices originating within the Earth’s surface and (2) identification of very high energy neutrino sources. The Ashra detector uses newly developed light collectors which realize both a 42°-diameter field-of-view and arcminute resolution. Therefore, it has superior angular resolution for imaging Cherenkov air showers. In this paper, we estimate the sensitivity of and cosmic-ray background resulting from application of the Ashra-1 Cherenkov τ shower observation method. Both data from a commissioning run and a long-term observation (with fully equipped trigger system and one light collector) are presented. Our estimates are based on a detailed Monte Carlo simulation which describes all relevant shower processes from neutrino interaction to Cherenkov photon detection produced by τ air showers. In addition, the potential to determine the arrival direction of Cherenkov showers is evaluated by using the maximum likelihood method. We conclude that the Ashra-1 detector is a unique probe into detection of very high energy neutrinos and their accelerators.     

Searching for tau neutrinos with Cherenkov telescopes

Cherenkov telescopes have the capability of detecting high energy tau neutrinos in the energy range of 1–1000 PeV by searching for very inclined showers. If a tau lepton, produced by a tau neutrino, escapes from the Earth or a mountain, it will decay and initiate a shower in the air which can be detected by an air shower fluorescence or Cherenkov telescope. In this paper, we present detailed Monte Carlo simulations of corresponding event rates for the VERITAS and two proposed Cherenkov Telescope Array sites: Meteor Crater and Yavapai Ranch, which use representative AGN neutrino flux models and take into account topographic conditions of the detector sites. The calculated neutrino sensitivities depend on the observation time and the shape of the energy spectrum, but in some cases are comparable or even better than corresponding neutrino sensitivities of the IceCube detector. For VERITAS and the considered Cherenkov Telescope Array sites the expected neutrino sensitivities are up to factor 3 higher than for the MAGIC site because of the presence of surrounding mountains.     

Simulations of reflected radio signals from cosmic ray induced air showers

We present the calculation of coherent radio pulses emitted by extensive air showers induced by ultra-high energy cosmic rays accounting for reflection on the Earth’s surface. Results have been obtained with a simulation program that calculates the contributions from shower particles after reflection at a surface plane. The properties of the radiation are discussed in detail emphasizing the effects of reflection. The shape of the frequency spectrum is shown to be closely related to the angle of the observer with respect to shower axis, becoming hardest in the Cherenkov direction. The intensity of the flux at a fixed observation angle is shown to scale with the square of the primary particle energy to very good accuracy indicating the coherent aspect of the emission. The simulation methods of this paper provide the foundations for energy reconstruction of experiments looking at the Earth from balloons and satellites. They can also be used in dedicated studies of existing and future experimental proposals.     

Thinned simulations of extremely energetic showers in dense media for radio applications

This paper presents our results for full 3-D simulations of very-high to ultra-high energy electromagnetic cascades – and the associated coherent Cherenkov radiation – as might be produced by high-energy neutrino interactions in dense media. Using “thinning” techniques, we develop an algorithm based on the existing “ZHS” code, and demonstrate that the new “ZHS-thinned” code can produce fast and accurate results for showers up to . Using ZHS-thinned, we develop new parameterisations for the radiation from showers in ice, salt, and the lunar regolith, with a separate treatment of the megaregolith (deep regolith). Our parameterisations include for the first time a method to simulate fluctuations in shower length induced by the LPM effect. Our results, which avoid the pit-falls of scaling simulations from lower energies, allow improved calculations of the detection probability for experiments searching for high-energy neutrinos using the radio technique.     

The ExaVolt Antenna: A large-aperture, balloon-embedded antenna for ultra-high energy particle detection

We describe the scientific motivation, experimental basis, design methodology, and simulated performance of the ExaVolt Antenna (EVA) mission, and planned ultra-high energy (UHE) particle observatory under development for NASA’s suborbital super-pressure balloon program in Antarctica. EVA will improve over ANITA’s integrated totals – the current state-of-the-art in UHE suborbital payloads – by 1–2 orders of magnitude in a single flight. The design is based on a novel application of toroidal reflector optics which utilizes a super-pressure balloon surface, along with a feed-array mounted on an inner membrane, to create an ultra-large radio antenna system with a synoptic view of the Antarctic ice sheet below it. Radio impulses arise via the Askaryan effect when UHE neutrinos interact within the ice, or via geosynchrotron emission when UHE cosmic rays interact in the atmosphere above the continent. EVA’s instantaneous antenna aperture is estimated to be several hundred m² for detection of these events within a 150–600 MHz band. For standard cosmogenic UHE neutrino models, EVA should detect of order 30 events per flight in the EeV energy regime. For UHE cosmic rays, of order 15,000 geosynchrotron events would be detected in total, several hundred above 10 EeV, and of order 60 above the GZK cutoff energy.     

Coherent Cherenkov radio pulses from hadronic showers up to EeV energies

The Cherenkov radio pulse emitted by hadronic showers of energies in the EeV range in ice is calculated for the first time using full three dimensional simulations of both shower development and the coherent radio pulse emitted as the excess charge develops in the shower. A Monte Carlo, ZHAireS, has been developed for this purpose combining the high energy hadronic interaction capabilities of AIRES, and the dense media propagation capabilities of TIERRAS, with the precise low energy tracking and specific algorithms developed to calculate the radio emission in ZHS. A thinning technique is implemented to allow the simulation of radio pulses induced by showers up to 10 EeV in ice. The code is validated comparing the results for electromagnetic and hadronic showers to those obtained with GEANT4 and ZHS codes. The contribution to the pulse of other shower particles in addition to electrons and positrons, mainly protons, pions and muons, is found to be below 3% for 10 PeV and above proton induced showers. The characteristics of hadronic showers and the corresponding Cherenkov frequency spectra are compared with those from purely electromagnetic showers. The dependence of the spectra on shower energy and high-energy hadronic model is addressed and parameterizations for the radio emission in hadronic showers in ice are given for practical applications.     

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.     

Cross calibration of telescope optical throughput efficiencies using reconstructed shower energies for the Cherenkov Telescope Array

For reliable event reconstruction of Imaging Atmospheric Cherenkov Telescopes (IACTs), calibration of the optical throughput efficiency is required. Within current facilities, this is achieved through the use of ring shaped images generated by muons. Here, a complementary approach is explored, achieving cross calibration of elements of IACT arrays through pairwise comparisons between telescopes, focussing on its applicability to the upcoming Cherenkov Telescope Array (CTA). Intercalibration of telescopes of a particular type using eventwise comparisons of shower image amplitudes has previously been demonstrated to recover the relative telescope optical responses. A method utilising the reconstructed energy as an alternative to image amplitude is presented, enabling cross calibration between telescopes of varying types within an IACT array. Monte Carlo studies for two plausible CTA layouts have shown that this calibration procedure recovers the relative telescope response efficiencies at the few per cent level.     

A Monte Carlo template based analysis for air-Cherenkov arrays

We present a high-performance event reconstruction algorithm: an Image Pixel-wise fit for Atmospheric Cherenkov Telescopes (ImPACT). The reconstruction algorithm is based around the likelihood fitting of camera pixel amplitudes to an expected image template. A maximum likelihood fit is performed to find the best-fit shower parameters. A related reconstruction algorithm has already been shown to provide significant improvements over traditional reconstruction for both the CAT and H.E.S.S. experiments. We demonstrate a significant improvement to the template generation step of the procedure, by the use of a full Monte Carlo air shower simulation in combination with a ray-tracing optics simulation to more accurately model the expected camera images. This reconstruction step is combined with an MVA-based background rejection.Examples are shown of the performance of the ImPACT analysis on both simulated and measured (from a strong VHE source) gamma-ray data from the H.E.S.S. array, demonstrating an improvement in sensitivity of more than a factor two in observation time over traditional image moments-fitting methods, with comparable performance to previous likelihood fitting analyses. ImPACT is a particularly promising approach for future large arrays such as the Cherenkov Telescope Array (CTA) due to its improved high-energy performance and suitability for arrays of mixed telescope types.     

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an Extensive Air Shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10¹⁷ and 10¹⁹ eV and zenith angles up to 65°. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte Carlo results showing how LTP functions from data are in good agreement with simulations.     

Improved energy resolution for VHE gamma ray astronomy with systems of Cherenkov telescopes

We present analysis techniques to improve the energy resolution of stereoscopic systems of imaging atmospheric Cherenkov telescopes, using the HEGRA telescope system as an example. The techniques include (i) the determination of the height of the shower maximum, which is then taken into account in the energy determination, and (ii) the determination of the location of the shower core with the additional constraint that the direction of the gamma rays is known a priori. This constraint can be applied for gamma ray point sources, and results in a significant improvement in the localization of the shower core, which translates into better energy resolution. Combining both techniques, the HEGRA telescopes reach an energy resolution between 9% and 12%, over the entire energy range from 1 TeV to almost 100 TeV. Options for further improvements of the energy resolution are discussed.     

Radio Cherenkov signals from the Moon: Neutrinos and cosmic rays

Neutrino production of radio Cherenkov signals in the Moon is the object of radio telescope observations. Depending on the energy range and detection parameters, the dominant contribution to the neutrino signal may come from interactions of the neutrino on the Moon facing the telescope, rather than neutrinos that have traversed a portion of the Moon. Using the approximate analytic expression of the effective lunar aperture from a recent paper by Gayley, Mutel and Jaeger, we evaluate the background from cosmic ray interactions in the lunar regolith. We also consider the modifications to the effective lunar aperture from generic non-standard model neutrino interactions. A background to neutrino signals are radio Cherenkov signals from cosmic ray interactions. For cosmogenic neutrino fluxes, neutrino signals will be difficult to observe because of low neutrino flux at the high energy end and large cosmic ray background in the lower energy range considered here. We show that lunar radio detection of neutrino interactions is best suited to constrain or measure neutrinos from astrophysical sources and probe non-standard neutrino-nucleon interactions such as microscopic black hole production.     

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.     

Performance of a scintillation detector array operated with LHAASO-KM2A electronics

A scintillation detector array composed of 115 detectors and covering an area of about 20000 m² was installed at the end of 2016 at the Yangbajing international cosmic ray observatory and has been taking data since then. The array is equipped with electronics from Large High Altitude Air Shower Observatory Square Kilometer Complex Array (LHAASO-KM2A) and, in turn, currently serves as the largest debugging and testing platform for the LHAASO-KM2A. Furthermore, the array was used to study the performance of a wide field-of-view air Cherenkov telescope by providing accurate information on the shower core, direction and energy, etc. This work is mainly dealing with the scintillation detector array. The experimental setup and the offline calibration are described in detail. Then, a thorough comparison between the data and Monte Carlo (MC) simulations is presented and a good agreement is obtained. With the even-odd method, the resolutions of the shower direction and core are measured. Finally, successful observations of the expected Moon’s and Sun’s shadows of cosmic rays (CRs) verify the measured angular resolution.     

The center of lateral iso-density contours for inclined cosmic air showers

The horizontal lateral density of a cosmic air shower with a non-zero zenith angle is asymmetric. The asymmetry consist of a stretching of the iso-density contours to ellipses and to a shift of the center of the elliptic contours with respect to the core of the shower. The shift is caused by atmospheric attenuation. The modeling of the attenuation results in an equation for the shift as a function of zenith angle and the size of the iso-density contours. A more accurate equation is obtained by investigating the shift in lateral densities of simulated showers. It is shown how the shift can be incorporated in an elliptic lateral density function. A linear approximation for the shift allows for an analytical solution for the shifted elliptic density. Its predictions for the polar variations of the density are compared with data of simulated showers.     

The air shower maximum probed by Cherenkov effects from radio emission

Radio detection of cosmic-ray-induced air showers has come to a flight the last decade. Along with the experimental efforts, several theoretical models were developed. The main radio-emission mechanisms are established to be the geomagnetic emission due to deflection of electrons and positrons in Earth’s magnetic field and the charge-excess emission due to a net electron excess in the air shower front. It was only recently shown that Cherenkov effects play an important role in the radio emission from air showers. In this article we show the importance of these effects to extract quantitatively the position of the shower maximum from the radio signal, which is a sensitive measure for the mass of the initial cosmic ray. We also show that the relative magnitude of the charge-excess and geomagnetic emission changes considerably at small observer distances where Cherenkov effects apply.     

Measurement of the cosmic ray composition at the knee with the SPASE-2/AMANDA-B10 detectors

The mass composition of high-energy cosmic rays at energies above 10¹⁵ eV can provide crucial information for the understanding of their origin. Air showers were measured simultaneously with the SPASE-2 air shower array and the AMANDA-B10 Cherenkov telescope at the South Pole. This combination has the advantage to sample almost all high-energy shower muons and is thus a new approach to the determination of the cosmic ray composition. The change in the cosmic ray mass composition was measured versus existing data from direct measurements at low energies. Our data show an increase of the mean log atomic mass lnA by about 0.8 between 500 TeV and 5 PeV. This trend of an increasing mass through the “knee” region is robust against a variety of systematic effects.     

Upgrading and testing the 3D reconstruction of gamma-ray air showers as observed with an array of Imaging Atmospheric Cherenkov Telescopes

Stereoscopic arrays of Imaging Atmospheric Cherenkov Telescopes allow to reconstruct gamma-ray-induced showers in three dimensions, which offers several advantages: direct access to the shower parameters in space and straightforward calorimetric measurement of the incident energy. In addition, correlations between the different images of the same shower are taken into account. An analysis method based on a simple 3D-model of electromagnetic showers was recently implemented in the framework of the H.E.S.S. experiment. In the present article, the method is completed by an additional quality criterion, which reduces the background contamination by a factor of about 2 in the case of extended sources, while keeping gamma-ray efficiency at a high level. On the other hand, the dramatic flares of the blazar PKS 2155-304 in July 2006, which provided H.E.S.S. data with an almost pure gamma-ray sample, offered the unique opportunity of a precision test of the 3D-reconstruction method as well as of the H.E.S.S. simulations used in its calibration. An agreement at a few percent level is found between data and simulations for the distributions of all 3D shower parameters.     

Upgrading and testing the 3D reconstruction of gamma-ray air showers as observed with an array of Imaging Atmospheric Cherenkov Telescopes

Stereoscopic arrays of Imaging Atmospheric Cherenkov Telescopes allow to reconstruct gamma-ray-induced showers in three dimensions, which offers several advantages: direct access to the shower parameters in space and straightforward calorimetric measurement of the incident energy. In addition, correlations between the different images of the same shower are taken into account. An analysis method based on a simple 3D-model of electromagnetic showers was recently implemented in the framework of the H.E.S.S. experiment. In the present article, the method is completed by an additional quality criterion, which reduces the background contamination by a factor of about 2 in the case of extended sources, while keeping gamma-ray efficiency at a high level. On the other hand, the dramatic flares of the blazar PKS 2155-304 in July 2006, which provided H.E.S.S. data with an almost pure gamma-ray sample, offered the unique opportunity of a precision test of the 3D-reconstruction method as well as of the H.E.S.S. simulations used in its calibration. An agreement at a few percent level is found between data and simulations for the distributions of all 3D shower parameters.