Part II: Compact Objects | Chapter 3

Supernovae

The explosive deaths of stars: core-collapse mechanisms, thermonuclear Type Ia detonations, light curve physics, and explosive nucleosynthesis

Overview

Supernovae are among the most energetic phenomena in the Universe, releasing up to $3 \times 10^{53}$ erg (mostly in neutrinos) and briefly outshining their entire host galaxy. They are classified into two fundamentally different physical mechanisms: core-collapse supernovae (Types II, Ib, Ic), which mark the death of massive stars ($M \gtrsim 8\,M_\odot$), and thermonuclear supernovae (Type Ia), which arise from the detonation of carbon-oxygen white dwarfs in binary systems.

In this chapter we derive the physics of iron core collapse and the neutrino-driven explosion mechanism, analyze the Chandrasekhar detonation model for Type Ia supernovae, compute light curves powered by radioactive $^{56}\text{Ni}$ decay, derive the Arnett rule, and examine the nucleosynthesis of heavy elements.

1. Core-Collapse Mechanism

Massive stars build up an iron core through successive nuclear burning stages (hydrogen, helium, carbon, neon, oxygen, silicon). Iron ($^{56}\text{Fe}$) has the highest binding energy per nucleon, so no further exothermic fusion is possible. When the iron core exceeds the Chandrasekhar mass, electron degeneracy pressure can no longer support it, and catastrophic collapse ensues.

1.1 Electron Capture and Photodisintegration

Two processes accelerate the collapse. First, as the core density exceeds $\sim 10^{13}$ kg m$^{-3}$, electron capture on protons and nuclei becomes energetically favorable:

$$e^- + p \to n + \nu_e$$

This neutronization removes the electrons that were providing degeneracy pressure. Second, at temperatures exceeding $\sim 5 \times 10^9$ K, photodisintegration of iron nuclei absorbs thermal energy:

$$\gamma + {}^{56}\text{Fe} \to 13\,{}^{4}\text{He} + 4\,n \qquad \Delta E = 124.4\;\text{MeV}$$

Each iron nucleus absorbs 124.4 MeV, converting thermal energy into nuclear binding energy and further reducing the pressure support. The collapse proceeds on a free-fall timescale:

$$\boxed{t_{\text{ff}} = \sqrt{\frac{3\pi}{32 G\rho}} \approx 0.04\left(\frac{\rho}{10^{13}\;\text{kg m}^{-3}}\right)^{-1/2}\;\text{s}}$$

1.2 Core Bounce and Shock Formation

The inner core collapses homologously (maintaining a self-similar velocity profile) until it reaches nuclear density ($\rho_{\text{nuc}} \approx 2.7 \times 10^{17}$ kg m$^{-3}$), where the strong nuclear force becomes repulsive and the equation of state stiffens dramatically. The inner core rebounds (the "core bounce"), launching a shock wave into the infalling outer core.

The initial shock energy is:

$$E_{\text{shock}} \sim \frac{GM_{\text{core}}^2}{R_{\text{core}}} \sim 3 \times 10^{53}\;\text{erg}$$

However, the shock loses energy to photodisintegration of infalling iron ($\sim 8.8$ MeV per nucleon for $\sim 0.5\,M_\odot$ of iron, totaling $\sim 5 \times 10^{51}$ erg) and to neutrino losses. The shock stalls at a radius of $\sim 100\text{--}200$ km within milliseconds.

1.3 Neutrino-Driven Explosion

The proto-neutron star radiates $\sim 3 \times 10^{53}$ erg in neutrinos over $\sim 10$ seconds. A small fraction ($\sim 5\text{--}10\%$) of this neutrino energy is absorbed behind the stalled shock via:

$$\nu_e + n \to p + e^-, \qquad \bar{\nu}_e + p \to n + e^+$$

This neutrino heating revives the shock. The condition for successful explosion (the critical luminosity condition) requires that the neutrino heating rate exceeds the rate of energy advection through the gain region:

$$\boxed{L_\nu^{\text{crit}} \propto (M_{\text{PNS}})^{5/3}\,\dot{M}^{3/5}}$$

where $\dot{M}$ is the mass accretion rate through the shock. Multi-dimensional effects (convection, the standing accretion shock instability — SASI) are crucial for successful explosions in modern simulations.

2. Type Ia Supernovae: Thermonuclear Detonation

Type Ia supernovae result from the thermonuclear explosion of a C/O white dwarf. Unlike core-collapse supernovae, no compact remnant is left behind — the star is completely disrupted.

2.1 Carbon Deflagration-to-Detonation

When carbon ignites at the center of a Chandrasekhar-mass white dwarf ($\rho_c \sim 2 \times 10^{12}$ kg m$^{-3}$,$T \sim 3 \times 10^8$ K), the burning initially proceeds as a subsonic deflagration. The laminar flame speed is determined by thermal conduction:

$$v_{\text{lam}} = \sqrt{\frac{\kappa_{\text{th}} \epsilon_{\text{nuc}}}{\rho c_p T}} \sim 10^{-2}\,c_s$$

where $\kappa_{\text{th}}$ is the thermal diffusivity, $\epsilon_{\text{nuc}}$is the nuclear energy generation rate, and $c_s$ is the sound speed. Rayleigh-Taylor instabilities wrinkle the flame surface, increasing the effective burning rate. The deflagration-to-detonation transition (DDT) occurs when turbulence creates conditions for a detonation wave.

2.2 Nuclear Energy Release

The total nuclear energy released by burning $\sim 1.4\,M_\odot$ of C/O to nuclear statistical equilibrium (predominantly $^{56}\text{Ni}$) is:

$$E_{\text{nuc}} = \frac{M_{\text{WD}}}{\bar{A}\,m_u}\,q_{\text{nuc}} \approx \frac{1.4\,M_\odot}{14\,m_u}\times 1\;\text{MeV} \approx 1.6 \times 10^{51}\;\text{erg}$$

This exceeds the gravitational binding energy of the white dwarf ($\sim 5 \times 10^{50}$ erg), so the entire star is unbound. The kinetic energy of the ejecta is typically $E_K \approx 1.0\text{--}1.3 \times 10^{51}$ erg, corresponding to expansion velocities of $\sim 10{,}000\text{--}15{,}000$ km/s.

3. Supernova Light Curves and the Arnett Rule

After the explosion, the luminosity of a supernova is powered by the radioactive decay of $^{56}\text{Ni}$:

$$^{56}\text{Ni} \xrightarrow{\tau = 8.8\,\text{d}} {}^{56}\text{Co} \xrightarrow{\tau = 111.3\,\text{d}} {}^{56}\text{Fe}$$

3.1 Diffusion Timescale

The photon diffusion time through the expanding ejecta of mass $M_{\text{ej}}$and velocity $v_{\text{ej}}$ is:

$$t_{\text{diff}} = \left(\frac{\kappa M_{\text{ej}}}{v_{\text{ej}} c}\right)^{1/2} \approx 20\left(\frac{\kappa}{0.1\;\text{cm}^2/\text{g}}\right)^{1/2}\left(\frac{M_{\text{ej}}}{M_\odot}\right)^{1/2}\left(\frac{v_{\text{ej}}}{10^4\;\text{km/s}}\right)^{-1/2}\;\text{days}$$

where $\kappa$ is the opacity (dominated by electron scattering and line opacity from iron-group elements).

3.2 The Arnett Rule

David Arnett (1982) showed that at the time of peak luminosity, the instantaneous radioactive energy deposition rate equals the observed luminosity:

$$\boxed{L_{\text{peak}} = \dot{Q}_{\text{Ni}}(t_{\text{peak}}) = M_{\text{Ni}}\,\epsilon_{\text{Ni}}\,e^{-t_{\text{peak}}/\tau_{\text{Ni}}}}$$

where $\epsilon_{\text{Ni}} = 3.9 \times 10^{10}$ erg g$^{-1}$ s$^{-1}$is the specific energy generation rate of $^{56}\text{Ni}$ decay. This provides a direct measurement of the $^{56}\text{Ni}$ mass from the peak luminosity. Typical Type Ia supernovae produce $\sim 0.6\,M_\odot$of $^{56}\text{Ni}$, while normal core-collapse supernovae produce $\sim 0.07\,M_\odot$.

3.3 Late-Time Radioactive Tail

At late times ($t \gg t_{\text{diff}}$), the ejecta become optically thin and the luminosity directly tracks the $^{56}\text{Co}$ decay:

$$L(t) \approx M_{\text{Ni}}\,\epsilon_{\text{Co}}\,e^{-t/\tau_{\text{Co}}}$$

This produces the characteristic linear decline (0.98 mag per 100 days) in the magnitude light curve, which is the signature of $^{56}\text{Co}$ radioactive decay.

4. Explosive Nucleosynthesis

Supernovae are the primary site for the synthesis of elements heavier than oxygen. The nucleosynthetic yields depend on the peak temperature and density reached by each mass shell during the explosion.

4.1 Nuclear Statistical Equilibrium

At temperatures above $\sim 5 \times 10^9$ K, nuclear reactions are fast enough to reach nuclear statistical equilibrium (NSE). In NSE, the abundance of each isotope $(Z, A)$ is determined by the Saha equation for nuclei:

$$Y(Z,A) = G(Z,A)\left(\frac{\rho}{m_u}\right)^{A-1}\frac{A^{3/2}}{2^A}\left(\frac{2\pi\hbar^2}{m_u k_B T}\right)^{3(A-1)/2} Y_p^Z Y_n^{A-Z}\,e^{B(Z,A)/(k_B T)}$$

where $B(Z,A)$ is the nuclear binding energy and $G$ is the partition function. At the electron fractions typical of the inner ejecta ($Y_e \approx 0.5$), the dominant product is $^{56}\text{Ni}$ (28 protons, 28 neutrons).

4.2 The r-Process

The rapid neutron-capture process (r-process) produces about half of the elements heavier than iron. It requires an extremely neutron-rich environment ($Y_e < 0.25$) with high neutron density ($n_n > 10^{20}$ cm$^{-3}$). The r-process path runs along nuclei far from the valley of stability, where neutron captures are faster than beta decays:

$$\tau_{\text{n-capture}} = \frac{1}{n_n \langle\sigma v\rangle} \ll \tau_\beta$$

The primary r-process site has been debated for decades. The 2017 detection of a kilonova (AT2017gfo) associated with the neutron star merger GW170817 confirmed that neutron star mergers are a major r-process site, producing $\sim 0.05\,M_\odot$ of r-process material including elements like gold, platinum, and uranium.

5. Supernova Remnant Evolution

After the explosion, the ejecta interact with the surrounding interstellar medium, forming a supernova remnant (SNR) that evolves through several distinct phases.

5.1 The Sedov-Taylor Phase

Once the swept-up ISM mass exceeds the ejecta mass, the remnant enters the adiabatic (Sedov-Taylor) phase. Dimensional analysis gives the expansion law:

$$\boxed{R(t) = \xi_0\left(\frac{E}{\rho_0}\right)^{1/5} t^{2/5}}$$

where $\xi_0 \approx 1.15$ is a dimensionless constant from the self-similar Sedov-Taylor solution, $E$ is the explosion energy, and $\rho_0$ is the ambient ISM density. The shock velocity decelerates as $v_s \propto t^{-3/5}$.

5.2 Transition to Radiative Phase

When the shock velocity drops below $\sim 200$ km/s (after about $3 \times 10^4$ years for typical parameters), the post-shock gas cools efficiently by line emission. A dense, cool shell forms behind the shock, and the remnant enters the radiative (snowplow) phase with $R \propto t^{2/7}$. Eventually the expansion velocity drops to the ISM sound speed and the remnant merges with the ISM after $\sim 10^6$ years.

Applications

Cosmological Distance Ladder

Type Ia supernovae, standardized through the width-luminosity relation (Phillips relation), serve as the premier distance indicators at cosmological distances. The discovery that distant Type Ia supernovae are dimmer than expected in a decelerating universe led to the revelation that the expansion of the Universe is accelerating, implying the existence of dark energy ($\Omega_\Lambda \approx 0.7$).

Chemical Evolution of Galaxies

Core-collapse supernovae primarily produce alpha elements (O, Ne, Mg, Si, S, Ca, Ti), while Type Ia supernovae produce iron-peak elements with a time delay of$\sim 1$ Gyr. The [$\alpha$/Fe] ratio in stellar atmospheres serves as a cosmic clock, tracing the chemical enrichment history of galaxies. The knee in the [$\alpha$/Fe] vs [Fe/H] relation marks the onset of Type Ia contributions.

Supernova Neutrino Detection

A Galactic core-collapse supernova would produce $\sim 10^{58}$ neutrinos, with detectors like Super-Kamiokande expecting $\sim 10{,}000$ events and Hyper-Kamiokande $\sim 100{,}000$ events. The neutrino signal arrives minutes to hours before the optical emission (which must diffuse through the stellar envelope), enabling the Supernova Early Warning System (SNEWS) to alert optical telescopes before the supernova becomes visible. The neutrino light curve encodes the nuclear equation of state, the formation of a neutron star or black hole, and possible exotic physics (neutrino mass hierarchy, axion emission). The next Galactic supernova will be a once-in-a-generation event for neutrino physics.

Pair-Instability Supernovae

Very massive stars ($140\text{--}260\,M_\odot$) with low metallicity may encounter the pair-instability: at core temperatures exceeding $10^9$ K, photons produce electron-positron pairs, reducing the radiation pressure support and triggering explosive oxygen and silicon burning. Unlike core-collapse supernovae, pair- instability supernovae completely disrupt the star, leaving no compact remnant and producing up to $50\,M_\odot$ of $^{56}\text{Ni}$. These events may have been common among the first generation of stars (Population III) and could be detectable by JWST at $z \sim 10\text{--}20$.

Historical Notes

Supernovae have been recorded throughout human history. The "guest star" of 1054 CE, observed by Chinese astronomers, produced the Crab Nebula and its central pulsar. Tycho Brahe's supernova of 1572 and Kepler's supernova of 1604 were the last Galactic supernovae visible to the naked eye. SN 1987A in the Large Magellanic Cloud was the nearest supernova since 1604 and provided the first detection of supernova neutrinos (by Kamiokande-II, IMB, and Baksan), confirming the neutrino-driven mechanism. The classification scheme (Type I vs Type II based on hydrogen lines, with sub-types Ia, Ib, Ic) was established by Rudolph Minkowski in 1941 and refined by Fritz Zwicky and subsequent workers.

The field of supernova cosmology began with the calibration of Type Ia supernovae as standard candles. Mark Phillips (1993) discovered the width-luminosity relation, and the teams led by Saul Perlmutter (Supernova Cosmology Project) and Brian Schmidt/Adam Riess (High-z Supernova Search Team) independently announced in 1998 that distant Type Ia supernovae were fainter than expected, implying that the expansion of the Universe is accelerating. This discovery of dark energy is among the most important findings in the history of physics, earning Perlmutter, Schmidt, and Riess the 2011 Nobel Prize.

The computational modeling of core-collapse supernovae has been one of the grand challenges of astrophysics for over 50 years. The first successful 3D simulations producing explosions with realistic energies and nucleosynthetic yields were achieved in the 2010s and 2020s, using codes such as FLASH, CHIMERA, Fornax, and PROMETHEUS. These simulations confirm that multi-dimensional effects — convection, the SASI, and neutrino-driven turbulence — are essential for the explosion mechanism. The neutrino signal from SN 1987A (24 events across three detectors over $\sim 10$ seconds) remains the only direct observation of core-collapse neutrinos, but it confirmed the basic theoretical picture of a proto-neutron star cooling through neutrino emission. The next Galactic supernova will be observed with detectors $\sim 1000$ times more sensitive, revolutionizing our understanding of the explosion mechanism.

Supernova Classification in Detail

The supernova classification scheme is based primarily on spectral features near maximum light, with secondary subdivision by light curve shape:

Type Ia: No hydrogen, strong Si II absorption at 6150 $\text{\AA}$. Thermonuclear detonation of C/O white dwarf. Standardizable candles with peak $M_B \approx -19.3$ mag.

Type Ib: No hydrogen, no Si II, strong He I lines. Core collapse of massive star that lost its hydrogen envelope (Wolf-Rayet progenitor or binary stripping).

Type Ic: No hydrogen, no Si II, no He I. Core collapse of star that lost both hydrogen and helium envelopes. Broad-lined Ic (Ic-BL) are associated with long-duration GRBs.

Type II-P: Hydrogen lines present, extended luminosity plateau ($\sim 100$ days) powered by hydrogen recombination wave propagating inward through the massive hydrogen envelope.

Type II-L: Hydrogen lines, linear decline after maximum (less massive hydrogen envelope than II-P).

Type IIn: Narrow hydrogen emission lines from interaction with dense circumstellar medium ejected in pre-supernova eruptions. Can be extremely luminous ($L > 10^{44}$ erg/s) when the CSM is massive.

Superluminous supernovae (SLSNe) have peak luminosities exceeding$7 \times 10^{43}$ erg/s ($M < -21$ mag) and are powered by either magnetar spin-down energy injection or massive CSM interaction. They are predominantly found in low-metallicity dwarf galaxies and may trace the most massive stellar explosions.

Time-Domain Supernova Science

Modern wide-field surveys (ZTF, ATLAS, and soon LSST) are discovering supernovae at unprecedented rates. The Zwicky Transient Facility detects $\sim 3{,}000$supernovae per year, while LSST is expected to find $\sim 100{,}000$per year. Rapid classification using machine learning algorithms enables prompt follow-up spectroscopy and multi-wavelength observations. Early detection (within hours of explosion) reveals the shock breakout — the moment the supernova shock wave reaches the stellar surface — which constrains the progenitor radius and the presence of circumstellar material. The growing sample of well-observed supernovae is refining the Phillips relation for Type Ia standardization and revealing new subclasses that challenge existing theoretical models.

Gravitational Waves from Core-Collapse Supernovae

Core-collapse supernovae produce gravitational waves through several mechanisms: core bounce ($f \sim 700$ Hz), convective overturn and SASI ($f \sim 100\text{--}300$ Hz), and proto-neutron star oscillations ($f \sim 1\text{--}2$ kHz). The gravitational wave energy is $\sim 10^{-8}\text{--}10^{-7}\,M_\odot c^2$, much less than from compact binary mergers but detectable for Galactic events. A supernova at $10$ kpc would produce strains of $h \sim 10^{-21}\text{--}10^{-20}$at LIGO/Virgo, providing direct information about the dynamics of core collapse and the formation of the compact remnant. The coincident detection of gravitational waves and neutrinos from the next Galactic supernova will be a landmark multi-messenger event, probing the deepest interior of the collapsing star.

Computational Exploration

The following simulation models supernova light curves powered by radioactive nickel decay, compares Type Ia and core-collapse events, computes the Sedov-Taylor blast wave evolution, and illustrates explosive nucleosynthesis yields.

Supernova Light Curves, Blast Waves, and Nucleosynthesis

Python

script.py230 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a0a')
for ax in axes.flat:
    ax.set_facecolor('#0f0f0f')
    ax.tick_params(colors='white')
    ax.xaxis.label.set_color('white')
    ax.yaxis.label.set_color('white')
    ax.title.set_color('white')
    for spine in ax.spines.values():
        spine.set_color('#333')

# Constants
M_sun = 1.989e33  # grams
L_sun = 3.828e33  # erg/s
c_light = 3e10  # cm/s
day = 86400.0  # seconds

# --- Plot 1: Type Ia and II light curves ---
ax1 = axes[0, 0]

def nickel_heating(t, M_Ni):
    """Energy deposition rate from Ni-56 and Co-56 decay (erg/s)"""
    tau_Ni = 8.8 * day
    tau_Co = 111.3 * day
    eps_Ni = 3.9e10  # erg/g/s
    eps_Co = 6.78e9  # erg/g/s
    # Ni-56 decay
    q_Ni = M_Ni * eps_Ni * np.exp(-t / tau_Ni)
    # Co-56 decay (produced from Ni-56)
    q_Co = M_Ni * eps_Co * (tau_Co / (tau_Co - tau_Ni)) * (np.exp(-t / tau_Co) - np.exp(-t / tau_Ni))
    return q_Ni + q_Co

def arnett_lightcurve(t, M_ej, v_ej, M_Ni, kappa=0.1):
    """Simplified Arnett (1982) model"""
    # Diffusion timescale
    t_m = np.sqrt(2 * kappa * M_ej / (v_ej * c_light))
    x = t / t_m
    L = np.zeros_like(t)
    dt_int = t[1] - t[0] if len(t) > 1 else day
    for i in range(len(t)):
        # Integrate heating with diffusion kernel
        t_prime = np.linspace(0, t[i], max(2, int(t[i]/dt_int) + 1))
        if len(t_prime) < 2:
            continue
        x_p = t_prime / t_m
        x_i = t[i] / t_m
        integrand = nickel_heating(t_prime, M_Ni) * np.exp(-(x_i**2 - x_p**2)) * 2 * x_p / t_m
        L[i] = np.trapezoid(integrand, t_prime)
    return L

t_days = np.linspace(1, 300, 300)
t_sec = t_days * day

# Type Ia: M_ej ~ 1.4 Msun, v ~ 10000 km/s, M_Ni ~ 0.6 Msun
L_Ia = arnett_lightcurve(t_sec, 1.4*M_sun, 1e9, 0.6*M_sun, kappa=0.1)
mag_Ia = -2.5 * np.log10(np.maximum(L_Ia, 1e30) / (3e43)) - 19.3  # Approximate absolute mag

# Type II-P: M_ej ~ 10 Msun, v ~ 4000 km/s, M_Ni ~ 0.07 Msun
L_IIP = arnett_lightcurve(t_sec, 10*M_sun, 4e8, 0.07*M_sun, kappa=0.3)
# Add plateau from recombination
T_recomb = 6000  # K
R_phot = 4e8 * t_sec  # photosphere expanding
L_plateau = 4 * np.pi * R_phot**2 * 5.67e-5 * T_recomb**4
# Plateau dominates for first ~100 days
plateau_mask = t_days < 100
L_IIP_total = np.copy(L_IIP)
L_IIP_total[plateau_mask] = np.maximum(L_IIP[plateau_mask], L_plateau[plateau_mask] * 0.3)

# Type Ic-BL (broad-lined, like GRB-SN): M_Ni ~ 0.4 Msun
L_IcBL = arnett_lightcurve(t_sec, 3*M_sun, 2e9, 0.4*M_sun, kappa=0.07)

ax1.semilogy(t_days, np.maximum(L_Ia, 1e37) / L_sun, color='#ff4444', linewidth=2.5, label='Type Ia')
ax1.semilogy(t_days, np.maximum(L_IIP_total, 1e37) / L_sun, color='#ffa500', linewidth=2.5, label='Type II-P')
ax1.semilogy(t_days, np.maximum(L_IcBL, 1e37) / L_sun, color='#44bbff', linewidth=2.5, label='Type Ic-BL')

# Radioactive decay slope
t_late = np.linspace(120, 300, 50)
L_Co_ref = 1e8 * np.exp(-t_late * day / (111.3 * day))
ax1.semilogy(t_late, L_Co_ref, 'w--', linewidth=1, alpha=0.4, label='Co-56 decay slope')

ax1.set_xlabel('Days since explosion', fontsize=12)
ax1.set_ylabel('Luminosity (L_sun)', fontsize=12)
ax1.set_title('Supernova Light Curves', fontsize=14, fontweight='bold')
ax1.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax1.set_xlim(0, 300)
ax1.set_ylim(1e5, 2e10)
ax1.grid(True, alpha=0.15)

# --- Plot 2: Sedov-Taylor blast wave ---
ax2 = axes[0, 1]

E_sn = 1e51  # erg
n_ism_values = [0.1, 1.0, 10.0]  # cm^-3
colors_ism = ['#44bbff', '#ffa500', '#ff4444']
m_H = 1.67e-24  # grams

t_yr = np.logspace(1, 6, 500)
t_s = t_yr * 3.156e7

for n_ism, col in zip(n_ism_values, colors_ism):
    rho_ism = n_ism * m_H
    xi_0 = 1.15
    R_ST = xi_0 * (E_sn / rho_ism)**0.2 * t_s**0.4  # cm
    R_pc = R_ST / 3.086e18

# Transition to snowplow at ~30,000 yr
    v_shock = 0.4 * R_ST / t_s
    t_rad = t_yr[np.argmin(np.abs(v_shock - 2e7))]  # v ~ 200 km/s

ax2.loglog(t_yr, R_pc, color=col, linewidth=2.5, label=f'n = {n_ism} cm^-3')
    # Mark radiative transition
    idx_rad = np.argmin(np.abs(v_shock - 2e7))
    if idx_rad > 0 and idx_rad < len(t_yr) - 1:
        ax2.plot(t_yr[idx_rad], R_pc[idx_rad], 'o', color=col, markersize=8, zorder=5)

# Reference slopes
t_ref = np.logspace(2, 4.5, 50)
R_ref = 0.3 * (t_ref / 1e3)**0.4
ax2.loglog(t_ref, R_ref, 'w--', linewidth=1, alpha=0.3, label=r'$R propto t^{2/5}$')

t_ref2 = np.logspace(4.5, 6, 50)
R_ref2 = R_ref[-1] * (t_ref2 / t_ref[-1])**(2.0/7.0)
ax2.loglog(t_ref2, R_ref2, 'w:', linewidth=1, alpha=0.3, label=r'$R propto t^{2/7}$')

ax2.set_xlabel('Time (years)', fontsize=12)
ax2.set_ylabel('Radius (pc)', fontsize=12)
ax2.set_title('Supernova Remnant Expansion', fontsize=14, fontweight='bold')
ax2.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax2.set_xlim(10, 1e6)
ax2.set_ylim(0.01, 100)
ax2.grid(True, alpha=0.15)

# --- Plot 3: Nucleosynthesis yields ---
ax3 = axes[1, 0]

elements = ['C', 'O', 'Ne', 'Mg', 'Si', 'S', 'Ca', 'Ti', 'Cr', 'Fe', 'Ni', 'Zn']
atomic_num = [6, 8, 10, 12, 14, 16, 20, 22, 24, 26, 28, 30]

# Yields in solar masses (approximate)
yields_CCSNe = [0.1, 1.5, 0.3, 0.15, 0.1, 0.04, 0.01, 0.002, 0.005, 0.07, 0.04, 0.001]
yields_Ia = [0.05, 0.14, 0.005, 0.01, 0.15, 0.08, 0.015, 0.003, 0.01, 0.6, 0.05, 0.001]

x_pos = np.arange(len(elements))
width = 0.35

bars1 = ax3.bar(x_pos - width/2, yields_CCSNe, width, color='#ffa500', alpha=0.8, label='Core-collapse (25 M_sun)')
bars2 = ax3.bar(x_pos + width/2, yields_Ia, width, color='#ff4444', alpha=0.8, label='Type Ia')

ax3.set_xlabel('Element', fontsize=12)
ax3.set_ylabel('Yield (M_sun)', fontsize=12)
ax3.set_title('Nucleosynthesis Yields per Event', fontsize=14, fontweight='bold')
ax3.set_xticks(x_pos)
ax3.set_xticklabels(elements, fontsize=10)
ax3.set_yscale('log')
ax3.set_ylim(5e-4, 5)
ax3.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax3.grid(True, alpha=0.15, axis='y')

# --- Plot 4: Peak Ni mass vs peak luminosity (Arnett rule) ---
ax4 = axes[1, 1]

M_Ni_range = np.linspace(0.01, 1.2, 100)  # solar masses
# Arnett rule: L_peak ~ eps_Ni * M_Ni * exp(-t_peak/tau_Ni)
# For Type Ia, t_peak ~ 17 days
t_peak_Ia = 17 * day
eps_Ni_rate = 3.9e10  # erg/g/s
tau_Ni = 8.8 * day

L_peak = M_Ni_range * M_sun * eps_Ni_rate * np.exp(-t_peak_Ia / tau_Ni)
M_bol = -2.5 * np.log10(L_peak / 3e35) + 5 * np.log10(10 / 10)
abs_mag = -2.5 * np.log10(L_peak / (3.02e35)) - 0.2

ax4.plot(M_Ni_range, abs_mag, color='#ffa500', linewidth=3, label='Arnett rule prediction')

# Observed Type Ia SNe
obs_Ni = [0.1, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1]
obs_mag = [-17.5, -18.5, -18.9, -19.1, -19.3, -19.45, -19.55, -19.65, -19.75, -19.8]
ax4.scatter(obs_Ni, obs_mag, color='#ffdd00', s=60, zorder=5, edgecolor='white', linewidth=0.5, label='Observed Type Ia')

# Phillips relation note
ax4.annotate('Brighter = more Ni\n= broader LC', xy=(0.9, -19.65), xytext=(0.4, -18.0),
            color='#ff6644', fontsize=10, arrowprops=dict(arrowstyle='->', color='#ff6644'))

ax4.set_xlabel('Ni-56 Mass (M_sun)', fontsize=12)
ax4.set_ylabel('Peak Absolute Magnitude', fontsize=12)
ax4.set_title('Arnett Rule: Ni Mass vs Peak Luminosity', fontsize=14, fontweight='bold')
ax4.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax4.invert_yaxis()
ax4.grid(True, alpha=0.15)

# Print summary
print("=" * 65)
print("SUPERNOVAE: KEY RESULTS")
print("=" * 65)

print("\nCore-collapse supernova energetics:")
print(f"  Gravitational energy: ~3e53 erg")
print(f"  Neutrino energy: ~3e53 erg (99%)")
print(f"  Kinetic energy: ~1e51 erg")
print(f"  Radiated energy: ~1e49 erg (0.01%)")

print("\nType Ia supernova:")
print(f"  Nuclear energy: ~1.6e51 erg")
print(f"  Kinetic energy: ~1.1e51 erg")
print(f"  Typical Ni-56 mass: 0.6 M_sun")

print("\nFree-fall timescales:")
for rho in [1e13, 1e14, 1e15]:
    t_ff = np.sqrt(3*np.pi/(32 * 6.674e-8 * rho))
    print(f"  rho = {rho:.0e} kg/m^3: t_ff = {t_ff:.3f} s")

print("\nSedov-Taylor remnant sizes (E=1e51 erg, n=1 cm^-3):")
for t_age in [100, 1000, 10000, 100000]:
    rho_0 = 1.67e-24  # g/cm^3
    R = 1.15 * (1e51/rho_0)**0.2 * (t_age*3.156e7)**0.4
    print(f"  t = {t_age:6d} yr: R = {R/3.086e18:.1f} pc")

print("\nRadioactive decay chains:")
print(f"  Ni-56 -> Co-56: tau = 8.8 d, eps = 3.9e10 erg/g/s")
print(f"  Co-56 -> Fe-56: tau = 111.3 d, eps = 6.78e9 erg/g/s")

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
print("\nPlot saved to output.png")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Chapter Summary

Core-collapse supernovae result from the gravitational collapse of iron cores in massive stars. Electron capture and photodisintegration trigger the collapse on a free-fall timescale $t_{\text{ff}} \sim 0.04$ s. Core bounce launches a shock that stalls and is revived by neutrino heating, with multi-dimensional instabilities playing a critical role.

Type Ia supernovae arise from thermonuclear detonation of C/O white dwarfs near the Chandrasekhar mass. The Arnett rule connects peak luminosity directly to the $^{56}\text{Ni}$ mass synthesized.

Supernovae are the primary engines of chemical evolution, with core-collapse events producing alpha elements and Type Ia events producing iron-peak elements. Supernova remnants expand through the Sedov-Taylor phase ($R \propto t^{2/5}$) before transitioning to the radiative snowplow phase.

← Neutron Stars Gamma-Ray Bursts →