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¹⁸ eV). Despite the flux of these particles being extremely low, the area of ∼3000 km² 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² 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.

Despite the uncovered association of a high-energy neutrino with the apparent flaring state of blazar TXS 0506+056 in 2017, the mechanisms leading to astrophysical particle acceleration and neutrino production are still uncertain. Recent studies found that when transparent to γ-rays, γ-flaring blazars do not have the opacity for protons to produce neutrinos. Here we present observational evidence for an alternative explanation, in which γ-ray emission is suppressed during efficient neutrino production. A large proton and target photon density helps produce neutrinos while temporarily suppressing the observable γ-emission due to a large γ γ opacity. We show that the Fermi-LAT γ-flux of blazar PKS 1502+106 was at a local minimum when IceCube recorded the coincident high-energy neutrino IC-190730A. Using data from the OVRO 40 m Telescope, we find that radio emission from PKS 1502+106 at the time period of the coincident neutrino IC-190730A was in a high state, in contrast to earlier time periods when radio and γ fluxes are correlated for both low and high states. This points to an active outflow that is γ-suppressed at the time of neutrino production. We find similar local γ-suppression in other blazars, including in MAGIC’s TeV flux of TXS 0506+056 and Fermi-LAT’s flux of blazar PKS B1424-418 at the time of coincident IceCube neutrino detections. Using temporary γ-suppression, neutrino–blazar coincidence searches could be substantially more sensitive than previously assumed, enabling the identification of the origin of IceCube’s diffuse neutrino flux possibly with already existing data.

The statistics of black holes and their masses strongly suggests that their mass distribution has a cut-off towards lower masses near 3 × 10⁶ M_⊙. This is consistent with a classical formation mechanism from the agglomeration of the first massive stars in the universe. However, when the masses of the stars approach 10⁶ M_⊙, the stars become unstable and collapse, possibly forming the first generation of cosmological black holes. Here, we speculate that the claimed detection of an isotropic radio background may constitute evidence of the formation of these first supermassive black holes, since their data are compatible in spectrum and intensity with synchrotron emission from the remnants. The model proposed fulfils all observational conditions for the background, in terms of single-source strength, number of sources, far-infrared and gamma-ray emission. The observed high-energy neutrino flux is consistent with our calculations in flux and spectrum. The proposal described in this paper may also explain the early formation and growth of massive bulge-less disc galaxies as derived from the massive, gaseous shell formed during the explosion prior to the formation of a supermassive black hole.

LOPES, the LOFAR prototype station, was an antenna array for cosmic-ray air showers operating from 2003 to 2013 within the KASCADE-Grande experiment. Meanwhile, the analysis is finished and the data of air-shower events measured by LOPES are available with open access in the KASCADE Cosmic Ray Data Center (KCDC). This article intends to provide a summary of the achievements, results, and lessons learned from LOPES. By digital, interferometric beamforming the detection of air showers became possible in the radio-loud environment of the Karlsruhe Institute of Technology (KIT). As a prototype experiment, LOPES tested several antenna types, array configurations and calibration techniques, and pioneered analysis methods for the reconstruction of the most important shower parameters, i.e., the arrival direction, the energy, and mass-dependent observables such as the position of the shower maximum. In addition to a review and update of previously published results, we also present new results based on end-to-end simulations including all known instrumental properties. For this, we applied the detector response to radio signals simulated with the CoREAS extension of CORSIKA, and analyzed them in the same way as measured data. Thus, we were able to study the detector performance more accurately than before, including some previously inaccessible features such as the impact of noise on the interferometric cross-correlation beam. These results led to several improvements, which are documented in this paper and can provide useful input for the design of future cosmic-ray experiments based on the digital radio-detection technique.

A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 2004 January 1 and 2020 December 31 are given for cosmic rays that have energies in the range 78–166 EeV. Details are also given on a further nine very energetic events that have been used in the calibration procedure adopted to determine the energy of each primary. A sky plot of the arrival directions of the most energetic particles is shown. No interpretations of the data are offered.