Magnetar Engines in Broad-lined Type Ic Supernovae and a Unified Picture for Magnetar-powered Stripped-envelope Supernovae

1. Problem Statement

Broad-lined Type Ic supernovae (SNe Ic-BL) are energetic transients often associated with long-duration gamma-ray bursts (lGRBs). While their high explosion energies (up to $\sim 10^{52}$ erg) suggest an engine beyond the standard neutrino-driven mechanism, the exact nature of this engine remains debated.

• The Energy Puzzle: Previous studies fitting SNe Ic-BL lightcurves with pure $^{56}$Ni decay models failed to reproduce both the peak luminosity and the late-time decline simultaneously. Hybrid models including a magnetar engine were proposed, but they often yielded inconsistent parameters (e.g., "energy puzzles" where inferred explosion energies exceeded the rotational energy of the inferred magnetar).
• Lack of Statistical Cohesion: There was no comprehensive statistical analysis of a large sample of SNe Ic-BL to constrain magnetar parameters (initial spin period $P_i$ and magnetic field $B_p$) and ejecta properties ($M_{ej}$, $M_{Ni}$) simultaneously. Furthermore, the physical connections between SNe Ic-BL, other magnetar-powered transients (SLSNe Ic, FBOTs), and ordinary stripped-envelope supernovae (SESNe) were not unified.

2. Methodology

The authors developed a semi-analytical model and performed a statistical analysis on the largest available sample of SNe Ic-BL.

• Model Framework:

  • A hybrid model combining magnetar spin-down energy injection and radioactive decay of $^{56}$Ni.
  • The magnetar is treated as a newborn millisecond magnetar losing rotational energy via magnetic dipole radiation.
  • The model accounts for the conversion of magnetar rotational energy into ejecta kinetic energy (mechanical heating) and thermal radiation, alongside adiabatic cooling and gamma-ray trapping.
  • Crucially, the initial kinetic energy ($E_{kin,i}$) is constrained by a linear relationship with ejecta mass derived from ordinary SESNe (neutrino-driven mechanism), rather than treating it as a free parameter. This addresses the "energy puzzle" by ensuring the magnetar's contribution is physically consistent with the total explosion energy.

• Data Collection:

  • A sample of 80 spectroscopically confirmed SNe Ic-BL was collected from literature (up to March 2026).
  • The sample includes 11 GRB-associated SNe (GRB-SNe) and 69 non-GRB SNe (Non-GRB-SNe).
  • Data requirements included multi-band lightcurves with coverage of the rising phase, peak, and late-time $^{56}$Ni tail.

• Fitting Procedure:

  • Used Markov Chain Monte Carlo (MCMC) methods (emcee) to fit multi-band lightcurves.
  • Parameters fitted: $P_i$, $B_p$, ejecta mass ($M_{ej}$), $^{56}$Ni mass ($M_{Ni}$), gamma-ray opacity ($\kappa_\gamma$), floor temperature ($T_{floor}$), and explosion time.
  • Model comparison was performed using the Akaike Information Criterion (AIC) against pure $^{56}$Ni models.

3. Key Results

A. Magnetar Engine Validation

• The magnetar + $^{56}$Ni model provides a significantly better fit to the data than pure $^{56}$Ni models (AIC differences $\Delta AIC \sim 10^4 - 10^5$).
• The model successfully reproduces the peak excess luminosity and the late-time decline for all 80 SNe.

B. Inferred Physical Parameters (Medians with 1$\sigma$)

• Initial Spin Period ($P_i$): $\sim 2.04^{+1.84}_{-0.96}$ ms.
• Magnetic Field ($B_p$): $\sim 3.96^{+3.28}_{-1.40} \times 10^{15}$ G.
• Ejecta Mass ($M_{ej}$): $\sim 2.30^{+1.48}_{-1.02} M_\odot$.
• $^{56}$Ni Mass ($M_{Ni}$): $\sim 0.18^{+0.14}_{-0.09} M_\odot$.
• $^{56}$Ni Fraction ($f_{Ni}$): $\sim 0.08^{+0.05}_{-0.04}$.

C. Statistical Correlations

• $M_{ej} - P_i$ Anti-correlation: A strong, statistically significant anti-correlation exists ($P_i \propto M_{ej}^{-0.95}$). This implies that progenitors with larger ejecta masses produce faster-spinning magnetars (or vice versa), likely due to angular momentum conservation in the collapsing core.
• $M_{Ni} - M_{ej}$ Correlation: A strong positive correlation exists, suggesting a consistent $^{56}$Ni production efficiency across the sample.

D. GRB-SNe vs. Non-GRB-SNe

• No Intrinsic Difference: Statistical tests (Energy Distance test) indicate no significant difference between the underlying populations of GRB-SNe and Non-GRB-SNe.
• Observational Bias: GRB-SNe appear slightly brighter and have slightly higher $M_{Ni}$ and lower $B_p$. The authors attribute this to observational bias: GRB-SNe are detected at cosmological distances where bright afterglows outshine fainter SNe, selecting only the intrinsically brightest events.

E. Unified Picture with Other Transients

The authors compared SNe Ic-BL with Superluminous SNe (SLSNe Ic) and Fast Blue Optical Transients (FBOTs):

• Universal $M_{ej} - P_i$ Correlation: A universal anti-correlation ($P_i \propto M_{ej}^{-0.49}$) spans SNe Ic-BL, SLSNe Ic, and FBOTs, suggesting a common physical origin.
• Distinct Parameter Spaces:
  • FBOTs: Low ejecta mass ($M_{ej} \lesssim 0.5 M_\odot$), slower spin ($P_i \sim 10$ ms).
  • SLSNe Ic: High ejecta mass ($M_{ej} \gtrsim 0.5 M_\odot$), slower spin ($P_i \sim 3$ ms), lower $B_p$ ($\lesssim 10^{15}$ G).
  • SNe Ic-BL: High ejecta mass ($M_{ej} \gtrsim 0.5 M_\odot$), fastest spin ($P_i \sim 2$ ms), highest $B_p$ ($\gtrsim 10^{15}$ G).
• Lightcurve Morphology: The ratio of diffusion timescale to spin-down timescale ($t_{diff}/t_{sd}$) dictates the observational class.
  • SLSNe Ic/FBOTs: $t_{diff} \sim t_{sd}$ (Radiatively powered).
  • SNe Ic-BL: $t_{diff} \gg t_{sd}$ (Mechanically powered; energy goes into kinetic energy, not peak luminosity).

4. Proposed Unified Progenitor Framework

The paper proposes a unified binary evolution scenario:

• Origin: All SESNe (ordinary and magnetar-powered) originate from helium stars in binary systems formed via stable Roche-lobe overflow (RLOF) or common-envelope evolution (CEE).
• Differentiation Mechanism: The presence of a rapidly rotating magnetar depends on tidal synchronization in close binaries.
  • SNe Ic-BL & SLSNe Ic: Originate from massive helium stars ($4-40 M_\odot$) in tight binaries ($P_{orb} \lesssim 1-2$ days). Tidal forces spin up the progenitor, leading to a fast-rotating magnetar upon collapse.
  • FBOTs: Originate from low-mass helium stars ($2-4 M_\odot$) in close binaries, potentially corresponding to magnetar-powered ultra-stripped SNe.
  • Ordinary SESNe: Originate from wider binaries where tidal spin-up is inefficient; the resulting neutron stars rotate slowly, and the explosion is primarily neutrino-driven (magnetar-aided).

5. Significance and Conclusion

• Resolution of the Energy Puzzle: By constraining the initial kinetic energy via the $M_{ej}$ relation, the authors demonstrate that magnetar engines can naturally explain the high explosion energies of SNe Ic-BL without requiring unphysical parameters.
• Unified Classification: The study establishes a physical classification scheme for SESNe based on the dominant energy source (neutrino-driven vs. magnetar-powered) and the efficiency of energy conversion (radiative vs. mechanical).
• GRB Connection: The results strongly support the hypothesis that lGRBs are powered by millisecond magnetars formed in SNe Ic-BL, with the jet formation likely driven by the ultra-high magnetic fields ($B_p \gtrsim 10^{15}$ G) and rapid rotation of the central engine.
• Progenitor Insight: The work links the diversity of transient observations directly to binary orbital parameters and progenitor masses, providing a coherent evolutionary path for massive stars leading to various explosion types.

This paper effectively unifies the phenomenology of SNe Ic-BL, SLSNe Ic, and FBOTs under a single magnetar engine framework, distinguishing them by progenitor mass, binary separation, and the resulting magnetar properties.