{"id":"1e4ce452-a093-4e9d-b664-eeb766b09ec4","arxiv_id":"1908.05119","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"IC 443 shows neutral iron K-alpha line enhancements in two molecular-cloud interaction regions, consistent with MeV protons leaking out of the supernova remnant.","lead":"We found bright clumps of a neutral iron X-ray line in two regions of the supernova remnant IC 443, right where the remnant is smashing into dense molecular clouds. The likely source is low-energy (MeV) protons accelerated by the remnant, which would fill a missing piece in cosmic-ray acceleration physics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Admitted factor-40 mismatch between the diffusion coefficient needed for Reg 2's size and the fitted escape model makes the claimed quantitative consistency with the CR escape model unsupported.","rationale":"The reader's weakest assumption is the same as the most load-bearing concern I identify: the Appendix itself documents that reproducing the observed angular size of Reg 2 requires a diffusion coefficient about 40 times larger than the fitted model value, and the paper responds with untested escape clauses rather than a quantitative treatment. This matters because the central argument for co-acceleration rests on the claimed consistency between the observed Fe I K-alpha intensity and the CR escape model prediction; if the transport parameters are free, that consistency is not a precise test. I do not see a reason to move the verdict, because the concern is addressable and the observational result may still stand, so the reader's CONDITIONAL verdict remains appropriate. Other issues, such as the a posteriori selection of regions from the same line map and the reliance on unpublished CO data, are real but secondary to the quantitative self-consistency of the proposed model.","tokens_in":10168,"tokens_out":9405,"duration_ms":103923,"concrete_test":"Recompute the Fe I K-alpha intensity prediction of the Appendix with a self-consistent treatment of the two escape clauses: (i) inject CRs over the full shock-MC contact region rather than a point source, and (ii) set the in-cloud diffusion coefficient to D' = 4e26 cm^2/s while keeping the fitted gamma-ray parameters fixed. If the predicted I6.4 keV integrated over the Reg 1 and Reg 2 extraction areas changes by more than the quoted +/- 0.1 photons s^-1 cm^-2 sr^-1, or if the gamma-ray fit must be redone, then the paper's consistency claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the Fe I K-alpha intensities in Regs 1 and 2 are explained by MeV protons leaking from the IC 443 shock into adjacent molecular clouds depends on the quantitative comparison with the CR escape model in Section 4.2 and the Appendix. That comparison has a self-consistency problem that the paper itself admits. From Equation A2 with the fitted chi = 0.011 and delta = 0.58, the diffusion coefficient at 10 MeV is D ~ 9e24 cm^2/s. Yet the 4-arcmin (about 2 pc) extent of Reg 2, combined with the 10-MeV ionization cooling time t_cool ~ 500 yr at nH = 730 cm^-3, requires D' >= R^2/(6 t_cool) ~ 4e26 cm^2/s, roughly 40 times larger. The Appendix responds that a different coefficient inside/outside the MC or an extended injection region would be allowed, but neither possibility is modeled or constrained. The observed-to-predicted intensity agreement (0.273 +/- 0.119 and 0.344 +/- 0.277 vs 0.29 +/- 0.10) is therefore not a test of the escape model as fitted; it is an order-of-magnitude match with two adjustable transport parameters. This weakens the co-acceleration conclusion, though it does not refute it, because an extended injection geometry could in principle resolve the discrepancy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports Suzaku XIS observations of IC 443 and constructs an Fe I K-alpha line intensity map. It identifies two blob-like enhancements, Reg 1 in the northwest and Reg 2 in the middle of the remnant, both associated with molecular clouds. Spectral fits yield line intensities of 0.273 +/- 0.119 and 0.344 +/- 0.277 photons s^-1 cm^-2 sr^-1 for Regs 1 and 2, with large equivalent widths. The authors argue that the line is produced by inner-shell ionization of neutral Fe by MeV protons leaking from the SNR shock into adjacent molecular clouds, and they compare the observed intensity with a prediction from a CR escape model fitted to gamma-ray data. They find consistency (predicted 0.29(+0.10/-0.11)) and conclude that MeV protons are accelerated together with GeV-TeV protons in the SNR.","tokens_in":10413,"tokens_out":3937,"duration_ms":40584,"significance":"If the detection and the quantitative comparison were robust, this paper would provide a rare direct probe of the long-missing MeV cosmic-ray population in an SNR, with implications for the co-acceleration of low- and high-energy protons. The analysis uses an observable (Fe I K-alpha line) that is independent of the gamma-ray data used to fit the CR model, which is a genuine strength, and the model is described in enough detail to be checked. The main caveats are the marginal significance of one of the two regions, the lack of a trials correction for source selection, and an admitted factor-of-40 discrepancy in the diffusion coefficient needed to explain the spatial extent of Reg 2. These issues do not refute the scenario but they currently prevent the paper from making a strong quantitative case for the escape model as fitted.","major_comments":[{"comment":"","section":"Appendix, Eq. (A2) and following paragraph; Section 4.2"},{"comment":"","section":"Sections 3.1 and 3.2, Table 2"},{"comment":"","section":"Section 4.2 and Appendix"}],"minor_comments":[{"comment":"","section":"Section 4.1"},{"comment":"","section":"Section 4.2"},{"comment":"","section":"Section 3.2"},{"comment":"","section":"Figure 1 and Appendix"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and addresses an important question, but the advertised quantitative consistency with the CR escape model is undermined by the factor-of-40 diffusion discrepancy acknowledged in the Appendix. I would encourage the editor to request a revision that either models the proposed remedies or tempers the central claim accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is the map: Nobukawa et al. have imaged the 6.4 keV Fe I K-alpha line across IC 443 with Suzaku and found two new bright blobs, one in the northwest (Reg 1) and one in the middle (Reg 2), both sitting on molecular clouds. That is a step beyond Hirayama et al. (2019), who only had the northeast. The equivalent widths (>1.2 keV and 0.7 keV) are a good argument for proton bombardment rather than electrons or X-ray irradiation, and the PWN irradiation estimate seems to kill that alternative cleanly. Credit where due: this is a solid observational addition to the low-energy cosmic-ray story.\n\nThe model comparison is also a legitimate cross-check. They fit the Makino et al. escape model to the gamma-ray spectrum, then predict the Fe line intensity without fitting the line data. The predicted 0.29 photons s^-1 cm^-2 sr^-1 sits comfortably inside the observed values for Reg 1 and Reg 2. That is not circular.\n\nThe soft spots are real, and the paper itself points to the biggest one. Reg 2 is only a 2.1 sigma detection, and both regions were selected from the same line map used to quote significances, so the numbers overstate confidence. The CO data are private communication, so the density and geometry are hard to check. And in the appendix they admit that the diffusion coefficient needed to explain the 4-arcmin size of Reg 2 is about forty times larger than the value in their fitted model. They wave at a different diffusion coefficient inside the cloud or an extended injection region, but they do not model either. So the abstract's \"consistent with the prediction\" is too strong. What they have is an order-of-magnitude match with a transport knob that can absorb the discrepancy.\n\nNone of this kills the paper. The line is probably real in Reg 1, the association with clouds is compelling, and the co-acceleration scenario is the natural reading. But the quantitative support is suggestive, not demonstrative. The right fix is a revision that softens the conclusion and either models the extended injection or clearly frames it as an unconstrained possibility.\n\nThis is a paper for SNR and cosmic-ray specialists, especially the LECR crowd. It deserves a serious referee, and a good referee could help the authors turn a plausible result into a convincing one. Send it out.","headline":"A genuinely new Fe K-alpha map of IC 443 with a plausible LECR origin, but the quantitative agreement with the escape model is weaker than the abstract claims.","tokens_in":11070,"tokens_out":2073,"would_cite":true,"duration_ms":21824,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two blob-like 6.4 keV Fe I K-alpha enhancements in IC 443 are attributed to MeV protons accelerated at the SNR shock and leaking into molecular clouds, with intensities matching a model that also produces GeV-TeV protons.","keywords":["supernova remnants","IC 443","Fe I K-alpha line","low-energy cosmic rays","MeV protons","molecular clouds","X-ray spectroscopy","cosmic-ray acceleration"],"falsifier":"Measure the spatial profile of the 6.4 keV line across Reg 2 with a higher-resolution X-ray instrument: if the escape model is right, the line brightness should decline at roughly the ionization-cooling distance of ~10 MeV protons (about 500 yr of travel at $n_H=730$ cm$^{-3}$), and the line should not extend smoothly far beyond the 4-arcmin blob; a detection of 6.4 keV emission substantially more extended than that, or a centroid that shifts with position toward the pulsar wind nebula, would require X-ray irradiation or a very different diffusion environment.","tokens_in":9915,"feed_emoji":"🔭","tokens_out":10338,"duration_ms":95206,"temperature":0.7,"pith_summary":"This paper reports the discovery of two bright blob-like enhancements of the 6.4 keV Fe I K-$\\alpha$ line in the supernova remnant IC 443, one in the northwest and one near the middle, both coinciding with molecular clouds that the remnant's shock is interacting with. The paper argues that the most plausible origin is inner-shell ionization of neutral iron by MeV protons accelerated at the SNR shock and then leaking into the molecular clouds, rather than emission from hot plasma or X-ray irradiation by the nearby pulsar wind nebula. The measured line intensities, $0.273 \\pm 0.119$ and $0.344 \\pm 0.277$ photons s$^{-1}$ cm$^{-2}$ sr$^{-1}$, agree with the intensity of $0.29^{+0.10}_{-0.11}$ photons s$^{-1}$ cm$^{-2}$ sr$^{-1}$ predicted by a cosmic-ray escape model when the accelerated proton spectrum is fixed by the GeV-TeV gamma-ray observations. If the interpretation is right, IC 443 co-accelerates MeV protons with the GeV and TeV protons whose pion-decay gamma-rays are already observed, providing a direct observational bridge between low-energy and high-energy cosmic-ray acceleration in supernova remnants.","feed_headline":"Two iron-line blobs tie IC 443 to MeV cosmic rays","feed_subtitle":"Measured 6.4 keV glow matches the prediction if the same shock accelerates MeV and GeV-TeV protons.","key_machinery":"The mechanism that carries the argument is inner-shell ionization of neutral iron: a MeV proton knocks out a K-shell electron from Fe, and the resulting vacancy is filled by an electron that emits the Fe I K-alpha photon at 6.4 keV. The paper separates proton bombardment from electron bombardment by the line's equivalent width, and it anchors the proton flux by fitting a cosmic-ray escape model to the observed gamma-ray spectrum. In that model, the momentum at which protons escape the SNR decreases as the shock expands, producing a break in the proton spectrum and hence a bend in the gamma-ray spectrum near ~1 GeV; the same model then predicts the 6.4 keV line intensity, and the measured value matches the prediction to within the quoted uncertainties.","core_discovery":"The paper's central claim is that the Fe I K-$\\alpha$ line in IC 443, enhanced in two blob-like regions where the supernova remnant meets molecular clouds, is produced by low-energy cosmic-ray protons in the MeV band. Neutral iron atoms in the clouds are inner-shell ionized by these protons and emit at 6.4 keV; the large equivalent widths, >1.2 keV in Reg 1 and 0.7$^{+0.9}_{-0.6}$ keV in Reg 2, match the proton-bombardment expectation rather than the 0.2--0.4 keV expected for electron bombardment. The authors rule out X-ray irradiation from the pulsar wind nebula 1SAX J0617.1+2221 because the required flux would be two orders of magnitude above its observed flux. Applying a cosmic-ray escape model in which the same shock accelerates protons across energies, with the proton spectrum fixed by GeV and very-high-energy gamma-ray data, predicts a line intensity of $0.29^{+0.10}_{-0.11}$ photons s$^{-1}$ cm$^{-2}$ sr$^{-1}$, consistent with both measured regions; no significant line is found in a third region used as reference.","pith_inferences":["A testable extension is to map the 6.4 keV line at higher angular resolution: if the escape model is right, the line should trace the shock-cloud contact surface and fade within the roughly 500-yr ionization-cooling distance of 10 MeV protons, whereas a smoother or more extended profile would point to a larger diffusion coefficient inside the cloud than the model assumes.","The roughly 40-fold gap between the diffusion coefficient needed to explain the 4-arcmin size of Reg 2 and the coefficient used in the model could be resolved if future instruments detect a brightness gradient across the line-emitting cloud; such a gradient would measure the effective diffusion coefficient in the molecular-cloud environment.","The same method could be applied to other supernova remnants with both gamma-ray spectra and Fe I K-alpha detections to test whether co-acceleration of MeV and GeV-TeV protons is generic or peculiar to IC 443-like remnants."],"forward_implications":["MeV cosmic-ray protons are accelerated together with GeV and TeV protons at the IC 443 shock, meaning the same acceleration mechanism populates the low-energy cosmic rays that ionize interstellar gas.","Fe I K-alpha line imaging around molecular-cloud-interacting SNRs can act as a spatial probe of sub-GeV cosmic-ray escape, complementing gamma-ray and H3+ ionization-rate measurements.","The gamma-ray spectrum of IC 443 requires a smaller diffusion coefficient than those inferred for W28 and W44, so middle-aged SNRs may differ in how they release low-energy cosmic rays.","Because Reg 2 coincides with a dense molecular core and a gamma-ray clump, the data place escaping low-energy and high-energy cosmic rays in the same clouds, allowing a joint test of escape models."],"supporting_citations":[{"why":"Supplies the cosmic-ray escape model that the paper applies to IC 443 and from which the predicted Fe I K-alpha intensity is derived.","marker":"Makino et al. (2019)"},{"why":"Earlier X-ray analysis that detected the Fe I K-alpha line in the northeast of IC 443 and established the low-energy cosmic-ray interpretation being extended here.","marker":"Hirayama et al. (2019)"},{"why":"Provides the equivalent-width values that distinguish electron bombardment (0.2-0.4 keV) from proton bombardment (>0.6 keV).","marker":"Dogiel et al. (2011)"},{"why":"Reports the flux and morphology of the pulsar wind nebula 1SAX J0617.1+2221 used to rule out X-ray irradiation of the clouds.","marker":"Bocchino & Bykov (2001)"},{"why":"Gives the Fe K-shell ionization cross sections whose energy dependence identifies ~10 MeV protons as the relevant ionizing particles.","marker":"Tatischeff et al. (2012)"},{"why":"Provides the diffusion model for cosmic rays escaping a supernova remnant, including the form of the diffusion coefficient used in the fitting.","marker":"Ohira et al. (2011)"},{"why":"Source of the ionization-cooling time estimate for 10 MeV protons used in the Appendix to compare the observed emission-region size with the predicted diffusion.","marker":"Mannheim & Schlickeiser (1994)"},{"why":"GeV gamma-ray observations of IC 443 that constrain the high-energy proton spectrum in the escape-model fit.","marker":"Ackermann et al. (2013)"},{"why":"Very-high-energy gamma-ray data used together with other gamma-ray measurements to fix the accelerated proton spectrum.","marker":"Acciari et al. (2009)"}],"fun_headline_variants":["IC 443's iron glow points to low-energy cosmic rays","MeV protons leave iron fingerprint in supernova remnant","Fe K-alpha bumps in IC 443 hint at MeV proton acceleration","New iron-line maps tie IC 443's cloud shocks to MeV cosmic rays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that enough ~10 MeV protons can travel from the shock into the molecular clouds and spread over the observed 4-arcmin region before they lose energy by ionizing gas, and the paper's own diffusion estimate falls short of this requirement by about a factor of 40, with the gap handled by invoking unmodelled cloud geometry or extended injection.","fun_headline_variants_meta":{"raw":{"variants":["IC 443's iron glow points to low-energy cosmic rays","MeV protons leave iron fingerprint in supernova remnant","Fe K-alpha bumps in IC 443 hint at MeV proton acceleration","New iron-line maps tie IC 443's cloud shocks to MeV cosmic rays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00049,"raw_usage":{"total_tokens":2398,"prompt_tokens":922,"completion_tokens":1476,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":1403}},"tokens_in":538,"tokens_out":1476,"duration_ms":10251,"temperature":1.0,"reasoning_tokens":1403,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:23:14.185812+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spatial profile of the 6.4 keV line across Reg 2 with a higher-resolution X-ray instrument: if the escape model is right, the line brightness should decline at roughly the ionization-cooling distance of ~10 MeV protons (about 500 yr of travel at $n_H=730$ cm$^{-3}$), and the line should not extend smoothly far beyond the 4-arcmin blob; a detection of 6.4 keV emission substantially more extended than that, or a centroid that shifts with position toward the pulsar wind nebula, would require X-ray irradiation or a very different diffusion environment.","supporting_citations":[{"cited_title":"K., Matsumoto, & H., Ohira, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the cosmic-ray escape model that the paper applies to IC 443 and from which the predicted Fe I K-alpha intensity is derived."},{"cited_title":"2011 PASJ, 63, 535","cited_arxiv_id":null,"evidence_quote":"Provides the equivalent-width values that distinguish electron bombardment (0.2-0.4 keV) from proton bombardment (>0.6 keV)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the flux and morphology of the pulsar wind nebula 1SAX J0617.1+2221 used to rule out X-ray irradiation of the clouds."},{"cited_title":"2012, A&A, 546, A88","cited_arxiv_id":null,"evidence_quote":"Gives the Fe K-shell ionization cross sections whose energy dependence identifies ~10 MeV protons as the relevant ionizing particles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the diffusion model for cosmic rays escaping a supernova remnant, including the form of the diffusion coefficient used in the fitting."},{"cited_title":"\\ 2013, Science, 339, 807","cited_arxiv_id":null,"evidence_quote":"GeV gamma-ray observations of IC 443 that constrain the high-energy proton spectrum in the escape-model fit."},{"cited_title":"A., et al","cited_arxiv_id":null,"evidence_quote":"Very-high-energy gamma-ray data used together with other gamma-ray measurements to fix the accelerated proton spectrum."}],"review_version":1}