-
Atmospheric monitoring is an integral part of the design of the Cherenkov
Telescope Array (CTA), as atmospheric conditions affect the observations by
Imaging Atmospheric Cherenkov Telescopes (IACT) in multiple ways. The variable
optical properties of the atmosphere are a major contribution to the systematic
uncertainty in the determination of the energy and flux of the gamma photons.
Both the development of the air-shower and the production of Cherenkov light
depend on the molecular profile of the atmosphere. Additionally, the rapidly
changing aerosol profile, affecting the transmission of the Cherenkov light,
needs to be monitored on short time scales. Establishing a procedure to select
targets based on current atmospheric conditions can increase the efficiency of
the use of the observation time. The knowledge of atmospheric properties of the
future CTA locations and their annual and short-term variations in advance is
essential so that the atmospheric calibration can be readily applied to first
scientific data. To this end, some devices are already installed at one or both
of the selected sites...
-
Our study exploits the Constant Intensity Cut principles applied
simultaneously to muonic and electromagnetic detectors of cosmic rays. We use
the fact that the ordering of events according to their signal sizes induced in
different types of surface detectors provides information about the mass
composition of primary cosmic-ray beam, with low sensitivity to details of
hadronic interactions. Composition analysis at knee energies is performed using
Monte Carlo simulations for extensive air showers having maxima located far
away from a hypothetical observatory. Another type of a hypothetical
observatory is adopted to examine composition of ultra-high energy primaries
which initiate vertical air showers with maxima observed near surface
detectors.
-
Indications of a discrepancy between simulations and data on the number of
muons in cosmic ray (CR) showers exist over a large span of energies. We focus
in particular on the excess of multi-muon bundles observed by the DELPHI
detector at LEP and on the excess in the muon number in general reported by the
Pierre Auger Observatory. Even though the primary CR energies relevant for
these experiments differ by orders of magnitude, we can find a single mechanism
which can simultaneously increase predicted muon counts for both, while not
violating constraints from accelerators or from the longitudinal shower
development as observed by the Pierre Auger Observatory. We present a brief
theoretical motivation and describe a practical implementation of such a model,
based on the addition of soft particles to interactions above a chosen energy
threshold. Results of an extensive set of simulations show the behavior of this
model in various parts of a simplified parameter space.
-
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.
-
We focus on the primary composition of cosmic rays with the highest energies
that cause extensive air showers in the Earth's atmosphere. A way of examining
the two lowest order moments of the sample distribution of the depth of shower
maximum is presented. The aim is to show that useful information about the
composition of the primary beam can be inferred with limited knowledge we have
about processes underlying these observations. In order to describe how the
moments of the depth of shower maximum depend on the type of primary particles
and their energies, we utilize a superposition model. Using the principle of
maximum entropy, we are able to determine what trends in the primary
composition are consistent with the input data, while relying on a limited
amount of information from shower physics. Some capabilities and limitations of
the proposed method are discussed. In order to achieve a realistic description
of the primary mass composition, we pay special attention to the choice of the
parameters of the superposition model. We present two examples that demonstrate
what consequences can be drawn for energy dependent changes in the primary
composition.
-
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 introduce a simple branching model for the development of hadronic showers
in the Earth's atmosphere. Based on this model, we show how the size of the
pionic component followed by muons can be estimated. Several aspects of the
subsequent muonic component are also discussed. We focus on the energy
evolution of the muon production depth. We also estimate the impact of the
primary particle mass on the size of the hadronic component. Even though a
precise calculation of the development of air showers must be left to complex
Monte Carlo simulations, the proposed model can reveal qualitative insight into
the air shower physics.
-
The Cherenkov Telescope Array (CTA) is the next generation of ground-based
very high energy gamma-ray instruments and is planned to be built on two sites
(one in each hemisphere) in the coming years, with full array operation
foreseen to begin 2020. The goal of performing high precision gamma-ray energy
measurements while maximizing the use of observation time demands detailed and
fast information about atmospheric conditions. Besides LIDARs designed to
monitor clouds and aerosol content of the atmosphere in the pointing direction
of the CTA telescopes, we propose to use the "FRAM" (F(/Ph)otometric Robotic
Atmospheric Monitor) device, which is a small robotic astronomical telescope
with a large field of view and a sensitive CCD camera that together ensure
precise atmospheric characterization over the complete field-of-view of the
CTA. FRAM will use stellar photometry to measure atmospheric extinction across
the field of view of the CTA without interfering with the observation (unlike
laser-based methods). This allows FRAM to operate with high temporal resolution
and provide both real-time data for on-the-fly scheduling decisions and an
offline database for calibration and selection of scientific data. The fast
robotic mount of the telescope supports quick observation of multiple fields
when the array is split and even a check of the conditions in the directions of
the upcoming observations is possible. The FRAM concept is built upon
experience gained with a similar device operated at the Pierre Auger
Observatory. A working prototype of FRAM proposed for CTA is being built in
Prague for extensive testing before deployment on site; first results and
experiences with this prototype are presented.
-
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.
-
Many atmospheric and climatic criteria have to be taken into account for the
selection of a suitable site for the next generation of imaging air-shower
Cherenkov telescopes, the "Cherenkov Telescope Array" CTA. Such data are not
available with sufficient precision, thus a comparison of the proposed sites
and final decision based on a comprehensive characterization is impossible.
Identical cross-calibrated instruments have been developed which allow for
precise comparison between sites, the cross-validation of existing data, and
the ground-validation of satellite data. The site characterization work package
of the CTA consortium opted to construct and deploy 9 copies of an autonomous
multi-purpose weather sensor, incorporating an infrared cloud sensor, a newly
developed sensor for measuring the light of the night sky, and an
All-Sky-Camera, the whole referred to as Autonomous Tool for Measuring
Observatory Site COnditions PrEcisely (ATMOSCOPE). We present here the hardware
that was combined into the ATMOSCOPE and characterize its performance.
-
The All Sky Camera (ASC) was developed as an universal device for a
monitoring of the night sky quality and night sky background measurement. ASC
system consists of an astronomical CCD camera, a fish eye lens, a control
computer and associated electronics. The measurement is carried out during
astronomical twilight. The analysis results are the cloud fraction (the
percentage of the sky covered by clouds), night sky brightness (in mag/arcsec2)
and light background in the field of view of the camera. The analysis of the
cloud fraction is based on the astrometry (comparison to catalogue positions)
of the observed stars.
-
Contributions of the Pierre Auger Collaboration to the 33rd International
Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
-
The All Sky Camera (ASC) was developed as a universal device for the
monitoring of the night sky quality. Eight ASCs are already installed and
measure night sky parameters at eight of the candidate sites of the Cherenkov
Telescope Array (CTA) gamma-ray observatory. The ACS system consists of an
astronomical CCD camera, a fish eye lens, a control computer and associated
electronics. The measurement is carried out during astronomical night. The
images are automatically taken every 5 minutes and automatically processed
using the control computer of the device. The analysis results are the cloud
fraction (the percentage of the sky covered by clouds) and night sky brightness
(in mag/arcsec$^{2}$)
-
Inspired by recent astrophysical observations of leptonic excesses measured
by satellite experiments, we consider the impact of some general models of the
dark sector on the muon production in extensive air showers. We present a
compact approximative expression for the bremsstrahlung of a massive photon
from an electron and use it within Monte Carlo simulations to estimate the
amount of weakly interacting photon-like massive particles that could be
produced in an extensive air shower. We find that the resulting muon production
is by many orders of magnitude below the average muon count in a shower and
thus unobservable.
-
Using measured events from the fluorescence detector of the Pierre Auger
Observatory, an unbiased distribution of the atmospheric slant depths where
showers reach their maxima has been obtained. Analyzing the tail of this
distribution the proton-air cross-section for particle production at
center-of-mass energies per nucleon of 57 TeV is determined to be
[505$\pm$22(stat)+28$-$36(syst)] mb. Systematic uncertainties in the analysis
arise from the limited knowledge of the primary mass composition, the need to
use shower simulations and the selection of events. For the purpose of making
comparisons with accelerator data we also calculate the inelastic and total
proton-proton cross-sections using an extended Glauber model.
Our study exploits the Constant Intensity Cut principles applied simultaneously to muonic and electromagnetic detectors of cosmic rays. We use the fact that the ordering of events according to their signal sizes induced in different types of surface detectors provides information about the mass composition of primary cosmic-ray beam, with low sensitivity to details of hadronic interactions. Composition analysis at knee energies is performed using Monte Carlo simulations for extensive air showers having maxima located far away from a hypothetical observatory. Another type of a hypothetical observatory is adopted to examine composition of ultra-high energy primaries which initiate vertical air showers with maxima observed near surface detectors.
Indications of a discrepancy between simulations and data on the number of muons in cosmic ray (CR) showers exist over a large span of energies. We focus in particular on the excess of multi-muon bundles observed by the DELPHI detector at LEP and on the excess in the muon number in general reported by the Pierre Auger Observatory. Even though the primary CR energies relevant for these experiments differ by orders of magnitude, we can find a single mechanism which can simultaneously increase predicted muon counts for both, while not violating constraints from accelerators or from the longitudinal shower development as observed by the Pierre Auger Observatory. We present a brief theoretical motivation and describe a practical implementation of such a model, based on the addition of soft particles to interactions above a chosen energy threshold. Results of an extensive set of simulations show the behavior of this model in various parts of a simplified parameter space.
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.
We focus on the primary composition of cosmic rays with the highest energies that cause extensive air showers in the Earth's atmosphere. A way of examining the two lowest order moments of the sample distribution of the depth of shower maximum is presented. The aim is to show that useful information about the composition of the primary beam can be inferred with limited knowledge we have about processes underlying these observations. In order to describe how the moments of the depth of shower maximum depend on the type of primary particles and their energies, we utilize a superposition model. Using the principle of maximum entropy, we are able to determine what trends in the primary composition are consistent with the input data, while relying on a limited amount of information from shower physics. Some capabilities and limitations of the proposed method are discussed. In order to achieve a realistic description of the primary mass composition, we pay special attention to the choice of the parameters of the superposition model. We present two examples that demonstrate what consequences can be drawn for energy dependent changes in the primary composition.
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 introduce a simple branching model for the development of hadronic showers in the Earth's atmosphere. Based on this model, we show how the size of the pionic component followed by muons can be estimated. Several aspects of the subsequent muonic component are also discussed. We focus on the energy evolution of the muon production depth. We also estimate the impact of the primary particle mass on the size of the hadronic component. Even though a precise calculation of the development of air showers must be left to complex Monte Carlo simulations, the proposed model can reveal qualitative insight into the air shower physics.
The Cherenkov Telescope Array (CTA) is the next generation of ground-based very high energy gamma-ray instruments and is planned to be built on two sites (one in each hemisphere) in the coming years, with full array operation foreseen to begin 2020. The goal of performing high precision gamma-ray energy measurements while maximizing the use of observation time demands detailed and fast information about atmospheric conditions. Besides LIDARs designed to monitor clouds and aerosol content of the atmosphere in the pointing direction of the CTA telescopes, we propose to use the "FRAM" (F(/Ph)otometric Robotic Atmospheric Monitor) device, which is a small robotic astronomical telescope with a large field of view and a sensitive CCD camera that together ensure precise atmospheric characterization over the complete field-of-view of the CTA. FRAM will use stellar photometry to measure atmospheric extinction across the field of view of the CTA without interfering with the observation (unlike laser-based methods). This allows FRAM to operate with high temporal resolution and provide both real-time data for on-the-fly scheduling decisions and an offline database for calibration and selection of scientific data. The fast robotic mount of the telescope supports quick observation of multiple fields when the array is split and even a check of the conditions in the directions of the upcoming observations is possible. The FRAM concept is built upon experience gained with a similar device operated at the Pierre Auger Observatory. A working prototype of FRAM proposed for CTA is being built in Prague for extensive testing before deployment on site; first results and experiences with this prototype are presented.
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.
Many atmospheric and climatic criteria have to be taken into account for the selection of a suitable site for the next generation of imaging air-shower Cherenkov telescopes, the "Cherenkov Telescope Array" CTA. Such data are not available with sufficient precision, thus a comparison of the proposed sites and final decision based on a comprehensive characterization is impossible. Identical cross-calibrated instruments have been developed which allow for precise comparison between sites, the cross-validation of existing data, and the ground-validation of satellite data. The site characterization work package of the CTA consortium opted to construct and deploy 9 copies of an autonomous multi-purpose weather sensor, incorporating an infrared cloud sensor, a newly developed sensor for measuring the light of the night sky, and an All-Sky-Camera, the whole referred to as Autonomous Tool for Measuring Observatory Site COnditions PrEcisely (ATMOSCOPE). We present here the hardware that was combined into the ATMOSCOPE and characterize its performance.
The All Sky Camera (ASC) was developed as an universal device for a monitoring of the night sky quality and night sky background measurement. ASC system consists of an astronomical CCD camera, a fish eye lens, a control computer and associated electronics. The measurement is carried out during astronomical twilight. The analysis results are the cloud fraction (the percentage of the sky covered by clouds), night sky brightness (in mag/arcsec2) and light background in the field of view of the camera. The analysis of the cloud fraction is based on the astrometry (comparison to catalogue positions) of the observed stars.
The All Sky Camera (ASC) was developed as a universal device for the monitoring of the night sky quality. Eight ASCs are already installed and measure night sky parameters at eight of the candidate sites of the Cherenkov Telescope Array (CTA) gamma-ray observatory. The ACS system consists of an astronomical CCD camera, a fish eye lens, a control computer and associated electronics. The measurement is carried out during astronomical night. The images are automatically taken every 5 minutes and automatically processed using the control computer of the device. The analysis results are the cloud fraction (the percentage of the sky covered by clouds) and night sky brightness (in mag/arcsec$^{2}$)
Inspired by recent astrophysical observations of leptonic excesses measured by satellite experiments, we consider the impact of some general models of the dark sector on the muon production in extensive air showers. We present a compact approximative expression for the bremsstrahlung of a massive photon from an electron and use it within Monte Carlo simulations to estimate the amount of weakly interacting photon-like massive particles that could be produced in an extensive air shower. We find that the resulting muon production is by many orders of magnitude below the average muon count in a shower and thus unobservable.
Using measured events from the fluorescence detector of the Pierre Auger Observatory, an unbiased distribution of the atmospheric slant depths where showers reach their maxima has been obtained. Analyzing the tail of this distribution the proton-air cross-section for particle production at center-of-mass energies per nucleon of 57 TeV is determined to be [505$\pm$22(stat)+28$-$36(syst)] mb. Systematic uncertainties in the analysis arise from the limited knowledge of the primary mass composition, the need to use shower simulations and the selection of events. For the purpose of making comparisons with accelerator data we also calculate the inelastic and total proton-proton cross-sections using an extended Glauber model.