pith. sign in

arxiv: 2606.03881 · v1 · pith:U7JL4PKXnew · submitted 2026-06-02 · 🌌 astro-ph.HE

Suppressed diffusion and gamma-ray emission from the Cygnus Bubble

Pith reviewed 2026-06-28 08:54 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords gamma-ray emissionCygnus OB2cosmic ray accelerationsuppressed diffusionLHAASOstar clustersPeV energiesparticle transport
0
0 comments X

The pith

Explaining the PeV gamma rays from the Cygnus Bubble with steady hadronic acceleration requires strongly suppressed diffusion over 150 parsecs.

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

The paper models the transport of non-thermal particles in the Cygnus OB2 region to test whether acceleration at the cluster wind termination shock or a central source can produce the diffuse gamma-ray emission observed by LHAASO up to PeV energies. It solves the transport equation numerically for steady and bursting scenarios, includes interactions with molecular clouds and penetrating Galactic cosmic rays, and computes the resulting gamma-ray maps and spectra. The calculations show that both the observed energy spectrum and the extended spatial morphology can be matched only when the diffusion coefficient is taken to be spatially dependent and much smaller than the typical Galactic value over a region at least 150 pc across. This result implies that conventional assumptions about cosmic-ray propagation break down in this environment.

Core claim

A spatially dependent Bohm diffusion coefficient is required to reproduce both the spectrum and morphology in the cluster wind scenario. Penetrating Galactic cosmic rays can contribute significantly to the gamma-ray emission above ∼300 TeV. A suppressed diffusion coefficient with respect to the Galactic average in a region extending to at least 150 pc from the cluster center is needed to reproduce the LHAASO morphology. The conclusion is that explaining both the spectrum and morphology of the ∼PeV emission with hadrons accelerated in a non-relativistic steady source requires extreme assumptions, with the possibility that some of the highest-energy gamma rays originate from sources behind the

What carries the argument

Numerical solution of the particle transport equation incorporating a spatially dependent diffusion coefficient to obtain the spatial and spectral distributions of non-thermal particles and their gamma-ray emission via pp interactions.

If this is right

  • Penetrating Galactic cosmic rays can contribute significantly to the gamma-ray emission above ∼300 TeV.
  • The LHAASO morphology requires suppressed diffusion extending at least 150 pc from the cluster center.
  • Some of the highest-energy gamma rays may originate from sources behind the Cygnus association rather than the cluster itself.

Where Pith is reading between the lines

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

  • If diffusion suppression on this scale is common, cosmic-ray propagation models for other star clusters would need revision.
  • Transient or relativistic acceleration sites may become favored explanations for PeV emission once steady non-relativistic sources are ruled out by morphology.
  • Targeted observations of molecular clouds at different distances from the center could test whether particle penetration matches the suppressed-diffusion prediction.

Load-bearing premise

The gamma-ray morphology observed by LHAASO is produced by the spatial distribution of particles whose diffusion is strongly suppressed over a region extending to at least 150 pc from the cluster center.

What would settle it

A direct measurement of the diffusion coefficient near Cygnus OB2 showing values close to the standard Galactic average, or gamma-ray imaging that reveals a morphology inconsistent with the suppressed-diffusion distribution.

Figures

Figures reproduced from arXiv: 2606.03881 by Ben Li, Elena Amato, Pasquale Blasi.

Figure 2
Figure 2. Figure 2: Particle spectrum at Rs and 3Rs (top), and total gamma￾ray emission from the Cygnus bubble (bottom), for uniform Bohm diffusion (dotted lines), non-uniform Bohm diffusion with η2 = 0 (dashed lines), and non-uniform Bohm diffusion with η2 = 0.04 (solid lines). the gamma-ray data points measured by Fermi-LAT, ARGO, HAWC, and LHAASO. For the Kolmogorov and Kraichnan dif￾fusion coefficients, the maximum energy… view at source ↗
Figure 3
Figure 3. Figure 3: The Galactic cosmic-ray proton spectrum as measured in [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Total gamma-ray flux from the Cygnus bubble (top) and gamma-ray morphologies in di [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: As in Figure [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Continuous injection from a central point source, assum [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Impulsive injection from a central point source, assuming [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
read the original abstract

Recent gamma-ray observations indicate that star clusters can be efficient particle accelerators. In particular, LHAASO has detected diffuse gamma-ray emission from Cygnus OB2 extending to $\gtrsim$ PeV energies, indicating that particles are accelerated to at least $\gtrsim$1 PeV. In this work, we study the gamma-ray emission from the Cygnus region assuming particle acceleration either at the termination shock of the cluster wind (WTS) or in an unspecified source at the bubble center, taken to be either steady or bursting. We numerically solve the transport equation for non-thermal particles in all scenarios and derive their spatial and spectral distributions throughout the bubble. We then calculate the gamma-ray emission from pp interactions, including the contribution from particles interacting with the surrounding molecular cloud, which may help explain the extended emission observed by LHAASO. We also include the penetration of Galactic cosmic rays (GCRs) and the resulting shock reacceleration. The predicted emission is compared with Fermi-LAT, HAWC and LHAASO observations. For three diffusion models, we find that a spatially dependent Bohm diffusion coefficient is required to reproduce both the spectrum and morphology in the cluster wind scenario. Penetrating GCRs can contribute significantly to the gamma-ray emission above $\sim$300 TeV. A suppressed diffusion coefficient with respect to the Galactic average in a region extending to at least 150 pc from the cluster center is needed to reproduce the LHAASO morphology. Our conclusion is that explaining both the spectrum and morphology of the $\sim$PeV emission with hadrons accelerated in a non-relativistic steady source requires extreme assumptions. We also speculate on the possibility that some of the highest-energy gamma rays may originate from sources behind the Cygnus association.

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 models non-thermal particle transport in the Cygnus Bubble, considering acceleration at the cluster wind termination shock or a central source (steady or bursting). The transport equation is solved numerically for three diffusion models, with gamma-ray emission computed from pp interactions including molecular-cloud targets and Galactic cosmic-ray penetration/reacceleration. Predicted spectra and morphologies are compared to Fermi-LAT, HAWC, and LHAASO data. The central claim is that reproducing both the ~PeV spectrum and the extended LHAASO morphology requires a spatially dependent Bohm diffusion coefficient suppressed relative to the Galactic average over a region extending to at least 150 pc; the authors conclude that hadronic emission from a non-relativistic steady source therefore demands extreme assumptions and speculate that some highest-energy photons may originate behind the association.

Significance. If robust, the result would underscore the tension between extended PeV-scale emission and conventional diffusion in star clusters, providing a detailed numerical framework that incorporates multiple source scenarios, cloud targets, and GCR contributions. The explicit inclusion of GCR penetration as a significant component above ~300 TeV is a constructive element. The work supplies a useful benchmark for similar regions even if the 'extreme assumptions' conclusion requires qualification.

major comments (2)
  1. [§3] §3 (diffusion models): the conclusion that 'extreme assumptions' are required rests on only three diffusion models. The abstract states that a spatially dependent Bohm coefficient is needed to match both spectrum and LHAASO morphology; without evidence that these three exhaust plausible turbulence spectra (e.g., Kolmogorov/Kraichnan) or radial suppression profiles, other prescriptions could potentially reproduce the observed morphology with milder or more localized suppression, directly affecting the load-bearing claim.
  2. [Abstract] Abstract and results section on LHAASO comparison: the suppression factor is adjusted to reproduce the LHAASO morphology (and spectrum), which introduces circularity; while the paper compares to independent observations, the central requirement of suppression over ≥150 pc is defined by the fit itself, weakening the assertion that the data independently demand extreme conditions.
minor comments (2)
  1. [§2] Notation for the diffusion coefficient D(r,E) should be defined once and used consistently when comparing the three models.
  2. [Figure 4] Figure captions for the morphology maps should explicitly state the energy range and the contribution breakdown (WTS vs. GCR) shown in each panel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and constructive comments. We address each major comment point by point below, providing the strongest honest defense of the manuscript while acknowledging where clarification or qualification is warranted. Revisions will be made to improve the presentation without altering the core scientific conclusions.

read point-by-point responses
  1. Referee: [§3] §3 (diffusion models): the conclusion that 'extreme assumptions' are required rests on only three diffusion models. The abstract states that a spatially dependent Bohm coefficient is needed to match both spectrum and LHAASO morphology; without evidence that these three exhaust plausible turbulence spectra (e.g., Kolmogorov/Kraichnan) or radial suppression profiles, other prescriptions could potentially reproduce the observed morphology with milder or more localized suppression, directly affecting the load-bearing claim.

    Authors: We appreciate the referee raising this point about model exhaustiveness. The three diffusion models were selected to bracket the range of plausible behaviors, with the spatially dependent Bohm case representing the minimal diffusion (strongest suppression) needed to confine PeV particles over the observed scales. Standard Kolmogorov or Kraichnan spectra typically yield larger diffusion coefficients at high rigidities than Bohm, which would necessitate comparable or greater suppression to reproduce the LHAASO morphology; thus, they would not relax the requirement for extreme conditions. We will add a short discussion in §3 explaining this rationale and noting the limitations of not exploring every possible radial profile. This qualifies the claim appropriately while preserving the result that strong, extended suppression is required in the scenarios considered. revision: partial

  2. Referee: [Abstract] Abstract and results section on LHAASO comparison: the suppression factor is adjusted to reproduce the LHAASO morphology (and spectrum), which introduces circularity; while the paper compares to independent observations, the central requirement of suppression over ≥150 pc is defined by the fit itself, weakening the assertion that the data independently demand extreme conditions.

    Authors: We disagree that the procedure introduces problematic circularity. Adjusting model parameters to match data is the standard method for determining the physical conditions implied by observations. The LHAASO morphology supplies an independent spatial constraint that cannot be satisfied by varying only the normalization or spectral index; only a diffusion coefficient suppressed relative to the Galactic average and extending to at least 150 pc reproduces the observed extent. The spectrum and morphology are compared jointly to multiple independent datasets (Fermi-LAT, HAWC, LHAASO), and the GCR contribution is treated separately. We will revise the abstract and results section to clarify this distinction between the fitting process and the conclusion that steady hadronic acceleration under conventional diffusion is disfavored. No change to the scientific interpretation is required. revision: partial

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper numerically solves the transport equation for non-thermal particles under three explicit diffusion models, computes gamma-ray emission from pp interactions (including molecular clouds and GCR penetration), and directly compares the output spectra and morphologies to independent Fermi-LAT, HAWC, and LHAASO observations. The statement that a suppressed, spatially dependent Bohm-like coefficient is required follows from the mismatch between standard models and data; this is an explicit modeling result, not a self-definitional loop, a fitted parameter renamed as a prediction, or any load-bearing self-citation. The derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on fitted diffusion parameters chosen to match observations and on standard assumptions about particle transport and hadronic emission; no new physical entities are introduced.

free parameters (2)
  • diffusion suppression factor
    Value chosen to reproduce LHAASO morphology; suppressed relative to Galactic average over region of at least 150 pc from center
  • normalization of accelerated particle spectrum
    Adjusted to match observed gamma-ray flux levels
axioms (2)
  • standard math The cosmic ray transport equation governs the spatial and spectral distribution of non-thermal particles
    Invoked when numerically solving for particle distributions throughout the bubble
  • domain assumption Gamma-ray emission is calculated from pp interactions of accelerated hadrons and penetrating GCRs
    Used to derive the predicted emission compared with Fermi-LAT, HAWC and LHAASO data

pith-pipeline@v0.9.1-grok · 5863 in / 1557 out tokens · 37698 ms · 2026-06-28T08:54:13.628366+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

55 extracted references · 54 canonical work pages · 11 internal anchors

  1. [1]

    , keywords =

    Galactic factories of cosmic-ray electrons and positrons. , keywords =. doi:10.1103/PhysRevD.103.083010 , archivePrefix =. 2010.11955 , primaryClass =

  2. [2]

    Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks

    Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks. , keywords =. doi:10.1103/PhysRevLett.114.085003 , archivePrefix =. 1412.0672 , primaryClass =

  3. [3]

    , keywords =

    The origin of very high-energy diffuse -ray emission: The case for galactic source cocoons. , keywords =. doi:10.1051/0004-6361/202554796 , archivePrefix =. 2503.14651 , primaryClass =

  4. [4]

    , keywords =

    Gamma rays from star clusters and implications for the origin of Galactic cosmic rays. , keywords =. doi:10.1051/0004-6361/202453017 , archivePrefix =. 2501.16097 , primaryClass =

  5. [5]

    , keywords =

    Understanding the TeV -ray emission surrounding the young massive star cluster Westerlund 1. , keywords =. doi:10.1051/0004-6361/202245444 , archivePrefix =. 2301.10496 , primaryClass =

  6. [6]

    , keywords =

    Multiple emission components in the Cygnus cocoon detected from Fermi-LAT observations. , keywords =. doi:10.1051/0004-6361/202245573 , archivePrefix =. 2301.04504 , primaryClass =

  7. [7]

    , keywords =

    Deciphering the gamma-ray emission in the Cygnus region. , keywords =. doi:10.1051/0004-6361/202555531 , archivePrefix =. 2508.21644 , primaryClass =

  8. [8]

    , keywords =

    Cosmic ray production in superbubbles. , keywords =. doi:10.1093/mnras/stac543 , archivePrefix =. 2201.07488 , primaryClass =

  9. [9]

    , keywords =

    Supernova remnants in super bubbles acting as cosmic ray accelerators. , keywords =. doi:10.1051/0004-6361/202554145 , archivePrefix =. 2505.02523 , primaryClass =

  10. [10]

    , keywords =

    A deep spectromorphological study of the -ray emission surrounding the young massive stellar cluster Westerlund 1. , keywords =. doi:10.1051/0004-6361/202244323 , archivePrefix =. 2207.10921 , primaryClass =

  11. [11]

    Monthly Notices of the Royal Astronomical Society , volume =

    Blasi, Pasquale and Morlino, Giovanni , title =. Monthly Notices of the Royal Astronomical Society , volume =. 2023 , month =. doi:10.1093/mnras/stad1662 , url =

  12. [12]

    , keywords =

    Particle acceleration in winds of star clusters. , keywords =. doi:10.1093/mnras/stab690 , archivePrefix =. 2102.09217 , primaryClass =

  13. [13]

    arXiv e-prints , keywords =

    Probing star clusters as cosmic ray factories. arXiv e-prints , keywords =. doi:10.48550/arXiv.2307.03477 , archivePrefix =. 2307.03477 , primaryClass =

  14. [14]

    , keywords =

    Gamma-Rays from Active Regions in the Galaxy - the Possible Contribution of Stellar Winds. , keywords =. doi:10.1007/BF00167503 , adsurl =

  15. [15]

    , keywords =

    Cosmic-ray acceleration at stellar wind terminal shocks. , keywords =. doi:10.1086/163652 , adsurl =

  16. [16]

    , keywords =

    Cosmic rays and gamma-rays from OB stars. , keywords =. doi:10.1086/159623 , adsurl =

  17. [17]

    Superbubbles and energetic particles in the Galaxy. I. Collective effects of particle acceleration. , keywords =. doi:10.1051/0004-6361:20041269 , archivePrefix =. astro-ph/0405531 , primaryClass =

  18. [18]

    Galactic cosmic rays after the AMS-02 observations

    Galactic cosmic rays after the AMS-02 observations. , keywords =. doi:10.1103/PhysRevD.99.103023 , archivePrefix =. 1904.10220 , primaryClass =

  19. [19]

    International Journal of Modern Physics D , keywords =

    The origin of Galactic cosmic rays: Challenges to the standard paradigm. International Journal of Modern Physics D , keywords =. doi:10.1142/S0218271819300222 , archivePrefix =. 1903.11584 , primaryClass =

  20. [20]

    , year = 1972, month = aug, volume =

    The Dynamical Effects of Stellar Mass Loss on Diffuse Nebulae. , year = 1972, month = aug, volume =

  21. [21]

    , keywords =

    Interstellar bubbles. , keywords =. doi:10.1086/181908 , adsurl =

  22. [22]

    Interstellar bubbles. II. Structure and evolution. , keywords =. doi:10.1086/155692 , adsurl =

  23. [23]

    Lack of thermal energy in superbubbles: hint of cosmic rays?

    Lack of thermal energy in superbubbles: hint of cosmic rays?. , keywords =. doi:10.1093/mnras/stx2427 , archivePrefix =. 1705.10448 , primaryClass =

  24. [24]

    , keywords =

    Hydrodynamic simulation of Cygnus OB2: the absence of a cluster wind termination shock. , keywords =. doi:10.1093/mnras/stae1627 , archivePrefix =. 2406.13589 , primaryClass =

  25. [25]

    Energy spectra of gamma-rays, electrons and neutrinos produced at proton-proton interactions in the very high energy regime

    Energy spectra of gamma rays, electrons, and neutrinos produced at proton-proton interactions in the very high energy regime. , keywords =. doi:10.1103/PhysRevD.74.034018 , archivePrefix =. astro-ph/0606058 , primaryClass =

  26. [26]

    Parametrization of gamma-ray production cross-sections for pp interactions in a broad proton energy range from the kinematic threshold to PeV energies

    Parametrization of gamma-ray production cross sections for p p interactions in a broad proton energy range from the kinematic threshold to PeV energies. , keywords =. doi:10.1103/PhysRevD.90.123014 , archivePrefix =. 1406.7369 , primaryClass =

  27. [27]

    , keywords =

    Pion Production Momentum Loss of Cosmic Ray Hadrons. , keywords =. doi:10.1088/0004-637X/802/2/114 , adsurl =

  28. [28]

    , keywords =

    Morphology of Gamma-Ray Halos around Middle-aged Pulsars: Influence of the Pulsar Proper Motion. , keywords =. doi:10.3847/1538-4357/ac235e , archivePrefix =. 2010.15731 , primaryClass =

  29. [29]

    , keywords =

    Slow Diffusion is Necessary to Explain the -Ray Pulsar Halos. , keywords =. doi:10.3847/1538-4357/ac8b8a , archivePrefix =. 2107.07395 , primaryClass =

  30. [30]

    , keywords =

    Prospects for the detection of very-high-energy pulsars with LHAASO and SWGO. , keywords =. doi:10.1093/mnras/stae1497 , archivePrefix =. 2407.00262 , primaryClass =

  31. [31]

    , keywords =

    Cygnus X-3 as a semi-hidden PeVatron. , keywords =. doi:10.1051/0004-6361/202554701 , archivePrefix =. 2503.11448 , primaryClass =

  32. [32]

    Fermi Large Area Telescope Fourth Source Catalog.Astron

    Fermi Large Area Telescope Fourth Source Catalog. , keywords =. doi:10.3847/1538-4365/ab6bcb , archivePrefix =. 1902.10045 , primaryClass =

  33. [33]

    Science , keywords =

    A Cocoon of Freshly Accelerated Cosmic Rays Detected by Fermi in the Cygnus Superbubble. Science , keywords =. doi:10.1126/science.1210311 , adsurl =

  34. [34]

    Massive Stars as Major Factories of Galactic Cosmic Rays

    Massive stars as major factories of Galactic cosmic rays. Nature Astronomy , keywords =. doi:10.1038/s41550-019-0724-0 , archivePrefix =. 1804.02331 , primaryClass =

  35. [35]

    Identification of the TeV Gamma-ray Source ARGO J2031+4157 with the Cygnus Cocoon

    Identification of the TeV Gamma-Ray Source ARGO J2031+4157 with the Cygnus Cocoon. , keywords =. doi:10.1088/0004-637X/790/2/152 , archivePrefix =. 1406.6436 , primaryClass =

  36. [36]

    Nature Astronomy , keywords =

    HAWC observations of the acceleration of very-high-energy cosmic rays in the Cygnus Cocoon. Nature Astronomy , keywords =. doi:10.1038/s41550-021-01318-y , archivePrefix =. 2103.06820 , primaryClass =

  37. [37]

    Science Bulletin , keywords =

    An ultrahigh-energy -ray bubble powered by a super PeVatron. Science Bulletin , keywords =. doi:10.1016/j.scib.2023.12.040 , archivePrefix =. 2310.10100 , primaryClass =

  38. [38]

    Cygnus X-3: A variable petaelectronvolt gamma-ray source

    Cygnus X-3: A variable petaelectronvolt gamma-ray source. arXiv e-prints , keywords =. doi:10.48550/arXiv.2512.16638 , archivePrefix =. 2512.16638 , primaryClass =

  39. [39]

    , year = 2021, month = jun, volume =

    Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 -ray Galactic sources. , year = 2021, month = jun, volume =. doi:10.1038/s41586-021-03498-z , adsurl =

  40. [40]

    , keywords =

    Solar Modulation of Galactic Cosmic-Ray Protons Based on a Modified Force-field Approach. , keywords =. doi:10.3847/1538-4357/ac1fe8 , adsurl =

  41. [41]

    Astroparticle Physics , keywords =

    The shape of the cosmic ray proton spectrum. Astroparticle Physics , keywords =. doi:10.1016/j.astropartphys.2020.102441 , archivePrefix =. 1911.01311 , primaryClass =

  42. [42]

    , keywords =

    Measurement of Ultra-High-Energy Diffuse Gamma-Ray Emission of the Galactic Plane from 10 TeV to 1 PeV with LHAASO-KM2A. , keywords =. doi:10.1103/PhysRevLett.131.151001 , archivePrefix =. 2305.05372 , primaryClass =

  43. [43]

    C., Stone, E

    Galactic Cosmic Rays in the Local Interstellar Medium: Voyager 1 Observations and Model Results. , keywords =. doi:10.3847/0004-637X/831/1/18 , adsurl =

  44. [44]

    C., Cummings, A

    Cosmic ray measurements from Voyager 2 as it crossed into interstellar space. Nature Astronomy , year = 2019, month = nov, volume =. doi:10.1038/s41550-019-0928-3 , adsurl =

  45. [45]

    New very local interstellar spectra for electrons, positrons, protons and light cosmic ray nuclei

    New Very Local Interstellar Spectra for Electrons, Positrons, Protons, and Light Cosmic Ray Nuclei. , keywords =. doi:10.3847/1538-4357/ab1e4a , archivePrefix =. 1902.10438 , primaryClass =

  46. [46]

    , keywords =

    Precision Measurement of the Proton Flux in Primary Cosmic Rays from Rigidity 1 GV to 1.8 TV with the Alpha Magnetic Spectrometer on the International Space Station. , keywords =. doi:10.1103/PhysRevLett.114.171103 , adsurl =

  47. [47]

    Science Advances , keywords =

    Measurement of the cosmic ray proton spectrum from 40 GeV to 100 TeV with the DAMPE satellite. Science Advances , keywords =. doi:10.1126/sciadv.aax3793 , archivePrefix =. 1909.12860 , primaryClass =

  48. [48]

    Proton and Helium Spectra from the CREAM-III Flight

    Proton and Helium Spectra from the CREAM-III Flight. , keywords =. doi:10.3847/1538-4357/aa68e4 , archivePrefix =. 1704.02512 , primaryClass =

  49. [49]

    Aharonian and Y.X

    Zhen Cao and F. Aharonian and Y.X. Bai and Y.W. Bao and D. Bastieri and X.J. Bi and Y.J. Bi and W. Bian and A.V. Bukevich and C.M. Cai and W.Y. Cao and Zhe Cao and J. Chang and J.F. Chang and A.M. Chen and E.S. Chen and G.H. Chen and H.X. Chen and Liang Chen and Long Chen and M.J. Chen and M.L. Chen and Q.H. Chen and S. Chen and S.H. Chen and S.Z. Chen an...

  50. [50]

    Aartsen, M

    Cosmic ray spectrum and composition from PeV to EeV using 3 years of data from IceTop and IceCube. , keywords =. doi:10.1103/PhysRevD.100.082002 , archivePrefix =. 1906.04317 , primaryClass =

  51. [51]

    Advances in Space Research , keywords =

    Modeling of GeV-TeV gamma-ray emission of Cygnus Cocoon. Advances in Space Research , keywords =. doi:10.1016/j.asr.2022.01.029 , archivePrefix =. 2203.09842 , primaryClass =

  52. [52]

    , keywords =

    On the Potential of Bright, Young Pulsars to Power Ultrahigh Gamma-Ray Sources. , keywords =. doi:10.3847/2041-8213/ac66cf , archivePrefix =. 2204.09440 , primaryClass =

  53. [53]

    Journal of High Energy Astrophysics , keywords =

    ASTRI Mini-Array core science at the Observatorio del Teide. Journal of High Energy Astrophysics , keywords =. doi:10.1016/j.jheap.2022.05.005 , archivePrefix =. 2208.03177 , primaryClass =

  54. [54]

    Astroparticle Physics , year = 2013, month = mar, volume =

    Introducing the CTA concept. Astroparticle Physics , year = 2013, month = mar, volume =. doi:10.1016/j.astropartphys.2013.01.007 , adsurl =

  55. [55]

    , keywords =

    Modeling Cosmic-ray Electron Spectra and Synchrotron Emission in the Multiphase Interstellar Medium. , keywords =. doi:10.3847/1538-4357/ade5ad , archivePrefix =. 2507.00142 , primaryClass =