pith. sign in

arxiv: 2603.04821 · v2 · pith:XOPUUCCDnew · submitted 2026-03-05 · 🌌 astro-ph.HE · astro-ph.IM

Resolving diffusion signatures in distant pulsar halos with current and future experiments

Yong-Jian Wei , En-Sheng Chen , Kun Fang , Xiao-Jun Bi This is my paper

Pith reviewed 2026-05-21 12:06 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IM
keywords pulsar halosgamma-ray observationscosmic-ray propagationmorphological discriminationLHAASOCTAdiffusion modelsinterstellar medium
0
0 comments X

The pith

Simulations show that CTA and LHAASO-KM2A can distinguish diffusion patterns in pulsar halos from simpler models through better resolution and higher photon counts.

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

Pulsar halos are extended gamma-ray regions around pulsars produced by diffusing cosmic rays, offering a direct look at how particles move through the interstellar medium near these sources. Confirming candidates has been hard because current instruments struggle to separate the specific shapes expected from diffusion from other possible spatial distributions. This paper runs simulations of mock observations for LHAASO-KM2A and the Cherenkov Telescope Array that include realistic detector performance. The work finds that collecting more photons and achieving sharper angular resolution both increase the ability to tell diffusion-based halos apart from simplified alternatives. If these capabilities hold in real data, the number of confirmed pulsar halos should rise and give clearer pictures of cosmic-ray behavior close to pulsars.

Core claim

Using mock observations with realistic instrumental responses, the ability of LHAASO-KM2A and CTA to distinguish diffusion-based halo morphologies from alternative simplified spatial models is significantly enhanced by increased photon statistics and improved angular resolution, with CTA gaining from its superior angular resolution and LHAASO-KM2A from its large effective area at the highest energies.

What carries the argument

Morphological discrimination between diffusion-based pulsar halo models and simplified alternative spatial models, evaluated through simulated observations that incorporate each experiment's instrumental responses.

If this is right

  • CTA observations will use superior angular resolution to strengthen morphological separation of halo candidates at lower energies.
  • LHAASO-KM2A will exploit its large effective area at the highest energies to increase sensitivity for halo morphology studies.
  • The total number of firmly identified pulsar halos is expected to grow beyond the current limited sample.
  • Confirmed halos will supply tighter constraints on cosmic-ray diffusion in the immediate vicinity of pulsars.

Where Pith is reading between the lines

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

  • These discrimination techniques could be applied to other extended high-energy sources to test diffusion versus alternative emission geometries.
  • Comparing simulated discrimination power against actual future data sets would provide a direct test of how well current diffusion models match observations.
  • Success in expanding the halo sample might allow statistical studies of how diffusion length scales with pulsar age or spin-down power.

Load-bearing premise

The mock observations with realistic instrumental responses accurately capture the differences between diffusion-based halo morphologies and the alternative simplified spatial models that would be present in real data.

What would settle it

If real observations of pulsar halo candidates with LHAASO-KM2A or CTA fail to show the expected statistical preference for diffusion morphologies over simplified models at the levels indicated by the simulations, the discrimination improvement would not hold.

Figures

Figures reproduced from arXiv: 2603.04821 by En-sheng Chen, Kun Fang, Xiao-Jun Bi, Yong-Jian Wei.

Figure 1
Figure 1. Figure 1: FIG. 1. Radial [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Discrimination curves for pulsar halos achievable with CTA-North and LHAASO-KM2A. [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
read the original abstract

Pulsar halos provide a unique probe of cosmic-ray propagation in the vicinity of pulsars and have important implications for our understanding of particle diffusion in the interstellar medium. However, the number of firmly identified pulsar halos remains limited. One of the main challenges is the difficulty in unambiguously confirming halo candidates through precise morphological measurements with current $\gamma$-ray observations. In this work, we investigate the prospects for identifying pulsar halo candidates through morphological discrimination using simulations of two advanced $\gamma$-ray experiments: LHAASO-KM2A and the Cherenkov Telescope Array (CTA). Using mock observations with realistic instrumental responses, we assess the ability of each experiment to distinguish diffusion-based halo morphologies from alternative simplified spatial models. Our results show that both increased photon statistics and improved angular resolution significantly enhance the power of morphological discrimination. In particular, CTA benefits from its superior angular resolution, while LHAASO-KM2A gains sensitivity from its large effective area at the highest energies. These results indicate that future $\gamma$-ray observations have the potential to expand the sample of pulsar halos and provide further insights into cosmic-ray transport around pulsars.

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 claims that forward simulations of mock gamma-ray observations of distant pulsar halos, incorporating realistic instrumental responses for LHAASO-KM2A and CTA, demonstrate that higher photon statistics and better angular resolution substantially improve the ability to morphologically discriminate diffusion-based halo profiles from simplified alternative spatial models such as Gaussians or power laws.

Significance. If the central results hold, the work would provide actionable guidance for expanding the sample of confirmed pulsar halos and constraining cosmic-ray diffusion in the ISM near pulsars. The forward-simulation approach with realistic responses is a strength, as it yields concrete, experiment-specific predictions rather than purely theoretical forecasts.

major comments (2)
  1. [§3] §3 (Mock Observation Setup): The central discrimination result relies on the chosen alternative models (Gaussian, power-law, etc.) spanning the range of plausible non-diffusion morphologies. The manuscript should demonstrate that the reported gains in discrimination power persist when the mocks incorporate energy-dependent diffusion coefficients or spatially varying ISM density fluctuations, as these effects could produce morphologies that overlap more with the diffusion case.
  2. [§4] §4 (Results on Morphological Discrimination): The quantitative improvement from photon statistics and angular resolution is load-bearing for the conclusions about CTA and LHAASO-KM2A. The paper must report the exact statistical metric (e.g., likelihood ratio or Bayes factor) used to quantify discrimination and show that it remains robust under reasonable variations in the diffusion index.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'realistic instrumental responses' would benefit from a one-sentence clarification of the key components included (e.g., PSF, energy resolution, background model) to aid readers who do not reach the methods section.
  2. [Figures] Figure captions: Several result figures would be clearer if they explicitly labeled the diffusion index value and the exact alternative model parameters used in each panel.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive review and positive assessment of the forward-simulation approach. We address the two major comments point by point below. Revisions have been made to improve clarity and robustness where feasible.

read point-by-point responses

  1. Referee: §3 (Mock Observation Setup): The central discrimination result relies on the chosen alternative models (Gaussian, power-law, etc.) spanning the range of plausible non-diffusion morphologies. The manuscript should demonstrate that the reported gains in discrimination power persist when the mocks incorporate energy-dependent diffusion coefficients or spatially varying ISM density fluctuations, as these effects could produce morphologies that overlap more with the diffusion case.

    Authors: We appreciate the referee's point that more physically motivated alternatives could reduce the separation from the diffusion morphology. Our baseline alternatives were chosen to match the simplified profiles routinely fitted to current data. A full re-simulation suite that self-consistently includes energy-dependent diffusion coefficients and spatially varying ISM fluctuations would require a substantial expansion of the computational campaign and is beyond the scope of the present study. In the revised manuscript we have added a dedicated paragraph in §3 that discusses this limitation, explains why the reported gains in discrimination power are expected to be conservative (because diffusion introduces energy-dependent morphology absent from static alternatives), and notes that future work should explore these extensions. No new quantitative results are claimed for the complex cases. revision: partial

  2. Referee: §4 (Results on Morphological Discrimination): The quantitative improvement from photon statistics and angular resolution is load-bearing for the conclusions about CTA and LHAASO-KM2A. The paper must report the exact statistical metric (e.g., likelihood ratio or Bayes factor) used to quantify discrimination and show that it remains robust under reasonable variations in the diffusion index.

    Authors: We agree that the precise statistical measure must be stated explicitly. The discrimination power was quantified via a likelihood-ratio test between the diffusion model and each alternative spatial template; the test statistic and its distribution under the null are now written out in §4. To demonstrate robustness, we have added an appendix that repeats the full analysis for diffusion indices varied by ±0.2 around the fiducial value. The improvement in separation significance with increased photon statistics and better angular resolution persists across this range, with only modest quantitative shifts. The corresponding figures and tables have been included in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: forward simulation of instrumental discrimination power

full rationale

The paper conducts forward modeling of mock gamma-ray observations for LHAASO-KM2A and CTA, folding in realistic instrumental responses to quantify the statistical ability to separate diffusion-based halo morphologies from simplified alternatives such as Gaussian or power-law profiles. This process generates new synthetic data and applies discrimination metrics; it does not derive any quantity that reduces by construction to a fitted parameter, self-citation, or input ansatz. The central results on gains from photon statistics and angular resolution emerge directly from the simulation outputs without re-labeling or self-referential closure. The analysis is self-contained against external benchmarks of instrument performance and does not invoke load-bearing uniqueness theorems or prior author work as the sole justification for its conclusions.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. Simulations presumably rely on standard instrument response functions and a diffusion model whose details are not provided.

pith-pipeline@v0.9.0 · 5741 in / 1059 out tokens · 57618 ms · 2026-05-21T12:06:39.571782+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

29 extracted references · 29 canonical work pages · 6 internal anchors

  1. [1]

    Therefore, we employ the AIC as a model selection metric, which is widely used in data analysis for comparing non-nested models

    AIC: Since the diffusion model and the alternative models are non-nested, the likeli- hood ratio test and Wilks’ theorem cannot be directly applied to assess the statistical significance of their differences. Therefore, we employ the AIC as a model selection metric, which is widely used in data analysis for comparing non-nested models. The AIC is defined ...

    work page

  2. [2]

    A. A. Abdoet al., Astrophys. J. Lett.664, L91 (2007), arXiv:0705.0707 [astro-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2007

  3. [3]

    P. A. Caraveo, G. F. Bignami, A. De Luca, S. Mereghetti, A. Pellizzoni, R. Mignani, A. Tur, and W. Becker, Science301, 1345 (2003)

    work page 2003

  4. [4]

    Geminga's puzzling pulsar wind nebula

    B. Posselt, G. Pavlov, P. Slane, R. Romani, N. Bucciantini, A. Bykov, O. Kargaltsev, M. Weis- skopf, and C. Y. Ng, Astrophys. J.835, 66 (2017), arXiv:1611.03496 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  5. [5]

    A. U. Abeysekaraet al.(HAWC), Science358, 911 (2017), arXiv:1711.06223 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  6. [6]

    Fang, Front

    K. Fang, Front. Astron. Space Sci.9, 1022100 (2022), arXiv:2209.13294 [astro-ph.HE]

    work page arXiv 2022

  7. [7]

    Liu, Int

    R.-Y. Liu, Int. J. Mod. Phys. A37, 2230011 (2022), arXiv:2207.04011 [astro-ph.HE]

    work page arXiv 2022

  8. [8]

    L´ opez-Coto, E

    R. L´ opez-Coto, E. de O˜ na Wilhelmi, F. Aharonian, E. Amato, and J. Hinton, Nature Astron. 6, 199 (2022), arXiv:2202.06899 [astro-ph.HE]

    work page arXiv 2022

  9. [9]

    Amato and S

    E. Amato and S. Recchia, Riv. Nuovo Cim.47, 399 (2024), arXiv:2409.00659 [astro-ph.HE]

    work page arXiv 2024

  10. [10]

    Albertet al.(HAWC), 3HWC: The Third HAWC Cat- alog of Very-High-Energy Gamma-ray Sources, Astro- phys

    A. Albertet al.(HAWC), Astrophys. J.905, 76 (2020), arXiv:2007.08582 [astro-ph.HE]

    work page arXiv 2020

  11. [11]

    Caoet al.(LHAASO), Astrophys

    Z. Caoet al.(LHAASO), Astrophys. J. Suppl.271, 25 (2024), arXiv:2305.17030 [astro-ph.HE]

    work page arXiv 2024

  12. [13]

    Design Concepts for the Cherenkov Telescope Array

    M. Actiset al.(CTA Consortium), Exper. Astron.32, 193 (2011), arXiv:1008.3703 [astro- ph.IM]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  13. [14]

    Schwefer, R

    G. Schwefer, R. Parsons, and J. Hinton, Astropart. Phys.163, 103008 (2024), arXiv:2406.17502 [astro-ph.HE]

    work page arXiv 2024

  14. [15]

    B. M. Gaensler and P. O. Slane, Ann. Rev. Astron. Astrophys.44, 17 (2006), arXiv:astro- ph/0601081

    work page arXiv 2006

  15. [16]

    Albertet al.(HAWC), Astrophys

    A. Albertet al.(HAWC), Astrophys. J.974, 246 (2024)

    work page 2024

  16. [17]

    B. S. Acharyaet al.(CTA Consortium),Science with the Cherenkov Telescope Array(WSP,

    work page

  17. [18]

    arXiv:1709.07997 [astro-ph.IM]

    work page internal anchor Pith review Pith/arXiv arXiv

  18. [19]

    Ctao instrument response functions - prod5 version v0.1,

    C. T. A. Observatory and C. T. A. Consortium, “Ctao instrument response functions - prod5 version v0.1,” (2021)

    work page 2021

  19. [20]

    Celli and G

    S. Celli and G. Peron, Astron. Astrophys.689, A258 (2024), arXiv:2403.03731 [astro-ph.HE]

    work page arXiv 2024

  20. [21]

    X. Bai, B. Y. Bi, X. J. Bi, Z. Cao, S. Z. Chen, Y. Chen, A. Chiavassa, X. H. Cui, Z. G. 14 Dai, D. della Volpe, T. Di Girolamo, G. Di Sciascio, Y. Z. Fan, J. Giacalone, Y. Q. Guo, H. H. He, T. L. He, M. Heller, D. Huang, Y. F. Huang, H. Jia, L. T. Ksenofontov, D. Leahy, F. Li, Z. Li, E. W. Liang, P. Lipari, R. Y. Liu, Y. Liu, S. Liu, X. Ma, O. Martineau-H...

    work page arXiv 2019

  21. [22]

    Caoet al.(LHAASO), Chin

    Z. Caoet al.(LHAASO), Chin. Phys. C48, 065001 (2024), arXiv:2401.01038 [astro-ph.IM]

    work page arXiv 2024

  22. [23]

    T. K. Gaisser, T. Stanev, and S. Tilav, Front. Phys. (Beijing)8, 748 (2013), arXiv:1303.3565 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  23. [24]

    Aharonianet al.(LHAASO), Extended Very- High-Energy Gamma-Ray Emission Surrounding PSR J0622+3749 Observed by LHAASO-KM2A, Phys

    F. Aharonianet al.(LHAASO), Phys. Rev. Lett.126, 241103 (2021), arXiv:2106.09396 [astro- ph.HE]

    work page arXiv 2021

  24. [25]

    Caoet al.(LHAASO), Sci

    Z. Caoet al.(LHAASO), Sci. China Phys. Mech. Astron.68, 279504 (2025), arXiv:2410.04425 [astro-ph.HE]

    work page arXiv 2025

  25. [26]

    Fang, S.-Q

    K. Fang, S.-Q. Xi, L.-Z. Bao, X.-J. Bi, and E.-S. Chen, Phys. Rev. D106, 123017 (2022), arXiv:2207.13533 [astro-ph.HE]

    work page arXiv 2022

  26. [27]

    Sabri, Y

    K. Sabri, Y. Gallant, J. Devin, and K. Feijen, arXiv e-prints , arXiv:2510.03183 (2025), arXiv:2510.03183 [astro-ph.HE]

    work page arXiv 2025

  27. [28]

    Albertet al.(HAWC), Astrophys

    A. Albertet al.(HAWC), Astrophys. J. Lett.944, L29 (2023), arXiv:2301.04646 [astro-ph.HE]

    work page arXiv 2023

  28. [29]

    Wach and A

    T. Wach and A. M. W. Mitchell, arXiv e-prints , arXiv:2510.02802 (2025), arXiv:2510.02802 [astro-ph.HE]

    work page arXiv 2025

  29. [30]

    De Sarkar, arXiv e-prints , arXiv:2601.21689 (2026), arXiv:2601.21689 [astro-ph.HE]

    A. De Sarkar, arXiv e-prints , arXiv:2601.21689 (2026), arXiv:2601.21689 [astro-ph.HE]. 15

    work page arXiv 2026