pith. sign in

arxiv: 1907.02037 · v1 · pith:5QREDWV7new · submitted 2019-07-03 · 🌌 astro-ph.HE

On the Gamma-Ray Nebula of Vela Pulsar. I. Very Slow Diffusion of Energetic Electrons within the TeV Nebula

Yiwei Bao , Siming Liu , Yang Chen This is my paper

Pith reviewed 2026-05-25 09:31 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords Vela pulsarpulsar wind nebulaVela XTeV gamma raysdiffusion coefficientenergetic electronsslow diffusionimpulsive injection
0
0 comments X

The pith

TeV electrons in the Vela X nebula diffuse more than a thousand times slower than in the interstellar medium.

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

The paper treats the electrons that produce the observed TeV gamma rays as having been released in one brief episode when the supernova remnant reverse shock compressed the pulsar wind nebula into its current cocoon shape. An exact solution of the equation governing particle spreading and energy loss then reproduces both the overall gamma-ray energy spectrum and the observed fall-off of brightness away from the center. This match fixes the diffusion coefficient at 1 times 10 to the 26 square centimeters per second for 10 TeV particles. The value lies more than three orders of magnitude below the rate measured in the surrounding galaxy and matches an earlier limit obtained for the Geminga nebula. The authors conclude that such slow movement of high-energy particles may be typical inside pulsar wind nebulae.

Core claim

By approximating the TeV-gamma-ray-emitting electrons in the Vela X PWN as injected impulsively at the moment the cocoon formed from the interaction between the SNR reverse shock and the PWN, and solving the diffusion-loss equation analytically, the broadband spectral energy distribution and surface brightness profile are reproduced simultaneously, determining the diffusion coefficient of TeV electrons and positrons to be 1 × 10^{26} cm² s^{-1} for 10 TeV particles. This coefficient is more than three orders of magnitude lower than that in the interstellar medium.

What carries the argument

Analytical solution to the diffusion-loss equation applied to electrons and positrons injected impulsively into the cocoon formed by SNR reverse shock interaction with the PWN.

If this is right

  • The size and spectrum of the TeV nebula directly constrain the diffusion rate of its parent particles.
  • The same slow diffusion rate accounts for the compact appearance of the gamma-ray emitting region.
  • Slow diffusion of high-energy particles is likely to occur in other pulsar wind nebulae as well.
  • The derived coefficient is consistent with the independent limit obtained for the Geminga TeV nebula.

Where Pith is reading between the lines

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

  • If slow diffusion is widespread in PWNe, models of how particles eventually escape into the interstellar medium and contribute to the cosmic-ray population would need to include a long residence time inside the nebula.
  • The result raises the possibility that ordered magnetic structures or enhanced turbulence inside PWNe are responsible for the suppression of particle motion.
  • Repeated application of the same impulsive-injection model to other well-observed PWNe could reveal whether the diffusion coefficient depends systematically on particle energy or on nebula age.

Load-bearing premise

The TeV-gamma-ray-emitting electrons are approximated as having been injected all at once when the cocoon formed due to the reverse shock interaction.

What would settle it

An independent estimate of the diffusion coefficient for 10 TeV particles inside Vela X, obtained for example from the spatial distribution of synchrotron X-rays, that differs by more than a factor of a few from 10^{26} cm² s^{-1}.

Figures

Figures reproduced from arXiv: 1907.02037 by Siming Liu, Yang Chen, Yiwei Bao.

Figure 1
Figure 1. Figure 1: Spatially integrated electron distribution function of annuli for α = 1.7. There are 10 annuli with a width of 12′ [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Fit to the γ-ray spectra of the inner and outer regions (left), and the normalized surface brightness profile (right). The H.E.S.S. data are taken from Abramowski et al. (2012), Fermi data are from Tibaldo et al. (2018). −           ◦        [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Spatial distribution of γ-ray indices. The data are taken from Abramowski et al. (2012) [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Dependence of the model fit on the diffusion coefficient for α = 1.7. The solid lines in the right panel represent the spectra of the inner region, and the dashed lines represent the spectra of the outer region. 4. DISCUSSION AND CONCLUSIONS The discovery of a high-energy spectral cutoff in combination with a recent detection of a very hard GeV spectrum indicate that the γ-ray emission of the TeV nebula of… view at source ↗
read the original abstract

High-energy particle transport in pulsar wind nebulae (PWNe) plays an essential role in explaining the characteristics revealed in multiwavelength observations. In this paper, the TeV-gamma-ray-emitting electrons in the Vela X PWN are approximated to be injected impulsively when the cocoon is formed due to the interaction between the SNR reverse shock and the PWN. By solving the diffusion-loss equation analytically, we reproduce the broadband spectral energy distribution and surface brightness profile simultaneously. The diffusion coefficient of TeV electrons and positrons, which is well constrained by the spectral and spatial properties of the TeV nebula, is thus determined to be $1 \times 10^{26}$\,cm$^{2}$\,s$^{-1}$ for 10\,TeV electrons and positrons. This coefficient is more than three orders of magnitude lower than that in the interstellar medium, in agreement with a constraint recently obtained from HAWC observations of a TeV nebula associated with the Geminga pulsar. These results suggest that slow diffusion of high-energy particles might be common in PWNe.

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 / 0 minor

Summary. The paper models the TeV gamma-ray nebula around the Vela pulsar by approximating the injection of energetic electrons as impulsive at the epoch when the SNR reverse shock forms the cocoon. An analytical solution to the diffusion-loss equation is used to simultaneously reproduce the broadband spectral energy distribution and the surface brightness profile, yielding a diffusion coefficient D = 1 × 10^{26} cm² s^{-1} at 10 TeV that is more than three orders of magnitude below typical interstellar-medium values. The result is presented as evidence that slow diffusion may be common in PWNe, consistent with recent HAWC constraints on the Geminga nebula.

Significance. If the central result holds, the work supplies a concrete, observationally anchored value for suppressed diffusion inside a PWN, with direct implications for cosmic-ray transport and multiwavelength modeling of such sources. The use of an analytical solution to the transport equation is a clear methodological strength, enabling a joint fit to both spectrum and morphology without requiring full numerical simulations.

major comments (1)
  1. [Model setup and abstract] The impulsive-injection assumption at cocoon formation (stated in the abstract) is load-bearing for the quoted value of D. Because the spatial and spectral evolution of the nebula depend on the injection time profile, a continuous or extended injection history tied to ongoing pulsar activity would generally require a different diffusion coefficient to match the same SED and surface-brightness data; the manuscript does not quantify this sensitivity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback. We address the single major comment point by point below.

read point-by-point responses

  1. Referee: The impulsive-injection assumption at cocoon formation (stated in the abstract) is load-bearing for the quoted value of D. Because the spatial and spectral evolution of the nebula depend on the injection time profile, a continuous or extended injection history tied to ongoing pulsar activity would generally require a different diffusion coefficient to match the same SED and surface-brightness data; the manuscript does not quantify this sensitivity.

    Authors: We agree that the impulsive-injection assumption is central to the quoted diffusion coefficient and that the manuscript does not quantify sensitivity to alternative injection histories. The approximation is adopted because cocoon formation marks the epoch of significant confinement by the reverse shock, after which the nebula's evolution is modeled via diffusion and losses; the short cooling timescale of TeV electrons (~few kyr) limits the contribution from much earlier injection. Nevertheless, a continuous injection profile tied to ongoing pulsar activity would indeed alter the required D. In the revised manuscript we will add a short discussion (in Section 4) comparing the impulsive case to a simple continuous-injection model and showing that the diffusion coefficient remains suppressed by at least two orders of magnitude relative to typical ISM values, preserving the main conclusion. revision: yes

Circularity Check

0 steps flagged

No circularity: D obtained via standard fitting of analytical diffusion-loss solution to multiwavelength data

full rationale

The paper assumes impulsive injection at cocoon formation, solves the diffusion-loss equation analytically, and adjusts the free parameter D to simultaneously match the observed TeV SED and surface-brightness profile. This is ordinary parameter inference from data under stated model assumptions, not a reduction of the claimed result to its inputs by construction. No self-citations appear in the provided text, no uniqueness theorems are invoked, and the derived D is compared to an independent HAWC constraint on Geminga. The result is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central numerical claim rests on one fitted parameter (the diffusion coefficient) and one structural modeling choice (impulsive injection at cocoon formation time). No new physical entities are introduced.

free parameters (1)
  • diffusion coefficient D at 10 TeV = 1e26 cm2/s
    Value chosen to simultaneously reproduce the broadband spectral energy distribution and the surface brightness profile of the TeV nebula.
axioms (1)
  • domain assumption TeV electrons and positrons are injected impulsively at the moment the SNR reverse shock forms the cocoon around the PWN
    This timing assumption allows the diffusion-loss equation to be solved analytically and directly determines the required diffusion coefficient.

pith-pipeline@v0.9.0 · 5729 in / 1422 out tokens · 38714 ms · 2026-05-25T09:31:20.478764+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

34 extracted references · 34 canonical work pages · 2 internal anchors

  1. [1]

    A., Ackermann, M., Ajello, M., et al

    Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJ, 713, 146

    work page 2010

  2. [2]

    2017, PhRvD, 95, 082007

    Abdollahi, S., Ackermann, M., Ajello, M., et al. 2017, PhRvD, 95, 082007

    work page 2017

  3. [3]

    U., Albert, A., Alfaro, R., et al

    Abeysekara, A. U., Albert, A., Alfaro, R., et al. 2017, Science, 358, 911

    work page 2017

  4. [4]

    2012, A&A, 548, A38

    Abramowski, A., Acero, F., Aharonian, F., et al. 2012, A&A, 548, A38

    work page 2012

  5. [5]

    2014, Physical Review Letters, 113, 121101

    Accardo, L., Aguilar, M., Aisa, D., et al. 2014, Physical Review Letters, 113, 121101

    work page 2014

  6. [6]

    2012, Physical Review Letters, 108, 011103

    Ackermann, M., Ajello, M., Allafort, A., et al. 2012, Physical Review Letters, 108, 011103

    work page 2012

  7. [7]

    C., Bazilevskaya, G

    Adriani, O., Barbarino, G. C., Bazilevskaya, G. A., et al. 2013, Physical Review Letters, 111, 081102

    work page 2013

  8. [8]

    2019, Physical Review Letters, 122, 041102

    Aguilar, M., Ali Cavasonza, L., Ambrosi, G., et al. 2019, Physical Review Letters, 122, 041102

    work page 2019

  9. [9]

    G., Bazer-Bachi, A

    Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 448, L43

    work page 2006

  10. [10]

    1995, Nature, 373, 587

    Aschenbach, B., Egger, R., & Tr¨ umper, J. 1995, Nature, 373, 587

    work page 1995

  11. [11]

    M., Aharonian, F

    Atoyan, A. M., Aharonian, F. A., & V¨ olk, H. J. 1995, PhRvD, 52, 3265

    work page 1995

  12. [12]

    2019, ApJ, accepted (Paper II)

    Bao, Y., & Chen, Y. 2019, ApJ, accepted (Paper II)

    work page 2019

  13. [13]

    Blasi, P., Amato, E., & Serpico, P. D. 2012, Physical Review Letters, 109, 061101 DAMPE Collaboration, Ambrosi, G., An, Q., et al. 2017, Nature, 552, 63 de Jager, O. C., Slane, P. O., & LaMassa, S. 2008, ApJL, 689, L125

    work page 2012

  14. [14]

    2018, ApJ, 863, 30

    Fang, K., Bi, X.-J., Yin, P.-F., & Yuan, Q. 2018, ApJ, 863, 30

    work page 2018

  15. [15]

    A., Bietenholz, M

    Frail, D. A., Bietenholz, M. F., & Markwardt, C. B. 1997, ApJ, 475, 224

    work page 1997

  16. [16]

    J., & Takahara, F

    Fujita, Y., Ohira, Y., Tanaka, S. J., & Takahara, F. 2009, ApJL, 707, L179

    work page 2009

  17. [17]

    2010, SF2A-2010: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, 313

    Rowell, G. 2010, SF2A-2010: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, 313

    work page 2010

  18. [18]

    M., & Slane, P

    Gaensler, B. M., & Slane, P. O. 2006, ARA&A, 44, 17

    work page 2006

  19. [19]

    A., Funk, S., Parsons, R

    Hinton, J. A., Funk, S., Parsons, R. D., & Ohm, S. 2011, ApJL, 743, L7 10 Bao et al

    work page 2011

  20. [20]

    Hooper, D., Blasi, P., & Serpico, P. D. 2009, JCAP, 1, 025

    work page 2009

  21. [21]

    2017, PhRvD, 96, 103013

    Hooper, D., Cholis, I., Linden, T., & Fang, K. 2017, PhRvD, 96, 103013

    work page 2017

  22. [22]

    2018, ApJ, 866, 143

    Huang, Z.-Q., Fang, K., Liu, R.-Y., & Wang, X.-Y. 2018, ApJ, 866, 143

    work page 2018

  23. [23]

    2010, MNRAS, 409, L35

    Li, H., & Chen, Y. 2010, MNRAS, 409, L35

    work page 2010

  24. [24]

    A., & Moskalenko, I

    Aharonian, F. A., & Moskalenko, I. V. 2013, ApJ, 768, 73

    work page 2013

  25. [25]

    B., & ¨Ogelman, H

    Markwardt, C. B., & ¨Ogelman, H. 1995, Nature, 375, 40

    work page 1995

  26. [26]

    2008, ApJ, 682, 400

    Orlando, E., & Bouchet, L. 2008, ApJ, 682, 400

    work page 2008

  27. [27]

    A local fading accelerator and the origin of TeV cosmic ray electrons

    Recchia, S., Gabici, S., Aharonian, F. A., & Vink, J. 2018, arXiv:1811.07551

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  28. [28]

    P., Pavlov, G

    Reynolds, S. P., Pavlov, G. G., Kargaltsev, O., et al. 2017, SSRv, 207, 175

    work page 2017

  29. [29]

    2019, MNRAS,

    Shi, Z.-D., & Liu, S. 2019, MNRAS,

    work page 2019

  30. [30]

    2018, ApJ, 865, 86

    Slane, P., Lovchinsky, I., Kolb, C., et al. 2018, ApJ, 865, 86

    work page 2018

  31. [31]

    Staszak, D., & VERITAS Collaboration 2015, 34th International Cosmic Ray Conference (ICRC2015), 34, 411

    work page 2015

  32. [32]

    2018, A&A, 617, A78

    Tibaldo, L., Zanin, R., Faggioli, G., et al. 2018, A&A, 617, A78

    work page 2018

  33. [33]

    GeV observations of the extended pulsar wind nebulae constrain the pulsar interpretations of the cosmic-ray positron excess

    Xi, S.-Q., Liu, R.-Y., Huang, Z.-Q., et al. 2018, arXiv:1810.10928

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  34. [34]

    2012, ApJ, 745, 140 10-1 100 101 102 Energy [TeV] 10-14 10-13 10-12 10-11 10-10 Flux [TeV−1 cm−2 s−1] SED of Inner and Outer Region Inner Region Outer Region

    Yan, H., Lazarian, A., & Schlickeiser, R. 2012, ApJ, 745, 140 10-1 100 101 102 Energy [TeV] 10-14 10-13 10-12 10-11 10-10 Flux [TeV−1 cm−2 s−1] SED of Inner and Outer Region Inner Region Outer Region

    work page 2012