Pith. sign in

REVIEW 3 major objections 4 minor 6 references

The radio re-brightening of the Type IIb SN 2001ig

T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2502.01740 v1 pith:5R4YUHZA submitted 2025-02-03 astro-ph.HE astro-ph.GA

Roberto Soria , Thomas D. Russell , Eli Wiston , Siying Cheng , Raffaella Margutti , Kovi Rose , Stuart Ryder , Giacomo Terreran This is my paper · ORCID
classification astro-ph.HEastro-ph.GA
keywords supernovae:generalISM:supernovaremnantsstars:massiveTypeIIbsupernovaeradiore-brighteningcircumstellarmediumSN2001ig
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

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.

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 (3)
  1. [Section 4.2, Eqs. (5)–(9)] 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.
  2. [Section 4.2 and Fig. 4] 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.
  3. [Section 4.2, final paragraph] 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.
minor comments (4)
  1. [Section 3, item (i); Section 5] 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.
  2. [Table 1] 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.
  3. [Equation (3)] 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.
  4. [Throughout] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

Radio re-brightening is an independent observational result; the CSM-density interpretation is a model fit under stated assumptions, not a circular derivation.

full rationale

The central claim (Section 3) is an observational comparison: 'there is strong re-brightening compared with the extrapolated power-law decline of 2001–2002. Today's radio flux density is two orders of magnitude higher than we expect assuming a simple power-law model.' This comparison is based on the early ATCA/VLA light curve (Ryder et al. 2004) and does not use the late-time model to define the baseline. The physical parameters (Rbrk, density enhancement, alpha_over, Mdot/v_w) are obtained by inverting the Chevalier (1998) synchrotron self-absorption scalings with stated assumptions (epsilon_B=epsilon_e=0.1, f=0.5, time-independent fractions) and are explicitly fit parameters: 'The fit parameters are: Rbrk, ρ0,wind, ρ0,over, and αover.' Inverting observed fluxes into a CSM density profile is inverse modeling, not a self-definitional reduction; no quantity is defined in terms of the quantity it is claimed to predict, and no fitted parameter is renamed as a prediction. The paper is transparent about the main degeneracies: 'there are still strong degeneracies between Rbrk and the normalizations of ρ0,wind and ρ0,over', and it states that the density models 'can best be interpreted as lower limits on the CSM density profile.' The claim that the flattening density index is 'a solid result regardless of the choices of Mej and ϵB' is overstrong because only the normalization of ϵB was varied, not its time evolution; however, this is a model-robustness caveat, not circularity. Citations to Ryder et al. (2004), Rose et al. (2024), and Chevalier (1998) are data sources or standard framework references, not a self-citation chain that forces the result. Overall, the re-brightening is self-contained against external benchmarks, and the derived CSM properties are honestly presented as model fits rather than independent predictions.

Assumptions & free parameters 11 free parameters · 8 assumptions · 0 invented entities

The central re-brightening measurement is model-independent, but every quantitative claim (density, mass-loss rate, break radius, density enhancement) is obtained through the assumed microphysical constants and the two-zone CSM model. The main driver of uncertainty is the assumed equipartition fractions, which set the absolute density scale.

free parameters (11)
  • epsilon_B = 0.1
    Assumed fraction of post-shock energy in the magnetic field; directly scales the inferred CSM density and mass-loss rate. Lower values would imply higher densities.
  • epsilon_e = 0.1
    Assumed fraction of post-shock energy in relativistic electrons; enters the equipartition scalings.
  • filling_factor_f = 0.5
    Assumed fraction of the spherical volume that emits; enters the radius and B-field scalings (Eqs. 1-2).
  • electron_index_p = 2.83 +/- 0.02
    Fitted with MCMC to the multiepoch radio spectra; sets the optically thin spectral slope.
  • Rbrk = ~3e17 cm (0.1 pc)
    Break radius where the shock encounters the overdensity; fitted in the MCMC model (Section 4.2, Eq. 10). Degenerate with density normalizations because of the 3-10 yr data gap.
  • rho0_wind = best-fit values in Appendix Figures 7-8 (cgs)
    CSM density normalization before the break; fitted to reproduce the early-time flux.
  • rho0_over = best-fit values in Appendix Figures 7-8 (cgs)
    CSM density normalization after the break; fitted to reproduce the late-time re-brightening.
  • alpha_over = 0.75
    Power-law index of the CSM density profile after the break; fitted. The best-fit value is the same for both ejecta-mass cases.
  • initial_ejecta_velocity_4Msun = 13,888 km/s
    For the 4 Msun ejecta model, the initial ejecta velocity is allowed to be a free parameter to match the early emission (Section 4.2).
  • ejecta_mass_Mej = 1 or 4 Msun
    Two bracketing values chosen from the literature; not fitted to the radio data, but changes the inferred shock dynamics and density profile.
  • explosion_energy_E = 1e51 erg
    Assumed explosion energy for the 1 Msun model. For the 4 Msun model the effective energy is encoded in the fitted initial velocity.
assumptions (8)
  • domain assumption Chevalier (1998) synchrotron self-absorption scalings relate peak flux, shock radius, and magnetic field (Eqs. 1-2 and 9).
    Standard radio-SN model invoked without independent verification for this object; all quantitative results rest on it.
  • domain assumption Thin-shell approximation for the shock dynamics (Eq. 5).
    Assumes the swept-up ejecta and CSM form a thin shell; standard but an approximation.
  • domain assumption Ejecta density profile follows a broken power law with index n = 10 (Eqs. 6-8).
    Assumed for massive stars with radiative envelopes; the exact index affects the deceleration rate.
  • domain assumption The early-time CSM follows a steady wind profile rho proportional to R^-2 (Section 4.2).
    Motivated by the fitted alpha = 1.96 +/- 0.01 to early data; used as the baseline from which the late-time excess is measured.
  • domain assumption The radio spectrum at each epoch is a broken power law with optically thick slope 5/2 and optically thin slope (1-p)/2 (Eq. 9).
    Standard synchrotron self-absorption model; assumes a single homogeneous emitting zone.
  • domain assumption The distance to NGC 7424 is 10.8 Mpc (Section 2).
    Adopted as the average of Hubble-flow and Tully-Fisher distances; affects the luminosity and linear scale.
  • domain assumption The explosion date is MJD 52246 (2001 Dec 3) (Section 2).
    Taken from Ryder et al. 2004 based on the early radio light curve; all ages and velocities depend on it.
  • domain assumption The 2013-2024 flux densities can be rescaled to reference frequencies using alpha = -1 (Figure 2 caption).
    The spectral index is measured to be about -1, so the rescaling is reasonable, but any deviation adds a small systematic error.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The radio re-brightening of the Type IIb SN 2001ig." pith.science (2026). https://pith.science/paper/5R4YUHZA

@misc{pith2026250201740,
  author       = {Pith},
  title        = {Pith review of: The radio re-brightening of the Type IIb SN 2001ig},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5R4YUHZA}},
  note         = {Machine review of arXiv:2502.01740}
}
read the original abstract

We study the late-time evolution of the compact Type IIb SN 2001ig in the spiral galaxy NGC 7424, with new and unpublished archival data from the Australia Telescope Compact Array and the Australian Square Kilometre Array Pathfinder. More than two decades after the SN explosion, its radio luminosity is showing a substantial re-brightening: it is now two orders of magnitude brighter than expected from the standard model of a shock expanding into a uniform circumstellar wind (i.e., with a density scaling as R^{-2}). This suggests that the SN ejecta have reached a denser shell, perhaps compressed by the fast wind of the Wolf-Rayet progenitor or expelled centuries before the final stellar collapse. We model the system parameters (circumstellar density profile, shock velocity, mass loss rate), finding that the denser layer was encountered when the shock reached a distance of ~0.1 pc; the mass-loss rate of the progenitor immediately before the explosion was Mdot/v_w ~ 10^{-7} Msun/yr/(km/s). We compare SN 2001ig with other SNe that have shown late-time re-brightenings, and highlight the opposite behaviour of some extended Type IIb SNe which show instead a late-time flux cut-off.

Figures

Figures reproduced from arXiv: 2502.01740 by the authors.

Figure 1
Figure 1. 9 GHz ATCA contours (cyan) of SN 2001ig observed on 2024 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Flux density of SN 2001ig at 2.4 GHz, 4.85 GHz and 8.55 GHz, based on ATCA observations, compared with the canonical evolution [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. ASKAP detections of SN 2001ig at late times, starting [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Inferred CSM density profiles as a function of shock radius, [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 6
Figure 6. Figure 6: Flux density evolution of the best-fitting spectral model for [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: Corner plot of the best-fitting model parameters and their uncertainties for the 1 [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]
Figure 8
Figure 8. Figure 8: As in Figure [PITH_FULL_IMAGE:figures/full_fig_p014_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

6 extracted references · 4 canonical work pages

  1. [3]

    7772 (December): 2

    IAU Circ. 7772 (December): 2. Podsiadlowski, Ph., J. J. L. Hsu, P. C. Joss, and R. R. Ross. 1993. The progenitor of supernova 1993J: a stripped supergiant in a binary system? Nature 364, no. 6437 (August): 509–511. https://doi.org/10.1038/364509a0. Pooley, David, J. Craig Wheeler, Jozsef Vinkó, Vikram V. Dwarkadas, Tamas Szalai, Jeffrey M. Silverman, Made...

    arXiv 1993

  2. [1018]

    arXiv: astro- ph/0512413 [astro-ph]

    https://doi.org/10.1086/507571. arXiv: astro- ph/0512413 [astro-ph]. Soria, Roberto, Siying Cheng, Manfred W. Pakull, Christian Motch, and Thomas D. Russell. 2024. A multiband look at ultraluminous X-ray sources in NGC 7424. MNRAS 529, no. 2 (April): 1169–1186. https: //doi.org/10.1093/mnras/stae551. arXiv: 2402.09512 [astro-ph.HE]. Soulain, A., F. Millou...

    work page arXiv 2024

  3. [1986]

    ApJ 301 (February): 790

    Radio Supernovae. ApJ 301 (February): 790. https://doi.org/10. 1086/163944. Weiler, Kurt W., Nino Panagia, Marcos J. Montes, and Richard A. Sramek

    work page

  4. [2002]

    ARA&A 40 (January): 387–438

    Radio Emission from Supernovae and Gamma-Ray Bursters. ARA&A 40 (January): 387–438. https://doi.org/10.1146/annurev. astro.40.060401.093744. Weiler, Kurt W., Nino Panagia, and Richard A. Sramek. 1990. Radio Emission from Supernovae. II. SN 1986J: A Different Kind of Type II.ApJ 364 (December): 611. https://doi.org/10.1086/169444. Weiler, Kurt W., Christop...

    arXiv 1990

  5. [2021]

    Radio observations of SN2004dk with VLITE confirm late-time re-brightening

    Radio Observations of SN2004dk with VLITE Confirm Late- time Rebrightening. ApJ 923, no. 1 (December): 32. https://doi.org/10. 3847/1538-4357/ac2154. arXiv: 2101.07348 [astro-ph.HE]. Brethauer, Daniel, Raffaella Margutti, Danny Milisavljevic, and Michael Bi- etenholz. 2020. Six years of luminous X-ray emission from the strongly interacting type-Ib SN 2014...

    work page Pith review arXiv

  6. [7424]

    Uncovering the Putative B-Star Binary Companion of the SN 1993J Progenitor

    IAU Circ. 7772 (December): 1. Filippenko, A. V., and R. Chornock. 2002. Supernovae 2001ig, 2002eg, 2002ew, 2002gc. IAU Circ. 7988 (October): 3. Filippenko, Alexei V. 1988. Supernova 1987K: Type II in Y outh, Type Ib in Old Age. AJ 96 (December): 1941. https://doi.org/10.1086/114940. . 1997. Optical Spectra of Supernovae.ARA&A 35 (January): 309–355. https:...

    work page Pith review arXiv 2002

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.