Radiative Supernova Bubbles

• Radiative supernova bubbles are regions of shock-heated, multi-phase gas created by massive star explosions with pronounced radiative cooling.
• They trace key evolutionary phases—from free expansion to momentum-conserving stages—that influence energy retention and feedback efficiency in the ISM.
• Detailed simulations and analytic scalings elucidate dependencies on mixing, density, and energy injection schemes critical for modeling ISM dynamics.

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 (Fierlinger et al., 2015, Diesing et al., 2024).

• 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 (Diesing et al., 2024).

This framework applies to both isolated, rapidly cooling remnants and the more complex hydrodynamics of overlapping events in superbubbles from clustered SNe (Yadav et al., 2016, Sharma et al., 2014).

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 (Fierlinger et al., 2015, Diesing et al., 2024).

• 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 (Fierlinger et al., 2015).
• 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}$ (Yadav et al., 2016). 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 (Sharma et al., 2014))	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 (Yadav et al., 2016).

• 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 (Yadav et al., 2016). 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 (Yadav et al., 2016, Sharma et al., 2014).
• 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 (Orr et al., 2021, Fielding et al., 2018).

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 (Milisavljevic et al., 2015).
• 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 (Akashi et al., 2020).
• 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 (Diesing et al., 2024).

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 (Fierlinger et al., 2015).
• 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 (Fielding et al., 2018, Orr et al., 2021).
• 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\%$ (Orr et al., 2021).
• 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 (Diesing et al., 2024).

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 (Yadav et al., 2016). 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 (Sharma et al., 2014).
• 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 (Diesing et al., 2024, Milisavljevic et al., 2015).

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, (Diesing et al., 2024)):

$$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 (Yadav et al., 2016):

$$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 (Orr et al., 2021).
• 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 (Fielding et al., 2018).

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.