pith. sign in

arxiv: 2507.10276 · v3 · pith:AQLKY27Cnew · submitted 2025-07-14 · 🌌 astro-ph.HE

Method for testing diffusive shock acceleration and diffusion propagation of 1-100 TeV cosmic electrons with multiwavelength observations of the Geminga halo and pulsar wind nebula

Weikang Gao , Li-Zhuo Bao , Kun Fang , En-sheng Chen , Siming Liu , HongBo Hu This is my paper

Pith reviewed 2026-05-22 00:50 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords diffusive shock accelerationdiffusion propagationGeminga halopulsar wind nebulacosmic ray electronsTeV gamma raysHAWCLHAASO
0
0 comments X

The pith

Geminga observations confirm diffusive shock acceleration and diffusion for cosmic electrons between 1 and 100 TeV.

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

This paper presents a method for testing the diffusive shock acceleration and diffusion propagation theories using observations of the Geminga pulsar wind nebula and its associated halo. By analyzing multiwavelength spectra and angular profiles from instruments like HAWC and Fermi-LAT, the authors compare data to model predictions in the 1-100 TeV range. The results indicate that these standard cosmic ray models are consistent with the experimental observations. This testing elevates validation of the theories from MeV to sub-PeV energies. However, current data limitations prevent high-precision checks on the energy dependence of the diffusion coefficient.

Core claim

A method is developed to test the diffusive shock acceleration and diffusion propagation model between one and several hundred TeV energies through the latest spectral and morphological data of the Geminga region from HAWC and Fermi-LAT. Our results show that the theories of diffusive shock acceleration and diffusion propagation are consistent with experimental observations. Future HAWC and LHAASO-KM2A observations will yield higher-precision results, and the confirmation of a rapidly increasing diffusion coefficient above 100 TeV would serve as important evidence supporting the diffusive shock acceleration and diffusion propagation theory.

What carries the argument

multiwavelength spectral and morphological fitting of Geminga halo and pulsar wind nebula data to predictions from diffusive shock acceleration combined with energy-dependent diffusion

If this is right

  • The standard cosmic ray model remains valid for electron acceleration and propagation at TeV energies.
  • A confirmed rapid rise in the diffusion coefficient above 100 TeV supplies direct support for diffusive shock acceleration at these energies.
  • The same comparison technique can be applied to other pulsar wind nebulae and halos for additional tests.
  • Higher-precision data will enable quantitative checks on the exact energy dependence of diffusion.

Where Pith is reading between the lines

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

  • If the diffusion increase holds, it may help explain why cosmic-ray spectra appear sharper at higher energies than simple propagation models predict.
  • The method offers a template for testing whether diffusion properties are universal across different regions of the interstellar medium.
  • Similar multiwavelength tests on other sources could distinguish between competing cosmic-ray origin scenarios at TeV scales.

Load-bearing premise

The wide energy bins in published morphological data still allow a consistency test between observed spectra and the energy dependence predicted by diffusive shock acceleration theory.

What would settle it

A future narrow-bin measurement of the angular profile or spectrum above 100 TeV showing no rapid increase in the diffusion coefficient would falsify the claimed consistency.

Figures

Figures reproduced from arXiv: 2507.10276 by En-sheng Chen, Hongbo Hu, Kun Fang, Li-Zhuo Bao, Siming Liu, Weikang Gao.

Figure 2
Figure 2. Figure 2: FIG. 2. Spectra fit of Geminga halo. The black data points [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The diffusion coefficients corresponding to the fitted parameters under the condition that [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Brightness profile of Geminga halo measued by HAWC in two energy bins. The top two figures are profiles under [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
read the original abstract

Diffusive shock acceleration and diffusion propagation are essential components of the standard cosmic ray model. These theories are based on extensive observations of high-energy solar processes, providing substantial direct evidence in the MeV energy range. Although the model is widely and successfully used to explain high-energy cosmic phenomena, direct validation has been elusive. The multi-wavelength spectra and angular profile measurements of the Geminga pulsar wind nebula and its pulsar halo, particularly the precise spectral observations by HAWC and LHAASO-KM2A in recent years, offer a rare opportunity to test these theories with cosmic rays energies between 1 TeV and several hundred TeV. These observations are expected to elevate the direct testing of theoretical models from multi-MeV to sub-PeV energies. In this work, a method is developed to test the diffusive shock acceleration and diffusion propagation model between one and several hundred TeV energies through the latest spectral and morphological data of the Geminga region from HAWC and Fermi-LAT. Our results show that the theories of diffusive shock acceleration and diffusion propagation are consistent with experimental observations. However, the published morphological data adopted rather wide energy bins and currently do not allow a high precision test of the inferred energy dependent diffusion coefficient by observed energy spectra with DSA theory. It is anticipated that future HAWC and LHAASO-KM2A observations will yield higher-precision results, and the confirmation of a rapidly increasing diffusion coefficient above 100 TeV would serve as important evidence supporting the diffusive shock acceleration and diffusion propagation theory. Similar tests would be both important and valuable for other models.

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

Summary. The paper develops a method to test diffusive shock acceleration (DSA) and diffusion propagation using multi-wavelength spectra and angular profiles of the Geminga pulsar wind nebula and halo, drawing on HAWC and Fermi-LAT data in the 1-100 TeV range. It reports consistency between the theories and observations while explicitly noting that wide energy bins in the published morphological data preclude a high-precision test of the inferred energy-dependent diffusion coefficient against DSA expectations. Future higher-resolution observations are expected to confirm a rapidly rising diffusion coefficient above 100 TeV.

Significance. If the central consistency result holds under improved data, the work supplies a concrete route to validate DSA and diffusion models at TeV-sub-PeV energies using a well-studied source, extending direct tests beyond the MeV regime. The approach is grounded in existing multi-messenger observations and highlights a falsifiable signature (rapid D(E) increase above 100 TeV), giving it moderate significance pending narrower energy bins.

major comments (1)
  1. [Abstract] Abstract: the central claim of consistency between DSA/diffusion propagation and the 1-100 TeV Geminga data rests on an inferred energy-dependent diffusion coefficient D(E); however, the text itself states that wide energy bins in HAWC morphological data prevent high-precision comparison of this D(E) with DSA predictions, so alternative functional forms remain compatible within the same bins and the support for the claim is correspondingly limited.
minor comments (1)
  1. Clarify in the methods or results section how the morphological fitting procedure extracts D(E) and whether the procedure remains stable when the energy-bin width is varied within the published uncertainties.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful review and constructive feedback on our manuscript. We address the single major comment below and have revised the abstract to better reflect the limitations of the current data while preserving the core contribution of the testing method.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim of consistency between DSA/diffusion propagation and the 1-100 TeV Geminga data rests on an inferred energy-dependent diffusion coefficient D(E); however, the text itself states that wide energy bins in HAWC morphological data prevent high-precision comparison of this D(E) with DSA predictions, so alternative functional forms remain compatible within the same bins and the support for the claim is correspondingly limited.

    Authors: We agree that the wide energy bins in the published HAWC morphological data limit the precision of any comparison between the inferred D(E) and DSA predictions, allowing alternative functional forms to remain compatible within uncertainties. This limitation is already stated explicitly in the original abstract and main text. To address the referee's concern, we have revised the abstract to replace the phrasing 'consistent with experimental observations' with 'broadly consistent with the DSA and diffusion propagation model within the uncertainties set by current energy binning,' while emphasizing that the primary advance is the development of a testable method and that narrower bins from future observations will enable higher-precision discrimination. This revision accurately reflects the evidential weight of the present data without overstating it. revision: yes

Circularity Check

0 steps flagged

No significant circularity; consistency claim rests on external data comparison

full rationale

The paper develops a method using HAWC and Fermi-LAT spectral and morphological observations of the Geminga region to test DSA and diffusion propagation models in the 1-100 TeV range. It reports that the theories are consistent with the data while explicitly acknowledging that wide energy bins in published morphological data prevent high-precision comparison of any inferred energy-dependent diffusion coefficient against DSA expectations. No equations or steps are quoted that reduce a claimed prediction to a fitted input by construction, nor is there load-bearing self-citation or ansatz smuggling; the test uses independent external observations rather than redefining the target result in terms of itself.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that Geminga multiwavelength data can test DSA and diffusion at TeV energies. No new entities are postulated. The energy-dependent diffusion coefficient is treated as inferable from data, which functions as a fitted element.

free parameters (1)
  • energy dependent diffusion coefficient
    Inferred from observed energy spectra to enable comparison with DSA theory; specific functional form or normalization not detailed in abstract.
axioms (1)
  • domain assumption Geminga pulsar wind nebula and halo observations are suitable for testing diffusive shock acceleration and diffusion propagation between 1 and several hundred TeV
    Invoked in the abstract as providing a rare opportunity to test the standard model at these energies.

pith-pipeline@v0.9.0 · 5851 in / 1419 out tokens · 71632 ms · 2026-05-22T00:50:56.520873+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

50 extracted references · 50 canonical work pages · 3 internal anchors

  1. [1]

    Amato, Particle acceleration in pulsars and pulsar wind nebulae, inFoundations of Cosmic Ray Astro- physics(IOS Press, 2024) pp

    E. Amato, Particle acceleration in pulsars and pulsar wind nebulae, inFoundations of Cosmic Ray Astro- physics(IOS Press, 2024) pp. 183–221

    work page 2024

  2. [2]

    L. O. Drury, An introduction to the theory of diffusive shock acceleration of energetic particles in tenuous plas- mas, Reports on Progress in Physics46, 973 (1983)

    work page 1983

  3. [3]

    Sironi and A

    L. Sironi and A. Spitkovsky, Acceleration of particles at the termination shock of a relativistic striped wind, The Astrophysical Journal741, 39 (2011)

    work page 2011

  4. [4]

    Hoshino, J

    M. Hoshino, J. Arons, Y. A. Gallant, and A. B. Langdon, Relativistic magnetosonic shock waves in synchrotron sources - Shock structure and nonthermal acceleration of positrons, Astrophys. J.390, 454 (1992)

    work page 1992

  5. [5]

    Klecker, R

    B. Klecker, R. Mewaldt, J. Bieber, A. Cummings, L. Drury, J. Giacalone, J. Jokipii, F. Jones, M. Krainev, M. Lee,et al., Anomalous cosmic rays, Space Science Re- 10 views83, 259 (1998)

    work page 1998

  6. [6]

    Stanev and Stanev,High energy cosmic rays, Vol

    T. Stanev and Stanev,High energy cosmic rays, Vol. 175 (Springer, 2010)

    work page 2010

  7. [7]

    Hillas, Can diffusive shock acceleration in supernova remnants account for high-energy galactic cosmic rays?, Journal of Physics G: Nuclear and Particle Physics31, R95 (2005)

    A. Hillas, Can diffusive shock acceleration in supernova remnants account for high-energy galactic cosmic rays?, Journal of Physics G: Nuclear and Particle Physics31, R95 (2005)

    work page 2005

  8. [8]

    S. P. Reynolds, Supernova remnants at high energy, Annu. Rev. Astron. Astrophys.46, 89 (2008)

    work page 2008

  9. [9]

    Florinski, S

    V. Florinski, S. Ferreira, and N. Pogorelov, Galactic cos- mic rays in the outer heliosphere: Theory and models, Space Science Reviews176, 147 (2013)

    work page 2013

  10. [10]

    J. J. Kroon, P. A. Becker, J. D. Finke, and C. D. Dermer, Electron acceleration in pulsar-wind termination shocks: an application to the crab nebula gamma-ray flares, The Astrophysical Journal833, 157 (2016)

    work page 2016

  11. [11]

    Florinski, G

    V. Florinski, G. Zank, and N. Pogorelov, Galactic cos- mic ray transport in the global heliosphere, Journal of Geophysical Research: Space Physics108(2003)

    work page 2003

  12. [12]

    K´ ota, Theory and modeling of galactic cosmic rays: Trends and prospects, Space Science Reviews176, 391 (2013)

    J. K´ ota, Theory and modeling of galactic cosmic rays: Trends and prospects, Space Science Reviews176, 391 (2013)

    work page 2013

  13. [13]

    A. W. Strong, I. V. Moskalenko, and V. S. Ptuskin, Cosmic-ray propagation and interactions in the galaxy, Annu. Rev. Nucl. Part. Sci.57, 285 (2007)

    work page 2007

  14. [14]

    Amato, The streaming instability: a review, Memorie della Societa Astronomica Italiana, v

    E. Amato, The streaming instability: a review, Memorie della Societa Astronomica Italiana, v. 82, p. 806 (2011) 82, 806 (2011)

    work page 2011

  15. [15]

    Fang, X.-J

    K. Fang, X.-J. Bi, and P.-F. Yin, Possible origin of the slow-diffusion region around geminga, Monthly Notices of the Royal Astronomical Society488, 4074 (2019)

    work page 2019

  16. [16]

    Albert, R

    A. Albert, R. Alfaro, C. Alvarez, J. Arteaga-Vel´ azquez, D. A. Rojas, H. A. Solares, R. Babu, E. Belmont-Moreno, A. Bernal, K. Caballero-Mora,et al., Precise measure- ments of tev halos around geminga and monogem pulsars with hawc, The Astrophysical Journal974, 246 (2024)

    work page 2024

  17. [17]

    Y. Guo, Y. Zhang, Q. Yuan, H. Hu, Z. Cao, F. Aha- ronian, Q. An, L. Bai, Y. Bai, Y. Bao,et al., Observa- tions of extended very-high-energy halos around geminga and monogem with the lhaaso-km2a, PoS ICRC2021964 (2021)

    work page 2021

  18. [18]

    R. N. Manchester, G. B. Hobbs, A. Teoh, and M. Hobbs, The australia telescope national facility pulsar catalogue, The Astronomical Journal129, 1993 (2005)

    work page 1993

  19. [19]

    A. A. Abdo, B. Allen, D. Berley, S. Casanova, C. Chen, D. Coyne, B. Dingus, R. Ellsworth, L. Fleysher, R. Fleysher,et al., Tev gamma-ray sources from a sur- vey of the galactic plane with milagro, The Astrophysical Journal664, L91 (2007)

    work page 2007

  20. [20]

    Abeysekara, A

    A. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. ´Alvarez, R. Arceo, J. Arteaga-Vel´ azquez, D. Avila Ro- jas, H. Ayala Solares, A. Barber,et al., Extended gamma- ray sources around pulsars constrain the origin of the positron flux at earth, Science358, 911 (2017)

    work page 2017

  21. [21]

    Di Mauro, S

    M. Di Mauro, S. Manconi, and F. Donato, Detection of aγ-ray halo around geminga with the fermi-lat data and implications for the positron flux, Physical Review D100, 123015 (2019)

    work page 2019

  22. [22]

    R.-Y. Liu, H. Yan, and H. Zhang, Understanding the multiwavelength observation of geminga’s tev halo: the role of anisotropic diffusion of particles, Physical review letters123, 221103 (2019)

    work page 2019

  23. [23]

    Evoli, T

    C. Evoli, T. Linden, and G. Morlino, Self-generated cosmic-ray confinement in tev halos: Implications for tev γ-ray emission and the positron excess, Physical Review D98, 063017 (2018)

    work page 2018

  24. [24]

    Caraveo, G

    P. Caraveo, G. Bignami, A. DeLuca, S. Mereghetti, A. Pellizzoni, R. Mignani, A. Tur, and W. Becker, Geminga’s tails: A pulsar bow shock probing the inter- stellar medium, Science301, 1345 (2003)

    work page 2003

  25. [25]

    G. G. Pavlov, S. Bhattacharyya, and V. E. Zavlin, New x-ray observations of the geminga pulsar wind nebula, The Astrophysical Journal715, 66 (2010)

    work page 2010

  26. [26]

    C. Hui, J. Lee, A. Kong, P. Tam, J. Takata, K. Cheng, and D. Ryu, Rapid x-ray variations of the geminga pulsar wind nebula, The Astrophysical Journal846, 116 (2017)

    work page 2017

  27. [27]

    Posselt, G

    B. Posselt, G. G. Pavlov, P. O. Slane, R. Romani, N. Buc- ciantini, A. M. Bykov, O. Kargaltsev, M. C. Weisskopf, and C.-Y. Ng, Geminga’s puzzling pulsar wind nebula, The Astrophysical Journal835, 66 (2017)

    work page 2017

  28. [28]

    Pellizzoni, F

    A. Pellizzoni, F. Govoni, P. Esposito, M. Murgia, and A. Possenti, Detection of continuum radio emission asso- ciated with geminga, Monthly Notices of the Royal As- tronomical Society: Letters416, L45 (2011)

    work page 2011

  29. [29]

    Giacinti, A

    G. Giacinti, A. Mitchell, R. L´ opez-Coto, V. Joshi, R. Par- sons, and J. Hinton, Halo fraction in tev-bright pul- sar wind nebulae, Astronomy & Astrophysics636, A113 (2020)

    work page 2020

  30. [30]

    Manconi, J

    S. Manconi, J. Woo, R.-Y. Shang, R. Krivonos, C. Tang, M. Di Mauro, F. Donato, K. Mori, and C. J. Hailey, Geminga’s pulsar halo: An x-ray view, Astronomy & As- trophysics689, A326 (2024)

    work page 2024

  31. [31]

    Caprioli, P

    D. Caprioli, P. Blasi, and E. Amato, On the escape of par- ticles from cosmic ray modified shocks, Monthly Notices of the Royal Astronomical Society396, 2065 (2009)

    work page 2065

  32. [32]

    Florinski and J

    V. Florinski and J. Jokipii, Cosmic-ray spectra at spher- ical termination shocks, The Astrophysical Journal591, 454 (2003)

    work page 2003

  33. [33]

    Fang, X.-J

    K. Fang, X.-J. Bi, S.-J. Lin, and Q. Yuan, Klein–nishina effect and the cosmic ray electron spectrum, Chinese Physics Letters38, 039801 (2021)

    work page 2021

  34. [34]

    G. R. Blumenthal and R. J. Gould, Bremsstrahlung, synchrotron radiation, and compton scattering of high- energy electrons traversing dilute gases, Reviews of mod- ern Physics42, 237 (1970)

    work page 1970

  35. [35]

    naima: a Python package for inference of relativistic particle energy distributions from observed nonthermal spectra

    V. Zabalza, naima: a python package for inference of relativistic particle energy distributions from observed nonthermal spectra, arXiv preprint arXiv:1509.03319 (2015)

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  36. [36]

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

    S.-Q. Xi, R.-Y. Liu, Z.-Q. Huang, K. Fang, and X.- Y. Wang, GeV observations of the extended pulsar wind nebulae constrain the pulsar interpretations of the cosmic-ray positron excess, Astrophys. J.878, 104 (2019), arXiv:1810.10928 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  37. [37]

    Yuan, S.-J

    Q. Yuan, S.-J. Lin, K. Fang, and X.-J. Bi, Propagation of cosmic rays in the ams-02 era, Physical Review D95, 083007 (2017)

    work page 2017

  38. [38]

    Foreman-Mackey, D

    D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Good- man, emcee: The mcmc hammer, Publications of the As- tronomical Society of the Pacific125, 306 (2013)

    work page 2013

  39. [39]

    J. Kirk, A. Guthmann, Y. Gallant, and A. Achterberg, Particle acceleration at ultrarelativistic shocks: an eigen- function method, The Astrophysical Journal542, 235 (2000)

    work page 2000

  40. [40]

    Z. Cao, D. della Volpe, S. Liu, X. Bi, Y. Chen, B. Piaz- zoli, L. Feng, H. Jia, Z. Li, X. Ma,et al., The large high altitude air shower observatory (lhaaso) science book 11 (2021 edition), arXiv preprint arXiv:1905.02773 (2019)

    work page arXiv 2021

  41. [41]

    Faherty, F

    J. Faherty, F. M. Walter, and J. Anderson, The trigono- metric parallax of the neutron star geminga, inIso- lated Neutron Stars: From the Surface to the Interior (Springer, 2007) pp. 225–230

    work page 2007

  42. [42]

    J. R. Knies, M. Sasaki, and P. P. Plucinsky, Suzaku obser- vations of the monogem ring and the origin of the gemini hαring, Monthly Notices of the Royal Astronomical So- ciety477, 4414 (2018)

    work page 2018

  43. [43]

    A. M. Bykov, E. Amato, A. E. Petrov, A. M. Kras- silchtchikov, and K. P. Levenfish, Pulsar wind nebu- lae with bow shocks: non-thermal radiation and cos- mic ray leptons, Space Sci. Rev.207, 235 (2017), arXiv:1705.00950 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  44. [44]

    Dempsey and P

    P. Dempsey and P. Duffy, Radiative losses and cut-offs of energetic particles at relativistic shocks, Monthly Notices of the Royal Astronomical Society378, 625 (2007)

    work page 2007

  45. [45]

    Heavens and K

    A. Heavens and K. Meisenheimer, Particle acceleration in extragalactic sources: the role of synchrotron losses in determining the spectrum, Monthly Notices of the Royal Astronomical Society225, 335 (1987)

    work page 1987

  46. [46]

    Sironi and B

    L. Sironi and B. Cerutti, Particle acceleration in pul- sar wind nebulae: Pic modelling, Modelling Pulsar Wind Nebulae , 247 (2017)

    work page 2017

  47. [47]

    P´ etri and Y

    J. P´ etri and Y. Lyubarsky, Magnetic reconnection at the termination shock of a striped pulsar wind, International Journal of Modern Physics D17, 1961 (2008)

    work page 1961

  48. [48]

    Xi, R.-Y

    S.-Q. Xi, R.-Y. Liu, Z.-Q. Huang, K. Fang, and X.-Y. Wang, Gev observations of the extended pulsar wind neb- ulae constrain the pulsar interpretations of the cosmic- ray positron excess, The Astrophysical Journal878, 104 (2019)

    work page 2019

  49. [49]

    P. A. Caraveo, G. F. Bignami, R. Mignani, and L. G. Taff, Parallax observations with the hubble space tele- scope yield the distance to geminga, The Astrophysical Journal461, L91 (1996)

    work page 1996

  50. [50]

    Fang, H.-B

    K. Fang, H.-B. Hu, X.-J. Bi, and E.-S. Chen, Effect of the magnetic field correlation length on the gamma-ray pul- sar halo morphology under anisotropic diffusion, Physical Review D108, 023017 (2023)

    work page 2023