pith. sign in

arxiv: 2603.07061 · v2 · pith:YRWTZ447new · submitted 2026-03-07 · 🌌 astro-ph.HE

Anisotropic Diffusion in Pulsar Halos: Interpreting the asymmetric morphology of Geminga and Monogem halos measured by HAWC

Si-Zhe Wu , Chao-Ming Li , Ruo-Yu Liu This is my paper

Pith reviewed 2026-05-25 07:06 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords pulsar halosanisotropic diffusioninterstellar magnetic fieldsGemingaMonogemAlfvénic Mach numbermagnetic coherence lengthHAWC
0
0 comments X

The pith

Modeling of Geminga and Monogem halos shows distinct mean magnetic field orientations but similar Alfvénic Mach numbers near 0.2.

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

The paper applies an anisotropic diffusion model to the asymmetric shapes of the Geminga and Monogem pulsar halos observed by HAWC. It concludes that the mean magnetic field directions differ in the two regions, meaning the pulsars sit in separate magnetic coherence zones, while the turbulence strength measured by the Alfvénic Mach number stays close to 0.2 in both. This leads to an estimated local magnetic field coherence length of about 100 parsecs. The work positions pulsar halo morphologies as a tool for probing interstellar magnetic turbulence. A sympathetic reader would see this as evidence that diffusion is not uniform but shaped by local field geometry.

Core claim

The asymmetric morphologies of the Geminga and Monogem halos arise from anisotropic diffusion of electrons and positrons in the interstellar medium, where the shape depends on the viewing angle of the mean magnetic field, the Alfvénic Mach number, and the pulsar distance. Fitting the HAWC data yields different magnetic field orientations for the two halos but comparable Mach numbers around 0.2, implying a coherence length of approximately 100 pc.

What carries the argument

Anisotropic diffusion model for pulsar halos, where diffusion is faster along the mean magnetic field and suppressed perpendicular to it, controlled by the Alfvénic Mach number.

If this is right

  • The two halos occupy different magnetic coherence regions.
  • Both regions have Alfvénic Mach numbers near 0.2.
  • The local magnetic coherence length is approximately 100 pc.
  • Pulsar halo shapes serve as a diagnostic for interstellar magnetic turbulence properties.

Where Pith is reading between the lines

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

  • High-resolution imaging of additional halos could trace magnetic field changes across the local interstellar medium.
  • The same modeling framework may apply to halos observed by LHAASO.
  • Results could refine cosmic-ray transport calculations that assume uniform diffusion.

Load-bearing premise

The observed asymmetric morphologies result from anisotropic diffusion whose shape is determined by the magnetic field viewing angle, Alfvénic Mach number, and pulsar distance.

What would settle it

A measurement showing that the Geminga and Monogem halos share the same mean magnetic field orientation would falsify the claim of separate coherence regions.

Figures

Figures reproduced from arXiv: 2603.07061 by Chao-Ming Li, Ruo-Yu Liu, Si-Zhe Wu.

Figure 1
Figure 1. Figure 1: Sketch figure for the definition of ϕ and χ. The location of B-field represents the projection of the mean mag￾netic field direction on the celestial sphere. Please refer to the appendix for the formula used to calculate the equatorial coordinates of the magnetic field direction. The injection term in Eq. (1), Q(Ee, t), is assumed to follow the spindown evolution of the pulsar, i.e., Q(Ee, t) = Q0 ( 1 + t … view at source ↗
Figure 2
Figure 2. Figure 2: Variation of the characteristic diffusion angle θd with azimuth ζ for Geminga and Monogem with different integration radius. The colored dashed line represents the characteristic diffusion angle at each azimuth, calculated from the median of the posterior distribution of the anisotropic model parameters. The solid colored line and the shaded region indicate the median and the 68% error range fitted for eac… view at source ↗
Figure 3
Figure 3. Figure 3: Spectra fitting to the spectra of Geminga (left panel) and Monogem (right panel), with blue and orange bands being uncertainty bands of 68% confidence intervals for the cases of rmax = 100 pc and 200 pc respectively. Gray bands show for comparison the results obtained by HAWC with the diffusion template (A. Albert et al. 2024). The black data points indicate the data used in the ηe fit, while the gray data… view at source ↗
Figure 4
Figure 4. Figure 4: The 90% confidence interval contour of the mag￾netic field direction distribution in the celestial sphere for Geminga and Monogem. The solid and dashed lines repre￾sent the results obtained with integration radii rmax of 100 pc and 200 pc, respectively. rmax = 100 pc and 200 pc around each pulsar. Notably, the two pulsars are spatially proximate, separated by ∼ 100 pc. Consequently, their emission regions … view at source ↗
Figure 5
Figure 5. Figure 5: The corner plots for Geminga (left) and Monogem (right) with integration radius rmax = 100 pc obtained using the MCMC method. Corner plots for the MCMC fitting in the case of rmax = 100 pc and rmax = 200 pc are shown in Figures 5 and 6, respectively. The fitting results for Monogem with rmax = 200 pc exhibit a bimodal distribution. For the peak with relatively smaller MA, θd in the third quadrant is obtain… view at source ↗
Figure 6
Figure 6. Figure 6: Same as [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
read the original abstract

Pulsar halos are produced by electrons and positrons diffusing in the interstellar medium around their parent pulsar wind nebulae. Recent observations by HAWC and LHAASO have revealed asymmetric morphologies in the halos surrounding Geminga and Monogem. The anisotropic diffusion model provides a natural explanation for such asymmetries, where the morphology is determined by the viewing angle of the mean magnetic field, the Alfv\'enic Mach number ($M_{\rm A}$), and the pulsar distance. In this work, we model the measured morphologies based on this framework and constrain the properties of interstellar magnetic turbulence. We find that the mean magnetic field orientations within the two halos are different, implying that they reside in different magnetic coherence regions, whereas the Alfv\'enic Mach numbers are relatively close ($M_{\rm A}\sim 0.2$). The results suggest a local magnetic field coherence length of approximately 100pc. Our study demonstrates that the morphology of pulsar halos serves as a powerful diagnostic tool for the properties of interstellar magnetic fields, highlighting the need for more accurate morphological measurements and sophisticated diffusion modeling in future studies.

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

0 major / 3 minor

Summary. The paper models the asymmetric morphologies of the Geminga and Monogem pulsar halos observed by HAWC as arising from anisotropic diffusion of electrons and positrons. The morphology depends on the viewing angle of the mean magnetic field, the Alfvénic Mach number M_A, and pulsar distance. Fitting yields different mean field orientations for the two halos (implying distinct magnetic coherence regions) but similar M_A ≈ 0.2, from which a local coherence length of ~100 pc is inferred. The work positions halo morphology as a diagnostic for interstellar magnetic turbulence.

Significance. If the fits are robust, the results supply an independent constraint on local ISM turbulence parameters (M_A and field orientation) that complements synchrotron and Faraday-rotation studies. The finding of comparable M_A despite differing orientations supports a picture of spatially varying but statistically similar turbulence regimes within ~100 pc of the Sun. This approach could be extended to additional HAWC/LHAASO halos for a statistical map of local magnetic coherence.

minor comments (3)
  1. [Abstract] The abstract states that morphologies are modeled and parameters constrained but supplies no reference to the diffusion tensor, likelihood function, or data-selection cuts; a one-sentence pointer to the relevant equation or section would clarify the fitting procedure for readers.
  2. [Results] The inference that the coherence length is ~100 pc appears to rest on the pulsars lying in different coherence regions; the manuscript should state the projected physical separation between Geminga and Monogem and the quantitative criterion used to convert orientation difference into a length scale.
  3. [Methods] Notation for the anisotropic diffusion coefficient (parallel vs. perpendicular) and the precise definition of M_A should be introduced with an equation in the methods section before the fitting results are presented.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the constructive and positive assessment of our manuscript, including the recognition of its potential to complement other probes of local ISM turbulence. The recommendation for minor revision is noted. No specific major comments were provided in the report, so we have no individual points requiring response or revision at this stage.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper applies an established anisotropic diffusion framework to fit observed HAWC morphologies of Geminga and Monogem halos, thereby constraining parameters (M_A ~0.2, differing field orientations, ~100 pc coherence length). This is ordinary parameter estimation from data; the abstract and claims contain no equations or steps in which a reported result is definitionally identical to a fitted input, a prediction is statistically forced by the fit itself, or a load-bearing premise reduces to a self-citation chain. The derivation remains self-contained against external morphological measurements.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the anisotropic diffusion model fully accounts for the observed asymmetry and that the three parameters (viewing angle, M_A, distance) are sufficient to determine morphology. No new entities are introduced. The fitted M_A and coherence length are results rather than free parameters supplied by the authors.

axioms (1)
  • domain assumption Asymmetric halo morphology is produced by anisotropic diffusion determined by viewing angle of mean B, Alfvénic Mach number, and pulsar distance.
    Stated in abstract paragraph 2 as the framework used to model the data.

pith-pipeline@v0.9.0 · 5746 in / 1366 out tokens · 31613 ms · 2026-05-25T07:06:09.369595+00:00 · methodology

Review history (2 revisions) →

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

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

  1. [1]

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

    Abeysekara, A. U., Albert, A., Alfaro, R., et al. 2017, Science, 358, 911, doi: 10.1126/science.aan4880

    work page doi:10.1126/science.aan4880 2017

  2. [2]

    X., et al

    Aharonian, F., An, Q., Axikegu, Bai, L. X., et al. 2021, Phys. Rev. Lett., 126, 241103, doi: 10.1103/PhysRevLett.126.241103

    work page doi:10.1103/physrevlett.126.241103 2021

  3. [3]

    C., et al

    Albert, A., Alfaro, R., Arteaga-Velázquez, J. C., et al. 2023, ApJL, 944, L29, doi: 10.3847/2041-8213/acb5ee

    work page doi:10.3847/2041-8213/acb5ee 2023

  4. [4]

    2024, The Astrophysical Journal, 974, 246, doi: 10.3847/1538-4357/ad738e

    Albert, A., Alfaro, R., Alvarez, C., et al. 2024, The Astrophysical Journal, 974, 246, doi: 10.3847/1538-4357/ad738e

    work page doi:10.3847/1538-4357/ad738e 2024

  5. [5]

    2024, Nuovo Cimento Rivista Serie, 47, 399, doi: 10.1007/s40766-024-00059-8

    Amato, E., & Recchia, S. 2024, Nuovo Cimento Rivista Serie, 47, 399, doi: 10.1007/s40766-024-00059-8

    work page doi:10.1007/s40766-024-00059-8 2024

  6. [6]

    2022, ApJ, 936, 183, doi: 10.3847/1538-4357/ac8b8a

    Bao, L.-Z., Fang, K., Bi, X.-J., & Wang, S.-H. 2022, ApJ, 936, 183, doi: 10.3847/1538-4357/ac8b8a

    work page doi:10.3847/1538-4357/ac8b8a 2022

  7. [7]

    2025a, Phys

    Bao, Y., Giacinti, G., Liu, R.-Y., Zhang, H.-M., & Chen, Y. 2025a, Phys. Rev. Lett., 135, 191002, doi: 10.1103/xqvb-lkrr

    work page doi:10.1103/xqvb-lkrr

  8. [8]

    2025b, Phys

    Bao, Y., Liu, R.-Y., Giacinti, G., Zhang, H.-M., & Chen, Y. 2025b, Phys. Rev. D, 112, 123017, doi: 10.1103/9461-gltv

    work page doi:10.1103/9461-gltv

  9. [9]

    2021, Nature, 594, 33

    Cao, Z., Aharonian, F., An, Q., et al. 2021, Nature, 594, 33

    work page 2021

  10. [10]

    2025, Science China

    Cao, Z., Aharonian, F., Axikegu, et al. 2025, Science China

    work page 2025

  11. [11]

    Physics, Mechanics, and Astronomy, 68, 279504, doi: 10.1007/s11433-024-2508-5

    work page doi:10.1007/s11433-024-2508-5

  12. [12]

    2023, in 38th International Cosmic Ray Conference (ICRC2023), Vol

    Chen, E., et al. 2023, in 38th International Cosmic Ray Conference (ICRC2023), Vol. ICRC2023, 613 9 De La Torre Luque, P., Fornieri, O., & Linden, T. 2022, PhRvD, 106, 123033, doi: 10.1103/PhysRevD.106.123033

    work page doi:10.1103/physrevd.106.123033 2023

  13. [13]

    2022, Frontiers in Astronomy and Space Sciences, 9, 1022100, doi: 10.3389/fspas.2022.1022100

    Fang, K. 2022, Frontiers in Astronomy and Space Sciences, 9, 1022100, doi: 10.3389/fspas.2022.1022100

    work page doi:10.3389/fspas.2022.1022100 2022

  14. [14]

    PHECT: A lightweight computation tool for pulsar halo emission

    Fang, K. 2025, arXiv e-prints, arXiv:2508.13667, doi: 10.48550/arXiv.2508.13667

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2508.13667 2025

  15. [15]

    2019, Mon

    Fang, K., Bi, X.-J., & Yin, P.-F. 2019, Mon. Not. R. Astron. Soc., 488, 4074, doi: 10.1093/mnras/stz1974

    work page doi:10.1093/mnras/stz1974 2019

  16. [16]

    2023, Phys

    Fang, K., Hu, H.-B., Bi, X.-J., & Chen, E.-S. 2023, Phys. Rev. D, 108, 023017, doi: 10.1103/PhysRevD.108.023017

    work page doi:10.1103/physrevd.108.023017 2023

  17. [17]

    W., Lang D., Goodman J., 2013, @doi [ ] 10.1086/670067 , http://adsabs.harvard.edu/abs/2013PASP..125..306F 125, 306

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, pasp, 125, 306, doi: 10.1086/670067 H. E. S. S. Collaboration, Aharonian, F., Ait Benkhali, F., et al. 2023, A&A, 673, A148, doi: 10.1051/0004-6361/202245776

    work page doi:10.1086/670067 2013

  18. [18]

    2017, Phys

    Hooper, D., Cholis, I., Linden, T., & Fang, K. 2017, Phys. Rev. D, 96, 103013, doi: 10.1103/PhysRevD.96.103013 Jóhannesson, G., Porter, T. A., & Moskalenko, I. V. 2019, ApJ, 879, 91, doi: 10.3847/1538-4357/ab258e

    work page doi:10.1103/physrevd.96.103013 2017

  19. [19]

    2026, Phys

    Li, C.-M., Liu, R.-Y., & Yan, H. 2026, Phys. Rev. D, 113, L021303, doi: 10.1103/vt3s-rbj1

    work page doi:10.1103/vt3s-rbj1 2026

  20. [20]

    2022, International Journal of Modern Physics A, 37, 2230011, doi: 10.1142/S0217751X22300113

    Liu, R.-Y. 2022, International Journal of Modern Physics A, 37, 2230011, doi: 10.1142/S0217751X22300113

    work page doi:10.1142/s0217751x22300113 2022

  21. [21]

    2019, Phys

    Liu, R.-Y., Yan, H., & Zhang, H. 2019, Phys. Rev. Lett., 123, 221103, doi: 10.1103/PhysRevLett.123.221103 López-Coto, R., de Oña Wilhelmi, E., Aharonian, F.,

    work page doi:10.1103/physrevlett.123.221103 2019

  22. [22]

    2022, Nature Astronomy, 6, 199, doi: 10.1038/s41550-021-01580-0 López-Coto, R., & Giacinti, G

    Amato, E., & Hinton, J. 2022, Nature Astronomy, 6, 199, doi: 10.1038/s41550-021-01580-0 López-Coto, R., & Giacinti, G. 2018, MNRAS, 479, 4526, doi: 10.1093/mnras/sty1821

    work page doi:10.1038/s41550-021-01580-0 2022

  23. [23]

    N., Hobbs G

    Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, The Astronomical Journal, 129, 1993, doi: 10.1086/428488

    work page internal anchor Pith review doi:10.1086/428488 2005

  24. [24]

    2023, Mon

    Pavaskar, P., Yan, H., & Cho, J. 2023, Mon. Not. R. Astron. Soc., 523, 1056, doi: 10.1093/mnras/stad1237

    work page doi:10.1093/mnras/stad1237 2023

  25. [25]

    2018, PhRvD, 97, 123008, doi: 10.1103/PhysRevD.97.123008

    Profumo, S., Reynoso-Cordova, J., Kaaz, N., & Silverman, M. 2018, PhRvD, 97, 123008, doi: 10.1103/PhysRevD.97.123008

    work page doi:10.1103/physrevd.97.123008 2018

  26. [26]

    A., et al

    Recchia, S., Di Mauro, M., Aharonian, F. A., et al. 2021, Physical Review D, 104, doi: 10.1103/physrevd.104.123017

    work page doi:10.1103/physrevd.104.123017 2021

  27. [27]

    2019, MNRAS, 484, 3491, doi: 10.1093/mnras/stz268

    Tang, X., & Piran, T. 2019, MNRAS, 484, 3491, doi: 10.1093/mnras/stz268

    work page doi:10.1093/mnras/stz268 2019

  28. [28]

    V., et al

    Trotta, R., Jóhannesson, G., Moskalenko, I. V., et al. 2011, The Astrophysical Journal, 729, 106, doi: 10.1088/0004-637X/729/2/106

    work page doi:10.1088/0004-637x/729/2/106 2011

  29. [29]

    2024, ApJ, 969, 9, doi: 10.3847/1538-4357/ad43e1

    Wu, Q.-Z., Li, C.-M., Liang, X.-H., Ge, C., & Liu, R.-Y. 2024, ApJ, 969, 9, doi: 10.3847/1538-4357/ad43e1

    work page doi:10.3847/1538-4357/ad43e1 2024

  30. [30]

    2008, The Astrophysical Journal, 673, 942, doi: 10.1086/524771

    Yan, H., & Lazarian, A. 2008, The Astrophysical Journal, 673, 942, doi: 10.1086/524771

    work page doi:10.1086/524771 2008

  31. [31]

    Z., & Wang, X.-Y

    Yan, K., Liu, R.-Y., Chen, S. Z., & Wang, X.-Y. 2022, The Astrophysical Journal, 935, 65, doi: 10.3847/1538-4357/ac7feb

    work page doi:10.3847/1538-4357/ac7feb 2022

  32. [32]

    2025, The Astrophysical Journal, 987, 19, doi: 10.3847/1538-4357/add6a4

    Yan, K., Wu, S., & Liu, R.-Y. 2025, The Astrophysical Journal, 987, 19, doi: 10.3847/1538-4357/add6a4

    work page doi:10.3847/1538-4357/add6a4 2025