A search for anisotropy in the arrival directions of ultra high energy cosmic rays recorded at the Pierre Auger Observatory

Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

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

Deep observation of the NGC1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse γ-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E ≥ 100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6º from the centre, have been detected as point-like VHE γ-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index Γ = 3.6 ± 0.2_(stat) ± 0.2_(syst) between 90 GeV and 1200 GeV. We do not detect any diffuse γ-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes α. For α = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be ≤ 1 − 2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to ≤ 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B_(0), that depend on the CR distribution. For α = 2.2, B_(0) ≥ 5 − 8 µG, which is below the ∼25 µG inferred from Faraday rotation measurements, whereas for α ≤ 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our γ-ray flux upper limits independently of the magnetic field strength.

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 x 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. (C) 2011 Elsevier B.V. All rights reserved.

Energy estimation of cosmic rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30-80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy-corrected for geometrical effects-is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The impact of the air-fluorescence yield on the reconstructed shower parameters of ultra-high energy cosmic rays.

An accurate knowledge of the fluorescence yield and its dependence on atmospheric properties such as pressure, temperature or humidity is essential to obtain a reliable measurement of the primary energy of cosmic rays in experiments using the fluorescence technique. In this work, several sets of fluorescence yield data (i.e. absolute value and quenching parameters) are described and compared. A simple procedure to study the effect of the assumed fluorescence yield on the reconstructed shower parameters (energy and shower maximum depth) as a function of the primary features has been developed. As an application, the effect of water vapor and temperature dependence of the collisional cross section on the fluorescence yield and its impact on the reconstruction of primary energy and shower maximum depth has been studied. Published by Elsevier B.V.

The Pierre Auger Observatory scaler mode for the study of solar activity modulation of galactic cosmic rays

Since data-taking began in January 2004, the Pierre Auger Observatory has been recording the count rates of low energy secondary cosmic ray particles for the self-calibration of the ground detectors of its surface detector array. After correcting for atmospheric effects, modulations of galactic cosmic rays due to solar activity and transient events are observed. Temporal variations related with the activity of the heliosphere can be determined with high accuracy due to the high total count rates. In this study, the available data are presented together with an analysis focused on the observation of Forbush decreases, where a strong correlation with neutron monitor data is found.

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, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th 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(-6)(+7))%, 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. (C) 2010 Elsevier B.V. All rights reserved.

Measurement of the energy spectrum of cosmic rays above 10(18) eV using the Pierre Auger Observatory

We report a measurement of the flux of cosmic rays with unprecedented precision and Statistics using the Pierre Auger Observatory Based on fluorescence observations in coincidence with at least one Surface detector we derive a spectrum for energies above 10(18) eV We also update the previously published energy spectrum obtained with the surface detector array The two spectra are combined addressing the systematic uncertainties and, in particular. the influence of the energy resolution on the spectral shape The spectrum can be described by a broken power law E-gamma with index gamma = 3 3 below the ankle which is measured at log(10)(E-ankle/eV) = 18 6 Above the ankle the spectrum is described by a power law with index 2 6 followed by a flux suppression, above about log(10)(E/eV) = 19 5, detected with high statistical significance (C) 2010 Elsevier B V All rights reserved

Measurement of the cosmic ray energy spectrum using hybrid events of the Pierre Auger Observatory

The energy spectrum of ultra-high energy cosmic rays above 10(18)eV is measured using the hybrid events collected by the Pierre Auger Observatory between November 2005 and September 2010. The large exposure of the Observatory allows the measurement of the main features of the energy spectrum with high statistics. Full Monte Carlo simulations of the extensive air showers (based on the CORSIKA code) and of the hybrid detector response are adopted here as an independent cross check of the standard analysis (Phys. Lett. B 685, 239 (2010)). The dependence on mass composition and other systematic uncertainties are discussed in detail and, in the full Monte Carlo approach, a region of confidence for flux measurements is defined when all the uncertainties are taken into account. An update is also reported of the energy spectrum obtained by combining the hybrid spectrum and that measured using the surface detector array.

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 degrees, detected at the Pierre Auger Observatory. the geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the similar to 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.

Measurement of the radiation energy in the radio signal of extensive air showers as a universal estimator of cosmic-ray energy.

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

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 theta)(max), sensitive to the mass composition of cosmic rays above 3 x 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 modeling that must be resolved before the mass composition can be inferred from (sec theta)(max).

Search for isotropic γ radiation in the cosmological window between 65 and 200 TeV

Electromagnetic energy injected into the universe above a few hundred TeV is expected to pile up as γ radiation in a relatively narrow energy interval below 100 TeV due to its interaction with the 2.7^°K background radiation. We present an upper limit (90% C.L.) on the ratio of primary γ to charged cosmic rays in the energy interval 65–160 TeV (80–200 TeV) of 10.3 • 10^−3 (7.8 • 10^−3). Data from the HEGRA cosmic-ray detector complex consisting of a wide angle Čerenkov array (AIROBICC) measuring the lateral distribution of air Čerenkov light and a scintillator array, were used with a novel method to discriminate γ-ray and hadron induced air showers. If the presently unmeasured universal far infrared background radiation is not too intense, the result rules out a topological-defect origin of ultrahigh energy cosmic rays for masses of the X particle released by the defects equal to or larger than about 10^16 GeV.

Search for very high energy gamma radiation from the radio bright region DR4 of the SNR G78.2+2.1

Data from the HEGRA air shower array are used to set an upper limit on the emission of gamma-radiation above 25 (18) TeV from the direction of the radio bright region DR4 within the SNR G78.2 + 2.1 of 2.5 (7.1). 10^-13 cm^-2 sec^-1. The shock front of SNR G78.2 + 2.1 probably recently overtook the molecular cloud Gong 8 which then acts as a target for the cosmic rays produced within the SNR, thus leading to the expectation of enhanced gamma-radiation. Using a model of Drury, Aharonian and Völk which assumes that SNRs are the sources of galactic cosmic rays via first order Fermi acceleration, we calculated a theoretical prediction for the gamma-ray flux from the DR4 region and compared it with our experimental flux limit. Our 'best estimate' value for the predicted flux lies a factor of about 18 above the upper limit for gamma-ray energies above 25 TeV. Possible reasons for this discrepancy are discussed.