Source author record

G. Golup

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

ResearcherUnclaimed source record
Topic enrichment pending

Catalog footprint

What is connected

97works
0topics
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

97 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

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

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2021arXiv

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

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

An All-Sky Search for Three Flavors of Neutrinos from Gamma-Ray Bursts with the IceCube Neutrino Observatory

We present the results and methodology of a search for neutrinos produced in the decay of charged pions created in interactions between protons and gamma-rays during the prompt emission of 807 gamma-ray bursts (GRBs) over the entire sky. This three-year search is the first in IceCube for shower-like Cherenkov light patterns from electron, muon, and tau neutrinos correlated with GRBs. We detect five low-significance events correlated with five GRBs. These events are consistent with the background expectation from atmospheric muons and neutrinos. The results of this search in combination with those of IceCube's four years of searches for track-like Cherenkov light patterns from muon neutrinos correlated with Northern-Hemisphere GRBs produce limits that tightly constrain current models of neutrino and ultra high energy cosmic ray production in GRB fireballs.

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

All-flavour Search for Neutrinos from Dark Matter Annihilations in the Milky Way with IceCube/DeepCore

We present the first IceCube search for a signal of dark matter annihilations in the Milky Way using all-flavour neutrino-induced particle cascades. The analysis focuses on the DeepCore sub-detector of IceCube, and uses the surrounding IceCube strings as a veto region in order to select starting events in the DeepCore volume. We use 329 live-days of data from IceCube operating in its 86-string configuration during 2011-2012. No neutrino excess is found, the final result being compatible with the background-only hypothesis. From this null result, we derive upper limits on the velocity-averaged self-annihilation cross-section, < σ_A v >, for dark matter candidate masses ranging from 30 GeV up to 10 TeV, assuming both a cuspy and a flat-cored dark matter halo profile. For dark matter masses between 200 GeV and 10 TeV, the results improve on all previous IceCube results on < σ_A v >, reaching a level of 10^{-23} cm^3 s^-1, depending on the annihilation channel assumed, for a cusped NFW profile. The analysis demonstrates that all-flavour searches are competitive with muon channel searches despite the intrinsically worse angular resolution of cascades compared to muon tracks in IceCube.

preprint2016arXiv

Anisotropy in Cosmic-Ray Arrival Directions in the Southern Hemisphere with Six Years of Data from the IceCube Detector

The IceCube Neutrino Observatory has accumulated a total of 318 billion cosmic-ray induced muon events between May 2009 and May 2015. This data set was used for a detailed analysis of the cosmic-ray arrival direction anisotropy in the TeV to PeV energy range. The observed global anisotropy features large regions of relative excess and deficit, with amplitudes on the order of $10^{-3}$ up to about 100 TeV. A decomposition of the arrival direction distribution into spherical harmonics shows that most of the power is contained in the low-multipole ($\ell\leq 4$) moments. However, higher multipole components are found to be statistically significant down to an angular scale of less than $10^{\circ}$, approaching the angular resolution of the detector. Above 100 TeV, a change in the morphology of the arrival direction distribution is observed, and the anisotropy is characterized by a wide relative deficit whose amplitude increases with primary energy up to at least 5\,PeV, the highest energies currently accessible to IceCube. No time dependence of the large- and small-scale structures is observed in the six-year period covered by this analysis. The high-statistics data set reveals more details on the properties of the anisotropy and is potentially able to shed light on the various physical processes that are responsible for the complex angular structure and energy evolution.

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

Characterization of the Atmospheric Muon Flux in IceCube

Muons produced in atmospheric cosmic ray showers account for the by far dominant part of the event yield in large-volume underground particle detectors. The IceCube detector, with an instrumented volume of about a cubic kilometer, has the potential to conduct unique investigations on atmospheric muons by exploiting the large collection area and the possibility to track particles over a long distance. Through detailed reconstruction of energy deposition along the tracks, the characteristics of muon bundles can be quantified, and individual particles of exceptionally high energy identified. The data can then be used to constrain the cosmic ray primary flux and the contribution to atmospheric lepton fluxes from prompt decays of short-lived hadrons. In this paper, techniques for the extraction of physical measurements from atmospheric muon events are described and first results are presented. The multiplicity spectrum of TeV muons in cosmic ray air showers for primaries in the energy range from the knee to the ankle is derived and found to be consistent with recent results from surface detectors. The single muon energy spectrum is determined up to PeV energies and shows a clear indication for the emergence of a distinct spectral component from prompt decays of short-lived hadrons. The magnitude of the prompt flux, which should include a substantial contribution from light vector meson di-muon decays, is consistent with current theoretical predictions.

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

Improved limits on dark matter annihilation in the Sun with the 79-string IceCube detector and implications for supersymmetry

We present an improved event-level likelihood formalism for including neutrino telescope data in global fits to new physics. We derive limits on spin-dependent dark matter-proton scattering by employing the new formalism in a re-analysis of data from the 79-string IceCube search for dark matter annihilation in the Sun, including explicit energy information for each event. The new analysis excludes a number of models in the weak-scale minimal supersymmetric standard model (MSSM) for the first time. This work is accompanied by the public release of the 79-string IceCube data, as well as an associated computer code for applying the new likelihood to arbitrary dark matter models.

preprint2016arXiv

Lowering IceCube's Energy Threshold for Point Source Searches in the Southern Sky

Observation of a point source of astrophysical neutrinos would be a "smoking gun" signature of a cosmic-ray accelerator. While IceCube has recently discovered a diffuse flux of astrophysical neutrinos, no localized point source has been observed. Previous IceCube searches for point sources in the southern sky were restricted by either an energy threshold above a few hundred TeV or poor neutrino angular resolution. Here we present a search for southern sky point sources with greatly improved sensitivities to neutrinos with energies below 100 TeV. By selecting charged-current $ν_μ$ interacting inside the detector, we reduce the atmospheric background while retaining efficiency for astrophysical neutrino-induced events reconstructed with sub-degree angular resolution. The new event sample covers three years of detector data and leads to a factor of ten improvement in sensitivity to point sources emitting below 100 TeV in the southern sky. No statistically significant evidence of point sources was found, and upper limits are set on neutrino emission from individual sources. A posteriori analysis of the highest-energy ~100 TeV starting event in the sample found that this event alone represents a $2.8σ$ deviation from the hypothesis that the data consists only of atmospheric background.

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2016arXiv

Search for Astrophysical Tau Neutrinos in Three Years of IceCube Data

The IceCube Neutrino Observatory has observed a diffuse flux of TeV-PeV astrophysical neutrinos at 5.7σ significance from an all-flavor search. The direct detection of tau neutrinos in this flux has yet to occur. Tau neutrinos become distinguishable from other flavors in IceCube at energies above a few hundred TeV, when the cascade from the tau neutrino charged current interaction becomes resolvable from the cascade from the tau lepton decay. This paper presents results from a dedicated search for tau neutrinos with energies between 214 TeV and 72 PeV. The analysis searches for IceCube optical sensors that observe two separate pulses in a single event - one from the tau neutrino interaction, and a second from the tau decay. This is the first IceCube tau neutrino search to be more sensitive to tau neutrinos than to any other neutrino flavor. No candidate events were observed in three years of IceCube data. For the first time, a differential upper limit on astrophysical tau neutrinos is derived around the PeV energy region, which is nearly three orders of magnitude lower in energy than previous limits from dedicated tau neutrino searches.

preprint2016arXiv

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

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

preprint2016arXiv

Search for ultrarelativistic magnetic monopoles with the Pierre Auger Observatory

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

preprint2016arXiv

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

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

preprint2016arXiv

The Pierre Auger Observatory Upgrade - Preliminary Design Report

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

preprint2015arXiv

A combined maximum-likelihood analysis of the high-energy astrophysical neutrino flux measured with IceCube

Evidence for an extraterrestrial flux of high-energy neutrinos has now been found in multiple searches with the IceCube detector. The first solid evidence was provided by a search for neutrino events with deposited energies $\gtrsim30$ TeV and interaction vertices inside the instrumented volume. Recent analyses suggest that the extraterrestrial flux extends to lower energies and is also visible with throughgoing, $ν_μ$-induced tracks from the Northern hemisphere. Here, we combine the results from six different IceCube searches for astrophysical neutrinos in a maximum-likelihood analysis. The combined event sample features high-statistics samples of shower-like and track-like events. The data are fit in up to three observables: energy, zenith angle and event topology. Assuming the astrophysical neutrino flux to be isotropic and to consist of equal flavors at Earth, the all-flavor spectrum with neutrino energies between 25 TeV and 2.8 PeV is well described by an unbroken power law with best-fit spectral index $-2.50\pm0.09$ and a flux at 100 TeV of $\left(6.7_{-1.2}^{+1.1}\right)\cdot10^{-18}\,\mathrm{GeV}^{-1}\mathrm{s}^{-1}\mathrm{sr}^{-1}\mathrm{cm}^{-2}$. Under the same assumptions, an unbroken power law with index $-2$ is disfavored with a significance of 3.8 $σ$ ($p=0.0066\%$) with respect to the best fit. This significance is reduced to 2.1 $σ$ ($p=1.7\%$) if instead we compare the best fit to a spectrum with index $-2$ that has an exponential cut-off at high energies. Allowing the electron neutrino flux to deviate from the other two flavors, we find a $ν_e$ fraction of $0.18\pm0.11$ at Earth. The sole production of electron neutrinos, which would be characteristic of neutron-decay dominated sources, is rejected with a significance of 3.6 $σ$ ($p=0.014\%$).

preprint2015arXiv

Atmospheric and Astrophysical Neutrinos above 1 TeV Interacting in IceCube

The IceCube Neutrino Observatory was designed primarily to search for high-energy (TeV--PeV) neutrinos produced in distant astrophysical objects. A search for $\gtrsim 100$~TeV neutrinos interacting inside the instrumented volume has recently provided evidence for an isotropic flux of such neutrinos. At lower energies, IceCube collects large numbers of neutrinos from the weak decays of mesons in cosmic-ray air showers. Here we present the results of a search for neutrino interactions inside IceCube's instrumented volume between 1~TeV and 1~PeV in 641 days of data taken from 2010--2012, lowering the energy threshold for neutrinos from the southern sky below 10 TeV for the first time, far below the threshold of the previous high-energy analysis. Astrophysical neutrinos remain the dominant component in the southern sky down to 10 TeV. From these data we derive new constraints on the diffuse astrophysical neutrino spectrum, $Φ_ν = 2.06^{+0.4}_{-0.3} \times 10^{-18} \left({E_ν}/{10^5 \,\, \rm{GeV}} \right)^{-2.46 \pm 0.12} {\rm {GeV^{-1} \, cm^{-2} \, sr^{-1} \, s^{-1}} } $, as well as the strongest upper limit yet on the flux of neutrinos from charmed-meson decay in the atmosphere, 1.52 times the benchmark theoretical prediction used in previous IceCube results at 90\% confidence.

preprint2015arXiv

Determining neutrino oscillation parameters from atmospheric muon neutrino disappearance with three years of IceCube DeepCore data

We present a measurement of neutrino oscillations via atmospheric muon neutrino disappearance with three years of data of the completed IceCube neutrino detector. DeepCore, a region of denser instrumentation, enables the detection and reconstruction of atmospheric muon neutrinos between 10 GeV and 100 GeV, where a strong disappearance signal is expected. The detector volume surrounding DeepCore is used as a veto region to suppress the atmospheric muon background. Neutrino events are selected where the detected Cherenkov photons of the secondary particles minimally scatter, and the neutrino energy and arrival direction are reconstructed. Both variables are used to obtain the neutrino oscillation parameters from the data, with the best fit given by $Δm^2_{32}=2.72^{+0.19}_{-0.20}\times 10^{-3}\,\mathrm{eV}^2$ and $\sin^2θ_{23} = 0.53^{+0.09}_{-0.12}$ (normal mass hierarchy assumed). The results are compatible and comparable in precision to those of dedicated oscillation experiments.

preprint2015arXiv

Development of a General Analysis and Unfolding Scheme and its Application to Measure the Energy Spectrum of Atmospheric Neutrinos with IceCube

We present the development and application of a generic analysis scheme for the measurement of neutrino spectra with the IceCube detector. This scheme is based on regularized unfolding, preceded by an event selection which uses a Minimum Redundancy Maximum Relevance algorithm to select the relevant variables and a Random Forest for the classification of events. The analysis has been developed using IceCube data from the 59-string configuration of the detector. 27,771 neutrino candidates were detected in 346 days of livetime. A rejection of 99.9999% of the atmospheric muon background is achieved. The energy spectrum of the atmospheric neutrino flux is obtained using the TRUEE unfolding program. The unfolded spectrum of atmospheric muon neutrinos covers an energy range from 100 GeV to 1 PeV. Compared to the previous measurement using the detector in the 40-string configuration, the analysis presented here, extends the upper end of the atmospheric neutrino spectrum by more than a factor of two, reaching an energy region that has not been previously accessed by spectral measurements.

preprint2015arXiv

Evidence for Astrophysical Muon Neutrinos from the Northern Sky with IceCube

Results from the IceCube Neutrino Observatory have recently provided compelling evidence for the existence of a high energy astrophysical neutrino flux utilizing a dominantly Southern Hemisphere dataset consisting primarily of nu_e and nu_tau charged current and neutral current (cascade) neutrino interactions. In the analysis presented here, a data sample of approximately 35,000 muon neutrinos from the Northern sky was extracted from data taken during 659.5 days of livetime recorded between May 2010 and May 2012. While this sample is composed primarily of neutrinos produced by cosmic ray interactions in the Earth's atmosphere, the highest energy events are inconsistent with a hypothesis of solely terrestrial origin at 3.7 sigma significance. These neutrinos can, however, be explained by an astrophysical flux per neutrino flavor at a level of Phi(E_nu) = 9.9^{+3.9}_{-3.4} times 10^{-19} GeV^{-1} cm^{-2} sr^{-1} s^{-1} ({E_nu / 100 TeV})^{-2}, consistent with IceCube's Southern Hemisphere dominated result. Additionally, a fit for an astrophysical flux with an arbitrary spectral index was performed. We find a spectral index of 2.2^{+0.2}_{-0.2}, which is also in good agreement with the Southern Hemisphere result.

preprint2015arXiv

First combined search for neutrino point-sources in the Southern Hemisphere with the ANTARES and IceCube neutrino telescopes

We present the results of searches for point-like sources of neutrinos based on the first combined analysis of data from both the ANTARES and IceCube neutrino telescopes. The combination of both detectors which differ in size and location forms a window in the Southern sky where the sensitivity to point sources improves by up to a factor of two compared to individual analyses. Using data recorded by ANTARES from 2007 to 2012, and by IceCube from 2008 to 2011, we search for sources of neutrino emission both across the Southern sky and from a pre-selected list of candidate objects. No significant excess over background has been found in these searches, and flux upper limits for the candidate sources are presented for $E^{-2.5}$ and $E^{-2}$ power-law spectra with different energy cut-offs.

preprint2015arXiv

Flavor Ratio of Astrophysical Neutrinos above 35 TeV in IceCube

A diffuse flux of astrophysical neutrinos above $100\,\mathrm{TeV}$ has been observed at the IceCube Neutrino Observatory. Here we extend this analysis to probe the astrophysical flux down to $35\,\mathrm{TeV}$ and analyze its flavor composition by classifying events as showers or tracks. Taking advantage of lower atmospheric backgrounds for shower-like events, we obtain a shower-biased sample containing 129 showers and 8 tracks collected in three years from 2010 to 2013. We demonstrate consistency with the $(f_e:f_μ:f_τ)_\oplus\approx(1:1:1)_\oplus$ flavor ratio at Earth commonly expected from the averaged oscillations of neutrinos produced by pion decay in distant astrophysical sources. Limits are placed on non-standard flavor compositions that cannot be produced by averaged neutrino oscillations but could arise in exotic physics scenarios. A maximally track-like composition of $(0:1:0)_\oplus$ is excluded at $3.3σ$, and a purely shower-like composition of $(1:0:0)_\oplus$ is excluded at $2.3σ$.

preprint2015arXiv

Measurement of the Atmospheric $ν_e$ Spectrum with IceCube

We present a measurement of the atmospheric $ν_e$ spectrum at energies between 0.1 TeV and 100 TeV using data from the first year of the complete IceCube detector. Atmospheric $ν_e$ originate mainly from the decays of kaons produced in cosmic-ray air showers. This analysis selects 1078 fully contained events in 332 days of livetime, then identifies those consistent with particle showers. A likelihood analysis with improved event selection extends our previous measurement of the conventional $ν_e$ fluxes to higher energies. The data constrain the conventional $ν_e$ flux to be $1.3^{+0.4}_{-0.3}$ times a baseline prediction from a Honda's calculation, including the knee of the cosmic-ray spectrum. A fit to the kaon contribution ($ξ$) to the neutrino flux finds a kaon component that is $ξ=1.3^{+0.5}_{-0.4}$ times the baseline value. The fitted/measured prompt neutrino flux from charmed hadron decays strongly depends on the assumed astrophysical flux and shape. If the astrophysical component follows a power law, the result for the prompt flux is $0.0^{+3.0}_{-0.0}$ times a calculated flux based on the work by Enberg, Reno and Sarcevic.

preprint2015arXiv

Search for Dark Matter Annihilation in the Galactic Center with IceCube-79

The Milky Way is expected to be embedded in a halo of dark matter particles, with the highest density in the central region, and decreasing density with the halo-centric radius. Dark matter might be indirectly detectable at Earth through a flux of stable particles generated in dark matter annihilations and peaked in the direction of the Galactic Center. We present a search for an excess flux of muon (anti-) neutrinos from dark matter annihilation in the Galactic Center using the cubic-kilometer-sized IceCube neutrino detector at the South Pole. There, the Galactic Center is always seen above the horizon. Thus, new and dedicated veto techniques against atmospheric muons are required to make the southern hemisphere accessible for IceCube. We used 319.7 live-days of data from IceCube operating in its 79-string configuration during 2010 and 2011. No neutrino excess was found and the final result is compatible with the background. We present upper limits on the self-annihilation cross-section, $\left<σ_\mathrm{A} v\right>$, for WIMP masses ranging from 30 GeV up to 10 TeV, assuming cuspy (NFW) and flat-cored (Burkert) dark matter halo profiles, reaching down to $\simeq 4 \cdot 10^{-24}$ cm$^3$ s$^{-1}$, and $\simeq 2.6 \cdot 10^{-23}$ cm$^3$ s$^{-1}$ for the $ν\overlineν$ channel, respectively.

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

Search for Prompt Neutrino Emission from Gamma-Ray Bursts with IceCube

We present constraints derived from a search of four years of IceCube data for a prompt neutrino flux from gamma-ray bursts (GRBs). A single low-significance neutrino, compatible with the atmospheric neutrino background, was found in coincidence with one of the 506 observed bursts. Although GRBs have been proposed as candidate sources for ultra-high energy cosmic rays, our limits on the neutrino flux disfavor much of the parameter space for the latest models. We also find that no more than $\sim1\%$ of the recently observed astrophysical neutrino flux consists of prompt emission from GRBs that are potentially observable by existing satellites.

preprint2015arXiv

Search for Transient Astrophysical Neutrino Emission with IceCube-DeepCore

We present the results of a search for astrophysical sources of brief transient neutrino emission using IceCube and DeepCore data acquired between May 15th 2012 and April 30th 2013. While the search methods employed in this analysis are similar to those used in previous IceCube point source searches, the data set being examined consists of a sample of predominantly sub-TeV muon neu- trinos from the Northern Sky (-5$^{\circ}$ < δ < 90$^{\circ}$ ) obtained through a novel event selection method. This search represents a first attempt by IceCube to identify astrophysical neutrino sources in this relatively unexplored energy range. The reconstructed direction and time of arrival of neutrino events is used to search for any significant self-correlation in the dataset. The data revealed no significant source of transient neutrino emission. This result has been used to construct limits at timescales ranging from roughly 1$\,$s to 10 days for generic soft-spectra transients. We also present limits on a specific model of neutrino emission from soft jets in core-collapse supernovae.

preprint2015arXiv

Searches for Relativistic Magnetic Monopoles in IceCube

Various extensions of the Standard Model motivate the existence of stable magnetic monopoles that could have been created during an early high-energy epoch of the Universe. These primordial magnetic monopoles would be gradually accelerated by cosmic magnetic fields and could reach high velocities that make them visible in Cherenkov detectors such as IceCube. Equivalently to electrically charged particles, magnetic monopoles produce direct and indirect Cherenkov light while traversing through matter at relativistic velocities. This paper describes searches for relativistic (v>0.76c) and mildly relativistic (v>0.51c) monopoles, each using one year of data taken in 2008/09 and 2011/12 respectively. No monopole candidate was detected. For a velocity above 0.51c the monopole flux is constrained down to a level of 1.55x10^-18 cm-2 s-1 sr-1. This is an improvement of almost two orders of magnitude over previous limits.

preprint2015arXiv

Searches for small-scale anisotropies from neutrino point sources with three years of IceCube data

Recently, IceCube found evidence for a diffuse signal of astrophysical neutrinos in an energy range of $60\,\mathrm{TeV}$ to the $\mathrm{PeV}$-scale. The origin of those events, being a key to understanding the origin of cosmic rays, is still an unsolved question. So far, analyses have not succeeded to resolve the diffuse signal into point-like sources. Searches including a maximum-likelihood-ratio test, based on the reconstructed directions and energies of the detected down- and up-going neutrino candidates, were also performed on IceCube data leading to the exclusion of bright point sources. In this paper, we present two methods to search for faint neutrino point sources in three years of IceCube data, taken between 2008 and 2011. The first method is an autocorrelation test, applied separately to the northern and southern sky. The second method is a multipole analysis, which expands the measured data in the northern hemisphere into spherical harmonics and uses the resulting expansion coefficients to separate signal from background. With both methods, the results are consistent with the background expectation with a slightly more sparse spatial distribution, corresponding to an underfluctuation. Depending on the assumed number of sources, the resulting upper limit on the flux per source in the northern hemisphere for an $E^{-2}$ energy spectrum ranges from $1.5 \cdot 10^{-8}\,\mathrm{GeV}/(\mathrm{cm}^2 \mathrm{s})$, in the case of one assumed source, to $4 \cdot 10^{-10} \,\mathrm{GeV}/(\mathrm{cm}^2 \mathrm{s})$, in the case of $3500$ assumed sources.

preprint2015arXiv

Searches for Time Dependent Neutrino Sources with IceCube Data from 2008 to 2012

In this paper searches for flaring astrophysical neutrino sources and sources with periodic emission with the IceCube neutrino telescope are presented. In contrast to time integrated searches, where steady emission is assumed, the analyses presented here look for a time dependent signal of neutrinos using the information from the neutrino arrival times to enhance the discovery potential. A search was performed for correlations between neutrino arrival times and directions as well as neutrino emission following time dependent lightcurves, sporadic emission or periodicities of candidate sources. These include active galactic nuclei, soft $γ$-ray repeaters, supernova remnants hosting pulsars, micro-quasars and X-ray binaries. The work presented here updates and extends previously published results to a longer period that covers four years of data from 2008 April 5 to 2012 May 16 including the first year of operation of the completed 86-string detector. The analyses did not find any significant time dependent point sources of neutrinos and the results were used to set upper limits on the neutrino flux from source candidates.

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.

preprint2014arXiv

Energy Reconstruction Methods in the IceCube Neutrino Telescope

Accurate measurement of neutrino energies is essential to many of the scientific goals of large-volume neutrino telescopes. The fundamental observable in such detectors is the Cherenkov light produced by the transit through a medium of charged particles created in neutrino interactions. The amount of light emitted is proportional to the deposited energy, which is approximately equal to the neutrino energy for $ν_e$ and $ν_μ$ charged-current interactions and can be used to set a lower bound on neutrino energies and to measure neutrino spectra statistically in other channels. Here we describe methods and performance of reconstructing charged-particle energies and topologies from the observed Cherenkov light yield, including techniques to measure the energies of uncontained muon tracks, achieving average uncertainties in electromagnetic-equivalent deposited energy of $\sim 15\%$ above 10 TeV.

preprint2014arXiv

IceCube Search for Dark Matter Annihilation in nearby Galaxies and Galaxy Clusters

We present the results of a first search for self-annihilating dark matter in nearby galaxies and galaxy clusters using a sample of high-energy neutrinos acquired in 339.8 days of live time during 2009/10 with the IceCube neutrino observatory in its 59-string configuration. The targets of interest include the Virgo and Coma galaxy clusters, the Andromeda galaxy, and several dwarf galaxies. We obtain upper limits on the cross section as a function of the weakly interacting massive particle mass between 300 GeV and 100 TeV for the annihilation into b bbar, W+W-, τ+τ-, μ+μ-, and ννbar. A limit derived for the Virgo cluster, when assuming a large effect from subhalos, challenges the weak interacting massive particle interpretation of a recently observed GeV positron excess in cosmic rays.

preprint2014arXiv

IceCube-Gen2: A Vision for the Future of Neutrino Astronomy in Antarctica

The recent observation by the IceCube neutrino observatory of an astrophysical flux of neutrinos represents the "first light" in the nascent field of neutrino astronomy. The observed diffuse neutrino flux seems to suggest a much larger level of hadronic activity in the non-thermal universe than previously thought and suggests a rich discovery potential for a larger neutrino observatory. This document presents a vision for an substantial expansion of the current IceCube detector, IceCube-Gen2, including the aim of instrumenting a $10\,\mathrm{km}^3$ volume of clear glacial ice at the South Pole to deliver substantial increases in the astrophysical neutrino sample for all flavors. A detector of this size would have a rich physics program with the goal to resolve the sources of these astrophysical neutrinos, discover GZK neutrinos, and be a leading observatory in future multi-messenger astronomy programs.

preprint2014arXiv

Multimessenger Search for Sources of Gravitational Waves and High-Energy Neutrinos: Results for Initial LIGO-Virgo and IceCube

We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find no significant coincident events, and use the search results to derive upper limits on the rate of joint sources for a range of source emission parameters. For the optimistic assumption of gravitational-wave emission energy of $10^{-2}$\,M$_\odot$c$^2$ at $\sim 150$\,Hz with $\sim 60$\,ms duration, and high-energy neutrino emission of $10^{51}$\,erg comparable to the isotropic gamma-ray energy of gamma-ray bursts, we limit the source rate below $1.6 \times 10^{-2}$\,Mpc$^{-3}$yr$^{-1}$. We also examine how combining information from gravitational waves and neutrinos will aid discovery in the advanced gravitational-wave detector era.

preprint2014arXiv

Multipole analysis of IceCube data to search for dark matter accumulated in the Galactic halo

Dark matter which is bound in the Galactic halo might self-annihilate and produce a flux of stable final state particles, e.g. high energy neutrinos. These neutrinos can be detected with IceCube, a cubic-kilometer sized Cherenkov detector. Given IceCube's large field of view, a characteristic anisotropy of the additional neutrino flux is expected. In this paper we describe a multipole method to search for such a large-scale anisotropy in IceCube data. This method uses the expansion coefficients of a multipole expansion of neutrino arrival directions and incorporates signal-specific weights for each expansion coefficient. We apply the technique to a high-purity muon neutrino sample from the Northern Hemisphere. The final result is compatible with the null-hypothesis. As no signal was observed, we present limits on the self-annihilation cross-section averaged over the relative velocity distribution $<σv>$ down to $1.9\cdot 10^{-23}\,\mathrm{cm}^3\mathrm{s}^{-1}$ for a dark matter particle mass of $700\,\mathrm{GeV}$ to $1000\,\mathrm{GeV}$ and direct annihilation into $ν\barν$. The resulting exclusion limits come close to exclusion limits from $γ$-ray experiments, that focus on the outer Galactic halo, for high dark matter masses of a few TeV and hard annihilation channels.

preprint2014arXiv

Observation of High-Energy Astrophysical Neutrinos in Three Years of IceCube Data

A search for high-energy neutrinos interacting within the IceCube detector between 2010 and 2012 provided the first evidence for a high-energy neutrino flux of extraterrestrial origin. Results from an analysis using the same methods with a third year (2012-2013) of data from the complete IceCube detector are consistent with the previously reported astrophysical flux in the 100 TeV - PeV range at the level of $10^{-8}\, \mathrm{GeV}\, \mathrm{cm}^{-2}\, \mathrm{s}^{-1}\, \mathrm{sr}^{-1}$ per flavor and reject a purely atmospheric explanation for the combined 3-year data at $5.7 σ$. The data are consistent with expectations for equal fluxes of all three neutrino flavors and with isotropic arrival directions, suggesting either numerous or spatially extended sources. The three-year dataset, with a livetime of 988 days, contains a total of 37 neutrino candidate events with deposited energies ranging from 30 to 2000 TeV. The 2000 TeV event is the highest-energy neutrino interaction ever observed.

preprint2014arXiv

Search for non-relativistic Magnetic Monopoles with IceCube

The IceCube Neutrino Observatory is a large Cherenkov detector instrumenting $1\,\mathrm{km}^3$ of Antarctic ice. The detector can be used to search for signatures of particle physics beyond the Standard Model. Here, we describe the search for non-relativistic, magnetic monopoles as remnants of the GUT (Grand Unified Theory) era shortly after the Big Bang. These monopoles may catalyze the decay of nucleons via the Rubakov-Callan effect with a cross section suggested to be in the range of $10^{-27}\,\mathrm{cm^2}$ to $10^{-21}\,\mathrm{cm^2}$. In IceCube, the Cherenkov light from nucleon decays along the monopole trajectory would produce a characteristic hit pattern. This paper presents the results of an analysis of first data taken from May 2011 until May 2012 with a dedicated slow-particle trigger for DeepCore, a subdetector of IceCube. A second analysis provides better sensitivity for the brightest non-relativistic monopoles using data taken from May 2009 until May 2010. In both analyses no monopole signal was observed. For catalysis cross sections of $10^{-22}\,(10^{-24})\,\mathrm{cm^2}$ the flux of non-relativistic GUT monopoles is constrained up to a level of $Φ_{90} \le 10^{-18}\,(10^{-17})\,\mathrm{cm^{-2}s^{-1}sr^{-1}}$ at a 90% confidence level, which is three orders of magnitude below the Parker bound. The limits assume a dominant decay of the proton into a positron and a neutral pion. These results improve the current best experimental limits by one to two orders of magnitude, for a wide range of assumed speeds and catalysis cross sections.

preprint2014arXiv

Searches for Extended and Point-like Neutrino Sources with Four Years of IceCube Data

We present results on searches for point-like sources of neutrinos using four years of IceCube data, including the first year of data from the completed 86-string detector. The total livetime of the combined dataset is 1,373 days. For an E$^{-2}$ spectrum the median sensitivity at 90\% C.L. is $\sim 10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$ for energies between 1 TeV$-$1 PeV in the northern sky and $\sim 10^{-11}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$ for energies between 100 TeV $-$ 100 PeV in the southern sky. The sensitivity has improved from both the additional year of data and the introduction of improved reconstructions compared to previous publications. In addition, we present the first results from an all-sky search for extended sources of neutrinos. We update results of searches for neutrino emission from stacked catalogs of sources, and test five new catalogs; two of Galactic supernova remnants and three of active galactic nuclei. In all cases, the data are compatible with the background-only hypothesis, and upper limits on the flux of muon neutrinos are reported for the sources considered.

preprint2014arXiv

The IceProd Framework: Distributed Data Processing for the IceCube Neutrino Observatory

IceCube is a one-gigaton instrument located at the geographic South Pole, designed to detect cosmic neutrinos, iden- tify the particle nature of dark matter, and study high-energy neutrinos themselves. Simulation of the IceCube detector and processing of data require a significant amount of computational resources. IceProd is a distributed management system based on Python, XML-RPC and GridFTP. It is driven by a central database in order to coordinate and admin- ister production of simulations and processing of data produced by the IceCube detector. IceProd runs as a separate layer on top of other middleware and can take advantage of a variety of computing resources, including grids and batch systems such as CREAM, Condor, and PBS. This is accomplished by a set of dedicated daemons that process job submission in a coordinated fashion through the use of middleware plugins that serve to abstract the details of job submission and job management from the framework.

preprint2013arXiv

Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector

We report on results of an all-sky search for high-energy neutrino events interacting within the IceCube neutrino detector conducted between May 2010 and May 2012. The search follows up on the previous detection of two PeV neutrino events, with improved sensitivity and extended energy coverage down to approximately 30 TeV. Twenty-six additional events were observed, substantially more than expected from atmospheric backgrounds. Combined, both searches reject a purely atmospheric origin for the twenty-eight events at the $4σ$ level. These twenty-eight events, which include the highest energy neutrinos ever observed, have flavors, directions, and energies inconsistent with those expected from the atmospheric muon and neutrino backgrounds. These properties are, however, consistent with generic predictions for an additional component of extraterrestrial origin.

preprint2013arXiv

First observation of PeV-energy neutrinos with IceCube

We report on the observation of two neutrino-induced events which have an estimated deposited energy in the IceCube detector of 1.04 $\pm$ 0.16 and 1.14 $\pm$ 0.17 PeV, respectively, the highest neutrino energies observed so far. These events are consistent with fully contained particle showers induced by neutral-current $ν_{e,μ,τ}$ ($\barν_{e,μ,τ}$) or charged-current $ν_{e}$ ($\barν_{e}$) interactions within the IceCube detector. The events were discovered in a search for ultra-high energy neutrinos using data corresponding to 615.9 days effective livetime. The expected number of atmospheric background is $0.082 \pm 0.004 \text{(stat)}^{+0.041}_{-0.057} \text{(syst)}$. The probability to observe two or more candidate events under the atmospheric background-only hypothesis is $2.9\times10^{-3}$ ($2.8σ$) taking into account the uncertainty on the expected number of background events. These two events could be a first indication of an astrophysical neutrino flux, the moderate significance, however, does not permit a definitive conclusion at this time.

preprint2013arXiv

Improvement in Fast Particle Track Reconstruction with Robust Statistics

The IceCube project has transformed one cubic kilometer of deep natural Antarctic ice into a Cherenkov detector. Muon neutrinos are detected and their direction inferred by mapping the light produced by the secondary muon track inside the volume instrumented with photomultipliers. Reconstructing the muon track from the observed light is challenging due to noise, light scattering in the ice medium, and the possibility of simultaneously having multiple muons inside the detector, resulting from the large flux of cosmic ray muons. This manuscript describes work on two problems: (1) the track reconstruction problem, in which, given a set of observations, the goal is to recover the track of a muon; and (2) the coincident event problem, which is to determine how many muons are active in the detector during a time window. Rather than solving these problems by developing more complex physical models that are applied at later stages of the analysis, our approach is to augment the detectors early reconstruction with data filters and robust statistical techniques. These can be implemented at the level of on-line reconstruction and, therefore, improve all subsequent reconstructions. Using the metric of median angular resolution, a standard metric for track reconstruction, we improve the accuracy in the initial reconstruction direction by 13%. We also present improvements in measuring the number of muons in coincident events: we can accurately determine the number of muons 98% of the time.

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

Observation of the cosmic-ray shadow of the Moon with IceCube

We report on the observation of a significant deficit of cosmic rays from the direction of the Moon with the IceCube detector. The study of this "Moon shadow" is used to characterize the angular resolution and absolute pointing capabilities of the detector. The detection is based on data taken in two periods before the completion of the detector: between April 2008 and May 2009, when IceCube operated in a partial configuration with 40 detector strings deployed in the South Pole ice, and between May 2009 and May 2010 when the detector operated with 59 strings. Using two independent analysis methods, the Moon shadow has been observed to high significance (> 6 sigma) in both detector configurations. The observed location of the shadow center is within 0.2 degrees of its expected position when geomagnetic deflection effects are taken into account. This measurement validates the directional reconstruction capabilities of IceCube.

preprint2013arXiv

Probing the origin of cosmic-rays with extremely high energy neutrinos using the IceCube Observatory

We have searched for extremely high energy neutrinos using data taken with the IceCube detector between May 2010 and May 2012. Two neutrino induced particle shower events with energies around 1 PeV were observed, as reported previously. In this work, we investigate whether these events could originate from cosmogenic neutrinos produced in the interactions of ultra-high energy cosmic-rays with ambient photons while propagating through intergalactic space. Exploiting IceCube's large exposure for extremely high energy neutrinos and the lack of observed events above 100 PeV, we can rule out the corresponding models at more than 90% confidence level. The model independent quasi-differential 90% CL upper limit, which amounts to $E^2 ϕ_{ν_e + ν_μ+ ν_τ} = 1.2 \times 10^{-7}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ at 1 EeV, provides the most stringent constraint in the energy range from 10 PeV to 10 EeV. Our observation disfavors strong cosmological evolution of the highest energy cosmic ray sources such as the Fanaroff-Riley type II class of radio galaxies.

preprint2013arXiv

Search for a diffuse flux of astrophysical muon neutrinos with the IceCube 59-string configuration

A search for high-energy neutrinos was performed using data collected by the IceCube Neutrino Observatory from May 2009 to May 2010, when the array was running in its 59-string configuration. The data sample was optimized to contain muon neutrino induced events with a background contamination of atmospheric muons of less than 1%. These data, which are dominated by atmospheric neutrinos, are analyzed with a global likelihood fit to search for possible contributions of prompt atmospheric and astrophysical neutrinos, neither of which have yet been identified. Such signals are expected to follow a harder energy spectrum than conventional atmospheric neutrinos. In addition, the zenith angle distribution differs for astrophysical and atmospheric signals. A global fit of the reconstructed energies and directions of observed events is performed, including possible neutrino flux contributions for an astrophysical signal and atmospheric backgrounds as well as systematic uncertainties of the experiment and theoretical predictions. The best fit yields an astrophysical signal flux for $ν_μ+ \barν_μ$ of $E^2 \cdot Φ(E) = 0.25 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$, and a zero prompt component. Although the sensitivity of this analysis for astrophysical neutrinos surpasses the Waxman and Bahcall upper bound, the experimental limit at 90% confidence level is a factor of 1.5 above at a flux of $E^2 \cdot Φ(E) = 1.44 \cdot 10^{-8} \textrm{GeV} \textrm{cm}^{-2} \textrm{s}^{-1} \textrm{sr}^{-1}$.

preprint2013arXiv

Search for neutrino-induced particle showers with IceCube-40

We report on the search for neutrino-induced particle-showers, so-called cascades, in the IceCube-40 detector. The data for this search was collected between April 2008 and May 2009 when the first 40 IceCube strings were deployed and operational. Three complementary searches were performed, each optimized for different energy regimes. The analysis with the lowest energy threshold (2 TeV) targeted atmospheric neutrinos. A total of 67 events were found, consistent with the expectation of 41 atmospheric muons and 30 atmospheric neutrino events. The two other analyses targeted a harder, astrophysical neutrino flux. The analysis with an intermediate threshold of 25 TeV lead to the observation of 14 cascade-like events, again consistent with the prediction of 3.0 atmospheric neutrino and 7.7 atmospheric muon events. We hence set an upper limit of $E^2 Φ_{lim} \leq 7.46\times10^{-8}\,\mathrm{GeV sr^{-1} s^{-1} cm^{-2}}$ (90% C.L.) on the diffuse flux from astrophysical neutrinos of all neutrino flavors, applicable to the energy range 25 TeV to 5 PeV, assuming an $E_ν^{-2}$ spectrum and a neutrino flavor ratio of 1:1:1 at the Earth. The third analysis utilized a larger and optimized sample of atmospheric muon background simulation, leading to a higher energy threshold of 100 TeV. Three events were found over a background prediction of 0.04 atmospheric muon events and 0.21 events from the flux of conventional and prompt atmospheric neutrinos. Including systematic errors this corresponds to a $2.7σ$ excess with respect to the background-only hypothesis. Our observation of neutrino event candidates above 100 TeV complements IceCube's recently observed evidence for high-energy astrophysical neutrinos.

preprint2013arXiv

Search for time-independent neutrino emission from astrophysical sources with 3 years of IceCube data

We present the results of a search for neutrino point sources using the IceCube data collected between April 2008 and May 2011 with three partially completed configurations of the detector: the 40-, 59- and 79-string configurations. The live-time of this data set are 1,040 days. An unbinned maximum likelihood ratio test was used to search for an excess of neutrinos above the atmospheric background at any given direction in the sky. By adding two more years of data with improved event selection and reconstruction techniques, the sensitivity was improved by a factor 3.5 or more with respect to the previously published results obtained with the 40-string configuration of IceCube. We performed an all-sky survey and a dedicated search using a catalog of \textit{a priori} selected objects observed by other telescopes. In both searches, the data are compatible with the background-only hypothesis. In the absence of evidence for a signal, we set upper limits on the flux of muon neutrinos. For an E$^{-2}$ neutrino spectrum, the observed limits are between 0.9 and $23.2\times 10^{-12}$ TeV$^{-1}$ cm$^{-2}$s$^{-1}$. We also report upper limits for neutrino emission from groups of sources which were selected according to theoretical models or observational parameters and analysed with a stacking approach.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Observation of Cosmic Ray Anisotropy with the IceTop Air Shower Array

We report on the observation of anisotropy in the arrival direction distribution of cosmic rays at PeV energies. The analysis is based on data taken between 2009 and 2012 with the IceTop air shower array at the South Pole. IceTop, an integral part of the IceCube detector, is sensitive to cosmic rays between 100 TeV and 1 EeV. With the current size of the IceTop data set, searches for anisotropy at the 10^-3 level can, for the first time, be extended to PeV energies. We divide the data set into two parts with median energies of 400 TeV and 2 PeV, respectively. In the low energy band, we observe a strong deficit with an angular size of about 30 degrees and an amplitude of (-1.58 +/- 0.46 (stat) +/- 0.52 (sys)) x 10^(-3) at a location consistent with previous observations of cosmic rays with the IceCube neutrino detector. The study of the high energy band shows that the anisotropy persists to PeV energies and increases in amplitude to (-3.11 +/- 0.38 (stat) +/- 0.96 (sys)) x 10^(-3).

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

The Rapid Atmospheric Monitoring System of the Pierre Auger Observatory

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray Showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

preprint2011arXiv

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

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

preprint2010arXiv

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

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

preprint2009arXiv

The Fluorescence Detector of the Pierre Auger Observatory

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

preprint2009arXiv

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

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