-
There has been a recent surge in interest in the detection of $\tau$
lepton-induced air showers from detectors at altitude. When a $\tau$ neutrino
($\nu_\tau$) enters the Earth it produces $\tau$ leptons as a result of nuclear
charged current interactions. In some cases, this process results in a $\tau$
lepton exiting the surface of the Earth, which can subsequently decay in the
atmosphere and produce an extensive air shower. These upward-going air showers
can be detected via fluorescence, optical Cherenkov, or geomagnetic radio
emission. Several experiments have been proposed to detect these signals. We
present a comprehensive simulation of the production of $\tau$ leptons by
$\nu_\tau$'s propagating through Earth to aid the design of future experiments.
These simulations for $\nu_\tau$'s and leptons in the energy range from
$10^{15}$ eV to $10^{21}$ eV treat the full range of incidence angles from
Earth-skimming to diametrically-traversing. Propagation of $\nu_\tau$'s and
leptons include the effects of rock and an ocean or ice layer of various
thicknesses. The interaction models include $\nu_\tau$ regeneration and account
for uncertainties in the Standard Model neutrino cross-section and in the
photo-nuclear contribution to the $\tau$ energy loss rate.
-
We investigate whether transition radiation from a particle shower crossing
the interface between Earth and air and induced by an Earth-skimming neutrino
can explain the upward event announced recently by the ANITA Collaboration.
While the properties of the observed signal can in principle be explained with
transition radiation, the flux necessary for a successful explanation is in
tension with the current best limits from the Pierre Auger Observatory, the
IceCube neutrino detector and the ANITA balloon. We also show that in this
scenario, the direction of the incoming neutrino is determined precisely to
within a few degrees, from the polarization properties of the observed events
and thanks to the Earth opacity to ultra high energy neutrinos.
-
Radio emission in atmospheric showers is currently interpreted in terms of
radiation due to the deviation of the charged particles in the magnetic field
of the Earth and to the charge excess (Askaryan effect). Each of these
mechanisms has a distinctive polarization. The complex signal patterns can be
qualitatively explained as the interference (superposition) of the fields
induced by each mechanism. In this work we explicitly and quantitatively test a
simple phenomenological model based on this idea. The model is constructed by
isolating each of the two components at the simulation level and by making use
of approximate symmetries for each of the contributions separately. The results
of the model are then checked against full ZHAireS Monte Carlo simulations of
the electric field calculated from first principles. We show that the simple
model describes radio emission at a few percent level in a wide range of
shower-observer geometries and on a shower-by-shower basis. As a consequence,
this approach provides a simple method to reduce the computing time needed to
accurately predict the electric field of radio pulses emitted from air showers,
with many practical applications in experimental situations of interest.
-
We measure the energy emitted by extensive air showers in the form of radio
emission in the frequency range from 30 to 80 MHz. Exploiting the accurate
energy scale of the Pierre Auger Observatory, we obtain a radiation energy of
15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV
arriving perpendicularly to a geomagnetic field of 0.24 G, scaling
quadratically with the cosmic-ray energy. A comparison with predictions from
state-of-the-art first-principle calculations shows agreement with our
measurement. The radiation energy provides direct access to the calorimetric
energy in the electromagnetic cascade of extensive air showers. Comparison with
our result thus allows the direct calibration of any cosmic-ray radio detector
against the well-established energy scale of the Pierre Auger Observatory.
-
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger
Observatory and is used to detect the radio emission of cosmic-ray air showers.
These observations are compared to the data of the surface detector stations of
the Observatory, which provide well-calibrated information on the cosmic-ray
energies and arrival directions. The response of the radio stations in the 30
to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of
the incoming electric field. For the latter, the energy deposit per area is
determined from the radio pulses at each observer position and is interpolated
using a two-dimensional function that takes into account signal asymmetries due
to interference between the geomagnetic and charge-excess emission components.
The spatial integral over the signal distribution gives a direct measurement of
the energy transferred from the primary cosmic ray into radio emission in the
AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air
shower arriving perpendicularly to the geomagnetic field. This radiation energy
-- corrected for geometrical effects -- is used as a cosmic-ray energy
estimator. Performing an absolute energy calibration against the
surface-detector information, we observe that this radio-energy estimator
scales quadratically with the cosmic-ray energy as expected for coherent
emission. We find an energy resolution of the radio reconstruction of 22% for
the data set and 17% for a high-quality subset containing only events with at
least five radio stations with signal.
-
Coherent radiation at radio frequencies from high-energy showers fully
contained in a dense radio-transparent medium - like ice, salt or regolith -
has been extensively investigated as a promising technique to search for
ultra-high energy (UHE) neutrinos. Additional emission in the form of
transition radiation may occur when a neutrino-induced shower produced close to
the Earth surface emerges from the ground into atmospheric air. We present the
first detailed evaluation of transition radiation from high-energy showers
crossing the boundary between two different media. We found that transition
radiation is sizable over a wide solid angle and coherent up to $\sim$ 1 GHz.
These properties encourage further work to evaluate the potential of a
large-aperture UHE neutrino experiment based on detection of transition
radiation.
-
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are
detectable with the Surface Detector array of the Pierre Auger Observatory. We
report here on searches through Auger data from 1 January 2004 until 20 June
2013. No neutrino candidates were found, yielding a limit to the diffuse flux
of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound
predictions. Neutrino identification is attempted using the broad
time-structure of the signals expected in the SD stations, and is efficiently
done for neutrinos of all flavors interacting in the atmosphere at large zenith
angles, as well as for "Earth-skimming" neutrino interactions in the case of
tau neutrinos. In this paper the searches for downward-going neutrinos in the
zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for
upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L.
single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an
$E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times
10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~
s^{-1}~ sr^{-1}}$.
-
A measurement of the cosmic-ray spectrum for energies exceeding
$4{\times}10^{18}$ eV is presented, which is based on the analysis of showers
with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger
Observatory between 1 January 2004 and 31 December 2013. The measured spectrum
confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$
eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with
index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by
a smooth suppression region. For the energy ($E_\text{s}$) at which the
spectral flux has fallen to one-half of its extrapolated value in the absence
of suppression, we find
$E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$
eV.
-
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.
-
Contributions of the Pierre Auger Collaboration to the 33rd International
Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
-
The detection of electromagnetic pulses from high energy showers is used as a
means to search for Ultra-High Energy cosmic ray and neutrino interactions. An
approximate formula has been obtained to numerically evaluate the radio pulse
emitted by a charged particle that instantaneously accelerates, moves at
constant speed along a straight track and halts again instantaneously. The
approximate solution is applied to the particle track after dividing it in
smaller subintervals. The resulting algorithm (often referred to as the ZHS
algorithm) is also the basis for most of the simulations of the electric field
produced in high energy showers in dense media. In this work, the
electromagnetic pulses as predicted with the ZHS algorithm are compared to
those obtained with an exact solution of the electric field produced by a
charged particle track. The precise conditions that must apply for the
algorithm to be valid are discussed and its accuracy is addressed. This
comparison is also made for electromagnetic showers in dense media. The ZHS
algorithm is shown to describe Cherenkov radiation and to be valid for most
situations of interest concerning detectors searching for Ultra-High Energy
neutrinos. The results of this work are also relevant for the simulation of
pulses emitted from air showers.
-
Using detailed Monte Carlo simulations we have characterized the features of
the radio emission of inclined air showers in the Ultra-High Frequency band
(300 MHz - 3 GHz). The Fourier-spectrum of the radiation is shown to have a
sizable intensity well into the GHz frequency range. The emission is mainly due
to transverse currents induced by the geomagnetic field and to the excess
charge produced by the Askaryan effect. At these frequencies only a
significantly reduced volume of the shower around the axis contributes
coherently to the signal observed on the ground. The size of the coherently
emitting volume depends on frequency, shower geometry and observer position,
and is interpreted in terms of the relative time delays. At ground level, the
maximum emission at high frequencies is concentrated in an elliptical ring-like
region around the intersection of a Cherenkov cone with its vertex at shower
maximum and the ground. The frequency spectrum of inclined showers when
observed at positions that view shower maximum in the Cherenkov direction, is
shown to be in broad agreement with the pulses detected by the Antarctic
Impulsive Transient Antenna (ANITA) experiment, making the interpretation that
they are due to Ultra-High Energy Cosmic Ray atmospheric showers consistent
with our simulations. These results are also of great importance for
experiments aiming to detect molecular bremsstrahlung radiation in the GHz
range as they present an important background for its detection.
-
We present predictions for the radio pulses emitted by extensive air showers
using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full
complexity of ultra-high energy cosmic-ray induced shower development in the
atmosphere, and allows the calculation of the electric field in both the time
and frequency domains. We do not presuppose any emission mechanism, but our
results are compatible with a superposition of geomagnetic and charge excess
radio emission effects. We investigate the polarization of the electric field
as well as the effects of the refractive index n and shower geometry on the
radio pulses. We show that geometry, coupled to the relativistic effects that
appear when using a realistic n>1, play a prominent role on the radio emission
of air showers.
-
An in-depth characterization of coherent radio Cherenkov pulses from particle
showers in dense dielectric media, referred to as the Askaryan effect, is
presented. The time-domain calculation developed in this article is based on a
form factor to account for the lateral dimensions of the shower. It is
computationally efficient and able to reproduce the results of detailed
particle shower simulations with high fidelity in most regions of practical
interest including Fresnel effects due to the longitudinal development of the
shower. In addition, an intuitive interpretation of the characteristics of the
Askaryan pulse is provided. We expect our approach to benefit the analysis of
radio pulses in experiments exploiting the radio technique.
-
The Cherenkov radio pulse emitted by hadronic showers in ice is calculated
for showers of energies in the EeV range. This is obtained with 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 and optimized 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 pions and muons, is found to be below 1%. 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.
-
The electric field of the Cherenkov radio pulse produced by a single charged
particle track in a dielectric medium is derived from first principles. An
algorithm is developed to obtain the pulse in the time domain for numerical
calculations. The algorithm is implemented in a Monte Carlo simulation of
electromagnetic showers in dense media (specifically designed for coherent
radio emission applications) as might be induced by interactions of ultra-high
energy neutrinos. The coherent Cherenkov radio emission produced by such
showers is obtained simultaneously both in the time and frequency domains. A
consistency check performed by Fourier-transforming the pulse in time and
comparing it to the frequency spectrum obtained directly in the simulations
yields, as expected, fully consistent results. The reversal of the time
structure inside the Cherenkov cone and the signs of the corresponding pulses
are addressed in detail. The results, besides testing algorithms used for
reference calculations in the frequency domain, shed new light into the
properties of the radio pulse in the time domain. The shape of the pulse in the
time domain is directly related to the depth development of the excess charge
in the shower and its width to the observation angle with respect to the
Cherenkov direction. This information can be of great practical importance for
interpreting actual data.
-
We study the frequency and angular dependences of Cherenkov radio pulses
originated by the excess of electrons in electromagnetic showers in different
dense media. We develop a simple model to relate the main characteristics of
the electric field spectrum to the properties of the shower such as
longitudinal and lateral development. This model allows us to establish the
scaling of the electric field spectrum with the properties of the medium such
as density, radiation length, Moliere radius, critical energy and refraction
index. We normalize the predictions of the scaling relations to the numerical
results obtained in our own developed GEANT4-based Monte Carlo simulation, and
we give a unified parameterization of the frequency spectrum and angular
distribution of the electric field in ice, salt, and the lunar regolith, in
terms of the relevant properties of the media. Our parameterizations are valid
for electromagnetic showers below the energy at which the
Landau-Pomeranchuk-Migdal effect starts to be relevant in these media. They
also provide an approximate estimate of radio emission in hadronic showers
induced by high energy cosmic rays or neutrinos.
-
The possibilities of detecting high energy neutrinos through inclined showers
produced in the atmosphere are addressed with an emphasis on the detection of
air showers by arrays of particle detectors. Rates of inclined showers produced
by both down-going neutrino interactions and by up-coming $\tau$ decays from
earth-skimming neutrinos as a function of shower energy are calculated with
analytical methods using two sample neutrino fluxes with different spectral
indices. The relative contributions from different flavors and charged, neutral
current and resonant interactions are compared for down-going neutrinos
interacting in the atmosphere. No detailed description of detectors is
attempted but rough energy thresholds are implemented to establish the ranges
of energies which are more suitable for neutrino detection through inclined
showers. Down-going and up-coming rates are compared.
-
In this article I review the main theoretical problems that are posed by the
highest energy end of the observed cosmic ray spectrum, stressing the
importance of establishing their composition in order to decide between
proposed scenarios. I then discuss the possibilities that are opened by the
detection of inclined showers with extensive air shower arrays. Recent progress
in modelling magnetic deviations for these showers has allowed the analysis of
inclined showers that were detected by the Haverah Park experiment. This
analysis disfavours models that predict a large proportion of photons in the
highest energy cosmic rays and open up new possibilities for future shower
array detectors particularly those, like the Pierre Auger Observatory, using
water Cherenkov detectors.
-
In this article I briefly discuss the characteristics of the Auger
observatories paying particular attention to the role of inclined showers, both
in the search for high energy neutrino interactions deep in the atmosphere and
as an alternative tool for the study of cosmic rays, particulartly their
composition.
-
Here an attempt is made at summarizing the presentations, most of which were
about the highest energy particles observed in nature. Particular attention is
paid to the solutions to the Ultra High Energy Cosmic Ray particles, to the new
and forthcoming data, to the new proposals for experiments and to the role of
primary composition, that were amongst the most discussed subtopics.
-
The role of Horizontal Showers induced by cosmic rays is discussed in detail.
A new approach to the calculation of the muon component in horizontal air
showers induced by protons, heavier nuclei or photons is presented which allows
a simple analytical evaluation of the muon density profiles at ground level.
The results of the first application of these results to horizontal air showers
detected at the Haverah Park Array by the recently started {\sl Leeds-Santiago
collaboration}, leading to important restrictions on composition at ultra high
energies, are reported.
-
We study the potential of the Pierre Auger detector for horizontal air
showers initiated by ultra high energy neutrino. Assuming some simple trigger
requirements we obtain measurable event rates for neutrino fluxes from AGN,
from topological defects and from the interactions of cosmic rays with the
microwave background.
-
The emission of particles from black holes created in the early Universe has
detectable astrophysical consequences. The most stringent bound on their
abundance has been obtained from the absence of a detectable diffuse flux of
100 MeV photons. Further scrutiny of these bounds is of interest as they, for
instance, rule out primordial black holes as a dark matter candidate. We here
point out that these bounds can, in principle, be improved by studying the
diffuse cosmic neutrino flux. Measurements of near-vertical atmospheric
neutrino fluxes in a region of low geomagnetic latitude can provide a
competitive bound. The most favorable energy to detect a possible diffuse flux
of primordial black hole origin is found to be a few MeV. We also show that
measurements of the diffuse $\nu _\tau$ flux is the most promising to improve
the existing bounds deduced from gamma-ray measurements. Neutrinos from
individual black hole explosions can be detected in the GeV-TeV energy region.
We find that the kilometer-scale detectors, recently proposed, are able to
establish competitive bounds.
There has been a recent surge in interest in the detection of $\tau$ lepton-induced air showers from detectors at altitude. When a $\tau$ neutrino ($\nu_\tau$) enters the Earth it produces $\tau$ leptons as a result of nuclear charged current interactions. In some cases, this process results in a $\tau$ lepton exiting the surface of the Earth, which can subsequently decay in the atmosphere and produce an extensive air shower. These upward-going air showers can be detected via fluorescence, optical Cherenkov, or geomagnetic radio emission. Several experiments have been proposed to detect these signals. We present a comprehensive simulation of the production of $\tau$ leptons by $\nu_\tau$'s propagating through Earth to aid the design of future experiments. These simulations for $\nu_\tau$'s and leptons in the energy range from $10^{15}$ eV to $10^{21}$ eV treat the full range of incidence angles from Earth-skimming to diametrically-traversing. Propagation of $\nu_\tau$'s and leptons include the effects of rock and an ocean or ice layer of various thicknesses. The interaction models include $\nu_\tau$ regeneration and account for uncertainties in the Standard Model neutrino cross-section and in the photo-nuclear contribution to the $\tau$ energy loss rate.
We investigate whether transition radiation from a particle shower crossing the interface between Earth and air and induced by an Earth-skimming neutrino can explain the upward event announced recently by the ANITA Collaboration. While the properties of the observed signal can in principle be explained with transition radiation, the flux necessary for a successful explanation is in tension with the current best limits from the Pierre Auger Observatory, the IceCube neutrino detector and the ANITA balloon. We also show that in this scenario, the direction of the incoming neutrino is determined precisely to within a few degrees, from the polarization properties of the observed events and thanks to the Earth opacity to ultra high energy neutrinos.
Radio emission in atmospheric showers is currently interpreted in terms of radiation due to the deviation of the charged particles in the magnetic field of the Earth and to the charge excess (Askaryan effect). Each of these mechanisms has a distinctive polarization. The complex signal patterns can be qualitatively explained as the interference (superposition) of the fields induced by each mechanism. In this work we explicitly and quantitatively test a simple phenomenological model based on this idea. The model is constructed by isolating each of the two components at the simulation level and by making use of approximate symmetries for each of the contributions separately. The results of the model are then checked against full ZHAireS Monte Carlo simulations of the electric field calculated from first principles. We show that the simple model describes radio emission at a few percent level in a wide range of shower-observer geometries and on a shower-by-shower basis. As a consequence, this approach provides a simple method to reduce the computing time needed to accurately predict the electric field of radio pulses emitted from air showers, with many practical applications in experimental situations of interest.
We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.
Coherent radiation at radio frequencies from high-energy showers fully contained in a dense radio-transparent medium - like ice, salt or regolith - has been extensively investigated as a promising technique to search for ultra-high energy (UHE) neutrinos. Additional emission in the form of transition radiation may occur when a neutrino-induced shower produced close to the Earth surface emerges from the ground into atmospheric air. We present the first detailed evaluation of transition radiation from high-energy showers crossing the boundary between two different media. We found that transition radiation is sizable over a wide solid angle and coherent up to $\sim$ 1 GHz. These properties encourage further work to evaluate the potential of a large-aperture UHE neutrino experiment based on detection of transition radiation.
Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
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.
The detection of electromagnetic pulses from high energy showers is used as a means to search for Ultra-High Energy cosmic ray and neutrino interactions. An approximate formula has been obtained to numerically evaluate the radio pulse emitted by a charged particle that instantaneously accelerates, moves at constant speed along a straight track and halts again instantaneously. The approximate solution is applied to the particle track after dividing it in smaller subintervals. The resulting algorithm (often referred to as the ZHS algorithm) is also the basis for most of the simulations of the electric field produced in high energy showers in dense media. In this work, the electromagnetic pulses as predicted with the ZHS algorithm are compared to those obtained with an exact solution of the electric field produced by a charged particle track. The precise conditions that must apply for the algorithm to be valid are discussed and its accuracy is addressed. This comparison is also made for electromagnetic showers in dense media. The ZHS algorithm is shown to describe Cherenkov radiation and to be valid for most situations of interest concerning detectors searching for Ultra-High Energy neutrinos. The results of this work are also relevant for the simulation of pulses emitted from air showers.
Using detailed Monte Carlo simulations we have characterized the features of the radio emission of inclined air showers in the Ultra-High Frequency band (300 MHz - 3 GHz). The Fourier-spectrum of the radiation is shown to have a sizable intensity well into the GHz frequency range. The emission is mainly due to transverse currents induced by the geomagnetic field and to the excess charge produced by the Askaryan effect. At these frequencies only a significantly reduced volume of the shower around the axis contributes coherently to the signal observed on the ground. The size of the coherently emitting volume depends on frequency, shower geometry and observer position, and is interpreted in terms of the relative time delays. At ground level, the maximum emission at high frequencies is concentrated in an elliptical ring-like region around the intersection of a Cherenkov cone with its vertex at shower maximum and the ground. The frequency spectrum of inclined showers when observed at positions that view shower maximum in the Cherenkov direction, is shown to be in broad agreement with the pulses detected by the Antarctic Impulsive Transient Antenna (ANITA) experiment, making the interpretation that they are due to Ultra-High Energy Cosmic Ray atmospheric showers consistent with our simulations. These results are also of great importance for experiments aiming to detect molecular bremsstrahlung radiation in the GHz range as they present an important background for its detection.
We present predictions for the radio pulses emitted by extensive air showers using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full complexity of ultra-high energy cosmic-ray induced shower development in the atmosphere, and allows the calculation of the electric field in both the time and frequency domains. We do not presuppose any emission mechanism, but our results are compatible with a superposition of geomagnetic and charge excess radio emission effects. We investigate the polarization of the electric field as well as the effects of the refractive index n and shower geometry on the radio pulses. We show that geometry, coupled to the relativistic effects that appear when using a realistic n>1, play a prominent role on the radio emission of air showers.
An in-depth characterization of coherent radio Cherenkov pulses from particle showers in dense dielectric media, referred to as the Askaryan effect, is presented. The time-domain calculation developed in this article is based on a form factor to account for the lateral dimensions of the shower. It is computationally efficient and able to reproduce the results of detailed particle shower simulations with high fidelity in most regions of practical interest including Fresnel effects due to the longitudinal development of the shower. In addition, an intuitive interpretation of the characteristics of the Askaryan pulse is provided. We expect our approach to benefit the analysis of radio pulses in experiments exploiting the radio technique.
The Cherenkov radio pulse emitted by hadronic showers in ice is calculated for showers of energies in the EeV range. This is obtained with 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 and optimized 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 pions and muons, is found to be below 1%. 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.
The electric field of the Cherenkov radio pulse produced by a single charged particle track in a dielectric medium is derived from first principles. An algorithm is developed to obtain the pulse in the time domain for numerical calculations. The algorithm is implemented in a Monte Carlo simulation of electromagnetic showers in dense media (specifically designed for coherent radio emission applications) as might be induced by interactions of ultra-high energy neutrinos. The coherent Cherenkov radio emission produced by such showers is obtained simultaneously both in the time and frequency domains. A consistency check performed by Fourier-transforming the pulse in time and comparing it to the frequency spectrum obtained directly in the simulations yields, as expected, fully consistent results. The reversal of the time structure inside the Cherenkov cone and the signs of the corresponding pulses are addressed in detail. The results, besides testing algorithms used for reference calculations in the frequency domain, shed new light into the properties of the radio pulse in the time domain. The shape of the pulse in the time domain is directly related to the depth development of the excess charge in the shower and its width to the observation angle with respect to the Cherenkov direction. This information can be of great practical importance for interpreting actual data.