pith. machine review for the scientific record. sign in

arxiv: 2604.23103 · v1 · submitted 2026-04-25 · 🌌 astro-ph.HE

Recognition: unknown

UHECR doublets and their conditional association with nearby radio galaxies

Victor Barbosa Martins
Authors on Pith no claims yet

Pith reviewed 2026-05-08 07:38 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords ultra-high-energy cosmic raysUHECR multipletsradio galaxiesFornax APierre Auger ObservatoryGalactic magnetic field backtrackingsource identification
0
0 comments X

The pith

Ultra-high-energy cosmic rays above 32 EeV form doublets that backtrack to nearby radio galaxies with 5.8 sigma significance after accounting for Galactic deflections.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a fixed-window search for pairs of ultra-high-energy cosmic rays arriving within 3 degrees and 15 days, applies a rigidity-based filter to suppress random overlaps, and backtracks the events through models of the Milky Way's magnetic field. When the arrival directions are compared to the ten closest bright radio galaxies, the matches reach 5.8 sigma post-trial, driven largely by a 4.5 sigma joint signal from eight events pointing near Fornax A. A sympathetic reader would care because the sources of the highest-energy particles have remained unidentified for decades; localizing them to specific, nearby objects would directly constrain the acceleration sites and the nuclear composition at injection.

Core claim

Applying the 3-degree, 15-day spatiotemporal search plus kinematic filter to 16 years of Pierre Auger data yields 28 UHECR doublets above 32 EeV. Backtracking these trajectories with three Galactic magnetic field models across eight nuclear species and testing against the ten nearest bright radio galaxies produces an overall conditional post-trial significance of 5.8 sigma, including 4.5 sigma for the joint association of eight multiplets with the Fornax A region. The results indicate that radio galaxies accelerate heavy nuclei (Z > 3) over long timescales, most likely inside the mildly relativistic backflows of their extended radio lobes, with the particles reaching Earth as independent had

What carries the argument

The spatiotemporal multiplet search method that applies fixed 3-degree spatial and 15-day temporal windows together with a kinematic filter to retain only high-rigidity particles before backtracking.

If this is right

  • Radio galaxies, especially Fornax A, function as long-term accelerators of heavy UHECRs (Z > 3).
  • These nuclei are observed at Earth primarily as independent secondary fragments above 32 EeV.
  • The associations survive after backtracking through multiple Galactic magnetic field models for a range of nuclear charges.
  • Time-integrated source searches are less sensitive than this fixed-window multiplet approach for the same exposure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the result holds, future larger datasets could test whether additional radio galaxies contribute once statistics improve.
  • Detailed modeling of particle fragmentation and deflection inside radio-lobe backflows could become a direct prediction to confront with composition data.
  • The method's reliance on conditional statistics invites similar time-dependent searches in other high-energy messengers to cross-check source candidates.
  • Confirmation would favor scenarios in which UHECR composition grows heavier with energy, consistent with some propagation models but open to direct test.

Load-bearing premise

The chosen 3-degree and 15-day windows plus the kinematic filter genuinely isolate physical associations instead of random alignments, and conditioning the significance on the individual best-fit associations does not inflate the reported 5.8 sigma value through post-selection bias.

What would settle it

Repeating the identical analysis on an independent UHECR dataset from Telescope Array or with a newer Galactic magnetic field model that produces no comparable post-trial significance for the same or neighboring radio-galaxy associations would falsify the central claim.

read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports a search for UHECR doublets in 16 years of Pierre Auger data above 32 EeV using a pre-defined 3° angular and 15-day temporal window plus a kinematic filter, identifying 28 multiplets. These are backtracked with three Galactic magnetic field models and eight nuclear species; the best-fit directions are then tested against the ten nearest bright radio galaxies, yielding a conditional post-trial significance of 5.8 σ overall (and 4.5 σ for the subset of eight multiplets associated with Fornax A). The authors interpret this as evidence that radio galaxies, especially Fornax A, accelerate heavy (Z > 3) UHECRs.

Significance. If the conditional significances survive a properly constructed null distribution that repeats the full best-fit selection step, the result would be a notable step toward source identification for UHECRs, strengthening the case for radio galaxies over other candidate classes. The analysis is data-driven and uses existing Auger events, but its impact hinges on whether the reported p-values correctly incorporate the optimization over the 24 backtracking combinations.

major comments (2)
  1. [Abstract and statistical methods] Abstract and statistical methods section: The quoted 5.8 σ (and 4.5 σ for Fornax A) post-trial significance is conditioned on selecting, for each of the 28 multiplets, the single best association out of 24 discrete backtracking combinations (3 GMF models × 8 nuclear species). The manuscript does not state whether the Monte Carlo realizations used to evaluate the post-trial p-value repeat this identical per-multiplet best-fit selection procedure. If the null distribution instead fixes the association choice or omits the optimization, the reported significance is inflated by post-selection bias and cannot be taken at face value. This directly affects the central claim.
  2. [Methods and results] Methods and results sections on multiplet selection: The kinematic filter and the choice of exactly ten radio galaxies are presented as fixed and pre-defined, yet the overall pipeline still involves an after-the-fact optimization over 24 backtracking options before the significance is computed. A concrete test is required: the null Monte Carlo must generate synthetic multiplets under the same window and filter, then apply the identical 24-option best-fit selection before counting how often the observed level of association with the ten galaxies (or with Fornax A) is exceeded. Without this explicit repetition, the conditional significance cannot be interpreted as a frequentist p-value.
minor comments (2)
  1. [Abstract and Methods] The abstract and main text should explicitly define the kinematic filter (rigidity cut, energy threshold, etc.) and state whether any of its parameters were tuned on the data.
  2. [Results] Figure captions and tables listing the 28 multiplets should include the chosen nuclear species and GMF model for each best-fit association so that the conditioning step is transparent.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review, which highlights an important point about the construction of the null distribution. We address the two major comments below and will revise the manuscript to make the statistical procedure fully explicit.

read point-by-point responses

  1. Referee: [Abstract and statistical methods] The quoted 5.8 σ post-trial significance is conditioned on selecting the single best association out of 24 backtracking combinations for each multiplet. The manuscript does not state whether the Monte Carlo realizations repeat this identical per-multiplet best-fit selection. If the null omits the optimization, the significance is inflated by post-selection bias.

    Authors: We confirm that the Monte Carlo realizations used to compute the post-trial p-value do repeat the full per-multiplet optimization over the 24 combinations (3 GMF models × 8 species) before testing associations with the radio galaxies. This was done to ensure the null distribution incorporates the same selection step that was applied to the data. We acknowledge that the manuscript text does not state this explicitly and will revise the statistical methods section to describe the procedure in full detail, including a brief pseudocode outline of the MC pipeline. revision: yes

  2. Referee: [Methods and results] The kinematic filter and choice of ten radio galaxies are fixed, yet the pipeline involves after-the-fact optimization over 24 backtracking options. A concrete test is required: the null MC must generate synthetic multiplets under the same window and filter, then apply the identical 24-option best-fit selection before counting associations.

    Authors: The null Monte Carlo does generate synthetic multiplets using the identical 3° × 15-day window and kinematic filter, followed by the full 24-combination backtracking and best-fit selection for each realization before evaluating the conditional association rate with the ten galaxies (or with Fornax A). This matches the data analysis exactly. We will add an explicit paragraph and a short algorithmic description in the methods section to document this step and remove any ambiguity about the frequentist construction of the p-value. revision: yes

Circularity Check

1 steps flagged

Conditional post-trial significance depends on data-driven best-fit selection over 24 backtracking options

specific steps
  1. fitted input called prediction [Abstract]
    "Testing the backtracked directions against the ten nearest bright radio galaxies yields an overall post-trial significance of 5.8 σ, conditioned on the individual best-fit associations. Specifically, we find a 4.5 σ conditional post-trial significance for the joint association of 8 of these multiplets with the Fornax A region alone."

    The significance is computed after selecting, for each of the 28 multiplets, the single best-fit association from 24 backtracking combinations (3 GMF models × 8 nuclear species). Because this choice is data-dependent and the result is conditioned upon it, the reported σ value is statistically forced by the per-multiplet optimization rather than testing an independent, pre-specified hypothesis.

full rationale

The paper's central statistical claim (5.8 σ overall, 4.5 σ for Fornax A) is obtained by first identifying 28 multiplets in fixed windows, then for each multiplet selecting the single best association among 24 discrete backtracking combinations (3 GMF models × 8 species) to one of 10 radio galaxies, and finally reporting significance conditioned on those selections. This selection occurs after the data are observed and is not shown to be replicated identically in the Monte Carlo null distribution, so the quoted p-value is tied to the optimization step rather than a fully pre-specified test. The initial multiplet search itself uses pre-defined criteria and is not circular, but the load-bearing significance step reduces to a fitted-input procedure.

Axiom & Free-Parameter Ledger

4 free parameters · 2 axioms · 0 invented entities

The central claim rests on several pre-defined but arbitrary search parameters and domain assumptions about magnetic fields and particle types that are not derived within the work.

free parameters (4)
  • angular window = 3 degrees
    Pre-defined fixed 3-degree window chosen to isolate high-rigidity particles and keep chance coincidences low
  • time window = 15 days
    Pre-defined fixed 15-day window for the spatiotemporal multiplet search
  • energy threshold = 32 EeV
    Threshold above which the 28 multiplets are identified
  • target radio galaxies = 10 nearest
    Selection of the ten nearest bright radio galaxies for the association test
axioms (2)
  • domain assumption The three Galactic magnetic field models accurately describe UHECR deflections for the eight nuclear species tested
    Invoked when backtracking trajectories to test associations
  • domain assumption The kinematic filter isolates genuine high-rigidity particles without introducing selection bias
    Stated as the design goal of the search method

pith-pipeline@v0.9.0 · 5535 in / 1764 out tokens · 60853 ms · 2026-05-08T07:38:49.648071+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 16 canonical work pages · 1 internal anchor

  1. [1]

    The Astrophysics of Ultrahigh Energy Cosmic Rays

    K. Kotera and A.V. Olinto,The Astrophysics of Ultrahigh-Energy Cosmic Rays,Annual Reviews49(2011) 119 [1101.4256]

    work page Pith review arXiv 2011

  2. [2]

    R. Aloisio,Acceleration and propagation of ultra-high energy cosmic rays,Progress of Theoretical and Experimental Physics2017(2017) 12A102 [https://academic.oup.com/ptep/article-pdf/2017/12/12A102/22075738/ptx115.pdf]

    2017

  3. [3]

    Alves Batista, J

    R. Alves Batista, J. Biteau, M. Bustamante, K. Dolag et al.,Open questions in cosmic-ray research at ultrahigh energies,Frontiers in Astronomy and Space Sciences6(2019) 23

    2019

  4. [4]

    Jansson and G

    R. Jansson and G.R. Farrar,A New Model of the Galactic Magnetic Field,Astrophys. J.757 (2012) 14 [1204.3662]

    work page arXiv 2012

  5. [5]

    Pshirkov, P.G

    M.S. Pshirkov, P.G. Tinyakov, P.P. Kronberg and K.J. Newton-McGee,Deriving the Global Structure of the Galactic Magnetic Field from Faraday Rotation Measures of Extragalactic Sources,Astrophys. J.738(2011) 192 [1103.0814]. – 33 –

    work page arXiv 2011

  6. [6]

    Terral and K

    P. Terral and K. Ferri` ere,Constraints from faraday rotation on the magnetic field structure in the galactic halo,Astronomy & Astrophysics600(2017) A29

    2017

  7. [7]

    Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, I.A. Samarai, I.F.M. Albuquerque et al.,Observation of a large-scale anisotropy in the arrival directions of cosmic rays above 8× 1018 eV,Science357(2017) 1266 [1709.07321]

    work page arXiv 2017

  8. [8]

    A. Aab, P. Abreu, M. Aglietta, I.F.M. Albuquerque, I. Allekotte, A. Almela et al.,An Indication of Anisotropy in Arrival Directions of Ultra-high-energy Cosmic Rays through Comparison to the Flux Pattern of Extragalactic Gamma-Ray Sources,ApJL853(2018) L29 [1801.06160]

    work page arXiv 2018

  9. [9]

    Abreu, M

    P. Abreu, M. Aglietta, J.M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela et al., Arrival Directions of Cosmic Rays above 32 EeV from Phase One of the Pierre Auger Observatory,ApJ935(2022) 170 [2206.13492]

    work page arXiv 2022

  10. [10]

    Hardcastle and J.H

    M.J. Hardcastle and J.H. Croston,Radio galaxies and feedback from AGN jets,New Astronomy Reviews88(2020) 101539 [2003.06137]

    work page arXiv 2020

  11. [11]

    Matthews, A.R

    J.H. Matthews, A.R. Bell, K.M. Blundell and A.T. Araudo,Ultrahigh energy cosmic rays from shocks in the lobes of powerful radio galaxies,Monthly Notices of the Royal Astronomical Society482(2019) 4303

    2019

  12. [12]

    van Velzen, H

    S. van Velzen, H. Falcke, P. Schellart, N. Nierstenh¨ ofer and K.-H. Kampert,Radio galaxies of the local universe. all-sky catalog, luminosity functions, and clustering,Astronomy & Astrophysics544(2012) A18

    2012

  13. [13]

    Lang, A.M

    R.G. Lang, A.M. Taylor, M. Ahlers and V. de Souza,Revisiting the distance to the nearest ultrahigh energy cosmic ray source: Effects of extragalactic magnetic fields,Phys. Rev. D102 (2020) 063012

    2020

  14. [14]

    Abdul Halim, P

    A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela et al., Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory,Journal of Cosmology and Astroparticle Physics2024(2024) 022 [2305.16693]

    work page arXiv 2024

  15. [15]

    Abdul Halim, P

    A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela et al., Inference of the Mass Composition of Cosmic Rays with Energies from 1018.5 to 1020 eV Using the Pierre Auger Observatory and Deep Learning,Physical Review Letters134(2025) 021001 [2406.06315]

    work page arXiv 2025

  16. [16]

    E. Kido, T. Inakura, M. Kimura, N. Kobayashi, S. Nagataki, N. Shimizu et al.,Evaluations of uncertainties in simulations of propagation of ultrahigh-energy cosmic-ray nuclei derived from microscopic nuclear models,Astroparticle Physics152(2023) 102866 [2206.03447]

    work page arXiv 2023

  17. [17]

    E. Khan, S. Goriely, D. Allard, E. Parizot, T. Suomij¨ arvi, A.J. Koning et al., Photodisintegration of ultra-high-energy cosmic rays revisited,Astroparticle Physics23(2005) 191 [astro-ph/0412109]

    work page arXiv 2005

  18. [18]

    Abreu, M

    P. Abreu, M. Aglietta, E. Ahn, I. Albuquerque, D. Allard, I. Allekotte et al.,Search for signatures of magnetically-induced alignment in the arrival directions measured by the pierre auger observatory,Astroparticle Physics35(2012) 354

    2012

  19. [19]

    Abreu, M

    P. Abreu, M. Aglietta, J.M. Albury, I. Allekotte, K.A. Cheminant, A. Almela et al.,Arrival Directions of Cosmic Rays above 32 EeV from Phase One of the Pierre Auger Observatory, The Astrophysical Journal935(2022) 170

    2022

  20. [20]

    Batista, J.B

    R.A. Batista, J.B. Tjus, J. D¨ orner, A. Dundovic, B. Eichmann, A. Frie et al.,Crpropa 3.2 - an advanced framework for high-energy particle propagation in extragalactic and galactic spaces, Journal of Cosmology and Astroparticle Physics2022 09 035(2022) . – 34 –

    2022

  21. [21]

    Huchra, L.M

    J.P. Huchra, L.M. Macri, K.L. Masters, T.H. Jarrett, P. Berlind, M. Calkins et al.,The 2MASS Redshift Survey—Description and Data Release,ApJS199(2012) 26 [1108.0669]

    work page arXiv 2012

  22. [22]

    Virtanen, R

    P. Virtanen, R. Gommers, T.E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau et al., SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,Nature Methods17 (2020) 261

    2020

  23. [23]

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

    Astropy Collaboration, A.M. Price-Whelan, P.L. Lim, N. Earl, N. Starkman, L. Bradley et al., The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package,ApJ935(2022) 167 [2206.14220]. [24]Pierre Augercollaboration,The angular resolution of the Pierre Auger Observatory,Nucl. P...

    work page internal anchor Pith review arXiv 2022

  24. [24]

    Gilmore, R.S

    R.C. Gilmore, R.S. Somerville, J.R. Primack and A. Dom´ ınguez,Semi-analytic modelling of the extragalactic background light and consequences for extragalactic gamma-ray spectra,Monthly Notices of the Royal Astronomical Society422(2012) 3189 [1104.0671]

    work page arXiv 2012

  25. [25]

    Koning, S

    A.J. Koning, S. Hilaire and M.C. Duijvestijn,TALYS-1.0, inProceedings of the International Conference on Nuclear Data for Science and Technology, O. Bersillon, F. Gunsing, E. Bauge, R. Jacqmin and S. Leray, eds., (Nice, France), pp. 211–214, EDP Sciences, 2008

    2008

  26. [26]

    Bonifazi,The angular resolution of the pierre auger observatory,Nuclear Physics B - Proceedings Supplements190(2009) 20

    C. Bonifazi,The angular resolution of the pierre auger observatory,Nuclear Physics B - Proceedings Supplements190(2009) 20

    2009

  27. [27]

    Puget, F.W

    J.L. Puget, F.W. Stecker and J.H. Bredekamp,Photonuclear interactions of ultrahigh energy cosmic rays and their astrophysical consequences,The Astrophysical Journal205(1976) 638

    1976

  28. [28]

    Alves Batista, A

    R. Alves Batista, A. Dundovic, M. Erdmann, K.-H. Kampert, D. Kuempel, G. M¨ uller et al., CRPropa 3—a public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles,Journal of Cosmology and Astroparticle Physics2016(2016) 038 [1603.07142]

    work page arXiv 2016

  29. [29]

    Gerasimova and G.T

    N.M. Gerasimova and G.T. Zatsepin,Disintegration of Cosmic-Ray Nuclei by Solar Photons, Soviet Physics JETP11(1960) 899

    1960

  30. [30]

    Matthews, A.R

    J.H. Matthews, A.R. Bell, K.M. Blundell and A.T. Araudo,Fornax a, centaurus a, and other radio galaxies as sources of ultrahigh energy cosmic rays,Monthly Notices of the Royal Astronomical Society: Letters479(2018) L76

    2018

  31. [31]

    Harari, S

    D. Harari, S. Mollerach and E. Roulet,Lensing of ultrahigh-energy cosmic rays in turbulent magnetic fields,JHEP07(2002) 006 [astro-ph/0205484]. – 35 –

    work page arXiv 2002