Pith. sign in

REVIEW 4 major objections 6 minor 14 references

Studying Cosmic-ray Interactions in Giant Molecular Clouds with the HAWC Gamma-ray Observatory

T0 review · 4 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read This HAWC search for TeV gamma-ray emission from three nearby molecular clouds finds no significant excess and places 95% upper limits above the expected hadronic cosmic-ray flux.

desk verdict Honest, well-scoped HAWC upper limits for three GMCs; the missing threshold-stability and systematics work keeps the conclusion provisional, but it's a useful calibration measurement. read the letter →

arxiv 1908.06073 v1 pith:WOMSDRV5 submitted 2019-08-16 astro-ph.HE

Hugo Alberto Ayala Solares (for the HAWC Collaboration) This is my paper · ORCID
classification astro-ph.HE
keywords cosmicraysgamma-rayastronomygiantmolecularcloudsTeVupperlimitsHAWCGalacticcosmic-rayseahadronicinteractionsdustopticaldepthtemplates
verification ladder T0 review T1 audit T2 compute T3 formal

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 tests whether the cosmic-ray flux measured near Earth is the same as the “sea” of Galactic cosmic rays by looking for the gamma rays those cosmic rays should produce when they collide with gas in three nearby giant molecular clouds: Aquila Rift, Hercules, and Taurus. Using about 1127 days of HAWC observatory data, the search finds no significant gamma-ray excess from any of the three clouds. The paper therefore reports 95% credible upper limits on the TeV gamma-ray flux in four energy bins between 1 and 100 TeV. Those upper limits sit above the flux expected from purely hadronic collisions of the local cosmic-ray spectrum with the cloud gas, which is below $10^{-11}\,\mathrm{TeV}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$ above 10 TeV. The result is consistent with the uniform cosmic-ray sea paradigm but does not confirm it; the paper is explicit that the results are preliminary.

What carries the argument

The spatial templates for the three clouds are built from Planck dust optical depth maps at 353 GHz, converted to molecular-hydrogen column density through $N_{\mathrm{H_2}} = \tau_D / (\tau_D/N_{\mathrm{H_2}})_{\mathrm{ref}}$ with a reference value of $1.18\times 10^{-26}\,\mathrm{cm}^2$, followed by opacity cuts ($5\times 10^{-5}$ for Taurus and Aquila, $2.5\times 10^{-5}$ for Hercules) and normalization to the mean column density. The gamma-ray expectation is computed by convolving the locally measured cosmic-ray spectrum with a pion-to-gamma production cross section, using both a column-density route (Eq. 3.1) and a mass-and-distance route (Eq. 3.2). A likelihood-ratio test statistic is maximized under a fixed power-law spectrum with index 2.75, and a Markov-Chain Monte Carlo with a uniform prior on normalization produces the 95% credible intervals when the test statistic stays below 25.

What would settle it

A TeV observation with substantially better sensitivity that resolves one of these clouds and measures a gamma-ray flux above the 95% upper limits reported in Table 3 would contradict the reported limits, while a measured flux significantly below the hadronic expectation would contradict the uniform cosmic-ray-sea paradigm the paper sets out to test.

Watch

Extended reading notes

Core claim

No significant TeV gamma-ray excess is observed from Aquila Rift, Hercules, or Taurus in HAWC data, so the analysis places 95% credible upper limits on the gamma-ray flux in four half-decade energy bins from 1 to 100 TeV. The expected gamma-ray flux from pure hadronic interactions of the local cosmic-ray “sea” with the passive molecular clouds is below $10^{-11}\,\mathrm{TeV}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$ above 10 TeV, and the measured upper limits lie above that expected range (see Table 3 and Figure 2). Thus the data cannot distinguish the hadronic expectation from background, and the paradigm that the locally measured cosmic-ray flux represents the Galactic sea remains viable but unproven at TeV energies.

Load-bearing premise

The result depends on the Planck dust optical depth maps and the chosen opacity thresholds tracing exactly the molecular gas that would emit gamma rays; if dust does not trace the gas, the spatial templates and cloud masses, and therefore the derived upper limits, no longer correspond to the clouds as claimed.

Editorial extensions

If this is right

  • The reported upper limits quantify the TeV gamma-ray flux from three high-latitude molecular clouds in four energy bins and can be compared directly with hadronic models of cosmic-ray interactions.
  • If the local cosmic-ray sea is indeed uniform, longer HAWC exposures or a more sensitive TeV instrument should eventually detect these clouds at or near the expected hadronic flux level.
  • The non-detection places no pressure on the standard paradigm, but it also does not verify it; the allowed flux range still includes both the hadronic expectation and zero.
  • The same template-and-likelihood pipeline can be applied to other molecular clouds outside the Galactic plane to build a larger sample of TeV upper limits.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: because the upper limits sit above the hadronic expectation, combining many high-latitude clouds could push the effective sensitivity down to the predicted flux and provide a sharper test of the uniform cosmic-ray sea.
  • Editorial inference: the strong dependence on Planck dust opacity and the chosen thresholds means that improved dust-to-gas calibration, or comparison with CO-based templates, would likely tighten or shift these limits; this is a testable extension rather than a stated result.
  • Editorial inference: if the local cosmic-ray sea is not uniform, high-latitude clouds at different distances and Galactic latitudes—such as these three—are natural probes of small-scale anisotropies in the TeV cosmic-ray density.
  • Editorial inference: the quoted expected flux assumes the locally measured spectrum extends unchanged to the clouds; any future direct measurement of a spectral break in the 1–100 TeV range would change the predicted flux and could make the current limits more or less constraining.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 6 minor

Summary. The paper reports a search for TeV gamma-ray emission from three giant molecular clouds (Taurus, Aquila Rift, Hercules) using 1127 days of HAWC data. The analysis builds spatial templates from Planck 353 GHz dust optical depth maps, assumes a power-law spectrum with fixed index alpha = 2.75 in four half-decade energy bins from 1 to 100 TeV, and uses a maximum-likelihood/MCMC procedure to derive 95% credible upper limits. No significant excess (TS < 25) is found in any cloud. The limits are compared with two model predictions for hadronic gamma-ray emission obtained by convolving the AMS-measured local cosmic-ray flux with the Kafexhiu et al. pion-production parametrization, using either the Planck column density or the cloud mass and distance. The authors conclude that the measured upper limits lie above the expected hadronic flux, so the local cosmic-ray 'sea' paradigm cannot be constrained by the present sensitivity.

Significance. The measurement addresses a question of genuine interest: whether the local cosmic-ray spectrum is representative of the Galactic cosmic-ray sea, using high-latitude giant molecular clouds as targets in the TeV band. HAWC's wide field of view and high duty cycle are well matched to these extended, low-surface-brightness sources. The upper-limit machinery is standard and internally consistent, and the comparison with expectation is anchored to an external cosmic-ray spectrum (AMS) and an independent pion-production parametrization, so there is no statistical circularity in the central comparison. The paper also makes the quasi-differential limits and cloud properties available for future joint analyses. Its significance is limited by the fact that the upper limits are above the expected fluxes, so the result is a null measurement with modest constraining power, and by the lack of a systematic treatment of the template definition and spectral model, which must be addressed before the 'upper limit above expectation' conclusion can be taken as robust.

major comments (4)
  1. [Sec. 2.1, Eq. (2.3), Fig. 1] The spatial templates are defined by hand-chosen Planck opacity thresholds (5e-5 for Taurus and Aquila, 2.5e-5 for Hercules) with no stability test and no comparison with an independent tracer such as CO. The same Planck-based map enters both the search template and the expected-flux models (Eqs. 3.1-3.2), so the threshold choice couples the observed upper limits to the predicted flux in an uncontrolled way. A lower threshold may add unrelated HI or dust cirrus, while a higher threshold may exclude diffuse molecular gas and lower the expected flux, making the conclusion that the limits are above the expectation easier to reach. Please include a threshold scan (e.g., varying tau_cut by factors of two) with the resulting limits and expected fluxes, or a CO-based cross-check for at least one cloud, and quote the corresponding systematic spread in Table 3.
  2. [Sec. 2, Eq. (2.2)] The spectral index of the assumed power law is fixed at alpha = 2.75, but the paper does not test the dependence of the quasi-differential upper limits on this choice. Because the normalization K is fitted with the spectral shape fixed inside each energy bin, a harder or softer spectrum within the plausible hadronic range (e.g., alpha = 2.3-3.0) can shift the derived K and therefore the limits in Table 3. Please provide a scan over alpha, or an argument that 2.75 is conservative in every energy bin, before the limits are used as a test of the cosmic-ray-sea model.
  3. [Sec. 3, Eqs. (3.1)-(3.2), Table 2, Fig. 2] The expected-model bands in Fig. 2 do not include the uncertainties in the dust-to-gas conversion factor (tau_D/NH2)_ref, the opacity-threshold choice, or the distance and mass determinations quoted in Table 2. The green band is only the 10th-90th percentile spread of the column-density distribution and does not capture the absolute normalization uncertainty that enters Eq. (3.1) through NH2 and Eq. (3.2) through M/D^2. A systematic error budget for the expected flux, for example a +/-30% variation of the conversion factor combined with the distance uncertainties, is needed to support the statement that the upper limits lie above the expected hadronic range.
  4. [Sec. 2, Table 3] The 95% credible intervals in Table 3 are statistical only. The paper does not discuss the dominant experimental systematics in the HAWC energy scale, angular resolution, background model, or the fhit binning, which for a TeV measurement can be comparable to the statistical uncertainties. A summary of the dominant systematic uncertainties is necessary for a journal-level upper-limit result, and the absence of such a discussion currently limits the robustness of the central claim.
minor comments (6)
  1. [Sec. 2.1] The sentence 'Then we normalize the map to the mean column density of the cloud times the size of the spatial bin (1/(⟨NH2⟩ dΩ))' is unclear; please define the normalization factor explicitly and state the units of the resulting template.
  2. [Sec. 2, Eq. (2.1)] The test statistic is missing the logarithm; it should read TS = 2 ln[L(S+B)/L(B)], or the symbols should be defined as log-likelihoods rather than likelihoods.
  3. [Table 1] The 'midpoint' values (1.77, 5.59, 17.7, 55.97) are geometric means of the bin edges; please label them as 'logarithmic midpoint' or 'geometric mean energy' to avoid confusion.
  4. [Table 2] The column labeled 'Size [sr]' is a solid angle, not a physical size; rename it to 'Solid angle Ω [sr]' for consistency with the text.
  5. [Figure 2] The caption says the band represents the uncertainty in the upper limit, while the text describes the green band as the 10th-90th percentile of the column density distribution; please clarify which quantity the bands represent.
  6. [References and text] Reference [14] is incomplete ('Schlafly, ApJ (2014)') and should be completed; also fix the typo 'has already being exploited' to 'has already been exploited'.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: the CR-sea expectation is anchored to external AMS and cross-section inputs, while the HAWC upper limits are data-derived; the only self-citations are calibration references that do not carry the conclusion.

full rationale

The central comparison is independent of the fitted HAWC parameters. The predicted hadronic flux is obtained by convolving the AMS local CR spectrum with the Kafexhiu et al. pion-production parameterization and scaling by cloud column density or mass/distance from Planck dust maps and Schlafly distances. The HAWC upper limits come from maximizing the likelihood in each fixed energy bin and taking the 95% credible interval; no parameter of the CR-sea model is fitted to HAWC data, and the expected flux is not derived from the fitted normalization K. The Planck-based template is an input to the upper-limit calculation, so the mask choices affect both sides of the comparison, but that is a systematic/correctness issue rather than a circular reduction. The self-citations to earlier HAWC calibration and analysis papers provide instrument response, event reconstruction, and limit-setting conventions; they do not supply the physical expectation against which the data are judged. No equation or fitted parameter is renamed as a prediction, and the uniqueness of the result is not imported from same-author prior work.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The analysis depends on standard HAWC calibration, Planck dust gas tracers, and literature values for cloud masses and distances. The only hand-set parameters are the spectral index and the opacity cuts; no new particles or physical entities are introduced.

free parameters (2)
  • Spectral index alpha in the gamma-ray template = 2.75 (fixed by hand)
    Chosen as the assumed photon index for the power-law spectral model in Eq. 2.2. The upper limits depend on this choice, which is not fitted or varied.
  • Planck dust opacity cut thresholds for cloud templates = 5e-5 (Taurus, Aquila), 2.5e-5 (Hercules)
    Hand-picked thresholds in Sec. 2.1 used to select high-density regions. They define the spatial template and therefore the normalization and the upper limits.
assumptions (5)
  • domain assumption The molecular clouds are passive targets: gamma-ray emission is dominated by pp interactions of the ambient cosmic-ray sea, with negligible contributions from local accelerators or leptonic processes.
    Used throughout Sec. 3 to compute expected fluxes from Eqs. 3.1 and 3.2. If this fails, the comparison to hadronic models is invalid.
  • domain assumption Dust optical depth at 353 GHz traces molecular hydrogen column density with the reference conversion (tau_D/NH2)_ref = 1.18e-26 cm2.
    Invoked in Eq. 2.3 to build both the spatial templates and the cloud masses. Dust-to-gas ratio variations or line-of-sight confusion would bias the results.
  • ad hoc to paper The fixed spectral index alpha = 2.75 is adequate for computing quasi-differential upper limits over 1 to 100 TeV.
    Used in Eq. 2.2. No variation or test of this choice is shown, so it is specific to this analysis.
  • domain assumption Distances and masses from Schlafly (2014) and the cited catalogs are correct.
    Table 2 inputs enter the expected flux models in Eqs. 3.1 and 3.2. Incorrect distances or masses would change the expected level and the interpretation.
  • domain assumption HAWC detector response, background estimation, and fractional PMT binning from Refs. [8,10,11] are valid.
    The analysis relies on the HAWC collaboration's established reconstruction and likelihood methods without re-derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Studying Cosmic-ray Interactions in Giant Molecular Clouds with the HAWC Gamma-ray Observatory." pith.science (2026). https://pith.science/paper/WOMSDRV5

@misc{pith2026190806073,
  author       = {Pith},
  title        = {Pith review of: Studying Cosmic-ray Interactions in Giant Molecular Clouds with the HAWC Gamma-ray Observatory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WOMSDRV5}},
  note         = {Machine review of arXiv:1908.06073}
}
abstract

The cosmic-ray flux in the Galaxy can be characterized by combining the knowledge of the distribution of gas in the Galaxy and the observation of gamma rays. We analyze the data from the HAWC Observatory to look for gamma rays in three galactic giant molecular clouds, that are outside the galactic plane ($|b|>5^{\circ}$). We can then test the paradigm that the measured local cosmic-ray flux is the same as the "sea" of Galactic cosmic rays. Due to its large field of view, and high duty cycle, HAWC is suitable to search for gamma rays from large structures in the TeV gamma-ray regime. We present here preliminary results from measurements of the Aquila Rift, Hercules and Taurus molecular clouds.

Figures

Figures reproduced from arXiv: 1908.06073 by the authors.

Figure 1
Figure 1. GMCs obtained using data from the Planck survey [12]. The clouds are Taurus (top-left corner), Aquila (top-right corner) and Hercules (bottom). The figures were obtained after applying the method men￾tioned in Sec. 2.1. The coordinate system shown is galactic coordinates and the units are after normalizing the map. Since we will need the mass of the cloud later for calculating the expectation of gamma rays, this is … view at source ↗
Figure 2
Figure 2. shows the upper limits together with the expected model in green for Eq. 3.1 and orange for Eq. 3.2. The green band corresponds to the range from the 10th to 90th percentiles of the column density distribution, while the line corresponds to the median. The credible interval upper limits can be found in [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 12 canonical work pages

  1. [1]

    Aguilar, D

    M. Aguilar, D. Aisa, B. Alpat, et al., Phys. Rev. Lett. 114 (2015) 171103

    work page 2015

  2. [2]

    Strong, I

    A. Strong, I. V . Moskalenko, and V . S. Ptuskin,Ann. Rev. Sci. 57 (2007) 285–327

    work page 2007

  3. [3]

    Issa and A

    M. Issa and A. W. Wolfendale, Nature 292 (1981) 430–433

    work page 1981

  4. [4]

    Aharonian, Space Science Reviews 99 (Oct, 2001) 187–196

    F. Aharonian, Space Science Reviews 99 (Oct, 2001) 187–196

    work page 2001

  5. [5]

    Casanova, F

    S. Casanova, F. A. Aharonian, Y . Fukui, et al.,Publications of the Astronomical Society of Japan 62 (2010) 769–777

    work page 2010

  6. [6]

    R. Yang, E. de Ona Wilhelmi, and F. Aharonian, A&A 566 (2014)

    work page 2014

  7. [7]

    Aharonian, G

    F. Aharonian, G. Peron, R. Yang, et al., arXiv e-prints (2018) arXiv:1811.12118

    arXiv 2018

  8. [8]

    HAWC Collaboration et al., ApJ 843 (2017) 39

    work page 2017

Show all 14 references
  1. [9]

    Vianello, R

    G. Vianello, R. J. Lauer, P. W. Younk, et al., arXiv e-prints (Aug, 2015) arXiv:1508.07479

    2015 arXiv

  2. [10]

    HAWC Collaboration et al., ApJ 842 (2017) 9

    2017

  3. [11]

    HAWC Collaboration et al., ApJ 853 (Feb, 2018) 154

    2018

  4. [12]

    6 Giant Molecular Clouds with HA WC Hugo Alberto Ayala Solares1

    Planck Collaboration et al., A&A 536 (2011) 16. 6 Giant Molecular Clouds with HA WC Hugo Alberto Ayala Solares1

    2011

  5. [13]

    Kafexhiu, F

    E. Kafexhiu, F. Aharonian, A. M. Taylor, and G. S. Vila, PRD 90 (Dec., 2014) 123014

    2014

  6. [14]

    Schlafly, ApJ (2014). 7

    2014

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.