Source author record

S. Buitink

S. Buitink 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

96works
7topics
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

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

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

Determining Atmospheric Electric Fields using MGMR3D

Cosmic-ray particles impinging on the atmosphere induce high-energy particle cascades in air, an Extensive Air Shower (EAS), emitting coherent radio emission. This emission is affected by the presence of strong electric fields during thunderstorm conditions. To reconstruct the atmospheric electric field from the measured radio footprint of the EAS we use an analytic model for the calculation of the radio emission, MGMR3D. In this work we make an extensive comparison between the results of a microscopic model for radio emission, CoREAS, to obtain an improved parametrization for MGMR3D in the presence of atmospheric electric fields, as well as confidence intervals. The approach to extract the electric field structure is applied successfully to an event with a complicated radio footprint measured by LOFAR during thunderstorm conditions. This shows that, with the improved parametrization, MGMR3D can be used to extract the structure of the atmospheric electric field.

preprint2022arXiv

In situ, broadband measurement of the radio frequency attenuation length at Summit Station, Greenland

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or $10^{17}$ electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_α$. We find an approximately linear dependence of $L_α$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_α\rangle = \big( (1154 \pm 121) - (0.81 \pm 0.14) (ν/$MHz$)\big)$ m for frequencies $ν\in [145 - 350]$ MHz.

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

Reconstructing the neutrino energy for in-ice radio detectors

Starting in summer 2021, the Radio Neutrino Observatory in Greenland (RNO-G) will search for astrophysical neutrinos at energies >10 PeV by detecting the radio emission from particle showers in the ice around Summit Station, Greenland. We present an extensive simulation study that shows how RNO-G will be able to measure the energy of such particle cascades, which will in turn be used to estimate the energy of the incoming neutrino that caused them. The location of the neutrino interaction is determined using the differences in arrival times between channels and the electric field of the radio signal is reconstructed using a novel approach based on Information Field Theory. Based on these properties, the shower energy can be estimated. We show that this method can achieve an uncertainty of 13% on the logarithm of the shower energy after modest quality cuts and estimate how this can constrain the energy of the neutrino. The method presented in this paper is applicable to all similar radio neutrino detectors, such as the proposed radio array of IceCube-Gen2.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

Triboelectric Backgrounds to radio-based UHE Neutrino Exeperiments

The proposed IceCube-Gen2 (ICG2) seeks to instrument ~500 sq. km of Antarctic ice near the geographic South Pole with radio antennas, in order to observe the highest energy (E>1 EeV) neutrinos in the Universe. To this end, ICG2 will use the impulsive radio-frequency (RF) signal produced by neutrino interactions in polar ice caps. In such experiments, rare single event candidates must be unambiguously separated from background; to date, signal identification strategies primarily reject thermal noise and anthropogenic backgrounds. Here, we consider the possibility that fake neutrino signals may also be naturally generated via the 'triboelectric effect'. This broadly includes any process in which force applied at a boundary layer results in displacement of surface charge, generating a potential difference ΔV. Wind blowing over granular surfaces such as snow can induce such a ΔV, with subsequent discharge. Discharges over nanosecond-timescales can then lead to RF emissions at characteristic MHz-GHz frequencies. We find that such backgrounds are evident in the several neutrino experiments considered, and are generally characterized by: a) a threshold wind velocity which likely depends on the experimental signal trigger threshold and layout; for the experiments considered herein, this value is typically O(10 m/s), b) frequency spectra generally shifted to the low-end of the frequency regime to which current radio experiments are typically sensitive (100-200 MHz), c) for the strongest background signals, an apparent preference for discharges from above-surface structures, although the presence of more isotropic, lower amplitude triboelectric discharges cannot be excluded.

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2017arXiv

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

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

preprint2016arXiv

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

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

preprint2016arXiv

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

Cosmic-ray energy spectrum and composition up to the ankle - the case for a second Galactic component

We have carried out a detailed study to understand the observed energy spectrum and composition of cosmic rays with energies up to ~10^18 eV. Our study shows that a single Galactic component with subsequent energy cut-offs in the individual spectra of different elements, optimised to explain the observed spectra below ~10^14 eV and the knee in the all-particle spectrum, cannot explain the observed all-particle spectrum above ~2x10^16 eV. We discuss two approaches for a second component of Galactic cosmic rays -- re-acceleration at a Galactic wind termination shock, and supernova explosions of Wolf-Rayet stars, and show that the latter scenario can explain almost all observed features in the all-particle spectrum and the composition up to ~10^18 eV, when combined with a canonical extra-galactic spectrum expected from strong radio galaxies or a source population with similar cosmological evolution. In this two-component Galactic model, the knee at ~ 3x10^15 eV and the second knee at ~10^17 eV in the all-particle spectrum are due to the cut-offs in the first and second components, respectively. We also discuss several variations of the extra-galactic component, from a minimal contribution to scenarios with a significant component below the ankle (at ~4x10^18 eV), and find that extra-galactic contributions in excess of regular source evolution are neither indicated nor in conflict with the existing data. Our main result is that the second Galactic component predicts a composition of Galactic cosmic rays at and above the second knee that largely consists of helium or a mixture of helium and CNO nuclei, with a weak or essentially vanishing iron fraction, in contrast to most common assumptions. This prediction is in agreement with new measurements from LOFAR and the Pierre Auger Observatory which indicate a strong light component and a rather low iron fraction between ~10^17 and 10^18 eV.

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

Measurement of the circular polarization in radio emission from extensive air showers confirms emission mechanisms

We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields.

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

Obscured flat spectrum radio AGN as sources of high-energy neutrinos

Active Galactic Nuclei (AGN) are believed to be one of the main source candidates for the high-energy (TeV-PeV) cosmic neutrino flux recently discovered by the IceCube neutrino observatory. Nevertheless, several correlation studies between AGN and the cosmic neutrinos detected by IceCube show no significance. Therefore, in this article we consider a specific sub-class of AGN for which an increased neutrino production is expected. This sub-class contains AGN for which their high-energy jet is pointing toward Earth. Furthermore, we impose the condition that the jet is obscured by gas or dust surrounding the AGN. A method is presented to determine the total column density of the obscuring medium, which is probed by determining the relative X-ray attenuation with respect to the radio flux as obtained from the AGN spectrum. The total column density allows us to probe the interaction of the jet with the surrounding matter which leads to additional neutrino production. Finally, starting from two different source catalogs, this method is applied to specify a sample of low redshift radio galaxies for which an increased neutrino production is expected.

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

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 high-energy neutrinos from dust obscured Blazars

The recent discovery of high-energy cosmic neutrinos by the IceCube neutrino observatory opens up a new field in physics, the field of neutrino astronomy. Using the IceCube neutrino detector we plan to search for high-energy neutrinos emitted from Active Galactic Nuclei (AGN), since AGN are believed to be one of the most promising sources of the most energetic cosmic rays and hence of high-energy neutrinos. We discuss a specific type of AGN which we plan to investigate in more detail with data obtained by the IceCube observatory. The main properties of the AGN category in which we are interested are given by a high-energy jet which is pointing in our line of sight defining a class of AGN, called Blazars, and in particular the ones that are obscured by surrounding dust. The jet-matter interaction is expected to give an increased high-energy neutrino production. The properties of this specific type of AGN are expected to give very distinct features in the electromagnetic spectrum, which are discussed in detail.

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.

preprint2016arXiv

Timing calibration and spectral cleaning of LOFAR time series data

We describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns.

preprint2015arXiv

A First Search for Cosmogenic Neutrinos with the ARIANNA Hexagonal Radio Array

The ARIANNA experiment seeks to observe the diffuse flux of neutrinos in the 10^8 - 10^10 GeV energy range using a grid of radio detectors at the surface of the Ross Ice Shelf of Antarctica. The detector measures the coherent Cherenkov radiation produced at radio frequencies, from about 100 MHz to 1 GHz, by charged particle showers generated by neutrino interactions in the ice. The ARIANNA Hexagonal Radio Array (HRA) is being constructed as a prototype for the full array. During the 2013-14 austral summer, three HRA stations collected radio data which was wirelessly transmitted off site in nearly real-time. The performance of these stations is described and a simple analysis to search for neutrino signals is presented. The analysis employs a set of three cuts that reject background triggers while preserving 90% of simulated cosmogenic neutrino triggers. No neutrino candidates are found in the data and a model-independent 90% confidence level Neyman upper limit is placed on the all flavor neutrino+antineutrino flux in a sliding decade-wide energy bin. The limit reaches a minimum of 1.9x10^-23 GeV^-1 cm^-2 s^-1 sr^-1 in the 10^8.5 - 10^9.5 GeV energy bin. Simulations of the performance of the full detector are also described. The sensitivity of the full ARIANNA experiment is presented and compared with current neutrino flux models.

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

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

Influence of Atmospheric Electric Fields on the Radio Emission from Extensive Air Showers

The atmospheric electric fields in thunderclouds have been shown to significantly modify the intensity and polarization patterns of the radio footprint of cosmic-ray-induced extensive air showers. Simulations indicated a very non-linear dependence of the signal strength in the frequency window of 30-80 MHz on the magnitude of the atmospheric electric field. In this work we present an explanation of this dependence based on Monte-Carlo simulations, supported by arguments based on electron dynamics in air showers and expressed in terms of a simplified model. We show that by extending the frequency window to lower frequencies additional sensitivity to the atmospheric electric field is obtained.

preprint2015arXiv

Livetime and sensitivity of the ARIANNA Hexagonal Radio Array

The ARIANNA collaboration completed the installation of the hexagonal radio array (HRA) in December 2014, serving as a pilot program for a planned high energy neutrino telescope located about 110 km south of McMurdo Station on the Ross Ice Shelf near the coast of Antarctica. The goal of ARIANNA is to measure both diffuse and point fluxes of astrophysical neutrinos at energies in excess of 1016 eV. Upgraded hardware has been installed during the 2014 deployment season and stations show a livetime of better than 90% between commissioning and austral sunset. Though designed to observe radio pulses from neutrino interactions originating within the ice below each detector, one station was modified to study the low-frequency environment and signals from above. We provide evidence that the HRA observed both continuous emission from the Galaxy and a transient solar burst. Preliminary work on modeling the (weak) Galactic signal confirm the absolute sensitivity of the HRA detector system.

preprint2015arXiv

Measurement of the cosmic-ray energy spectrum above $10^{16}$ eV with the LOFAR Radboud Air Shower Array

The energy reconstruction of extensive air showers measured with the LOFAR Radboud Air Shower Array (LORA) is presented in detail. LORA is a particle detector array located in the center of the LOFAR radio telescope in the Netherlands. The aim of this work is to provide an accurate and independent energy measurement for the air showers measured through their radio signal with the LOFAR antennas. The energy reconstruction is performed using a parameterized relation between the measured shower size and the cosmic-ray energy obtained from air shower simulations. In order to illustrate the capabilities of LORA, the all-particle cosmic-ray energy spectrum has been reconstructed, assuming that cosmic rays are composed only of protons or iron nuclei in the energy range between $\sim2\times10^{16}$ and $2\times10^{18}$ eV. The results are compatible with literature values and a changing mass composition in the transition region from a galactic to an extragalactic origin of cosmic rays.

preprint2015arXiv

Performance of the ARIANNA Hexagonal Radio Array

Installation of the ARIANNA Hexagonal Radio Array (HRA) on the Ross Ice Shelf of Antarctica has been completed. This detector serves as a pilot program to the ARIANNA neutrino telescope, which aims to measure the diffuse flux of very high energy neutrinos by observing the radio pulse generated by neutrino-induced charged particle showers in the ice. All HRA stations ran reliably and took data during the entire 2014-2015 austral summer season. A new radio signal direction reconstruction procedure is described, and is observed to have a resolution better than a degree. The reconstruction is used in a preliminary search for potential neutrino candidate events in the data from one of the newly installed detector stations. Three cuts are used to separate radio backgrounds from neutrino signals. The cuts are found to filter out all data recorded by the station during the season while preserving 85.4% of simulated neutrino events that trigger the station. This efficiency is similar to that found in analyses of previous HRA data taking seasons.

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

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

preprint2015arXiv

Search for dark matter annihilations in the Sun with the 79-string IceCube detector

We have performed a search for muon neutrinos from dark matter annihilation in the center of the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are consistent with the expected background from atmospheric muons and neutrinos. Upper limits are therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV.

preprint2015arXiv

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

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

preprint2015arXiv

The lunar Askaryan technique: a technical roadmap

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

preprint2015arXiv

The radio emission pattern of air showers as measured with LOFAR - a tool for the reconstruction of the energy and the shower maximum

The pattern of the radio emission of air showers is finely sampled with the Low-Frequency ARray (LOFAR). A set of 382 measured air showers is used to test a fast, analytic parameterization of the distribution of pulse powers. Using this parameterization we are able to reconstruct the shower axis and give estimators for the energy of the air shower as well as the distance to the shower maximum.

preprint2014arXiv

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

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

preprint2014arXiv

LORA: A scintillator array for LOFAR to measure extensive air showers

The measurement of the radio emission from extensive air showers, induced by high-energy cosmic rays is one of the key science projects of the LOFAR radio telescope. The LOfar Radboud air shower Array (LORA) has been installed in the core of LOFAR in the Netherlands. The main purpose of LORA is to measure the properties of air showers and to trigger the read-out of the LOFAR radio antennas to register extensive air showers. The experimental set-up of the array of scintillation detectors and its performance are described.

preprint2014arXiv

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

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

preprint2014arXiv

Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

preprint2014arXiv

Polarized radio emission from extensive air showers measured with LOFAR

We present LOFAR measurements of radio emission from extensive air showers. We find that this emission is strongly polarized, with a median degree of polarization of nearly $99\%$, and that the angle between the polarization direction of the electric field and the Lorentz force acting on the particles, depends on the observer location in the shower plane. This can be understood as a superposition of the radially polarized charge-excess emission mechanism, first proposed by Askaryan and the geomagnetic emission mechanism proposed by Kahn and Lerche. We calculate the relative strengths of both contributions, as quantified by the charge-excess fraction, for $163$ individual air showers. We find that the measured charge-excess fraction is higher for air showers arriving from closer to the zenith. Furthermore, the measured charge-excess fraction also increases with increasing observer distance from the air shower symmetry axis. The measured values range from $(3.3\pm 1.0)\%$ for very inclined air showers at $25\, \mathrm{m}$ to $(20.3\pm 1.3)\%$ for almost vertical showers at $225\, \mathrm{m}$. Both dependencies are in qualitative agreement with theoretical predictions.

preprint2014arXiv

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

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

preprint2014arXiv

The shape of the radio wavefront of extensive air showers as measured with LOFAR

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

preprint2013arXiv

An Absence of Neutrinos Associated with Cosmic Ray Acceleration in Gamma-Ray Bursts

Gamma-Ray Bursts (GRBs) have been proposed as a leading candidate for acceleration of ultra high-energy cosmic rays, which would be accompanied by emission of TeV neutrinos produced in proton-photon interactions during acceleration in the GRB fireball. Two analyses using data from two years of the IceCube detector produced no evidence for this neutrino emission, placing strong constraints on models of neutrino and cosmic-ray production in these sources.

preprint2013arXiv

Detecting cosmic rays with the LOFAR radio telescope

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

preprint2013arXiv

Measurement of South Pole ice transparency with the IceCube LED calibration system

The IceCube Neutrino Observatory, approximately 1 km^3 in size, is now complete with 86 strings deployed in the Antarctic ice. IceCube detects the Cherenkov radiation emitted by charged particles passing through or created in the ice. To realize the full potential of the detector, the properties of light propagation in the ice in and around the detector must be well understood. This report presents a new method of fitting the model of light propagation in the ice to a data set of in-situ light source events collected with IceCube. The resulting set of derived parameters, namely the measured values of scattering and absorption coefficients vs. depth, is presented and a comparison of IceCube data with simulations based on the new model is shown.

preprint2013arXiv

Measurement of the Atmospheric $ν_e$ flux in IceCube

We report the first measurement of the atmospheric electron neutrino flux in the energy range between approximately 80 GeV and 6 TeV, using data recorded during the first year of operation of IceCube's DeepCore low energy extension. Techniques to identify neutrinos interacting within the DeepCore volume and veto muons originating outside the detector are demonstrated. A sample of 1029 events is observed in 281 days of data, of which 496 $\pm$ 66(stat.) $\pm$ 88(syst.) are estimated to be cascade events, including both electron neutrino and neutral current events. The rest of the sample includes residual backgrounds due to atmospheric muons and charged current interactions of atmospheric muon neutrinos. The flux of the atmospheric electron neutrinos is consistent with models of atmospheric neutrinos in this energy range. This constitutes the first observation of electron neutrinos and neutral current interactions in a very large volume neutrino telescope optimized for the TeV energy range.

preprint2013arXiv

Search for Neutrinos from Annihilating Dark Matter in the Direction of the Galactic Center with the 40-String IceCube Neutrino Observatory

A search for muon neutrinos from dark matter annihilations in the Galactic Center region has been performed with the 40-string configuration of the IceCube Neutrino Observatory using data collected in 367 days of live-time starting in April 2008. The observed fluxes were consistent with the atmospheric background expectations. Upper limits on the self-annihilation cross-section are obtained for dark matter particle masses ranging from 100 GeV to 10 TeV. In the case of decaying dark matter, lower limits on the lifetime have been determined for masses between 200 GeV and 20 TeV.

preprint2013arXiv

Searches for high-energy neutrino emission in the Galaxy with the combined IceCube-AMANDA detector

We report on searches for neutrino sources at energies above 200 GeV in the Northern sky of the galactic plane, using the data collected by the South Pole neutrino telescopes IceCube and AMANDA. The galactic region considered here includes the Local Arm towards the Cygnus region and our closest approach to the Perseus Arm. The data have been collected between 2007 and 2009 when AMANDA was an integrated part of IceCube, which was still under construction and operated with 22-strings (2007-8) and 40-strings (2008-9) of optical modules deployed in the ice. By combining the larger IceCube detector with the lower energy threshold of the more compact AMANDA detector, we obtain an improved sensitivity at energies below $\sim$10 TeV with respect to previous searches. The analyses presented here are: a scan for point sources within the galactic plane; a search optimized for multiple and extended sources in the Cygnus region, which might be below the sensitivity of the point source scan; and studies of seven pre-selected neutrino source candidates. For one of them, Cygnus X-3, a time-dependent search for neutrinos in coincidence with observed radio and X-ray flares has been performed. No evidence of a signal is found, and upper limits are reported for each of the searches. We investigate neutrino spectra proportional to E$^{-2}$ and E$^{-3}$ to cover the entire range of possible spectra. The soft E$^{-3}$ spectrum results in an energy distribution similar to a source with cut-off below $\sim$50 TeV. For the considered region of the galactic plane, the 90% confidence level muon neutrino flux upper limits are in the range E$^3$dN/dE$\sim 5.4 - 19.5 \times 10^{-11} \rm{TeV^{2} cm^{-2} s^{-1}}$ for point-like neutrino sources in the energy region [180.0 GeV - 20.5 TeV]. These represent the most stringent upper limits for soft-spectra neutrino sources within the Galaxy reported to date.

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.

preprint2012arXiv

A Search for UHE Tau Neutrinos with IceCube

The first dedicated search for ultra-high energy (UHE) tau neutrinos of astrophysical origin was performed using the IceCube detector in its 22-string configuration with an instrumented volume of roughly 0.25 km^3. The search also had sensitivity to UHE electron and muon neutrinos. After application of all selection criteria to approximately 200 live-days of data, we expect a background of 0.60 +/- 0.19 (stat.) $^{+0.56}_{-0.58}$ (syst.) events and observe three events, which after inspection emerge as being compatible with background but are kept in the final sample. Therefore, we set an upper limit on neutrinos of all-flavors from UHE astrophysical sources at 90% CL of $E^{2} Φ(ν_{x}) < 16.3 * 10^-8 GeV cm^-2 sr^-1 s^-1 over an estimated primary neutrino energy range of 340 TeV to 200 PeV.

preprint2012arXiv

All-particle cosmic ray energy spectrum measured with 26 IceTop stations

We report on a measurement of the cosmic ray energy spectrum with the IceTop air shower array, the surface component of the IceCube Neutrino Observatory at the South Pole. The data used in this analysis were taken between June and October, 2007, with 26 surface stations operational at that time, corresponding to about one third of the final array. The fiducial area used in this analysis was 0.122 km^2. The analysis investigated the energy spectrum from 1 to 100 PeV measured for three different zenith angle ranges between 0° and 46°. Because of the isotropy of cosmic rays in this energy range the spectra from all zenith angle intervals have to agree. The cosmic-ray energy spectrum was determined under different assumptions on the primary mass composition. Good agreement of spectra in the three zenith angle ranges was found for the assumption of pure proton and a simple two-component model. For zenith angles θ < 30°, where the mass dependence is smallest, the knee in the cosmic ray energy spectrum was observed between 3.5 and 4.32 PeV, depending on composition assumption. Spectral indices above the knee range from -3.08 to -3.11 depending on primary mass composition assumption. Moreover, an indication of a flattening of the spectrum above 22 PeV were observed.

preprint2012arXiv

An improved method for measuring muon energy using the truncated mean of dE/dx

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the muon energy. Because high-energy muons (E_mu > 1 TeV) lose energy randomly, the spread in dE/dx values is quite large, leading to a typical energy resolution of 0.29 in log10(E_mu) for a muon observed over a 1 km path length in the IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the muon energy. This method results in an energy resolution of 0.22 in log10(E_mu), which gives a 26% improvement. This technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors.

preprint2012arXiv

Cosmic Ray Composition and Energy Spectrum from 1-30 PeV Using the 40-String Configuration of IceTop and IceCube

The mass composition of high energy cosmic rays depends on their production, acceleration, and propagation. The study of cosmic ray composition can therefore reveal hints of the origin of these particles. At the South Pole, the IceCube Neutrino Observatory is capable of measuring two components of cosmic ray air showers in coincidence: the electromagnetic component at high altitude (2835 m) using the IceTop surface array, and the muonic component above ~1 TeV using the IceCube array. This unique detector arrangement provides an opportunity for precision measurements of the cosmic ray energy spectrum and composition in the region of the knee and beyond. We present the results of a neural network analysis technique to study the cosmic ray composition and the energy spectrum from 1 PeV to 30 PeV using data recorded using the 40-string/40-station configuration of the IceCube Neutrino Observatory.

preprint2012arXiv

IceTop: The surface component of IceCube

IceTop, the surface component of the IceCube Neutrino Observatory at the South Pole, is an air shower array with an area of 1 km2. The detector allows a detailed exploration of the mass composition of primary cosmic rays in the energy range from about 100 TeV to 1 EeV by exploiting the correlation between the shower energy measured in IceTop and the energy deposited by muons in the deep ice. In this paper we report on the technical design, construction and installation, the trigger and data acquisition systems as well as the software framework for calibration, reconstruction and simulation. Finally the first experience from commissioning and operating the detector and the performance as an air shower detector will be discussed.

preprint2012arXiv

Lateral Distribution of Muons in IceCube Cosmic Ray Events

In cosmic ray air showers, the muon lateral separation from the center of the shower is a measure of the transverse momentum that the muon parent acquired in the cosmic ray interaction. IceCube has observed cosmic ray interactions that produce muons laterally separated by up to 400 m from the shower core, a factor of 6 larger distance than previous measurements. These muons originate in high pT (> 2 GeV/c) interactions from the incident cosmic ray, or high-energy secondary interactions. The separation distribution shows a transition to a power law at large values, indicating the presence of a hard pT component that can be described by perturbative quantum chromodynamics. However, the rates and the zenith angle distributions of these events are not well reproduced with the cosmic ray models tested here, even those that include charm interactions. This discrepancy may be explained by a larger fraction of kaons and charmed particles than is currently incorporated in the simulations.

preprint2012arXiv

Search for Galactic PeV Gamma Rays with the IceCube Neutrino Observatory

Gamma-ray induced air showers are notable for their lack of muons, compared to hadronic showers. Hence, air shower arrays with large underground muon detectors can select a sample greatly enriched in photon showers by rejecting showers containing muons. IceCube is sensitive to muons with energies above ~500 GeV at the surface, which provides an efficient veto system for hadronic air showers with energies above 1 PeV. One year of data from the 40-string IceCube configuration was used to perform a search for point sources and a Galactic diffuse signal. No sources were found, resulting in a 90% C.L. upper limit on the ratio of gamma rays to cosmic rays of 1.2 x 10^(-3)for the flux coming from the Galactic Plane region (-80 deg < l < -30 deg; -10 deg < b < 5 deg) in the energy range 1.2 - 6.0 PeV. In the same energy range, point source fluxes with E^(-2) spectra have been excluded at a level of (E/TeV)^2 dΦ/dE ~ 10^(-12)-10^(-11) cm^2/s/TeV depending on source declination. The complete IceCube detector will have a better sensitivity, due to the larger detector size, improved reconstruction and vetoing techniques. Preliminary data from the nearly-final IceCube detector configuration has been used to estimate the 5 year sensitivity of the full detector. It is found to be more than an order of magnitude better, allowing the search for PeV extensions of known TeV gamma-ray emitters.

preprint2012arXiv

Search for Relativistic Magnetic Monopoles with IceCube

We present the first results in the search for relativistic magnetic monopoles with the IceCube detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume of 1 km$^{3}$. This analysis searches data taken on the partially completed detector during 2007 when roughly 0.2 km$^{3}$ of ice was instrumented. The lack of candidate events leads to an upper limit on the flux of relativistic magnetic monopoles of $Φ_{\mathrm{90%C.L.}}\sim 3\e{-18}\fluxunits$ for $β\geq0.8$. This is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost $γ$ below $10^{7}$. This result is then interpreted for a wide range of mass and kinetic energy values.

preprint2012arXiv

Searches for periodic neutrino emission from binary systems with 22 and 40 strings of IceCube

In this paper we present the results of searches for periodic neutrino emission from a catalog of binary systems. Such modulation, observed in the photon flux, would be caused by the geometry of these systems. In the analysis, the period is fixed by these photon observations, while the phase and duration of the neutrino emission are treated as free parameters to be fit with the data. If the emission occurs during ~20% or less of the total period, this analysis achieves better sensitivity than a time-integrated analysis. We use the IceCube data taken from May 31, 2007 to April 5, 2008 with its 22-string configuration, and from April 5, 2008 to May 20, 2009 with its 40-string configuration. No evidence for neutrino emission is found, with the strongest excess occurring for Cygnus X-3 at 2.1 sigma significance after accounting for trials. Neutrino flux upper limits for both periodic and time-integrated emission are provided.

preprint2012arXiv

Use of event-level neutrino telescope data in global fits for theories of new physics

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

preprint2011arXiv

A Search for a Diffuse Flux of Astrophysical Muon Neutrinos with the IceCube 40-String Detector

The IceCube Neutrino Observatory is a 1 km$^{3}$ detector currently taking data at the South Pole. One of the main strategies used to look for astrophysical neutrinos with IceCube is the search for a diffuse flux of high-energy neutrinos from unresolved sources. A hard energy spectrum of neutrinos from isotropically distributed astrophysical sources could manifest itself as a detectable signal that may be differentiated from the atmospheric neutrino background by spectral measurement. This analysis uses data from the IceCube detector collected in its half completed configuration which operated between April 2008 and May 2009 to search for a diffuse flux of astrophysical muon neutrinos. A total of 12,877 upward going candidate neutrino events have been selected for this analysis. No evidence for a diffuse flux of astrophysical muon neutrinos was found in the data set leading to a 90 percent C.L. upper limit on the normalization of an $E^{-2}$ astrophysical $ν_μ$ flux of $8.9 \times 10^{-9} \ \mathrm{GeV \ cm^{-2} \ s^{-1} \ sr^{-1}}$. The analysis is sensitive in the energy range between $35 \ \mathrm{TeV} - 7 \ \mathrm{PeV}$. The 12,877 candidate neutrino events are consistent with atmospheric muon neutrinos measured from 332 GeV to 84 TeV and no evidence for a prompt component to the atmospheric neutrino spectrum is found.

preprint2011arXiv

Background studies for acoustic neutrino detection at the South Pole

The detection of acoustic signals from ultra-high energy neutrino interactions is a promising method to measure the tiny flux of cosmogenic neutrinos expected on Earth. The energy threshold for this process depends strongly on the absolute noise level in the target material. The South Pole Acoustic Test Setup (SPATS), deployed in the upper part of four boreholes of the IceCube Neutrino Observatory, has monitored the noise in Antarctic ice at the geographic South Pole for more than two years down to 500 m depth. The noise is very stable and Gaussian distributed. Lacking an in-situ calibration up to now, laboratory measurements have been used to estimate the absolute noise level in the 10 to 50 kHz frequency range to be smaller than 20 mPa. Using a threshold trigger, sensors of the South Pole Acoustic Test Setup registered acoustic pulse-like events in the IceCube detector volume and its vicinity. Acoustic signals from refreezing IceCube holes and from anthropogenic sources have been used to localize acoustic events. Monte Carlo simulations of sound propagating from the established sources to the SPATS sensors have allowed to check corresponding model expectations. An upper limit on the neutrino flux at energies $E_ν> 10^{11}$ GeV is derived from acoustic data taken over eight months.

preprint2011arXiv

Constraints on high-energy neutrino emission from SN 2008D

SN 2008D, a core collapse supernova at a distance of 27 Mpc, was serendipitously discovered by the Swift satellite through an associated X-ray flash. Core collapse supernovae have been observed in association with long gamma-ray bursts and X-ray flashes and a physical connection is widely assumed. This connection could imply that some core collapse supernovae possess mildly relativistic jets in which high-energy neutrinos are produced through proton-proton collisions. The predicted neutrino spectra would be detectable by Cherenkov neutrino detectors like IceCube. A search for a neutrino signal in temporal and spatial correlation with the observed X-ray flash of SN 2008D was conducted using data taken in 2007-2008 with 22 strings of the IceCube detector. Events were selected based on a boosted decision tree classifier trained with simulated signal and experimental background data. The classifier was optimized to the position and a "soft jet" neutrino spectrum assumed for SN 2008D. Using three search windows placed around the X-ray peak, emission time scales from 100 - 10000 s were probed. No events passing the cuts were observed in agreement with the signal expectation of 0.13 events. Upper limits on the muon neutrino flux from core collapse supernovae were derived for different emission time scales and the principal model parameters were constrained.

preprint2011arXiv

Constraints on the Extremely-high Energy Cosmic Neutrino Flux with the IceCube 2008-2009 Data

We report on a search for extremely-high energy neutrinos with energies greater than $10^6$ GeV using the data taken with the IceCube detector at the South Pole. The data was collected between April 2008 and May 2009 with the half completed IceCube array. The absence of signal candidate events in the sample of 333.5 days of livetime significantly improves model independent limit from previous searches and allows to place a limit on the diffuse flux of cosmic neutrinos with an $E^{-2}$ spectrum in the energy range $2.0 \times 10^{6}$ $-$ $6.3 \times 10^{9}$ GeV to a level of $E^2 ϕ\leq 3.6 \times 10^{-8}$ ${\rm GeV cm^{-2} sec^{-1}sr^{-1}}$.

preprint2011arXiv

First search for atmospheric and extraterrestrial neutrino-induced cascades with the IceCube detector

We report on the first search for atmospheric and for diffuse astrophysical neutrino-induced showers (cascades) in the IceCube detector using 257 days of data collected in the year 2007-2008 with 22 strings active. A total of 14 events with energies above 16 TeV remained after event selections in the diffuse analysis, with an expected total background contribution of $8.3\pm 3.6$. At 90% confidence we set an upper limit of $E^2Φ_{90%CL}<3.6\times10^{-7} GeV \cdot cm^{-2} \cdot s^{-1}\cdot sr^{-1} $ on the diffuse flux of neutrinos of all flavors in the energy range between 24 TeV and 6.6 PeV assuming that $Φ\propto E^{-2}$ and that the flavor composition of the $ν_e : ν_μ: ν_τ$ flux is $1 : 1 : 1$ at the Earth. The atmospheric neutrino analysis was optimized for lower energies. A total of 12 events were observed with energies above 5 TeV. The observed number of events is consistent with the expected background, within the uncertainties.

preprint2011arXiv

Limits on Neutrino Emission from Gamma-Ray Bursts with the 40 String IceCube Detector

IceCube has become the first neutrino telescope with a sensitivity below the TeV neutrino flux predicted from gamma-ray bursts if GRBs are responsible for the observed cosmic-ray flux above $10^{18}$ eV. Two separate analyses using the half-complete IceCube detector, one a dedicated search for neutrinos from $p γ$-interactions in the prompt phase of the GRB fireball, and the other a generic search for any neutrino emission from these sources over a wide range of energies and emission times, produced no evidence for neutrino emission, excluding prevailing models at 90% confidence.

preprint2011arXiv

Multi-year search for dark matter annihilations in the Sun with the AMANDA-II and IceCube detectors

A search for an excess of muon-neutrinos from dark matter annihilations in the Sun has been performed with the AMANDA-II neutrino telescope using data collected in 812 days of livetime between 2001 and 2006 and 149 days of livetime collected with the AMANDA-II and the 40-string configuration of IceCube during 2008 and early 2009. No excess over the expected atmospheric neutrino background has been observed. We combine these results with the previously published IceCube limits obtained with data taken during 2007 to obtain a total livetime of 1065 days. We provide an upper limit at 90% confidence level on the annihilation rate of captured neutralinos in the Sun, as well as the corresponding muon flux limit at the Earth, both as functions of the neutralino mass in the range 50 GeV-5000 GeV. We also derive a limit on the neutralino-proton spin-dependent and spin-independent cross section. The limits presented here improve the previous results obtained by the collaboration between a factor of two and five, as well as extending the neutralino masses probed down to 50 GeV. The spin-dependent cross section limits are the most stringent so far for neutralino masses above 200 GeV, and well below direct search results in the mass range from 50 GeV to 5 TeV.

preprint2011arXiv

Neutrino analysis of the September 2010 Crab Nebula flare and time-integrated constraints on neutrino emission from the Crab using IceCube

We present the results for a search of high-energy muon neutrinos with the IceCube detector in coincidence with the Crab nebula flare reported on September 2010 by various experiments. Due to the unusual flaring state of the otherwise steady source we performed a prompt analysis of the 79-string configuration data to search for neutrinos that might be emitted along with the observed gamma-rays. We performed two different and complementary data selections of neutrino events in the time window of 10 days around the flare. One event selection is optimized for discovery of E^-2 neutrino spectrum typical of 1st order Fermi acceleration. A similar event selection has also been applied to the 40-string data to derive the time-integrated limits to the neutrino emission from the Crab. The other event selection was optimized for discovery of neutrino spectra with softer spectral index and TeV energy cut-offs as observed for various galactic sources in gamma-rays. The 90% CL best upper limits on the Crab flux during the 10 day flare are 4.73 x 10^-11 cm-2 s-1 TeV-1 for an E^-2 neutrino spectrum and 2.50 x 10^-10 cm-2 s-1 TeV-1 for a softer neutrino spectra of E-2.7, as indicated by Fermi measurements during the flare. IceCube has also set a time-integrated limit on the neutrino emission of the Crab using 375.5 days of livetime of the 40-string configuration data. This limit is compared to existing models of neutrino production from the Crab and its impact on astrophysical parameters is discussed. The most optimistic predictions of some models are already rejected by the IceCube neutrino telescope with more than 90% CL.

preprint2011arXiv

Observation of an Anisotropy in the Galactic Cosmic Ray arrival direction at 400 TeV with IceCube

In this paper we report the first observation in the Southern hemisphere of an energy dependence in the Galactic cosmic ray anisotropy up to a few hundred TeV. This measurement was performed using cosmic ray induced muons recorded by the partially deployed IceCube observatory between May 2009 and May 2010. The data include a total of 33$\times 10^{9}$ muon events with a median angular resolution of $\sim3^{\circ}$ degrees. A sky map of the relative intensity in arrival direction over the Southern celestial sky is presented for cosmic ray median energies of 20 and 400 TeV. The same large-scale anisotropy observed at median energies around 20 TeV is not present at 400 TeV. Instead, the high energy skymap shows a different anisotropy structure including a deficit with a post-trial significance of -6.3$σ$. This anisotropy reveals a new feature of the Galactic cosmic ray distribution, which must be incorporated into theories of the origin and propagation of cosmic rays.

preprint2011arXiv

Observation of Anisotropy in the Arrival Directions of Galactic Cosmic Rays at Multiple Angular Scales with IceCube

Between May 2009 and May 2010, the IceCube neutrino detector at the South Pole recorded 32 billion muons generated in air showers produced by cosmic rays with a median energy of 20 TeV. With a data set of this size, it is possible to probe the southern sky for per-mille anisotropy on all angular scales in the arrival direction distribution of cosmic rays. Applying a power spectrum analysis to the relative intensity map of the cosmic ray flux in the southern hemisphere, we show that the arrival direction distribution is not isotropic, but shows significant structure on several angular scales. In addition to previously reported large-scale structure in the form of a strong dipole and quadrupole, the data show small-scale structure on scales between 15 degrees and 30 degrees. The skymap exhibits several localized regions of significant excess and deficit in cosmic ray intensity. The relative intensity of the smaller-scale structures is about a factor of 5 weaker than that of the dipole and quadrupole structure. The most significant structure, an excess localized at right ascension 122.4 degrees and declination -47.4 degrees, extends over at least 20 degrees in right ascension and has a post-trials significance of 5.3 sigma. The origin of this anisotropy is still unknown.

preprint2011arXiv

Search for Dark Matter from the Galactic Halo with the IceCube Neutrino Observatory

Self-annihilating or decaying dark matter in the Galactic halo might produce high energy neutrinos detectable with neutrino telescopes. We have conducted a search for such a signal using 276 days of data from the IceCube 22-string configuration detector acquired during 2007 and 2008. The effect of halo model choice in the extracted limit is reduced by performing a search that considers the outer halo region and not the Galactic Center. We constrain any large scale neutrino anisotropy and are able to set a limit on the dark matter self-annihilation cross section of <σ_{A}v> \simeq 10^{-22} cm^3/s for WIMP masses above 1 TeV, assuming a monochromatic neutrino line spectrum.

preprint2011arXiv

Searching for soft relativistic jets in Core-collapse Supernovae with the IceCube Optical Follow-up Program

Context. Transient neutrino sources such as Gamma-Ray Bursts (GRBs) and Supernovae (SNe) are hypothesized to emit bursts of high-energy neutrinos on a time-scale of \lesssim 100 s. While GRB neutrinos would be produced in high relativistic jets, core-collapse SNe might host soft-relativistic jets, which become stalled in the outer layers of the progenitor star leading to an efficient production of high-energy neutrinos. Aims. To increase the sensitivity to these neutrinos and identify their sources, a low-threshold optical follow-up program for neutrino multiplets detected with the IceCube observatory has been implemented. Methods. If a neutrino multiplet, i.e. two or more neutrinos from the same direction within 100 s, is found by IceCube a trigger is sent to the Robotic Optical Transient Search Experiment, ROTSE. The 4 ROTSE telescopes immediately start an observation program of the corresponding region of the sky in order to detect an optical counterpart to the neutrino events. Results. No statistically significant excess in the rate of neutrino multiplets has been observed and furthermore no coincidence with an optical counterpart was found. Conclusion. The search allows, for the first time, to set stringent limits on current models predicting a high-energy neutrino flux from soft relativistic hadronic jets in core-collapse SNe. We conclude that a sub-population of SNe with typical Lorentz boost factor and jet energy of 10 and 3\times10^{51} erg, respectively, does not exceed 4.2% at 90% confidence.

preprint2011arXiv

Time-Dependent Searches for Point Sources of Neutrinos with the 40-String and 22-String Configurations of IceCube

This paper presents searches for flaring sources of neutrinos using the IceCube neutrino telescope. For the first time, a search is performed over the entire parameter space of energy, direction and time looking for neutrino flares of 20 microseconds to a year duration from astrophysical sources among the atmospheric neutrino and muon backgrounds. Searches which integrate over time are less sensitive to flares because they are affected by a larger background of atmospheric neutrinos and muons that can be reduced by the time constraint. Flaring sources considered here, such as active galactic nuclei, soft gamma-ray repeaters and gamma-ray bursts, are promising candidate neutrino emitters. We used mainly data taken between April 5, 2008 and May 20, 2009 by a partially completed configuration of IceCube with 40 strings. For the presented searches an unbinned maximum likelihood method is used with a time-dependent prior to test several different source hypotheses. An "untriggered" search covers any possible time-dependent emission from sources not correlated to any other observation using other astrophysical messengers such as photons. Moreover, a similar time scan is performed for a predefined catalogue of sources that exhibit intense photon flares. Searches triggered by multi-wavelength information on flares from blazars and soft gamma-ray repeaters are performed using the 40 string data and also the data taken by the previous configuration of 22 strings in operation between May 31, 2007 and April 5, 2008. Flares for which extensive and continuous monitoring is available from Fermi-LAT and SWIFT and flares detected by imaging Cherenkov telescopes with shorter time-scale monitoring are considered. The results from all searches are compatible with a fluctuation of the background.

preprint2010arXiv

A new limit on the Ultra-High-Energy Cosmic-Ray flux with the Westerbork Synthesis Radio Telescope

A particle cascade (shower) in a dielectric, for example as initiated by an ultra-high energy cosmic ray, will have an excess of electrons which will emit coherent Čerenkov radiation, known as the Askaryan effect. In this work we study the case in which such a particle shower occurs in a medium just below its surface. We show, for the first time, that the radiation transmitted through the surface is independent of the depth of the shower below the surface when observed from far away, apart from trivial absorption effects. As a direct application we use the recent results of the NuMoon project, where a limit on the neutrino flux for energies above $10^{22}$\,eV was set using the Westerbork Synthesis Radio Telescope by measuring pulsed radio emission from the Moon, to set a limit on the flux of ultra-high-energy cosmic rays.

preprint2010arXiv

Calibration and Characterization of the IceCube Photomultiplier Tube

Over 5,000 PMTs are being deployed at the South Pole to compose the IceCube neutrino observatory. Many are placed deep in the ice to detect Cherenkov light emitted by the products of high-energy neutrino interactions, and others are frozen into tanks on the surface to detect particles from atmospheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 made by Hamamatsu Photonics. This paper describes the laboratory characterization and calibration of these PMTs before deployment. PMTs were illuminated with pulses ranging from single photons to saturation level. Parameterizations are given for the single photoelectron charge spectrum and the saturation behavior. Time resolution, late pulses and afterpulses are characterized. Because the PMTs are relatively large, the cathode sensitivity uniformity was measured. The absolute photon detection efficiency was calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured characteristics are discussed in the context of their relevance to IceCube event reconstruction and simulation efforts.

preprint2010arXiv

Constraints on the flux of Ultra-High Energy neutrinos from WSRT observations

Ultra-high energy (UHE) neutrinos and cosmic rays initiate particle cascades underneath the Moon's surface. These cascades have a negative charge excess and radiate Cherenkov radio emission in a process known as the Askaryan effect. The optimal frequency window for observation of these pulses with radio telescopes on the Earth is around 150 MHz. By observing the Moon with the Westerbork Synthesis Radio Telescope array we are able to set a new limit on the UHE neutrino flux. The PuMa II backend is used to monitor the Moon in 4 frequency bands between 113 and 175 MHz with a sampling frequency of 40 MHz. The narrowband radio interference is digitally filtered out and the dispersive effect of the Earth's ionosphere is compensated for. A trigger system is implemented to search for short pulses. By inserting simulated pulses in the raw data, the detection efficiency for pulses of various strength is calculated. With 47.6 hours of observation time, we are able to set a limit on the UHE neutrino flux. This new limit is an order of magnitude lower than existing limits. In the near future, the digital radio array LOFAR will be used to achieve an even lower limit.

preprint2010arXiv

Measurement of the atmospheric neutrino energy spectrum from 100 GeV to 400 TeV with IceCube

A measurement of the atmospheric muon neutrino energy spectrum from 100 GeV to 400 TeV was performed using a data sample of about 18,000 up-going atmospheric muon neutrino events in IceCube. Boosted decision trees were used for event selection to reject mis-reconstructed atmospheric muons and obtain a sample of up-going muon neutrino events. Background contamination in the final event sample is less than one percent. This is the first measurement of atmospheric neutrinos up to 400 TeV, and is fundamental to understanding the impact of this neutrino background on astrophysical neutrino observations with IceCube. The measured spectrum is consistent with predictions for the atmospheric muon neutrino plus muon antineutrino flux.

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

Search for a Lorentz-violating sidereal signal with atmospheric neutrinos in IceCube

A search for sidereal modulation in the flux of atmospheric muon neutrinos in IceCube was performed. Such a signal could be an indication of Lorentz-violating physics. Neutrino oscillation models, derivable from extensions to the Standard Model, allow for neutrino oscillations that depend on the neutrino's direction of propagation. No such direction-dependent variation was found. A discrete Fourier transform method was used to constrain the Lorentz and CPT-violating coefficients in one of these models. Due to the unique high energy reach of IceCube, it was possible to improve constraints on certain Lorentz-violating oscillations by three orders of magnitude with respect to limits set by other experiments.

preprint2010arXiv

Simulation of radio emission from air showers in atmospheric electric fields

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

preprint2010arXiv

The first search for extremely-high energy cosmogenic neutrinos with the IceCube Neutrino Observatory

We report on the results of the search for extremely-high energy (EHE) neutrinos with energies above $10^7$ GeV obtained with the partially ($\sim$30%) constructed IceCube in 2007. From the absence of signal events in the sample of 242.1 days of effective livetime, we derive a 90% C.L. model independent differential upper limit based on the number of signal events per energy decade at $E^2 ϕ_{ν_e+ν_μ+ν_τ}\simeq 1.4 \times 10^{-6}$ GeV cm$^{-2}$ sec$^{-1}$ sr$^{-1}$ for neutrinos in the energy range from $3\times10^7$ to $3\times10^9$ GeV.

preprint2009arXiv

Air Shower Measurements with LOFAR

Air showers from cosmic rays emit short, intense radio pulses. LOFAR is a new radio telescope, that is being built in the Netherlands and Europe. Designed primarily as a radio interferometer, the core of LOFAR will have a high density of radio antennas, which will be extremely well calibrated. This makes LOFAR a unique tool for the study of the radio properties of single air showers. Triggering on the radio emission from air showers means detecting a short radio pulse and discriminating real events from radio interference. At LOFAR we plan to search for pulses in the digital data stream - either from single antennas or from already beam-formed data - and calculate several parameters characterizing the pulse shape to pick out real events in a second stage. In addition, we will have a small scintillator array to test and confirm the performance of the radio only trigger.

preprint2009arXiv

Improved flux limits for neutrinos with energies above 10$^{22}$ eV from observations with the Westerbork Synthesis Radio Telescope

Particle cascades initiated by ultra-high energy (UHE) neutrinos in the lunar regolith will emit an electromagnetic pulse with a time duration of the order of nano seconds through a process known as the Askaryan effect. It has been shown that in an observing window around 150 MHz there is a maximum chance for detecting this radiation with radio telescopes commonly used in astronomy. In 50 hours of observation time with the Westerbork Synthesis Radio Telescope array we have set a new limit on the flux of neutrinos, summed over all flavors, with energies in excess of $4\times10^{22}$ eV.

preprint2009arXiv

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

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

preprint2009arXiv

Monte Carlo simulations of air showers in atmospheric electric fields

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