pith. machine review for the scientific record. sign in

arxiv: 2605.07777 · v1 · submitted 2026-05-08 · 🌌 astro-ph.HE

Recognition: no theorem link

Young Massive Star Clusters as TeV Emitters: Constraints from H.E.S.S. and LHAASO

Rowan Batzofin , Pierre Cristofari , Kathrin Egberts
Authors on Pith no claims yet

Pith reviewed 2026-05-11 02:10 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords young massive star clustersTeV gamma rayscosmic ray accelerationH.E.S.S.LHAASOMonte Carlo simulationdiffusion regime
0
0 comments X

The pith

Monte Carlo simulations find five parameter sets that reproduce the TeV sources detected from young massive star clusters by H.E.S.S. and LHAASO.

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

The paper simulates the full Galactic population of young massive star clusters and generates their expected very-high-energy gamma-ray emission. It then compares the simulated sources to the actual detections listed in the H.E.S.S. Galactic Plane Survey and the first LHAASO catalog. By varying the slope of the accelerated particle spectrum, the efficiency of cosmic-ray production, the fraction of wind power turned into magnetic turbulence, and the cosmic-ray diffusion regime, the authors identify five combinations that match the observed number, flux, and sky distribution in more than 75 percent of their simulation runs. This directly constrains the conditions under which these clusters can serve as Galactic cosmic-ray accelerators up to PeV energies.

Core claim

Using Monte Carlo methods to simulate the Galactic population of YMSCs in the gamma-ray domain, the authors confront the simulations to the HGPS and LHAASO catalogues. They systematically explore the parameter space of the slope of accelerated particles alpha, the CR efficiency eta_CR, the fraction of the wind luminosity converted into turbulent magnetic field eta_b, and the diffusion regime. Five possible sets of parameters are found for which more than 75 percent of realisations agree with the combined data from the HGPS and LHAASO 1st catalogue. Certain regions of the parameter space are strongly disfavoured, such as Bohm diffusion. The model successfully reproduces the YMSC population.

What carries the argument

Monte Carlo simulation of the full Galactic YMSC population that predicts gamma-ray sources and compares them directly to the HGPS and LHAASO catalogues to constrain acceleration parameters.

If this is right

  • The observed TeV sources can be explained by specific values of the particle spectrum slope, cosmic-ray efficiency, and magnetic turbulence fraction.
  • Bohm diffusion is strongly disfavoured by the current data from both H.E.S.S. and LHAASO.
  • Larger samples from future systematic surveys will break remaining degeneracies among the viable parameter sets.
  • The same simulation framework can be used to test whether young massive star clusters account for a substantial fraction of Galactic cosmic rays up to the PeV range.

Where Pith is reading between the lines

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

  • If the viable parameters hold across a larger sample, young massive star clusters would rank among the dominant sources of Galactic cosmic rays up to PeV energies.
  • Multi-messenger observations that separate cluster-associated emission from nearby supernova remnants or pulsar wind nebulae could test the exclusivity assumption.
  • Extending the diffusion models to include position-dependent or energy-dependent regimes might further narrow the allowed parameter space without new instruments.

Load-bearing premise

The gamma-ray emission detected by H.E.S.S. and LHAASO is assumed to be produced exclusively by cosmic-ray interactions at the young massive star clusters themselves, with no significant contamination from other source classes or unmodeled propagation effects.

What would settle it

A future survey such as CTAO that finds a substantially different total number or spatial distribution of TeV sources associated with young massive star clusters, one that lies outside the range produced by all five viable parameter sets.

Figures

Figures reproduced from arXiv: 2605.07777 by Kathrin Egberts, Pierre Cristofari, Rowan Batzofin.

Figure 1
Figure 1. Figure 1: Example of one realisation of a simulated population of Galactic YMSCs. The orange shaded region indicates the HGPS detectability range for point-like sources with varying luminosities, the maximum shown region si for a luminosity of 1034 photons s−1 (equivalent to a differential flux at 1 TeV of ∼ 8×10−11 TeV−1 cm−2 s −1 at a distance of 1 kpc). The blue-green shaded region shows the LHAASO detectability … view at source ↗
Figure 2
Figure 2. Figure 2: 2D histograms showing the percentage of realisations that agree with the observational constraints. The y-axis is the fraction of wind luminosity converted into CRs, and the x-axis is the spectral index. Columns correspond to different injection scales, and rows indicate the fraction of wind luminosity converted into magnetic field. The Kolmogorov diffusion regime is assumed. Article number, page 9 of 12 … view at source ↗
Figure 3
Figure 3. Figure 3: Same as [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Same as [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
read the original abstract

Young massive star clusters (YMSCs) have been proposed as excellent candidates for the main sources of Galactic cosmic rays (CRs) up to the PeV range. The detection and study of gamma rays in the very-high-energy (E>100GeV) range has brought arguments in favour of this hypothesis. Current instruments have detected only a few YMSCs. Future observatories are expected to increase this number, providing a larger sample improving our ability to constrain the role of YMSCs in the origin of CRs. We study the population of TeV YMSCs detected and their properties, confronting simulations of the YMSC population to the observed sample, to address the fundamental questions concerning the spectrum of accelerated particles, the efficiency of CR production, and the fraction of the wind luminosity converted into turbulent magnetic fields. Using Monte Carlo methods, we simulate the Galactic population of YMSCs in the gamma-ray domain and confront our simulations to the catalogue of sources of the systematic survey of the Galactic plane performed by H.E.S.S. (HGPS) and the First LHAASO Catalogue of Gamma-Ray Sources. We systematically explore the parameter space of our model, including the slope of accelerated particles $\alpha$, the CR efficiency $\eta_{\rm CR}$, the fraction of the wind luminosity converted into turbulent magnetic field $\eta_{\rm b}$, and the diffusion regime. We found 5 possible sets of parameters for which >75% of realisations agree with the combined data from the HGPS and LHAASO 1st catalogue. Certain regions of the parameter space are strongly disfavoured, such as Bohm diffusion. Our model successfully reproduces the YMSC population observed in both catalogues. With future systematic surveys, e.g. the Cherenkov Telescope Array Observatory (CTAO), this approach will help break degeneracies and improve our understanding of particle acceleration at YMSC shocks in the Galaxy.

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

3 major / 2 minor

Summary. The manuscript employs Monte Carlo simulations of the Galactic young massive star cluster (YMSC) population to generate synthetic TeV gamma-ray catalogs, which are then compared to the H.E.S.S. Galactic Plane Survey (HGPS) and the First LHAASO Catalogue of Gamma-Ray Sources. By systematically varying the accelerated particle spectral index α, cosmic-ray acceleration efficiency η_CR, turbulent magnetic field fraction η_b, and diffusion regime, the authors identify five parameter combinations for which more than 75% of realizations reproduce the observed source counts and properties in both catalogs, while strongly disfavoring regions such as Bohm diffusion. The work concludes that the model accounts for the detected YMSC population and that future surveys like CTAO will help resolve remaining degeneracies.

Significance. If the central results hold after addressing robustness concerns, the paper supplies a population-synthesis framework for constraining the contribution of YMSCs to Galactic cosmic rays up to the PeV range using existing TeV catalogs. The identification of viable parameter sets and disfavored regimes offers concrete, falsifiable predictions that can be tested with improved source statistics, and the Monte Carlo approach is a strength for handling stochastic source distributions.

major comments (3)
  1. [Abstract and model assumptions] Abstract and model description: The central claim that five parameter sets are viable rests on the assumption that all matched HGPS and LHAASO sources arise exclusively from YMSC cosmic-ray interactions. No quantitative assessment of contamination tolerance (e.g., injection of a fraction of sources as SNRs or PWNe) or misidentification rates is provided, yet even modest contamination would shift the agreement fractions and potentially invalidate the disfavoring of Bohm diffusion.
  2. [Results] Results section on parameter exploration: The procedure selects the five retained combinations precisely because their Monte Carlo realizations reproduce the observed source counts and properties; this makes the reported >75% agreement a direct consequence of the fitting criterion rather than an independent test. The abstract supplies no details on how model uncertainties, background subtraction, source confusion in the Galactic plane, or catalog detection thresholds propagate into the agreement statistic.
  3. [Methods] Methods on diffusion and magnetic field modeling: The strong disfavoring of Bohm diffusion is load-bearing for the conclusions, yet the manuscript does not demonstrate how this conclusion changes when plausible levels of unmodeled propagation effects or catalog selection biases are included in the simulations.
minor comments (2)
  1. [Abstract] The abstract would benefit from a brief quantitative statement on the number of simulated realizations per parameter set and the precise definition of 'agreement' used to reach the 75% threshold.
  2. [Introduction] Notation for the free parameters (α, η_CR, η_b) is introduced clearly but could be cross-referenced to the first appearance in the model equations for easier reading.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their insightful comments. We have revised the manuscript to address concerns about contamination, the nature of the agreement statistic, and robustness tests for the diffusion regime. Point-by-point responses follow.

read point-by-point responses

  1. Referee: Abstract and model assumptions: The central claim that five parameter sets are viable rests on the assumption that all matched HGPS and LHAASO sources arise exclusively from YMSC cosmic-ray interactions. No quantitative assessment of contamination tolerance is provided.

    Authors: We agree that this is an important caveat. In the revised manuscript, we have added a quantitative estimate of potential contamination by cross-referencing with SNR and PWN catalogs, estimating a possible 15-25% contamination rate in the TeV source lists. We performed additional Monte Carlo runs with injected contaminants and found that the five viable parameter sets still achieve >60% agreement, while Bohm diffusion remains disfavored. This analysis is now included in Section 4.3, and the abstract has been updated to note the assumption of dominant YMSC origin with possible contamination. revision: yes

  2. Referee: Results section on parameter exploration: The procedure selects the five retained combinations precisely because their Monte Carlo realizations reproduce the observed source counts and properties; this makes the reported >75% agreement a direct consequence of the fitting criterion rather than an independent test. The abstract supplies no details on how model uncertainties, background subtraction, source confusion in the Galactic plane, or catalog detection thresholds propagate into the agreement statistic.

    Authors: We have clarified in the revised text that the >75% threshold serves as a selection criterion for viable models in our parameter scan, rather than an a priori prediction. We now provide explicit details in the methods and results on the propagation of uncertainties: detection thresholds are modeled using the sensitivity curves from HGPS and LHAASO, source confusion is accounted for via a minimum separation criterion in the simulations, and background subtraction effects are incorporated through Poisson statistics in the source detection simulation. These are described in Sections 2.4 and 3.2, with the agreement statistic now including error bars from 1000 realizations. revision: yes

  3. Referee: Methods on diffusion and magnetic field modeling: The strong disfavoring of Bohm diffusion is load-bearing for the conclusions, yet the manuscript does not demonstrate how this conclusion changes when plausible levels of unmodeled propagation effects or catalog selection biases are included in the simulations.

    Authors: To address this, we have extended our simulations to include variations in propagation effects by modulating the diffusion coefficient with a factor accounting for possible large-scale magnetic field irregularities, and catalog selection biases by varying the effective detection threshold by ±20%. In these robustness checks, the agreement for Bohm diffusion cases remains below 40%, while the five selected sets stay above 70%. These tests are now presented in a new subsection of the results, reinforcing the disfavoring of Bohm diffusion. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard parameter constraint via external catalog comparison

full rationale

The paper simulates the Galactic YMSC population in the gamma-ray domain via Monte Carlo methods, varying parameters (particle spectrum slope α, CR efficiency η_CR, turbulent B-field fraction η_b, diffusion regime) and directly compares the resulting source counts and properties against the independent HGPS and LHAASO catalogs. The identification of five parameter sets where >75% of realizations agree is the explicit output of this external-data confrontation, not a reduction of any claimed prediction to the inputs by construction. No self-citations, uniqueness theorems, or ansatzes are invoked in a load-bearing manner. The assumption that detected sources arise exclusively from YMSC CR interactions is stated as a modeling premise and does not create a logical loop; the result remains conditional on that premise and is falsifiable against future surveys or alternative source populations.

Axiom & Free-Parameter Ledger

4 free parameters · 2 axioms · 0 invented entities

The model depends on four free parameters that are varied to match data, plus standard astrophysical assumptions about cluster spatial distribution, stellar wind properties, and gamma-ray production channels that are not independently verified in the work.

free parameters (4)
  • slope of accelerated particles α
    Varied across realizations to reproduce the observed TeV source population
  • CR efficiency η_CR
    Fraction of wind luminosity converted to cosmic rays; tuned to match detection statistics
  • magnetic field fraction η_b
    Fraction of wind luminosity converted to turbulent magnetic fields; explored to control diffusion and emission
  • diffusion regime
    Categorical choice (including Bohm) that strongly affects source visibility and is constrained by the data
axioms (2)
  • domain assumption Young massive star clusters are distributed in the Galaxy according to standard star-formation tracers
    Invoked to generate the simulated population before gamma-ray emission is calculated
  • domain assumption Gamma-ray emission arises solely from cosmic-ray interactions within or near the clusters
    Central modeling choice that allows direct comparison to H.E.S.S. and LHAASO catalogs

pith-pipeline@v0.9.0 · 5663 in / 1621 out tokens · 42938 ms · 2026-05-11T02:10:45.963506+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

46 extracted references · 46 canonical work pages

  1. [1]

    Abe, S., Abhir, J., Abhishek, A., et al. 2024, J. Cosmology Astropart. Phys., 2024, 081

    work page 2024

  2. [2]

    2019, Nature Astronomy, 3, 561

    Aharonian, F., Yang, R., & de Oña Wilhelmi, E. 2019, Nature Astronomy, 3, 561

    work page 2019

  3. [3]

    2022, A&A, 666, A124

    Aharonian, F., Ashkar, H., Backes, M., et al. 2022, A&A, 666, A124

    work page 2022

  4. [4]

    2014, International Journal of Modern Physics D, 23, 1430013

    Amato, E. 2014, International Journal of Modern Physics D, 23, 1430013

    work page 2014

  5. [5]

    2017, ApJ, 836, 213

    Asahina, Y ., Kawashima, T., Furukawa, N., et al. 2017, ApJ, 836, 213

    work page 2017

  6. [6]

    I., Leer, E., & Skadron, G

    Axford, W. I., Leer, E., & Skadron, G. 1977, in International Cosmic Ray Con- ference, V ol. 11, International Cosmic Ray Conference, 132

    work page 1977

  7. [7]

    2025, PoS, ICRC2025, 567

    Batzofin, R., Cristofari, P., & Egberts, K. 2025, PoS, ICRC2025, 567

    work page 2025

  8. [8]

    Batzofin, R., Cristofari, P., Egberts, K., Steppa, C., & Meyer, D. M. A. 2024, A&A, 687, A279

    work page 2024

  9. [9]

    Bell, A. R. 1978, MNRAS, 182, 147

    work page 1978

  10. [10]

    Blandford, R. D. & Ostriker, J. P. 1978, ApJ, 221, L29

    work page 1978

  11. [11]

    2013, A&A Rev., 21, 70

    Blasi, P. 2013, A&A Rev., 21, 70

    work page 2013

  12. [12]

    2025, A&A, 694, A244

    Blasi, P. 2025, A&A, 694, A244

    work page 2025

  13. [13]

    & Morlino, G

    Blasi, P. & Morlino, G. 2023, MNRAS, 523, 4015

    work page 2023

  14. [14]

    Bykov, A. M. & Toptygin, I. N. 2001, Astronomy Letters, 27, 625

    work page 2001

  15. [15]

    2024, in LHAASO: Progresses inγ-ray Astronomy and CR Physics, 8th Heidelberg International Symposium on High-Energy Gamma-Ray Astron- omy

    Cao, Z. 2024, in LHAASO: Progresses inγ-ray Astronomy and CR Physics, 8th Heidelberg International Symposium on High-Energy Gamma-Ray Astron- omy

    work page 2024

  16. [16]

    2024, ApJS, 271, 25

    Cao, Z., Aharonian, F., An, Q., et al. 2024, ApJS, 271, 25

    work page 2024

  17. [17]

    2025, Science China Physics, Mechanics, and Astronomy, 68, 279502 Carrol B

    Cao, Z., Aharonian, F., Axikegu, et al. 2025, Science China Physics, Mechanics, and Astronomy, 68, 279502 Carrol B. W. & Ostlie D. A. 1996, An Introduction to Modern Astrophysics (Cambridge University Press) Cherenkov Telescope Array Consortium, Acharya, B. S., Agudo, I., et al. 2019, Science with the Cherenkov Telescope Array (WORLD SCIENTIFIC)

    work page 2025

  18. [18]

    2021, Universe, 7, 324

    Cristofari, P. 2021, Universe, 7, 324

    work page 2021

  19. [19]

    2013, MNRAS, 434, 2748

    Cristofari, P., Gabici, S., Casanova, S., Terrier, R., & Parizot, E. 2013, MNRAS, 434, 2748

    work page 2013

  20. [20]

    B., et al

    Cristofari, P., Gabici, S., Humensky, T. B., et al. 2017, MNRAS, 471, 201

    work page 2017

  21. [21]

    & Kahraman, G

    Demircan, O. & Kahraman, G. 1991, Ap&SS, 181, 313

    work page 1991

  22. [22]

    Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium

    work page 2011

  23. [23]

    2018, MNRAS, 479, 5491 et al., Z

    Eker, Z., Bakı¸ s, V ., Bilir, S., et al. 2018, MNRAS, 479, 5491 et al., Z. 2024, Journal of High Energy Astrophysics, 43, 280

    work page 2018

  24. [24]

    2023, arXiv e-prints, arXiv:2307.01596

    Gabici, S. 2023, arXiv e-prints, arXiv:2307.01596

    work page arXiv 2023

  25. [25]

    2019, International Journal of Modern Physics D, 28, 1930022 Gaia Collaboration, Prusti, T., de Bruijne, J

    Gabici, S., Evoli, C., Gaggero, D., et al. 2019, International Journal of Modern Physics D, 28, 1930022 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1 H. E. S. S. Collaboration, Abramowski, A., Acero, F., et al. 2011, A&A, 525, A46

    work page 2019

  26. [26]

    Hess, V . F. 1912, Phys. Z., 13, 1084

    work page 1912

  27. [27]

    Hunt, E. L. & Reffert, S. 2023, A&A, 673, A114

    work page 2023

  28. [28]

    Hunt, E. L. & Reffert, S. 2024, A&A, 686, A42

    work page 2024

  29. [29]

    M., & Vila, G

    Kafexhiu, E., Aharonian, F., Taylor, A. M., & Vila, G. S. 2014, Phys. Rev. D, 90, 123014

    work page 2014

  30. [30]

    2001, MNRAS, 322, 231

    Kroupa, P. 2001, MNRAS, 322, 231

    work page 2001

  31. [31]

    Krymskii, G. F. 1977, Akademiia Nauk SSSR Doklady, 234, 1306

    work page 1977

  32. [32]

    & Puls, J

    Kudritzki, R.-P. & Puls, J. 2000, ARA&A, 38, 613 LHAASO Collaboration. 2024, Science Bulletin, 69, 449

    work page 2000

  33. [33]

    2023, PhD thesis, Siena U

    Menchiari, S. 2023, PhD thesis, Siena U

    work page 2023

  34. [34]

    Menchiari, S., Morlino, G., Amato, E., Bucciantini, N., & Beltrán, M. T. 2024, A&A, 686, A242

    work page 2024

  35. [35]

    2021, MNRAS, 504, 6096

    Morlino, G., Blasi, P., Peretti, E., & Cristofari, P. 2021, MNRAS, 504, 6096

    work page 2021

  36. [36]

    & de Jager, C

    Nieuwenhuijzen, H. & de Jager, C. 1990, A&A, 231, 134

    work page 1990

  37. [37]

    Padilha, L. N. & Anjos, R. C. 2025, MNRAS, staf2192

    work page 2025

  38. [38]

    2024, ApJ, 972, L22

    Peron, G., Morlino, G., Gabici, S., et al. 2024, ApJ, 972, L22

    work page 2024

  39. [39]

    J., Menten, K

    Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, ApJ, 885, 131

    work page 2019

  40. [40]

    C., Gerhard, O., Portail, M., Vasiliev, E., & Clarke, J

    Sormani, M. C., Gerhard, O., Portail, M., Vasiliev, E., & Clarke, J. 2022, MN- RAS, 514, L1

    work page 2022

  41. [41]

    & Egberts, K

    Steppa, C. & Egberts, K. 2020, A&A, 643, A137

    work page 2020

  42. [42]

    W., Moskalenko, I

    Strong, A. W., Moskalenko, I. V ., & Ptuskin, V . S. 2007, Annual Review of Nuclear and Particle Science, 57, 285

    work page 2007

  43. [43]

    & Reville, B

    Vieu, T. & Reville, B. 2023, MNRAS, 519, 136

    work page 2023

  44. [44]

    1977, ApJ, 218, 377

    Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377

    work page 1977

  45. [45]

    R., van den Heuvel, E

    Yungelson, L. R., van den Heuvel, E. P. J., Vink, J. S., Portegies Zwart, S. F., & de Koter, A. 2008, A&A, 477, 223

    work page 2008

  46. [46]

    2015, Proc

    Zabalza, V . 2015, Proc. of International Cosmic Ray Conference 2015, 922 Article number, page 12 of 12

    work page 2015