{"id":"46cc8fc1-0285-4154-ba0f-3b5decf5fbcd","arxiv_id":"2608.10451","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":8,"one_line_summary":"First JWST mid-infrared spectra of a normal Type Ia supernova from maximum light reveal strongly stratified ejecta with stable nickel innermost, cobalt at intermediate radii, and argon outside, plus resolved ~800 km/s substructure in [Ca IV] emission.","lead":"JWST captured the first mid-infrared spectrum of a normal Type Ia supernova at maximum brightness, plus follow-up spectra at 23 and 84 days. The data show a layered explosion with stable nickel innermost, radioactive cobalt above it, and argon in the outer ejecta, along with small-scale clumping in calcium emission.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The [CoIII] flat-top interpreted as a central Co deficit may instead be collisional de-excitation in the dense core; the paper's own critical-density argument predicts this, so the headline stratification is not uniquely established.","rationale":"The reader's weakest assumption (optically thin, symmetric line inversion) points at the same machinery, but my concern is a specific, physical failure mode that the paper itself introduces in Section 3.1 and never closes: the low critical density of [CoIII]. If the inner core is still collisionally de-excited at +84 d, the flat top is a line-formation artifact, not a Co hole. This directly affects the most novel physical conclusion, the strong Ni/Co/Ar stratification, and the quantitative OVL and v50 metrics derived from it. The data quality, public release, and honest caveats elsewhere (e.g., the acknowledged Mg II identification uncertainty) support accepting the observational sequence, but the Co-hole inference needs a dedicated density/critical-density check before the stratification claim can be taken as established. I therefore recommend CONDITIONAL rather than unconditional ACCEPT: add the calculation (or reframe the Co inner edge as a line-formation effect) and the claim is secure.","tokens_in":29943,"tokens_out":14354,"duration_ms":152231,"concrete_test":"Compute n_e(v) and n_crit([CoIII] 11.88 µm), n_crit([NiIII] 11.00 µm) from the CMFGEN N100 model used in Section 5 at ~104 days, and overlay the v<2000 km/s region. If n_e(v<2000) is between n_crit([CoIII]) and n_crit([NiIII]) at +84 d, collisional quenching of [CoIII] is viable and the inferred Co hole is not established. The decisive forward test: rerun the 104-day CMFGEN calculation with the Co abundance kept constant below 2000 km/s and check whether the synthetic [CoIII] line develops a flat-topped core from density/critical-density effects alone; if it does, the stratification claim must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In Section 4.1, the flat-topped [CoIII] 11.88 µm core (v≲2000 km/s) is converted via Eq. (1) into zero emissivity there, and the paper concludes that radioactive Co is absent from the inner ~2000 km/s. This is load-bearing: it is the unique piece behind the 'Co absent within ~2000 km/s' part of the stratification claim and the Ni–Co OVL=0.56 value. The interpretation is not uniquely forced. Section 3.1 already invokes the ~5–6× lower critical density of [CoIII] relative to [NiIII] to explain the weak Co lines at +1 d: in a dense core, collisional de-excitation suppresses [CoIII] emission even where Co is abundant. The same mechanism operating at +84 d in the innermost, still-dense ejecta produces exactly a flat-topped profile: the line turns on first in lower-density outer regions and remains quenched inside, yielding the observed delayed emergence and central plateau. The paper rules out an opacity origin for the flat top (no continuum, sharp [ArIII]) but never rules out collisional de-excitation. If this is happening, Eq. (1) maps a line-formation effect to zero emissivity, and v50(Co), the 19% Co-flux-interior-to-Ni quantile, and OVL are all biased. A density threshold predicts the same morphology, so the observation does not uniquely demonstrate a central Co hole.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"SN 2025rbs is a normal but fast-declining Type Ia supernova in NGC 7331 observed with JWST NIRSpec/MIRI at +1, +23, and +84 days after B-band maximum, together with contemporaneous ground-based optical and NIR spectroscopy. The paper presents the first maximum-light MIR spectrum and the earliest MIR spectroscopic sequence of a normal SN Ia. It identifies the emerging permitted and forbidden lines, tracks the photospheric-to-nebular transition in the MIR, and uses line inversion at +84 days to infer strongly stratified ejecta: stable Ni concentrated at the lowest velocities, radioactive Co at intermediate velocities with an apparent deficit inside ~2000 km/s, and Ar in an outer shell. The paper also reports reproducible small-scale substructure in [Ca IV] 3.21 µm (inter-epoch correlation r=0.92, characteristic clump scale ~800 km/s) and finds that CMFGEN models underpredict the candidate MIR Mg II features despite approximately reproducing the NIR Mg II 1.0927 µm line.","tokens_in":30216,"tokens_out":14531,"duration_ms":138924,"significance":"The dataset is genuinely novel and will be a reference point for early-time MIR spectroscopy of SN Ia. The paper is careful in several respects: bootstrap error propagation, Monte Carlo and block-bootstrap checks of the [Ca IV] substructure, cross-checks between [Ca IV] and [Ca V], and explicit caveats about model dependence. The stratification and clumping results, if established, would provide genuinely new constraints on explosion models and on mixing in thermonuclear supernovae. The main caveat is that the 'Co absent inside ~2000 km/s' conclusion rests on interpreting the flat-topped [Co III] profile as a compositional hole rather than a line-formation effect; a quantitative test is required before this central claim can be considered secure.","major_comments":[{"comment":"The flat-topped [Co III] 11.88 µm profile is interpreted in §4.1 as the absence of Co emission for v ≲ 2000 km/s, and this inference drives v50(Co), the 19% Co-flux-interior-to-Ni quantile, and the Ni–Co OVL. However, §3.1 invokes the ~5–6× lower critical density of [Co III] relative to [Ni III] to explain the weakness of Co lines at +1 day. The same mechanism can produce a flat-topped profile at +84 days if the innermost ejecta remain near or above the [Co III] critical density: the line would turn on preferentially in lower-density outer layers, and the observed profile would show a central plateau even where Co is abundant. The paper rules out an opacity origin (no continuum, fully nebular [Ar III], peaked Ni lines), but it does not rule out this collisional line-formation effect. Please estimate the electron density profile in the inner ejecta at +84 days relative to n_crit([Co III]), or show quantitatively from the +23 and +84 day profiles that the flat-top inner edge is stationary in velocity (as expected for a compositional hole) rather than shrinking inward as density declines (as expected for a critical-density threshold). Without such a test, the 'Co absent within ~2000 km/s' statement is not uniquely established.","section":"§3.1 and §4.1"},{"comment":"The bootstrap uncertainties reported in Table 2 propagate spectral noise under the assumption that the flat-topped [Co III] plateau corresponds to exactly zero emissivity. They do not include the systematic uncertainty in that zero-emissivity assumption: Eq. (1) assigns j(v)=0 wherever dF/dv=0, so any residual plateau flux from fitting uncertainty, weak blending, or a line-formation effect is mapped to a complete central hole. Please quantify how v50(Co), the 19% flux fraction, and Ni–Co OVL change if the emissivity inside v_inner is allowed to be a small finite fraction (e.g., 3–10%) of the peak emissivity, or if the flat top is instead modeled as a density-dependent suppression. This robustness check is important because the stratification claim depends precisely on the zero-emissivity mapping of the plateau.","section":"§4.1 and Table 2"}],"minor_comments":[{"comment":"The text refers to the 'Phillip relations'; this should be the 'Phillips relations.'","section":"§2.3"},{"comment":"The term 'root-mean-square (RSM)' is a typo; it should be 'root-mean-square (RMS).'","section":"§4.2"},{"comment":"The phrase 'χ2/dof = 0.30 for two degrees of freedom' is ambiguous; please state the χ² value and the number of degrees of freedom separately.","section":"§4.2"},{"comment":"The critical-density comparison between [Co III] and [Ni III] would benefit from citing the atomic data source and stating the adopted collision strengths or references, since this argument is now central to the interpretation in §4.1.","section":"§3.1"},{"comment":"The text calls the [Co III] 11.88 µm line 'essentially isolated,' but Figure 4 marks contamination from [Co II], [S IV], and [Ni II] in parts of the profile; please clarify which wavelength ranges were used for the fit and inversion.","section":"§4.1 and Figure 4"},{"comment":"The caption mentions 'ACF and PSD scales,' but the text only defines the autocorrelation function; please either describe the power spectral density analysis or remove the PSD reference.","section":"Figure 6"}],"recommendation":"major_revision","confidential_remarks":"To the editor: this is a strong observational paper and I support publication after the authors address the [Co III] flat-top ambiguity. The reader's report recommends acceptance, but I regard the central stratification claim as currently not uniquely established because the paper's own critical-density discussion in §3.1 provides a plausible alternative mechanism that is not ruled out in §4.1. The requested quantitative density/evolution test is feasible within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper gives you the first maximum-light MIR spectrum and the earliest MIR spectroscopic sequence of a normal SN Ia, from JWST NIRSpec/MIRI at +1, +23, +84 d. That alone is a real advance. The data reduction is careful, the line identifications are guided by CMFGEN models, and the Ca substructure analysis is the most thorough I have seen on this topic: the +23/+84 d residual correlation (r≈0.92), the cross-check with [Ca V] (r≈0.66), and the ACF-to-clump-scale Monte Carlo calibration all hold together. The Mg II discrepancy between MIR and NIR is a useful, concrete challenge for radiative-transfer codes. The paper deserves a serious referee and will be widely cited as a benchmark dataset.\n\nThe soft spot is the Co stratification. The flat-topped [Co III] 11.88 µm profile is read as zero emissivity inside ~2000 km s−1, leading to the claim that radioactive Co is absent from the core. The paper rules out an opacity origin by pointing to the absence of continuum and the sharp [Ar III] edges. But collisional de-excitation is not opacity. In Section 3.1 the authors themselves invoke the ~5–6× lower critical density of [Co III] to explain its weakness at +1 d. At +84 d the innermost ejecta are still the densest part of the outflow, and nothing in the paper demonstrates that they have fallen below the [Co III] critical density. If the core is still collisionally quenched, Eq. (1) maps a line-formation effect into an emissivity hole, and the derived v50(Co), the 19% Co-flux-interior-to-Ni quantile, and the Ni–Co OVL are all biased. A density threshold produces exactly the observed flat top, so the observation does not uniquely establish a central Co deficit. This is a load-bearing part of the stratification story, not a footnote.\n\nThe other assumptions—optically thin MIR at +84 d, spherical symmetry in the inversion—are stated clearly, and the paper hedges the Mg II identification appropriately. The asymmetry of [Ar III] suggests the symmetry assumption is approximate at best.\n\nWho is this for? Anyone doing SN Ia observations, MIR spectroscopy, or radiative-transfer modeling of these explosions. It should go to peer review; the data and the Ca clump result are strong, and the Co interpretation needs either an additional test (e.g., density/ionization estimate for the core) or a more cautious claim.","headline":"Genuinely new JWST MIR data (first maximum-light spectrum of a normal SN Ia) and a careful, cross-checked analysis—but the central claim of a Co-free inner core is undercut by a collisional de-excitation alternative the paper does not rule out.","tokens_in":31150,"tokens_out":3630,"would_cite":true,"duration_ms":36901,"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":"The first maximum-light mid-infrared spectrum of a normal Type Ia supernova shows its ejecta were still strongly layered by day 84, with stable nickel deepest, radioactive cobalt above it, and argon forming an outer shell.","keywords":["Type Ia supernovae","mid-infrared spectroscopy","JWST observations","SN 2025rbs","ejecta stratification","nucleosynthesis","nebular spectra","clumpy ejecta"],"falsifier":"Observe SN 2025rbs again in the MIR at, say, +150 days or later. If the [Co III] 11.88 μm line fills in at its center or narrows as the density drops, then the flat-topped profile at +84 days was an opacity effect rather than a real absence of cobalt inside about 2000 km/s, and the stratified model would need revision. The clump interpretation would be falsified if higher-resolution [Ca IV] and [Ca V] spectra taken at a third epoch no longer show the same residuals at the same velocities, or if the residuals track a continuum or instrumental artifact rather than a fixed clump population.","tokens_in":29682,"feed_emoji":"🔭","tokens_out":10256,"duration_ms":85419,"temperature":0.7,"pith_summary":"This paper presents the first maximum-light mid-infrared spectrum and the earliest mid-infrared spectroscopic sequence of a normal Type Ia supernova, using JWST observations of SN 2025rbs at +1, +23, and +84 days past B-band maximum. It claims that these spectra trace the MIR ejecta from a continuum-plus-lines photospheric phase into a fully nebular phase, and that the +84 day forbidden lines directly map a strongly stratified interior: stable nickel sits at the lowest velocities, radioactive cobalt lies above it but is absent inside about 2000 km/s, and argon forms an outer shell. The paper also reports small ripples in the calcium line [Ca IV] 3.21 μm with a characteristic scale near 800 km/s, repeated at two epochs, which it reads as clumping or patchy ionization in the calcium-emitting regions. The claim matters because MIR spectroscopy near maximum opens a simultaneous view of the freshly emerging inner ejecta and the fading outer burning products, giving explosion models a new set of layers to reproduce.","feed_headline":"First maximum-light mid-IR spectra of a Type Ia reveal layered ejecta","feed_subtitle":"Three mid-infrared epochs map nickel, cobalt, and argon to separate velocity shells, with calcium clumping at ~800 km/s.","key_machinery":"The load-bearing machinery is the one-dimensional line-inversion identity $j(v)\\propto (1/v)\\,dF_\\nu/dv$, which turns a Doppler-broadened line profile into a radial emissivity for each ion; the paper applies it to the +84 day [Ni III], [Co III], and [Ar III] lines and then volume-weights the results, $p(v)\\propto j(v) v^2$, to define half-flux velocities and overlap coefficients. For the calcium substructure, the machinery is an autocorrelation function of the high-pass-filtered residuals, calibrated by Monte Carlo forward models of a clump-filled spherical shell to recover an intrinsic clump FWHM and a filling factor. The radiative-transfer comparisons are made with the CMFGEN non-local-thermodynamic-equilibrium code, used differentially (subtracting a single ion from the calculation) to identify which transitions dominate each observed feature.","core_discovery":"At +84 days, the MIR spectrum of SN 2025rbs is fully nebular, and its isolated forbidden lines yield, through a line-inversion formula $j(v)\\propto (1/v)\\,dF/dv$, velocity-resolved emissivities for different elements. Half of the [Ni III] 11.00 μm flux comes from below $5.0\\times10^3$ km/s, half of [Co III] 11.88 μm from below $7.0\\times10^3$ km/s, and half of [Ar III] 8.99 μm from below $11.6\\times10^3$ km/s; pairwise overlaps fall from 0.56 (Ni–Co) to 0.28 (Co–Ar) to 0.08 (Ni–Ar). The flat-topped Co profile requires the radioactive cobalt to be missing inside roughly 2000 km/s, implying the center is dominated by stable iron-group material. On a smaller scale, the [Ca IV] 3.21 μm profile shows 3–4% fluctuations with correlation lengths 530 and 590 km/s at +23 and +84 days, a strong correlation ($r=0.92$) between epochs, and an inferred intrinsic clump scale of 780 km/s; the fainter [Ca V] 4.16 μm line shows a correlated residual ($r=0.66$), supporting a physical origin. Independently, the paper identifies MIR features at 4.6, 9.6, and 13.0 μm as Mg II transitions and shows that radiative-transfer models reproduce the NIR Mg II 1.0927 μm line while underpredicting these MIR lines by a large factor, pointing to transition-dependent excitation rather than a simple abundance deficit.","pith_inferences":["This reader infers that the OVL(Ni, Co) overlap coefficient could be measured across a sample of normal Type Ia supernovae; if it varies systematically with decline rate or peak luminosity, it would give an empirical axis for separating delayed-detonation from sub-Chandrasekhar double-detonation models.","This reader infers that the ~800 km/s calcium clump scale may be a fossil of the deflagration plume scale, and because homologous expansion preserves velocity, future high-signal-to-noise [Ca IV] spectra of many events can test whether this scale is universal or depends on explosion properties.","This reader infers that the paper's Mg II comparison suggests a sharper test: a model with the correct Mg ionization balance should reproduce the ratio of MIR to NIR Mg II flux without rescaling, and observing the same four lines in another fast-declining supernova would show whether the underprediction is event-specific or a property of the model grids.","This reader infers that confirming the [Ca V]/[Ca IV] correlation with better signal to noise would help distinguish patchy ionization from abundance clumping, because [Ca V] requires a higher ionization energy and would sample different spatial regions if ionization varies."],"forward_implications":["In the MIR, SN 2025rbs is fully nebular by +84 days even though its optical and near-infrared spectra are still transitional, so MIR observations can certify the onset of the nebular phase earlier than shorter-wavelength data.","The separation of Ni, Co, and Ar into distinct radial layers, with essentially no Ar inside the Ni 90%-flux velocity and no radioactive Co inside about 2000 km/s, implies limited macroscopic mixing in this explosion.","If the Ca substructure is physical, it places roughly 5100 clumps of characteristic size 780 km/s filling about 10% of the adopted emitting shell, and the persistent two-epoch pattern means the clumps remain fixed in velocity during expansion.","The Mg II discrepancy identifies a specific model failure: the MIR recombination-cascade transitions and the NIR photoexcited 1.0927 μm line require different ionization and excitation treatments, so a successful model must reproduce both simultaneously.","Observations beginning at maximum light are necessary to catch the transient MIR Mg II features and the early emergence of forbidden [Ni III] and [Ar III], which vanish or become blended within weeks."],"supporting_citations":[{"why":"Source of the line-inversion formula $j(v)\\propto (1/v)\\,dF/dv$ used to convert each nebular line profile into an emissivity.","marker":"C. Fransson & R. A. Chevalier 1989"},{"why":"Provides the fitting, bootstrap, and line-inversion methodology, and the comparison [Ca IV] profiles of SN 2022aaiq and SN 2024gy that help rule out fixed instrumental patterns.","marker":"L. A. Kwok et al. 2026"},{"why":"Supplies the CMFGEN atomic data, differential single-ion line lists, and the relativistic-tilt discussion used for line identification and profile analysis.","marker":"S. Blondin et al. 2023"},{"why":"Supplies the N100 delayed-detonation ejecta structure used for the baseline radiative-transfer models.","marker":"I. R. Seitenzahl et al. 2013"},{"why":"Supplies the MERGER 1109 violent-merger ejecta structure used to test whether higher Mg abundance and lower ionization explain the MIR Mg II features.","marker":"R. Pakmor et al. 2012"},{"why":"Provides the statistical clump framework and filling-factor equation from which the clump scale, volume filling factor, and clump count are derived.","marker":"N. N. Chugai 1994"},{"why":"Established the set of MIR forbidden lines expected from Type Ia supernovae, against which the +84 day identifications are checked.","marker":"C. L. Gerardy et al. 2007"},{"why":"Obtained the only earlier SN Ia MIR spectral sequence (SN 2014J, +37 to +117 days), the baseline this paper extends to maximum light.","marker":"C. M. Telesco et al. 2015"},{"why":"The CMFGEN code that generates the model spectra used throughout Section 5 for line identification.","marker":"D. J. Hillier & L. Dessart 2012"}],"fun_headline_variants":["First max-light mid-IR spectra dissect Type Ia layers","JWST mid-IR maps stable nickel core, cobalt ring, argon shell","JWST mid-IR unveils stratified ejecta in Type Ia supernova","First MIR sequence maps nickel to argon shells in SN Ia"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The stratification and clump measurements assume the +84 day mid-infrared lines are optically thin emitters viewed in a spherically symmetric outflow, so that a line profile can be inverted directly into a radial emissivity; significant residual opacity or an asymmetric geometry would shift the inferred velocity shells and clump sizes.","fun_headline_variants_meta":{"raw":{"variants":["First max-light mid-IR spectra dissect Type Ia layers","JWST mid-IR maps stable nickel core, cobalt ring, argon shell","JWST mid-IR unveils stratified ejecta in Type Ia supernova","First MIR sequence maps nickel to argon shells in SN Ia"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001032,"raw_usage":{"total_tokens":4509,"prompt_tokens":1272,"completion_tokens":3237,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":888,"completion_tokens_details":{"reasoning_tokens":3163}},"tokens_in":888,"tokens_out":3237,"duration_ms":19763,"temperature":1.0,"reasoning_tokens":3163,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:20:12.062135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe SN 2025rbs again in the MIR at, say, +150 days or later. If the [Co III] 11.88 μm line fills in at its center or narrows as the density drops, then the flat-topped profile at +84 days was an opacity effect rather than a real absence of cobalt inside about 2000 km/s, and the stratified model would need revision. The clump interpretation would be falsified if higher-resolution [Ca IV] and [Ca V] spectra taken at a third epoch no longer show the same residuals at the same velocities, or if the residuals track a continuum or instrumental artifact rather than a fixed clump population.","supporting_citations":[{"cited_title":"A., Liu, C., Jha, S","cited_arxiv_id":null,"evidence_quote":"Provides the fitting, bootstrap, and line-inversion methodology, and the comparison [Ca IV] profiles of SN 2022aaiq and SN 2024gy that help rule out fixed instrumental patterns."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the statistical clump framework and filling-factor equation from which the clump scale, volume filling factor, and clump count are derived."}],"review_version":1}