{"id":"97dd17b3-3ae0-4039-ae74-2be0b76e5ac5","arxiv_id":"2502.01740","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"SN 2001ig has re-brightened in the radio more than 20 years after the explosion, likely because its shock hit a dense circumstellar shell.","lead":"Radio observations show the Type IIb supernova SN 2001ig is glowing much more brightly two decades after exploding, reversing the expected fade. This re-brightening likely means the blast wave has hit a dense shell of gas shed by the star in its final years.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Radio re-brightening is robust, but the 'denser CSM shell' interpretation depends on the untested assumption that epsilon_B and epsilon_e are constant; time-evolving microphysical fractions could mimic the flux rise without any density enhancement.","rationale":"The paper's observational discovery is compelling: the late-time ATCA and ASKAP detections are many epochs, multiple frequencies, and the re-brightening is large (two orders of magnitude) relative to a well-sampled early power-law decline with a stable optically thin spectral index. I agree with the reader's conditional assessment. The weakest point is not the photometry but the physical interpretation. The reader identified the Chevalier SSA scalings and the assumed microphysical fractions as load-bearing. I sharpen that concern: the decisive assumption is not just the values of epsilon_B and epsilon_e but their constancy in time. The observed flux in the optically thin regime is determined by the shock radius, magnetic field (proportional to sqrt(epsilon_B rho v_s^2)), and non-thermal electron energy density. If epsilon_B or epsilon_e increase with time, the radio flux can rise without any change in the CSM density profile. The paper explicitly assumes these fractions do not evolve (Section 4.2), and the 3-10 year data gap prevents the light curve from ruling out a gradual microphysical evolution. A time-varying microphysics model with a standard wind could in principle reproduce the re-brightening, and the paper does not test this. The proposed test—allowing epsilon_B(t) and epsilon_e(t) as power laws with no density break and comparing model evidence—directly addresses this degeneracy. Unless such a model is strongly disfavoured, the claim that the radio re-brightening is caused by a denser CSM shell remains plausible but unproven. The paper's transparency about degeneracies and lower limits supports keeping the verdict CONDITIONAL rather than moving to REJECT: the observations and the qualitative inference of an enhanced interaction region are likely correct, but the quantitative parameters (R_brk ≈ 0.1 pc, mass-loss rate, density enhancement) should be treated as model-dependent. Thus no change to the reader's verdict.","tokens_in":22197,"tokens_out":11442,"duration_ms":111561,"concrete_test":"Re-run the MCMC fitting of Section 4.2 on the full 2001-2024 light curves with a pure wind CSM (rho ∝ R^-2, no break) and allow the microphysical fractions to vary as power laws in time, epsilon_B(t) = epsilon_B0 (t/t0)^q and epsilon_e(t) = epsilon_e0 (t/t0)^r, with all other assumptions (Mej, n, f) unchanged. Compare the Bayesian evidence (e.g., Δln Z) of this model against the two-zone density-break model. If the time-varying microphysics model fits comparably well and yields reasonable q and r values, the density-shell interpretation is degenerate and the paper's central explanation is not uniquely supported; if the density-break model is strongly preferred, the concern is laid to rest.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The re-brightening itself is on solid observational ground: multiple ATCA and ASKAP epochs show fluxes roughly two orders of magnitude above the early power-law extrapolation, with a stable spectral index alpha ≈ -1 (Section 3). The load-bearing point is the interpretation that the shock has entered a denser CSM region. The model maps observed fluxes to shock radius, magnetic field, and CSM density using fixed microphysical fractions epsilon_B = epsilon_e = 0.1 (Section 4.1, Eqs. 1-3), and Section 4.2 states that 'these fractions do not evolve with time'. If epsilon_B or epsilon_e vary with time (e.g., due to changing shock velocity, obliquity, or magnetic-field amplification efficiency), the synchrotron flux can rise while the CSM density follows a normal R^-2 wind, so the inferred density enhancement and the break radius R_brk ≈ 3e17 cm would not be required. The paper's caveat that lower epsilon_B implies higher densities addresses only the normalization degeneracy, not this time-variability degeneracy. The absence of radio observations between roughly 3 and 10 years (the gap highlighted in Figure 4) means a smoothly growing efficiency could be absorbed into the model without being noticed. Therefore, the central physical claim rests on the time-independence of the microphysical parameters, which is the least secure link in the argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports new ATCA (2.1, 5.5, 9 GHz) and ASKAP (0.89–1.37 GHz) observations of SN 2001ig at ages 12–22 years, showing that its radio flux has re-brightened by about two orders of magnitude relative to the power-law decline measured in 2001–2003. The spectral index has remained around α ≈ –1 during the available late-time epochs. Using the Chevalier synchrotron self-absorption formalism with a two-zone CSM density profile and a thin-shell dynamical model, the authors infer that the shock encountered a density enhancement at Rbrk ≈ 0.1 pc, with a shallower density profile (α_over = 0.75) beyond the break, and estimate a pre-explosion mass-loss rate Mdot/v_w ≈ 1×10^-7 Msun/yr/(km/s). They interpret the re-brightening as evidence of interaction with a dense shell and compare SN 2001ig with other late-time re-brightening SNe.","tokens_in":22437,"tokens_out":8555,"duration_ms":76494,"significance":"The observational result is significant: SN 2001ig is one of the closest and best-monitored Type IIb SNe, and this is the first report of a late-time radio re-brightening in this object. The multi-frequency detections at ATCA and ASKAP, the consistency of the spectral index, and the high signal-to-noise of the 2024 detections make the re-brightening itself a solid empirical result. The paper also makes good use of archival data and is transparent about degeneracies in the model, particularly the gap in coverage and the normalization degeneracy with ε_B. However, the physical interpretation—a denser CSM shell—rests on the assumed time-independence of the microphysical parameters, which is not tested. If that assumption fails, the inferred density enhancement and break radius are not unique. The comparison with other re-brightening SNe is useful but also inherits the same model dependence.","major_comments":[{"comment":"The modelling assumes that ε_B and ε_e are constant in time, and this assumption is load-bearing for the claim that the re-brightening is caused by a denser CSM region. A gradual increase of ε_B (or of the electron acceleration efficiency) over 20 years, with no change in the R^-2 wind density, would raise the synchrotron flux in the same gradual way, and the 3–10 yr observational gap (Fig. 4) means such a trend could go undetected. The paper's caveat that lower ε_B implies higher CSM densities (Section 4.2) addresses only the absolute normalization, not this time-dependence. Please add a physical justification for time-constant microphysics, an explicit test with time-varying ε_B, or a substantial softening of the physical conclusion in the abstract and Section 5.","section":"Section 4.2, Eqs. (5)–(9)"},{"comment":"The text in Section 4.2 acknowledges 'strong degeneracies between Rbrk and the normalizations of ρ0,wind and ρ0,over', but the abstract and Section 5 quote a single value, Rbrk ≈ 0.1 pc, without reporting the posterior distribution. The corner plots (Figs. 7–8) are not summarized numerically. Please report the median and credible intervals for Rbrk, ρ0,over/ρ0,wind, and α_over for both the 1 M⊙ and 4 M⊙ cases, and present the density-break radius as a range in the abstract and conclusions.","section":"Section 4.2 and Fig. 4"},{"comment":"The statement that 'the inferred trend of a flattening density index at late times is a solid result regardless of the choices of Mej and ε_B' is too strong. The density-index flattening is derived from the same constant-ε_B mapping as the other quantities, so it inherits the time-variability degeneracy described above; a time-dependent ε_B could flatten the inferred density profile without any actual change in the CSM slope. Please qualify this claim or support it with a test.","section":"Section 4.2, final paragraph"}],"minor_comments":[{"comment":"The phrase 'in 2004 April' appears twice where the most recent 5.5 GHz measurement is quoted; the flux density of 3.3 mJy corresponds to the 2024 April 24 observation, so the year should be 2024 in both places.","section":"Section 3, item (i); Section 5"},{"comment":"The last row of Table 1 contains a stray colon in '2024-04-24:17:14:00'; it should read '2024-04-24 17:14:00' for consistency with the other rows.","section":"Table 1"},{"comment":"The notation 'B_p^2/1 G' is ambiguous; the numerical evaluation corresponds to (B_p / 1 G)^2, so the equation should be typeset accordingly to avoid confusion about the units.","section":"Equation (3)"},{"comment":"Several author names and words contain stray spaces from LaTeX source rendering, e.g., 'W olf-Rayet', 'Y oon', 'V oronkov'; these should be corrected in the final version.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The observational discovery is solid and well-presented, and the authors are appropriately cautious about the model degeneracies they do discuss. The main gap is that the time-independence of the microphysical parameters is the least secure link in the chain from data to the physical conclusion, yet it is not addressed beyond a normalization caveat. A revision that either defends the constancy of ε_B, tests a time-varying model, or substantially tempers the abstract and Section 5 claims would make the paper suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, the short version: this paper is worth a serious referee. The radio re-brightening of SN 2001ig is a genuine new result—two orders of magnitude above the early power-law decline, seen in ATCA and ASKAP data at multiple epochs and frequencies, with the spectral index stuck near -1. That is a robust observational statement and it makes the object a nearby anchor for late-time CSM interaction in compact Type IIb SNe. The paper also does the comparative work properly, placing 2001ig against SN 1987A, 2014C, 2003bg, and the VLASS/VAST samples, and it is transparent about its modeling assumptions. The mass-loss rate and break radius are derived, not measured, and the paper says so. Credit where due: the data reduction looks careful, including a weather-decorrelation correction on the 2024 February point, and the corner plots show real MCMC output rather than hand-waving. The modeling follows Chevalier's SSA formalism with fixed epsilon_B = epsilon_e = 0.1, f = 0.5, thin shell, and a two-zone CSM. The soft spot, as the stress-test note says, is the assumption that the microphysical fractions don't evolve with time. If epsilon_B grows over the decades, some part of the flux rise could be mimicked without any true density enhancement. The 3-10 year gap in coverage (Figure 4) means a slow increase in efficiency could hide in the model. I don't think this kills the paper's central qualitative claim—the flux rise is large and the spectral index stays put, which is more naturally read as a genuine CSM over-density—but it does mean the quoted break radius and mass-loss-rate error bars are optimistic. The paper actually acknowledges the Rbrk degeneracy up front, which is good, though it stops short of folding the epsilon_B time-variability into the systematic budget. The density profile inferred at late times is a fitted result, not an independent measurement; readers should treat 'denser shell at 0.1 pc' as a plausible model-dependent interpretation, not an observed fact. That is exactly the kind of distinction a good referee can sharpen. My recommendation: send it out. A serious referee will push on the time-constant-epsilon assumption, the gap-driven degeneracy, and the lack of formal error bars on Mdot and Rbrk, but the observational discovery is solid and the paper is already honest about many of its own limitations. It deserves space in the literature and will get cited. I'd bring it to reading group to talk about how much of the late-time rebrightening population is real CSM structure versus evolving microphysics.","headline":"Solid new observation of late-time radio re-brightening in SN 2001ig; the physical interpretation is plausible but rests on untested time-constant microphysics, so the paper is worth refereeing with the model claims treated as provisional.","tokens_in":23093,"tokens_out":1981,"would_cite":true,"duration_ms":18797,"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":"SN 2001ig's radio emission is re-brightening two decades after explosion, to flux densities about a hundred times above the 2001–2002 power-law decline, which the authors attribute to the shock reaching a dense shell at roughly $0.1$ pc.","keywords":["supernovae: general","ISM: supernova remnants","stars: massive","Type IIb supernovae","radio re-brightening","circumstellar medium","SN 2001ig"],"falsifier":"Very long baseline interferometry at 5–9 GHz can settle the central claim: at 10.8 Mpc, a shock at $0.1$ pc would subtend about 2 milliarcseconds, so resolving the remnant at that size would confirm the radius, while a measured radius a factor of two away from $0.1$ pc would falsify the model's mapping between flux and radius.","tokens_in":21921,"feed_emoji":"📡","tokens_out":13922,"duration_ms":111297,"temperature":0.7,"pith_summary":"This paper reports that the Type IIb supernova SN 2001ig, in the galaxy NGC 7424, is re-brightening in the radio more than two decades after it exploded: its measured flux density is now about two orders of magnitude higher than the power-law decline seen in 2001–2002 would predict. The authors interpret this as the blast wave running into a denser circumstellar shell, and they model the encounter as a break in the density profile at a radius of roughly $0.1$ pc, where the surrounding medium stops following the standard $R^{-2}$ wind and instead declines as $R^{-0.75}$. From the same modelling they derive a pre-explosion mass-loss rate of $\\dot{M}/v_w \\sim 10^{-7}\\,M_\\odot\\,\\text{yr}^{-1}\\,\\text{km}^{-1}\\,\\text{s}$. If correct, the result turns a single nearby supernova into a probe of the final centuries of mass loss from a stripped, compact progenitor star, and it sharpens the distinction between compact and extended Type IIb events.","feed_headline":"SN 2001ig's radio glow is now 100 times brighter than expected","feed_subtitle":"New ATCA and ASKAP data suggest the blast wave hit a dense shell at about 0.1 pc.","key_machinery":"The quantitative argument rests on the Chevalier (1998) synchrotron self-absorption formalism, which turns a measured peak flux density into a shock radius, magnetic field, and circumstellar density. With energy equipartition fractions $\\epsilon_B = \\epsilon_e = 0.1$ and a volume filling factor $f = 0.5$, the paper uses the peak flux at 18.8 GHz to fix the early shock radius and velocity, then integrates the thin-shell momentum equation with ejecta and CSM density profiles to evolve the shock. The resulting radio spectra are fit to all epochs with an MCMC that returns the break radius $R_{\\rm brk}$ and the post-break density index $\\alpha_{\\rm over}$. This machinery is what allows measured flux densities to be translated into a physical picture of a dense shell at $0.1$ pc.","core_discovery":"The paper's central claim is that the late-time radio emission of SN 2001ig is genuinely re-brightening rather than continuing the power-law decline established in 2001–2002, with today's flux density two orders of magnitude above the extrapolation. The re-brightening started no later than an age of about 4000 days, may have begun as early as 1000 days, and is still ongoing; the spectral index has stayed at $\\alpha \\approx -1$, which the authors take as evidence that the same optically thin synchrotron mechanism is at work. To explain it, they construct a two-zone circumstellar medium in which the shock first moves through a standard wind ($\\rho \\propto R^{-2}$) and then, at $R_{\\rm brk} \\approx 3\\times10^{17}$ cm, encounters a region with a shallower density slope ($\\rho \\propto R^{-0.75}$), producing an overdensity of about an order of magnitude at the current epoch relative to the extrapolated wind. Both a $1\\,M_\\odot$ and a $4\\,M_\\odot$ ejecta-mass model reproduce the observations, and the inferred immediate pre-explosion mass-loss rate is $\\dot{M}/v_w \\simeq 1.3\\times10^{-7}\\,M_\\odot\\,\\text{yr}^{-1}\\,\\text{km}^{-1}\\,\\text{s}$. The authors explicitly note that the density normalisation depends on the assumed microphysical energy fractions and should be read as a lower limit.","pith_inferences":["A lower $\\epsilon_B$ than the assumed 0.1, as favoured by some recent shock-acceleration models, would push the inferred CSM densities upward by about an order of magnitude; that would make the shell mass large enough that a binary common-envelope or eruptive mass-loss origin becomes more attractive than the minimal wind-compression picture the paper adopts.","The early-time roughly 150-day radio modulation and the late-time shell may share a single origin: if the early modulation is a spiral density pattern from a wind-wind collision in a binary system, the same pattern swept outward could assemble the shell at $0.1$ pc, and a single spiral-density model fit to both early and late epochs would test this directly.","A systematic radio survey of Type IIb SNe at ages of 10–30 years could map the apparent dichotomy between re-brightening and flux cutoff; if the two behaviours turn out to be a continuum rather than two classes, the progenitor-radius interpretation would need revision.","X-ray follow-up of SN 2001ig over the next few years could corroborate the shell-interaction picture independently of the radio: a shock running into denser material should produce enhanced thermal and inverse-Compton X-ray emission."],"forward_implications":["If the re-brightening is caused by a dense shell at about $0.1$ pc, continued radio monitoring should show the flux eventually peaking and then declining as the shock exits the overdense region, with the turnover timescale set by the shell thickness.","The inferred mass-loss rate of $\\dot{M}/v_w \\sim 10^{-7}\\,M_\\odot\\,\\text{yr}^{-1}\\,\\text{km}^{-1}\\,\\text{s}$ connects the explosion to a specific mass-loss episode centuries before collapse, placing a direct constraint on late-stage stellar evolution models.","SN 2001ig becomes a clean example of the compact-progenitor path for Type IIb SNe, standing in contrast to extended-progenitor events like SN 1993J whose radio flux drops sharply at late times.","If the dense shell is hydrogen-rich, optical recombination lines such as H-alpha may reappear or brighten in the coming years, offering an independent test of the shell's composition.","The post-break density index $\\alpha_{\\rm over} \\approx 0.75$ implies the shock will keep producing significant radio emission for an extended period, making SN 2001ig a promising target for long-term radio monitoring."],"supporting_citations":[{"why":"Supplies the synchrotron self-absorption equations that convert peak flux density into shock radius and magnetic field, the basis of both the mass-loss estimate and the density profile fits.","marker":"Chevalier 1998"},{"why":"Provides the 2001–2003 radio light curve and the ~150-day modulation; its power-law decline is the baseline from which the re-brightening is measured.","marker":"Ryder et al. 2004"},{"why":"Defines the canonical radio-SN model of rise and power-law decline used to extrapolate the early-time light curve to late times.","marker":"Kurt W. Weiler et al. 2002"},{"why":"Gives the thin-shell momentum equation and CSM interaction framework used to evolve the shock radius and velocity.","marker":"Chevalier and Fransson 2017"},{"why":"Provides the SN 2014C case of a stripped SN hitting a dense shell, the main comparative template and physical scenario for the shell interaction.","marker":"Margutti et al. 2017"},{"why":"Supplies the ASKAP sample of late-time radio re-brightenings and the method used to estimate mass-loss rates from peak flux.","marker":"Rose et al. 2024"},{"why":"Provides the VLASS sample of late-time re-brightening SNe and the critique of the off-axis jet interpretation used to argue against that scenario here.","marker":"Stroh et al. 2021"},{"why":"Establishes the compact-progenitor interpretation and the He II detection that hinted at ongoing CSM interaction, motivating the late-time study.","marker":"Ryder et al. 2018"}],"fun_headline_variants":["SN 2001ig radio re-brightens 100x above expected","Supernova 2001ig flares again in radio after 20 years","Radio re-brightening of SN 2001ig points to dense shell","SN 2001ig's radio glow jumps two orders of magnitude","Late-time radio re-brightening of SN 2001ig from dense shell"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fractions of post-shock energy going into magnetic fields and relativistic electrons are both 10 percent and that the emitting region fills half the sphere; if those assumed fractions are wrong, the inferred densities, shell radius, and mass-loss rate all shift, and the paper itself notes that lower $\\epsilon_B$ would imply higher densities.","fun_headline_variants_meta":{"raw":{"variants":["SN 2001ig radio re-brightens 100x above expected","Supernova 2001ig flares again in radio after 20 years","Radio re-brightening of SN 2001ig points to dense shell","SN 2001ig's radio glow jumps two orders of magnitude","Late-time radio re-brightening of SN 2001ig from dense shell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000911,"raw_usage":{"total_tokens":3989,"prompt_tokens":1095,"completion_tokens":2894,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":711,"completion_tokens_details":{"reasoning_tokens":2795}},"tokens_in":711,"tokens_out":2894,"duration_ms":19661,"temperature":1.0,"reasoning_tokens":2795,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T14:38:22.514960+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Very long baseline interferometry at 5–9 GHz can settle the central claim: at 10.8 Mpc, a shock at $0.1$ pc would subtend about 2 milliarcseconds, so resolving the remnant at that size would confirm the radius, while a measured radius a factor of two away from $0.1$ pc would falsify the model's mapping between flux and radius.","supporting_citations":[],"review_version":1}