Expansion and photon diffusion reshape supernova energy
Relate an illustrative ejecta energy to a diffusion light-curve scale without equating luminosity and explosion power.
On this page
The question for today
What you’ll understand
- Distinguish ejecta kinetic energy from internal radiation energy and heating.
- Estimate an expanding-ejecta diffusion timescale with stated geometric conventions.
- Explain why a light-curve peak need not mark optical depth one.
A supernova shock can deposit substantial internal energy and accelerate ejecta, but a later optical light curve does not simply report that initial shock power. Expansion does work, photons require time to escape, and radioactive or central-engine heating can add energy after the explosion. A transparent one-zone calculation separates these budgets and produces a useful timescale. Its numerical coefficient is a model estimate, rather than a precision measurement of ejecta mass.
Specify homologous expansion and separate energy reservoirs
Consider spherical ejecta of fixed mass M with outer radius R≈vt, once the initial stellar radius is negligible. Homologous expansion means each fluid element maintains approximately constant velocity proportional to radius. For a uniform-density homologous sphere, v is the outer-edge speed and the kinetic energy is 3Mv²/10, obtained by integrating one half velocity squared over mass. Using Mv²/2 instead would describe a different single-speed shell approximation. Declaring the velocity convention prevents an apparently small coefficient from silently changing inferred energy.
Let Eᵢₙₜ be the ejecta’s trapped radiation energy. In a radiation-dominated one-zone model its pressure is Eᵢₙₜ/(3V), so the first law includes expansion work, deposited heating and escaping bolometric luminosity. Heating means energy actually thermalized in the ejecta; radioactive decay energy escaping as gamma rays is not automatically deposited. A shock’s initial internal energy can become kinetic energy through pressure work, and at late freely expanding times the remaining internal energy can be much smaller than the kinetic reservoir.
Diffusion competes with the decreasing expansion time
Use a constant grey opacity κₒₚ in area per mass and uniform density as an illustrative closure. The center-to-edge optical depth is κₒₚρR. A random-walk escape estimate is that optical depth times R/c; its coefficient depends on the density profile and boundary treatment. Because radius grows while density falls as radius cubed, the instantaneous diffusion estimate decreases as one over time. Equating it to the expansion age defines a characteristic light-curve width, with a square-root dependence on mass and opacity.
At that crossing, the optical depth is approximately c/v, not one. For nonrelativistic ejecta this can still be much larger than unity: photons escape efficiently relative to expansion while the ejecta remain optically thick. The actual time of maximum luminosity also depends on the heating history, ionization, recombination and density profile. The familiar statement that luminosity at peak tracks instantaneous deposited heating belongs to particular diffusion models and is approximate for real events. It is not a general identity equating peak radiation with initial kinetic energy.
A spectrum requires more than a diffusion clock
A bolometric luminosity describes total radiated power, whereas an observed band samples a spectrum shaped by wavelength-dependent opacity and thermalization. A blackbody photospheric estimate can illustrate colors but does not reproduce a line-rich supernova spectrum. Recombination can alter the effective opacity and produce a plateau, while circumstellar interaction can supply continuing shock power. Radioactive heating, magnetar heating and interaction may yield partially degenerate optical curves. Distinguishing them requires velocities, spectra, temporal coverage and consistent deposited-energy accounting, rather than one fitted width.
The heat lab isolates diffusion at fixed geometry. Reset and set both controls to mode one and diffusivity 0.25; Step twice to time one, predicting amplitude exp(−0.25)≈0.779. Reset and set mode one and diffusivity one; Step twice, predicting exp(−1)≈0.368. The instrument’s κ is a diffusivity, not the opacity κₒₚ: physical radiative diffusivity scales inversely with opacity at fixed density. This periodic mode has no expansion, radioactive heating, photosphere or emergent supernova light curve.
Work it through
Estimate ejecta energy and a diffusion width
Take hypothetical uniform-density homologous ejecta with M=2×10³⁰ kg, outer speed v=10⁷ m/s and κₒₚ=0.01 m²/kg. Use c=3×10⁸ m/s. Find Eₖ, the characteristic t* and the optical depth at t*.
The uniform homologous kinetic energy is (3/10)Mv²=0.3×2×10³⁰×10¹⁴=6×10⁴³ J. The outer-speed convention determines the coefficient.
The estimated timescale squared is 3κₒₚM/(4πcv)=6×10²⁸/(12π×10¹⁵)≈1.59×10¹² s². Hence t*≈1.26×10⁶ s, or 14.6 days.
At this timescale τ≈c/v=30. The radius is approximately 1.26×10¹³ m, and τR/c≈1.26×10⁶ s, checking the crossing condition.
The illustrative model gives Eₖ=6×10⁴³ J and t*≈14.6 days while τ≈30; the width is not the time of optical thinning.
The living notebook
Explore the model
Interactive lab · heat
Keep diffusivity distinct from opacity
Reset and set BOTH controls n=1,κ=0.25; Step twice to t=1: predict amplitude ≈0.779. Reset and set n=1,κ=1; Step twice: predict ≈0.368. This κ is diffusivity, whereas radiative diffusivity decreases as opacity rises at fixed density. No expanding ejecta or heating is modeled.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA diffusion instrument illustrates transport but does not supply a supernova opacity or peak luminosity.
Pause & reason
Check your understanding
Your turn
Try it yourself
A radiation-dominated homologously expanding parcel has trapped internal energy 10⁴³ J at radius R₀. Over a later interval neglect both heating and radiative escape and let the radius grow to 10R₀. Find the remaining internal energy and the energy transferred through expansion work.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
With only adiabatic radiation work, Eᵢₙₜ scales as V⁻¹ᐟ³ and hence as R⁻¹.
Compare with the solution
The internal energy falls by a factor of ten to 10⁴² J. Expansion work removes 10⁴³−10⁴²=9×10⁴² J from that reservoir. No energy has been radiated under the stated zero-escape assumption. This limiting calculation explains why the initial shock heat is not all available for late luminosity; maintaining exactly free expansion assumes the added kinetic energy is dynamically negligible at the stage considered.
Keep the idea
Take it with you
- Kinetic energy, trapped radiation and deposited heating are different budgets.
- An expanding diffusion width depends on geometry, opacity and velocity conventions.
- Efficient photon escape relative to expansion can occur at optical depth much greater than one.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.