pith. the verified trust layer for science. sign in

arxiv: 2508.13667 · v2 · submitted 2025-08-19 · 🌌 astro-ph.HE · astro-ph.IM

PHECT: A lightweight computation tool for pulsar halo emission

Pith reviewed 2026-05-18 22:59 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IM
keywords pulsar halosgamma-ray emissioncosmic-ray propagationtransport modelingfinite-volume methodsinterstellar mediumnumerical tools
0
0 comments X p. Extension

The pith

PHECT is a lightweight tool that models pulsar halo gamma-ray emission with multiple transport models and stable finite-volume methods on irregular grids.

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

The paper presents PHECT, a software package for computing the gamma-ray output of pulsar halos produced when electrons and positrons scatter off background photons in the surrounding interstellar medium. It moves beyond the common assumption of simple homogeneous diffusion by letting users select from several more complex transport descriptions. All settings are supplied in a YAML file, and the underlying solver uses finite-volume discretizations that remain stable even when the spatial grid is uneven or the diffusion coefficient changes abruptly. The design targets the coming era of higher-resolution observations and the growing catalog of detected halos, providing a practical way to test how particles actually propagate on scales of tens of parsecs.

Core claim

PHECT enables modeling of pulsar halo emission by incorporating multiple transport models extending beyond standard diffusion and adopting finite-volume discretizations that remain stable on non-uniform grids and in the presence of discontinuous diffusion coefficients, with all computations driven by a YAML configuration file and no manual code changes required.

What carries the argument

Finite-volume discretization of the particle transport equations, which maintains numerical stability on non-uniform spatial grids and across jumps in the diffusion coefficient.

If this is right

  • Users can switch among different possible origins of pulsar halos simply by editing the configuration file.
  • The numerical method supports flexible, non-uniform grids that concentrate resolution where needed near the pulsar.
  • Discontinuous diffusion coefficients can be handled without solver breakdown, allowing more realistic medium properties.
  • The lightweight, YAML-driven interface makes the tool practical for analyzing the growing sample of observed halos with higher-precision data.

Where Pith is reading between the lines

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

  • The same stable discretization could be applied to other localized cosmic-ray sources that involve sharp changes in propagation conditions.
  • Integration with larger-scale Galactic propagation codes could test whether halo-scale data are consistent with global models.
  • Public release of the YAML interface may encourage community comparisons of transport assumptions against the same observational targets.

Load-bearing premise

The finite-volume scheme stays stable and produces reliable results when applied to the transport models and grid types used for pulsar halo calculations.

What would settle it

A side-by-side comparison in which PHECT output for standard diffusion on a uniform grid deviates substantially from an analytic solution or from an independent established code would show the discretization is unsuitable.

Figures

Figures reproduced from arXiv: 2508.13667 by Kun Fang.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic illustrations of two of the possible slow-diffusion interpretations. The star [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. One-dimensional [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Two-dimensional [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Average of [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
read the original abstract

$\gamma$-ray pulsar halos, most likely formed by inverse Compton scattering of electrons and positrons propagating in the pulsar-surrounding interstellar medium with background photons, serve as an ideal probe for Galactic cosmic-ray propagation on small scales (typically tens of parsecs). While the associated electron and positron propagation is often modeled using homogeneous and isotropic diffusion, termed here as standard diffusion, the actual transport process is expected to be more complex. This work introduces the Pulsar Halo Emission Computation Tool (PHECT), a lightweight software designed for modeling pulsar halo emission. PHECT incorporates multiple transport models extending beyond standard diffusion, accounting for different possible origins of pulsar halos. Users can conduct necessary computations simply by configuring a YAML file without manual code edits. Furthermore, the tool adopts finite-volume discretizations that remain stable on non-uniform grids and in the presence of discontinuous diffusion coefficients. PHECT is ready for the increasingly precise observational data and the rapidly growing sample of pulsar halos.

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

1 major / 2 minor

Summary. The manuscript introduces PHECT, a lightweight software tool for modeling γ-ray emission from pulsar halos formed by inverse Compton scattering of electrons and positrons. It supports multiple transport models that extend beyond the standard homogeneous isotropic diffusion approximation, employs finite-volume discretizations asserted to remain stable on non-uniform grids and with discontinuous diffusion coefficients, and is configured entirely through a YAML file without requiring code modifications. The tool is presented as ready for use with increasingly precise observational data and the growing sample of detected pulsar halos.

Significance. If the numerical implementation delivers the claimed stability and flexibility, PHECT would provide a practical, accessible resource for exploring non-standard cosmic-ray transport scenarios on tens-of-parsec scales. The YAML-driven configuration and support for multiple models address usability needs in the field as HAWC, LHAASO, and similar instruments deliver higher-quality halo spectra.

major comments (1)
  1. [Numerical methods section] Section describing the numerical methods (likely §3 or equivalent): the central claim that the finite-volume discretizations 'remain stable on non-uniform grids and in the presence of discontinuous diffusion coefficients' is presented without any benchmark tests, convergence studies, error analysis, or comparisons against known solutions. This demonstration is load-bearing for the tool's asserted suitability for pulsar-halo modeling and must be supplied before the stability assertion can be evaluated.
minor comments (2)
  1. [Abstract and §1] The abstract and introduction would benefit from a concise list or table of the specific transport models implemented beyond standard diffusion.
  2. [User interface / configuration section] An example YAML configuration file (or excerpt) should be included to illustrate the claimed 'no manual code edits' workflow.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the need for explicit validation of the numerical methods. We address the major comment below and will revise the manuscript to incorporate the requested demonstrations.

read point-by-point responses
  1. Referee: [Numerical methods section] Section describing the numerical methods (likely §3 or equivalent): the central claim that the finite-volume discretizations 'remain stable on non-uniform grids and in the presence of discontinuous diffusion coefficients' is presented without any benchmark tests, convergence studies, error analysis, or comparisons against known solutions. This demonstration is load-bearing for the tool's asserted suitability for pulsar-halo modeling and must be supplied before the stability assertion can be evaluated.

    Authors: We agree that the stability claims for the finite-volume discretizations require explicit numerical validation to be fully substantiated. The manuscript introduces PHECT as a practical tool and emphasizes its configuration flexibility, but we acknowledge that the absence of benchmark tests, convergence studies, and error analysis weakens the presentation of the numerical core. In the revised manuscript we will expand the numerical methods section with a new subsection containing benchmark tests against known analytical solutions for the diffusion equation. These will cover uniform and non-uniform grids as well as cases with discontinuous diffusion coefficients, together with convergence rates and L2-error analyses. The added material will directly support the tool's claimed suitability for pulsar-halo modeling while preserving the lightweight character of the code. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The manuscript introduces PHECT as a configurable software tool implementing multiple transport models and finite-volume discretizations for pulsar halo emission modeling. No derivation chain, first-principles predictions, or fitted parameters are presented that reduce by construction to the paper's own inputs or self-citations. The work is self-contained as a description of code architecture and usage via YAML configuration, with stability claims tied directly to the chosen numerical methods rather than any circular redefinition or imported uniqueness theorem.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard domain assumptions about particle transport and emission processes rather than introducing new free parameters or entities.

axioms (1)
  • domain assumption Pulsar halos are formed by inverse Compton scattering of electrons and positrons propagating in the interstellar medium with background photons.
    This physical picture underpins the emission modeling performed by the tool.

pith-pipeline@v0.9.0 · 8769 in / 1056 out tokens · 75926 ms · 2026-05-18T22:59:10.552534+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

81 extracted references · 81 canonical work pages · 20 internal anchors

  1. [1]

    For electrons with energies Ee >∼ 1 GeV, the primary energy-loss processes during prop- agation are synchrotron radiation in magnetic fields and ICS with background photons

    Common framework The electron propagation in the ISM is typically described by the diffusion-loss equation, written as: ∂N ∂t = ∇ · (D∇N) + ∂(bN) ∂Ee + Q , (1) where Ee is the electron kinetic energy, N = N(Ee, r, t) is the energy differential number density of electrons, b = b(Ee) ≡ | dEe/dt| is the absolute energy-loss rate, Q = Q(Ee, r, t) is the sourc...

  2. [2]

    In addition to the physically motivated models discussed in Sec

    Model-specific modules The differences between the models in PHECT primarily manifest in the procedures for solving the electron propagation equation. In addition to the physically motivated models discussed in Sec. II, PHECT also includes phenomenological models such as normal diffusion and spherically symmetric two-zone diffusion. These models are valua...

  3. [3]

    (1) using the Green’s function method is given by N(Ee, r, t) = Z R3 d3r0 Z dEe,0 Z dt0 Q(Ee,0, r0, t0) G(Ee, r, t; Ee,0, r0, t0)

    normal and normEx models The general solution to Eq. (1) using the Green’s function method is given by N(Ee, r, t) = Z R3 d3r0 Z dEe,0 Z dt0 Q(Ee,0, r0, t0) G(Ee, r, t; Ee,0, r0, t0) . (B1) The normal and normEx models are based on homogeneous and isotropic diffusion, corre- sponding to the Green’s function: G(Ee, r, t; Ee,0, r0, t0) = 1 b(Ee)(4πλ)3/2 exp...

  4. [4]

    Given that Drr = M 4 ADzz, if we define a new coordinate z′ such that z′ = M 2 Az, the diffusion equation becomes isotropic in the r − z′ coordinate system

    aniso model Anisotropic diffusion can be solved in cylindrical coordinates through a coordinate trans- formation [56]. Given that Drr = M 4 ADzz, if we define a new coordinate z′ such that z′ = M 2 Az, the diffusion equation becomes isotropic in the r − z′ coordinate system. Based on Eqs. (B1) and (B2), it is straightforward to derive its solution in cyli...

  5. [5]

    [72]): G(Ee, r, t; Ee,0, r0, t0) = ρ(α) 3 (|r − r0|λ−1/α) b(Ee)λ3/α δ(t − t0 − τ)H(τ) , (B14) where ρ(α) 3 (r) = 1 2π2r Z ∞ 0 exp (−kα) sin(kr)kdk

    superdiff model For the fractional diffusion equation, we can extend the Green’s function as follows (see, e.g., Ref. [72]): G(Ee, r, t; Ee,0, r0, t0) = ρ(α) 3 (|r − r0|λ−1/α) b(Ee)λ3/α δ(t − t0 − τ)H(τ) , (B14) where ρ(α) 3 (r) = 1 2π2r Z ∞ 0 exp (−kα) sin(kr)kdk . (B15) The solution for the superdiff model is then N(Ee, r, ts) = Z ts tmin dt0 qt(t0) qE(...

  6. [6]

    III A 2, J¨ utnner introduced a relativistic correction to the Maxwell–Boltzmann distribution

    normJ model As mentioned in Sec. III A 2, J¨ utnner introduced a relativistic correction to the Maxwell–Boltzmann distribution. Due to the formal similarity between the solution kernel of the diffusion equa- tion and the Maxwell–Boltzmann distribution, Ref. [49] employed an analogous method to address the superluminal problem in the diffusion equation. Fo...

  7. [7]

    Spherically symmetric diffusion We divide the space into spherical shells with inner and outer radii of rj−1/2 and rj+1/2, respectively, where j = 0 , 1, 2, ..., J9, referred to as control volumes (CVs). The principle of FVM is to integrate the equation over each CV and transform the volume integral into the flux difference across the interfaces, ensuring...

  8. [8]

    The derivation process is similar to that in Appendix C 1, and we present the results directly below

    Cylindrically symmetric diffusion Within the framework of the operator splitting method, we discretize the 2D radially symmetric diffusion in the r direction and the 1D diffusion in the z direction separately. The derivation process is similar to that in Appendix C 1, and we present the results directly below. The semi-discrete scheme for radially symmetr...

  9. [9]

    Fang, Front

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

  10. [10]

    Liu, Int

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

  11. [11]

    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]

  12. [12]

    Amato and S

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

  13. [13]

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

  14. [14]

    Extended gamma-ray sources around pulsars constrain the origin of the positron flux at Earth

    A. Abeysekara et al. (HAWC), Science 358, 911 (2017), arXiv:1711.06223 [astro-ph.HE]

  15. [15]

    Using HAWC to Discover Invisible Pulsars

    T. Linden, K. Auchettl, J. Bramante, I. Cholis, K. Fang, D. Hooper, T. Karwal, and S. W. Li, Phys. Rev. D 96, 103016 (2017), arXiv:1703.09704 [astro-ph.HE]

  16. [16]

    Fang and X.-J

    K. Fang and X.-J. Bi, Phys. Rev. D 105, 103007 (2022), arXiv:2203.01546 [astro-ph.HE]

  17. [17]

    Albert et al

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

  18. [18]

    Aharonian et al

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

  19. [19]

    Fang, S.-Q

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

  20. [20]

    Cao et al

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

  21. [21]

    Albert et al

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

  22. [22]

    Albert et al

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

  23. [23]

    Kumar, arXiv e-prints , arXiv:2508.01934 (2025), arXiv:2508.01934 [astro-ph.HE]

    S. Kumar, arXiv e-prints , arXiv:2508.01934 (2025), arXiv:2508.01934 [astro-ph.HE]

  24. [24]

    Di Mauro, S

    M. Di Mauro, S. Manconi, and F. Donato, Phys. Rev. D 101, 103035 (2020), arXiv:1908.03216 [astro-ph.HE]

  25. [25]

    Albert et al., Phys

    A. Albert et al., Phys. Rev. Lett. 134, 171005 (2025), arXiv:2505.00175 [astro-ph.HE]

  26. [26]

    John and T

    I. John and T. Linden, arXiv e-prints , arXiv:2503.17442 (2025), arXiv:2503.17442 [astro- ph.HE]. 34

  27. [27]

    Fang, X.-J

    K. Fang, X.-J. Bi, and P.-F. Yin, Astrophys. J. 884, 124 (2019), arXiv:1906.08542 [astro- ph.HE]

  28. [28]

    Martin, A

    P. Martin, A. Marcowith, and L. Tibaldo, Astron. Astrophys. 665, A132 (2022), arXiv:2206.11803 [astro-ph.HE]

  29. [29]

    Bourguinat, C

    L.-M. Bourguinat, C. Evoli, P. Martin, and S. Recchia, arXiv e-prints , arXiv:2507.01495 (2025), arXiv:2507.01495 [astro-ph.HE]

  30. [30]

    Self-Generated Cosmic-Ray Confinement in TeV Halos: Implications for TeV Gamma-Ray Emission and the Positron Excess

    C. Evoli, T. Linden, and G. Morlino, Phys. Rev. D 98, 063017 (2018), arXiv:1807.09263 [astro-ph.HE]

  31. [31]

    Fang, X.-J

    K. Fang, X.-J. Bi, and P.-F. Yin, Mon. Not. Roy. Astron. Soc. 488, 4074 (2019), arXiv:1903.06421 [astro-ph.HE]

  32. [32]

    R.-Y. Liu, H. Yan, and H. Zhang, Phys. Rev. Lett. 123, 221103 (2019), arXiv:1904.11536 [astro-ph.HE]

  33. [33]

    Recchia, M

    S. Recchia, M. Di Mauro, F. A. Aharonian, L. Orusa, F. Donato, S. Gabici, and S. Manconi, Phys. Rev. D 104, 123017 (2021), arXiv:2106.02275 [astro-ph.HE]

  34. [34]

    Mukhopadhyay and T

    P. Mukhopadhyay and T. Linden, Phys. Rev. D 105, 123008 (2022), arXiv:2111.01143 [astro- ph.HE]

  35. [35]

    De La Torre Luque, O

    P. De La Torre Luque, O. Fornieri, and T. Linden, Phys. Rev. D 106, 123033 (2022), arXiv:2205.08544 [astro-ph.HE]

  36. [36]

    HAWC Observations Strongly Favor Pulsar Interpretations of the Cosmic-Ray Positron Excess

    D. Hooper, I. Cholis, T. Linden, and K. Fang, Phys. Rev. D 96, 103013 (2017), arXiv:1702.08436 [astro-ph.HE]

  37. [37]

    Two-zone diffusion of electrons and positrons from Geminga explains the positron anomaly

    K. Fang, X.-J. Bi, P.-F. Yin, and Q. Yuan, Astrophys. J. 863, 30 (2018), arXiv:1803.02640 [astro-ph.HE]

  38. [38]

    Pulsar TeV Halos Explain the TeV Excess Observed by Milagro

    T. Linden and B. J. Buckman, Phys. Rev. Lett. 120, 121101 (2018), arXiv:1707.01905 [astro- ph.HE]

  39. [39]

    Yan, R.-Y

    K. Yan, R.-Y. Liu, R. Zhang, C.-M. Li, Q. Yuan, and X.-Y. Wang, Nature Astron. 8, 628 (2024), arXiv:2307.12363 [astro-ph.HE]

  40. [40]

    Guo and Q

    Y.-Y. Guo and Q. Yuan, JHEAp 43, 227 (2024), arXiv:2407.10705 [astro-ph.HE]

  41. [41]

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

  42. [42]

    Giacinti, A

    G. Giacinti, A. Mitchell, R. L´ opez-Coto, V. Joshi, R. Parsons, and J. Hinton, Astron. Astro- phys. 636, A113 (2020), arXiv:1907.12121 [astro-ph.HE]. 35

  43. [43]

    Propagation of cosmic rays in the AMS-02 era

    Q. Yuan, S.-J. Lin, K. Fang, and X.-J. Bi, Phys. Rev. D 95, 083007 (2017), arXiv:1701.06149 [astro-ph.HE]

  44. [44]

    B. B. Wang, G. P. Zank, L. L. Zhao, and L. Adhikari, Astrophys. J. 932, 65 (2022)

  45. [45]

    Yao et al., Astrophys

    J. Yao et al., Astrophys. J. 939, 75 (2022), arXiv:2209.14059 [astro-ph.HE]

  46. [46]

    L.-Z. Bao, K. Fang, X.-J. Bi, and S.-H. Wang, Astrophys. J.936, 183 (2022), arXiv:2107.07395 [astro-ph.HE]

  47. [47]

    Barreto-Mota, E

    L. Barreto-Mota, E. M. de Gouveia Dal Pino, S. Xu, and A. Lazarian, arXiv e-prints , arXiv:2405.12146 (2024), arXiv:2405.12146 [astro-ph.HE]

  48. [48]

    Y. Xiao, J. Zhang, and S. Xu, arXiv e-prints , arXiv:2506.15031 (2025), arXiv:2506.15031 [astro-ph.HE]

  49. [49]

    Fang, X.-J

    K. Fang, X.-J. Bi, S.-J. Lin, and Q. Yuan, Chin. Phys. Lett. 38, 039801 (2021), arXiv:2007.15601 [astro-ph.HE]

  50. [50]

    J. M. Migliazzo, B. M. Gaensler, D. C. Backer, B. W. Stappers, E. van der Swaluw, and R. G. Strom, Astrophys. J. Lett. 567, L141 (2002), arXiv:astro-ph/0202063

  51. [51]

    S. B. Popov and R. Turolla, Astrophys. Space Sci. 341, 457 (2012), arXiv:1204.0632 [astro- ph.HE]

  52. [52]

    Suzuki, A

    H. Suzuki, A. Bamba, and S. Shibata, Astrophys. J. 914, 103 (2021), arXiv:2104.10052 [astro- ph.HE]

  53. [53]

    R. N. Manchester, G. B. Hobbs, A. Teoh, and M. Hobbs, Astron. J. 129, 1993 (2005), arXiv:astro-ph/0412641

  54. [54]

    Radiative losses and cut-offs of energetic particles at relativistic shocks

    P. Dempsey and P. Duffy, Mon. Not. Roy. Astron. Soc. 378, 625 (2007), arXiv:0704.0168 [astro-ph]

  55. [55]

    T. A. Porter, G. Johannesson, and I. V. Moskalenko, Astrophys. J. Supp. 262, 30 (2022), arXiv:2112.12745 [astro-ph.HE]

  56. [56]

    V. L. Ginzburg and S. I. Syrovatskii, The Origin of Cosmic Rays (Macmillan, New York, 1964)

  57. [57]

    Superluminal problem in diffusion of relativistic particles and its phenomenological solution

    R. Aloisio, V. Berezinsky, and A. Gazizov, Astrophys. J. 693, 1275 (2009), arXiv:0805.1867 [astro-ph]

  58. [58]

    Lv, X.-J

    X.-J. Lv, X.-J. Bi, K. Fang, P.-F. Yin, and M.-J. Zhao, Astrophys. J. 982, 85 (2025), arXiv:2410.19396 [astro-ph.HE]

  59. [59]

    Kawanaka and R

    N. Kawanaka and R. Takahashi, Phys. Rev. D 111, 123016 (2025), arXiv:2410.14829 [astro- 36 ph.HE]

  60. [60]

    Measuring the Local Diffusion Coefficient with H.E.S.S. Observations of Very High-Energy Electrons

    D. Hooper and T. Linden, Phys. Rev. D 98, 083009 (2018), arXiv:1711.07482 [astro-ph.HE]

  61. [61]

    Constraining the slow-diffusion zone size and electron injection spectral index for the Geminga pulsar halo

    K. Fang, Phys. Rev. D 109, 043041 (2024), arXiv:2310.16594 [astro-ph.HE]

  62. [62]

    S. M. Osipov, A. M. Bykov, A. E. Petrov, and V. I. Romansky, J. Phys. Conf. Ser. 1697, 012009 (2020)

  63. [63]

    Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients

    C. Evoli, D. Gaggero, A. Vittino, G. Di Bernardo, M. Di Mauro, A. Ligorini, P. Ullio, and D. Grasso, JCAP 2, 015, arXiv:1607.07886 [astro-ph.HE]

  64. [64]

    Fang, H.-B

    K. Fang, H.-B. Hu, X.-J. Bi, and E.-S. Chen, Phys. Rev. D 108, 023017 (2023), arXiv:2304.04570 [astro-ph.HE]

  65. [65]

    J. Xia, X. Bi, K. Fang, and S. Liu, Astrophys. J. 978, 162 (2025), arXiv:2409.01653 [astro- ph.HE]

  66. [66]

    S.-H. Wang, K. Fang, X.-J. Bi, and P.-F. Yin, Phys. Rev. D 103, 063035 (2021), arXiv:2101.01438 [astro-ph.HE]

  67. [67]

    Fang, S.-Q

    K. Fang, S.-Q. Xi, and X.-J. Bi, Phys. Rev. D 104, 103024 (2021), arXiv:2107.02140 [astro- ph.HE]

  68. [68]

    Y. Bao, G. Giacinti, R.-Y. Liu, H.-M. Zhang, and Y. Chen, arXiv e-prints , arXiv:2407.02478 (2024), arXiv:2407.02478 [astro-ph.HE]

  69. [69]

    K. Yan, S. Wu, and R.-Y. Liu, Astrophys. J. 987, 19 (2025), arXiv:2507.08526 [astro-ph.HE]

  70. [70]

    Di Mauro, S

    M. Di Mauro, S. Manconi, and F. Donato, Phys. Rev. D 100, 123015 (2019), [Erratum: Phys.Rev.D 104, 089903 (2021)], arXiv:1903.05647 [astro-ph.HE]

  71. [71]

    Zhang, R.-Y

    Y. Zhang, R.-Y. Liu, S. Z. Chen, and X.-Y. Wang, Astrophys. J. 922, 130 (2021)

  72. [72]

    Cosmic-Ray Propagation in Light of Recent Observation of Geminga

    G. Johannesson, T. A. Porter, and I. V. Moskalenko, Astrophys. J. 879, 91 (2019), arXiv:1903.05509 [astro-ph.HE]

  73. [73]

    S. Niu, Q. Yuan, S.-N. Zhang, L. Lei, L. Ji, and Y.-Z. Fan, arXiv e-prints , arXiv:2501.17046 (2025), arXiv:2501.17046 [astro-ph.HE]

  74. [74]

    Khokhriakova, W

    A. Khokhriakova, W. Becker, P. Predehl, J. S. Sanders, M. Freyberg, and A. Schwope, Res. Notes AAS 9, 154 (2025), arXiv:2507.02750 [astro-ph.HE]

  75. [75]

    R.-Y. Liu, C. Ge, X.-N. Sun, and X.-Y. Wang, Astrophys. J.875, 149 (2019), arXiv:1904.11438 [astro-ph.HE]

  76. [76]

    Khokhriakova, W

    A. Khokhriakova, W. Becker, G. Ponti, M. Sasaki, B. Li, and R. Y. Liu, Astron. Astrophys. 683, A180 (2024), arXiv:2310.10454 [astro-ph.HE]. 37

  77. [77]

    Manconi, J

    S. Manconi, J. Woo, R.-Y. Shang, R. Krivonos, C. Tang, M. Di Mauro, F. Donato, K. Mori, and C. J. Hailey, Astron. Astrophys. 689, A326 (2024), arXiv:2403.10902 [astro-ph.HE]

  78. [78]

    John and T

    I. John and T. Linden, Phys. Rev. D 107, 103021 (2023), arXiv:2206.04699 [astro-ph.HE]

  79. [79]

    J. Xia, X. Lv, K. Fang, and S. Liu, Phys. Rev. D 111, 123048 (2025), arXiv:2503.15052 [astro-ph.HE]

  80. [80]

    A. A. Lagutin, Y. A. Nikulin, and V. V. Uchaikin, Nuclear Physics B Proceedings Supplements 97, 267 (2001)

Showing first 80 references.