{"id":"64721d10-517b-4966-896d-507e038a50d0","arxiv_id":"2501.12613","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"HAWC detects IC 443 up to ~30 TeV without a spectral cutoff and finds a new extended gamma-ray source, HAWC J0615+2213, that may be a TeV halo.","lead":"The HAWC observatory measured gamma rays from the supernova remnant IC 443 and reports a point source with a spectrum that shows no cutoff up to about 30 TeV, along with a new extended source named HAWC J0615+2213. This is relevant to the question of whether supernova remnants can accelerate cosmic rays to PeV energies, and whether a middle-aged pulsar is powering a new TeV halo.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The point-source association with IC 443 is not secured: the best-fit HAWC position is 0.14 deg from the pulsar wind nebula CXOU J061705.3+222212, and the paper itself concedes PWN emission cannot be discarded.","rationale":"The reader's weakest-assumption identifies exactly the assumption I find most load-bearing: the identification of the HAWC point source with IC 443. The paper's central astrophysical claim, that IC 443 remains a viable hadronic PeVatron candidate, requires this association. The data, as presented in Table 1 and Fig. 1, place the point source closer to the pulsar wind nebula than to the SNR catalog position, and Sec. 4.1 explicitly warns that PWN emission cannot be discarded. A PWN origin would make the source leptonic, so the hadronic interpretation would no longer follow. This is not an internal inconsistency: the likelihood analysis and spectral measurements are plausible, and the raw detection of a point source is not in question. Rather, the physical interpretation is underdetermined by the angular resolution, so the paper's own caveat should be elevated to a central limitation. The reader's conditional verdict already reflects this, so I do not recommend a change in verdict, but the concrete test above would either strengthen the association or require the paper to present the source as positionally ambiguous. A secondary overreach, the use of the flat posterior on the proton cutoff to suggest PeVatron capability, is real but less fundamental than the association problem.","tokens_in":20695,"tokens_out":7617,"duration_ms":83614,"concrete_test":"Perform a binned likelihood fit of the HAWC data (fhit and energy schemes) with the GDE and extended source fixed to the best-fit model, then compare three point-source models: (a) source fixed at the Green-catalog IC 443 position, (b) source fixed at CXOU J061705.3+222212, and (c) source position free. Compute the TS between models (a)/(b) and (c), and also compute TS for adding a second point source at the alternative position. If the PWN-fixed model is within ΔTS < 9 of the free-position model, or if adding a PWN component gives ΔTS > 16, the paper should describe the source as positionally ambiguous and retract or qualify the hadronic PeVatron interpretation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's hadronic PeVatron conclusion depends on identifying the point source with the SNR IC 443. However, Table 1 places the point source at (94.42, 22.35), 0.26 deg from the Green-catalog IC 443 position and 0.14 deg from the pulsar wind nebula CXOU J061705.3+222212. The paper states in Sec 4.1 that 'emission from the pulsar wind nebula should not be discarded' and that HAWC cannot distinguish where the emission originates. Because a PWN origin would be leptonic, the high-energy (up to ~30 TeV) spectrum would not be a hadronic SNR signature, and the claim that IC 443 remains a viable hadronic PeVatron candidate would lose its basis. The 'continuation' of the Fermi/VERITAS/MAGIC SED is suggestive but not decisive, since the HAWC observation is a single power law over a broad region and the morphology is unresolved. This is the load-bearing assumption: if the source is CXOU J061705.3+222212 rather than the SNR shell, the main astrophysical interpretation changes, while the raw detection remains valid.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":20867,"tokens_out":8411,"duration_ms":79467,"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":[{"comment":"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.","section":"§3.2 (Table 1), §4.1"},{"comment":"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.","section":"§4.1"},{"comment":"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.","section":"Abstract; §3.3; Appendix B"}],"minor_comments":[{"comment":"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'.","section":"Footnote 48"},{"comment":"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.","section":"Figure 1 caption"},{"comment":"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.","section":"§4.2.3, Eq. (9)"},{"comment":"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.","section":"§4.2.3"},{"comment":"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.","section":"§4.2.1"},{"comment":"Adding the systematic uncertainties for the energy-scheme parameters would make the cross-check more informative; as printed, the table shows only statistical uncertainties.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent observational study with robust detections of two gamma-ray sources, transparent model selection, and careful treatment of systematics. My main concern is that the abstract and conclusions overstate the association with IC 443 and the absence of a spectral cutoff; these issues are fixable through reframing and additional likelihood or positional tests. I would not reject the manuscript, because the measurements themselves are valuable and the authors are transparent about the PWN caveat in §4.1. The editor may wish to ensure that the revision addresses the association quantitatively and softens the PeVatron wording in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a careful, transparent HAWC analysis of the IC 443 region, and the new extended source HAWC J0615+2213 (TS 88.8, hard index -2.49) is a solid detection. But the point source the authors associate with IC 443 sits 0.26 deg from the SNR center and 0.14 deg from the pulsar wind nebula CXOU J061705.3+222212, and the paper itself concedes that PWN emission cannot be discarded. That ambiguity is load-bearing for the PeVatron story.\n\nWhat's genuinely new: the HAWC point-source spectrum as a simple power law to ~30 TeV with no sign of a cutoff, and the extended source with its tentative halo interpretation. The analysis is well executed—maximum likelihood with HAL/threeML, explicit model-selection thresholds, systematic uncertainties from detector responses and GDE modeling, and a comparison between binning schemes. The SED points and the G189.6 upper limits are useful. The authors are also honest about what they can't do: they explicitly defer a full halo analysis and note that HAWC's angular resolution cannot separate where the point-source emission originates.\n\nThe soft spots are in the interpretation, not the measurement. First, the association of the point source with IC 443 rather than the PWN is not established. The 'continuation' of the Fermi/VERITAS/MAGIC SED is suggestive but not decisive, since the HAWC source is a single power law over a broad region and the morphology is unresolved. The argument that leptonic emission can't reach 30 TeV is based on the SNR environment (dense clouds, ~10 µG field); a PWN has a different environment, and the Klein-Nishina argument the paper uses is itself sensitive to the target photon field. If the point source is the PWN, the hadronic PeVatron interpretation loses its basis. Second, the flat posterior on the proton cutoff energy is not evidence for a high cutoff—it just means the data don't constrain it. The phrase 'strong candidate for a PeVatron' overstates what a flat posterior shows. Third, LHAASO has already reported an extended source in this region; the HAWC result is a new measurement and a two-component decomposition, not a wholly new discovery, though the authors handle this by testing the LHAASO model against their data.\n\nNone of this sinks the paper as an observational contribution. The extended source detection is statistically strong and the systematics work is honest. The issue is that the headline interpretation should be framed as conditional: if the point source is the remnant, it remains a PeVatron candidate; if it's the PWN, the case is much weaker. The authors already say most of this in Section 4.1; they just don't carry it into the abstract/conclusions with the same weight.\n\nI'd send this to a competent referee. The analysis deserves review and the community will want the extended source measurement in the literature. For my own work, I'd cite the extended source detection, not the PeVatron claim.","headline":"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.","tokens_in":22043,"tokens_out":3909,"would_cite":true,"duration_ms":39414,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["IC 443","supernova remnants","PeVatron","gamma-ray astronomy","TeV halo","pulsar wind nebula","HAWC observatory","cosmic-ray acceleration"],"falsifier":"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.","tokens_in":20433,"feed_emoji":"🌌","tokens_out":9430,"duration_ms":85604,"temperature":0.7,"pith_summary":"Using 2966 days of data from the HAWC observatory, this paper asks whether the supernova remnant IC 443 accelerates cosmic rays to PeV energies. It reports a point source that it associates with IC 443, with a simple power-law spectrum of index $-3.14\\pm0.18$ and no evidence for a cutoff up to about 30 TeV. Because electrons in this environment cannot easily reach the energies needed to produce such gamma rays, the authors argue that the tail is best explained by hadronic cosmic rays and that IC 443 remains a viable PeVatron candidate. The same analysis uncovers a new extended source, HAWC J0615+2213, with a harder spectrum (index $-2.49\\pm0.08$), which the paper proposes may be a TeV halo powered by the middle-aged pulsar B0611+22. The point source sits 0.26 degrees from the catalog position of IC 443, and the paper itself notes that emission from the pulsar wind nebula should not be discarded.","feed_headline":"IC 443's gamma-ray spectrum shows no cutoff up to 30 TeV","feed_subtitle":"Featureless power law keeps hadronic cosmic-ray acceleration on the table; a new extended source may be a TeV halo.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Establishes the hadronic pion-decay origin of IC 443's gamma-ray emission at GeV energies and supplies the gas-density assumption used in the spectral modeling.","marker":"M. Ackermann et al. 2013"},{"why":"Provides the GeV spectral data included in the combined pion-decay fit.","marker":"A. A. Abdo et al. 2010"},{"why":"Provides TeV data that the HAWC point-source spectrum is compared against.","marker":"V. A. Acciari et al. 2009"},{"why":"Provides earlier TeV data used in the multi-wavelength spectral fit.","marker":"J. Albert et al. 2007"},{"why":"Describes the HAWC detector performance, energy reconstruction, and binning schemes used for the analysis.","marker":"A. Albert et al. 2024"},{"why":"Supplies the background-estimation method on which the maximum-likelihood analysis relies.","marker":"A. U. Abeysekara et al. 2019"},{"why":"Presents the LHAASO catalog source in the IC 443 region that the paper compares with its own model.","marker":"Z. Cao et al. 2024"},{"why":"Provides the spectral-fitting code used to compare pion-decay and inverse-Compton models.","marker":"V. Zabalza 2015"},{"why":"Gives the diffusion-coefficient estimate used to test whether the extended source is a TeV halo.","marker":"A. M. Atoyan et al. 1995"}],"fun_headline_variants":["HAWC finds no spectral cutoff for IC 443 up to 30 TeV","New HAWC source HAWC J0615+2213 may be TeV halo","IC 443's gamma rays: no cutoff, new neighbor source","PeVatron candidate IC 443 shows featureless spectrum","HAWC data: IC 443's spectrum extends to 30 TeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["HAWC finds no spectral cutoff for IC 443 up to 30 TeV","New HAWC source HAWC J0615+2213 may be TeV halo","IC 443's gamma rays: no cutoff, new neighbor source","PeVatron candidate IC 443 shows featureless spectrum","HAWC data: IC 443's spectrum extends to 30 TeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1753,"prompt_tokens":1096,"completion_tokens":657,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":712,"completion_tokens_details":{"reasoning_tokens":558}},"tokens_in":712,"tokens_out":657,"duration_ms":6644,"temperature":1.0,"reasoning_tokens":558,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:59:13.778780+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"A., Ackermann , M., Ajello , M., et al","cited_arxiv_id":null,"evidence_quote":"Provides the GeV spectral data included in the combined pion-decay fit."},{"cited_title":"A., Aliu , E., Arlen , T., et al","cited_arxiv_id":null,"evidence_quote":"Provides TeV data that the HAWC point-source spectrum is compared against."},{"cited_title":"2007, ApJL, 664, L87","cited_arxiv_id":null,"evidence_quote":"Provides earlier TeV data used in the multi-wavelength spectral fit."},{"cited_title":"2024, ApJ, 972, 144","cited_arxiv_id":null,"evidence_quote":"Describes the HAWC detector performance, energy reconstruction, and binning schemes used for the analysis."},{"cited_title":"U., Albert , A., Alfaro , R., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the background-estimation method on which the maximum-likelihood analysis relies."},{"cited_title":"2024, ApJS, 271, 25","cited_arxiv_id":null,"evidence_quote":"Presents the LHAASO catalog source in the IC 443 region that the paper compares with its own model."},{"cited_title":"M., Aharonian , F","cited_arxiv_id":null,"evidence_quote":"Gives the diffusion-coefficient estimate used to test whether the extended source is a TeV halo."}],"review_version":1}