Demonstrating Agreement between Radio and Fluorescence Measurements of the Depth of Maximum of Extensive Air Showers at the Pierre Auger Observatory

We show, for the first time, radio measurements of the depth of shower maximum (X_{max}) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence dataset, and between a subset of air showers observed simultaneously with both radio and fluorescence techniques, a measurement setup unique to the Pierre Auger Observatory. Furthermore, we show radio X_{max} resolution as a function of energy and demonstrate the ability to make competitive high-resolution X_{max} measurements with even a sparse radio array. With this, we show that the radio technique is capable of cosmic-ray mass composition studies, both at Auger and at other experiments.

Limits to Gauge Coupling in the Dark Sector Set by the Nonobservation of Instanton-Induced Decay of Super-Heavy Dark Matter in the Pierre Auger Observatory Data

Instantons, which are nonperturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decaying in the Galactic halo. These particles could have been produced during the post-inflationary epoch and match the relic abundance of dark matter inferred today. The nonobservation of the signatures searched for allows us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: α_{X}≲0.09, for 10^{9}≲M_{X}/GeV<10^{19}. Conversely, we obtain that, for instance, a reduced coupling constant α_{X}=0.09 excludes masses M_{X}≳3×10^{13}  GeV. In the context of dark matter production from gravitational interactions alone, we illustrate how these bounds are complementary to those obtained on the Hubble rate at the end of inflation from the nonobservation of tensor modes in the cosmological microwave background.

The Relationship of Lightning Radio Pulse Amplitudes and Source Altitudes as Observed by LOFAR

When a lightning flash is propagating in the atmosphere it is known that especially the negative leaders emit a large number of very high frequency (VHF) radio pulses. It is thought that this is due to streamer activity at the tip of the growing negative leader. In this work, we have investigated the dependence of the strength of this VHF emission on the altitude of such emission for two lightning flashes as observed by the Low Frequency ARray (LOFAR) radio telescope. We find for these two flashes that the extracted amplitude distributions are consistent with a power-law, and that the amplitude of the radio emissions decreases very strongly with source altitude, by more than a factor of 2 from 1 km altitude up to 5 km altitude. In addition, we do not find any dependence on the extracted power-law with altitude, and that the extracted power-law slope has an average around 3, for both flashes.

Measurement of the Fluctuations in the Number of Muons in Extensive Air Showers with the Pierre Auger Observatory

We present the first measurement of the fluctuations in the number of muons in extensive air showers produced by ultrahigh energy cosmic rays. We find that the measured fluctuations are in good agreement with predictions from air shower simulations. This observation provides new insights into the origin of the previously reported deficit of muons in air shower simulations and constrains models of hadronic interactions at ultrahigh energies. Our measurement is compatible with the muon deficit originating from small deviations in the predictions from hadronic interaction models of particle production that accumulate as the showers develop.

Features of the Energy Spectrum of Cosmic Rays above 2.5×10^{18} eV Using the Pierre Auger Observatory

We report a measurement of the energy spectrum of cosmic rays above 2.5×10^{18}  eV based on 215 030 events. New results are presented: at about 1.3×10^{19}  eV, the spectral index changes from 2.51±0.03(stat)±0.05(syst) to 3.05±0.05(stat)±0.10(syst), evolving to 5.1±0.3(stat)±0.1(syst) beyond 5×10^{19}  eV, while no significant dependence of spectral features on the declination is seen in the accessible range. These features of the spectrum can be reproduced in models with energy-dependent mass composition. The energy density in cosmic rays above 5×10^{18}  eV is [5.66±0.03(stat)±1.40(syst)]×10^{53}  erg Mpc^{-3}.

Determining Electric Fields in Thunderclouds With the Radiotelescope LOFAR

An analysis is presented of electric fields in thunderclouds using a recently proposed method based on measuring radio emission from extensive air shower events during thunderstorm conditions. This method can be regarded as a tomography of thunderclouds using cosmic rays as probes. The data cover the period from December 2011 till August 2014. We have developed an improved fitting procedure to be able to analyze the data. Our measurements show evidence for the main negative-charge layer near the -10° isotherm. This we have seen for a winter as well as for a summer cloud where multiple events pass through the same cloud and also the vertical component of the electric field could be reconstructed. On the day of measurement of some cosmic-ray events showing evidence for strong fields, no lightning activity was detected within 100 km distance. For the winter events, the top heights were between 5 and 6 km, while in the summer, typical top heights of 9 km were seen. Large horizontal components in excess of 70 kV/m of the electric fields are observed in the middle and top layers.

Radio Emission Reveals Inner Meter-Scale Structure of Negative Lightning Leader Steps

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.

Needle-like structures discovered on positively charged lightning branches

Lightning is a dangerous yet poorly understood natural phenomenon. Lightning forms a network of plasma channels propagating away from the initiation point with both positively and negatively charged ends-called positive and negative leaders¹. Negative leaders propagate in discrete steps, emitting copious radio pulses in the 30-300-megahertz frequency band^2-8 that can be remotely sensed and imaged with high spatial and temporal resolution^9-11. Positive leaders propagate more continuously and thus emit very little high-frequency radiation¹². Radio emission from positive leaders has nevertheless been mapped^13-15, and exhibits a pattern that is different from that of negative leaders^11-13,16,17. Furthermore, it has been inferred that positive leaders can become transiently disconnected from negative leaders^{9,12,16,18-20}, which may lead to current pulses that both reconnect positive leaders to negative leaders^{11,16,17,20-22} and cause multiple cloud-to-ground lightning events¹. The disconnection process is thought to be due to negative differential resistance¹⁸, but this does not explain why the disconnections form primarily on positive leaders²², or why the current in cloud-to-ground lightning never goes to zero²³. Indeed, it is still not understood how positive leaders emit radio-frequency radiation or why they behave differently from negative leaders. Here we report three-dimensional radio interferometric observations of lightning over the Netherlands with unprecedented spatiotemporal resolution. We find small plasma structures-which we call 'needles'-that are the dominant source of radio emission from the positive leaders. These structures appear to drain charge from the leader, and are probably the reason why positive leaders disconnect from negative ones, and why cloud-to-ground lightning connects to the ground multiple times.

Recent results from the KASCADE-Grande data analysis

KASCADE, together with its extension KASCADE-Grande measured individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The data collection was fully completed at the end of 2013 and the experiment was dismantled. However, the data analysis is still in progress. Recently, we published a new result on upper limits to the flux of ultra-high energy gamma rays, which set constraints on some fundamental astrophysical models. We also use the data to investigate the validity of the new hadronic interactions models like SIBYLL version 2.3c or EPOS-LHC. In addition, we updated and improved the webbased platform of the KASCADE Cosmic Ray Data Centre (KCDC), where now the data from KASCADE and KASCADE-Grande of more than 20 years measurements is available, including corresponding Monte-Carlo simulated events based on three different hadronic interaction models. In this contribution, recent results from KASCADE-Grande and the update of KCDC is briefly discussed.

Study of themuon content of high-energy air showers with KASCADE-Grande

In this work, we report measurements on the muon content (Eth > 230 MeV) of extensive air showers (EAS) induced by cosmic rays with primary energy from 10 PeV up to 1 EeV performed with the KASCADE-Grande experiment. The measurements are confronted with SIBYLL 2.3. The results are focused on the dependence of the total muon number and the lateral density distribution of muons in EAS on the zenith angle and the total number of charged particles in the shower. We also present updated results of a detailed study of the attenuation length of shower muons, which reveal a deviation between the measured data and the predictions of the post-LHC hadronic interaction models SIBYLL 2.3, QGSJET-II-04 and EPOS-LHC.

First analysis of inclined air showers detected by Tunka-Rex

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array for the detection of radio emission from cosmic-ray air showers in the frequency band of 30 to 80 MHz and for primary energies above 100 PeV. The standard analysis of Tunka-Rex includes events with zenith angle of up to 50?. This cut is determined by the efficiency of the external trigger. However, due to the air-shower footprint increasing with zenith angle and due to the more efficient generation of radio emission (the magnetic field in the Tunka valley is almost vertical), there are a number of ultra-high-energy inclined events detected by Tunka-Rex. In this work we present a first analysis of a subset of inclined events detected by Tunka-Rex. We estimate the energies of the selected events and test the efficiency of Tunka-Rex antennas for detection of inclined air showers.

Summary of the main results of the KASCADE and KASCADE-Grande experiments

The KASCADE and KASCADE-Grande experiments operated in KIT-Campus North, Karlsruhe (Germany) from 1993 to 2012. The two experiments studied primary cosmic rays in the energy range from 1014 eV to 1018 eV, investigating the change of slope of the spectrum detected at 2 - 4 × 1015 eV, the so called knee. We briefly review the performance of the experiments and then the main results obtained in the operation of both experiments: the test of hadronic interaction models, the all particle primary spectrum, the elemental composition of primary cosmic rays (with the first claim of a knee-like feature of the heavy primaries spectrum) and the search for large scale anisotropies.

Signal recognition and background suppression by matched filters and neural networks for Tunka-Rex

The Tunka Radio Extension (Tunka-Rex) is a digital antenna array, which measures radio emission of the cosmic-ray air-showers in the frequency band of 30-80 MHz. Tunka-Rex is co-located with the TAIGA experiment in Siberia and consists of 63 antennas, 57 of them are in a densely instrumented area of about 1 km2. In the present workwe discuss the improvements of the signal reconstruction applied for Tunka-Rex. At the first stage we implemented matched filtering using averaged signals as template. The simulation study has shown that matched filtering allows one to decrease the threshold of signal detection and increase its purity. However, the maximum performanceof matched filtering is achievable only in case of white noise, while in reality the noise is not fully random due to different reasons. To recognize hidden features of the noise and treat them, we decided to use convolutional neural network with autoencoder architecture. Taking the recorded trace as an input, the autoencoder returns denoised traces, i.e. removes all signal-unrelated amplitudes. We present the comparison between the standard method of signal reconstruction, matched filtering and the autoencoder, and discuss the prospects of application of neural networks for lowering the threshold of digital antenna arrays for cosmic-ray detection.

Present status and prospects of the Tunka Radio Extension

The Tunka Radio Extension (Tunka-Rex) is a digital radio array operating in the frequency band of 30-80 MHz and detecting radio emission from air-showers produced by cosmic rays with energies above 100 PeV. The experimentis installed at the site of the TAIGA (Tunka Advanced Instrument for cosmic rays and Gamma Astronomy) observatory and performs joint measurements with the co-located particle and air-Cherenkov detectors in passive mode receiving a trigger from the latter. Tunka-Rex collects data since 2012, and during the last five years went throughseveral upgrades. As a result the density of the antenna field was increased by three times since its commission. In this contribution we present the latest results of Tunka-Rex experiment, particularly an updated analysis and efficiency study, which have been applied to the measurement of the mean shower maximum as a function of energy for cosmic rays of energies up to EeV. The future plans are also discussed: investigations towards an energy spectrum of cosmic rays with Tunka-Rex and their mass composition using a combination of Tunka-Rex data with muon measurements by the particle detector Tunka-Grande.

TAIGA - a hybrid array for high energy gamma astronomy and cosmic ray physics

The physics motivations and advantages of the new TAIGA (Tunka Advanced Instrument for cosmic ray physics and Gamma Astronomy) detector are presented. TAIGA aims at gamma-ray astronomy at energies from a few TeV to several PeV, as well as cosmic ray physics from 100 TeV to several EeV. For the energy range 30 – 200 TeV the sensitivity of 10 km2 area TAIGA array for the detection of local sources is expected to be 5 × 10-14 erg cm-2 sec-1 for 300 h of observations. Reconstruction of the given EAS energy, incoming direction and its core position, based on the timing TAIGA-HiSCORE data, allows one to increase a distance between the IACTs up to 600-1000 m. The low investments together with the high sensitivity for energies ≥ 30-50 TeV make this pioneering technique very attractive for exploring the galactic PeVatrons and cosmic rays. At present the TAIGA first stage has been constructed in Tunka valley, 50 km West from the Lake Baikal. The first experimental results of the TAIGA first stage are presented.

Tests of the SIBYLL 2.3 high-energy hadronic interaction model using the KASCADE-Grande muon data

The KASCADE-Grande observatory was a ground-based air shower array dedicated to study the energy and composition of cosmic rays in the energy interval E = 1 PeV –1 EeV. The experiment consisted of different detector systems which allowed the simultaneous measurement of distinct components of air showers (EAS), such as the muon content. In this contribution, we study the total muon number and the lateral density distribution of muons in EAS detected by KASCADE-Grande as a function of the zenith angle and the total number of charged particles. The attenuation length of the muon content of EAS is also measured. The results are compared with the predictions of the SIBYLL 2.3 high-energy hadronic interaction model.

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.

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8±0.7(stat)±6.7(syst)  MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principles calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

Accelerator Measurements of Magnetically Induced Radio Emission from Particle Cascades with Applications to Cosmic-Ray Air Showers

For 50 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 ultrahigh-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.

Search for patterns by combining cosmic-ray energy and arrival directions at the Pierre Auger Observatory

Energy-dependent patterns in the arrival directions of cosmic rays are searched for using data of the Pierre Auger Observatory. We investigate local regions around the highest-energy cosmic rays with [Formula: see text] eV by analyzing cosmic rays with energies above [Formula: see text] eV arriving within an angular separation of approximately 15[Formula: see text]. We characterize the energy distributions inside these regions by two independent methods, one searching for angular dependence of energy-energy correlations and one searching for collimation of energy along the local system of principal axes of the energy distribution. No significant patterns are found with this analysis. The comparison of these measurements with astrophysical scenarios can therefore be used to obtain constraints on related model parameters such as strength of cosmic-ray deflection and density of point sources.

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 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 polarized emission process in the GHz frequency range. The measurements show that microwave radiation offers a new means of studying air showers at E≥10^{17}  eV.

Measurement of the proton-air cross section at √s=57 TeV with the Pierre Auger Observatory

We report a measurement of the proton-air cross section for particle production at the center-of-mass energy per nucleon of 57 TeV. This is derived from the distribution of the depths of shower maxima observed with the Pierre Auger Observatory: systematic uncertainties are studied in detail. Analyzing the tail of the distribution of the shower maxima, a proton-air cross section of [505±22(stat)(-36)(+28)(syst)] mb is found.

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×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×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.

Measurement of the depth of maximum of extensive air showers above 10{18} eV

We describe the measurement of the depth of maximum, X{max}, of the longitudinal development of air showers induced by cosmic rays. Almost 4000 events above 10;{18} eV observed by the fluorescence detector of the Pierre Auger Observatory in coincidence with at least one surface detector station are selected for the analysis. The average shower maximum was found to evolve with energy at a rate of (106{-21}{+35}) g/cm{2}/decade below 10{18.24+/-0.05} eV, and (24+/-3) g/cm{2}/decade above this energy. The measured shower-to-shower fluctuations decrease from about 55 to 26 g/cm{2}. The interpretation of these results in terms of the cosmic ray mass composition is briefly discussed.

Observation of the suppression of the flux of cosmic rays above 4 x 10 (19) eV

The energy spectrum of cosmic rays above 2.5 x 10;{18} eV, derived from 20,000 events recorded at the Pierre Auger Observatory, is described. The spectral index gamma of the particle flux, J proportional, variantE;{-gamma}, at energies between 4 x 10;{18} eV and 4 x 10;{19} eV is 2.69+/-0.02(stat)+/-0.06(syst), steepening to 4.2+/-0.4(stat)+/-0.06(syst) at higher energies. The hypothesis of a single power law is rejected with a significance greater than 6 standard deviations. The data are consistent with the prediction by Greisen and by Zatsepin and Kuz'min.

Upper limit on the diffuse flux of ultrahigh energy tau neutrinos from the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory is sensitive to Earth-skimming tau neutrinos that interact in Earth's crust. Tau leptons from nu(tau) charged-current interactions can emerge and decay in the atmosphere to produce a nearly horizontal shower with a significant electromagnetic component. The data collected between 1 January 2004 and 31 August 2007 are used to place an upper limit on the diffuse flux of nu(tau) at EeV energies. Assuming an E(nu)(-2) differential energy spectrum the limit set at 90% C.L. is E(nu)(2)dN(nu)(tau)/dE(nu)<1.3 x 10(-7) GeV cm(-2) s(-1) sr(-1) in the energy range 2 x 10(17) eV< E(nu)< 2 x 10(19) eV.

Detection and imaging of atmospheric radio flashes from cosmic ray air showers

The nature of ultrahigh-energy cosmic rays (UHECRs) at energies >10(20) eV remains a mystery. They are likely to be of extragalactic origin, but should be absorbed within approximately 50 Mpc through interactions with the cosmic microwave background. As there are no sufficiently powerful accelerators within this distance from the Galaxy, explanations for UHECRs range from unusual astrophysical sources to exotic string physics. Also unclear is whether UHECRs consist of protons, heavy nuclei, neutrinos or gamma-rays. To resolve these questions, larger detectors with higher duty cycles and which combine multiple detection techniques are needed. Radio emission from UHECRs, on the other hand, is unaffected by attenuation, has a high duty cycle, gives calorimetric measurements and provides high directional accuracy. Here we report the detection of radio flashes from cosmic-ray air showers using low-cost digital radio receivers. We show that the radiation can be understood in terms of the geosynchrotron effect. Our results show that it should be possible to determine the nature and composition of UHECRs with combined radio and particle detectors, and to detect the ultrahigh-energy neutrinos expected from flavour mixing.