Previous detections of individual astrophysical sources of neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017, we detected a high-energy neutrino, IceCube-170922A, with an energy of ~290 tera-electron volts. Its arrival direction was consistent with the location of a known γ-ray blazar, TXS 0506+056, observed to be in a flaring state. An extensive multiwavelength campaign followed, ranging from radio frequencies to γ-rays. These observations characterize the variability and energetics of the blazar and include the detection of TXS 0506+056 in very-high-energy γ-rays. This observation of a neutrino in spatial coincidence with a γ-ray-emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

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

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 ≡ 1018 eV). Despite the flux of these particles being extremely low, the area of ∼3000 km2 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 km2 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.

More than half a century after the discovery of ultra-high energy cosmic rays (UHECRs), their origin is still an open question. The study of anisotropies in the arrival directions of such particles is an essential ingredient to solve this puzzle. We update our previous analysis of large-scale anisotropies observed by the Pierre Auger Observatory using the latest data collected before the AugerPrime upgrade. We select events with zenith angles up to 80 degrees, implying a sky coverage of 85%, and energies above 4 EeV, for which the surface detector of the Observatory is fully efficient. Dipolar and quadrupolar amplitudes are evaluated through a combined Fourier analysis of the event count rate in right ascension and azimuth. The analysis is performed in three energy bins with boundaries at 4, 8, 16 and 32 EeV and two additional cumulative bins with energies above 8 and 32 EeV. The most significant signal is a dipolar modulation in right ascension for energies above 8 EeV, as previously reported, with statistical significance of 6.6 . Additionally, we report the measurements of the angular power spectrum for the same energy bins with the same dataset.

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 1017 eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 1018 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).