Wind Breakout in Supernovae and Stellar Winds

• Wind Breakout is the phenomenon where a supernova shock emerges through an optically thick circumstellar wind rather than the stellar surface, delaying and reshaping the breakout flash.
• Analytic models, radiation-hydrodynamics simulations, and Monte Carlo calculations are used to capture shock structure, diffusion timescales, and spectral evolution in the wind breakout process.
• Applications include diagnosing recent mass-loss episodes, probing circumstellar density profiles, and extending the concept to galactic superbubbles and compact-object merger ejecta.

Wind breakout most commonly denotes supernova shock breakout through an optically thick circumstellar wind rather than at the stellar photosphere. In that setting, a radiation‑mediated shock launched in the progenitor envelope continues into a wind with $\rho(r)\propto r^{-2}$, and the observable breakout occurs where the optical depth ahead of the shock falls to order $c/v_{\rm sh}$, delaying and reshaping the emergent flash relative to bare‑star breakout (Balberg et al., 2011, Svirski et al., 2014). In current literature, the term also appears in related but distinct contexts: relativistic breakout from Wolf–Rayet winds, early collisionless shocks and neutrino precursors, superbubble breakout from galactic discs into large‑scale winds, and magnetar‑wind‑driven shock breakout from merger ejecta (Giacinti et al., 2015, Granot et al., 2017, Roy et al., 2013, Wu et al., 2022).

1. Definition, geometry, and breakout condition

In the canonical stellar‑wind formulation, the circumstellar medium is a steady spherical wind,

$\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$

so the line‑of‑sight optical depth is

$$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$

Shock breakout in such a wind occurs when the optical depth ahead of the shock satisfies

$$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$

which gives a breakout radius

$$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$

This is the basic distinction between stellar‑surface breakout and wind breakout: the relevant emitting radius is $r_{\rm br}>r_\ast$, not the hydrostatic stellar radius (Giacinti et al., 2015).

For observed X‑ray flashes, this changes the characteristic timescale from the light‑crossing time $R_\ast/c$ to the dynamical time $R_{\rm bo}/v_{\rm bo}$. In the thick‑wind interpretation of SN 2008D, for example, the breakout condition is satisfied at $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$, whereas the adopted Wolf–Rayet radius is $c/v_{\rm sh}$0; breakout therefore occurs slightly above the stellar surface but still in the wind, and the rise time is set by $c/v_{\rm sh}$1 s rather than $c/v_{\rm sh}$2 s (Svirski et al., 2014).

Balberg & Loeb formulated the same problem in terms of observables of the breakout flash. For a wind $c/v_{\rm sh}$3, the diffusion time through the unshocked wind is

$c/v_{\rm sh}$4

and breakout is identified by $c/v_{\rm sh}$5. With total X‑ray energy $c/v_{\rm sh}$6 and diffusion time $c/v_{\rm sh}$7, one can infer a breakout radius

$c/v_{\rm sh}$8

which makes wind breakout a directly diagnostic probe of the dense inner circumstellar medium (Balberg et al., 2011).

2. Shock structure: radiation mediation, collisionless transition, and early shock formation

The standard picture assumes that a radiation‑dominated shock (RDS) propagates through the star and survives until breakout, where a collisionless shock (CS) forms. Giacinti and Bell showed that this picture is incomplete for moderately optically thick winds. Their central argument is that the RDS may fail soon after leaving the stellar core because the radiation field dilutes as $c/v_{\rm sh}$9 and cannot keep successive upstream shells co‑moving. For a steady $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$0 wind, they obtain an approximate condition

$\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$1

under which a collisionless shock forms deep inside the wind, at $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$2, well before classical breakout (Giacinti et al., 2015).

For red‑supergiant‑like parameters with $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$3 cm, $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$4, $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$5, and $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$6, Giacinti and Bell find $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$7, but the CS appears already at $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$8. Their 1D radiation‑hydrodynamics simulations show a sharp velocity and temperature discontinuity at $\rho(r)=\frac{\dot M}{4\pi r^2 u_{\rm w}},$9 embedded in a decaying radiation‑mediated structure; the CS persists and strengthens as it propagates outward (Giacinti et al., 2015).

This revision has immediate non‑thermal consequences. Once the CS exists, diffusive shock acceleration can begin before photons escape. Assuming Bohm diffusion and a surface magnetic field $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$0 G, the proton acceleration time is

$$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$1

shorter than the time available before breakout. The resulting $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$2 and $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$3 interactions produce TeV neutrinos that may precede the first breakout photons by

$$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$4

for red supergiants with $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$5 and $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$6, and $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$7 minutes for similar Wolf–Rayet progenitors (Giacinti et al., 2015).

The same work also delineates the opposite regime: in very thick winds, such as Type IIn systems with $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$8, the inequality above cannot be satisfied with realistic $$\tau(r)=\kappa\int_r^\infty \rho(r')dr'=\kappa \frac{\dot M}{4\pi u_{\rm w}}\frac{1}{r}.$$9, so the RDS survives to the conventional breakout depth and the classical sequence remains applicable (Giacinti et al., 2015).

3. Radiation physics, spectra, and temporal evolution

A major development in the theory of wind breakout is the recognition that the post‑breakout collisionless shock and the residual soft‑photon bath deposited by the earlier radiation‑mediated phase can be tightly coupled. Svirski and Nakar showed that in Wolf–Rayet explosions through thick winds, hot shocked electrons with $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$0 keV are cooled efficiently by inverse Compton scattering of soft photons injected from the stellar envelope. They parameterize the soft‑photon luminosity as

$$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$1

and find $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$2 for Wolf–Rayet stars. When

$$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$3

multiple inverse Compton scatterings produce a nearly flat spectrum,

$$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$4

between $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$5 and $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$6 (Svirski et al., 2014).

This framework naturally reproduces the unusually flat early X‑ray spectrum of SN 2008D. Using $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$7 s and $$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$8, Svirski and Nakar infer

$$\tau(r_{\rm br})\sim \frac{c}{u_{\rm s}}=\beta_{\rm s}^{-1},$$9

that is,

$$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$0

with $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$1 and $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$2. In their model the bolometric luminosity decays as

$$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$3

while $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$4, and the X‑ray band can steepen to $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$5 once the wind becomes optically thin enough that multiple inverse Compton scatterings control the X‑ray production (Svirski et al., 2014).

A complementary semi‑analytic treatment of fast, sub‑relativistic radiation‑mediated shocks predicts that when $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$6, photon escape steepens the shock, increases compression, and boosts downstream photon generation, so the characteristic breakout temperature drops by roughly an order of magnitude while the flux is still rising. In that picture the time‑integrated breakout spectrum is broad and non‑thermal even though the instantaneous emission is locally thermal, and the system is constrained by a closure among $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$7, $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$8, and $$r_{\rm br}\approx \frac{\kappa \dot M \beta_{\rm s}}{4\pi u_{\rm w}}.$$9 (Ioka et al., 2018).

Monte‑Carlo calculations of finite radiation‑mediated shocks with radiative losses sharpen this picture further. Ito, Levinson, and Nakar find that the $r_{\rm br}>r_\ast$0 peak depends strongly on shock velocity, from $r_{\rm br}>r_\ast$1 keV for $r_{\rm br}>r_\ast$2 to $r_{\rm br}>r_\ast$3 keV for $r_{\rm br}>r_\ast$4, and shifts to lower energies as losses increase. They also find a nearly flat low‑energy component,

$r_{\rm br}>r_\ast$5

below the peak, extending down to the absorption break and implying much brighter optical/UV emission than expected from a pure Wien extrapolation (Ito et al., 2020).

4. Observational exemplars and diagnostic use

SN 2008D/XRT 080109 is the best‑studied candidate for shock breakout through a thick Wolf–Rayet wind. In the Svirski–Nakar interpretation, its $r_{\rm br}>r_\ast$6 s rise, $r_{\rm br}>r_\ast$7 s luminous phase, flat $r_{\rm br}>r_\ast$8–$r_{\rm br}>r_\ast$9 keV spectrum, and later steepening are explained by breakout at $R_\ast/c$0 in a wind with $R_\ast/c$1, requiring $R_\ast/c$2 at small radii and implying an enhanced mass‑loss episode beginning approximately $R_\ast/c$3 days or less before core collapse (Svirski et al., 2014).

A different Monte‑Carlo treatment, which includes bulk Comptonization and aspherical shock geometry but assumes a steady spherical wind, reproduces the observed light curve and spectrum with a spherical reference velocity $R_\ast/c$4, breakout radius $R_\ast/c$5 cm, and wind normalization $R_\ast/c$6, corresponding to

$R_\ast/c$7

for $R_\ast/c$8. In that model, even an aspherical explosion breaks out simultaneously in all directions in a steady spherical $R_\ast/c$9 wind because $R_{\rm bo}/v_{\rm bo}$0 and therefore $R_{\rm bo}/v_{\rm bo}$1 is angle‑independent (Ohtani et al., 2018).

PTF 09uj provides a slower, ultraviolet‑selected example. Its GALEX light curve rose on a timescale of a few days to $R_{\rm bo}/v_{\rm bo}$2, and analytic modeling of a dense $R_{\rm bo}/v_{\rm bo}$3 circumstellar wind gives

$R_{\rm bo}/v_{\rm bo}$4

$R_{\rm bo}/v_{\rm bo}$5

for $R_{\rm bo}/v_{\rm bo}$6. The inferred dense wind must have been ejected over a few years, making PTF 09uj a prototype of shock breakout in a compact, optically thick circumstellar outflow rather than at the stellar photosphere (Ofek et al., 2010).

These examples illustrate a recurring diagnostic logic. Light‑curve timescales primarily constrain $R_{\rm bo}/v_{\rm bo}$7 and $R_{\rm bo}/v_{\rm bo}$8; spectral hardness constrains $R_{\rm bo}/v_{\rm bo}$9; and the comparison between inner and outer wind densities can reveal recent, enhanced pre‑explosion mass loss (Svirski et al., 2014, Ofek et al., 2010).

5. Relativistic, aspherical, and multi‑messenger extensions

Wind breakout is not restricted to sub‑relativistic shocks. Levinson and Nakar constructed an analytic model for relativistic radiation‑mediated shock breakout from a stellar wind in which pair creation self‑adjusts the opacity. In that regime, breakout occurs when the total optical depth ahead of the shock is

$R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$0

provided the local shock Lorentz factor is still $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$1. The observable scalings are then

$R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$2

and

$R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$3

with the closure relation

$R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$4

For typical parameters they find about $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$5 ergs radiated in MeV gamma‑rays over durations from a fraction of a second to an hour, potentially detectable to $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$6–$R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$7 Mpc (Granot et al., 2017).

At lower relativistic or fast‑Newtonian speeds, Monte‑Carlo simulations of finite radiation‑mediated shocks support a more nuanced statement about spectral formation. The emission around the $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$8 peak is close to a Wien spectrum, while the broad low‑energy component below the peak is nearly flat in $R_{\rm bo}\approx 6\times 10^{11}\,\mathrm{cm}$9. A plausible implication is that some apparently power‑law time‑integrated spectra mix distinct phases: a radiation‑mediated breakout component during the rise and a later collisionless‑shock component during the decay (Ito et al., 2020).

The multi‑messenger extension arises from the early collisionless‑shock scenario discussed above. If a CS forms well before breakout in a moderately thick wind, TeV neutrinos can escape while the radiation field remains trapped, providing a precursor signal that directly probes the otherwise invisible wind interior (Giacinti et al., 2015).

6. Broader astrophysical usages of the term

Outside stellar wind breakout from supernovae, the phrase also denotes breakout of superbubbles from galactic discs. In a stratified ISM, superbubbles driven by clustered supernovae break out when they reach the gas scale height and remain supersonic, allowing hot interior gas to vent into the halo. Roy et al. identified two thresholds in energy injection surface density: $c/v_{\rm sh}$00 for gentle breakout with Mach number $c/v_{\rm sh}$01–$c/v_{\rm sh}$02, and

$c/v_{\rm sh}$03

for vigorous breakout with Mach number $c/v_{\rm sh}$04–$c/v_{\rm sh}$05, corresponding to $c/v_{\rm sh}$06. They also found that thermal instability, not Rayleigh–Taylor instability, first fragments the shell for typical disc parameters (Roy et al., 2013).

Fielding and collaborators extended this clustered‑supernova picture with 3D stratified simulations. They find that superbubbles from clustered supernovae can break out well before the cluster’s supernovae cease, that post‑breakout radiative losses are reduced dramatically, and that a large fraction,

$c/v_{\rm sh}$07

of the injected supernova energy vents into the halo. Their winds can reach mass outflow rates of order the star‑formation rate, while the momentum flux remains only of order the direct supernova momentum input because hot gas vents before doing substantial work on the surroundings. They further argue that the onset of cluster‑driven superbubble breakout corresponds to a threshold $c/v_{\rm sh}$08 (Fielding et al., 2018).

A different extension appears in compact‑object merger theory. If a double‑neutron‑star merger leaves a long‑lived magnetar, its relativistic wind can drive a shock through anisotropic merger ejecta and produce ultraviolet/soft‑X‑ray shock breakout emission. In this context, breakout occurs at different times in different directions because the ejecta structure is strongly anisotropic, so the bolometric SBO light curve is broadened and can, in principle, be used to probe the angular structure of the ejecta (Wu et al., 2022).

Wind breakout therefore denotes not a single mechanism but a family of breakout phenomena in which an embedded shock emerges from an optically thick outflow, atmosphere, or shell. In supernovae, the dominant contemporary meaning remains shock breakout through a dense circumstellar wind; yet the same underlying idea—radiation trapped behind a shock and released only after escape conditions are satisfied—connects stellar explosions, galactic superbubbles, and compact‑merger ejecta (Balberg et al., 2011, Roy et al., 2013, Wu et al., 2022).