UHECR doublets and their conditional association with nearby radio galaxies

TL;DR

This study identifies significant associations between ultra-high-energy cosmic ray doublets and nearby radio galaxies, especially Fornax A, suggesting these galaxies as long-term accelerators of heavy UHECRs based on 16 years of Pierre Auger Observatory data.

Contribution

Introduces a novel spatiotemporal multiplet search method and demonstrates a strong statistical link between UHECRs and radio galaxies, highlighting Fornax A as a key source.

Findings

28 UHECR multiplets identified above 32 EeV

Post-trial significance of 5.8 sigma for association with radio galaxies

4.5 sigma significance for association with Fornax A

Abstract

The origin of ultra-high-energy cosmic rays (UHECRs) remains a fundamental question in astroparticle physics. While localized 3 $\sigma$ correlations with active galactic nuclei and starburst galaxies have been reported using time-integrated analyses, we propose and implement a spatiotemporal multiplet search method utilizing a pre-defined fixed window of 3 degrees and 15 days, a kinematic filter designed to isolate high-rigidity particles and keep chance coincidences low. Applying this method to 16 years of Pierre Auger Observatory data, we identify 28 UHECR multiplets (doublets) above 32 EeV. We backtrack these trajectories using three Galactic magnetic field models across eight distinct nuclear species. Testing the backtracked directions against the ten nearest bright radio galaxies yields an overall post-trial significance of 5.8 $\sigma$, conditioned on the individual best-fit associations. Specifically, we find a 4.5 $\sigma$ conditional post-trial significance for the joint association of 8 of these multiplets with the Fornax A region alone. These results point to radio galaxies, with a strong contribution from Fornax A, as long-term accelerators of heavy (Z > 3) UHECRs, possibly within the mildly relativistic backflows of their extended radio lobes, detected at Earth primarily as independent secondary fragments above 32 EeV.