[Paper Review] Supernova Explosions in Winds and Bubbles, with Applications to SN 1987A

Vikram V. Dwarkadas|arXiv (Cornell University)|Dec 21, 2006

Gamma-ray bursts and supernovae5 citations

TL;DR

This paper models supernova explosions in circumstellar wind-blown bubbles, focusing on SN 1987A’s shock evolution and emission. It uses hydrodynamic simulations to show that the shock's interaction with the equatorial ring causes a luminosity bump in X-rays and a nonlinear radio rise, though a constant particle energy density parameter best fits observed radio growth, suggesting limitations in spherical symmetry or injection models.

ABSTRACT

Massive stars can significantly modify the surrounding medium during their lifetime. When the stars explode as supernovae, the resulting shock wave expands within this modified medium and not within the interstellar medium. We explore the evolution of the medium around massive stars, and the expansion of the shock wave within this medium. We then apply these results to understanding the expansion of the shock wave in the ambient medium surrounding SN 1987A, and the evolution of the radio and X-ray emission in this case.

Motivation & Objective

• To understand how supernova shock waves evolve in circumstellar media shaped by massive stars' pre-explosion winds.
• To explain the observed radio and X-ray light curves of SN 1987A by modeling shock propagation through a structured, wind-blown bubble environment.
• To assess the limitations of spherically symmetric hydrodynamic models in reproducing the observed emission evolution, particularly the linear radio rise and X-ray bump.
• To identify key physical parameters—such as magnetic field, particle energy density, and shock speed—that govern emission signatures.
• To lay the groundwork for future 3D models including ionization, HII regions, and aspherical structures like the equatorial ring.

Proposed method

• Uses self-similar hydrodynamic solutions for wind-blown bubbles (W77 model) to describe the pre-supernova circumstellar medium, with shock radius scaling as R_sh ∝ t^{3/5}.
• Applies energy and momentum conservation to compute the shock wave's expansion through the bubble, including the reverse shock and interaction with the dense equatorial ring.
• Models radio emission using synchrotron radiation from relativistic electrons, assuming a power-law injection and scaling the emissivity with thermal pressure and magnetic field.
• Computes hard X-ray emission via the CHIANTI atomic code, assuming electron temperature ≈ 0.02 × ion temperature, and compares with observed spectra.
• Fits model parameters (e.g., K, B) to observed light curves, finding that constant K or B yields best agreement with the observed linear radio rise.
• Assesses model limitations, including spherical symmetry, constant temperature ratios, and missing soft X-ray contributions from shock-protrusion interactions.

Experimental results

Research questions

• RQ1 How does the shock wave from a core-collapse supernova evolve when it expands into a pre-existing wind-blown bubble?
• RQ2 Why does the radio luminosity of SN 1987A increase linearly over time, contrary to the model’s prediction of quadratic or faster growth?
• RQ3 What causes the observed bump in hard X-ray luminosity upon shock-ring interaction, and why is it not reproduced in the spherically symmetric model?
• RQ4 How do assumptions about particle injection and magnetic field evolution affect the fit to observed radio and X-ray light curves?
• RQ5 To what extent do hydrodynamic asymmetries and non-spherical structures like the equatorial ring influence the emission signatures?

Key findings

• The shock radius in the wind-blown bubble evolves as R_sh ∝ t^{3/5}, consistent with the W77 self-similar solution, with a predicted radius of ~48.8 pc after 10 Myr of wind activity.
• The model predicts a radio luminosity that increases quadratically or faster with time, but the observed data show a linear rise, indicating a mismatch in the assumed particle injection or energy density scaling.
• Best-fit parameters suggest that either the particle energy density parameter K or the magnetic field B must be approximately constant over time to match the observed linear radio growth, implying K ∝ P and B ∝ P^{0.5} for pressure P.
• The hard X-ray luminosity model shows a significant bump upon shock-ring interaction, which is not seen in observations, suggesting that the reverse shock and ring interaction are not well captured in the spherically symmetric approximation.
• The model provides a reasonable fit to the X-ray data before shock-ring interaction, but fails to reproduce the post-interaction behavior, highlighting the need for 3D hydrodynamic and ionization modeling.
• The electron-to-ion temperature ratio is assumed constant at 0.02, but this may vary, and future models should account for variable cooling and non-equilibrium ionization effects.