Pith. sign in

REVIEW 3 major objections 4 minor 24 references

The Milky Way's TeV gamma-ray glow extends far beyond the narrow galactic plane, according to a new analysis of eight years of HAWC data.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

HAWC's Pass 5 analysis of 8 years of data yields new TeV diffuse galactic emission profiles with tails beyond 3 degrees latitude and a galactic-center excess relative to standard cosmic-ray propagation models.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection Useful preliminary HAWC profiles of TeV diffuse emission; new Pass 5 data and a fair LHAASO comparison, but the 'significant' claims need numbers and the source-subtraction systematics are only partially bracketed. the 3 major comments →

arxiv 2509.03189 v1 pith:EAE4FYHR submitted 2025-09-03 astro-ph.HE

Dissecting the Diffuse Emission of the Galaxy with the HAWC Observatory

Georg Schwefer (for the HAWC Collaboration) This is my paper
classification astro-ph.HE
keywords Galactic diffuse gamma-ray emissionTeV gamma rayscosmic-ray propagationHAWCGalactic plane surveysource subtractionLHAASO comparisoninterstellar gas
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

The paper reports that the Milky Way's diffuse TeV gamma-ray emission—produced when cosmic rays collide with interstellar gas—is not confined to the thin galactic plane. After fitting and subtracting a catalog of sources, significant residual emission remains all along the plane and extends in broad, non-gaussian tails out to at least 10 degrees of galactic latitude. The shape of those tails matches expectations for cosmic rays interacting with atomic gas that rises hundreds of parsecs from the plane. Toward the galactic center, the measured flux exceeds predictions from conventional propagation models, hinting at extra cosmic-ray sources, a radial gradient, or unresolved faint sources. An independent analysis that mimics a recently published LHAASO-style mask recovers consistent spectra, lending confidence to the result.

Core claim

Using eight years of reprocessed Pass 5 data from a high-altitude water-Cherenkov observatory, the authors subtract a fitted source model—built with the automated workflow planned for a future Galactic Plane Survey—from maps of the TeV sky. The leftover emission is then fitted with three templates for the diffuse emission. The central finding is that this residual diffuse emission is significant across the Galactic plane and shows non-gaussian latitudinal tails out to |b| = 10 degrees, consistent with gas-correlated emission. The overall flux is higher than predicted by conventional models such as CRINGE, with the excess strongest toward the inner Galaxy. Spectral fits in inner and outer gal

What carries the argument

The source-subtraction workflow is the central mechanism: candidate sources are identified from significance maps at two correlation radii, merged if overlapping, and fit with generalized-Gaussian spatial models and log-parabola spectra, then subtracted from the data. The remaining emission is modeled with three diffuse templates—the CRINGE hadronic model with rescaled CO-to-H2 conversion, a Fermi-tuned variable-diffusion model, and a uniform-cosmic-ray model—and the spread across six fits (two source models, three templates) defines the systematic uncertainty bands.

Load-bearing premise

The fitted source catalog fully accounts for every compact and extended object in the field, so that whatever remains after subtraction is genuinely diffuse emission rather than misattributed source flux.

What would settle it

A future instrument with better angular resolution could test whether the apparent diffuse emission at |b| > 3 degrees resolves into known or new point sources; if so, the gas-correlated diffuse claim weakens. Concretely, running this source-subtraction pipeline on mock skymaps with a known diffuse component plus realistic unresolved-source populations would reveal whether the recovered latitudinal tails are biased by source leakage.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Diffuse TeV emission extends at least 10 degrees from the Galactic plane, with non-gaussian tails whose shape matches cosmic-ray interaction with interstellar gas.
  • The excess toward the Galactic center is larger than conventional models predict, implying either a more centrally peaked cosmic-ray source distribution, a radial gradient in accelerator properties, non-standard propagation, or unresolved faint sources.
  • The HAWC spectra in the LHAASO-defined inner and outer windows are consistent with LHAASO results where the analysis methods are matched, supporting the robustness of the diffuse-emission measurement.
  • The presence of the |b| > 3 degree tails is independent of the source model assumed, strengthening the claim that the extended emission is genuine.
  • The latitudinal and longitudinal flux profiles provide a new benchmark for TeV Galactic diffuse models, comparable in scope to what Fermi-LAT established at GeV energies.

Where Pith is reading between the lines

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

  • If the inner-Galaxy excess is partly unresolved sources, a deeper survey with higher angular resolution should resolve a portion of the reported diffuse flux into point sources, shifting the excess downward.
  • The latitudinal tails could be used to constrain the scale height of cosmic rays and gas in the inner Galaxy when combined with atomic and molecular gas maps, offering a direct probe of propagation perpendicular to the disk.
  • The same source-subtraction pipeline could be applied to simulated skymaps with a known diffuse component to calibrate the systematic bias introduced by source-diffuse confusion, providing a quantitative correction for future analyses.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This conference paper presents an 8-year HAWC (Pass 5) analysis of TeV Galactic diffuse gamma-ray emission. The authors subtract a catalog of point-like/extended sources using a CTAO-style detection and fitting algorithm, then fit the residual emission with three Galactic diffuse emission templates (CRINGE with rescaled X_CO, IEM-varmin-rescaled, and Uniform CR). They report longitudinal and latitudinal flux profiles with statistical, source-model, and template uncertainties, and find significant remaining emission throughout the Galactic plane, non-Gaussian latitude tails out to |b|=10 deg, and an excess toward the Galactic center relative to the CRINGE prediction. They also reproduce the LHAASO analysis in inner and outer windows and report broadly consistent spectra, with a ~50% higher inner-window normalization attributed to mask differences.

Significance. If the claims hold, this would be an important, independent TeV measurement of the Galactic diffuse emission from HAWC, especially the latitude tails beyond 3 deg and the inner-Galaxy excess. The analysis is more careful than typical conference proceedings: it uses state-of-the-art data processing (Pass 5), an open-source tool (Gammapy), and explicitly treats three classes of uncertainty (statistical, source model, template). The cross-check against LHAASO with an equivalent analysis method is a useful consistency test. However, the central scientific claims rest on qualitative significance assertions and on a source-subtraction systematic that is only partially bracketed, so the current manuscript is not yet suitable as a definitive measurement.

major comments (3)
  1. [Section 3.1, Figs. 2 and 3] The abstract and Section 3.1 assert 'significant remaining emission throughout the Galactic plane' and 'significant emission beyond 3 degree latitude', but no test statistic, p-value, or significance map for the residual after source subtraction is provided. The lower panels of Fig. 1 are residual significance maps, but they are not per-longitude-bin or per-latitude-bin significances and do not quantify the tail claim. Please define 'significant' and report a quantitative measure (e.g., TS or p-value) for the >3 deg emission and for the inner-Galaxy excess relative to the models.
  2. [Section 2.2 and 3.1] The source-model systematic is estimated from only two selection thresholds (all 3-sigma sources vs. sources with >=4 sigma or catalog association). This does not bracket the dominant degeneracy: the source-detection map already contains the diffuse emission, so compact diffuse structures can be absorbed as 'sources', while imperfectly modeled bright sources (as acknowledged for Crab and Geminga) can inject residual flux on degree scales. The claim that 'the presence of the tails ... is independent of the source model assumed' is therefore not established. A cross-check using a mask-based source exclusion (as in the LHAASO comparison) or an explicit unresolved-source population model would materially strengthen the paper.
  3. [Section 3.1, Fig. 3] The latitude-tail claim is based on a single longitude window, 20 deg < l < 60 deg, yet the abstract and conclusions state it as a general property of the Galactic diffuse emission. The longitude profile is integrated over |b|<10 deg and cannot confirm the tail beyond 3 deg elsewhere. Either restrict the claim to the measured range or show similar profiles for additional longitude bins (e.g., inner Galaxy and Cygnus), where source confusion and exposure differ significantly.
minor comments (4)
  1. [Section 1] Typos: 'singificant' should be 'significant'; 'Miky Way' should be 'Milky Way'.
  2. [Fig. 1 caption and Section 2.2] The phrase 'multiplied in each longitude bin by the best-fit normalization' is ambiguous. Clarify whether the normalization is an independent fit parameter per bin or a smooth function, and how this relates to the profile shown in Fig. 2.
  3. [Tables 1 and 2] In Table 2, for the single power-law fit in the outer window, the columns gamma2 and E_br should be omitted or marked with an em dash rather than left blank, to avoid visually implying a broken power law.
  4. [Section 3.2] The 50% higher inner-window normalization is attributed to mask differences, but no quantitative test is shown. If the [8] mask is not available, a statement that this is a plausible explanation rather than a demonstrated one would be more precise.

Circularity Check

0 steps flagged

No significant circularity: the central diffuse-emission measurement is a source-subtracted residual compared with independent templates and an external LHAASO cross-check; the main caveats are source/diffuse degeneracy and unquantified significance, not circular construction.

full rationale

The derivation chain is not circular in the sense of fitting a parameter and renaming it a prediction, nor does it reduce to a self-citation theorem. The source model is subtracted first, and the remaining GDE flux is then fitted with three external template models (CRINGE, Uniform CRs, IEM-varmin-rescaled). The headline claims—significant remaining emission throughout the Galactic plane and non-Gaussian tails beyond |b|>3°—are statements about the residual data after source subtraction, not about quantities encoded in the fit. The CRINGE model is cited from the same author's prior work (ref. [13]), but it is only one of three comparison templates and is not used as a uniqueness argument or as the sole support for any claim; the LHAASO-consistency analysis provides an independent benchmark. The paper explicitly acknowledges the real systematic degeneracy between large, low-surface-brightness sources and diffuse emission, and imperfect source fits near Crab and Geminga; these are honest limitations and correctness risks, not circularity. The per-longitude-bin normalizations mentioned in the Figure 1 caption are used to display residuals and construct profiles; however, the shape comparisons in the latitudinal profile and the external LHAASO comparison do not reduce to this fit. The word 'significant' is not quantified, but that is a statistical-reporting issue rather than a circularity issue. Overall, the central result retains independent content and is not forced by definition or by the author's own prior models.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central results rest on the instrument response and background model (Pass 5), the template models for the diffuse emission, and the completeness of the source-subtraction algorithm. The fitted source and template parameters are nuisance parameters; the key numbers in the paper are the residual profiles and their comparisons.

free parameters (3)
  • Source model parameters (spatial and spectral) for all fitted sources = not tabulated; fit per candidate
    The CTAO-style source detection and fitting algorithm assigns each candidate a generalized Gaussian spatial model and log-parabola spectrum; these parameters are fitted to the same data from which the residual diffuse emission is derived, so the residual depends on their values.
  • Diffuse template normalization per longitude/latitude region = not tabulated; best-fit normalization per bin
    For each of the three GDE models (CRINGE rescaled XCO, IEM-varmin rescaled, Uniform CR), the model is scaled to the residual data in each region/bin; the lower panel of Figure 1 uses the best-fit normalization per longitude bin.
  • Source detection thresholds (correlation radii 0.2 and 0.4 deg, 3 sigma threshold, 0.2 deg separation, 0.5 deg radius di = chosen by hand
    These algorithmic choices determine which sources are subtracted and thus affect the residual; the source-model uncertainty is estimated by comparing two threshold variants.
axioms (4)
  • domain assumption The HAWC Pass 5 event reconstruction, energy estimation (neural network), and gamma-hadron separation are correct within stated resolutions.
    Section 2.1: cuts on PSF containment, energy bias, and resolution rely on the Pass 5 instrument response being accurate.
  • domain assumption The hadronic background model after masking |b| < 10 deg correctly estimates cosmic-ray background; no significant diffuse gamma emission is absorbed by the background fit.
    Section 2.1: 'We mask out the Galactic plane... to avoid any diffuse emission being absorbed in the background.' If the mask is insufficient, the residual diffuse flux would be biased.
  • domain assumption The three GDE templates (CRINGE rescaled XCO, IEM-varmin-rescaled, Uniform CR) provide valid spatial and spectral shapes for the TeV diffuse emission, including gas distributions from external surveys.
    Section 2.2: the residual data are fitted with these templates; the 'diffuse model uncertainty' band only spans these three models and does not cover other gas maps or CR propagation scenarios.
  • domain assumption The CTAO-style source detection and fitting algorithm identifies essentially all non-diffuse gamma-ray sources contributing in the analysis region.
    Section 2.2 and 3.1: residuals after source subtraction are interpreted as GDE; the paper acknowledges imperfect fits near Crab, Geminga, Cygnus, and inner Galaxy, so a violation would directly bias the profiles.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Dissecting the Diffuse Emission of the Galaxy with the HAWC Observatory." pith.science (2026). https://pith.science/paper/EAE4FYHR

@misc{pith2026250903189,
  author       = {Pith},
  title        = {Pith review of: Dissecting the Diffuse Emission of the Galaxy with the HAWC Observatory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EAE4FYHR}},
  note         = {Machine review of arXiv:2509.03189}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Galactic diffuse gamma-ray emission is produced by the interaction of high-energy cosmic rays propagating through the Milky Way with interstellar gas and radiation fields. Its measurement can provide crucial insights into the acceleration and transport of cosmic rays throughout our Galaxy. Here, we present a new analysis of the TeV Galactic diffuse gamma-ray emission using 8 years of HAWC data. This data was processed with the updated Pass 5 processing, enhancing the sensitivity and resolution of the instrument. For the analysis, we make use of Gammapy, an open-source package for gamma-ray astronomy, and recent models of the Galactic diffuse emission at TeV energies. After subtracting the emission from sources using an algorithm akin to that developed for the foreseen CTAO Galactic plane survey, we find significant remaining emission throughout the Galactic plane. We show the latitudinal and longitudinal flux profiles of the emission in multiple parts of the galaxy taking into account various sources of uncertainty and compare to existing models. We find significant emission beyond 3 degree latitude, consistent in shape with the prediction for the interaction of cosmic rays with the interstellar gas. We also demonstrate that our results are consistent with recent LHAASO results when equivalent analysis methods are used.

Figures

Figures reproduced from arXiv: 2509.03189 by Georg Schwefer (for the HAWC Collaboration).

Figure 1
Figure 1. Figure 1: Significance maps of our dataset with a correlation radius of 0.4 ◦ between −10◦ < 𝑏 < 10◦ and 0 ◦ < 𝑙 < 240◦ . In the upper panel, we show the raw significance map including all sources and diffuse emission. The middle panel shows the residual significance after subtraction of the source model. The lower panel shows the residual significance after further subtraction of the CRINGE diffuse emission model, … view at source ↗
Figure 2
Figure 2. Figure 2: Longitudinal profile in integrated flux between 1 and 100 TeV of the residual Galactic emission after subtraction of the source model between −10◦ < 𝑏 < 10◦ and 0 ◦ < 𝑙 < 240◦ . We show three different uncertainties associated with the measurement as well as three different models for the GDE of the Milky Way described in the text. −10 −5 0 5 10 Galactic Latitude in deg 0.00 0.25 0.50 0.75 1.00 1.25 1.50 R… view at source ↗
Figure 4
Figure 4. Figure 4: Spectrum of the Galactic diffuse gamma-ray emission in the inner (left) and outer (right) Galactic window obtained from the analysis following the LHAASO methods. We also show the corresponding LHAASO results from [8]. The respective best-fit parameters and uncertainties are given in Tables 1 and 2. 4. Summary and Conclusions We have presented an updated and expanded analysis of the TeV Galactic diffuse ga… view at source ↗

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

24 extracted references · 8 canonical work pages · 2 internal anchors

  1. [1]

    Grenier et al

    I.A. Grenier et al. ARA&A53 (2015) 199

    work page 2015

  2. [2]

    Tibaldo et al.Universe7 (2021) 141 [2103.16423]

    L. Tibaldo et al.Universe7 (2021) 141 [2103.16423]

    Pith/arXiv arXiv 2021

  3. [3]

    Fermi-LATcollaboration Astrophys. J. Suppl.223(2016) 26 [1602.07246]

    Pith/arXiv arXiv 2016

  4. [4]

    H.E.S.S.collaboration Phys. Rev. D90(2014) 122007 [1411.7568]

    Pith/arXiv arXiv 2014

  5. [5]

    IceCubecollaboration Science380 (2023) 1338

    work page 2023

  6. [6]

    Tibet ASgamma collaborationPhys. Rev. Lett.126 (2021) 141101 [2104.05181]

    Pith/arXiv arXiv 2021

  7. [7]

    LHAASO collaborationPhys. Rev. Lett.131 (2023) 151001 [2305.05372]

    arXiv 2023

  8. [8]

    LHAASO collaborationPhys. Rev. Lett.134 (2025) 081002 [2411.16021]

    arXiv 2025

  9. [9]

    J.961(2024) 104 [2310.09117]

    HA WCcollaborationAstrophys. J.961(2024) 104 [2310.09117]

    work page arXiv 2024

  10. [10]

    He et al.Astrophys

    J. He et al.Astrophys. J.980 (2025) 17 [2501.02720]

    Pith/arXiv arXiv 2025

  11. [11]

    De La Torre Luque et al.2502.18268

    P. De La Torre Luque et al.2502.18268

    Pith/arXiv arXiv

  12. [12]

    Castro et al.2506.06593

    L.E.E. Castro et al.2506.06593

    Pith/arXiv arXiv

  13. [13]

    Schwefer et al.Astrophys

    G. Schwefer et al.Astrophys. J.949 (2023) 16 [2211.15607]

    Pith/arXiv arXiv 2023

  14. [14]

    De La Torre Luque et al.Astron

    P. De La Torre Luque et al.Astron. Astrophys.672(2023) A58 [2203.15759]

    Pith/arXiv arXiv 2023

  15. [15]

    Performance of the HAWC Observatory and TeV Gamma-Ray Measurements of the Crab Nebula with Improved Extensive Air Shower Reconstruction Algorithms

    HA WCcollaboration Astrophys. J.972(2024) 144 [2405.06050]

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  16. [16]

    J.881(2019) 134 [1905.12518]

    HA WCcollaboration Astrophys. J.881(2019) 134 [1905.12518]

    Pith/arXiv arXiv 2019

  17. [17]

    Validation of standardized data formats and tools for ground-level particle-based gamma-ray observatories

    HA WCcollaboration Astron. Astrophys.667 (2022) A36 [2203.05937]

    work page internal anchor Pith review Pith/arXiv arXiv 2022

  18. [18]

    Astrophys.678 (2023) A157 [2308.13584]

    Gammapycollaboration Astron. Astrophys.678 (2023) A157 [2308.13584]

    Pith/arXiv arXiv 2023

  19. [19]

    Acero et al

    F. Acero et al. Jan., 2025. 10.5281/zenodo.14760974

    doi:10.5281/zenodo.14760974 2025

  20. [20]

    CTA ConsortiumcollaborationJCAP10(2024) 081 [2310.02828]

    Pith/arXiv arXiv 2024

  21. [21]

    Pagliaroli et al.Universe9 (2023) 381

    G. Pagliaroli et al.Universe9 (2023) 381

    work page 2023

  22. [22]

    Kalberla and J

    P.M.W. Kalberla and J. KerpARA&A 47(2009) 27

    work page 2009

  23. [23]

    Bull.69(2024) 449 [2310.10100]

    LHAASO collaborationSci. Bull.69(2024) 449 [2310.10100]

    arXiv 2024

  24. [24]

    J.974(2024) 246

    HA WCcollaboration Astrophys. J.974(2024) 246. 8 Dissecting the Diffuse Emission of the Galaxy with the HAWC Observatory Georg Schwefer Full Author List: HAWC Collaboration R.Alfaro 1,C.Alvarez 2,A.Andrés 3,E.Anita-Rangel 3,M.Araya 4,J.C.Arteaga-Velázquez 5,D.AvilaRojas 3,H.A.AyalaSolares 6, R. Babu7, P. Bangale8, E. Belmont-Moreno1, A. Bernal3, K.S. Caba...

    work page 2024

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.