REVIEW 3 major objections 6 minor 1 cited by
Study of the IC 443 region with the HAWC observatory
T0 review · 3 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Using 2966 days of HAWC data, this paper reports that IC 443's gamma-ray spectrum extends to about 30 TeV without a cutoff, keeping the supernova remnant a viable PeVatron candidate, and identifies a new extended source near pulsar…
desk verdict Careful HAWC analysis of IC 443 with a solid extended-source detection, but the PeVatron interpretation rests on a point-source association the authors themselves concede is not secure. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument is carried by a maximum-likelihood source-search pipeline applied to HAWC data, which first adds point sources one by one when the test-statistic improvement exceeds 25, then tests each point source as a Gaussian extended source, and finally tests spectral curvature by comparing a simple power law with a log-parabola and an exponential cutoff. The spectral analyses use a power-law form $\Phi(E) = \Phi_0 (E/E_{\rm piv})^\alpha$ with pivot energy 2.3 TeV, and the hadronic interpretation is evaluated with a Markov Chain Monte Carlo spectral-fitting routine by comparing a broken power law and an exponential cutoff for the parent proton population. The decisive quantities are the test-statistic differences between models, the Bayesian Information Criterion between the broken power-law and cutoff proton spectra, and the flat posterior of the maximum proton energy, which together indicate that the data prefer a continuation of the spectrum over a cutoff.
What would settle it
A future observation with angular resolution good enough to localize the >10 TeV emission to the pulsar wind nebula rather than to the remnant shell would break the association with IC 443; alternatively, detecting a spectral cutoff below about 30 TeV in the point source would contradict the paper's central no-cutoff claim.
Extended reading notes
Core claim
On its own terms, the paper's central discovery is that the gamma-ray spectrum of the region around IC 443 continues as a featureless power law to roughly 30 TeV, with the likelihood analysis finding no significant curvature or cutoff. A maximum-likelihood fit to 2966 days of HAWC data yields a point source at $\alpha = 94.42^\circ$, $\delta = 22.35^\circ$ with spectral index $-3.14\pm0.18$, consistent with earlier TeV measurements, and the paper associates this point source with IC 443. In a combined spectral fit that includes prior GeV and TeV measurements plus the new HAWC data, a pion-decay model with a broken power-law proton spectrum is preferred over one with an exponential cutoff, and the posterior for the maximum proton energy stays flat out to 1 PeV, meaning the data do not require a cutoff. Alongside the point source the analysis finds an extended component, HAWC J0615+2213, with index $-2.49\pm0.08$; the extended emission is not well described by cosmic-ray illumination of known gas, and its centroid lies 0.29 degrees from the pulsar B0611+22, leading the authors to propose a TeV halo interpretation while deferring a full model to a later paper.
Load-bearing premise
The load-bearing assumption is that the point source belongs to the supernova remnant IC 443, even though its best-fit position sits 0.26 degrees from the remnant's catalog position and closer to the pulsar wind nebula; the paper itself concedes that the pulsar wind nebula interpretation should not be discarded.
Editorial extensions
If this is right
- If the result holds, IC 443 remains a viable hadronic PeVatron, with protons plausibly accelerated beyond the ~65 TeV limit of standard diffusive shock acceleration.
- The extended source HAWC J0615+2213, if a TeV halo, adds a ~90-kyr-old pulsar to the known halo population and implies a diffusion coefficient below the interstellar value by about a factor of 45.
- The non-detection of the neighbouring remnant G189.6+03.3, despite a hard GeV counterpart, means its flux must fall steeply above 1 TeV.
- Both sources are adequately described by simple power laws; the likelihood does not prefer a log-parabola or exponential-cutoff form for either.
Reading between the lines
- A testable extension is a joint spectral-spatial search with a more sensitive instrument above 10 TeV: if the no-cutoff tail belongs to the remnant, the point-source flux should continue to fall as a power law beyond 30 TeV.
- If the point source is actually dominated by the pulsar wind nebula, the hadronic interpretation is not destroyed, because inverse-Compton emission from the nebula is suppressed at the highest energies by the Klein-Nishina effect, so some hadronic contribution would still be required.
- The flat posterior on the maximum proton energy implies the current data cannot distinguish a cutoff at 65 TeV from one near 1 PeV; only observations that push the measured spectrum above ~50 TeV can make that discrimination.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a maximum-likelihood analysis of 2966 days of HAWC data toward the IC 443 region. The authors find a point source at (α = 94.42°, δ = 22.35°) that they associate with IC 443, with a power-law spectrum of index −3.14 ± 0.18 and no significant spectral curvature, and an extended source named HAWC J0615+2213 with index −2.49 ± 0.08 and Gaussian width 1.05°. They interpret the point source as gamma-ray emission from IC 443, argue that leptonic models cannot explain emission up to ~30 TeV, and conclude that IC 443 remains a viable hadronic PeVatron candidate. For the extended source, they discuss cosmic-ray illumination of interstellar gas, unresolved faint sources, and a possible TeV halo powered by pulsar B0611+22, deferring a detailed analysis to a future publication. The paper also reports upper limits on G189.6+03.3 and compares its model with LHAASO's catalog result for the region.
Significance. If the point-source association with IC 443 is correct, the work extends the known gamma-ray spectrum of this supernova remnant to ~30 TeV and supports hadronic acceleration above the ~65 TeV diffusive-shock-acceleration limit, which is directly relevant to Galactic PeVatron searches. The newly reported extended source HAWC J0615+2213 is also potentially valuable as a TeV-halo candidate around pulsar B0611+22. The analysis has notable strengths: it uses a transparent model-selection algorithm with explicit TS thresholds, compares two independent binning schemes, includes detector and modeling systematics, and shows residual significance maps consistent with background. However, the significance of the paper is tempered by two load-bearing weaknesses: the point-source identification with IC 443 is not quantitatively secured, and the claimed absence of a cutoff above ~30 TeV rests on sparse high-energy data. These issues mainly affect the astrophysical interpretation rather than the raw detections.
major comments (3)
- [§3.2 (Table 1), §4.1] The point-source association with IC 443 is not secured. The best-fit position is 0.26° from the Green-catalog position of IC 443 and 0.14° from the pulsar wind nebula CXOU J061705.3+222212, while the combined statistical and systematic position uncertainty in Table 1 is of order 0.1°. The paper itself states in §4.1 that 'emission from the pulsar wind nebula should not be discarded' and that HAWC cannot distinguish where the emission originates. Because the hadronic-PeVatron conclusion in §4.1 and §5 depends on this identification, the manuscript should present the association as a working hypothesis and support it with a quantitative test, for example a likelihood comparison of models with the point source placed at the IC 443 position, at the PWN position, or left free, together with a chance-coincidence estimate given the source density in the region.
- [§4.1] The leptonic-exclusion argument is environment-specific and does not apply if the point-source emission originates in the pulsar wind nebula. The electron maximum-energy bounds of roughly 3–10 TeV cited in §4.1 are derived for the SNR shock environment and do not constrain electrons accelerated in a young pulsar wind nebula. A PWN origin would be leptonic and could plausibly produce a hard, power-law-like gamma-ray spectrum; the Klein-Nishina remark in §4.1 does not by itself rule out emission at the highest observed energies. Consequently, the sentence concluding that 'hadronic processes are a more plausible explanation' is not established for the point-source emission as long as the PWN alternative remains viable.
- [Abstract; §3.3; Appendix B] The claim that the spectrum extends to ~30 TeV with no cutoff is stronger than the data support. In Table 3, the highest SED point with a statistically significant TS is at 15.4 TeV (TS = 11.1); the 38.2 TeV bin is an upper limit, and the 1σ energy range [0.3, 30] TeV is a derived detection interval rather than a direct detection of 30 TeV photons. The exponential-cutoff model is not preferred (Table 9), but this test has little sensitivity to a cutoff near or above 30 TeV given the sparse high-energy points. The conclusion should be framed as 'no cutoff is required by the data' rather than as a positive statement that the spectrum extends to ~30 TeV without a cutoff.
minor comments (6)
- [Footnote 48] The pulsar discovery is cited as 'Davies, J. and Lyne, A. and Seiradakis, J. (1972)'; please use standard author-year formatting, and correct the typo 'used to measured the distance'.
- [Figure 1 caption] The X-ray compact object is written as 'CXOU J061705.3+222122' in the caption but as 'CXOU J061705.3+222212' in the text and Table 1; the spelling should be unified.
- [§4.2.3, Eq. (9)] Equation (9) gives an upper limit on the diffusion coefficient, but the following sentence says 'We find that the diffusion coefficient to be 2.21 × 10^28 cm^2 s^-1'; please use 'upper limit' consistently in the text.
- [§4.2.3] The sentence 'All estimates are below ∼0.1 cm−3, the energy density of the interstellar medium' should have units of eV cm^-3 for the energy density.
- [§4.2.1] The sentence 'This would mean that there would still be some emission left that is not described by these faint sources' is unclear; please rephrase to indicate what fraction of the extended-source flux remains unexplained after accounting for unresolved sources.
- [Table 2] Adding the systematic uncertainties for the energy-scheme parameters would make the cross-check more informative; as printed, the table shows only statistical uncertainties.
Circularity Check
No significant circularity: the spectral measurements, cutoff search, and halo diffusion estimate are direct fits to HAWC data with external cross-checks.
full rationale
This is an observational measurement paper rather than a derivation that reduces to its inputs. The two central results—the point-source power-law index (-3.14 ± 0.18, no significant cutoff up to ~30 TeV) and the extended source HAWC J0615+2213 with index (-2.49 ± 0.08)—are obtained from maximum-likelihood fits to HAWC data using explicit model-comparison thresholds (TS > 25, TS > 16), and the SED points in Tables 3 and 4 are produced by fixing all model parameters except the normalization in each energy bin. No fitted parameter is renamed as a prediction of the same quantity. The Naima hadronic modeling takes the measured SED as input, uses externally motivated parameters (e.g., n = 20 cm^-3 from Ackermann et al. 2013, distances from parallax), and the flat posterior on the maximum proton energy is a data-driven result, not an assumed input. The TeV-halo diffusion coefficient D <= 2.21e28 cm2/s follows from the measured Gaussian width (theta = 1.05 deg), the externally measured pulsar distance (d0 = 3.55 kpc), and the standard cooling-time formula (Eq. 7); it is not fitted within the same model. Self-citations (A. Albert et al. 2023; A. U. Abeysekara et al. 2022; A. Albert et al. 2024) are methodological references for the analysis framework, energy reconstruction, and source-searching algorithm, and none is used to forbid alternative models or to supply a uniqueness theorem. The HAWC results are also cross-checked against external LHAASO catalog measurements in Appendix C, so the paper is self-contained against an independent benchmark. The association of the point source with IC 443 rather than the nearby pulsar wind nebula is an astrophysical interpretation, and the paper explicitly concedes that 'emission from the pulsar wind nebula should not be discarded'; this is a source-association caveat, not a circular derivation. No circular step could be identified with quoted text and a specific reduction.
Assumptions & free parameters
free parameters (6)
- Point source normalization Φ0 =
5.9e-14 TeV^-1 cm^-2 s^-1 at 2.3 TeV
- Point source spectral index α =
-3.14 ± 0.18
- Extended source normalization Φ0 =
3.18e-13 TeV^-1 cm^-2 s^-1 at 2.3 TeV
- Extended source spectral index α =
-2.49 ± 0.08
- Extended source Gaussian width σ =
1.05 deg
- GDE scale factor =
2.62 ± 1.20
assumptions (5)
- domain assumption The HAWC gamma/hadron selection and detector response are correct as described in previous HAWC papers.
- domain assumption The HERMES galactic diffuse emission template accurately represents the gamma-ray background in the IC 443 region.
- domain assumption The point source found at (94.42, 22.35) is the same object as IC 443 despite a 0.26 degree offset from the catalog position.
- domain assumption The extended source lies at the distance of the pulsar B0611+22, 3.55 kpc, for the halo calculations.
- domain assumption The interstellar medium diffusion coefficient is 1e30 cm^2/s and the electron cooling time is 12 kyr at 100 TeV.
Cite this review
Pith. "Pith review of Study of the IC 443 region with the HAWC observatory." pith.science (2026). https://pith.science/paper/GSF43FJJ
@misc{pith2026250112613,
author = {Pith},
title = {Pith review of: Study of the IC 443 region with the HAWC observatory},
year = {2026},
howpublished = {\url{https://pith.science/paper/GSF43FJJ}},
note = {Machine review of arXiv:2501.12613}
}
abstract
Supernova remnants are one potential source class considered a PeVatron (i.e. capable of accelerating cosmic rays above PeV energies). The shock fronts produced after the explosion of the supernova are ideal regions for particle acceleration. IC 443 is a supernova remnant that has been studied extensively at different wavelengths. Using 2966 days of gamma-ray data from the HAWC observatory, we study the emission of IC 443 with the objective of finding signatures of cosmic-ray acceleration at the PeV scale. Using a maximum likelihood method, we find a point source located at ($\alpha$=94.42$^{\circ}$, $\delta$=22.35$^{\circ}$) that we associate with IC 443. The measured spectrum is a simple power law with an index of $-3.14\pm$0.18, which is consistent with previous TeV observations. Although we cannot confirm that IC 443 is a hadronic PeVatron, we do not find any sign that the spectrum has a cut off at tens of TeV energies, with the spectrum extending to $\sim$30 TeV. Furthermore, we also find a new extended component in the region whose emission is described by a simple power law with an index of $-2.49\pm$0.08 and which we call HAWC J0615+2213. While we show evidence that this new source might be a new TeV halo, we defer a detailed analysis of this new source to another publication.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 1 Pith paper
-
Gamma rays as a signature of r-process producing supernovae: remnants and future Galactic explosions
Using a magnetorotational supernova model, the paper predicts that gamma-ray lines from tin-126 and iron-60 in Galactic remnants, and from tellurium, iodine, and antimony isotopes in a future nearby supernova, would b...
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
-
[3]
A., Ackermann , M., Ajello , M., et al
Abdo , A. A., Ackermann , M., Ajello , M., et al. 2009, ApJL, 706, L1
work page 2009
-
[4]
A., Ackermann , M., Ajello , M., et al
Abdo , A. A., Ackermann , M., Ajello , M., et al. 2010, ApJ, 712, 459
work page 2010
-
[5]
U., Albert , A., Alfaro , R., et al
Abeysekara , A. U., Albert , A., Alfaro , R., et al. 2019, ApJ, 881, 134
work page 2019
-
[6]
U., Albert , A., Alfaro , R., et al
Abeysekara , A. U., Albert , A., Alfaro , R., et al. 2017, ApJ, 843, 39
work page 2017
-
[7]
U., Albert , A., Alfaro , R., et al
Abeysekara , A. U., Albert , A., Alfaro , R., et al. 2022, in 37th International Cosmic Ray Conference, 828
work page 2022
-
[8]
U., Albert, A., Alfaro, R., et al
Abeysekara, A. U., Albert, A., Alfaro, R., et al. 2023, Nuclear Instruments and Methods in Physics Research Section A, 1052, 168253
work page 2023
Show all 64 references
-
[9]
A., et al
Abramowski , A., Aharonian , F., Benkhali , F. A., et al. 2016, , 531, 476
2016
-
[10]
A., Aliu , E., Arlen , T., et al
Acciari , V. A., Aliu , E., Arlen , T., et al. 2009, ApJL, 698, L133
2009
-
[11]
2013, Science, 339, 807
Ackermann, M., Ajello, M., Allafort, A., et al. 2013, Science, 339, 807
2013
-
[12]
2024, ApJ, 972, 144
Albert , A., Alfaro , R., Alvarez , C., et al. 2024, ApJ, 972, 144
2024
-
[13]
2021, , 914, 106
Albert , A., Alfaro , R., Alvarez , C., et al. 2021, , 914, 106
2021
-
[14]
2020, ApJ, 905, 76
Albert , A., Alfaro , R., Alvarez , C., & HAWC Collaboration . 2020, ApJ, 905, 76
2020
-
[15]
2023, , 954, 205
Albert , A., Alvarez , C., Avila Rojas , D., et al. 2023, , 954, 205
2023
-
[16]
2007, ApJL, 664, L87
Albert , J., Aliu , E., Anderhub , H., et al. 2007, ApJL, 664, L87
2007
-
[17]
2017, , 472, 51
Ambrocio-Cruz , P., Rosado , M., de la Fuente , E., Silva , R., & Blanco-Pi \ n on , A. 2017, , 472, 51
2017
-
[18]
M., Aharonian , F
Atoyan , A. M., Aharonian , F. A., & V \"o lk , H. J. 1995, , 52, 3265
1995
-
[19]
G., Ellison , D
Baring , M. G., Ellison , D. C., Reynolds , S. P., Grenier , I. A., & Goret , P. 1999, ApJ, 513, 311
1999
-
[20]
Bell, A. R. 2004, MNRAS, 353, 550
2004
-
[21]
R., Schure, K
Bell, A. R., Schure, K. M., Reville, B., & Giacinti, G. 2013, Monthly Notices of the Royal Astronomical Society, 431, 415
2013
-
[22]
& Becker , W
Camilloni , F. & Becker , W. 2023, AAP, 680, A83
2023
-
[23]
2024, ApJS, 271, 25
Cao , Z., Aharonian , F., An , Q., Axikegu , & (The Lhaaso Collaboration) . 2024, ApJS, 271, 25
2024
-
[24]
2009, MNRAS, 395, 895
Caprioli, D., Blasi, P., Amato, E., & Vietri, M. 2009, MNRAS, 395, 895
2009
-
[25]
Castelletti , G., Dubner , G., Clarke , T., & Kassim , N. E. 2011, AAP, 534, A21
2011
-
[26]
Chevalier , R. A. 1999, ApJ, 511, 798
1999
-
[27]
2020, Astroparticle Physics, 123, 102492
Cristofari, P., Blasi, P., & Amato, E. 2020, Astroparticle Physics, 123, 102492
2020
-
[28]
M., Ungerechts , H., Cohen , R
Dame , T. M., Ungerechts , H., Cohen , R. S., et al. 1987, ApJ, 322, 706
1987
-
[29]
and Lyne, A
Davies, J. and Lyne, A. and Seiradakis, J. 1972, Nature, 240, 229
1972
-
[30]
T., Goss , W
Deller , A. T., Goss , W. M., Brisken , W. F., et al. 2019, , 875, 100
2019
-
[31]
Drury, L. O. 1983, Reports on Progress in Physics, 46, 973
1983
-
[32]
2021, AAP, 653, A18
Dundovic , A., Evoli , C., Gaggero , D., & Grasso , D. 2021, AAP, 653, A18
2021
-
[33]
Erickson , W. C. & Mahoney , M. J. 1985, ApJ, 290, 596
1985
-
[34]
2017, JCAP, 2017, 015
Evoli , C., Gaggero , D., Vittino , A., et al. 2017, JCAP, 2017, 015
2017
-
[35]
Fesen , R. A. 1984, ApJ, 281, 658
1984
-
[36]
M., Chatterjee , S., Slane , P
Gaensler , B. M., Chatterjee , S., Slane , P. O., et al. 2006, ApJ, 648, 1037
2006
-
[37]
Giacinti , G., Mitchell , A. M. W., L \'o pez-Coto , R., et al. 2020, , 636, A113
2020
-
[38]
Green , D. A. 2019, Journal of Astrophysics and Astronomy, 40, 36
2019
-
[39]
2014, PRL, 113, 155005
Guo, F., Li, H., Daughton, W., & Liu, Y.-H. 2014, PRL, 113, 155005
2014
-
[40]
2016, AAP, 594, A116
HI4PI Collaboration , Ben Bekhti , N., Fl \"o er , L., et al. 2016, AAP, 594, A116
2016
-
[41]
Lagage, P. O. & Cesarsky, C. J. 1983, Astronomy and Astrophysics, 118, 223
1983
-
[42]
S., et al
Lee , J.-J., Koo , B.-C., Yun , M. S., et al. 2008, , 135, 796
2008
-
[43]
S., et al
Lee, J.-J., Koo, B.-C., Yun, M. S., et al. 2008, The Astronomical Journal, 135, 796
2008
-
[44]
2022, ApJ, 927, 226
Li , J., Jiang , B., & Zhao , H. 2022, ApJ, 927, 226
2022
-
[45]
2017, Phys
Linden, T., Auchettl, K., Bramante, J., et al. 2017, Phys. Rev. D, 96, 103016
2017
-
[46]
M., Clearfield , C
Olbert , C. M., Clearfield , C. R., Williams , N. E., Keohane , J. W., & Frail , D. A. 2001, , 554, L205
2001
-
[47]
2012, Space Science Reviews, 173, 535
Petrosian, V. 2012, Space Science Reviews, 173, 535
2012
-
[48]
2011, A&A, 536, 16
Planck Collaboration et al. 2011, A&A, 536, 16
2011
-
[49]
R., et al
Rajwade , K., Seymour , A., Lorimer , D. R., et al. 2016, , 462, 2518
2016
-
[50]
M., et al
Seta , M., Hasegawa , T., Dame , T. M., et al. 1998, , 505, 286
1998
-
[51]
2023, JCAP, 2023, 027
Sharma , P., Ou , Z., Henry-Cadrot , C., Dubos , C., & Suomij \"a rvi , T. 2023, JCAP, 2023, 027
2023
-
[52]
2008, , 485, 777
Troja , E., Bocchino , F., Miceli , M., & Reale , F. 2008, , 485, 777
2008
-
[53]
2008, The Astrophysical Journal, 689, L125
Troja, E., Bocchino, F., Miceli, M., Reale, F., & Dubner, G. 2008, The Astrophysical Journal, 689, L125
2008
-
[54]
2006, , 649, 258
Troja , E., Bocchino , F., & Reale , F. 2006, , 649, 258
2006
-
[55]
2021, AAP, 649, A14
Ustamujic , S., Orlando , S., Greco , E., et al. 2021, AAP, 649, A14
2021
-
[56]
J., Younk , P., et al
Vianello , G., Lauer , R. J., Younk , P., et al. 2015, arXiv e-prints, arXiv:1507.08343
2015 arXiv
-
[57]
2022, Frascati Phys.Ser., 74, 153
Vink , J. 2022, Frascati Phys.Ser., 74, 153
2022
-
[58]
V., Balser , D
Wenger , T. V., Balser , D. S., Anderson , L. D., & Bania , T. M. 2018, , 856, 52
2018
-
[59]
1893, AN, [3130] 131, 157
Wolf , M. 1893, AN, [3130] 131, 157
-
[60]
2009, The Astrophysical Journal, 705, L6
Yamaguchi, H., Ozawa, M., Koyama, K., et al. 2009, The Astrophysical Journal, 705, L6
2009
-
[61]
M., Manchester , R
Yao , J. M., Manchester , R. N., & Wang , N. 2017, , 835, 29
2017
-
[62]
2015, in International Cosmic Ray Conference, Vol
Zabalza , V. 2015, in International Cosmic Ray Conference, Vol. 34, 34th International Cosmic Ray Conference (ICRC2015), 922
2015
-
[63]
& Yan, H
Zhang, B. & Yan, H. 2011, ApJ, 726, 90
2011
-
[64]
2018, ApJ, 859, 141
Zhang , S., Tang , X., Zhang , X., et al. 2018, ApJ, 859, 141
2018
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.