{"id":"c9f56690-db63-4d05-81c6-a48cc3cbd29e","arxiv_id":"2505.06609","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A unified MCMC fit of 32 radio supernovae shows higher mass-loss for stripped-envelope progenitors and attributes apparent shock-acceleration anomalies to a dense inner circumstellar medium.","lead":"Researchers fit a standard radio emission model to 32 supernovae and found that stripped-envelope supernovae have about ten times higher inferred mass-loss rates than normal hydrogen-rich ones. The study also suggests that surprisingly strong magnetic field efficiencies and flat ejecta slopes are artifacts of a missing dense shell of gas close to the star.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass-loss gap is driven by the 100x adopted wind-velocity ratio, not by the fitted CSM-density scale (Eq. 4, Table 5).","rationale":"I agree with the reader that the paper is a careful systematic study and that CONDITIONAL is the right verdict. However, I identify a different weakest link. The reader's strongest_claim is the mass-loss difference, and the reader's rationale mentions wind-velocity assumptions, but their weakest_assumption focuses on the inner-CSM contamination of late-time parameters. In my reading, the more direct threat to the headline is the conversion from A~ to Mdot. Because Mdot is proportional to vw and the posterior medians of A~ are within 0.3 dex of each other (Table 5), the factor ~50 difference in the abstract is essentially the chosen vw ratio. The dense inner CSM explanation (Section 5.2) is an interpretation of the early-phase discrepancy; while plausible, it concerns the secondary claim. The mass-loss-difference claim is quantitative and central, and is the one most in need of a sensitivity test. The concrete check I propose would settle it: if the ratio of mass-loss posteriors is robust to a factor of ~3–10 change in the assumed SESN wind speed, the claim stands; otherwise it should be reframed as a difference in Mdot/vw (i.e., A~), with the physical mass-loss difference conditional on stellar wind physics. This does not move the verdict away from CONDITIONAL, but it sharpens the required revision: the abstract and summary should either state the wind-velocity assumption prominently or present the A~ (normalized mass-loss) result as the primary data-driven finding.","tokens_in":32899,"tokens_out":14826,"duration_ms":146068,"concrete_test":"Regenerate the mass-loss posteriors from the MCMC chains using the same A~ and s samples but with wind velocity treated as a nuisance parameter: draw vw from log-normal priors representing red supergiants (median 10 km/s) and stripped stars (median 1000 km/s, plus a conservative median of 100 km/s). Compute the 95% credible interval for the ratio Mdot_SESN/Mdot_II. If the interval includes 10 or below, the order-of-magnitude claim is not supported by the data. A simpler sensitivity check: recompute Figure 4 with vw(SESN)=300 km/s while keeping vw(II)=10 km/s; if the two mass-loss contours then overlap, the central claim is an artifact of the assumed velocity ratio rather than a robust inference.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract, §8) that SESN progenitors have mass-loss rates an order of magnitude above SN II rests on converting the fitted CSM density scale A~ into Mdot via Mdot = 4π vw × 5×10^11 g cm^-1 × A~ (Eq. 4). Across the successful fits, the median log A~ is 1.75 for SNe II and 1.47 for SESNe (Table 5): the data-constrained quantity A~ = Mdot/vw is essentially the same (slightly higher for SNe II). The reported factor ~50–100 in Mdot therefore comes almost entirely from the adopted wind velocities, 10 km/s for SNe II and 1000 km/s for SESNe (§5.1), not from the radio light curves. With only a handful of SNe II (N≈3–4) and 1σ scatter of ~1.5 dex in log A~, the difference in A~ is not statistically significant. Section 5.1 acknowledges that 'the choice of the CSM velocity involves uncertainty,' but this limitation is not propagated into the abstract claim or tested with alternative wind priors. As worded, the headline overstates what the radio data alone establish.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a systematic MCMC analysis of 32 radio-selected supernovae with clear light-curve peaks, fitting a standard analytical radio SN model to 27 objects. The fitted parameters are the CSM density scale A~, electron acceleration efficiency eps_e, magnetic field amplification efficiency eps_B, electron spectral index p, CSM density slope s, and outer ejecta density slope n. The authors report that stripped-envelope SN (SESN) progenitors have mass-loss rates about an order of magnitude higher than SN II progenitors, that shock acceleration efficiencies are below 10^-2 with equipartition allowed, and that the outer ejecta density slope tends to be shallow (n ~ 5). When early-phase (first 30 days) data are included, the fits prefer extremely high eps_B, which the authors interpret as evidence for a dense CSM component near the progenitor that is not captured by the single-component model.","tokens_in":33156,"tokens_out":4746,"duration_ms":49334,"significance":"If the central mass-loss claim held, this would be an important systematic constraint on massive-star mass loss and on the progenitors of stripped-envelope supernovae. The paper's strengths include the use of a uniform model framework across a relatively large sample, a transparent MCMC setup, and a robustness test with fsh = 4 that does not qualitatively change the fitted microphysical parameters. The paper is also candid about several limitations. However, the headline mass-loss difference is not actually driven by the fitted CSM density scale but by the assumed wind velocities, and the statistical significance of the group-level difference is not established. The dense-CSM interpretation also has a partially circular structure. These issues affect the paper's main conclusions and require revision before the claims can be taken at face value.","major_comments":[{"comment":"The abstract and Section 8 claim that SESN progenitors have mass-loss rates an order of magnitude above those of SNe II, but the fitted CSM density scale A~ is nearly the same for the two groups (median log A~ = 1.75 +/- 1.53 for SNe II vs 1.47 +/- 1.08 for SESNe, Table 5). Since Mdot = 4*pi*v_w * 5e11 g/cm * A~ (Eq. 4), the factor of ~100 in Mdot comes almost entirely from the assumed wind velocities (10 km/s for SNe II vs 1000 km/s for SESNe), not from the radio data. Section 5.1 acknowledges the uncertainty in the CSM velocity, but this limitation is not propagated into the abstract or summary claims, and no alternative wind-velocity priors are tested. As written, the headline overstates what the radio data alone establish.","section":"Section 5.1, Eq. (4), Table 5"},{"comment":"With only four SNe II in the successful-fit sample (SN 1987A, 2004dj, 2012aw, 2016X) and 1-sigma scatter of about 1.5 dex in log A~, the difference in A~ between SNe II and SESNe is not statistically significant. The manuscript should include a formal significance assessment (e.g., posterior overlap or a two-sample test) or explicitly state that the current sample cannot distinguish the group-level CSM density scales. As written, the 'order of magnitude' statement in the abstract is not supported by the fitted density scales alone.","section":"Table 5, Section 5.1"},{"comment":"The dense-CSM interpretation is introduced after excluding the first 30 days of data because the model cannot treat inhomogeneous CSM, and the same exclusion is then presented as evidence that a dense inner CSM is the missing setup. This has a circular structure. The manuscript does mention alternative explanations (genuinely large eps_B, time-dependent microphysics), but the abstract wording 'we identify the missing setup as a dense CSM' is stronger than the evidence. A direct test with a two-component or hydrodynamical CSM model, or at least an explicit statement that the dense-CSM interpretation is one of several mutually consistent possibilities, is needed to avoid circular reasoning.","section":"Sections 3.3 and 5.2"},{"comment":"The exclusion of early-phase data does not remove the dynamical influence of a dense inner CSM on the shock trajectory at later times. The manuscript itself states that the shock evolution is first influenced by the inner CSM and that subsequent hydrodynamical properties would be affected, and that this effect is not included in the model. If the inner CSM modifies the shock before day 30, the fitted late-time parameters, including A~ and s, could be biased. The paper should quantify this effect (for example, by comparing with the two-component calculations of Matsuoka et al. 2025) or explicitly state that the mass-loss and density-slope results are conditional on the unperturbed single-component CSM assumption.","section":"Section 3.3"}],"minor_comments":[{"comment":"The units in the mass-loss rate ranges are incomplete: '10^-7 <= Mdot < 10^-4 M_sun' and 'Mdot ~ 10^-3 M_sun' should read M_sun yr^-1 throughout.","section":"Section 8, first bullet"},{"comment":"The final column header 'References of chi^2_red < 5' appears garbled; it should be relabeled (e.g., 'Good fit (chi^2_red < 5)') and the checkmarks defined clearly in the caption.","section":"Table 1"},{"comment":"The imposed lower limit sigma_obs,i = 0.1 F_obs,i is applied uniformly, but its effect on the derived posterior widths and on the reduced chi-square threshold is not discussed. A brief justification or a sensitivity test would clarify how much of the quoted credible intervals depends on this assumption.","section":"Section 3.2"},{"comment":"The phrase 'A~ only in the inner CSM' refers to a localized enhancement that is not a fitted parameter of the model; the text should clarify that this is an inferred physical interpretation, not a quantity directly constrained by the MCMC.","section":"Section 5.2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a careful look. This is the first time anyone has run the same radio SN model, with the same priors and fitting procedure, across a decent sample of objects with clear radio peaks. That alone is a real step forward: earlier work was a pile of individual fits with different fixed parameters and different assumptions, so cross-object comparisons were shaky. The authors fit 27 of 32 objects, show corner plots, test the fsh = 4 systematic, impose a sensible flux-calibration error floor, and are transparent about where their posteriors hug the prior boundary (n ~ 5, especially). The microphysical conclusions — epsilon_e and epsilon_B around 1e-3 to 1e-2, p > 2.5, equipartition not excluded — are modest but honestly presented.\n\nThe soft spots are real but not fatal. The big one is flagged in the stress-test note, and I think it lands. The fitted quantity is log A~, the CSM density scale, and the median values for SNe II (1.75) and SESNe (1.47) are basically the same, even slightly higher for SNe II. The order-of-magnitude gap in Mdot comes almost entirely from multiplying by v_w = 10 km/s for SNe II and 1000 km/s for SESNe. Equation (4) is just Mdot = 4π v_w × 5e11 × A~. So the abstract's claim that SESN progenitors have Mdot ~ 1e-3 vs < 1e-4 for SNe II is not a statement the radio light curves themselves establish; it is a statement about the assumed wind velocities. Section 5.1 does acknowledge the uncertainty, but the abstract and summary don't, and there is no test with alternative wind priors. With only 3-4 SNe II and ~1.5 dex scatter on A~, the statistical case is weak.\n\nSecond, the dense-CSM interpretation is indirect. The authors exclude the first 30 days because the single-power-law model cannot handle an inhomogeneous inner CSM; then they show that including those points drives epsilon_B to ~1; then they argue the missing ingredient is a dense confined CSM. That is a reasonable narrative, consistent with flash spectroscopy and SN 2020oi modeling, but it is not a detection. They also note in Section 3.3 that the inner CSM can change the shock trajectory and thus bias the late-time fitted parameters, but they do not quantify this. Third, no code or data release, which matters for a paper whose main product is a set of fitted distributions. The n ~ 5 pile-up and the s < 2 tendency are flagged by the authors themselves, so I won't count those heavily.\n\nBottom line: this is a solid subfield contribution that deserves a serious referee, but the mass-loss-rate comparison needs to be reframed as conditional on wind velocities, or tested against a range of wind assumptions. I would send it to review and ask for that revision before publication.","headline":"A genuinely useful first systematic MCMC fit of 32 radio SNe, but the headline mass-loss gap between SNe II and SESNe is mostly inherited from the assumed wind velocities, not from the radio data.","tokens_in":33727,"tokens_out":1659,"would_cite":true,"duration_ms":19783,"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":"A uniform MCMC analysis of 27 radio supernovae finds that stripped-envelope progenitors lose mass at roughly $10^{-3}\\,M_\\odot\\,{\\rm yr}^{-1}$, an order of magnitude more than Type II progenitors, and attributes several anomalous fit…","keywords":["radio supernovae","circumstellar medium","mass-loss history","shock acceleration","magnetic field amplification","Markov chain Monte Carlo","stripped-envelope supernovae","confined circumstellar material"],"falsifier":"Observe a newly discovered stripped-envelope supernova with radio monitoring starting within days of explosion; the confined-CSM claim predicts that the early light curve will require a density enhancement near $10^{15}$ cm that is much higher than the extrapolated late-time wind and that a two-zone model will fit the data with $\\epsilon_B < 0.01$. If a two-zone fit still demands $\\epsilon_B\\sim 1$, or if no early absorption excess appears, the artifact interpretation fails.","tokens_in":32670,"feed_emoji":"📡","tokens_out":6115,"duration_ms":62314,"temperature":0.7,"pith_summary":"This paper argues that a uniform six-parameter radio-supernova model, fitted with MCMC to 27 well-observed radio supernovae, reveals two statistical facts. First, stripped-envelope supernova progenitors lose mass at about $10^{-3}\\,M_\\odot\\,{\\rm yr}^{-1}$ in the final years before explosion, an order of magnitude faster than Type II progenitors at $<10^{-4}\\,M_\\odot\\,{\\rm yr}^{-1}$. Second, several puzzling fit results—near-unity magnetic-field amplification efficiency and an unusually shallow outer ejecta slope—disappear if a dense circumstellar shell within roughly $10^{15}$ cm of the progenitor is added to the picture. The authors conclude that such a confined dense shell is common among radio-bright supernovae, and that the standard single-power-law circumstellar medium is incomplete.","feed_headline":"A dense shell surrounds many exploding stars","feed_subtitle":"MCMC fits to 27 radio supernovae show that stripped-envelope progenitors lost mass ten times faster than SNe II.","key_machinery":"The workhorse is the standard synchrotron self-absorbed radio supernova model: a homologously expanding ejecta with outer density slope $n$ drives a thin-shell shock into a power-law CSM with density $\\rho = D r^{-s}$, and the radio luminosity is computed from the relativistic electron spectrum and amplified magnetic field, with parameters $\\boldsymbol{\\theta}=(\\log \\tilde{A}_*, \\log\\epsilon_e, \\log\\epsilon_B, p, s, n)$. The MCMC sampler maps the posterior, and the key diagnostic is a degeneracy between the CSM density scale $\\tilde{A}_*$ and the magnetic efficiency $\\epsilon_B$: both parameters dim optically thick emission and brighten optically thin emission in similar ways, so a model lacking a dense inner shell misattributes the early-time excess to $\\epsilon_B\\sim 1$. Excluding the first 30 days restores moderate $\\epsilon_B$, and the paper reads the resulting discrepancy as evidence for a dense confined CSM at $r\\lesssim 10^{15}$ cm.","core_discovery":"The paper reports that, under one analytical radio emission model applied across 32 supernovae, the inferred mass-loss rate of stripped-envelope progenitors clusters around $10^{-3}\\,M_\\odot\\,{\\rm yr}^{-1}$, an order of magnitude higher than for Type II supernovae, while the electron acceleration efficiency and magnetic field amplification efficiency are both below $10^{-2}$ and their equipartition is not ruled out. When the same fits include the first 30 days of data, many objects demand $\\epsilon_B \\sim 1$, a value far above what particle-in-cell shock simulations produce. The paper interprets this as a model artifact: the early radio data trace a dense circumstellar medium in the immediate vicinity of the progenitor that is not smoothly connected to the outer wind, so the standard model compensates by inflating $\\epsilon_B$ and lowering the CSM density scale. It also finds outer ejecta slopes as small as $n\\sim 5$, below classical shock-breakout predictions, and links this to the same confined CSM flattening the outer ejecta structure.","pith_inferences":["A direct test of the confined-CSM interpretation: early multi-band radio follow-up within days of explosion should reveal free-free absorption or a density excess near $10^{15}$ cm in objects like those studied here, and a two-zone model with such a shell should fit the early data without requiring $\\epsilon_B\\sim 1$.","The $\\tilde{A}_*$-$\\epsilon_B$ degeneracy implies that published single-zone radio fits that fix $\\epsilon_e=\\epsilon_B=0.1$ carry a hidden systematic bias in inferred CSM density, so future work should fit density and microphysics jointly rather than fixing either.","The confined shell proposed here may be the same phenomenon seen in flash-spectroscopy supernovae and in some stripped-envelope events, suggesting a common pre-explosion mass-ejection mechanism across apparently different transient classes.","A hydrodynamic simulation of a shock crossing a dense inner shell would provide a decisive check: if it reproduces the early rise with $\\epsilon_B<0.01$, the artifact interpretation is confirmed; if not, the large $\\epsilon_B$ from early-included fits would reflect genuine shock microphysics."],"forward_implications":["Mass-loss histories of stripped-envelope progenitors should be revised upward by roughly an order of magnitude relative to equipartition-based estimates, with consequences for final stellar evolution and pre-supernova activity.","Radio-bright supernovae should commonly show signatures of dense confined CSM, including flash-spectroscopy features and early-time radio absorption.","Shock acceleration in non-relativistic supernova shocks need not be highly efficient: efficiencies below $10^{-2}$ with a viable near-equipartition ratio between electrons and magnetic field remain consistent with the data.","The outer ejecta of core-collapse supernovae interacting with a confined shell can be shallower than classical predictions, which affects how ejecta mass and kinetic energy are inferred from radio light curves.","Single-power-law fits that include early-time radio data should be interpreted with caution, because the near-unity $\\epsilon_B$ they produce is a warning sign of missing density structure rather than a reliable microphysics measurement."],"supporting_citations":[{"why":"Provides the standard review and formalism for radio emission from supernova-circumstellar medium interaction on which the model is built.","marker":"Chevalier & Fransson 2017"},{"why":"Supplies the thin-shell self-similar shock solution and synchrotron emission framework used to compute the light curves.","marker":"Chevalier 1982a"},{"why":"Derives the electron cooling treatment and synchrotron self-absorption spectrum that sets the model's spectral shape.","marker":"Fransson & Björnsson 1998"},{"why":"Supplies the comprehensive radio supernova catalog from which the paper selects its 32-object golden sample and provides the comparison for mass-loss estimates.","marker":"Bietenholz et al. 2021"},{"why":"Provides the emcee MCMC sampler used for all parameter estimation and posterior sampling.","marker":"Foreman-Mackey et al. 2013"},{"why":"Presents the SN 2020oi case study showing an inhomogeneous CSM near $10^{15}$ cm, the direct observational precedent for the missing dense inner shell.","marker":"Maeda et al. 2021"},{"why":"Shows that an inhomogeneous inner CSM affects the subsequent shock evolution, supporting the paper's caveat that early inner structure can bias late-time fitted parameters.","marker":"Matsuoka et al. 2025"},{"why":"Provides the classical theoretical prediction of the steep outer ejecta slope that the paper's $n\\sim 5$ result is compared against.","marker":"Matzner & McKee 1999"},{"why":"Provides representative ejecta masses and kinetic energies adopted for stripped-envelope supernovae lacking individual estimates.","marker":"Lyman et al. 2016"},{"why":"Particle-in-cell simulations that bound magnetic field amplification efficiency at values far below $\\epsilon_B\\sim 1$, making the early-included fit result physically implausible.","marker":"Sironi & Spitkovsky 2011"}],"fun_headline_variants":["Dense shells may surround many radio supernovae","Stripped-envelope supernovae lose mass ten times faster","MCMC on 27 supernovae uncovers dense circumstellar shells","Early radio data suggest dense shells around exploding stars","Supernova radio fits trace a dense inner shell"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that any dense shell near the star affects only the first thirty days, so that the later radio light curve can still be modeled with a single smooth wind and a shock trajectory unaffected by the inner structure.","fun_headline_variants_meta":{"raw":{"variants":["Dense shells may surround many radio supernovae","Stripped-envelope supernovae lose mass ten times faster","MCMC on 27 supernovae uncovers dense circumstellar shells","Early radio data suggest dense shells around exploding stars","Supernova radio fits trace a dense inner shell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1824,"prompt_tokens":1044,"completion_tokens":780,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":699}},"tokens_in":660,"tokens_out":780,"duration_ms":8023,"temperature":1.0,"reasoning_tokens":699,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:38:03.474414+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a newly discovered stripped-envelope supernova with radio monitoring starting within days of explosion; the confined-CSM claim predicts that the early light curve will require a density enhancement near $10^{15}$ cm that is much higher than the extrapolated late-time wind and that a two-zone model will fit the data with $\\epsilon_B < 0.01$. If a two-zone fit still demands $\\epsilon_B\\sim 1$, or if no early absorption excess appears, the artifact interpretation fails.","supporting_citations":[],"review_version":1}