Pith. sign in

REVIEW 5 major objections 5 minor 1 cited by

Hint at an axion-like particle from GRB 221009A

T0 review · 5 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read The paper argues that conventional physics cannot explain LHAASO's detection of GRB 221009A photons above 10 TeV, and that photon-axion-like-particle oscillations with ALP mass around 10^-10 eV and coupling around 4x10^-12 GeV^-1 resolve…

desk verdict An update of the authors' own ALP explanation for GRB 221009A, but the central 'should not have observed' claim is not actually established because no expected-count calculation is shown. read the letter →

arxiv 2412.21175 v1 pith:5HHEWO4H submitted 2024-12-30 astro-ph.HE hep-ph

Giorgio Galanti , Lara Nava , Marco Roncadelli , Fabrizio Tavecchio , Giacomo Bonnoli This is my paper · ORCID
classification astro-ph.HEhep-ph
keywords gamma-rayburstsGRB221009Aaxion-likeparticlesphoton-ALPoscillationsextragalacticbackgroundlightveryhighenergygammaraysLHAASO
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

The paper sets out to show that the LHAASO detection of GRB 221009A at energies between 13 and 18 TeV is incompatible with the current best model of the extragalactic background light, which predicts survival probabilities of about $3x10^{-6}$ at 15 TeV and $1x10^{-8}$ at 18 TeV. It then argues that photon-axion-like-particle oscillations in the magnetic fields of the host galaxy, intergalactic space, and the Milky Way restore the observed flux, with ALP mass around $10^{-10}$ eV and coupling around $4x10^{-12}$ $GeV^{-1}$. If the argument holds, the Universe is more transparent to very-high-energy gamma rays than conventionally believed, and GRB 221009A becomes the strongest current hint of axion-like particles.

What carries the argument

The key object is the Lagrangian coupling Laγγ = gaγγ E·B a, which lets a photon oscillate into an axion-like particle in the presence of an external magnetic field. The argument works by computing the photon-ALP transfer matrix through each magnetized medium — the GRB host galaxy (modelled as a spiral or a starburst), extragalactic space (with field strength 1 nG or below $10^{-15}$ G), and the Milky Way — and combining these with the EBL absorption baseline of the adopted model to obtain the final energy-dependent survival probability PALP. QED vacuum polarization and photon dispersion on the cosmic microwave background are included. This quantity, evaluated at the benchmark ALP parameters, is what reconciles the observed >10 TeV flux with the LHAASO spectrum.

What would settle it

A measurement that would settle the claim is an independent determination of the extragalactic background light that gives a survival probability at 15 TeV larger than about $3x10^{-6}$, which would remove the crisis; or, on the positive side, the detection of the predicted oscillatory spectral pattern above 10 TeV in GRB 221009A-like events, which would support the ALP interpretation.

Watch

Extended reading notes

Core claim

On its own terms, the paper's central discovery is a quantitative contradiction: under the state-of-the-art extragalactic background light model used in the paper, conventional physics gives essentially zero chance of observing a 15-18 TeV photon from a source at z=0.151, yet the LHAASO Collaboration saw such photons from GRB 221009A. The paper finds that including photon-ALP oscillations lifts the photon survival probability enough to match the observed spectrum using ALP parameters ma ~ $10^{-10}$ eV and gaγγ ~ $4x10^{-12}$ $GeV^{-1}$, values that respect the most stringent current bounds. It therefore claims that the GRB provides the strongest hint at ALP existence to date, reinforcing two earlier indications from blazars.

Load-bearing premise

The load-bearing premise is that the state-of-the-art EBL model used in the paper correctly describes how strongly very-high-energy photons are absorbed; if the true extragalactic background light attenuates photons less than this model says, the discrepancy shrinks and the ALP explanation is no longer necessary.

Editorial extensions

If this is right

  • If correct, the standard EBL absorption prediction fails for GRB 221009A, meaning the Universe at z ≈ 0.15 is more transparent to >10 TeV photons than the adopted EBL model implies.
  • The benchmark ALP mass and coupling lie inside current experimental and astrophysical bounds, so the proposed solution does not require new parameter space beyond what is already allowed.
  • The same ALP parameters connect GRB 221009A to two earlier blazar hints, making a consistent ALP interpretation across three independent source classes.
  • Among the new-physics options considered, Lorentz invariance violation cannot explain the detection under current limits, leaving photon-ALP oscillations as the only workable explanation in the paper's comparison.

Reading between the lines

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

  • A natural next test is to look for the same enhanced transparency in a population of high-redshift GRBs: if several >10 TeV detections accumulate, the ALP interpretation would gain statistical weight, whereas a single conventional cutoff would strain it.
  • Because the entire crisis rests on the adopted EBL model, an independent EBL determination with lower optical depth near 10-20 TeV would dissolve the tension without invoking new particles; the ALP claim should therefore be read as conditional on the EBL model's accuracy.
  • The same energy-dependent survival probability predicts spectral irregularities (oscillatory features) in the >10 TeV band; future very-high-energy gamma-ray observatories with fine spectral bins could look for these wiggles as a distinct ALP signature.
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

5 major / 5 minor

Summary. This manuscript argues that the LHAASO detection of GRB 221009A at z=0.151 with photons above 10 TeV is incompatible with the Saldana-Lopez et al. extragalactic background light (EBL) model, and that photon-ALP oscillations with ma ~ 1e-10 eV and gaγγ ~ 4e-12 GeV^-1 reconcile the observation. The authors quote EBL survival probabilities P_CP ~ 3e-6 at 15 TeV and ~1e-8 at 18 TeV, assert that the EBL-deabsorbed spectrum has no cutoff up to 7 TeV and can be extended to ~20 TeV, and then claim that conventional physics cannot explain the detection, whereas an ALP model from their prior PRL (ref 5) can. The paper reports survival probabilities with ALPs in a figure and concludes that GRB 221009A provides the strongest hint at ALP existence to date. Most quantitative details are delegated to refs 5 and 6.

Significance. If the central claim is correct, the paper would provide a significant astrophysical hint for ALPs with parameters compatible with dark matter candidates, and would imply that EBL attenuation at very high energies is overestimated by the Saldana-Lopez model. The manuscript builds on a state-of-the-art EBL model and a published ALP oscillation framework, and explicitly compares with existing ALP bounds, which are strengths. However, the core "should not have observed" claim is not yet quantitatively established because no expected-count calculation is presented, and the ALP parameter choice is partly optimized. The conclusion is therefore conditional and needs additional statistical support before the "strongest hint" statement can be sustained.

major comments (5)
  1. [Section 2 (right panel of Fig. 1) and Abstract] The claim that LHAASO "should not have observed" photons above 10 TeV is not supported by the survival probability alone. The paper quotes P_CP ≈ 3e-6 at 15 TeV and ≈ 1e-8 at 18 TeV, but never computes the expected number of detected events, which requires the emitted spectrum dN/dE, the LHAASO effective area A_eff(E), and the exposure T: N = T ∫ A_eff(E) (dN/dE) P(E) dE for E > 10 TeV. Without this calculation, a small survival probability can still be consistent with a detection if the intrinsic flux is large enough or if the spectrum hardens above 7 TeV. The statement that the EBL-deabsorbed LHAASO spectrum "shows no cutoff" up to 7 TeV does not logically assure an unbroken power-law extension to ~20 TeV. This missing statistical comparison is load-bearing because the entire ALP motivation depends on the alleged impossibility of the conventional detection.
  2. [Section 2, first paragraph] The tension with conventional physics is defined solely by the Saldana-Lopez (SL) EBL model, but the paper generalizes to "the Universe should be more transparent than currently believed" without testing robustness to EBL model uncertainties. The SL model is one of several recent EBL reconstructions, and survival probabilities vary among models; the authors should show that a representative set of current EBL models yields a similar suppression, or quantify the systematic spread. As written, the claimed problem is conditional on SL, whereas the ALP explanation is presented as an unconditional resolution. A robustness check against at least one or two other modern EBL models would strengthen the central argument considerably.
  3. [Section 2, Fig. 1 left and central panels] The benchmark ALP parameters (ma = 1e-10 eV, gaγγ = 4e-12 GeV^-1) are chosen both to maximize PALP at 15 TeV and to match two previous ALP hints (refs 15 and 17). The paper does not provide a likelihood or chi-square comparison between the conventional model and the ALP model against the observed LHAASO spectrum, nor does it show the sensitivity of the conclusion to variations of ma and gaγγ within the allowed band. Consequently, the statement that the ALP model succeeds "without the need of contrived choices of parameter values" is not quantitatively demonstrated. I recommend adding a parameter scan over the allowed region and a model-selection statistic (e.g., delta-chi-square or AIC) to support the claimed hint.
  4. [Section 2, host galaxy and extragalactic magnetic field assumptions] The robustness claim "we observe no substantial modification of our results in the two cases concerning the value assumed by Bext" is not supported by any figure, table, or quantitative statement in this manuscript. Likewise, the choice of a starburst host galaxy is asserted without a comparison to the spiral-host case. Since the ALP survival probability depends critically on the magnetic field configurations along the line of sight, these alternatives need to be shown or a specific reference to where they are computed must be supplied.
  5. [Section 2, method and Fig. 1] The central quantitative results—P_CP and PALP as functions of energy—are quoted from refs 5 and 6 without showing the transfer-matrix equations or the input parameters (magnetic field profiles, domain lengths, electron densities, ALP parameters) used to compute the PALP curves in Fig. 1. This prevents the reader from verifying the claimed factors, especially the ALP-induced enhancement above 10 TeV. If the manuscript is intended as a companion letter, it should state this explicitly and list the exact external parameters; otherwise an appendix with the calculation is needed.
minor comments (5)
  1. [Section 3, final paragraph] The text says the ALP parameters are within the most stringent bounds (refs 34 and 35), but Section 2 lists CAST (ref 33) as one of the most reliable bounds. Please clarify why CAST is omitted from the final statement or include it.
  2. [Abstract and Section 1] The phrase "the Universe should be more transparent than currently believed" is not a precise scientific statement; consider rephrasing to "the Saldana-Lopez EBL model predicts a lower transparency than inferred from the detection."
  3. [Section 2, benchmark parameters] Because the benchmark ALP parameters are taken from the authors' own previous hints (refs 15 and 17), the reader should be told explicitly whether the GRB analysis is a postdiction or a prediction with parameters fixed before the LHAASO data.
  4. [Section 2, right panel of Fig. 1] The caption should describe the LHAASO sensitivity curve and the spectral deabsorption procedure more completely, including the energy binning and the assumed intrinsic spectral model, so that the comparison in the right panel is reproducible.
  5. [Section 3, last sentence] The reference to "recent data released by the LHAASO Collaboration concerning the GRB 221009A spectrum above 10 TeV (ref 42)" is vague; specify which data set is used and how it differs from the earlier Science paper (ref 36).

Circularity Check

2 steps flagged · score 5.0 of 10

Partial circularity: ALP benchmark parameters are optimized at the LHAASO event energy and anchored to the authors' own previous hints, though the GRB data and EBL model are external.

  1. fitted input called prediction [Section 2, Figure 1 discussion (parameter choice and central panel)]
    "The left panel of Fig. 1 suggests us to take ma ≃ (10−11 − 10−7) eV and gaγγ ≃ (3 − 5) × 10−12 GeV−1 in order to maximize PALP around E = 15 TeV within all current most stringent ALP bounds."

    The benchmark parameters are explicitly varied to maximize the photon survival probability at E = 15 TeV, which the same section identifies as the 'benchmark energy for the LHAASO event'. The resulting enhanced P_ALP at 15 TeV is then presented as explaining the LHAASO observation above 10 TeV. Thus the target quantity, the survival probability at the observed highest-energy event, is optimized with respect to the model parameters before being used as evidence. The demonstration is therefore partly fitted to the datapoint it claims to explain.

  2. self citation load bearing [Section 2 and Section 3 (benchmark values and 'two other indications')]
    "we assume ma = 10−10 eV and gaγγ = 4 × 10−12 GeV−1 as benchmark values according to two previous indications at the ALP existence15,17 ... This fact regarding GRB 221009A represents a strong hint at the ALP existence, which adds to two other indications coming from blazars."

    The two 'previous indications' used to set the benchmark parameters are refs 15 and 17, whose author lists overlap with the present paper (Tavecchio, Roncadelli, Galanti, Bonnoli, De Angelis, Bignami). The conclusion that GRB 221009A 'adds to two other indications' therefore leans on the same prior self-citations that supply the parameter values. This is not fully load-bearing for the independent GRB calculation, because the ALP formalism and the LHAASO/SL data are external, but it makes the 'strongest hint to date' framing partly self-referential.

full rationale

The derivation chain is not globally circular: the LHAASO detection, the SL EBL attenuation, and the standard photon-ALP oscillation formalism are all independent inputs. The central comparison P_ALP versus P_CP is not an identity, and the model produces a full spectrum rather than a single fitted point. However, the paper's own text shows that the benchmark ALP parameters are selected in order to maximize P_ALP at E = 15 TeV, the energy of the LHAASO event the model is then used to explain; the survival-probability 'prediction' is therefore partly optimized against the target. The benchmark values are also justified by the authors' own previous ALP hints (refs 15,17), so the 'strongest hint to date' framing is partly self-referential. The absence of an expected-photon-count calculation for the 'should not have observed' claim is a substantive correctness/completeness issue, but it is not circularity because no fitted parameter is recycled into that claim. On balance, the independent external data and physics prevent a high score, but the parameter optimization and self-citation chain keep the paper from being fully self-contained.

Assumptions & free parameters 2 free parameters · 4 assumptions · 1 invented entities

The central ALP explanation rests on two chosen parameters (ma, gaγγ) taken from previous hints, an assumed EBL model, and an assumed intrinsic source spectrum. No direct detection of ALPs exists, making the evidence indirect.

free parameters (2)
  • ALP mass ma = 10^-10 eV (benchmark), range 10^-11 to 10^-7 eV
    Chosen to maximize photon survival probability around 15 TeV and aligned with previous ALP hints (refs 15,17). Not fitted to the GRB data directly, but selected from prior analyses and to match the observed energy.
  • ALP-photon coupling gaγγ = 4x10^-12 GeV^-1 (benchmark), range 3 to 5x10^-12 GeV^-1
    Selected from previous ALP hints and within the CAST bound; the range maximizes the photon survival probability at 15 TeV.
assumptions (4)
  • domain assumption Axion-like particles exist and interact with photons via the Lagrangian in Eq. (1).
    The entire solution assumes ALPs are real and that the coupling gaγγ has the specified strength. This is the hypothesis under test, not an established fact.
  • domain assumption The Saldana-Lopez et al. EBL model accurately represents the extragalactic background light.
    The conventional-physics survival probabilities (P_CP) are computed using this model; if it overestimates absorption, the tension with conventional physics disappears.
  • domain assumption The intrinsic gamma-ray spectrum of GRB 221009A continues as a power law without a cutoff between 7 TeV and 18 TeV.
    The paper extends the LHAASO-measured spectrum (no cutoff up to 7 TeV) to higher energies, which is necessary to claim a detection anomaly.
  • domain assumption The magnetic field models for the host galaxy, intergalactic medium, and Milky Way (refs 24,25,26,27,31) are adequate.
    The ALP conversion probability depends on these field configurations, and the paper adopts literature models without independent verification.
invented entities (1)
  • Axion-like particle (ALP) independent evidence
    purpose: To avoid EBL absorption by converting photons into non-interacting ALPs in magnetic fields and back into photons near Earth.
    ALPs are a theoretical construct with experimental bounds (e.g., CAST) but no direct detection. The paper's predicted mass and coupling range provides a falsifiable target for future direct searches, so there is an independent handle.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hint at an axion-like particle from GRB 221009A." pith.science (2026). https://pith.science/paper/5HHEWO4H

@misc{pith2026241221175,
  author       = {Pith},
  title        = {Pith review of: Hint at an axion-like particle from GRB 221009A},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5HHEWO4H}},
  note         = {Machine review of arXiv:2412.21175}
}
abstract

The detection by the LHAASO Collaboration of the gamma-ray burst GRB 221009A at redshift $z = 0.151$ with energies up to $(13-18) \, \rm TeV$ challenges conventional physics. Photons emitted with energies above $10 \, \rm TeV$ at this redshift can hardly be observed on Earth due to their interaction with the extragalactic background light (EBL). We show that indeed the LHAASO Collaboration should not have observed photons with energies above $10 \, \rm TeV$ if the state-of-the-art EBL model by Saldana-Lopez et al. is taken into account. A problem therefore arises: the Universe should be more transparent than currently believed. We also show that the issue is solved if we introduce the interaction of photons with axion-like particles (ALPs). ALPs are predicted by String Theory, are among the best candidates for dark matter and can produce spectral and polarization effects on astrophysical sources in the presence of external magnetic fields. In particular, for GRB 221009A, photon-ALP oscillations occur within the crossed magnetized media, i.e. the host galaxy, the extragalactic space, the Milky Way, partially reducing the EBL absorption to a level that explains the LHAASO detection of GRB 221009A and its observed spectrum without the need of contrived choices of parameter values, which are instead compulsory within proposed emission models within conventional physics. This fact regarding GRB 221009A represents a strong hint at the ALP existence, which adds to two other indications coming from blazars, a class of active galactic nuclei.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Extreme value distribution for gamma-ray-burst prompt data -- How unexpected was the BOAT event?

    astro-ph.HE 2025-07 conditional novelty 6.0 of 10

    Using generalized extreme value fits to 30-day maxima from Fermi-GBM and BATSE, the fluence and peak flux of GRB 221009A are extreme outliers, with median return periods of about 1000 years and 140 years respectively.

Reference graph

Works this paper leans on

43 extracted references · 37 canonical work pages · cited by 1 Pith paper

  1. [1]

    32677 (2022)

    LHAASO Collaboration, GCN Circular n. 32677 (2022)

    work page 2022

  2. [2]

    32648 (2022)

    Stargate Collaboration, GCN Circular n. 32648 (2022)

    work page 2022

  3. [3]

    Carpet-2 Collaboration, ATel #15669 (2022)

    work page 2022

  4. [4]

    Dwek and F

    E. Dwek and F. Krennrich, Astropart. Phys. 43, 112 (2013)

    work page 2013

  5. [5]

    Galanti et al

    G. Galanti et al. , Phys. Rev. Lett. 131, 251001 (2023)

    work page 2023

  6. [6]

    Saldana-Lopez et al

    A. Saldana-Lopez et al. , MNRAS 507, 5144 (2021)

    work page 2021

  7. [7]

    Galanti and M

    G. Galanti and M. Roncadelli, Universe 8, 253 (2022)

    2022

  8. [8]

    Ringwald, Phys

    A. Ringwald, Phys. Dark Univ. 1, 116 (2012)

    work page 2012

Show all 43 references
  1. [9]

    Heisenberg and H

    W. Heisenberg and H. Euler, Z. Phys. 98, 714 (1936)

    1936

  2. [10]

    Dobrynina, A

    A. Dobrynina, A. Kartavtsev and G. Raffelt, Phys. Rev. D 91, 083003 (2015)

    2015

  3. [11]

    Maiani, R

    L. Maiani, R. Petronzio and E. Zavattini, Phys. Lett. B 175, 359 (1986)

    1986

  4. [12]

    G. G. Raffelt and L. Stodolsky, Phys. Rev. D 37, 1237 (1988)

    1988

  5. [13]

    De Angelis, M

    A. De Angelis, M. Roncadelli and O. Mansutti, Phys. Rev. D 76, 121301 (2007)

    2007

  6. [14]

    Simet, D

    M. Simet, D. Hooper and P. D. Serpico, Phys. Rev. D 77, 063001 (2008)

    2008

  7. [15]

    Tavecchio, M

    F. Tavecchio, M. Roncadelli, G. Galanti and G. Bonnoli, Phys. Rev. D 86, 085036 (2012)

    2012

  8. [16]

    Galanti, F

    G. Galanti, F. Tavecchio and M. Landoni, MNRAS 491, 5268 (2020)

    2020

  9. [17]

    Galanti, M

    G. Galanti, M. Roncadelli, A. De Angelis and G. F. Bignami, MNRAS 493, 1553 (2020)

    2020

  10. [18]

    Bassan, A

    N. Bassan, A. Mirizzi and M. Roncadelli, JCAP 05, 010 (2010)

    2010

  11. [19]

    Galanti, Phys

    G. Galanti, Phys. Rev. D 105, 083022 (2022)

    2022

  12. [20]

    Galanti, Phys

    G. Galanti, Phys. Rev. D 107, 043006 (2023)

    2023

  13. [21]

    Galanti, M

    G. Galanti, M. Roncadelli, F. Tavecchio and E. Costa, Phys. Rev. D 107, 103007 (2023)

    2023

  14. [22]

    Galanti, M

    G. Galanti, M. Roncadelli and F. Tavecchio, Phys. Rev. D 108, 083017 (2023)

    2023

  15. [23]

    A. J. Levan et al ., Astrophys. J. 946, L28 (2023)

    2023

  16. [24]

    Beck, Astron

    R. Beck, Astron. Astrophys. Rev. 24, 4 (2016)

    2016

  17. [25]

    Lopez-Rodriguez et al

    E. Lopez-Rodriguez et al. , Astrophys. J. 914, 24 (2021)

    2021

  18. [26]

    Galanti and M

    G. Galanti and M. Roncadelli, Phys. Rev. D 98, 043018 (2018)

    2018

  19. [27]

    Galanti and M

    G. Galanti and M. Roncadelli, J. High Energy Astrophys. 20, 1 (2018)

    2018

  20. [28]

    P. P. Kronberg, H. Lesch and U. Hopp, Astrophys. J. 511, 56 (1999)

    1999

  21. [29]

    Furlanetto and A

    S. Furlanetto and A. Loeb, Astrophys. J. 556, 619 (2001)

    2001

  22. [30]

    M. S. Pshirkov, P. G. Tinyakov and F. R. Urban, Phys. Rev. Lett. 116, 191302 (2016)

    2016

  23. [31]

    Jansson and G

    R. Jansson and G. R. Farrar, Astrophys. J. 757, 14 (2012)

    2012

  24. [32]

    J. M. Yao, R. N. Manchester and N. Wang, Astrophys. J. 835, 29 (2017)

    2017

  25. [33]

    Anastassopoulos et al

    V. Anastassopoulos et al. [CAST Collaboration], Nature Physics 13, 584 (2017)

    2017

  26. [34]

    Sisk-Reyn´ eset al

    J. Sisk-Reyn´ eset al. , MNRAS 510, 1264 (2022)

    2022

  27. [35]

    Dessert, D

    C. Dessert, D. Dunsky and B. R. Safdi, Phys. Rev. D 105, 103034 (2022)

    2022

  28. [36]

    LHAASO Collaboration, Science 380, 1390 (2023)

    2023

  29. [37]

    Galanti, M

    G. Galanti, M. Roncadelli and F. Tavecchio, arXiv:2211.06935 (2022)

    2022 arXiv

  30. [38]

    M. M. Gonzalez et al. , Astrophys. J. 944 178 (2023)

    2023

  31. [39]

    Das and S

    S. Das and S. Razzaque, Astron. Astrophys. 670, 12 (2023)

    2023

  32. [40]

    Z.-C. Zhao, Y. Zhou and S. Wang, EPJC 83, 92 (2023)

    2023

  33. [41]

    Sahu et al

    S. Sahu et al. , Astrophys. J. 942 L30 (2023)

    2023

  34. [42]

    LHAASO Collaboration, Science Advances 9, eadj2778 (2023)

    2023

  35. [43]

    Arias et al

    P. Arias et al. , JCAP 06, 013 (2012)

    2012

Pith tools

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