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.

Cosmic rays are atomic nuclei arriving from outer space that reach the highest energies observed in nature. Clues to their origin come from studying the distribution of their arrival directions. Using 3 × 10 cosmic rays with energies above 8 × 10 electron volts, recorded with the Pierre Auger Observatory from a total exposure of 76,800 km sr year, we determined the existence of anisotropy in arrival directions. The anisotropy, detected at more than a 5.2σ level of significance, can be described by a dipole with an amplitude of [Formula: see text] percent toward right ascension α = 100 ± 10 degrees and declination δ = [Formula: see text] degrees That direction indicates an extragalactic origin for these ultrahigh-energy particles.

Neutrinos with energies above 1017 eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming τ neutrinos with nearly tangential trajectories relative to the Earth. No neutrino candidates were found in ∼ 14.7 years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The 90% C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an Eν−2 spectrum in the energy range 1.0 × 1017 eV –2.5 × 1019 eV is E2 dNν/dEν < 4.4 × 10−9 GeV cm−2 s−1 sr−1, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

We report a measurement of the energy spectrum of cosmic rays for energies above 2.5 × 10 18 eV based on 215,030 events recorded with zenith angles below 60°. A key feature of the work is that the estimates of the energies are independent of assumptions about the unknown hadronic physics or of the primary mass composition. The measurement is the most precise made hitherto with the accumulated exposure being so large that the measurements of the flux are dominated by systematic uncertainties except at energies above 5 × 10 19 eV. The principal conclusions are (1) The flattening of the spectrum near 5 × 10 18 eV, the so-called "ankle," is confirmed. (2) The steepening of the spectrum at around 5 × 10 19 eV is confirmed. (3) A new feature has been identified in the spectrum: in the region above the ankle the spectral index γ of the particle flux (∝ E −γ) changes from 2.51 AE 0.03 ðstatÞ AE 0.05 ðsystÞ to 3.05 AE 0.05 ðstatÞ AE 0.10 ðsystÞ before changing sharply to 5.1 AE 0.3 ðstatÞ AE 0.1 ðsystÞ above 5 × 10 19 eV. (4) No evidence for any dependence of the spectrum on declination has been found other than a mild excess from the Southern Hemisphere that is consistent with the anisotropy observed above 8 × 10 18 eV.

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

4We present a new method for probing the hadronic interaction models at ultra-high 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, firstly, inconsistencies between our data and predictions over a greater energy range and with substantially more events than in previous studies. Secondly, by calibrating the new parameter with fluorescence measurements from observations made at the Auger Observatory, we can infer the depth of shower maximum X max for a sample of over 81,000 events extending from 0.3 EeV to over 100 EeV. Above 30 EeV, the sample is nearly fourteen times larger than currently available from fluorescence measurements and extending the covered energy range by half a decade. The energy dependence of X max 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. 96.50.sb, 96.50.sd, 98.70.Sa

A search for ultra-high energy photons with energies above 1 EeV is performed using nine years of data collected by the Pierre Auger Observatory in hybrid operation mode. An unprecedented separation power between photon and hadron primaries is achieved by combining measurements of the longitudinal air-shower development with the particle content at ground measured by the fluorescence and surface detectors, respectively. Only three photon candidates at energies 1 − 2 EeV are found, which is compatible with the expected hadroninduced background. Upper limits on the integral flux of ultra-high energy photons of 0.027, 0.009, 0.008, 0.008 and 0.007 km −2 sr −1 yr −1 are derived at 95% C.L. for energy thresholds of 1, 2, 3, 5 and 10 EeV. These limits bound the fractions of photons in the all-particle integral flux below 0.1%, 0.15%, 0.33%, 0.85% and 2.7%. For the first time the photon fraction at EeV energies is constrained at the sub-percent level. The improved limits are below the flux of diffuse photons predicted by some astrophysical scenarios for cosmogenic photon production. The new results rule-out the early top-down models − in which ultra-high energy cosmic rays are produced by, e.g., the decay of super-massive particles − and challenge the most recent super-heavy dark matter models.

Context. Nearly a dozen star-forming galaxies have been detected in γ-rays by the Fermi observatory in the last decade. A remarkable property of this sample is the quasi-linear relation between the γ-ray luminosity and the star formation rate, which was obtained assuming that the latter is well traced by the infrared luminosity of the galaxies. The non-linearity of this relation has not been fully explained yet. Aims. We aim to determine the biases derived from the use of the infrared luminosity as a proxy for the star formation rate and to shed light on the more fundamental relation between the latter and the γ-ray luminosity. We expect to quantify and explain some trends observed in this relation. Methods. We compiled a near-homogeneous set of distances, ultraviolet, optical, infrared, and γ-ray fluxes from the literature for all known γ-ray emitting, star-forming galaxies. From these data, we computed the infrared and γ-ray luminosities, and star formation rates. We determined the best-fitting relation between the latter two, and we describe the trend using simple, population-orientated models for cosmic-ray transport and cooling. Results. We find that the γ-ray luminosity–star formation rate relation obtained from infrared luminosities is biased to shallower slopes. The actual relation is steeper than previous estimates, having a power-law index of 1.35 ± 0.05, in contrast to 1.23 ± 0.06. Conclusions. The unbiased γ-ray luminosity–star formation rate relation can be explained at high star formation rates by assuming that the cosmic-ray cooling region is kiloparsec-sized and pervaded by mild to fast winds. Combined with previous results about the scaling of wind velocity with star formation rate, our work provides support to advection as the dominant cosmic-ray escape mechanism in galaxies with low star formation rates.