pith. sign in

arxiv: 2604.11271 · v1 · submitted 2026-04-13 · 🌌 astro-ph.HE · astro-ph.GA

A Comparative Study of TeV Gamma-Ray Sources with Various Objects

Xin Zhou , Ji Yang , Yang Su , Xuepeng Chen , Yang Chen , Yan Sun , Qing-Zeng Yan , Shaobo Zhang This is my paper

Pith reviewed 2026-05-10 15:31 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords TeV gamma-ray sourcesLHAASOpulsar wind nebulaesupernova remnantsmicroquasarsmolecular cloudssource associationsgamma-ray emission
0
0 comments X

The pith

LHAASO TeV gamma-ray sources show real overlaps with about 20 percent of pulsar wind nebulae and supernova remnants.

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

The authors apply a randomization-adjusted overlap method to test whether LHAASO TeV sources line up with known catalogs of objects more often than chance would allow. They find clear excesses for pulsar wind nebulae, supernova remnants, and microquasars, while OB associations match random expectations. This points to the first two classes as major sites of particle acceleration to TeV energies. The work further shows that pulsar wind nebulae sitting near molecular clouds are brighter and that the gamma-ray positions shift toward those clouds. Roughly 60 percent of pulsar wind nebulae appear bright across both LHAASO energy bands, but only 10 percent of supernova remnants and microquasars do.

Core claim

The paper claims that the RAOC method yields association proportions of 0.19 with SNRs, 0.20 with PWNe, and 0.027 with microquasars for LHAASO sources. PWNe linked to molecular clouds show stronger emission and a systematic positional offset toward the clouds. About 60 percent of PWNe are gamma-ray bright in both WCDA and KM2A ranges, while the fraction is roughly 10 percent for SNRs and microquasars. Shell-type SNRs associate at about 0.1, and 70 percent of the two LHAASO components share a common origin.

What carries the argument

The Randomization-Adjusted Overlap Correlation (RAOC) method, which measures how much observed positional overlaps exceed those in randomized versions of the object catalogs.

Load-bearing premise

The randomization procedure must accurately reproduce the real spatial clustering and selection biases present in the SNR, PWN, and microquasar catalogs.

What would settle it

A new high-resolution TeV survey that measures no average positional offset between PWN-associated sources and their nearby molecular clouds would falsify the claimed role of clouds in enhancing emission.

Figures

Figures reproduced from arXiv: 2604.11271 by Ji Yang, Qing-Zeng Yan, Shaobo Zhang, Xin Zhou, Xuepeng Chen, Yang Chen, Yang Su, Yan Sun.

Figure 1
Figure 1. Figure 1: Left: Example case of two types of samples (two triangles and three circles) dis [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Top: l-b map of the locations of LHAASO TeV gamma-ray sources (red circles) and known SNRs (excluding the PWN sample; black pluses). An example of randomly distributed SNRs (blue crosses) and the adaptive refined mesh grid used to generate them are also shown (see Section 3 for details). The bottom left and middle panels show the l and b distributions of the SNRs (black) and of an example of randomly distr… view at source ↗
Figure 3
Figure 3. Figure 3: 12CO (J=1–0) intensity maps of the candidates of the associated MCs around PWNe G18.00−0.69, Eel, G21.88−0.10, G23.5+0.1, G25.24−0.19, G32.64+0.53, G36.01+0.06, and G63.7+1.1 (see Zhou et al. 2023, for details). The intensities of the 12CO emission are at least 3σ. The extents of the PWNe are indicated by red circles. The surrounding shell-like structures are outlined by black solid circles, except for PWN… view at source ↗
Figure 4
Figure 4. Figure 4: Same as Figure 3, but for PWNe CTB 87, G75.23 [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
read the original abstract

We investigate the relationships between LHAASO TeV gamma-ray sources and various kinds of objects, including pulsar wind nebulae (PWNe), supernova remnants (SNRs), HII regions, microquasars, and OB associations. We propose a Randomization-Adjusted Overlap Correlation (RAOC) method to statistically assess association probabilities and evaluate association proportions across catalogs. The results reveal statistically significant overlaps between LHAASO sources and SNRs, PWNe, and microquasars, supporting their role as important contributors to TeV gamma-ray emission. The estimated association proportions of LHAASO sources are 0.19$\pm$0.08 with SNRs, 0.20$\pm$0.04 with PWNe, and 0.027$\pm$0.008 with microquasars. The proportion of the gamma-ray sources associated with the subsample of shell-type SNRs is ~0.1. While HII regions also show potential association, particularly with the KM2A component, their large self-overlap ratio complicates precise estimation. In contrast, OB associations exhibit a high probability of chance coincidence, suggesting their limited contribution to TeV gamma-ray emission. Our analysis of TeV gamma-ray emission capabilities shows that ~60% of PWNe are gamma-ray bright in both the WCDA and KM2A energy ranges. For SNRs and microquasars, the TeV gamma-ray bright fraction is ~10%. The subsample of PWNe associated with molecular clouds (MCs) shows enhanced gamma-ray emission. Furthermore, positional analysis reveals a systematic offset of the gamma-ray sources overlapping with PWNe toward the associated MCs. These findings imply a role for MCs in PWN gamma-ray production. Additionally, self-correlation analysis indicates that about 70% of the WCDA and KM2A gamma-ray components share a common origin. The study also identifies selection effects in existing SNR catalogs and notes clustering among approximately 30% of HII regions within larger star-forming regions.

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 paper introduces the Randomization-Adjusted Overlap Correlation (RAOC) method to quantify associations between LHAASO TeV gamma-ray sources and catalogs of PWNe, SNRs, HII regions, microquasars, and OB associations. It reports statistically significant overlaps yielding association proportions of 0.19±0.08 (SNRs), 0.20±0.04 (PWNe), and 0.027±0.008 (microquasars), with ~60% of PWNe and ~10% of SNRs/microquasars being TeV-bright; it further notes enhanced emission and positional offsets toward molecular clouds for PWNe, complications from HII self-overlap, high chance-coincidence for OB associations, ~70% common origin between WCDA and KM2A components, and selection effects in SNR catalogs.

Significance. If the RAOC null distributions are shown to faithfully reproduce the observed spatial clustering, latitude distributions, and incompleteness of the input catalogs, the quantitative association fractions and the MC-PWN offset analysis would provide useful constraints on the Galactic TeV source population and the role of molecular clouds in particle acceleration. The self-correlation result between WCDA and KM2A components offers an internal consistency check that strengthens the data interpretation.

major comments (2)
  1. [Methods (RAOC)] Methods section on RAOC: the randomization procedure is described only at a high level ('statistically assess association probabilities') without specifying the kernel (e.g., whether Monte-Carlo positions are drawn from the observed Galactic density field, whether catalog angular sizes or flux limits are preserved, or the number of realizations). Because the reported proportions and significance statements rest entirely on the resulting null distribution, this omission is load-bearing.
  2. [Results (association proportions)] Results, association-proportion paragraph: the quoted uncertainties (0.19±0.08 for SNRs, 0.20±0.04 for PWNe) are presented as if they fully incorporate both Poisson counting and RAOC variance, yet no table or equation shows how the two contributions are combined or how catalog incompleteness propagates. This directly affects the claim of 'statistically significant overlaps'.
minor comments (2)
  1. [Abstract and Results] The abstract states that HII regions have a 'large self-overlap ratio' that 'complicates precise estimation,' but neither the value of this ratio nor the exact correction applied in RAOC is given in the text or tables.
  2. [Results (PWN-MC offsets)] The positional-offset analysis toward MCs for PWNe is presented without quoted uncertainties on the mean offset or a control sample of non-associated sources; adding these would clarify the robustness of the claimed systematic shift.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful comments, which have helped us improve the clarity and rigor of our manuscript. We address each major comment below and have made revisions to the paper as indicated.

read point-by-point responses

  1. Referee: [Methods (RAOC)] Methods section on RAOC: the randomization procedure is described only at a high level ('statistically assess association probabilities') without specifying the kernel (e.g., whether Monte-Carlo positions are drawn from the observed Galactic density field, whether catalog angular sizes or flux limits are preserved, or the number of realizations). Because the reported proportions and significance statements rest entirely on the resulting null distribution, this omission is load-bearing.

    Authors: We agree that the description of the RAOC method in the original manuscript was at a high level and lacked sufficient detail on the randomization kernel. This is a valid point, and we have revised the Methods section to provide a complete specification. Specifically, we now describe that the Monte Carlo positions are sampled from the observed Galactic longitude and latitude distributions of each catalog to preserve the spatial clustering and incompleteness properties. Catalog angular sizes are preserved by assigning the same extent to randomized objects, and flux limits are maintained by using the same selection criteria. We have also included a new figure in the supplementary material demonstrating that the null distributions accurately reproduce the observed latitude distributions and clustering statistics of the input catalogs. revision: yes

  2. Referee: [Results (association proportions)] Results, association-proportion paragraph: the quoted uncertainties (0.19±0.08 for SNRs, 0.20±0.04 for PWNe) are presented as if they fully incorporate both Poisson counting and RAOC variance, yet no table or equation shows how the two contributions are combined or how catalog incompleteness propagates. This directly affects the claim of 'statistically significant overlaps'.

    Authors: The referee correctly identifies that the uncertainty estimation was not fully transparent. The reported uncertainties primarily reflect the standard deviation from the RAOC null distributions, which incorporate the effects of spatial distributions and catalog properties. However, we did not explicitly detail the combination with Poisson counting statistics or the propagation of incompleteness. In the revised manuscript, we have added a dedicated paragraph and an equation in the Results section that explains the error budget: the total uncertainty is calculated as the root-sum-square of the RAOC variance and the Poisson uncertainty. We also discuss how incompleteness in the SNR and PWN catalogs is mitigated by using subsamples with known completeness and note that this affects the association proportions at the level of the quoted uncertainties. A table has been added to break down the uncertainty components for each association proportion. revision: yes

Circularity Check

0 steps flagged

No circularity: RAOC yields measured proportions from external catalog overlaps

full rationale

The paper defines RAOC as a randomization procedure applied to observed positions in independent catalogs (LHAASO sources vs. SNRs/PWNe/microquasars). The reported association proportions (0.19±0.08, 0.20±0.04, 0.027±0.008) and significance statements are direct numerical outputs of that procedure, not inputs or self-referential definitions. No equation reduces a claimed prediction to a fitted parameter by construction, no uniqueness theorem is imported from prior self-work, and no ansatz is smuggled via citation. The derivation chain is self-contained against the input catalogs and does not collapse to tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard catalog cross-matching assumptions and the validity of the RAOC randomization procedure; no new physical entities are introduced.

axioms (1)
  • domain assumption Spatial distributions of catalog objects can be adequately modeled by randomization for chance-coincidence estimation
    Invoked in the definition of RAOC to subtract random overlaps

pith-pipeline@v0.9.0 · 5697 in / 1239 out tokens · 39726 ms · 2026-05-10T15:31:12.650643+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

18 extracted references · 18 canonical work pages

  1. [1]

    D., Bania, T

    Anderson, L. D., Bania, T. M., Balser, D. S., et al. 2014, ApJS, 212, 1

    work page 2014

  2. [2]

    2013, A&A Rev., 21, 70

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

    work page 2013

  3. [3]

    X., et al

    Cao, Z., Aharonian, F., Bai, Y. X., et al. 2025a, arXiv e-prints, arXiv:2510.26112

    work page arXiv

  4. [4]

    A., de Grijs, R., Glushkova, E

    Chemel, A. A., de Grijs, R., Glushkova, E. V., & Dambis, A. K. 2022, MNRAS, 515, 4359

    work page 2022

  5. [5]

    2012, Advances in Space Research, 49, 1313

    Ferrand, G., & Safi-Harb, S. 2012, Advances in Space Research, 49, 1313

    work page 2012

  6. [6]

    Green, D. A. 2019, Journal of Astrophysics and Astronomy, 40, 36

    work page 2019

  7. [7]

    J., Rowell, G., Hofmann, W., et al

    Hampton, E. J., Rowell, G., Hofmann, W., et al. 2016, Journal of High Energy Astrophysics, 11, 1

    work page 2016

  8. [8]

    2025, ApJ, 990, 142

    Ji, S., Wang, Z., Zheng, D., & Zheng, J. 2025, ApJ, 990, 142

    work page 2025

  9. [9]

    2017, ApJ, 851, 37

    Liu, B., Chen, Y., Zhang, X., et al. 2017, ApJ, 851, 37

    work page 2017

  10. [10]

    F., & Rodríguez, L

    Mirabel, I. F., & Rodríguez, L. F. 1999, ARA&A, 37, 409

    work page 1999

  11. [11]

    Mitchell, A. M. W., & Celli, S. 2024, Journal of High Energy Astrophysics, 44, 340

    work page 2024

  12. [12]

    A., & McClintock, J

    Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49

    work page 2006

  13. [13]

    2025, ApJ, 980, 98

    Sun, J., Chen, Y., Bao, Y., Zhang, X., & Zhou, X. 2025, ApJ, 980, 98

    work page 2025

  14. [14]

    2026, ApJS, 282, 65

    Yang, J., Yan, Q.-Z., Su, Y., et al. 2026, ApJS, 282, 65

    work page 2026

  15. [15]

    2022, ApJ, 931, 128

    Zhang, X., Chen, Y., Zheng, F.-x., et al. 2022, ApJ, 931, 128

    work page 2022

  16. [16]

    2023, MNRAS, 521, 1931

    Zhong, W.-J., Zhang, X., Chen, Y., & Zhang, Q.-Q. 2023, MNRAS, 521, 1931

    work page 2023

  17. [17]

    2024, A&A, 683, A107

    Zhou, X., Su, Y., Yang, J., Chen, Y., & Jiang, Z. 2024, A&A, 683, A107

    work page 2024

  18. [18]

    2023, ApJS, 268, 61

    Zhou, X., Su, Y., Yang, J., et al. 2023, ApJS, 268, 61

    work page 2023