---
title: Planar Shock Breakout Cooling
url: https://www.emergentmind.com/topics/planar-shock-breakout-cooling
type: topic
---

# Planar Shock Breakout Cooling

Searching arXiv for the cited shock breakout papers to ground the article in the literature.
arxiv_search.query({"search_query":"all:\"Non-relativistic radiation mediated shock breakouts\" planar shock breakout cooling Sapir Katz Waxman", "start": 0, "max_results": 10})
arxiv_search.query({"search_query":"all:\"Demystifying shock breakout spectra\" OR all:\"Shock breakout theory\" OR all:\"Shock Breakout in 3-Dimensional Red Supergiant Envelopes\"", "start": 0, "max_results": 10})
arxiv_search.query({"search_query":"id:1103.5276 OR id:1103.5075 OR id:2412.06734 OR id:2206.04134 OR id:1607.01293", "start": 0, "max_results": 10})
Planar shock breakout cooling is the early radiative phase that follows the emergence of a non-relativistic radiation-mediated shock from a stellar surface or, more generally, from the edge of an optically thick medium whose breakout layer is geometrically thin compared with its radius. In this regime the local flow is treated as plane-parallel, breakout occurs when the optical depth ahead of the shock satisfies $\tau \sim c/v_{\rm sh}$, and the post-breakout bolometric luminosity enters a characteristic cooling decline that is asymptotically $t^{-4/3}$ under the standard diffusion, constant-opacity, radiation-dominated assumptions. The modern formulation is anchored by the exact planar breakout solutions of Sapir, Katz, and Waxman and their spherical projection to observable supernova light curves, with subsequent extensions to spectral classification, three-dimensional structure, circumstellar breakout, and opacity-dependent spectral formation [1103.5075] [1103.5276] [1607.01293] [2412.06734].

## 1. Physical definition and regime of validity

Planar shock breakout cooling refers to the phase in which the emitting layer has expanded by much less than the characteristic radius of the system, so curvature is negligible and each emitting patch can be approximated as a slab. In the standard stellar-surface problem, the preshock density near the surface is taken to follow a power law, $\rho \propto x^n$, the flow is non-relativistic, radiation pressure dominates, transport is handled by diffusion with constant opacity, and breakout is triggered once the remaining optical depth is comparable to the radiation-mediated shock width, $\tau \sim c/v_{\rm sh}$. The approximation is accurate when the scale height of the outer envelope is much smaller than the stellar radius and when the planar breakout duration is much shorter than the time on which spherical expansion becomes dynamically important [1103.5075] [1103.5276].

In the spherical projection used for observable light curves, four approximations are central: spherical symmetry in the global geometry, constant Thomson scattering opacity $\kappa$, a steady-state planar angular intensity at the surface, and a small scale-height condition. The surface angular distribution is taken from the steady-state Thomson-scattering solution of Chandrasekhar and is well approximated by
\[
h(\mu)\approx a_I+b_I\mu,\qquad a_I\simeq0.85,\quad b_I\simeq1.725,
\]
with normalization $\int_0^1 h(\mu)\mu\,d\mu=1$. These assumptions are explicitly intended for progenitors whose outer breakout layer is thin compared with $R$ [1103.5276].

A recurring distinction in the literature is between the planar phase and the later spherical cooling-envelope phase. In the planar phase, the breakout shell remains the luminosity-determining shell and the luminosity falls steeply; after the ejecta expands by a factor of order unity, the flow becomes spherical, the luminosity shell moves inward in mass coordinate, and the temporal decline is shallower. The transition is therefore dynamical as well as geometric, rather than merely a change of notation [1004.2496] [1607.01293].

## 2. Universal planar dynamics and bolometric scalings

The exact non-relativistic planar solution is formulated in terms of the breakout shock velocity $v_0$, the preshock density at breakout $\rho_0$, and the opacity $\kappa$. The natural scales are
\[
x_0=\frac{c}{\kappa\rho_0v_0},\qquad
t_0=\frac{x_0}{v_0}=\frac{c}{\kappa\rho_0v_0^2},
\]
and, in the notation of the exact planar problem,
\[
\mathcal{L}_0=\rho_0v_0^3,\qquad
\mathcal{E}_0=\kappa^{-1}\beta_0 c^2,
\]
with $v_0=\beta_0 c$. In these variables the solution is universal for fixed density index $n$: apart from weak $n$-dependence, the dimensionless hydrodynamic and radiative profiles do not depend on the absolute physical scale of the progenitor [1103.5075].

The universal planar luminosity rises before formal breakout because radiation diffuses ahead of the shock, peaks slightly before the Sakurai breakout time, and then relaxes to a cooling tail. The robust asymptotic result is
\[
\mathcal{L}(t)\propto t^{-4/3}\qquad (n>0),
\]
while the special case $n=0$ gives $\mathcal{L}(t)\propto t^{-9/8}$. The light-curve shape depends only weakly on $n$: for $n=1$ to $10$, the luminosity differs by less than about $25\%$ over the interval containing most of the energy. For the commonly used cases, the peak properties are explicitly tabulated: for $n=3$, $\mathcal{L}_{\rm peak}/\mathcal{L}_0=0.72$, $t_{\rm peak}/t_0=-1.25$, and $\Delta t_{\rm FWHM}/t_0=1.61$; for $n=3/2$, $\mathcal{L}_{\rm peak}/\mathcal{L}_0=0.77$, $t_{\rm peak}/t_0=-0.78$, and $\Delta t_{\rm FWHM}/t_0=1.49$ [1103.5075].

A central exact relation links the escaping flux to the acceleration of the surface:
\[
\frac{\partial v(\tau=0,t)}{\partial t}=\frac{\kappa}{c}\mathcal{L}(t),\qquad
v(\tau=0,t)=\frac{\kappa}{c}\int_{-\infty}^{t}\mathcal{L}(t')\,dt' .
\]
As a result, the asymptotic surface velocity is determined by the radiative fluence. In the spherical projection this gives two especially robust outputs:
\[
E_{\rm BO}=8.0\pi R^2\kappa^{-1}cv_0,\qquad
v_{\max}=2.0\,v_0,
\]
both accurate to about $10\%$ and largely insensitive to the detailed density profile, moderate asymmetry in arrival times, or modest changes in the surface angular intensity [1103.5075] [1103.5276].

## 3. Projection to observable light curves

The observable bolometric breakout signal is not the local planar flux itself but the surface emission convolved with both angular intensity and geometric time delay across the stellar disk. In the standard formulation,
\[
L_{\rm obs}(t)=\int_0^1 h(\mu)\,L\!\left(t-\frac{R(1-\mu)}{c}\right)\mu\,d\mu,
\]
where $\mu=\cos\theta$, $h(\mu)$ is the limb-darkened intensity profile, and the delay $R(1-\mu)/c$ accounts for finite light-travel time across the visible hemisphere. With the universal planar light curve written as $L(t)=4\pi R^2\mathcal{L}(t)$ and $\mathcal{L}(t)=\rho_0v_0^3\,\tilde{\mathcal{L}}(t/t_0,n)$, the observable bolometric history is determined by $R$, $v_0$, $\rho_0$, and $\kappa$, plus only a weak dependence on $n$ [1103.5276].

The planar approximation remains valid until spherical expansion appreciably reduces the optical depth of the emitting layers. The transition time is estimated as
\[
t_{\rm sph}\approx \frac{R}{4v_0},
\]
so the exact projected planar light curves are applicable for $t\ll t_{\rm sph}$. A particularly important regime is
\[
t_0=\frac{c}{\kappa\rho_0v_0^2}\ll \frac{R}{c},
\]
which is valid for most progenitors with roughly $R\lesssim10^{14}\,\mathrm{cm}$. In that limit the intrinsic planar pulse is brief, the observed light curve is dominated by light-travel-time smearing, and analytic expressions simplify substantially [1103.5276].

For $4v_0\ll c$, relevant to large red supergiants, the late planar tail persists observationally over
\[
\max(\Delta t_{\rm asym},R/c)\lesssim t\lesssim t_{\rm sph},
\]
with
\[
L\propto t^{-4/3}.
\]
In the same regime the luminosity can be inverted to estimate the progenitor radius,
\[
R\approx 2\times10^{13}
\left(\frac{L}{10^{43}\,\mathrm{erg\,s^{-1}}}\right)^{2/5}
\left(\frac{t}{1\,\mathrm{hr}}\right)^{8/15}\mathrm{cm}.
\]
Moderately asymmetric explosions enter through a spread in shock arrival times, $\Delta t_{\rm asym}$, which modifies the early observed light curve but does not change $E_{\rm BO}$ or $v_{\max}$ [1103.5276].

## 4. Thermalization, spectral formation, and the planar spectrum

The bolometric planar solution does not by itself fix the observed temperature or spectral shape. A separate thermalization problem determines whether the radiation remains in equilibrium. In the analytic treatment of Nakar and Sari, the luminosity shell is defined by $t_d\sim t$, the color shell is where the observed photon energy is set, and the thermal coupling parameter
\[
\eta \equiv \frac{n_{BB}}{\min\{t,t_d\}\,\dot n_{ph,ff}(T_{BB})}
\]
distinguishes equilibrium from non-equilibrium emission. If $\eta<1$, the radiation is thermalized; if $\eta>1$, photon production is insufficient and the observed temperature exceeds the blackbody temperature. During the planar phase, when the breakout shell itself controls the luminosity, the bolometric luminosity declines as
\[
L_{\rm obs}\approx \frac{E_0}{t_0}\left(\frac{t}{t_0}\right)^{-4/3},
\]
while the observed temperature can evolve quite differently depending on whether the breakout shell is initially thermalized [1004.2496].

The spectral classification of the planar phase can be organized by three timescales: the breakout-shell diffusion time
\[
t_{\rm bo}\approx \frac{c}{\kappa \rho_{\rm bo} v_{\rm bo}^2},
\]
the light-crossing time
\[
t_{\rm lc}=\frac{R}{c},
\]
and the thermal-equilibrium reveal time $t_{\rm eq}$. Because $t_{\rm eq}\ge t_{\rm bo}$, there are five allowed orderings and therefore five breakout scenarios. At any given time the spectrum is one of five types: a smeared blackbody, a smeared free-free spectrum, a mixed smeared free-free plus blackbody spectrum, an unsmeared self-absorbed free-free spectrum, or a blackbody. Before $t_{\rm eq}$ the spectrum is a Comptonized free-free spectrum; after $t_{\rm eq}$ a true blackbody component from thermalized material becomes visible. Once $t>\max(t_{\rm bo},t_{\rm lc})$, the flash is over and the later planar phase follows the familiar cooling scalings $T_{\rm obs}\propto t^{-1/3}$ and $L\propto t^{-4/3}$ if thermalized layers dominate [2412.06734].

Planar spectral calculations with local Compton equilibrium and bremsstrahlung photon production show that the surface temperature is determined primarily by $v_0$ and $\rho_0$ and only weakly by $n$. For hydrogen-helium envelopes, the peak surface temperature is fit by
\[
T_{\rm peak}\approx 9.44\exp{[12.63(v_0/c)^{1/2}]}\ {\rm eV},
\]
with more complete fits including a weak density dependence. The time-integrated emitted spectrum is a particularly robust prediction: it depends on $\mathcal{T}_{\rm peak}$ and $v_0$ alone, is relatively insensitive to light-travel-time smearing and slight deviations from spherical symmetry, and peaks at
\[
\nu_{\rm peak}=3\mathcal{T}_{\rm peak}
\]
in the paper’s notation. This robustness is one reason breakout fluence spectra are often treated as cleaner diagnostics than instantaneous spectra [1304.6428].

## 5. Departures from the idealized planar picture

The idealized planar solution is local. Observable breakout signals can depart from the textbook form because different surface patches need not reach breakout simultaneously. In the analytic spherical projection, a spread in shock arrival times $\Delta t_{\rm asym}$ changes the early-time shape but leaves $E_{\rm BO}$ and $v_{\max}$ essentially unchanged. This already implies that the early light curve is more sensitive to angular and temporal smearing than the integrated energetics are [1103.5276].

Three-dimensional radiation-hydrodynamic calculations of red-supergiant envelopes make this point explicit. In the Athena++ calculations, the outer envelope contains a low-density halo and large-scale density fluctuations, so breakout occurs at lower densities than in one-dimensional models and at different radii at different times. Local diffusion times are $t_{\rm diff}\approx1$–$2.2$ hr, but the dominant global timescale is the shock traversal across the corrugated breakout surface,
\[
t_{\rm cross}\approx \frac{\Delta R}{v_{\rm sh}},
\]
with $\Delta R\sim80\,R_\odot$ and $t_{\rm cross}\sim3$–$5$ hr. Measured breakout durations are $\sim3.0$–$4.4$ hr in 3D, versus $\sim0.45$–$0.55$ hr in 1D, and the longer duration lowers the predicted luminosity by a factor of $3$–$10$ to $L_{\rm bol}\sim10^{44}\,\mathrm{erg\ s^{-1}}$. The post-breakout decline of each local patch remains approximately planar, $L_{\rm bol}\propto t^{-4/3}$, but the observed light curve is a superposition of many patches at different stages of breakout and cooling. A common misconception is therefore that the observed rise time directly measures $R/c$; in these 3D models, it does not [2206.04134].

A complementary 3D core-collapse simulation of a $17\,M_\odot$ Type II-P progenitor finds shock breakout at $\sim92{,}000$ s in the southern hemisphere, $\sim118{,}000$ s in the equatorial direction, and $\sim170{,}000$ s in the northern hemisphere, so the angular spread in breakout time is almost a full day. The paper explicitly states that this would smear out the initial breakout flash, but it also explicitly states that the shock-breakout light curves themselves are not calculated before homology is reached. The early post-breakout evolution is instead treated hydrodynamically, with the thermal energy assumed to decline as $1/t$ in nearly adiabatic homologous expansion. This sharpens the distinction between hydrodynamic breakout structure and the dedicated planar radiative-cooling calculation [2411.03434].

Circumstellar material can change the problem qualitatively. If the optical depth of the CSM satisfies $\tau_{\rm CSM}\gtrsim c/V_{\rm bo,*}$, breakout occurs in the CSM rather than at the stellar surface, at a larger radius and with durations that can extend to days. In a dense shell with a sharp outer edge, the early evolution is still planar because the shocked layer is thin, but radiative losses during the planar phase become more important than in the stellar-envelope case. The analytic CSM solution retains the classic $L\propto t^{-4/3}$ limit when the shell is optically thick, but if $\tau_0$ is only marginally above $c/v_{\rm sh}$ the decline steepens and the total luminosity can show an intermediate-time behavior roughly $L\propto t^{-1}$ because multiple diffusion eigenmodes contribute [1607.01293] [2107.04048].

## 6. Opacity revisions, observational diagnostics, and interpretive use

A major revision of the classical spectral picture comes from frequency-dependent opacity with bound-free and bound-bound contributions from heavy elements. In fast Newtonian planar breakouts with $v/c\sim0.2$, adding TOPS opacity rather than assuming a fully ionized free-free medium increases photon production, helps the radiation maintain LTE to higher velocities, and can reduce the emission temperature by half and even an order of magnitude in the planar shock-breakout problem. In the representative case $\beta_{\rm bo}=0.1$ and $\rho_{\rm bo}=10^{-9}\,\mathrm{g\,cm^{-3}}$, the characteristic breakout temperature is around $\sim500$ eV in the free-free-only case and around $\sim150$ eV with TOPS opacity. For red-supergiant envelope breakout, the paper argues that the SED will very likely remain in LTE without stellar wind, implying much weaker X-ray emission than earlier simplified predictions [2509.16996].

Observationally, planar shock breakout cooling is valuable because the early signal depends mainly on a small set of breakout parameters. The rise time of the bolometric luminosity measures $t_{\rm bo}$, the total duration measures $\max(t_{\rm bo},t_{\rm lc})$, the appearance of a blackbody component constrains $t_{\rm eq}$, the total radiated energy constrains $E_{\rm bo}$, and the low-frequency turnover probes the self-absorption frequency. In scenarios with $t_{\rm lc}>t_{\rm bo}$, measuring the rise time, pulse duration, and total radiated energy allows direct estimates of $R$, $\rho_{\rm bo}$, and $v_{\rm bo}$. UV coverage is especially important because that is often where the transition from free-free to thermal emission is most diagnostic; facilities specifically highlighted for this purpose are ULTRASAT, UVEX, and Einstein Probe [2412.06734].

The practical astrophysical role of planar shock breakout cooling is therefore twofold. First, it provides the local building block for realistic breakout calculations in stars, winds, and compact dense shells. Second, it supplies diagnostic relations that are unusually insensitive to uncertain details: the bolometric light curve is weakly dependent on the density profile, the integrated breakout energy and maximum ejecta velocity are robust, and the late planar $t^{-4/3}$ decline can constrain radius when geometric smearing and strong asphericity are under control. At the same time, the literature shows that these same diagnostics can be degraded or reinterpreted by asymmetry, three-dimensional envelope structure, dense CSM, and realistic high-ionization opacity. Planar shock breakout cooling is thus best understood not as a complete description of early supernova light, but as the universal local regime from which more global and more realistic breakout phenomenology is constructed [1103.5075] [1103.5276] [2206.04134] [2509.16996].

Source: https://www.emergentmind.com/topics/planar-shock-breakout-cooling