Abstract
A promising energy range to look for angular correlations between cosmic rays of extragalactic origin and their sources is at the highest energies, above a few tens of EeV (1 EeV ≡ 10
18
...eV). Despite the flux of these particles being extremely low, the area of ∼3000 km
2
covered at the Pierre Auger Observatory, and the 17 yr data-taking period of the
Phase
1
of its operations, have enabled us to measure the arrival directions of more than 2600 ultra-high-energy cosmic rays above 32 EeV. We publish this data set, the largest available at such energies from an integrated exposure of 122,000 km
2
sr yr, and search it for anisotropies over the 3.4
π
steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scales, with ∼15° Gaussian spread or ∼25° top-hat radius, is obtained at the 4
σ
significance level for cosmic-ray energies above ∼40 EeV.
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 Anti-Coincidence 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 500 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.
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.
Using the data taken at the Pierre Auger Observatory between December 2004 and December 2012, we have examined the implications of the distributions of depths of atmospheric shower maximum (X ...sub(max)), using a hybrid technique, for composition and hadronic interaction models. We do this by fitting the distributions with predictions from a variety of hadronic interaction models for variations in the composition of the primary cosmic rays and examining the quality of the fit. Regardless of what interaction model is assumed, we find that our data are not well described by a mix of protons and iron nuclei over most of the energy range. Acceptable fits can be obtained when intermediate masses are included, and when this is done consistent results for the proton and iron-nuclei contributions can be found using the available models. We observe a strong energy dependence of the resulting proton fractions, and find no support from any of the models for a significant contribution from iron nuclei. However, we also observe a significant disagreement between the models with respect to the relative contributions of the intermediate components.
Abstract
In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 10
17
eV. We use an astrophysical model ...with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 10
18
eV (the so-called “ankle” feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources.
To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).
Abstract
A new analysis of the data set from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, ...which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above
with zenith angles up to 80° recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard
Fermi
-LAT sources (2FHL) and starburst galaxies from a sample that was examined with
Fermi
-LAT. Flux-limited samples, which include all types of galaxies from the
Swift
-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects, and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0
σ
, the highest value of the test statistic being for energies above
. The three alternative models are favored against isotropy with 2.7
σ
–3.2
σ
significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.
We present a new method for probing the hadronic interaction models at ultrahigh energy and extracting details about mass composition. This is done using the time profiles of the signals recorded ...with the water-Cherenkov detectors of the Pierre Auger Observatory. The profiles arise from a mix of the muon and electromagnetic components of air showers. Using the risetimes of the recorded signals, we define a new parameter, which we use to compare our observations with predictions from simulations. We find, first, inconsistencies between our data and predictions over a greater energy range and with substantially more events than in previous studies. Second, by calibrating the new parameter with fluorescence measurements from observations made at the Auger Observatory, we can infer the depth of shower maximum Xmax for a sample of over 81,000 events extending from 0.3 to over 100 EeV. Above 30 EeV, the sample is nearly 14 times larger than what is currently available from fluorescence measurements and extending the covered energy range by half a decade. The energy dependence of ⟨Xmax⟩ is compared to simulations and interpreted in terms of the mean of the logarithmic mass. We find good agreement with previous work and extend the measurement of the mean depth of shower maximum to greater energies than before, reducing significantly the statistical uncertainty associated with the inferences about mass composition.
We present the first measurement of the fluctuations in the number of muons in extensive air showers produced by ultrahigh energy cosmic rays. We find that the measured fluctuations are in good ...agreement with predictions from air shower simulations. This observation provides new insights into the origin of the previously reported deficit of muons in air shower simulations and constrains models of hadronic interactions at ultrahigh energies. Our measurement is compatible with the muon deficit originating from small deviations in the predictions from hadronic interaction models of particle production that accumulate as the showers develop.
The Pierre Auger collaboration reports new results bearing on the composition of cosmic rays. The muon number of air showers, created by cosmic rays and measured by the collaboration, is intriguingly ...at odds with all theoretical models, posing a challenge to our current understanding of the mass composition of cosmic rays. We present the first hybrid measurement of the average muon number in air showers at ultrahigh energies, initiated by cosmic rays with zenith angles between 62degrees and 80degrees. The measurement is based on 174 hybrid events recorded simultaneously with the surface detector array and the fluorescence detector of the Pierre Auger Observatory. The muon number for each shower is derived by scaling a simulated reference profile of the lateral muon density distribution at the ground until it fits the data. A 10 super(19) eV shower with a zenith angle of 67degrees, which arrives at the surface detector array at an altitude of 1450 m above sea level, contains on average (2.68 + or - 0.04 + or - 0.48(sys)) x 10 super(7) muons with energies larger than 0.3 GeV. The logarithmic gain d ln N sub( mu )/d ln E of muons with increasing energy between 4 X 10 super(18) eV and 5 x 10 super(19) eV is measured to be (1.029 + or - 0.024 + or - 0.030(sys)).
We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins 4, 8 EeV and E ≥ 8 EeV, the most ...significant signal is a dipolar modulation in R.A. at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index β = 0.79 0.19. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally, we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.