Source author record

T. Huege

T. Huege appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

124works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

124 published item(s)

preprint2026arXiv

Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observatory

We present a novel approach for assessing the muon content of air showers with large zenith angles on a combined analysis of their radio emission and particle footprint. We use the radiation energy reconstructed by the Auger Engineering Radio Array (AERA) as an energy estimator and determine the muon number independently with the water-Cherenkov detector array of the Pierre Auger Observatory, deployed on a 1500 m grid. We focus our analysis on air showers with primary energy above 4 EeV to ensure full detection efficiency. Over approximately ten years of accumulated data, we identify a set of 40 high-quality events that are used in the analysis. The estimated muon contents in data are compatible with those for iron primaries as predicted by current-generation hadronic interaction models. This result can be interpreted as a deficit of muons in simulations as a lighter mass composition has been established from Xmax measurements. This muon deficit was already observed in previous analyses of the Auger Collaboration and is confirmed using hybrid events that include radio measurements for the first time.

preprint2023arXiv

Results from recent analysis of KASCADE-Grande data

KASCADE and its extension array of KASCADE-Grande were devoted to measure individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The experiment has substantially contributed to investigate the energy spectrum and mass composition of cosmic rays in the transition region from galactic to extragalactic origin of cosmic rays as well as to quantify the characteristics of hadronic interaction models in the air shower development through validity tests using the multi-detector information from KASCADE-Grande. Although the data accumulation was completed in 2013, data analysis is still continuing. Recently, we investigated the reliability of the new hadronic interactions models of the SIBYLL version 2.3d only with the energy spectra from the KASCADE-Grande data. The evolution of the muon content of high energy air showers in the atmosphere is studied as well, using EPOS-LHC, SIBYLL 2.3, QGSJET-II-04 and SIBYLL 2.3c. In this talk, recent results from KASCADE-Grande and the update of the KASCADE Cosmic Ray Data Centre (KCDC) will be discussed.

preprint2022arXiv

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV.

preprint2022arXiv

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models.

preprint2022arXiv

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

The LOPES experiment was a radio interferometer built at the existing air shower array KASCADE-Grande in Karlsruhe, Germany. The last configuration of LOPES was called LOPES 3D and consisted of ten tripole antennas. Each of these antennas consisted of three crossed dipoles east-west, north-south, and vertically aligned. With this, LOPES 3D had the unique possibility to study the benefits of measurements with vertically aligned antennas in the environment of the well understood and calibrated particle detector array KASCADE-Grande. The measurements with three spatially coincident antennas allows a redundant reconstruction of the electric field vector. Several methods to exploit the redundancy were developed and tested. Furthermore, for the first time in LOPES, the background noise could be studied polarization- and direction dependent. With LOPES 3D it could be demonstrated that radio detection reaches a higher efficiency for inclined showers when including measurements with vertically aligned antennas and that the vertical component gets more important for the measurement of inclined showers. In this contribution we discuss a weighting scheme for the best combination of three redundant reconstructed electric field vectors. Furthermore, we discuss the influence of these weighting schemes on the ability to reconstruct air showers using the radio method. We show an estimate of the radio efficiency for inclined showers with focus on the benefits of measurements with vertically aligned antennas and we present the direction dependent noise in the different polarizations.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2022arXiv

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $δ_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $δ_{γ,0} > -10^{-21}$, $δ_{γ,1} > -10^{-40}$ eV$^{-1}$ and $δ_{γ,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $δ_{\mathrm{had},0} < 10^{-19}$, $δ_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $δ_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$σ$ CL.

preprint2022arXiv

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.

preprint2022arXiv

The SLAC T-510 experiment for radio emission from particle showers: detailed simulation study and interpretation

Over the last several decades, radio detection of air showers has been widely used to detect ultra-high-energy cosmic rays. We developed an experiment under controlled laboratory conditions at SLAC with which we measured the radio-frequency radiation from a charged particle shower produced by bunches of electrons as primaries with known energy. The shower took place in a target made of High Density Polyethylene located in a strong magnetic field. The experiment was designed so that Askaryan and magnetically-induced components of the radio emission could be measured independently. At the same time, we performed a detailed simulation of this experiment to predict the radio signal using two microscopic formalisms, endpoint and ZHS. In this paper, we present the simulation scheme and make a comparison with data characteristics such as linearity with magnetic field and amplitude. The simulations agree with the measurements within uncertainties and present a good description of the data. In particular, reflections within the target that accounted for the largest systematic uncertainties are addressed. The prediction of the amplitude of Askaryan emission agrees with measurements to within 5% for the endpoint formalism and 11% for the ZHS formalism. The amplitudes of magnetically-induced emission agree to within 5% for the endpoint formalism and less than 1% for the ZHS formalism. The agreement of the absolute scale of emission gives confidence in state-of-the-art air shower simulations which are based on the applied formalisms.

preprint2021arXiv

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector.

preprint2021arXiv

Final results of the LOPES radio interferometer for cosmic-ray air showers

LOPES, the LOFAR prototype station, was an antenna array for cosmic-ray air showers operating from 2003 - 2013 within the KASCADE-Grande experiment. Meanwhile, the analysis is finished and the data of air-shower events measured by LOPES are available with open access in the KASCADE Cosmic Ray Data Center (KCDC). This article intends to provide a summary of the achievements, results, and lessons learned from LOPES. By digital, interferometric beamforming the detection of air showers became possible in the radio-loud environment of the Karlsruhe Institute of Technology (KIT). As a prototype experiment, LOPES tested several antenna types, array configurations and calibration techniques, and pioneered analysis methods for the reconstruction of the most important shower parameters, i.e., the arrival direction, the energy, and mass-dependent observables such as the position of the shower maximum. In addition to a review and update of previously published results, we also present new results based on end-to-end simulations including all known instrumental properties. For this, we applied the detector response to radio signals simulated with the CoREAS extension of CORSIKA, and analyzed them in the same way as measured data. Thus, we were able to study the detector performance more accurately than before, including some previously inaccessible features such as the impact of noise on the interferometric cross-correlation beam. These results led to several improvements, which are documented in this paper and can provide useful input for the design of future cosmic-ray experiments based on the digital radio-detection technique.

preprint2020arXiv

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$σ$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $α_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

preprint2020arXiv

Radio emission from negative lightning leader steps reveals inner meter-scale structure

We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ μ$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models.

preprint2020arXiv

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere.

preprint2020arXiv

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

preprint2020arXiv

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

preprint2020arXiv

The SKA Particle Array Prototype: The First Particle Detector at the Murchison Radio-astronomy Observatory

We report on the design, deployment, and first results from a scintillation detector deployed at the Murchison Radio-astronomy Observatory (MRO). The detector is a prototype for a larger array -- the Square Kilometre Array Particle Array (SKAPA) -- planned to allow the radio-detection of cosmic rays with the Murchison Widefield Array and the low-frequency component of the Square Kilometre Array. The prototype design has been driven by stringent limits on radio emissions at the MRO, and to ensure survivability in a desert environment. Using data taken from Nov.\ 2018 to Feb.\ 2019, we characterize the detector response while accounting for the effects of temperature fluctuations, and calibrate the sensitivity of the prototype detector to through-going muons. This verifies the feasibility of cosmic ray detection at the MRO. We then estimate the required parameters of a planned array of eight such detectors to be used to trigger radio observations by the Murchison Widefield Array.

preprint2019arXiv

Combined sensitivity to the neutrino mass ordering with JUNO, the IceCube Upgrade, and PINGU

The ordering of the neutrino mass eigenstates is one of the fundamental open questions in neutrino physics. While current-generation neutrino oscillation experiments are able to produce moderate indications on this ordering, upcoming experiments of the next generation aim to provide conclusive evidence. In this paper we study the combined performance of the two future multi-purpose neutrino oscillation experiments JUNO and the IceCube Upgrade, which employ two very distinct and complementary routes towards the neutrino mass ordering. The approach pursued by the $20\,\mathrm{kt}$ medium-baseline reactor neutrino experiment JUNO consists of a careful investigation of the energy spectrum of oscillated $\barν_e$ produced by ten nuclear reactor cores. The IceCube Upgrade, on the other hand, which consists of seven additional densely instrumented strings deployed in the center of IceCube DeepCore, will observe large numbers of atmospheric neutrinos that have undergone oscillations affected by Earth matter. In a joint fit with both approaches, tension occurs between their preferred mass-squared differences $ Δm_{31}^{2}=m_{3}^{2}-m_{1}^{2} $ within the wrong mass ordering. In the case of JUNO and the IceCube Upgrade, this allows to exclude the wrong ordering at $>5σ$ on a timescale of 3--7 years --- even under circumstances that are unfavorable to the experiments' individual sensitivities. For PINGU, a 26-string detector array designed as a potential low-energy extension to IceCube, the inverted ordering could be excluded within 1.5 years (3 years for the normal ordering) in a joint analysis.

preprint2019arXiv

Design and Performance of the first IceAct Demonstrator at the South Pole

In this paper we describe the first results of a compact imaging air-Cherenkov telescope, IceAct, operating in coincidence with the IceCube Neutrino Observatory (IceCube) at the geographic South Pole. An array of IceAct telescopes (referred to as the IceAct project) is under consideration as part of the IceCube-Gen2 extension to IceCube. Surface detectors in general will be a powerful tool in IceCube-Gen2 for distinguishing astrophysical neutrinos from the dominant backgrounds of cosmic-ray induced atmospheric muons and neutrinos: the IceTop array is already in place as part of IceCube, but has a high energy threshold. Although the duty cycle will be lower for the IceAct telescopes than the present IceTop tanks, the IceAct telescopes may prove to be more effective at lowering the detection threshold for air showers. Additionally, small imaging air-Cherenkov telescopes in combination with IceTop, the deep IceCube detector or other future detector systems might improve measurements of the composition of the cosmic ray energy spectrum. In this paper we present measurements of a first 7-pixel imaging air Cherenkov telescope demonstrator, proving the capability of this technology to measure air showers at the South Pole in coincidence with IceTop and the deep IceCube detector.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

preprint2016arXiv

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

preprint2016arXiv

Accelerator measurements of magnetically-induced radio emission from particle cascades with applications to cosmic-ray air showers

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

preprint2016arXiv

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec θ)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec θ)_\mathrm{max}$.

preprint2016arXiv

Energy and Flux Measurements of Ultra-High Energy Cosmic Rays Observed During the First ANITA Flight

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.

preprint2016arXiv

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

preprint2016arXiv

Nanosecond-level time synchronization of autonomous radio detector stations for extensive air showers

To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used for cross-checks that indeed we reach nanosecond-scale timing accuracy by this correction. First, we operate a "beacon transmitter" which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.

preprint2016arXiv

Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory

Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to extend its range of detection and to directly measure the muon content of the particle showers. It consists of an infill of surface water-Cherenkov detectors accompanied by buried scintillator detectors used for muon counting. The main objectives of the AMIGA engineering array, referred to as the Unitary Cell, are to identify and resolve all engineering issues as well as to understand the muon-number counting uncertainties related to the design of the detector. The mechanical design, fabrication and deployment processes of the muon counters of the Unitary Cell are described in this document. These muon counters modules comprise sealed PVC casings containing plastic scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier tubes, and acquisition electronics. The modules are buried approximately 2.25 m below ground level in order to minimize contamination from electromagnetic shower particles. The mechanical setup, which allows access to the electronics for maintenance, is also described in addition to tests of the modules' response and integrity. The completed Unitary Cell has measured a number of air showers of which a first analysis of a sample event is included here.

preprint2016arXiv

Radio detection of cosmic ray air showers in the digital era

In 1965 it was discovered that cosmic ray air showers emit impulsive radio signals at frequencies below 100 MHz. After a period of intense research in the 1960s and 1970s, however, interest in the detection technique faded almost completely. With the availability of powerful digital signal processing techniques, new attempts at measuring cosmic ray air showers via their radio emission were started at the beginning of the new millennium. Starting with modest, small-scale digital prototype setups, the field has evolved, matured and grown very significantly in the past decade. Today's second-generation digital radio detection experiments consist of up to hundreds of radio antennas or cover areas of up to 17 km$^{2}$. We understand the physics of the radio emission in extensive air showers in detail and have developed analysis strategies to accurately derive from radio signals parameters which are related to the astrophysics of the primary cosmic ray particles, in particular their energy, arrival direction and estimators for their mass. In parallel to these successes, limitations inherent in the physics of the radio signals have also become increasingly clear. In this article, we review the progress of the past decade and the current state of the field, discuss the current paradigm of the radio emission physics and present the experimental evidence supporting it. Finally, we discuss the potential for future applications of the radio detection technique to advance the field of cosmic ray physics.

preprint2016arXiv

Radio measurements of the energy and the depth of the shower maximum of cosmic-ray air showers by Tunka-Rex

We reconstructed the energy and the position of the shower maximum of air showers with energies $E \gtrsim 100 $PeV applying a method using radio measurements performed with Tunka-Rex. An event-to-event comparison to air-Cherenkov measurements of the same air showers with the Tunka-133 photomultiplier array confirms that the radio reconstruction works reliably. The Tunka-Rex reconstruction methods and absolute scales have been tuned on CoREAS simulations and yield energy and $X_{\mathrm{max}}$ values consistent with the Tunka-133 measurements. The results of two independent measurement seasons agree within statistical uncertainties, which gives additional confidence in the radio reconstruction. The energy precision of Tunka-Rex is comparable to the Tunka-133 precision of $15 %$, and exhibits a $20 %$ uncertainty on the absolute scale dominated by the amplitude calibration of the antennas. For $X_{\mathrm{max}}$, this is the first direct experimental correlation of radio measurements with a different, established method. At the moment, the $X_{\mathrm{max}}$ resolution of Tunka-Rex is approximately $40 $g/cm$^2$. This resolution can probably be improved by deploying additional antennas and by further development of the reconstruction methods, since the present analysis does not yet reveal any principle limitations.

preprint2016arXiv

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

preprint2016arXiv

Search for ultrarelativistic magnetic monopoles with the Pierre Auger Observatory

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $γ=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $γ=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

preprint2016arXiv

Testing Hadronic Interactions at Ultrahigh Energies with Air Showers Measured by the Pierre Auger Observatory

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

preprint2016arXiv

The Pierre Auger Observatory Upgrade - Preliminary Design Report

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

preprint2015arXiv

Calibrating the absolute amplitude scale for air showers measured at LOFAR

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

preprint2015arXiv

Calibration of the absolute amplitude scale of the Tunka Radio Extension (ICRC 2015)

The Tunka Radio Extension (Tunka-Rex) is an array of 44 radio antenna stations, distributed over 3 km$^{2}$, constituting a radio detector for air showers with an energy threshold around 10$^{17}$ eV. It is an extension to Tunka-133, an air-Cherenkov detector in Siberia, which is used as an external trigger for Tunka-Rex and provides a reliable reconstruction of energy and shower maximum. Each antenna station consists of two perpendicularly aligned active antennas, called SALLAs. An antenna calibration of the SALLA with a commercial reference source enables us to reconstruct the detected radio signal on an absolute scale. Since the same reference source was used for the calibration of LOPES and, in a calibration campaign in 2014, also for LOFAR, these three experiments now have a consistent calibration and, therefore, absolute scale. This was a key ingredient to resolve a longer standing contradiction between measurements of two calibrated experiments. We will present how the calibration was performed and compare radio measurements of air showers from Tunka-Rex to model calculations with the radio simulation code CoREAS, confirming it within the scale uncertainty of the calibration of 18%.

preprint2015arXiv

Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations

Cosmic ray air showers emit radio pulses at MHz frequencies, which can be measured with radio antenna arrays - like LOPES at the Karlsruhe Institute of Technology in Germany. To improve the understanding of the radio emission, we test theoretical descriptions with measured data. The observables used for these tests are the absolute amplitude of the radio signal, and the shape of the radio lateral distribution. We compare lateral distributions of more than 500 LOPES events with two recent and public Monte Carlo simulation codes, REAS 3.11 and CoREAS (v 1.0). The absolute radio amplitudes predicted by REAS 3.11 are in good agreement with the LOPES measurements. The amplitudes predicted by CoREAS are lower by a factor of two, and marginally compatible with the LOPES measurements within the systematic scale uncertainties. In contrast to any previous versions of REAS, REAS 3.11 and CoREAS now reproduce the shape of the measured lateral distributions correctly. This reflects a remarkable progress compared to the situation a few years ago, and it seems that the main processes for the radio emission of air showers are now understood: The emission is mainly due to the geomagnetic deflection of the electrons and positrons in the shower. Less important but not negligible is the Askaryan effect (net charge variation). Moreover, we confirm that the refractive index of the air plays an important role, since it changes the coherence conditions for the emission: Only the new simulations including the refractive index can reproduce rising lateral distributions which we observe in a few LOPES events. Finally, we show that the lateral distribution is sensitive to the energy and the mass of the primary cosmic ray particles.

preprint2015arXiv

Digital radio detection of cosmic rays: achievements, status and perspectives

Over the past decade, radio detection of cosmic rays has matured from small-scale prototype experiments to installations spanning several km$^2$ with more than a hundred antennas. The physics of the radio signal is well understood and simulations and measurements are in good agreement. We have learned how to extract important cosmic ray parameters such as the geometry of the air shower and the energy of the primary particle from the radio signal, and have developed very promising approaches to also determine the mass of the primary particles. At the same time, limitations have become increasingly clear. I review the progress made in the past decade and provide a personal view on further potential for future development.

preprint2015arXiv

High-precision measurements of extensive air showers with the SKA

As of 2023, the Square Kilometre Array will constitute the world's largest radio telescope, offering unprecedented capabilities for a diverse science programme in radio astronomy. At the same time, the SKA will be ideally suited to detect extensive air showers initiated by cosmic rays in the Earth's atmosphere via their radio emission. With its very dense and uniform antenna spacing in a fiducial area of one km$^2$ and its large bandwidth of 50-350 MHz, the low-frequency part of the SKA will provide very precise measurements of individual cosmic ray air showers. These precision measurements will allow detailed studies of the mass composition of cosmic rays in the energy region of transition from a Galactic to an extragalactic origin. Also, the SKA will facilitate three-dimensional "tomography" of the electromagnetic cascades of air showers, allowing the study of particle interactions at energies beyond the reach of the LHC. Finally, studies of possible connections between air showers and lightning initiation can be taken to a new level with the SKA. We discuss the science potential of air shower detection with the SKA and report on the technical requirements and project status.

preprint2015arXiv

Improved absolute calibration of LOPES measurements and its impact on the comparison with REAS 3.11 and CoREAS simulations

LOPES was a digital antenna array detecting the radio emission of cosmic-ray air showers. The calibration of the absolute amplitude scale of the measurements was done using an external, commercial reference source, which emits a frequency comb with defined amplitudes. Recently, we obtained improved reference values by the manufacturer of the reference source, which significantly changed the absolute calibration of LOPES. We reanalyzed previously published LOPES measurements, studying the impact of the changed calibration. The main effect is an overall decrease of the LOPES amplitude scale by a factor of $2.6 \pm 0.2$, affecting all previously published values for measurements of the electric-field strength. This results in a major change in the conclusion of the paper 'Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations' published in Astroparticle Physics 50-52 (2013) 76-91: With the revised calibration, LOPES measurements now are compatible with CoREAS simulations, but in tension with REAS 3.11 simulations. Since CoREAS is the latest version of the simulation code incorporating the current state of knowledge on the radio emission of air showers, this new result indicates that the absolute amplitude prediction of current simulations now is in agreement with experimental data.

preprint2015arXiv

Investigation of the radio wavefront of air showers with LOPES measurements and CoREAS simulations (ARENA 2014)

We investigated the radio wavefront of cosmic-ray air showers with LOPES measurements and CoREAS simulations: the wavefront is of approximately hyperbolic shape and its steepness is sensitive to the shower maximum. For this study we used 316 events with an energy above 0.1 EeV and zenith angles below $45^\circ$ measured by the LOPES experiment. LOPES was a digital radio interferometer consisting of up to 30 antennas on an area of approximately 200 m x 200 m at an altitude of 110 m above sea level. Triggered by KASCADE-Grande, LOPES measured the radio emission between 43 and 74 MHz, and our analysis might strictly hold only for such conditions. Moreover, we used CoREAS simulations made for each event, which show much clearer results than the measurements suffering from high background. A detailed description of our result is available in our recent paper published in JCAP09(2014)025. The present proceeding contains a summary and focuses on some additional aspects, e.g., the asymmetry of the wavefront: According to the CoREAS simulations the wavefront is slightly asymmetric, but on a much weaker level than the lateral distribution of the radio amplitude.

preprint2015arXiv

KCDC - The KASCADE Cosmic-ray Data Centre

KCDC, the KASCADE Cosmic-ray Data Centre, is a web portal, where data of astroparticle physics experiments will be made available for the interested public. The KASCADE experiment, financed by public money, was a large-area detector for the measurement of high-energy cosmic rays via the detection of air showers. KASCADE and its extension KASCADE-Grande stopped finally the active data acquisition of all its components including the radio EAS experiment LOPES end of 2012 after more than 20 years of data taking. In a first release, with KCDC we provide to the public the measured and reconstructed parameters of more than 160 million air showers. In addition, KCDC provides the conceptional design, how the data can be treated and processed so that they are also usable outside the community of experts in the research field. Detailed educational examples make a use also possible for high-school students and early stage researchers.

preprint2015arXiv

Measurement of cosmic-ray air showers with the Tunka Radio Extension (Tunka-Rex)

Tunka-Rex is a radio detector for cosmic-ray air showers in Siberia, triggered by Tunka-133, a co-located air-Cherenkov detector. The main goal of Tunka-Rex is the cross-calibration of the two detectors by measuring the air-Cherenkov light and the radio signal emitted by the same air showers. This way we can explore the precision of the radio-detection technique, especially for the reconstruction of the primary energy and the depth of the shower maximum. The latter is sensitive to the mass of the primary cosmic-ray particles. In this paper we describe the detector setup and explain how electronics and antennas have been calibrated. The analysis of data of the first season proves the detection of cosmic-ray air showers and therefore, the functionality of the detector. We confirm the expected dependence of the detection threshold on the geomagnetic angle and the correlation between the energy of the primary cosmic-ray particle and the radio amplitude. Furthermore, we compare reconstructed amplitudes of radio pulses with predictions from CoREAS simulations, finding agreement within the uncertainties.

preprint2015arXiv

Precision measurements of cosmic ray air showers with the SKA

Supplemented with suitable buffering techniques, the low-frequency part of the SKA can be used as an ultra-precise detector for cosmic-ray air showers at very high energies. This would enable a wealth of scientific applications: the physics of the transition from Galactic to extragalactic cosmic rays could be probed with very high precision mass measurements, hadronic interactions could be studied up to energies well beyond the reach of man-made particle accelerators, air shower tomography could be performed with very high spatial resolution exploiting the large instantaneous bandwidth and very uniform instantaneous $u$-$v$ coverage of SKA1-LOW, and the physics of thunderstorms and possible connections between cosmic rays and lightning initiation could be studied in unprecedented levels of detail. In this article, we describe the potential of the SKA as an air shower radio detector from the perspective of existing radio detection efforts and discuss the associated technical requirements.

preprint2015arXiv

Revised absolute amplitude calibration of the LOPES experiment

One of the main aims of the LOPES experiment was the evaluation of the absolute amplitude of the radio signal of air showers. This is of special interest since the radio technique offers the possibility for an independent and highly precise determination of the energy scale of cosmic rays on the basis of signal predictions from Monte Carlo simulations. For the calibration of the amplitude measured by LOPES we used an external source. Previous comparisons of LOPES measurements and simulations of the radio signal amplitude predicted by CoREAS revealed a discrepancy of the order of a factor of two. A re-measurement of the reference calibration source, now performed for the free field, was recently performed by the manufacturer. The updated calibration values lead to a lowering of the reconstructed electric field measured by LOPES by a factor of $2.6 \pm 0.2$ and therefore to a significantly better agreement with CoREAS simulations. We discuss the updated calibration and its impact on the LOPES analysis results.

preprint2015arXiv

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

preprint2015arXiv

The lunar Askaryan technique: a technical roadmap

The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely large detector of ultra-high-energy particles. The high time resolution required to detect these pulses, which have a duration of around a nanosecond, puts this technique in a regime quite different from other forms of radio astronomy, with a unique set of associated technical challenges which have been addressed in a series of experiments by various groups. Implementing the methods and techniques developed by these groups for detecting lunar Askaryan pulses will be important for a future experiment with the Square Kilometre Array (SKA), which is expected to have sufficient sensitivity to allow the first positive detection using this technique. Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and the exclusion of spurious events from radio-frequency interference. We review the progress in each of these areas, and consider the further progress expected for future application with the SKA.

preprint2015arXiv

The Tunka Radio Extension: reconstruction of energy and shower maximum of the first year data (ICRC 2015)

Since its commissioning in autumn 2012, Tunka-Rex, the radio extension of the air-Cherenkov detector Tunka-133, performed three years of air shower measurements. Currently the detector consists of 44 antennas connected to air-Cherenkov and scintillator detectors, respectively, placed in the Tunka valley, Siberia. Triggered by these detectors, Tunka-Rex measures the radio signal up to EeV-scale air-showers. This configuration provides a unique possibility for cross-calibration between air-Cherenkov, radio and particle techniques. We present reconstruction methods for the energy and the shower maximum developed with CoREAS simulations, which allow for a precision competitive with the air-Cherenkov technique. We apply these methods to data acquired by Tunka-Rex in the first year which we use for cross-calibration, and we compare the results with the reconstruction of the energy and the shower maximum by Tunka-133, which provides also a reconstruction for the shower core used for the radio reconstruction. Our methods have shown that the atmospheric depth of the shower maximum ($X_\mathrm{max}$) can be reconstructed with a precision of better than 40 g/cm$^2$ for high quality events, in some cases even when only three antenna stations have signal. The energy precision is comparable with the air-Cherenkov precision of 15%. Soon the results will be checked with the independent data of the second year.

preprint2015arXiv

Tunka-Rex: the Cost-Effective Radio Extension of the Tunka Air-Shower Observatory

Tunka-Rex is the radio extension of the Tunka cosmic-ray observatory in Siberia close to Lake Baikal. Since October 2012 Tunka-Rex measures the radio signal of air-showers in coincidence with the non-imaging air-Cherenkov array Tunka-133. Furthermore, this year additional antennas will go into operation triggered by the new scintillator array Tunka-Grande measuring the secondary electrons and muons of air showers. Tunka-Rex is a demonstrator for how economic an antenna array can be without losing significant performance: we have decided for simple and robust SALLA antennas, and we share the existing DAQ running in slave mode with the PMT detectors and the scintillators, respectively. This means that Tunka-Rex is triggered externally, and does not need its own infrastructure and DAQ for hybrid measurements. By this, the performance and the added value of the supplementary radio measurements can be studied, in particular, the precision for the reconstructed energy and the shower maximum in the energy range of approximately $10^{17}-10^{18}\,$eV. Here we show first results on the energy reconstruction indicating that radio measurements can compete with air-Cherenkov measurements in precision. Moreover, we discuss future plans for Tunka-Rex.

preprint2014arXiv

A method for high precision reconstruction of air shower Xmax using two-dimensional radio intensity profiles

The mass composition of cosmic rays contains important clues about their origin. Accurate measurements are needed to resolve long-standing issues such as the transition from Galactic to extragalactic origin, and the nature of the cutoff observed at the highest energies. Composition can be studied by measuring the atmospheric depth of the shower maximum Xmax of air showers generated by high-energy cosmic rays hitting the Earth's atmosphere. We present a new method to reconstruct Xmax based on radio measurements. The radio emission mechanism of air showers is a complex process that creates an asymmetric intensity pattern on the ground. The shape of this pattern strongly depends on the longitudinal development of the shower. We reconstruct Xmax by fitting two-dimensional intensity profiles, simulated with CoREAS, to data from the LOFAR radio telescope. In the dense LOFAR core, air showers are detected by hundreds of antennas simultaneously. The simulations fit the data very well, indicating that the radiation mechanism is now well-understood. The typical uncertainty on the reconstruction of Xmax for LOFAR showers is 17 g/cm^2.

preprint2014arXiv

A search for point sources of EeV photons

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85° to +20°, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

preprint2014arXiv

A Targeted Search for Point Sources of EeV Neutrons

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

preprint2014arXiv

First Experimental Characterization of Microwave Emission from Cosmic Ray Air Showers

We report the first direct measurement of the overall characteristics of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers in the GHz frequency range. Microwave signals have been detected for more than 30 showers with energies above 3*10^16 eV. The observations presented in this Letter are consistent with a mainly forward-directed and polarised emission process in the GHz frequency range. The measurements show that microwave radiation offers a new means of studying air showers at energies above 10^17 eV.

preprint2014arXiv

Lunar detection of ultra-high-energy cosmic rays and neutrinos with the Square Kilometre Array

The origin of the most energetic particles in nature, the ultra-high-energy (UHE) cosmic rays, is still a mystery. Only the most energetic of these have sufficiently small angular deflections to be used for directional studies, and their flux is so low that even the 3,000 km^2 Pierre Auger detector registers only about 30 cosmic rays per year of these energies. A method to provide an even larger aperture is to use the lunar Askaryan technique, in which ground-based radio telescopes search for the nanosecond radio flashes produced when a cosmic ray interacts with the Moon's surface. The technique is also sensitive to UHE neutrinos, which may be produced in the decays of topological defects from the early universe. Observations with existing radio telescopes have shown that this technique is technically feasible, and established the required procedure: the radio signal should be searched for pulses in real time, compensating for ionospheric dispersion and filtering out local radio interference, and candidate events stored for later analysis. For the Square Kilometre Array (SKA), this requires the formation of multiple tied-array beams, with high time resolution, covering the Moon, with either SKA1-LOW or SKA1-MID. With its large collecting area and broad bandwidth, the SKA will be able to detect the known flux of UHE cosmic rays using the visible lunar surface - millions of square km - as the detector, providing sufficient detections of these extremely rare particles to address the mystery of their origin.

preprint2014arXiv

Reconstruction of inclined air showers detected with the Pierre Auger Observatory

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

preprint2014arXiv

Reconstruction of the energy and depth of maximum of cosmic-ray air-showers from LOPES radio measurements

LOPES is a digital radio interferometer located at Karlsruhe Institute of Technology (KIT), Germany, which measures radio emission from extensive air showers at MHz frequencies in coincidence with KASCADE-Grande. In this article, we explore a method (slope method) which leverages the slope of the measured radio lateral distribution to reconstruct crucial attributes of primary cosmic rays. First, we present an investigation of the method on the basis of pure simulations. Second, we directly apply the slope method to LOPES measurements. Applying the slope method to simulations, we obtain uncertainties on the reconstruction of energy and depth of shower maximum Xmax of 13% and 50 g/cm^2, respectively. Applying it to LOPES measurements, we are able to reconstruct energy and Xmax of individual events with upper limits on the precision of 20-25% for the primary energy and 95 g/cm^2 for Xmax, despite strong human-made noise at the LOPES site.

preprint2014arXiv

Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of $10^{19}$ eV at the Pierre Auger Observatory and the Telescope Array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

preprint2014arXiv

The wavefront of the radio signal emitted by cosmic ray air showers

Analyzing measurements of the LOPES antenna array together with corresponding CoREAS simulations for more than 300 measured events with energy above $10^{17}\,$eV and zenith angles smaller than $45^\circ$, we find that the radio wavefront of cosmic-ray air showers is of approximately hyperbolic shape. The simulations predict a slightly steeper wavefront towards East than towards West, but this asymmetry is negligible against the measurement uncertainties of LOPES. At axis distances $\gtrsim 50\,$m, the wavefront can be approximated by a simple cone. According to the simulations, the cone angle is clearly correlated with the shower maximum. Thus, we confirm earlier predictions that arrival time measurements can be used to study the longitudinal shower development, but now using a realistic wavefront. Moreover, we show that the hyperbolic wavefront is compatible with our measurement, and we present several experimental indications that the cone angle is indeed sensitive to the shower development. Consequently, the wavefront can be used to statistically study the primary composition of ultra-high energy cosmic rays. At LOPES, the experimentally achieved precision for the shower maximum is limited by measurement uncertainties to approximately $140\,$g/cm$^2$. But the simulations indicate that under better conditions this method might yield an accuracy for the atmospheric depth of the shower maximum, $X_\mathrm{max}$, better than $30\,$g/cm$^2$. This would be competitive with the established air-fluorescence and air-Cherenkov techniques, where the radio technique offers the advantage of a significantly higher duty-cycle. Finally, the hyperbolic wavefront can be used to reconstruct the shower geometry more accurately, which potentially allows a better reconstruction of all other shower parameters, too.

preprint2013arXiv

Ankle-like Feature in the Energy Spectrum of Light Elements of Cosmic Rays Observed with KASCADE-Grande

Recent results of the KASCADE-Grande experiment provided evidence for a mild knee-like structure in the all-particle spectrum of cosmic rays at $E = 10^{16.92 \pm 0.10} \, \mathrm{eV}$, which was found to be due to a steepening in the flux of heavy primary particles. The spectrum of the combined components of light and intermediate masses was found to be compatible with a single power law in the energy range from $10^{16.3} \, \mathrm{eV}$ to $10^{18} \, \mathrm{eV}$. In this paper, we present an update of this analysis by using data with increased statistics, originating both from a larger data set including more recent measurements and by using a larger fiducial area. In addition, optimized selection criteria for enhancing light primaries are applied. We find a spectral feature for light elements, namely a hardening at $E = 10^{17.08 \pm 0.08} \, \mathrm{eV}$ with a change of the power law index from $-3.25 \pm 0.05$ to $-2.79 \pm 0.08$.

preprint2013arXiv

Comparison of LOPES data and CoREAS simulations using a full detector simulation (ICRC2013)

The LOPES experiment at the Karlsruhe Institute of Technology, Germany, has been measuring radio emission of air showers for almost 10 years. For a better understanding of the emission process a detailed comparison of data with simulations is necessary. This is possible using a newly developed detector simulation including all LOPES detector components. After propagating a simulated event through this full detector simulation a standard LOPES like event file is written. LOPES data and CoREAS simulations can then be treated equally and the same analysis software can be applied to both. This gives the opportunity to compare data and simulations directly. Furthermore, the standard analysis software can be used with simulations which provide the possibility to check the accuracy regarding reconstruction of air shower parameters. We point out the advantages and present first results using such a full LOPES detector simulation. A comparison of LOPES data and the Monte Carlo code CoREAS based on an analysis using this detector simulation is shown.

preprint2013arXiv

Cosmic Ray Measurements with LOPES: Status and Recent Results (ARENA 2012)

LOPES is a digital antenna array at the Karlsruhe Institute of Technology, Germany, for cosmic-ray air-shower measurements. Triggered by the co-located KASCADE-Grande air-shower array, LOPES detects the radio emission of air showers via digital radio interferometry. We summarize the status of LOPES and recent results. In particular, we present an update on the reconstruction of the primary-particle properties based on almost 500 events above 100 PeV. With LOPES, the arrival direction can be reconstructed with a precision of at least 0.65°, and the energy with a precision of at least 20 %, which, however, does not include systematic uncertainties on the absolute energy scale. For many particle and astrophysics questions the reconstruction of the atmospheric depth of the shower maximum, Xmax, is important, since it yields information on the type of the primary particle and its interaction with the atmosphere. Recently, we found experimental evidence that the slope of the radio lateral distribution is indeed sensitive to the longitudinal development of the air shower, but unfortunately, the Xmax precision at LOPES is limited by the high level of anthropogenic radio background. Nevertheless, the developed methods can be transferred to next generation experiments with lower background, which should provide an Xmax precision competitive to other detection technologies.

preprint2013arXiv

Energy spectra of KASCADE-Grande based on shower size measurements and different hadronic interaction models

KASCADE-Grande is dedicated for investigations of cosmic-ray air showers in the primary energy range from 10 PeV to 1 EeV. The multi-detector system allows us to reconstruct charged particles, electron and muon numbers for individual air showers with high accuracies. Based on the shower size ($N_{ch}$) spectra of the charged particle component, the all-particle energy spectrum of cosmic rays is reconstructed, where attenuation effects in the atmosphere are corrected by applying the constant intensity cut method. The energy calibration is performed by using CORSIKA simulations with high-energy interaction models QGSJET-II-2, QGSJET-II-4, EPOS 1.99 and SIBYLL 2.1, where FLUKA has been used as low-energy interaction model for all cases. In the different hadronic models, different abundances for shower particles are predicted. Such model differences in the observables will be compared and discussed in this contribution. Furthermore, by using data with increasing statistics, the updated energy spectra by means of different interaction models will be presented.

preprint2013arXiv

High-energy cosmic rays measured with KASCADE-Grande

The detection of high-energy cosmic rays above a few hundred TeV is realized by the observation of extensive air-showers. By using the multi-detector setup of KASCADE-Grande, energy spectrum, elemental composition, and anisotropies of high-energy cosmic rays in the energy range from below the knee up to 2 EeV are investigated. In addition, the large high-quality data set permits distinct tests of the validity of hadronic interaction models used in interpreting air-shower measurements. After more than 16 years, the KASCADE-Grande experiment terminated measurements end of 2012. This contribution will give an overview of the main results of the data analysis achieved so far, and will report about the status of KCDC, the KASCADE Cosmic-ray Data Center, where via a web-based interface the data will be made available for the interested public.

preprint2013arXiv

Investigation on the energy and mass composition of cosmic rays using LOPES radio data

The sensitivity to the mass composition as well as the reconstruction of the energy of the primary particle are explored here by leveraging the features of the radio lateral distribution function. For the purpose of this analysis, a set of events measured with the LOPES experiment is reproduced with the latest CoREAS radio simulation code. Based on simulation predictions, a method which exploits the slope of the radio lateral distribution function is developed (Slope Method) and directly applied on measurements. As a result, the possibility to reconstruct both the energy and the depth of the shower maximum of the cosmic ray air shower using radio data and achieving relatively small uncertainties is presented.

preprint2013arXiv

KASCADE-Grande measurements of energy spectra for elemental groups of cosmic rays

The KASCADE-Grande experiment, located at KIT-Karlsruhe, Germany, consists of a large scintillator array for measurements of charged particles, N_ch, and of an array of shielded scintillation counters used for muon counting, N_mu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to 1000 PeV, thereby enabling the verification of a knee in the iron spectrum expected at approximately 100 PeV. Exploring the composition in this energy range is of fundamental importance for understanding the transition from galactic to extragalactic cosmic rays. Following earlier studies of elemental spectra reconstructed in the knee energy range from KASCADE data, we have now extended these measurements to beyond 100 PeV. By analysing the two-dimensional shower size spectrum N_ch vs. N_mu, we reconstruct the energy spectra of different mass groups by means of unfolding methods. The procedure and its results, giving evidence for a knee-like structure in the spectrum of iron nuclei, will be presented.

preprint2013arXiv

KASCADE-Grande measurements of energy spectra for elemental groups of cosmic rays

The KASCADE-Grande air shower experiment [W. Apel, et al. (KASCADE-Grande collaboration), Nucl. Instrum. Methods A 620 (2010) 202] consists of, among others, a large scintillator array for measurements of charged particles, Nch, and of an array of shielded scintillation counters used for muon counting, Nmu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to about 2000 PeV, where exploring the composition is of fundamental importance for understanding the transition from galactic to extragalactic origin of cosmic rays. Following earlier studies of the all-particle and the elemental spectra reconstructed in the knee energy range from KASCADE data [T. Antoni, et al. (KASCADE collaboration), Astropart. Phys. 24 (2005) 1], we have now extended these measurements to beyond 200 PeV. By analysing the two-dimensional shower size spectrum Nch vs. Nmu for nearly vertical events, we reconstruct the energy spectra of different mass groups by means of unfolding methods over an energy range where the detector is fully efficient. The procedure and its results, which are derived based on the hadronic interaction model QGSJET-II-02 and which yield a strong indication for a dominance of heavy mass groups in the covered energy range and for a knee-like structure in the iron spectrum at around 80 PeV, are presented. This confirms and further refines the results obtained by other analyses of KASCADE-Grande data, which already gave evidence for a knee-like structure in the heavy component of cosmic rays at about 80 PeV [W. Apel, et al. (KASCADE-Grande collaboration), Phys. Rev. Lett. 107 (2011) 171104].

preprint2013arXiv

LOPES 3D - vectorial measurements of radio emission from cosmic ray induced air showers

LOPES 3D is able to measure all three components of the electric field vector of the radio emission from air showers. This allows a better comparison with emission models. The measurement of the vertical component increases the sensitivity to inclined showers. By measuring all three components of the electric field vector LOPES 3D demonstrates by how much the reconstruction accuracy of primary cosmic ray parameters increases. Thus LOPES 3D evaluates the usefulness of vectorial measurements for large scale applications.

preprint2013arXiv

LOPES 3D reconfiguration and first measurements

The Radio detection technique of high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower due to their deflection in the Earth's magnetic field. The LOPES experiment at Karlsruhe Institute of Technology, Germany with its simple dipoles made major contributions to the revival of this technique. LOPES is working in the frequency range from 40 to 80 MHz and was reconfigured several times to improve and further develop the radio detection technique. In the current setup LOPES consists of 10 tripole antennas which measure the complete electric field vector of the radio emission from cosmic rays. LOPES is the first experiment measuring all three vectorial components at once and thereby gaining the full information about the electric field vector and not only a two-dimensional projection. Such a setup including also measurements of the vertical electric field component is expected to increase the sensitivity to inclined showers and help to advance the understanding of the emission mechanism. We present the reconfiguration and calibration procedure of LOPES 3D and discuss first measurements.

preprint2013arXiv

LOPES-3D, an antenna array for full signal detection of air-shower radio emission

To better understand the radio signal emitted by extensive air-showers and to further develop the radio detection technique of high-energy cosmic rays, the LOPES experiment was reconfigured to LOPES-3D. LOPES-3D is able to measure all three vectorial components of the electric field of radio emission from cosmic ray air showers. The additional measurement of the vertical component ought to increase the reconstruction accuracy of primary cosmic ray parameters like direction and energy, provides an improved sensitivity to inclined showers, and will help to validate simulation of the emission mechanisms in the atmosphere. LOPES-3D will evaluate the feasibility of vectorial measurements for large scale applications. In order to measure all three electric field components directly, a tailor-made antenna type (tripoles) was deployed. The change of the antenna type necessitated new pre-amplifiers and an overall recalibration. The reconfiguration and the recalibration procedure are presented and the operationality of LOPES-3D is demonstrated.

preprint2013arXiv

Mass sensitivity in the radio lateral distribution function

Measuring the mass composition of ultra-high energy cosmic rays is one of the main tasks in the cosmic rays field. Here we are exploring the composition signature in the coherent electromagnetic emission from extensive air showers, detected in the MHz frequency range. One of the experiments that successfully detects radio events in the frequency band of 40-80 MHz is the LOPES experiment at KIT. It is a digital interferometric antenna array and has the important advantage of taking data in coincidence with the particle detector array KASCADE-Grande. A possible method to look at the composition signature in the radio data, predicted by simulations, concerns the radio lateral distribution function, since its slope is strongly correlated with Xmax. Recent comparison between REAS3 simulations and LOPES data showed a significantly improved agreement in the lateral distribution function and for this reason an analysis on a possible LOPES mass signature through the slope method is promising. Trying to reproduce a realistic case, proton and iron showers are simulated with REAS3 using the LOPES selection information as input parameters. The obtained radio lateral distribution slope is analyzed in detail. The lateral slope method to look at the composition signature in the radio data is shown here and a possible signature of mass composition in the LOPES data is discussed.

preprint2013arXiv

Observation of Polarised Microwave Emission from Cosmic Ray Air Showers

We report on the first direct measurement of the basic features of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers. Microwave signals have been detected for more than 30 showers with energies above $3\times10^{16}$\,eV. The observations presented in this Letter are consistent with a mainly forward-beamed, coherent and polarised emission process in the GHz frequency range. An isotropic, unpolarised radiation is disfavoured as the dominant emission model. The measurements show that microwave radiation offers a new means of studying air showers at very high energy.

preprint2013arXiv

Reconstructing energy and Xmax of cosmic ray air showers using the radio lateral distribution measured with LOPES

The LOPES experiment, a digital radio interferometer located at KIT (Karlsruhe Institute of Technology), obtained remarkable results for the detection of radio emission from extensive air showers at MHz frequencies. Features of the radio lateral distribution function (LDF) measured by LOPES are explored in this work for a precise reconstruction of two fundamental air shower parameters: the primary energy and the shower Xmax. The method presented here has been developed on (REAS3-)simulations, and is applied to LOPES measurements. Despite the high human-made noise at the LOPES site, it is possible to reconstruct both the energy and Xmax for individual events. On the one hand, the energy resolution is promising and comparable to the one of the co-located KASCADE-Grande experiment. On the other hand, Xmax values are reconstructed with the LOPES measurements with a resolution of 90 g/cm2 . A precision on Xmax better than 30 g/cm2 is predicted and achievable in a region with a lower human-made noise level.

preprint2013arXiv

Simulating radio emission from air showers with CoREAS

CoREAS is a Monte Carlo code for the simulation of radio emission from extensive air showers. It implements the endpoint formalism for the calculation of electromagnetic radiation directly in CORSIKA. As such, it is parameter-free, makes no assumptions on the emission mechanism for the radio signals, and takes into account the complete complexity of the electron and positron distributions as simulated by CORSIKA. In this article, we illustrate the capabilities of CoREAS with simulations carried out in different frequency ranges from tens of MHz up to GHz frequencies, and describe in particular the emission characteristics at high frequencies due to Cherenkov effects arising from the varying refractive index of the atmosphere.

preprint2013arXiv

Study of the muon content of very high-energy EAS measured with the KASCADE-Grande observatory

The KASCADE-Grande detector is an air-shower array devoted to the study of primary cosmic rays with very high-energies (E = 10^{16} - 10^{18} eV). The instrument is composed of different particle detector systems suitable for the detailed study of the properties of Extensive Air Showers (EAS) developed by cosmic rays in the atmosphere. Among the EAS observables studied with the detector, the charged number of particles, the muon content (at different energy thresholds), and the number of electrons are found. By comparing the measurements of these air-shower parameters with the expectations from MC simulations, different hadronic interaction models can be tested at the high-energy regime with the KASCADE-Grande experiment. In this work, the results of a study on the evolution of the muon content of EAS with zenith angle, performed with the KASCADE-Grande instrument, is presented. Measurements are compared with predictions from MC simulations based on the QGSJET II, QGSJET II-04, SIBYLL 2.1 and EPOS 1.99 hadronic interaction models. A mismatch between experiment and simulations is observed. A similar problem is found for the evolution of the lateral distribution function of muons in the atmosphere.

preprint2013arXiv

The <lnA> study in the primary energy range 10^{16} - 10^{17} eV with the Muon Tracking Detector in the KASCADE-Grande experiment

The KASCADE-Grande Muon Tracking Detector enables with high accuracy the measurement of directions of EAS muons with energy above 0.8 GeV and up to 700 m distance from the shower centre. Reconstructed muon tracks are used to investigate muon pseudorapidity (eta) distributions. These distributions are nearly identical to the pseudorapidity distributions of their parent mesons produced in hadronic interactions. Comparison of the eta distributions from measured and simulated showers can be used to test the quality of the high energy hadronic interaction models. In this context a comparison of the QGSJet-II-2 and QGSJet-II-4 model will be shown. The pseudorapidity distributions reflect the longitudinal development of EAS and, as such, are sensitive to the mass of the cosmic rays primary particles. With various parameters of the eta distribution, obtained from the MTD data, it is possible to calculate the mean logarithmic mass of CRs. The results of the <lnA> analysis in the primary energy range 10^{16} eV - 10^{17} eV with the 1st quartile (Q1) of eta distribution will be presented.

preprint2013arXiv

The energy-spectrum of light primaries in the range from 10^{16.6} to 10^{18.2} eV

Data of the Grande extension of the KASCADE experiment allows us to study extensive air showers induced by primary cosmic rays with energies above 10^{16} eV. The energy of an event is estimated in terms of the number of charged particles (Nch ) and the number of muons (Nμ ) measured at an altitude of 110 m a.s.l. While a combination of the two numbers is used for the energy, the ratio defines the primary mass (group). The spectrum of the combined light and medium mass components, recently measured with KASCADE-Grande, was found to be compatible with both a single power-law and a broken power-law in the energy range between 10^{16.3} and 10^{18} eV. In this contribution we will present the investigation of possible structures in the spectrum of light primaries with increased statistics both from a larger data set including more recent measurements and by using a larger fiducial area than in the previous study. With the better statistical accuracy and with optimized selection criteria for enhancing light primaries we have found evidence for a hardening (ankle) of the spectrum of the light component at an energy of 10^{17.08+-0.08} eV.

preprint2013arXiv

The Tunka Radio Extension (Tunka-Rex): Status and First Results (ICRC 2013)

Tunka-Rex is a new radio antenna array which extends the Tunka experiment in Siberia close to lake Baikal. It consists of 20 antennas on an area of 1 km^2 which measure the radio emission of high-energy air showers. Tunka-Rex is triggered by the photomultiplier array of Tunka measuring air-Cherenkov light of air showers in the energy range from about 10 PeV to 1 EeV. This configuration allows for the worldwide first hybrid measurements of the radio and air-Cherenkov signal for the same events: an ideal situation to perform a cross-calibration between both methods. Consequently, the main goal of Tunka-Rex is to determine the achievable energy and Xmax precision of radio measurements by comparing them to the reconstruction of the air-Cherenkov measurements. Tunka-Rex started operation in autumn 2012, and already detected air-shower events. In this paper we present the status of Tunka-Rex and first results which indicate that Tunka-Rex measures indeed the radio emission by air showers and that is is sensitive to their energy.

preprint2013arXiv

The Tunka-Rex antenna station (ICRC 2013)

Tunka-Rex is the radio extension of Tunka-133, a 1 km^2 air-Cherenkov Detector for air showers in Siberia. Tunka-Rex began operation on October 8th 2012 with 20 radio antennas. Its main goals are to explore the possible precision of the radio detection technique in determination of primary energy and mass. Each radio antenna station consists of two perpendicular aligned active SALLA antennas, which receive the radio signal from air showers. The preamplified radio signal is transmitted to local cluster centers of the Tunka-133 DAQ, where it is filtered, amplified and digitized. To reconstruct the radio signal it is crucial to understand how it is affected in each of these steps. Thus, we have studied the combined response of the antenna, with its directional pattern and the analog electronics chain, consisting of a Low-Noise Amplifier and a filter amplifier. We discuss the hardware setup of Tunka-Rex and how a description of its response is obtained. Furthermore, we estimate systematic uncertainties on the reconstructed radio signal due to hardware effects (e.g., slight variations of the electronics properties). Finally, we present background measurements with the actual Tunka-Rex antennas.

preprint2013arXiv

Theory and simulations of air shower radio emission

A precise understanding of the radio emission from extensive air showers is of fundamental importance for the design of cosmic ray radio detectors as well as the analysis and interpretation of their data. In recent years, tremendous progress has been made in the understanding of the emission physics both in macroscopic and microscopic frameworks. A consistent picture has emerged: the emission stems mainly from time-varying transverse currents and a time-varying charge excess; in addition, Cherenkov-like compression of the emission due to the refractive index gradient in the atmosphere can lead to time-compression of the emitted pulses and thus high-frequency contributions in the signal. In this article, I discuss the evolution of the modelling in recent years, present the emission physics as it is understood today, and conclude with a description and comparison of the models currently being actively developed.

preprint2013arXiv

Thunderstorm Observations by Air-Shower Radio Antenna Arrays

Relativistic, charged particles present in extensive air showers lead to a coherent emission of radio pulses which are measured to identify the shower initiating high-energy cosmic rays. Especially during thunderstorms, there are additional strong electric fields in the atmosphere, which can lead to further multiplication and acceleration of the charged particles and thus have influence on the form and strength of the radio emission. For a reliable energy reconstruction of the primary cosmic ray by means of the measured radio signal it is very important to understand how electric fields affect the radio emission. In addition, lightning strikes are a prominent source of broadband radio emissions that are visible over very long distances. This, on the one hand, causes difficulties in the detection of the much lower signal of the air shower. On the other hand the recorded signals can be used to study features of the lightning development. The detection of cosmic rays via the radio emission and the influence of strong electric fields on this detection technique is investigated with the LOPES experiment in Karlsruhe, Germany. The important question if a lightning is initiated by the high electron density given at the maximum of a high-energy cosmic-ray air shower is also investigated, but could not be answered by LOPES. But, these investigations exhibit the capabilities of EAS radio antenna arrays for lightning studies. We report about the studies of LOPES measured radio signals of air showers taken during thunderstorms and give a short outlook to new measurements dedicated to search for correlations of lightning and cosmic rays.

preprint2013arXiv

Tunka-Rex: Status and Results of the First Measurements (RICAP 2013)

Tunka-Rex is the radio extension of Tunka-133 located in Siberia close to Lake Baikal. The latter is a photomultiplier array registering air-Cherenkov light from air showers induced by cosmic-ray particles with initial energies of approximately $10^{16}$ to $10^{18}$ eV. Tunka-Rex extends this detector with 25 antennas spread over an area of 1 km$^2$. It is triggered externally by Tunka-133, and detects the radio emission of the same air showers. The combination of an air-Cherenkov and a radio detector provides a facility for hybrid measurements and cross-calibration between the two techniques. The main goal of Tunka-Rex is to determine the precision of the reconstruction of air-shower parameters using the radio detection technique. It started operation in autumn 2012. We present the overall concept of Tunka-Rex, the current status of the array and first analysis results.

preprint2013arXiv

Vectorial Radio Interferometry with LOPES 3D

One successful detection technique for high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower. The LOPES experiment at Karlsruhe Institute of Technology, Germany, has made major contributions to the evolution of this technique. LOPES was reconfigured several times to improve and further develop the radio detection technique. In the latest setup LOPES consisted of 10 tripole antennas. With this, LOPES 3D was the first cosmic ray experiment measuring all three vectorial field components at once and thereby gaining the full information about the electric field vector. We present an analysis based on the data taken with special focus on the benefits of a direct measurement of the vertical polarization component. We demonstrate that by measuring all polarization components the detection and reconstruction efficiency is increased and noisy single channel data can be reconstructed by utilising the information from the other two channels of one antenna station.

preprint2012arXiv

A search for ultra-high energy neutrinos in highly inclined events at the Pierre Auger Observatory

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

preprint2012arXiv

Antennas for the Detection of Radio Emission Pulses from Cosmic-Ray induced Air Showers at the Pierre Auger Observatory

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

preprint2012arXiv

Description of Atmospheric Conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS)

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malargüe and averaged monthly models, the utility of the GDAS data is shown.

preprint2012arXiv

Experimental evidence for the sensitivity of the air-shower radio signal to the longitudinal shower development

We observe a correlation between the slope of radio lateral distributions, and the mean muon pseudorapidity of 59 individual cosmic-ray-air-shower events. The radio lateral distributions are measured with LOPES, a digital radio interferometer co-located with the multi-detector-air-shower array KASCADE-Grande, which includes a muon-tracking detector. The result proves experimentally that radio measurements are sensitive to the longitudinal development of cosmic-ray air-showers. This is one of the main prerequisites for using radio arrays for ultra-high-energy particle physics and astrophysics.

preprint2012arXiv

Radar reflection off extensive air showers

We investigate the possibility of detecting extensive air showers by the radar technique. Considering a bistatic radar system and different shower geometries, we simulate reflection of radio waves off the static plasma produced by the shower in the air. Using the Thomson cross-section for radio wave reflection, we obtain the time evolution of the signal received by the antennas. The frequency upshift of the radar echo and the power received are studied to verify the feasibility of the radar detection technique.

preprint2012arXiv

The Rapid Atmospheric Monitoring System of the Pierre Auger Observatory

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

preprint2012arXiv

The spectrum of high-energy cosmic rays measured with KASCADE-Grande

The energy spectrum of cosmic rays between 10**16 eV and 10**18 eV, derived from measurements of the shower size (total number of charged particles) and the total muon number of extensive air showers by the KASCADE-Grande experiment, is described. The resulting all-particle energy spectrum exhibits strong hints for a hardening of the spectrum at approximately 2x10**16 eV and a significant steepening at c. 8x10**16 eV. These observations challenge the view that the spectrum is a single power law between knee and ankle. Possible scenarios generating such features are discussed in terms of astrophysical processes that may explain the transition region from galactic to extragalactic origin of cosmic rays.

preprint2011arXiv

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

preprint2011arXiv

Anisotropy and chemical composition of ultra-high energy cosmic rays using arrival directions measured by the Pierre Auger Observatory

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

preprint2011arXiv

First results of the CROME experiment

It is expected that a radio signal in the microwave range is produced in the atmosphere due to molecular bremsstrahlung initiated by extensive air showers. The CROME (Cosmic-Ray Observation via Microwave Emission) experiment was built to search for this microwave signal. Radiation from the atmosphere is monitored in the extended C band (3.4--4.2 GHz) in coincidence with showers detected by the KASCADE-Grande experiment. The detector setup consists of several parabolic antennas and fast read-out electronics. The sensitivity of the detector has been measured with different methods. First results after half a year of data taking are presented.

preprint2011arXiv

Kneelike structure in the spectrum of the heavy component of cosmic rays observed with KASCADE-Grande

We report the observation of a steepening in the cosmic ray energy spectrum of heavy primary particles at about $8 \times 10^{16}$ eV. This structure is also seen in the all-particle energy spectrum, but is less significant. Whereas the `knee' of the cosmic ray spectrum at $3-5 \times 10^{15}$ eV was assigned to light primary masses by the KASCADE experiment, the new structure found by the KASCADE-Grande experiment is caused by heavy primaries. The result is obtained by independent measurements of the charged particle and muon components of the secondary particles of extensive air showers in the primary energy range of $10^{16}$ to $10^{18}$ eV. The data are analyzed on a single-event basis taking into account also the correlation of the two observables.

preprint2011arXiv

Search for First Harmonic Modulation in the Right Ascension Distribution of Cosmic Rays Detected at the Pierre Auger Observatory

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.

preprint2011arXiv

Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

preprint2011arXiv

The effect of the geomagnetic field on cosmic ray energy estimates and large scale anisotropy searches on data from the Pierre Auger Observatory

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

preprint2011arXiv

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^{17} and 10^{19} eV and zenith angles up to 65 degs. 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.

preprint2011arXiv

Trigger and Aperture of the Surface Detector Array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

preprint2010arXiv

Longitudinal EAS-Development Studies in the Air-Shower Experiment KASCADE-Grande

A large area (128 m^2) Muon Tracking Detector (MTD), located within the KASCADE experiment, has been built with the aim to identify muons (E_mu > 0.8 GeV) and their directions in extensive air showers by track measurements under more than 18 r.l. shielding. The orientation of the muon track with respect to the shower axis is expressed in terms of the radial- and tangential angles. By means of triangulation the muon production height H_mu is determined. By means of H_mu, a transition from light to heavy cosmic ray primary particle with increasing shower energy Eo from 1-10 PeV is observed. Muon pseudorapidity distributions for the first interactions above 15 km are studied and compared to Monte Carlo simulations.

preprint2010arXiv

New method for the time calibration of an interferometric radio antenna array

Digital radio antenna arrays, like LOPES (LOFAR PrototypE Station), detect high-energy cosmic rays via the radio emission from atmospheric extensive air showers. LOPES is an array of dipole antennas placed within and triggered by the KASCADE-Grande experiment on site of the Karlsruhe Institute of Technology, Germany. The antennas are digitally combined to build a radio interferometer by forming a beam into the air shower arrival direction which allows measurements even at low signal-to-noise ratios in individual antennas. This technique requires a precise time calibration. A combination of several calibration steps is used to achieve the necessary timing accuracy of about 1 ns. The group delays of the setup are measured, the frequency dependence of these delays (dispersion) is corrected in the subsequent data analysis, and variations of the delays with time are monitored. We use a transmitting reference antenna, a beacon, which continuously emits sine waves at known frequencies. Variations of the relative delays between the antennas can be detected and corrected for at each recorded event by measuring the phases at the beacon frequencies.

preprint2010arXiv

On noise treatment in radio measurements of cosmic ray air showers

Precise measurements of the radio emission by cosmic ray air showers require an adequate treatment of noise. Unlike to usual experiments in particle physics, where noise always adds to the signal, radio noise can in principle decrease or increase the signal if it interferes by chance destructively or constructively. Consequently, noise cannot simply be subtracted from the signal, and its influence on amplitude and time measurement of radio pulses must be studied with care. First, noise has to be determined consistently with the definition of the radio signal which typically is the maximum field strength of the radio pulse. Second, the average impact of noise on radio pulse measurements at individual antennas is studied for LOPES. It is shown that a correct treatment of noise is especially important at low signal-to-noise ratios: noise can be the dominant source of uncertainty for pulse height and time measurements, and it can systematically flatten the slope of lateral distributions. The presented method can also be transfered to other experiments in radio and acoustic detection of cosmic rays and neutrinos.

preprint2010arXiv

Prospects for determining air shower characteristics through geosynchrotron emission arrival times

Using simulations of geosynchrotron radiation from extensive air showers, we present a relation between the shape of the geosynchrotron radiation front and the distance of the observer to the maximum of the air shower. By analyzing the relative arrival times of radio pulses at several radio antennas in an air shower array, this relation may be employed to estimate the depth of maximum of an extensive air shower if its impact position is known, allowing an estimate for the primary particle's species. Vice versa, the relation provides an estimate for the impact position of the shower's core if an external estimate of the depth of maximum is available. In realistic circumstances, the method delivers reconstruction uncertainties down to 30 g/cm^2 when the distance to the shower core does not exceed 7 km. The method requires that the arrival direction is known with high precision.

preprint2010arXiv

Sensitivity of KASCADE-Grande data to hadronic interaction models

KASCADE-Grande is a large detector array dedicated for studies of high-energy cosmic rays in the primary energy range from 100 TeV to 1 EeV. The multi-detector concept of the experimental set-up offers the possibility to measure simultaneously various observables related to the electromagnetic, muonic, and hadronic air shower components. The experimental data are compared to predictions of CORSIKA simulations using high-energy hadronic interaction models (e.g. QGSJET or EPOS), as well as low-energy interaction models (e.g. FLUKA or GHEISHA). This contribution will summarize the results of such investigations. In particular, the validity of the new EPOS version 1.99 for EAS with energy around 100 PeV will be discussed.

preprint2010arXiv

Simulation of radio emission from air showers in atmospheric electric fields

We study the effect of atmospheric electric fields on the radio pulse emitted by cosmic ray air showers. Under fair weather conditions the dominant part of the radio emission is driven by the geomagnetic field. When the shower charges are accelerated and deflected in an electric field additional radiation is emitted. We simulate this effect with the Monte Carlo code REAS2, using CORSIKA-simulated showers as input. In both codes a routine has been implemented that treats the effect of the electric field on the shower particles. We find that the radio pulse is significantly altered in background fields of the order of ~100 V/cm and higher. Practically, this means that air showers passing through thunderstorms emit radio pulses that are not a reliable measure for the shower energy. Under other weather circumstances significant electric field effects are expected to occur rarely, but nimbostratus clouds can harbor fields that are large enough. In general, the contribution of the electric field to the radio pulse has polarization properties that are different from the geomagnetic pulse. In order to filter out radio pulses that have been affected by electric field effects, radio air shower experiments should keep weather information and perform full polarization measurements of the radio signal.

preprint2010arXiv

The KASCADE-Grande experiment: measurements of the all-particle energy spectrum of cosmic rays

The all-particle energy spectrum as measured by the KASCADE-Grande experiment for E = 10^{16} - 10^{18} eV is presented within the framework of the QGSJET II/FLUKA hadronic interaction models. Three different methods were applied based on the muon size and the total number of charged particles individually and in combination. From the study it is found that the spectrum cannot be completely described by a smooth power law due to the presence of characteristic features.

preprint2010arXiv

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the Véron-Cetty and Véron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

preprint2009arXiv

Lateral Distribution of the Radio Signal in Extensive Air Showers Measured with LOPES

The antenna array LOPES is set up at the location of the KASCADE-Grande extensive air shower experiment in Karlsruhe, Germany and aims to measure and investigate radio pulses from Extensive Air Showers. The coincident measurements allow us to reconstruct the electric field strength at observation level in dependence of general EAS parameters. In the present work, the lateral distribution of the radio signal in air showers is studied in detail. It is found that the lateral distributions of the electric field strengths in individual EAS can be described by an exponential function. For about 20% of the events a flattening towards the shower axis is observed, preferentially for showers with large inclination angle. The estimated scale parameters R0 describing the slope of the lateral profiles range between 100 and 200 m. No evidence for a direct correlation of R0 with shower parameters like azimuth angle, geomagnetic angle, or primary energy can be found. This indicates that the lateral profile is an intrinsic property of the radio emission during the shower development which makes the radio detection technique suitable for large scale applications.

preprint2009arXiv

Monte Carlo simulations of air showers in atmospheric electric fields

The development of cosmic ray air showers can be influenced by atmospheric electric fields. Under fair weather conditions these fields are small, but the strong fields inside thunderstorms can have a significant effect on the electromagnetic component of a shower. Understanding this effect is particularly important for radio detection of air showers, since the radio emission is produced by the shower electrons and positrons. We perform Monte Carlo simulations to calculate the effects of different electric field configurations on the shower development. We find that the electric field becomes important for values of the order of 1 kV/cm. Not only can the energy distribution of electrons and positrons change significantly for such field strengths, it is also possible that runaway electron breakdown occurs at high altitudes, which is an important effect in lightning initiation.

preprint2009arXiv

Simulations and theory of radio emission from cosmic ray air showers

In the last few years, interest in radio detection of cosmic ray air showers has risen continuously. By now, large-scale application of the radio technique is under investigation in the framework of LOFAR and the Pierre Auger Observatory. The experimental efforts are accompanied by new approaches to describe and model the underlying radiation mechanisms, to lay the foundation for an interpretation of the experimental data. In this article, I review the current radio emission theory and simulations, with slight focus on the geosynchrotron model and its predictions for the information content of the radio signals.

preprint2009arXiv

The Fluorescence Detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

preprint2009arXiv

Universality of electron-positron distributions in extensive air showers

Using a large set of simulated extensive air showers, we investigate universality features of electron and positron distributions in very-high-energy cosmic-ray air showers. Most particle distributions depend only on the depth of the shower maximum and the number of particles in the cascade at this depth. We provide multi-dimensional parameterizations for the electron-positron distributions in terms of particle energy, vertical and horizontal momentum angle, lateral distance, and time distribution of the shower front. These parameterizations can be used to obtain realistic electron-positron distributions in extensive air showers for data analysis and simulations of Cherenkov radiation, fluorescence signal, and radio emission.

preprint2009arXiv

Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

preprint2007arXiv

A radio air shower surface detector as an extension for IceCube and IceTop

The IceCube neutrino detector is built into the Antarctic ice sheet at the South Pole to measure high energy neutrinos. For this, 4800 photomultiplier tubes (PMTs) are being deployed at depths between 1450 and 2450 meters into the ice to measure neutrino induced charged particles like muons. IceTop is a surface air shower detector consisting of 160 Cherenkov ice tanks located on top of IceCube. To extend IceTop, a radio air shower detector could be built to significantly increase the sensitivity at higher shower energies and for inclined showers. As air showers induced by cosmic rays are a major part of the muonic background in IceCube, IceTop is not only an air shower detector, but also a veto to reduce the background in IceCube. Air showers are detectable by radio signals with a radio surface detector. The major emission process is the coherent synchrotron radiation emitted by e+ e- shower particles in the Earths magnetic field (geosynchrotron effect). Simulations of the expected radio signals of air showers are shown. The sensitivity and the energy threshold of different antenna field configurations are estimated.