---
title: Radiative Supernova Bubbles
url: https://www.emergentmind.com/topics/radiative-supernovae-bubbles
type: topic
---

# Radiative Supernova Bubbles

A radiative supernova bubble is a region of multi-phase, shock-heated gas that inflates following the explosive death of a massive star, where the post-shock plasma and swept-up shell undergo substantial radiative cooling, dramatically reducing the fraction of supernova energy retained in kinetic or thermal form. These bubbles play an essential role in the energetics and structure of the interstellar medium (ISM), determining the efficacy of stellar feedback, driving turbulence, and mediating galactic wind launching. Their properties, consequences, and observable manifestations are controlled by the interplay between hydrodynamics, radiative cooling, turbulent mixing at interfaces, and the clustering and temporal overlap of supernova events.

## 1. Phases of Radiative Supernova Bubble Evolution

The dynamical evolution of a supernova bubble spans a sequence of phases, with analytic and simulation results delineating the limits set by radiative losses:

- **Free-Expansion Phase**: The ejecta expand nearly ballistically until the swept-up ISM mass equals the ejecta mass. No significant radiative losses occur.
- **Sedov–Taylor (adiabatic) Phase**: The remnant behaves as a self-similar blastwave, with negligible cooling in the hot interior. The shock radius and velocity evolve as
  \[
  R_{\rm ST}(t) = \xi \left(\frac{E_{\rm SN}}{\rho_0}\right)^{1/5} t^{2/5}, \quad v_{\rm ST}(t) = \frac{2}{5} \frac{R_{\rm ST}}{t}
  \]
  with $\xi \simeq 1.15$ and $\rho_0$ the ambient density [1511.05151], [2404.15396].

- **Pressure-Driven Snowplow (PDS)**: Radiative cooling becomes efficient, particularly at the contact discontinuity (CD) and in a thin, dense shell. The shell is driven outward by the overpressure of the interior, with
  \[
  R_{\rm PDS}(t) \propto t^{2/7},\quad v_{\rm PDS}(t) \propto t^{-5/7}
  \]
  Cooling at this phase dominates energy loss, with significant mixing at the CD controlling the degree of radiative efficiency.

- **Momentum-Conserving Snowplow (MCS)**: When the interior pressure falls to ambient ISM levels, the shell coasts on its accumulated momentum, with
  \[
  R_{\rm MCS}(t) \propto t^{1/4},\quad v_{\rm MCS}(t) \propto t^{-3/4}
  \]

- **Shell Formation and Collapse**: A cold, dense shell forms behind the forward shock once $t_{\rm cool}$ drops below the dynamical time. Compression ratios can exceed 100, producing strong density and magnetic enhancements [2404.15396].

This framework applies to both isolated, rapidly cooling remnants and the more complex hydrodynamics of overlapping events in superbubbles from clustered SNe [1603.00815], [1402.6695].

## 2. Radiative Energy Losses, Mixing, and Retained Feedback

Radiative cooling in supernova bubbles is governed by
\[
\dot E_{\rm rad} = n^2 \Lambda(T) \Delta V
\]
where $n$ is density, $\Lambda(T)$ is the cooling function (typically $10^{-22}$ to $10^{-21}\,\mathrm{erg\,cm^3\,s^{-1}}$ for $10^5$–$10^7$ K gas), and $\Delta V$ the local volume [1511.05151], [2404.15396]. 

- **Mixing Scale Sensitivity:** The magnitude of radiative loss is acutely sensitive to the effective mixing length at the CD—controlled numerically by cell size or physically by turbulent diffusion. Finer resolution or imposed thresholds on cooling (e.g., no cooling below ambient density) can substantially increase the kinetic efficiency ($\epsilon_k$), while explicit conductive/turbulent mixing reduces it by smearing temperature and density gradients.
- **Kinetic Energy Retention:** In high-density ($n\sim100\,\mathrm{cm}^{-3}$) environments with isolated SNe, only $\sim$0.1% of input energy is retained in cold-shell kinetic form. Including pre-existing wind cavities (e.g., from a $60\,M_\odot$ progenitor star with $E_{\rm wind}=2.34 \times 10^{51}$ erg) raises this to $\sim1.5\%$, but still $\gtrsim98\%$ is lost to cooling when $E_{\rm feedback} \sim 3.34 \times 10^{51}$ erg [1511.05151].

- **Superbubble Context:** In clustered scenarios, mechanical efficiency ($\eta_{\rm mech}$) can reach $5$–$10\%$ for ISM density $n_{g0}\sim1\,\mathrm{cm}^{-3}$, with a scaling $\eta_{\rm mech}\propto n_{g0}^{-2/3}$ [1603.00815]. However, most input energy is still radiated, and only larger clusters ($N_{\rm OB}\gtrsim10^4$) sustain overpressured bubbles and steady winds.

| Regime            | Energy Retention           | Dominant Cooling Site              |
|-------------------|---------------------------|------------------------------------|
| Isolated SN       | $\sim$0.1%–1%             | Shell, contact discontinuity       |
| Single SN+Wind    | $\sim$1%–few%             | CD, shell (depends on mixing)      |
| Superbubble, $N>10^4$ | 5–10% (can reach 40% per [1402.6695]) | Shell, CD; losses reduce post-breakout |

## 3. Dynamics of Clustered Supernovae and Superbubble Formation

Superbubbles form as sequential SNe inject energy and mass into a common, hot, low-density bubble. Modeling and simulation reveal:

- **Energy Injection:** For $N_{\rm OB}$ SNe over $\tau_{\rm OB}\sim30$ Myr, the mechanical luminosity is
  \[
  L_w = \frac{N_{\rm OB} E_{\rm SN}}{\tau_{\rm OB}}
  \]
  The bubble's dynamical evolution follows $R(t)\propto t^{3/5}$ in the adiabatic limit, but this is suppressed by $(\eta_{\rm mech})^{1/5}$ due to cooling [1603.00815].

- **Mechanical Efficiency and Density Dependence:** As ambient density increases, mechanical efficiency falls; at $n_{g0}=10\,\mathrm{cm}^{-3}$, $\eta_{\rm mech}\sim 1.1\%$ at $10$ Myr [1603.00815]. Efficiency also rises with numerical resolution but does not converge without explicit physical diffusion.

- **Breakout and Galactic Winds:** For a superbubble to drive a sustained galactic wind (steady Chevalier & Clegg, CC85, solution), a critical $N_{\rm OB}\gtrsim10^4$ is required. Below this, sequential internal shocks dominate rather than a smooth wind [1603.00815], [1402.6695].

- **Fragmentation and Outflows:** Analytic models and observations of HI shells in local galaxies support the finding that most superbubbles stall/fragment within the ISM near or at the gas scale height unless in high-surface-density regions (nuclear starbursts or high-$z$ systems), enabling breakout and wind launching [2109.14626], [1807.08758].

## 4. Morphology, Non-Spherical Features, and Observational Diagnostics

Realistic supernova bubbles exhibit rich, often highly aspherical morphology. Key insights include:

- **Bubble-Like Interiors**: 3D maps of Cas A reveal multiple cavities (e.g., a $\sim$3 ly and a $\sim$1.5 ly bubble), connected to main-shell rings. These voids are the cross-sections of internal bubbles generated by turbulent mixing and radioactive $^{56}$Ni plumes [1501.07283]. 

- **Instabilities:** Rayleigh–Taylor and Kelvin–Helmholtz instabilities, seeded by neutrino-driven convection and SASI, drive large-scale mixing and preserve the "Swiss-cheese" topology of multiple bubbles.

- **Jet-Driven Bubbles:** 3D RHD simulations show that late, opposed jets ($t\sim50$–$100$ d) can carve anisotropic, hot, low-density bubbles, leading to earlier polar breakout of the photosphere and line-of-sight-dependent light curves [2006.01717].

- **Radiative Shell Observability:** The formation of a cold, dense shell at the radiative stage can be traced by nonthermal emission (radio, γ-rays) and, in principle, neutral hydrogen, but HI detection remains challenging [2404.15396]. 

| Morphology           | Diagnostic Method              | Key Observation          |
|----------------------|-------------------------------|--------------------------|
| Bubbles/cavities     | Near-IR [S III], [Fe II] lines| Cas A 3D interior, rings |
| Polar lobes          | Hydrodynamic simulation       | Polar-early, blue excess |
| Dense shell          | Nonthermal emission           | Radio, γ-ray brightening |

## 5. Feedback Consequences and Implications for ISM and Galaxy Evolution

Radiative supernova bubbles, while inefficient as engines for direct mechanical feedback, have substantial implications for structure formation and galaxy evolution:

- **Turbulence and Cloud Disruption:** Even when only $1$–$10\%$ of the input energy is retained, this is sufficient to regulate cloud-scale turbulence and drive GMC dispersal over Myr timescales [1511.05151].
- **Wind Loading and Galaxy Outflows:** The launching of galactic-scale winds depends sensitively on the ability of clustered SNe to reach the ISM scale height before cooling stalls expansion. Breakout leads to reduced radiative losses and efficient venting of energy and metals [1807.08758], [2109.14626].
- **Momentum Partitioning:** Local systems typically retain $<10\%$ of SN feedback momentum in the dense ISM, with the rest driving circumgalactic medium (CGM) flows. High-$z$ galaxies may retain up to $50\%$ [2109.14626].
- **Nonthermal Shell Signatures:** At the radiative stage, the shell enhances local density and magnetic fields, producing a rapid ($\sim100\times$) brightening in radio and γ-rays, detectable with current and upcoming instruments (e.g., CTA). The absence of X-ray emission from old remnants is predicted due to rapid synchrotron cooling [2404.15396].

## 6. Numerical Treatment, Convergence, and Physical Uncertainties

Accurately modeling radiative bubbles necessitates careful treatment of both microphysics and numerical resolution:

- **Mixing Layer Resolution:** The cooling layer at the shell or CD interface is extremely thin, so without explicit diffusion (thermal conduction, viscosity), numerical results for mechanical efficiency do not converge even at $\delta L\sim1$ pc [1603.00815]. This amplifies feedback if under-resolved.
- **Feedback Injection Schemes:** Numerical prescriptions depositing SN energy over large spatial regions (“thermal kernels”) drastically reduce efficiency due to immediate cooling; kinetic schemes or localized (few-pc) injection are more robust [1402.6695]. 
- **Physical Ingredients:** Inclusion of magnetic fields thickens shells and alters radiative loss fractions, while conduction evaporates shell mass, reshaping bubble interiors.
- **Observational Tests:** Multi-wavelength (radio, IR, γ-ray) monitoring can directly test the shell formation paradigm, as well as differentiate radiative bubbles from dense-clump interactions [2404.15396], [1501.07283].

## 7. Scaling Relations, Parameter Dependencies, and Empirical Comparisons

Core equations and scalings appear throughout the literature:

- **Shell Cooling Time and Radius (Cioffi et al. 1988, [2404.15396]):**
  \[
  t_{\rm rad} \approx 1.4 \times 10^4 E_{51}^{3/14} n_0^{-4/7}\, {\rm yr}; \quad R_{\rm rad} \approx 14 E_{51}^{2/7} n_0^{-3/7}\,{\rm pc}
  \]
- **Superbubble Radius Including Radiative Losses ([1603.00815]):**
  \[
  R(t) \simeq 58\,{\rm pc}\, \eta_{\rm mech,-1}^{1/5} \left(\frac{E_{\rm SN,51}\,N_{\rm OB,2}}{\tau_{\rm OB,30} n_{g0}}\right)^{1/5} t_{\rm Myr}^{3/5}
  \]
- **Observed HI Shells:** The break between stalling and breakout shells in the $R$–$v$ plane matches theoretical fragmentation criteria; most observed large shells in local disks are consistent with stalling at the ISM scale height, with breakout restricted to circumnuclear or high-$z$ starburst regions [2109.14626].
- **Energetic and Momentum Loading:** In stratified disc simulations, the post-breakout energy loading is $\eta_E\sim0.2$–$0.6$, with momentum loading $\eta_P\sim1$, and mass loading $\eta_M\sim0.5$–$1$ for escaping hot gas [1807.08758].

Collectively, these results establish radiative supernova bubbles as key, but radiatively inefficient, mediators of ISM structure, turbulence, feedback, and wind driving, with their macroscopic effects controlled by detailed microphysical and environmental parameters.

Source: https://www.emergentmind.com/topics/radiative-supernovae-bubbles