Pith. sign in

REVIEW 4 major objections 4 minor 35 references

The Large Magellanic Cloud as a source of the highest energy cosmic rays

T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that the highest-energy cosmic rays recorded by the Pierre Auger Observatory are clustered on the sky at the 1% significance level, in a direction close to the Large Magellanic Cloud, which it proposes as a possible source.

desk verdict A transparent but statistically under-powered speculation: the LMC source idea is fresh, but the 1% clustering claim is an uncorrected scan minimum and won't survive trial correction. read the letter →

arxiv 2507.16378 v1 pith:MDPSNDUV submitted 2025-07-22 astro-ph.HE

Roger Clay , Tadeusz Wibig This is my paper · ORCID
classification astro-ph.HE
keywords ultra-high-energycosmicraysPierreAugerObservatoryLargeMagellanicCloudCentaurusAarrival-directionanisotropyfactorialmomentsHillascriterioncosmic-raysources
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper asks whether the most energetic cosmic rays ever recorded, the top 100 events from the Pierre Auger Observatory above 100 EeV, carry any directional memory of where they were born. It finds that the most energetic events are not uniformly scattered: the five most energetic cluster in a patch of sky about 30 degrees across, with a random-chance probability of about 1%. That patch lies closer, in angular terms, to the Large Magellanic Cloud than to Centaurus A, the active galaxy usually discussed for this energy range. The paper argues that the LMC, despite being a dwarf galaxy, has the magnetic-field structures and proximity needed to accelerate particles to these energies. If correct, the result would mean ultra-high-energy cosmic-ray astronomy is possible and that nearby dwarf galaxies, not just active galactic nuclei, deserve serious consideration as sources.

What carries the argument

The carrying object is the factorial cumulant $F_k(\Delta\varphi)$, defined in equation (1), which counts how many groups of $k$ arrival directions have all mutual angular separations smaller than a cluster diameter $\Delta\varphi$, normalized by the number expected from a randomized sky. For orders $k=2$ and $k=3$, this is the statistic whose minimum over scanned angular diameters yields the reported 1% chance probability. The random background is built by randomizing right ascension while preserving the observatory's declination-dependent acceptance, and the same machinery is used to test robustness by varying the number of selected most-energetic events from 3 to 14 and by dropping the top five events.

What would settle it

Recompute the factorial-moment probabilities on the published Auger catalog with a full trial-factor correction: for every group size from 3 to 14 and every angular diameter, record the minimum chance probability, then compare that minimum against the distribution of minima from many randomized realizations. If the corrected global probability rises above about 10%, the clustering is consistent with chance and the LMC source hypothesis loses its basis. A second, independent falsifier is the next data release: if additional highest-energy events scatter outside the roughly 30-degree patch, the observed cluster is a statistical fluctuation.

Watch

Extended reading notes

Core claim

Using factorial cumulant moments of orders two through four computed from angular separations among the Auger events, the paper finds that the five most energetic showers deviate from a randomized sky at about the 1% probability level. Three of the five most energetic events, and four of the ten highest, lie within only about 30 degrees of one another. When the five most energetic events are replaced by the next five, the clustering disappears, showing the effect is driven by the extreme tail. The paper then notes that the cluster direction is closest to the Large Magellanic Cloud, whose 30 Doradus region has measured magnetic fields above 300 microgauss over tens of parsecs, satisfying the Hillas confinement criterion for iron nuclei at these rigidities. Back-tracking through Milky Way magnetic-field models gives deflections of about 15 degrees in that sky region, enough to connect the observed cluster to the LMC.

Load-bearing premise

The reported 1% significance depends on the assumption that the randomized sky used as the null hypothesis faithfully represents chance and that scanning over many group sizes and angular separations does not inflate the probability.

Editorial extensions

If this is right

  • If the 1% clustering is real, ultra-high-energy cosmic rays retain directional information at around 100 EeV, opening the possibility of charged-particle astronomy with existing observatories.
  • A source direction closer to the LMC than to Centaurus A would shift searches from active galactic nuclei to a nearby dwarf galaxy and its star-forming region 30 Doradus.
  • Because the LMC is only about 50 kiloparsecs away, particles from it would suffer essentially no GZK attenuation, naturally explaining why the most extreme tail of the energy spectrum may be dominated by a very local accelerator.
  • Future Auger data releases, including post-upgrade observations, can test the claim directly: additional highest-energy events should continue to fall in the same roughly 30-degree patch if the clustering is real.
  • The absence of clustering among events 6 through 10 implies the anisotropy lives in the extreme-energy tail, so larger future exposures should sharpen the signal rather than dilute it.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the clustering survives a proper trial-factor correction, the same factorial-moment scan applied to the Northern-hemisphere Telescope Array catalog could reveal whether the anisotropy is global or a southern-sky phenomenon, helping distinguish a physical source from an observatory-related effect.
  • The paper's 15-degree deflection estimate is rigidity dependent, so a testable extension would be to check whether the angular spread of the cluster shrinks for events with higher magnetic rigidity, as expected if the deflections are caused by the Galactic magnetic field.
  • The LMC hypothesis predicts that the excess direction should be spatially correlated with the 30 Doradus region and that its energy spectrum should continue to higher energies with little attenuation, offering a concrete multi-messenger test with gamma-ray and neutrino observatories.
  • The main unresolved question is statistical: recomputing the clustering probability while accounting for the scan over group sizes and angular diameters would settle whether the claimed 1% significance survives or dissolves into a larger chance probability.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. This paper analyzes the published Pierre Auger catalog of the 100 highest-energy cosmic-ray events recorded over 17 years. Using factorial-cumulant statistics, it reports that the five most energetic events are clustered at a chance probability of order 1%, with a minimum around 0.5% when scanning over group sizes and angular diameters, and notes that four of the ten highest-energy events lie within about 30 degrees of one another. Since the Large Magellanic Cloud is angularly close to the cluster and closer in angle than Centaurus A, the paper examines the LMC's viability as a source: its magnetic fields and Hillas criterion for heavy nuclei, its proximity, and its lack of GZK attenuation. The paper concludes that the LMC may be a source of these highest-energy events.

Significance. If the clustering is real, the paper points to a potentially important anisotropy in the highest-energy Auger events and identifies a nearby, relatively underexplored source candidate. Its strengths are that it uses public Auger data, makes an explicit a-priori choice for the top-five sample, tests sensitivity by discarding those events and by varying the selected number of events, and appropriately frames the LMC suggestion as an examination rather than a definitive detection. The analysis is not circular: the cluster is found first and the LMC is then considered as a candidate. However, the central statistical claim is not yet established because the reported probability is an uncorrected minimum over a multi-dimensional scan, so the manuscript's main conclusion is conditional on a global trial-correction test.

major comments (4)
  1. [Section II.A and Fig. 3] The p-values quoted as evidence for clustering are minima over angular diameter for each n and over n = 3,...,14. The text explicitly says these values 'should be taken with care', yet the abstract and Section II state the result as non-random at the 1% level. A Bonferroni-like correction with the number of scanned diameters and subset sizes would move p_min ≈ 0.005 to values of order 0.1 or larger, so the global chance probability may be fully consistent with isotropy. The authors should report a single global Monte Carlo p-value that includes all scans (diameters, n, and any post-hoc choice to highlight the cluster), or remove the 1% claim from the abstract.
  2. [Section II.A (null hypothesis)] The randomization of the null hypothesis is described in one sentence: 'randomize the angle of Right Ascension and the choice of declination, which is defined by the acceptance of a certain apparatus for a specific energy range.' It is not specified whether this reproduces the full Auger exposure map in both right ascension and declination, nor how the energies of the selected events enter the analysis. If the null is not the actual Auger exposure, the chance probabilities are biased. The authors should specify or implement an exposure-aware isotropic Monte Carlo.
  3. [Section III and IV (source assignment)] The assignment of the cluster to the LMC is post hoc. The paper does not enumerate candidate astrophysical objects within the approximately 15-degree uncertainty region (Cen A, the LMC, and other nearby galaxies) and does not account for the fact that the LMC was chosen after viewing the sky map and after the known Cen A association. A corrected significance for the specific LMC association, for example via a catalog scan or an independent dataset, is needed before the 'closer in angular terms' argument can carry the source claim.
  4. [Section III (magnetic deflection)] The 15-degree backtracking estimate from Unger and Farrar is used to favor the LMC over Cen A, but the authors themselves note that the modeled deflection directions are inconsistent and that the real deflections could be much lower. Because the angular separation between the cluster, the LMC, and Cen A is comparable to this deflection scale, the source attribution depends strongly on an uncertain Galactic magnetic-field model. The paper should propagate this uncertainty into the final claim, for example by presenting the cluster-source comparison under both zero-deflection and modeled-deflection assumptions.
minor comments (4)
  1. [Abstract and Section IV] The abstract says 'non-random at the 1% level' while Section IV says the clustering is 'also compatible with a random expectation at the 1% level'; these phrasings should be reconciled to avoid implying a discovery-level significance.
  2. [Fig. 2 and Eq. (1)] The Fig. 2 caption should explain how the shaded chance-probability bands are computed, including the simulation setup, the number of trials, and whether the Auger exposure is modeled; the 'norm.' factor in Eq. (1) should also be defined explicitly.
  3. [References] Reference [12] appears unrelated to the Milky Way magnetic-field statement in the Introduction and should be replaced or justified.
  4. [General] There are several small typographical issues: 'Pbservatory' in reference [20], 'SOFIA' capitalization, and the URL in reference [29] is for a news item rather than a citable data release.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity; the clustering is an empirical scan of published Auger data and the LMC hypothesis is a post-hoc candidate comparison using independent external inputs.

full rationale

The paper's central claim is an observed directional clustering among the five (and ten) highest-energy Auger events, quantified by factorial moments against an isotropically randomized background. This is not a derivation of a result from an input that already contains it: the cluster locations are taken from the public Auger catalog, the null distribution is generated by randomizing right ascension and sampling declinations from the detector acceptance, and no parameter is fitted from the cluster and then renamed as a prediction. The choice of the LMC as a candidate source is made after inspecting the sky map because it lies close in angle to the cluster; that is post-hoc selection, which is a statistical/interpretation concern, not circular reasoning. The plausibility arguments use independent external data: SOFIA magnetic-field measurements of 30 Doradus, Auger composition measurements, Hillas criteria, and the Unger-Farrar galactic magnetic-field backtracking. Cited prior work by one of the authors ([10], on the universal cosmic-ray spectrum and composition) is peripheral to the clustering analysis and is not load-bearing. The concern that the p_min(n) values in Fig. 3 are minima over diameter and event-count scans and may need a trial-factor correction is a substantive statistical criticism of the reported 1% significance, but it does not make the argument circular; the paper itself cautions that the probability values 'should be taken with care.' No circular step can be exhibited by reduction of an equation to its own input, so the appropriate finding is no significant circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new entities and fits no free parameters. Its conclusions rest on the public Auger catalog, on a randomized-background null model, on the Hillas criterion with measured magnetic fields, and on the assumed Galactic magnetic field deflection model.

assumptions (5)
  • standard math Factorial moments Fk as defined in Eq. (1) are a valid statistic for angular clustering, and the normalization 'norm.' correctly estimates the expected number of random k-tuples.
    Invoked in Section II to define the clustering measure and the comparison with random expectations.
  • domain assumption Randomizing Right Ascension and declination according to the Auger acceptance produces an unbiased uncorrelated background distribution.
    Section II.A: 'randomize the angle of Right Ascension and the choice of declination, which is defined by the acceptance'.
  • domain assumption The Galactic magnetic field backtracking map of Unger and Farrar (2025) gives deflections of order 15 degrees in the cluster region, sufficient to cover the LMC direction but not necessarily Centaurus A.
    Section III, citing [35]; this is the link between observed directions and the proposed LMC source.
  • domain assumption At the highest energies the cosmic ray composition is dominated by heavy nuclei.
    Section III, citing [33]; required for 30 Doradus to satisfy the Hillas criterion for iron rather than protons.
  • domain assumption The Hillas criterion and turbulent acceleration in magnetically dominated turbulence are sufficient to establish that 30 Doradus could in principle accelerate particles to 100+ EeV.
    Section III, citing [32] and [15]; the SOFIA measured B-field and size are used in a rule-of-thumb estimate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Large Magellanic Cloud as a source of the highest energy cosmic rays." pith.science (2026). https://pith.science/paper/MDPSNDUV

@misc{pith2026250716378,
  author       = {Pith},
  title        = {Pith review of: The Large Magellanic Cloud as a source of the highest energy cosmic rays},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MDPSNDUV}},
  note         = {Machine review of arXiv:2507.16378}
}
read the original abstract

The Pierre Auger Observatory has published properties of the 100 highest-energy cosmic ray events (to energies above 100 EeV) which it recorded over a 17 year period. We have examined the directional properties of these events and have taken particular note of the most energetically extreme events. We find that the most energetic events have directions which are grouped in a non-random way at the 1\% level. There is an apparent clustering in a limited region of the sky. Close to that direction is found Centaurus A, which has long been considered as a source of such particles, but we also note that a close-by dwarf galaxy, the Large Magellanic Cloud (LMC) is closer in angular terms. We examine the possibility that the LMC might be a source of observed cosmic rays at the highest energies.

Figures

Figures reproduced from arXiv: 2507.16378 by the authors.

Figure 1
Figure 1. FIG. 1. Sky map in celestial co-ordinates showing direction [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Factorial cumulants or second (a), third (b), fourth [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The probability that, among Auger [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 33 canonical work pages

  1. [1]

    K. W. Wo´ zniak (for the CREDO Collaboration), Detection of cosmic-ray ensembles with CREDO, EPJ Web Conf. 208, 15006 (2019)

    work page 2019

  2. [2]

    Abramowski et al

    A. Abramowski et al. (H.E.S.S. Collaboration), Discove ry of γ-ray emission from the extra- galactic pulsar wind nebula N 157B with H.E.S.S., A&A 545, L2 (2012). 10

    work page 2012

  3. [3]

    Abramowski et al

    A. Abramowski et al. (H.E.S.S. Collaboration), The exce ptionally powerful TeV γ-ray emit- ters in the Large Magellanic Cloud, Science 347, 406 (2015)

    work page 2015

  4. [4]

    A. A. Halim et al. (Pierre Auger Collaboration), Constra ining models for the origin of ultra- high-energy cosmic rays with a novel combined analysis of ar rival directions, spectrum, and composition data measured at the Pierre Auger Observatory, Journal of Cosmology and As- troparticle Physics 2024 (01), 022

    work page 2024

  5. [5]

    Aab et al

    A. Aab et al. (Pierre Auger Collaboration), The Pierre Au ger Cosmic Ray Observatory, Nu- clear Instruments and Methods in Physics Research A 798, 172 (2015)

    work page 2015

  6. [6]

    A. A. Halim et al. (Pierre Auger Collaboration), The Pier re Auger Observatory open data, Eur. Phys. J. C 85, 70 (2025)

    work page 2025

  7. [7]

    Coleman et al

    A. Coleman et al. (TA Collaboration), Ultra high energy c osmic rays the intersection of the cosmic and energy frontiers, Astroparticle Physics 149, 102819 (2023)

    work page 2023

  8. [8]

    Abu-Zayyad, T., et al., The surface detector array of the telescope array experiment, Nuclear Instruments and Methods in Physics Research A 689, 87 (2012)

    work page 2012

Show all 35 references
  1. [9]

    et al., New air fluorescence detectors employe d in the telescope array experiment, Nuclear Instruments and Methods in Physics Research A 676, 54 (2012)

    Tokuno, H. et al., New air fluorescence detectors employe d in the telescope array experiment, Nuclear Instruments and Methods in Physics Research A 676, 54 (2012)

    2012

  2. [10]

    Wibig and A

    T. Wibig and A. W. Wolfendale, Universal cosmic rays ene rgy spectrum and the mass com- position at the ’ankle’ and above (2017), arXiv:1706.05872 [astro-ph.HE]

    2017 arXiv

  3. [11]

    S. A. Mao et al., Magnetic field structure of the Large Mag ellanic Cloud from Faraday rota- tion measures of diffuse polarized emission, The Astrophysic al Journal 759, 25 (2012)

    2012

  4. [12]

    Unger and G

    M. Unger and G. Farrar, (In)Feasability of Studying Ult ra-High-Energy Cosmic Rays with Smartphones (2015)

    2015

  5. [13]

    Fermi, On the origin of the cosmic radiation, Phys

    E. Fermi, On the origin of the cosmic radiation, Phys. Re v. 75, 1169 (1949)

    1949

  6. [14]

    A. M. Hillas, Can diffusive shock acceleration in superno va remnants account for high-energy galactic cosmic rays?, Journal of Physics G: Nuclear and Par ticle Physics 31, R95 (2005)

    2005

  7. [15]

    Comisso, G

    L. Comisso, G. R. Farrar, and M. S. Muzio, Ultra-high-en ergy cosmic rays accelerated by magnetically dominated turbulence, The Astrophysical Jou rnal Letters 977, L18 (2024)

    2024

  8. [16]

    Greisen, End to the cosmic-ray spectrum?, Phys

    K. Greisen, End to the cosmic-ray spectrum?, Phys. Rev. Lett. 16, 748 (1966)

    1966

  9. [17]

    G. T. Zatsepin and V. A. Kuz’min, Upper Limit of the Spect rum of Cosmic Rays, ZhETF Pisma Redaktsiiu 4, 114 (1966). 11

    1966

  10. [18]

    A. A. Halim et al. (The Pierre Auger Collaboration), A ca talog of the highest-energy cos- mic rays recorded during Phase I of operation of the Pierre Au ger Observatory (2022), arXiv:2211.16020 [astro-ph.HE]

    2022 arXiv

  11. [19]

    Abreu et al

    P. Abreu et al. (Pierre Auger Collaboration), Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter, Astropartic le Physics 34, 314 (2010)

    2010

  12. [20]

    A. A. Halim et al. (The Pierre Auger Collaboration), The flux of ultra-high-energy cos- mic rays along the supergalactic plane measured at the Pierr e Auger Pbservatory (2024), arXiv:2407.06874 [astro-ph.HE]

    2024

  13. [21]

    Alasalı and Y

    T. Alasalı and Y. Ortakcı, Clustering techniques in dat a mining: A survey of methods chal- lenges, and applications, Computer Science 9, 32 (2024)

    2024

  14. [22]

    Bialas and R

    A. Bialas and R. Peschanski, Moments of rapidity distri butions as a measure of short-range fluctuations in high-energy collisions, Nuclear Physics B 273, 703 (1986)

    1986

  15. [23]

    Bialas and R

    A. Bialas and R. Peschanski, Intermittency in multipar ticle production at high energy, Nuclear Physics B 308, 857 (1988)

    1988

  16. [24]

    P. Lipa, P. Carruthers, H. Eggers, and B. Buschbeck, The correlation integral as probe of multiparticle correlations, Physics Letters B 285, 300 (1992)

    1992

  17. [25]

    De Wolf, I

    E. De Wolf, I. Dremin, and W. Kittel, Scaling laws for den sity correlations and fluctuations in multiparticle dynamics, Physics Reports 270, 1 (1996)

    1996

  18. [26]

    Kittel, Correlations and fluctuations in high-energ y collisions, in Particle Production Span- ning MeV and TeV Energies (Springer Netherlands, Dordrecht, 2000) Chap

    W. Kittel, Correlations and fluctuations in high-energ y collisions, in Particle Production Span- ning MeV and TeV Energies (Springer Netherlands, Dordrecht, 2000) Chap. 6, pp. 157–1 82

    2000

  19. [27]

    A. A. Abdo et al., Observations of the Large Magellanic C loud with Fermi, A&A 512, A7 (2010)

    2010

  20. [28]

    and et al., The SRG/eROSITA all-sky survey - first x-ray catalogues and data release of the western galactic hemisphere, A&A 682, A34 (2024)

    Merloni, A. and et al., The SRG/eROSITA all-sky survey - first x-ray catalogues and data release of the western galactic hemisphere, A&A 682, A34 (2024)

    2024

  21. [29]

    Sky & Telescope, The first all-sky X-ray map to be release d in 30 years reveals new wonders of the hot and energetic universe, https://skyandtelescop e.org/astronomy-news/first-all-sky- map-erosita (2020)

    2020

  22. [30]

    J. W. Cronin, Cosmic rays: the most energetic particles in the universe, Rev. Mod. Phys. 71, S165 (1999)

    1999

  23. [31]

    L. N. Tram, et al., SOFIA observations of 30 Doradus. II. Magnetic fields and large-scale gas kinematics, The Astrophysical Journal 946, 8 (2023). 12

    2023

  24. [32]

    A. M. Hillas, The origin of ultra-high-energy cosmic ra ys, Annual Review of Astronomy and Astrophysics 22, 425 (1984)

    1984

  25. [33]

    A. A. Halim et al. (The Pierre Auger Collaboration), Inf erence of the mass composition of cosmic rays with energies from 10 18.5 to 10 20 eV using the Pierre Auger Observatory and deep learning, Phys. Rev. Lett. 134, 021001 (2025)

    2025

  26. [34]

    A. A. Halim et al. (The Pierre Auger Collaboration), Lar ge-scale Cosmic-ray Anisotropies with 19 yr of Data from the Pierre Auger Observatory, The Astr ophysical Journal 976, 48 (2024), arXiv:2408.05292 [astro-ph.HE]

    2024 arXiv

  27. [35]

    Unger and G

    M. Unger and G. R. Farrar, The Galactic magnetic field and UHECR deflections (2025), arXiv:2502.15876 [astro-ph.HE]. 13

    2025 arXiv

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.