{"id":"e27bef50-738d-40ca-aea4-aab99dea3b1c","arxiv_id":"2506.14991","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"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 be detectable by COSI and next-generation MeV instruments, giving a direct…","lead":"This paper predicts specific gamma-ray lines that could reveal whether a supernova explosion made heavy r-process elements, using a magnetorotational supernova model. It identifies tin-126 in old remnants and tellurium/iodine lines in a future nearby explosion as detectable by upcoming MeV gamma-ray telescopes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Remnant detectability maps ignore the Doppler broadening the paper itself computes, comparing prompt narrow-line photon counts to line sensitivities; this overstates the 126Sb and 60Co remnant claims and is the most load-bearing weakness.","rationale":"The reader's weakest_assumption is the 0.3 Msun infalling-blob re-ejection that doubles the 60Fe yield. That is a real modeling choice, but it affects the 60Co remnant claims rather than the paper's headline 126Sn/126Sb signature. The more load-bearing problem for the central claim is the inconsistency between computing Doppler broadening and then ignoring it in the remnant sensitivity comparison: the right panels of Fig. 6 and the maps in Figs. 4-5 compare prompt, unbroadened line fluxes to narrow-line sensitivities. Since the line broadening is computed with the same ejecta velocity used throughout, this is an internal inconsistency rather than a disagreement with outside consensus. It is also easily testable: recomputing with broadened profiles and the tabulated sensitivities could move Vela Junior below COSI threshold and reduce the number of next-generation detections, although the qualitative conclusion that 126Sn is a distinctive r-process signature would likely survive for the nearest remnants. I therefore keep the reader's CONDITIONAL verdict unchanged: the paper is a useful feasibility study whose quantitative remnant lists need revision before being used for observation planning.","tokens_in":34319,"tokens_out":6092,"duration_ms":65666,"concrete_test":"Recompute the Fig. 4 and Fig. 5 detectability maps and the Fig. 6 right panels using the Doppler-broadened line profile from Eq. B9, integrated over each instrument's energy resolution (equivalently, require the total line flux to exceed sqrt(N) times the narrow-line sensitivity, with N = line FWHM / detector energy-resolution FWHM), and use the Table 4 sensitivity values. If Vela Junior/COSI or any next-generation remnant in the list falls below threshold, the remnant-based claims must be scaled back accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper computes Doppler broadening with vej = 0.015c in Sec. 2.4 and Eq. B9, and Fig. 6 left panels show the resulting broadened remnant spectra (line widths of roughly 2 percent around 666 keV and 1.17/1.33 MeV). Yet the remnant detectability maps in Figs. 4 and 5 and the right panels of Fig. 6 compare the total prompt photon flux in each line against published narrow-line sensitivities, implicitly assuming an unresolved line confined to one detector resolution bin. For a line broadened to N independent resolution elements, the total flux needed for a fixed significance grows roughly as sqrt(N); with N = 3-5 for COSI and next-generation detectors, sources sitting near threshold (Vela Junior with COSI, and possibly Crab, IC443, or others) could drop below detectability. Also, the COSI 126Sb sensitivity adopted in Sec. 3.1 (4e-6 ph/cm2/s) disagrees with the tabulated 4.2e-7 ph/cm2/s in Table 4, so the input sensitivities are internally inconsistent as well. Both effects bias the remnant observability claims upward, and they affect the primary 126Sn/126Sb signature, not just the 60Fe yield assumption. The qualitative conclusion that 126Sb is a distinctive r-process signature may survive, but the specific list of detectable remnants is not secure until the comparison is redone with the actual broadened line profiles.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the 35OC-RS magnetorotational supernova (MR-SN) model of Reichert et al. (2021) together with the PRISM network to compute gamma-ray spectra from radioactive decays in the ejecta, for both a future Galactic supernova and known supernova remnants. It predicts that second-peak r-process isotopes such as 126Sb (from 126Sn) and neutron-rich isotopes such as 60Co (from 60Fe) would be detectable in roughly eleven to twelve known Galactic remnants with COSI or next-generation MeV instruments, and that a future Galactic MR-SN would show r-process lines at 10 days (132Te, 131I, 132I), 100 days (103Ru, 106Rh, 95Nb), and 6 years (125Sb, 194Ir). The authors compare their predicted fluxes with published line and continuum sensitivities for INTEGRAL/SPI, COSI, AMEGO, and GRAMS, and argue that a detection of 126Sb would be a distinctive signature of r-process nucleosynthesis in a supernova.","tokens_in":34522,"tokens_out":4533,"duration_ms":43231,"significance":"If the predictions are correct, the paper offers a concrete, isotope-specific observational program for identifying r-process production in core-collapse supernovae, with falsifiable predictions for upcoming MeV instruments. The study has clear strengths: it builds on a published MR-SN simulation with tracer-based nucleosynthesis, validates the resulting abundance pattern against external metal-poor star data (Fig. 2b), explicitly propagates nuclear-mass uncertainties, and includes a detailed treatment of absorption and Doppler broadening in Appendix B. The qualitative conclusion that 126Sb and other second-peak lines are distinctive r-process signatures is well motivated. However, as detailed below, the specific remnant detectability claims are not yet secure because the comparison with line sensitivities is made inconsistently with the Doppler-broadened spectra the paper itself computes, and because a central 60Fe yield rests on an unverified infall-re-ejection assumption.","major_comments":[{"comment":"The remnant detectability maps compare the total prompt photon flux in each line against published narrow-line sensitivities, even though the paper itself computes Doppler broadening of these lines (Eq. B9, Sec. 2.4, left panels of Fig. 6). A line broadened to roughly 2% energy width, as shown for 126Sb and 60Co, will be spread over several independent detector resolution bins; for a fixed significance the required total flux increases approximately as the square root of the number of resolution elements. Repeating the comparison with the broadened spectra would move sources near threshold (e.g., Vela Junior with COSI, and possibly Crab, IC443, or others) below detectability. The specific list of detectable remnants is therefore not secure until the maps in Figs. 4 and 5 are recomputed with the broadened line profiles and an explicit spectral extraction significance.","section":"Sec. 2.4, Fig. 6, Figs. 4-5"},{"comment":"The adopted COSI line sensitivity for the 666/695 keV 126Sb lines is internally inconsistent: Sec. 3.1 uses Sl = 4e-6 ph/cm2/s, while Table 4 lists 4.2e-7 ph/cm2/s (with two years of observation time) for the same lines. This factor-of-ten discrepancy directly changes the COSI detection range quoted in Fig. 4 and the claim that Vela Junior and Vela SNR are above COSI sensitivity. Please state which value is correct, normalize observation times consistently, and re-evaluate the affected remnant conclusions.","section":"Sec. 3.1 and Table 4"},{"comment":"The central remnant claim for 60Co (twelve remnants, Fig. 5 top) relies on the 60Fe yield of 8.3e-3 M_sun, which is not a simulated result but an estimate obtained by assuming that the 0.3 M_sun infalling blob (Fig. 1, right panel) is re-ejected in the jet and follows the same nucleosynthesis as the existing jet ejecta except for the main r-process. The footnote to Table 5 states that the base simulation gives 4e-3 M_sun. If the blob instead accretes onto the protoneutron star, the 60Co fluxes roughly halve, removing the COSI detection of Vela Junior and some next-generation detections. The paper should either present this as a model-dependent range with the lower-yield case shown explicitly, or support the infall-re-ejection assumption with a dynamical argument or simulation.","section":"Sec. 2.1 and Appendix D, Table 5 footnote"},{"comment":"The quantitative remnant counts are based on single adopted values from often wide published ranges (e.g., IC443 age 3000-30000 yr adopted as 30000 yr; Vela Junior distance 0.5-1 kpc adopted as 0.7 kpc; G11.2-0.3 age 1400-2400 yr adopted as 2400 yr). Because the detectability boundary in Figs. 4 and 5 is steep in distance and, for 60Co, in age, these point choices directly affect whether individual remnants appear above threshold. The paper should show the remnant positions with error bars or as ranges, and state how many remnants remain detectable under the extreme ends of the adopted intervals.","section":"Table 1 and Figs. 4-5"}],"minor_comments":[{"comment":"The top panel caption reads 'MN-SN model'; this should be 'MR-SN model'.","section":"Fig. 5 caption"},{"comment":"The text refers to 'Jellyfish SNR (IC441)'; the standard name used elsewhere in the paper is IC443.","section":"Sec. 3.1"},{"comment":"The text states that 44Ti has a half-life of 60 days; its half-life is 59.1 years, and this typo should be corrected to avoid confusing 44Ti with 44Sc.","section":"Sec. 4"},{"comment":"The reference to 'Tables C and C in Appendix C' should be to the numbered tables (Tables 3 and 4) for the compiled continuum and line sensitivities.","section":"Sec. 2.5"},{"comment":"The sentence referring to 'SN 1987' should read 'SN 1987A'.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read of 2506.14991. It's the most concrete attempt I've seen to turn MR-SN r-process predictions into a MeV gamma-ray observing strategy for COSI and successors. The 126Sn/126Sb remnant signature is genuinely new, and the paper gives a prioritized list of remnants and epochs for a future Galactic event. The pipeline is standard—network plus spectra plus published sensitivities—but the authors do it carefully, compare their abundance pattern to metal-poor star data, and are transparent about nuclear uncertainties and the limits of the 2D simulation.\n\nThe soft spots are real but mostly localized. The load-bearing one is that the remnant detectability maps (Figs. 4, 5, and right panels of Fig. 6) compare the total prompt line flux to narrow-line sensitivities, even though the paper itself computes Doppler broadening of about 2% (Sec. 2.4). A line broadened over several resolution elements needs roughly sqrt(N) more flux for the same significance. Sources sitting near threshold—Vela Junior for COSI, maybe Crab, IC443, or others for next-gen—could drop below detectability. The qualitative point that 126Sb would be a distinctive r-process signature probably survives, but the twelve-remnant list is not secure until the comparison is redone with the actual broadened profiles.\n\nSecond: the 60Fe yield story. The 8.3e-3 Msun adopted in the text comes from adding 0.3 Msun of infalling material to the 4e-3 Msun from Reichert et al. 2021, under the assumption that the blob is re-ejected and follows jet nucleosynthesis. That's a modeling choice, not a simulated result, and the dozen-remnant 60Co claim leans on it. The paper does flag the 100% ejection assumption as an upper limit, and the table gives the range, so it's not hidden—but it should be treated as a systematic uncertainty.\n\nThird, and minor in comparison: the COSI and GRAMS line sensitivities used in Sections 3.1-3.2 don't match the appendix tables by factors of 6-10. Some are conservative, some optimistic, so the net direction isn't uniform, but it needs to be reconciled.\n\nThe future-supernova prompt spectra (10 days, 100 days, 6 years) are more robust, since they use continuum sensitivities and do include the reprocessed/broadened spectra. That part of the paper is solid and will be genuinely useful for planning.\n\nVerdict: send it to a serious referee. The core idea is worth engaging, and the flaws are fixable. I'd want the remnant comparison redone with broadened profiles and the sensitivity tables cleaned up before I'd trust the target list, but the paper earns a full review, not a desk reject.","headline":"A concrete, useful target list for COSI, but the remnant detectability maps ignore the Doppler broadening the paper itself computes, so the specific targets aren't yet secure.","tokens_in":35208,"tokens_out":3592,"would_cite":true,"duration_ms":34854,"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":"Gamma-ray lines could prove supernovae make r-process elements","keywords":["r-process nucleosynthesis","magnetorotational supernovae","MeV gamma-ray astronomy","supernova remnants","126Sn and 126Sb gamma-ray lines","60Fe and 60Co","COSI","next-generation MeV detectors"],"falsifier":"Observe the Vela region with COSI or a next-generation MeV telescope: the paper predicts 126Sb lines at 666 and 695 keV and 60Co lines at 1173 and 1332 keV with fluxes above SPI's existing 60Fe upper limit of about $1.1\\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$, so a targeted search that fails to find these lines in Vela Junior or Vela SNR would directly contradict the model's central prediction. A second, independent check would be to measure a remnant's 60Fe mass from the 1332 keV line and find it closer to standard CCSN yields (about $3.7\\times10^{-5}$ solar masses) than to the MR-SN yield ($8.3\\times10^{-3}$ solar masses).","tokens_in":33998,"feed_emoji":"🔭","tokens_out":6256,"duration_ms":56185,"temperature":0.7,"pith_summary":"The paper argues that gamma-ray observations can settle whether core-collapse supernovae produce r-process elements. Using a two-dimensional magnetorotational supernova model that ejects neutron-rich jet material, it computes the MeV gamma-ray spectra from radioactive decays over timescales of days to millenia. It claims that the 666 and 695 keV lines of 126Sb, fed by long-lived 126Sn, would be detectable in about eleven known supernova remnants and are a distinctive signature of r-process production. It also claims that enhanced 60Fe, seen as 60Co lines at 1173 and 1332 keV, would indicate the neutron-rich conditions favorable to the r-process, and that a future Galactic MR-SN would show second-peak and weak r-process lines at 10 days, 100 days, and 6 years. If detected, these lines would be direct evidence of r-process nucleosynthesis in supernovae; non-detection would constrain the yields.","feed_headline":"Gamma-ray lines could prove supernovae make r-process elements","feed_subtitle":"COSI and next-generation MeV telescopes could spot 126Sb and 60Co in eleven known remnants.","key_machinery":"The load-bearing object is the 35OC-RS magnetorotational supernova model, a two-dimensional simulation with neutrino transport whose jet ejects neutron-rich material. Nucleosynthesis along tracer trajectories and shock-heated outflow is evolved with the PRISM reaction network; the gamma-ray spectra are built from beta-decay flows using ENDF/B-VIII.0 and ENSDF decay intensities. The emitted photons are then propagated through an expanding, mixed ejecta with Doppler broadening and absorption, and a geometric estimate of positron annihilation adds the 511 keV and positronium contributions. Observability is set by published line and continuum sensitivities of INTEGRAL/SPI, COSI, and next-generation instruments such as GRAMS and AMEGO, scaled to one-year or burst observing times. This chain converts a hydrodynamical nucleosynthesis model into concrete predictions of detectable lines in specific remnants.","core_discovery":"The central claim is that an observation of 126Sb gamma rays in a supernova remnant stands out as a signature of an r-process-producing supernova. For the 35OC-RS magnetorotational supernova model, the paper predicts $2.5\\times10^{-4}$ solar masses of 126Sb, enough to place roughly eleven known remnants above the line sensitivities of COSI or next-generation MeV telescopes at their adopted ages and distances. It further predicts that the same neutron-rich jet boosts the 60Fe yield to $8.3\\times10^{-3}$ solar masses, making the 60Co 1173 and 1332 keV lines detectable in all twelve remnants considered and distinguishing MR-SNe from standard neutrino-driven supernovae at distances beyond a few kiloparsecs. For a future Galactic explosion, the calculation shows that second-peak r-process isotopes (132Te, 131I, 132I) appear above backgrounds at ten days, weak r-process isotopes (103Ru, 106Rh, 95Nb) at one hundred days, and 125Sb with the third-peak isotope 194Ir at about six years. The paper presents this as a concrete observational test of whether any core-collapse supernova channel makes r-process elements.","pith_inferences":["One could apply the same flux calculation to other proposed r-process supernova channels, such as collapsars; any neutron-rich outflow that makes second-peak nuclei should produce similar 126Sb and 125Sb lines, so the detection strategy is broader than magnetorotational supernovae alone.","A longer-duration magnetohydrodynamic simulation that follows the infalling blob through re-ejection would directly test the 60Fe yield; until then, the twelve-remnant 60Co count should be read with that modeling choice in mind.","Existing archival data from SPI may already be able to search for the 126Sb 666 keV line in Vela or Vela Junior, since the paper notes SPI's Vela 60Fe upper limit is an order of magnitude below the predicted flux; a re-analysis could provide a near-term test before COSI launches.","If both 126Sb and 60Co are seen together in the same remnant, the ratio of their fluxes would measure the neutron richness of the ejected jet material and could be compared with abundance patterns in metal-poor stars."],"forward_implications":["A positive detection of the 126Sb lines in any listed remnant would be direct evidence that the r-process reached the second peak in that supernova.","A non-detection at the predicted fluxes would place quantitative upper limits on r-process and 60Fe yields, constraining MR-SN models and possibly ruling out the magnetorotational origin for that remnant.","Detecting 60Co beyond about 4–5 kpc with next-generation instruments, or beyond 0.3 kpc with COSI, would favor an MR-SN over a standard neutrino-driven explosion, while a modest 60Co abundance in a nearby remnant would favor a standard CCSN.","Combined measurements of 126Sb and 60Co in one remnant would discriminate among explosion mechanisms: substantial 60Co without 126Sb suggests a neutron-rich but second-peak-failing outflow, and 126Sb without 60Co would motivate exotic supernova models.","A future Galactic MR-SN would show a time-ordered sequence of r-process lines (second peak at 10 days, weak r-process at 100 days, 125Sb and 194Ir at 6 years) that could measure the weak-to-second-peak ratio and how far the neutron-capture chain extended."],"supporting_citations":[{"why":"Supplies the 35OC-RS MR-SN model and its tracer trajectories, the foundation of the nucleosynthesis calculation.","marker":"Reichert et al. (2021)"},{"why":"Provides the two-dimensional MR-SN simulation setup from which the model is drawn.","marker":"Obergaulinger & Aloy (2017)"},{"why":"Provides the PRISM reaction network used to evolve abundances from the ejecta thermodynamics.","marker":"Sprouse et al. (2021)"},{"why":"Gives the gamma-ray spectral and radiative-transfer method used to compute prompt and processed spectra.","marker":"Wang et al. (2020)"},{"why":"First highlighted long-lived r-process isotopes as gamma-ray candidates in CCSNe and motivates the remnant search.","marker":"Qian et al. (1998)"},{"why":"Supplies the standard CCSN 60Fe yield that defines the non-MR-SN comparison in Fig. 5.","marker":"Limongi & Chieffi (2018)"},{"why":"Provides the SPI upper limit on 60Fe in the Vela region that the model's Vela prediction is compared against.","marker":"Wang et al. (2007)"},{"why":"Provides COSI line and continuum sensitivities used to define detectability.","marker":"Tomsick et al. (2023)"},{"why":"Provides GRAMS sensitivities used as the next-generation instrument benchmark.","marker":"Aramaki et al. (2020)"}],"fun_headline_variants":["126Sb gamma rays unmask r-process supernovae","MeV telescopes can find r-process in supernova remnants","Tin-126 signal proves supernovae make heavy elements","Future supernova bursts reveal r-process isotopes","Gamma-ray signatures trace r-process in explosions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim depends on the assumption that roughly 0.3 solar masses of infalling material seen in the simulation will be re-ejected by the jet and follow the same nucleosynthesis as the existing jet ejecta; this choice nearly doubles the 60Fe yield and is what puts twelve remnants above the 60Co detection threshold.","fun_headline_variants_meta":{"raw":{"variants":["126Sb gamma rays unmask r-process supernovae","MeV telescopes can find r-process in supernova remnants","Tin-126 signal proves supernovae make heavy elements","Future supernova bursts reveal r-process isotopes","Gamma-ray signatures trace r-process in explosions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1584,"prompt_tokens":1146,"completion_tokens":438,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":762,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":762,"tokens_out":438,"duration_ms":4667,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:46:36.815686+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the Vela region with COSI or a next-generation MeV telescope: the paper predicts 126Sb lines at 666 and 695 keV and 60Co lines at 1173 and 1332 keV with fluxes above SPI's existing 60Fe upper limit of about $1.1\\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$, so a targeted search that fails to find these lines in Vela Junior or Vela SNR would directly contradict the model's central prediction. A second, independent check would be to measure a remnant's 60Fe mass from the 1332 keV line and find it closer to standard CCSN yields (about $3.7\\times10^{-5}$ solar masses) than to the MR-SN yield ($8.3\\times10^{-3}$ solar masses).","supporting_citations":[],"review_version":2}