Source author record

F. G. Schröder

F. G. Schröder 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

86works
4topics
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

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

preprint2025arXiv

All-sky neutrino point-source search with IceCube combined track and cascade data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 04/06/2008 to 05/23/2022 to search the whole sky for neutrino sources and, using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing-power and sensitivity in the northern sky, with all-sky cascades, with good energy-resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube which targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5$σ$ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068's emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

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

Density of GeV muons in air showers measured with IceTop

We present a measurement of the density of GeV muons in near-vertical air showers using three years of data recorded by the IceTop array at the South Pole. Depending on the shower size, the muon densities have been measured at lateral distances between 200 m and 1000 m. From these lateral distributions, we derive the muon densities as functions of energy at reference distances of 600 m and 800 m for primary energies between 2.5 PeV and 40 PeV and between 9 PeV and 120 PeV, respectively. The muon densities are determined using, as a baseline, the hadronic interaction model Sibyll 2.1 together with various composition models. The measurements are consistent with the predicted muon densities within these baseline interaction and composition models. The measured muon densities have also been compared to simulations using the post-LHC models EPOS-LHC and QGSJet-II.04. The result of this comparison is that the post-LHC models together with any given composition model yield higher muon densities than observed. This is in contrast to the observations above 1 EeV where all model simulations yield for any mass composition lower muon densities than the measured ones. The post-LHC models in general feature higher muon densities so that the agreement with experimental data at the highest energies is improved but the muon densities are not correct in the energy range between 2.5 PeV and about 100 PeV.

preprint2022arXiv

First all-flavor search for transient neutrino emission using 3-years of IceCube DeepCore data

Since the discovery of a flux of high-energy astrophysical neutrinos, searches for their origins have focused primarily at TeV-PeV energies. Compared to sub-TeV searches, high-energy searches benefit from an increase in the neutrino cross section, improved angular resolution on the neutrino direction, and a reduced background from atmospheric neutrinos and muons. However, the focus on high energy does not preclude the existence of sub-TeV neutrino emission where IceCube retains sensitivity. Here we present the first all-flavor search from IceCube for transient emission of low-energy neutrinos, focusing on the energy region of 5.6-100 GeV using three years of data obtained with the IceCube-DeepCore detector. We find no evidence of transient neutrino emission in the data, thus leading to a constraint on the volumetric rate of astrophysical transient sources in the range of $\sim 705-2301\, \text{Gpc}^{-3}\, \text{yr}^{-1}$ for sources following a subphotospheric energy spectrum with a mean energy of 100 GeV and a bolometric energy of $10^{52}$ erg.

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 Relativistic Magnetic Monopoles with Eight Years of IceCube Data

We present an all-sky 90\% confidence level upper limit on the cosmic flux of relativistic magnetic monopoles using 2886 days of IceCube data. The analysis was optimized for monopole speeds between 0.750$c$ and 0.995$c$, without any explicit restriction on the monopole mass. We constrain the flux of relativistic cosmic magnetic monopoles to a level below $2.0\times 10^{-19} {\textrm{cm}}^{-2} {\textrm{s}}^{-1} {\textrm{sr}}^{-1}$ over the majority of the targeted speed range. This result constitutes the most strict upper limit to date for magnetic monopoles above the Cherenkov threshold and up to $β\sim 0.995$ and fills the gap between existing limits on the cosmic flux of non-relativistic and ultrarelativistic magnetic monopoles

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.

preprint2021arXiv

A search for neutrino emission from cores of Active Galactic Nuclei

The sources of the majority of the high-energy astrophysical neutrinos observed with the IceCube neutrino telescope at the South Pole are unknown. So far, only a gamma-ray blazar was compellingly associated with the emission of high-energy neutrinos. In addition, several studies suggest that the neutrino emission from the gamma-ray blazar population only accounts for a small fraction of the total astrophysical neutrino flux. In this work we probe the production of high-energy neutrinos in the cores of Active Galactic Nuclei (AGN), induced by accelerated cosmic rays in the accretion disk region. We present a likelihood analysis based on eight years of IceCube data, searching for a cumulative neutrino signal from three AGN samples created for this work. The neutrino emission is assumed to be proportional to the accretion disk luminosity estimated from the soft X-ray flux. Next to the observed soft X-ray flux, the objects for the three samples have been selected based on their radio emission and infrared color properties. For the largest sample in this search, an excess of high-energy neutrino events with respect to an isotropic background of atmospheric and astrophysical neutrinos is found, corresponding to a post-trial significance of 2.60$σ$. Assuming a power-law spectrum, the best-fit spectral index is $2.03^{+0.14}_{-0.11}$, consistent with expectations from particle acceleration in astrophysical sources. If interpreted as a genuine signal with the assumptions of a proportionality of X-ray and neutrino fluxes and a model for the sub-threshold flux distribution, this observation implies that at 100 TeV, 27$\%$ - 100$\%$ of the observed neutrinos arise from particle acceleration in the core of AGN.

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

Detection of a particle shower at the Glashow resonance with IceCube

The Glashow resonance describes the resonant formation of a $W^-$ boson during the interaction of a high-energy electron antineutrino with an electron, peaking at an antineutrino energy of 6.3 petaelectronvolts (PeV) in the rest frame of the electron. Whereas this energy scale is out of reach for currently operating and future planned particle accelerators, natural astrophysical phenomena are expected to produce antineutrinos with energies beyond the PeV scale. Here we report the detection by the IceCube neutrino observatory of a cascade of high-energy particles (a particle shower) consistent with being created at the Glashow resonance. A shower with an energy of $6.05 \pm 0.72$ PeV (determined from Cherenkov radiation in the Antarctic Ice Sheet) was measured. Features consistent with the production of secondary muons in the particle shower indicate the hadronic decay of a resonant $W^-$ boson, confirm that the source is astrophysical and provide improved directional localization. The evidence of the Glashow resonance suggests the presence of electron antineutrinos in the astrophysical flux, while also providing further validation of the standard model of particle physics. Its unique signature indicates a method of distinguishing neutrinos from antineutrinos, thus providing a way to identify astronomical accelerators that produce neutrinos via hadronuclear or photohadronic interactions, with or without strong magnetic fields. As such, knowledge of both the flavour (that is, electron, muon or tau neutrinos) and charge (neutrino or antineutrino) will facilitate the advancement of neutrino astronomy.

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.

preprint2021arXiv

IceCube Data for Neutrino Point-Source Searches Years 2008-2018

IceCube has performed several all-sky searches for point-like neutrino sources using track-like events, including a recent time-integrated analysis using 10 years of IceCube data. This paper accompanies the public data release of these neutrino candidates detected by IceCube between April 6, 2008 and July 8, 2018. The selection includes through-going tracks, primarily due to muon neutrino candidates, that reach the detector from all directions, as well as neutrino track events that start within the instrumented volume. An updated selection and reconstruction for data taken after April 2012 slightly improves the sensitivity of the sample. While more than 80% of the sample overlaps between the old and new versions, differing events can lead to changes relative to the previous 7 year event selection. An a posteriori estimate of the significance of the 2014-2015 TXS flare is reported with an explanation of observed discrepancies with previous results. This public data release, which includes 10 years of data and binned detector response functions for muon neutrino signal events, shows improved sensitivity in generic time-integrated point source analyses and should be preferred over previous releases.

preprint2021arXiv

Improved Characterization of the Astrophysical Muon-Neutrino Flux with 9.5 Years of IceCube Data

We present a measurement of the high-energy astrophysical muon-neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of ~650k neutrino-induced muon tracks from the Northern celestial hemisphere, corresponding to 9.5 years of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fluxes has been implemented based on recent models. The best-fit single power-law parameterization for the astrophysical energy spectrum results in a normalization of $ϕ_{\mathrm{@100TeV}}^{ν_μ+\barν_μ} = 1.44_{-0.26}^{+0.25} \times 10^{-18}\,\mathrm{GeV}^{-1}\mathrm{cm}^{-2}\mathrm{s}^{-1}\mathrm{sr}^{-1}$ and a spectral index $γ_{\mathrm{SPL}} = 2.37_{-0.09}^{+0.09}$, constrained in the energy range from $15\,\mathrm{TeV}$ to $5\,\mathrm{PeV}$. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class specific flux predictions from the literature and a model-independent spectral unfolding. The data is well consistent with a single power law hypothesis, however, spectra with softening above one PeV are statistically more favorable at a two sigma level.

preprint2021arXiv

Search for High-Energy Neutrinos from Ultra-Luminous Infrared Galaxies with IceCube

Ultra-luminous infrared galaxies (ULIRGs) have infrared luminosities $L_{\mathrm{IR}} \geq 10^{12} L_{\odot}$, making them the most luminous objects in the infrared sky. These dusty objects are generally powered by starbursts with star-formation rates that exceed $100~ M_{\odot}~ \mathrm{yr}^{-1}$, possibly combined with a contribution from an active galactic nucleus. Such environments make ULIRGs plausible sources of astrophysical high-energy neutrinos, which can be observed by the IceCube Neutrino Observatory at the South Pole. We present a stacking search for high-energy neutrinos from a representative sample of 75 ULIRGs with redshift $z \leq 0.13$ using 7.5 years of IceCube data. The results are consistent with a background-only observation, yielding upper limits on the neutrino flux from these 75 ULIRGs. For an unbroken $E^{-2.5}$ power-law spectrum, we report an upper limit on the stacked flux $Φ_{ν_μ+ \barν_μ}^{90\%} = 3.24 \times 10^{-14}~ \mathrm{TeV^{-1}~ cm^{-2}~ s^{-1}}~ (E/10~ \mathrm{TeV})^{-2.5}$ at 90% confidence level. In addition, we constrain the contribution of the ULIRG source population to the observed diffuse astrophysical neutrino flux as well as model predictions.

preprint2021arXiv

Searches for neutrinos from cosmic-ray interactions in the Sun using seven years of IceCube data

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SA$ν$s) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SA$ν$s also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. Here, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SA$ν$ flux is observed. After inclusion of systematic uncertainties, we set a 90\% upper limit of $1.02^{+0.20}_{-0.18}\cdot10^{-13}$~$\mathrm{GeV^{-1}cm^{-2}s^{-1}}$ at 1 TeV.

preprint2020arXiv

A search for IceCube events in the direction of ANITA neutrino candidates

During the first three flights of the Antarctic Impulsive Transient Antenna (ANITA) experiment, the collaboration detected several neutrino candidates. Two of these candidate events were consistent with an ultra-high-energy up-going air shower and compatible with a tau neutrino interpretation. A third neutrino candidate event was detected in a search for Askaryan radiation in the Antarctic ice, although it is also consistent with the background expectation. The inferred emergence angle of the first two events is in tension with IceCube and ANITA limits on isotropic cosmogenic neutrino fluxes. Here, we test the hypothesis that these events are astrophysical in origin, possibly caused by a point source in the reconstructed direction. Given that any ultra-high-energy tau neutrino flux traversing the Earth should be accompanied by a secondary flux in the TeV-PeV range, we search for these secondary counterparts in seven years of IceCube data using three complementary approaches. In the absence of any significant detection, we set upper limits on the neutrino flux from potential point sources. We compare these limits to ANITA's sensitivity in the same direction and show that an astrophysical explanation of these anomalous events under standard model assumptions is severely constrained regardless of source spectrum.

preprint2020arXiv

A Search for MeV to TeV Neutrinos from Fast Radio Bursts with IceCube

We present two searches for IceCube neutrino events coincident with 28 fast radio bursts (FRBs) and one repeating FRB. The first improves upon a previous IceCube analysis -- searching for spatial and temporal correlation of events with FRBs at energies greater than roughly 50 GeV -- by increasing the effective area by an order of magnitude. The second is a search for temporal correlation of MeV neutrino events with FRBs. No significant correlation is found in either search, therefore, we set upper limits on the time-integrated neutrino flux emitted by FRBs for a range of emission timescales less than one day. These are the first limits on FRB neutrino emission at the MeV scale, and the limits set at higher energies are an order-of-magnitude improvement over those set by any neutrino telescope.

preprint2020arXiv

ANTARES and IceCube Combined Search for Neutrino Point-like and Extended Sources in the Southern Sky

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

preprint2020arXiv

Characteristics of the diffuse astrophysical electron and tau neutrino flux with six years of IceCube high energy cascade data

We report on the first measurement of the astrophysical neutrino flux using particle showers (cascades) in IceCube data from 2010 -- 2015. Assuming standard oscillations, the astrophysical neutrinos in this dedicated cascade sample are dominated ($\sim 90 \%$) by electron and tau flavors. The flux, observed in the sensitive energy range from $16\,\mathrm{TeV}$ to $2.6\,\mathrm{PeV}$, is consistent with a single power-law model as expected from Fermi-type acceleration of high energy particles at astrophysical sources. We find the flux spectral index to be $γ=2.53\pm0.07$ and a flux normalization for each neutrino flavor of $ϕ_{astro} = 1.66^{+0.25}_{-0.27}$ at $E_{0} = 100\, \mathrm{TeV}$, in agreement with IceCube's complementary muon neutrino results and with all-neutrino flavor fit results. In the measured energy range we reject spectral indices $γ\leq2.28$ at $\ge3σ$ significance level. Due to high neutrino energy resolution and low atmospheric neutrino backgrounds, this analysis provides the most detailed characterization of the neutrino flux at energies below $\sim100\,{\rm{TeV}}$ compared to previous IceCube results. Results from fits assuming more complex neutrino flux models suggest a flux softening at high energies and a flux hardening at low energies (p-value $\ge 0.06$). The sizable and smooth flux measured below $\sim 100\,{\rm{TeV}}$ remains a puzzle. In order to not violate the isotropic diffuse gamma-ray background as measured by the Fermi-LAT, it suggests the existence of astrophysical neutrino sources characterized by dense environments which are opaque to gamma-rays.

preprint2020arXiv

Constraints on Neutrino Emission from Nearby Galaxies Using the 2MASS Redshift Survey and IceCube

The distribution of galaxies within the local universe is characterized by anisotropic features. Observatories searching for the production sites of astrophysical neutrinos can take advantage of these features to establish directional correlations between a neutrino dataset and overdensities in the galaxy distribution in the sky. The results of two correlation searches between a seven-year time-integrated neutrino dataset from the IceCube Neutrino Observatory, and the 2MASS Redshift Survey (2MRS) catalog are presented here. The first analysis searches for neutrinos produced via interactions between diffuse intergalactic Ultra-High Energy Cosmic Rays (UHECRs) and the matter contained within galaxies. The second analysis searches for low-luminosity sources within the local universe, which would produce subthreshold multiplets in the IceCube dataset that directionally correlate with galaxy distribution. No significant correlations were observed in either analyses. Constraints are presented on the flux of neutrinos originating within the local universe through diffuse intergalactic UHECR interactions, as well as on the density of standard candle sources of neutrinos at low luminosities.

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

IceCube Search for High-Energy Neutrino Emission from TeV Pulsar Wind Nebulae

Pulsar wind nebulae (PWNe) are the main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit nonthermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 yr of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

preprint2020arXiv

IceCube Search for Neutrinos Coincident with Compact Binary Mergers from LIGO-Virgo's First Gravitational-Wave Transient Catalog

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational wave (GW) detectors during their first and second observing runs. We present results from two searches targeting emission coincident with the sky localization of each gravitational wave event within a 1000 second time window centered around the reported merger time. One search uses a model-independent unbinned maximum likelihood analysis, which uses neutrino data from IceCube to search for point-like neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 second window for each of the 11 GW events. These limits range from 0.02-0.7 $\mathrm{GeV~cm^{-2}}$. We also set limits on the total isotropic equivalent energy, $E_{\mathrm{iso}}$, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 $\times$ 10$^{51}$ - 1.8 $\times$ 10$^{55}$ erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

preprint2020arXiv

Measurements of the Time-Dependent Cosmic-Ray Sun Shadow with Seven Years of IceCube Data -- Comparison with the Solar Cycle and Magnetic Field Models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3σ$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

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

Search for PeV Gamma-Ray Emission from the Southern Hemisphere with 5 Years of Data from the IceCube Observatory

The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a point-like source of PeV gamma rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several un-binned maximum likelihood searches for PeV gamma rays in the Southern Hemisphere using 5 years of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers, and provides excellent sensitivity to gamma rays between $\sim$0.6 PeV and 100 PeV. Our measurements of point-like and diffuse Galactic emission of PeV gamma rays are consistent with background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic Plane at 2 PeV to $2.61 \times 10^{-19}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ at 90% confidence, assuming an E$^{-3}$ spectrum, and we estimate 90% upper limits on point-like emission at 2 PeV between 10$^{-21}$ - 10$^{-20}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ for an E$^{-2}$ spectrum, depending on declination. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by H.E.S.S., and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma rays correlated with neutrinos from IceCube's high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

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.

preprint2019arXiv

A Search for Neutrino Point-Source Populations in 7 Years of IceCube Data with Neutrino-count Statistics

The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

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.

preprint2019arXiv

Efficient propagation of systematic uncertainties from calibration to analysis with the SnowStorm method in IceCube

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. This method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.

preprint2019arXiv

Search for Sources of Astrophysical Neutrinos Using Seven Years of IceCube Cascade Events

Low background searches for astrophysical neutrino sources anywhere in the sky can be performed using cascade events induced by neutrinos of all flavors interacting in IceCube with energies as low as ~1 TeV. Previously, we showed that even with just two years of data, the resulting sensitivity to sources in the southern sky is competitive with IceCube and ANTARES analyses using muon tracks induced by charge current muon neutrino interactions - especially if the neutrino emission follows a soft energy spectrum or originates from an extended angular region. Here, we extend that work by adding five more years of data, significantly improving the cascade angular resolution, and including tests for point-like or diffuse Galactic emission to which this dataset is particularly well-suited. For many of the signal candidates considered, this analysis is the most sensitive of any experiment. No significant clustering was observed, and thus many of the resulting constraints are the most stringent to date. In this paper we will describe the improvements introduced in this analysis and discuss our results in the context of other recent work in neutrino astronomy.

preprint2019arXiv

Time-integrated Neutrino Source Searches with 10 years of IceCube Data

This paper presents the results from point-like neutrino source searches using ten years of IceCube data collected between Apr.~6, 2008 and Jul.~10, 2018. We evaluate the significance of an astrophysical signal from a point-like source looking for an excess of clustered neutrino events with energies typically above $\sim1\,$TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the Northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of $2.9\,σ$ after accounting for statistical trials. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the Northern catalog are inconsistent with background at 3.3$\,σ$ significance. These results, all based on searches for a cumulative neutrino signal integrated over the ten years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.

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

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

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

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

Reconstruction of air-shower parameters for large-scale radio detectors using the lateral distribution

We investigate features of the lateral distribution function (LDF) of the radio signal emitted by cosmic ray air-showers with primary energies $> 0.1$~EeV and its connection to air-shower parameters such as energy and shower maximum using CoREAS simulations made for the configuration of the Tunka-Rex antenna array. Taking into account all significant contributions to the total radio emission, such as by the geomagnetic effect, the charge excess, and the atmospheric refraction we parameterize the radio LDF. This parameterization is two-dimensional and has several free parameters. The large number of free parameters is not suitable for experiments of sparse arrays operating at low SNR (signal-to-noise ratios). Thus, exploiting symmetries, we decrease the number of free parameters and reduce the LDF to a simple one-dimensional function. The remaining parameters can be fit with a small number of points, i.e. as few as the signal from three antennas above detection threshold. Finally, we present a method for the reconstruction of air-shower parameters, in particular, energy and $X_{\mathrm{max}}$ (shower maximum), which can be reached with a theoretical accuracy of better than 15\% and 30~g/cm$^2$, respectively.

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

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

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

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

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

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

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: a Radio Antenna Array for the Tunka Experiment (ARENA 2012)

Tunka-Rex, the Tunka radio extension, is an array of 20 antennas at the Tunka experiment close to Lake Baikal in Siberia. It started operation on 08 October 2012. The antennas are connected directly to the data acquisition of the Tunka main detector, a 1 square-km large array of 133 non-imaging photomultipliers observing the Cherenkov light of air showers in dark and clear nights. This allows to cross-calibrate the radio signal with the air-Cherenkov signal of the same air showers - in particular with respect to the energy and the atmospheric depth of the shower maximum, Xmax. Consequently, we can test whether in rural regions with low radio background the practically achievable radio precision comes close to the precision of the established fluorescence and air-Cherenkov techniques. At a mid-term perspective, due to its higher duty-cycle, Tunka-Rex can enhance the effective observing time of Tunka by an order of magnitude, at least in the interesting energy range above 100 PeV. Moreover, Tunka-Rex is very cost-effective, e.g., by using economic Short Aperiodic Loaded Loop Antennas (SALLAs). Thus, the results of Tunka-Rex and the comparison to other sophisticated radio arrays will provide crucial input for future large-scale cosmic-ray observatories, for which measurement precision as well as costs per area have to be optimized. In this paper we shortly describe the Tunka-Rex setup and discuss the technical and scientific goals of Tunka-Rex.

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

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.

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

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.