Explosive Nucleosynthesis

An iron core cannot burn — fusion past the iron peak costs energy rather than releasing it — so a massive star that has climbed through hydrogen, helium, carbon, neon, oxygen and silicon burning collapses. When the collapse rebounds, or when a carbon–oxygen white dwarf ignites a runaway, a shock tears outward through layers sitting in relatively cool hydrostatic equilibrium. In the shocked gas the temperature spikes to several billion kelvin for barely a second or two before adiabatic expansion cools it back down. In that brief window nuclear reactions run explosively — hence explosive nucleosynthesis — and then freeze, so the composition that flies away is not an equilibrium state but a frozen photograph of an interrupted march toward one. The outcome is decided by a race between two clocks: the nuclear reaction time and the hydrodynamic expansion time. Where reactions win, the gas reaches nuclear statistical equilibrium (NSE) and its abundances depend only on temperature, density and the electron fraction, collapsing onto the most tightly bound nuclei. Where expansion wins, the pre-shock composition is quenched unchanged. Because peak temperature falls with radius, the ejecta are sorted into concentric zones — complete silicon burning to iron-peak nuclei near the center, then incomplete silicon burning, explosive oxygen/neon burning, and unprocessed outer layers. The signature product is radioactive $^{56}\mathrm{Ni}$ , whose decay to $^{56}\mathrm{Fe}$ lights the supernova for months and whose mass sets the peak brightness — a clean thread from nuclear physics to the observed sky.

Key takeaways

Explosive nucleosynthesis is the burning when a supernova shock or thermonuclear runaway flashes stellar matter to $T\gtrsim 2\times10^{9}\,\mathrm{K}$ for seconds — far faster than hydrostatic.
Yields are set by a race between the nuclear reaction time $\tau_{\mathrm{nuc}}$ and the hydrodynamic time $\tau_{\mathrm{HD}}\simeq 446/\sqrt{\rho}\,$ s; the shorter clock wins.
Peak temperature falls with radius: complete Si-burning ( $T_9\gtrsim5$ ) makes $^{56}\mathrm{Ni}$ and the iron peak; incomplete Si-burning and explosive O/Ne/C-burning make the $\alpha$ -elements Si, S, Ar, Ca.
The hottest zones reach nuclear statistical equilibrium (NSE), a Saha-like state fixed only by $T$ , $\rho$ and $Y_e$ , favouring doubly magic $^{56}\mathrm{Ni}$ ( $Z=N=28$ ) at $Y_e\approx0.5$ .
The chain $^{56}\mathrm{Ni}\to{}^{56}\mathrm{Co}\to{}^{56}\mathrm{Fe}$ (half-lives 6.08 d and 77.2 d) powers the light curve; by Arnett's rule the peak luminosity is proportional to the $^{56}\mathrm{Ni}$ mass.
Key idea: explosive burning is a frozen photograph of an interrupted march toward nuclear statistical equilibrium — what survives is set by whether each reaction beats the seconds-long expansion, and the peak temperature at each radius selects the zone (complete Si-burning to $^{56}\mathrm{Ni}$ , incomplete Si-burning, or explosive O/Ne/C-burning) that sets the local yield.

The core insight

Explosive burning is a frozen photograph of an interrupted march toward nuclear statistical equilibrium: what survives is set by whether each reaction's timescale beats the seconds-long expansion, and the peak temperature reached at each radius selects which zone — complete Si-burning to $^{56}\mathrm{Ni}$ , incomplete Si-burning, or explosive O/Ne/C-burning — produces the local yield.

Where this lives today

Type Ia supernovae, whose $^{56}\mathrm{Ni}$ -set peak luminosities make them standardizable candles, are the distance-ladder rungs that revealed dark energy — and the puzzle of which white-dwarf systems explode is attacked with early light curves and nebular spectroscopy that read off the nickel distribution. Gamma-ray astronomy watches the decay directly: NuSTAR mapped the $^{44}\mathrm{Ti}$ from $\alpha$ -rich freeze-out across Cassiopeia A and SN 1987A, revealing asymmetric explosions, while INTEGRAL caught the $^{56}\mathrm{Co}$ 847 and 1238 keV lines from the Type Ia SN 2014J. Chemical-evolution surveys (APOGEE, GALAH, Gaia) trace the decline of $[\alpha/\mathrm{Fe}]$ with $[\mathrm{Fe}/\mathrm{H}]$ that encodes delayed Type Ia iron against prompt core-collapse $\alpha$ -elements, while short-lived radionuclides ( $^{44}\mathrm{Ti}$ , $^{60}\mathrm{Fe}$ ) in meteorites let cosmochemists test burning-zone yields atom by atom. The open frontier is the core-collapse mechanism itself — three-dimensional simulations still struggle to reproduce the observed $\sim0.05$ – $0.1\,M_\odot$ of ejected nickel.

The physical picture

Picture a wall of heat — the supernova shock — sweeping outward through the layered star, flashing the gas at each depth to a peak temperature for a second or two before expansion quenches it. The hotter the flash, the more completely the local nuclei rearrange toward the most bound configuration: the deepest, hottest shells relax all the way to iron-group nuclei, shallower shells reach only silicon-to-calcium, and the outermost shells barely change. One explosion thus writes a whole periodic-table's worth of products, indexed by peak temperature. The first figure is that stratification.

Figure

Peak shock temperature sorts the ejecta into concentric burning zones: complete silicon burning at T₉ ≳ 5 reaches nuclear statistical equilibrium and makes iron-peak ⁵⁶Ni, incomplete silicon burning at 4 ≲ T₉ ≲ 5 and explosive oxygen/neon burning at 3 ≲ T₉ ≲ 4 make silicon through calcium, and below T₉ ≈ 2 the outer layers are ejected unprocessed

Think of NSE like a chemical equilibrium: at the highest temperatures photodisintegration shatters nuclei into free protons, neutrons and $\alpha$ particles, and as the gas cools these reassemble into the most stable available nucleus — with equal protons and neutrons, $^{56}\mathrm{Ni}$ . But if cooling outruns the slow triple-alpha reassembly, free helium is stranded in the $\alpha$ -rich freeze-out. The event is radioactive at heart: weak decays are far too slow to build stable $^{56}\mathrm{Fe}$ directly, so the explosion makes $^{56}\mathrm{Ni}$ and lets it decay later, converting binding energy into the glow that unfolds over weeks. The second figure is that glow.

Figure

The radioactive light curve powered by ⁵⁶Ni → ⁵⁶Co → ⁵⁶Fe: full-trapping energy deposition falls steeply while ⁵⁶Ni with half-life 6.08 d dominates, then settles onto the ⁵⁶Co tail with half-life 77.2 d whose logarithmic slope is 0.0098 mag per day

The core equations

The clock the reactions must beat is the hydrodynamic expansion time. For shocked, self-gravitating matter this is essentially the free-fall time; the density-scaled form below is the conventional e-folding time of the density during expansion, the yardstick for every reaction.

\tau_{\mathrm{HD}}\simeq \frac{446}{\sqrt{\rho\,[\mathrm{g\,cm^{-3}}]}}\ \mathrm{s} \;=\; \frac{1}{\sqrt{24\pi G\rho}} \qquad (1)

where: $\tau_{\mathrm{HD}}$ — expansion/cooling time of the shocked matter; $\rho$ — local mass density; $G$ — Newton's constant. Inserting $G=6.674\times10^{-8}\,\mathrm{cgs}$ into $(24\pi G\rho)^{-1/2}$ gives $445.8\approx446$ with $\rho$ in $\mathrm{g\,cm^{-3}}$ . Reading it: at $\rho\sim10^{7}\,\mathrm{g\,cm^{-3}}$ the hot conditions last only $\approx0.14\,\mathrm{s}$ , so a reaction must beat a tenth of a second to matter. Each nucleus's fate is set by the ratio of its reaction time to that expansion time. Because the thermally averaged rate is ferociously temperature-sensitive (the Gamow peak), the ratio flips from "wins" to "loses" over a narrow temperature band.

\tau_{\mathrm{nuc}}=\frac{1}{n_j\langle\sigma v\rangle}\qquad\text{compared with}\qquad \tau_{\mathrm{HD}} \qquad (2)

where: $\tau_{\mathrm{nuc}}$ — timescale for a nuclear reaction; $n_j$ — number density of the reacting partner; $\langle\sigma v\rangle$ — thermally averaged cross section times relative velocity. Reading it: $\tau_{\mathrm{nuc}}\ll\tau_{\mathrm{HD}}$ gives equilibrium (NSE or quasi-equilibrium); $\tau_{\mathrm{nuc}}\gg\tau_{\mathrm{HD}}$ quenches the reaction and freezes the composition. To follow a fluid element as it cools we need its trajectory. Radiation pressure dominates the shocked gas, so the adiabat $P\propto T^4$ , $P\propto\rho^{4/3}$ fixes how temperature tracks density and gives the exponential cooling used for freeze-out.

T\propto\rho^{1/3}\quad\Longrightarrow\quad \rho(t)=\rho_0\,e^{-t/\tau},\ \ T(t)=T_0\,e^{-t/3\tau} \qquad (3)

where: $T,\rho$ — temperature and density after shock passage; $\tau\sim\tau_{\mathrm{HD}}$ — expansion timescale; $\rho_0,T_0$ — immediate post-shock values. Reading it: temperature falls three times more slowly than density, the standard Fowler–Hoyle prescription for tracking freeze-out. In the zones that reach equilibrium, the composition follows a Saha-like law: balancing each reaction against its inverse forces every nucleus's chemical potential to a sum of proton and neutron potentials, and the abundances collapse onto a function of three variables.

Y_i=\frac{G_i}{2^{A_i}}\,A_i^{3/2}\,(\rho N_A)^{A_i-1}\left(\frac{2\pi\hbar^2}{m_u k_B T}\right)^{\frac{3}{2}(A_i-1)}Y_p^{Z_i}Y_n^{N_i}\,e^{B_i/k_B T} \qquad (4)

where: $Y_i$ — molar abundance of nucleus $i$ ( $A_i,Z_i,N_i$ ); $G_i$ — partition function; $B_i$ — binding energy; $Y_p,Y_n$ — free proton and neutron molar abundances; $N_A$ — Avogadro's number; $m_u$ — atomic mass unit. Reading it: the factor $e^{B_i/k_B T}$ rewards binding energy, so at fixed $T,\rho,Y_e$ the equilibrium piles up on the most bound accessible nucleus. Which iron-group nucleus wins depends not on temperature but on the neutron-to-proton balance, encoded in the electron fraction — conserved because the weak interactions that change it are frozen out.

Y_e=\sum_i \frac{Z_i}{A_i}X_i \qquad (5)

where: $Y_e$ — electrons (equivalently protons) per baryon; $X_i$ — mass fraction of species $i$ ; $Z_i,A_i$ — its charge and mass number. Reading it: $Y_e=0.5$ selects symmetric $Z=N$ $^{56}\mathrm{Ni}$ ; a small excess ( $Y_e<0.5$ ) shifts the peak to $^{54}\mathrm{Fe}$ and $^{58}\mathrm{Ni}$ . The premier product does not survive as itself: it decays in two steps that are the engine of the light curve.

^{56}\mathrm{Ni}\ \xrightarrow{\;6.08\,\mathrm{d}\;}\ {}^{56}\mathrm{Co}\ \xrightarrow{\;77.2\,\mathrm{d}\;}\ {}^{56}\mathrm{Fe} \qquad (6)

where: $^{56}\mathrm{Ni}$ — doubly magic even–even nucleus ( $Z=N=28$ ), dominant NSE product at $Y_e=0.5$ ; half-lives 6.08 d and 77.2 d belong to $^{56}\mathrm{Ni}$ and $^{56}\mathrm{Co}$ . Reading it: the two decays deposit gamma rays and positrons that thermalize, and the observed luminosity tracks that deposition. Under full trapping the luminosity equals the decay heating summed over the ejecta — turning a light curve into nuclear physics.

L(t)=\int\!\big(\lambda_{\mathrm{Ni}}N_{\mathrm{Ni}}Q_{\mathrm{Ni}}+\lambda_{\mathrm{Co}}N_{\mathrm{Co}}Q_{\mathrm{Co}}\big)\,\mathrm{d}V \qquad (7)

where: $L(t)$ — luminosity powered by decay (full trapping); $\lambda=\ln 2/t_{1/2}$ — decay constant; $N$ — number of radioactive nuclei; $Q$ — energy deposited per decay. Reading it: since $N\propto$ the nickel mass, so is the luminosity — the seed of Arnett's rule.

Derivation AWhy nickel-56 dominates complete silicon burning

The goal is to show, straight from the NSE Saha law of equation (4), that a mass shell reaching complete silicon burning at $Y_e\approx0.5$ piles its equilibrium composition onto $^{56}\mathrm{Ni}$ rather than any other iron-group nucleus or the ultimately more bound $^{56}\mathrm{Fe}$ . The route is to isolate the two composition levers — the binding-energy exponential and the free-nucleon abundances raised to the proton and neutron powers — then follow the temperature down until the equilibrium re-forms and freezes.

Step 1. Strip equation (4) to its two composition-dependent factors, $Y_i\propto Y_p^{Z_i}Y_n^{N_i}\,e^{B_i/k_B T}$ . Why: at fixed $T$ , $\rho$ and $Y_e$ the partition function and prefactor vary weakly across neighbouring iron-group nuclei, so these two levers pick the winner.

Step 2. Evaluate the exponential. Binding energy per nucleon is maximal at the iron peak, so $e^{B_i/k_B T}$ drives the equilibrium there; among $Z=N$ nuclei, doubly magic $^{56}\mathrm{Ni}$ ( $Z=N=28$ ) carries the largest binding energy. Why: the exponential is the strongest term, making NSE a statement about which nucleus is most bound, not most abundant to start with.

Step 3. Impose the electron fraction. The weak interactions that change $Y_e$ are frozen, so the gas keeps its inherited $Y_e\approx0.5$ , i.e. $Y_p\approx Y_n$ ; then $Y_p^{Z}Y_n^{N}$ at fixed $A$ is largest when $Z=N$ . Why: equal free-nucleon abundances steer the equilibrium onto symmetric $Z=N$ nuclei rather than the neutron-rich $^{56}\mathrm{Fe}$ .

Step 4. Track the temperature down. Above $T_9\approx7$ photodisintegration dissolves even $^{56}\mathrm{Ni}$ into $\alpha$ particles and free nucleons; cooling through $T_9\sim5$ – $6$ re-forms $^{56}\mathrm{Ni}$ , which then freezes out. Why: the product is $^{56}\mathrm{Ni}$ , not $^{56}\mathrm{Fe}$ , because the slow weak decays that would make iron act only long after freeze-out.

Go deeper · Derivation B · 1 minThe alpha-rich freeze-out condition

Derivation BThe alpha-rich freeze-out condition

The goal is to derive why low-density, high-entropy complete silicon burning does not freeze out as pure iron-group material but strands a substantial mass fraction of free $\alpha$ particles — the $\alpha$ -rich freeze-out that leaves helium and $^{44}\mathrm{Ti}$ in the ejecta. The argument is a timescale comparison: the rate that would reassemble free $\alpha$ particles into carbon and beyond against the rate at which the gas expands and cools.

Step 1. Identify the bottleneck: reassembling free $\alpha$ particles must pass through triple-alpha, $3\,{}^{4}\mathrm{He}\to{}^{12}\mathrm{C}$ , since no stable nucleus at mass 5 or 8 lets it proceed two particles at a time. Why: this one three-body step throttles the whole reassembly funnel, so its rate decides whether the $\alpha$ particles are recaptured at all.

Step 2. Write its density dependence, $\tau_{3\alpha}^{-1}\propto\rho^{2}Y_\alpha^{3}\,\lambda_{3\alpha}(T)$ . Why: the $\rho^{2}$ signature of a three-body reaction is the crux — halving the density quarters the reassembly rate, while the expansion time lengthens only as $\rho^{-1/2}$ .

Step 3. Compare the clocks. Freeze-out sets in when $\tau_{3\alpha}\gtrsim\tau_{\mathrm{HD}}=446/\sqrt{\rho}\,$ s: at high density triple-alpha keeps pace and sweeps the $\alpha$ particles into iron-group nuclei (normal freeze-out); at low density it cannot. Why: the mismatch between $\rho^{2}$ and $\rho^{-1/2}$ makes dropping the density a one-way ticket to the $\alpha$ -rich regime.

Step 4. Read off the signature: enhanced $^{4}\mathrm{He}$ , extra $\alpha$ -chain nuclei, and altered iron-group ratios — most usefully enhanced $^{44}\mathrm{Ti}$ . Why: the leftover $\alpha$ particles reshape the yields in a way later gamma-ray telescopes can see.

Go deeper · Derivation C · 1 minPeak brightness set by the nickel-56 mass, Arnett's rule

Derivation CPeak brightness set by the nickel-56 mass, Arnett's rule

The goal is to connect the microscopic yield — the mass of $^{56}\mathrm{Ni}$ a shell synthesizes — to a supernova's single most observable property, its peak luminosity. The bridge is Arnett's rule: for a homologously expanding, radiation-dominated ejecta the luminosity at maximum equals the instantaneous radioactive heating rate, so the peak reads out the nickel mass directly.

Step 1. Write the heating rate from equation (7) in the trapped limit, $\dot E=\lambda_{\mathrm{Ni}}N_{\mathrm{Ni}}Q_{\mathrm{Ni}}+\lambda_{\mathrm{Co}}N_{\mathrm{Co}}Q_{\mathrm{Co}}$ , with $N_{\mathrm{Ni}}(0)=M_{\mathrm{Ni}}/(56\,m_u)$ . Why: every term is linear in the number of radioactive nuclei, so the whole heating rate scales with the nickel mass.

Step 2. Apply Arnett's rule. The ejecta expand homologously, so the diffusion time falls as the star grows and the light curve rises until leakage matches the shrinking internal store — at that instant $L_{\mathrm{peak}}\approx\dot E(t_{\mathrm{peak}})$ , and since $\dot E\propto M_{\mathrm{Ni}}$ , so is the peak. Why: the equality at maximum removes the unknown diffusion history, leaving luminosity proportional to nickel mass.

Step 3. Anchor with numbers. A normal Type Ia makes $M_{\mathrm{Ni}}\approx0.6\,M_\odot$ , giving initial deposition $\sim5\times10^{43}\,\mathrm{erg\,s^{-1}}$ and peak $L\sim10^{43}\,\mathrm{erg\,s^{-1}}$ near $t\approx18\,$ d; a core-collapse event makes only $M_{\mathrm{Ni}}\approx0.05$ – $0.1\,M_\odot$ and is fainter. Why: this quantitative link makes Type Ia supernovae calibratable standard candles.

Limiting cases

The regime is set by the peak shock temperature $T_9\equiv T/10^{9}\,\mathrm{K}$ at a given mass shell, which falls outward:

Case	Result	Meaning
Complete Si-burning / NSE	$T_9\gtrsim 5$	Full NSE; yield dominated by iron-group, principally $^{56}\mathrm{Ni}$ at $Y_e\approx0.5$ . At low density an $\alpha$ -rich freeze-out enhances $^{4}\mathrm{He}$ and $\alpha$ -chain nuclei.
Incomplete Si-burning	$4\lesssim T_9\lesssim 5$	Silicon only partly consumed; a quasi-equilibrium of Si- and Fe-groups. Products: $^{28}\mathrm{Si}$ , $^{32}\mathrm{S}$ , $^{40}\mathrm{Ca}$ with some $^{56}\mathrm{Ni}$ .
Explosive O- and Ne-burning	$3\lesssim T_9\lesssim 4$	Oxygen and neon burn to the $\alpha$ nuclei $^{28}\mathrm{Si}$ , $^{32}\mathrm{S}$ , $^{36}\mathrm{Ar}$ , $^{40}\mathrm{Ca}$ ; little iron-group.
Explosive C-burning / unprocessed	$T_9\lesssim 2$	Carbon burning yields $^{20}\mathrm{Ne}$ , $^{23}\mathrm{Na}$ , $^{24}\mathrm{Mg}$ ; below this the shock is too weak to alter the composition, which is simply ejected.

Common trapThe ejecta are not just the equilibrium NSE composition

It is tempting to compute the NSE abundances at the peak temperature and call that the yield. That is wrong except in the very hottest, densest zones. Expansion freezes the reactions out of equilibrium, and the freeze-out composition can differ sharply from the peak- $T$ NSE state. The clearest example is the $\alpha$ -rich freeze-out: the three-body triple-alpha reaction ( $\propto\rho^{2}$ ) cannot keep pace during rapid low-density cooling, so a surplus of free $\alpha$ particles pushes the yields toward $^{4}\mathrm{He}$ , $\alpha$ -chain nuclei and $^{44}\mathrm{Ti}$ rather than pure $^{56}\mathrm{Ni}$ . A second, independent trap is to think $T$ and $\rho$ alone fix the iron-group product — they do not; the electron fraction $Y_e$ does. Even a slight neutron excess ( $Y_e<0.5$ ) replaces $^{56}\mathrm{Ni}$ with neutron-richer $^{58}\mathrm{Ni}$ or $^{54}\mathrm{Fe}$ , which is why deep, electron-captured zones near a collapsing core give different iron-group ratios than the symmetric white-dwarf burning of a Type Ia.

Worked examples

1 · The two clocks in a silicon-burning shell (fully worked)

A mass shell is shock-heated to $T_9=5$ at a density $\rho=10^{7}\,\mathrm{g\,cm^{-3}}$ . Estimate how long the hot conditions last, and contrast it with a lower-density shell at $\rho=10^{5}\,\mathrm{g\,cm^{-3}}$ . What does the comparison say about complete versus $\alpha$ -rich freeze-out?

Step 1. Expansion clock at the dense shell from equation (1): $\tau_{\mathrm{HD}}=446/\sqrt{10^{7}}=\mathbf{0.141\ s}$ . The gas stays hot for roughly this long before cooling quenches it.

Step 2. At $T_9=5$ the silicon-burning reactions have $\tau_{\mathrm{nuc}}\ll10^{-3}\,$ s $\ll\tau_{\mathrm{HD}}$ , so the composition reaches NSE and, at $Y_e\approx0.5$ , relaxes onto $^{56}\mathrm{Ni}$ — iron-peak material.

Step 3. Drop the density to $\rho=10^{5}\,\mathrm{g\,cm^{-3}}$ : $\tau_{\mathrm{HD}}=446/\sqrt{10^{5}}=\mathbf{1.41\ s}$ , ten times longer. But triple-alpha reassembly scales as $\rho^{2}$ , so lowering $\rho$ by $10^{2}$ slows it by $10^{4}$ , overwhelming the tenfold gain in time.

Step 4. The reassembly loses the race, $\tau_{3\alpha}\gtrsim\tau_{\mathrm{HD}}$ , and free $\alpha$ particles are stranded. The dense shell (0.14 s) undergoes a normal freeze-out to $^{56}\mathrm{Ni}$ ; the rarefied shell (1.4 s) an $\alpha$ -rich freeze-out enriched in $^{4}\mathrm{He}$ and $^{44}\mathrm{Ti}$ — density, not temperature, is the switch.

2 · Peak luminosity from the nickel mass (fully worked)

A normal Type Ia supernova synthesizes $M_{\mathrm{Ni}}=0.6\,M_\odot$ of $^{56}\mathrm{Ni}$ . Using the $^{56}\mathrm{Ni}$ decay ( $t_{1/2}=6.08\,$ d, deposited energy $Q_{\mathrm{Ni}}\approx1.72\,$ MeV per decay), estimate the radioactive deposition rate just after the explosion and compare with the observed peak luminosity.

Step 1. Count the nuclei: $N_{\mathrm{Ni}}=M_{\mathrm{Ni}}/(56\,m_u)=(0.6\times1.989\times10^{33})/(56\times1.6605\times10^{-24})=\mathbf{1.28\times10^{55}}$ .

Step 2. Decay constant: $\lambda_{\mathrm{Ni}}=\ln2/t_{1/2}=0.693/(6.08\times86400\,\mathrm{s})=\mathbf{1.32\times10^{-6}\ s^{-1}}$ .

Step 3. Deposition at $t=0$ : $\dot E=\lambda_{\mathrm{Ni}}N_{\mathrm{Ni}}Q_{\mathrm{Ni}}=(1.32\times10^{-6})(1.28\times10^{55})(1.72\times1.602\times10^{-6})=\mathbf{4.7\times10^{43}\ erg\,s^{-1}}$ , a specific rate $\varepsilon_{\mathrm{Ni}}\approx3.9\times10^{10}\,\mathrm{erg\,s^{-1}\,g^{-1}}$ .

Step 4. By Arnett's rule the peak luminosity, near $t\approx18\,$ d after some $^{56}\mathrm{Ni}$ has decayed and photons diffuse out, is $L_{\mathrm{peak}}\sim10^{43}\,\mathrm{erg\,s^{-1}}$ — a few tenths of the $t=0$ deposition. A tenth of a solar mass more or less of $^{56}\mathrm{Ni}$ moves the peak brightness proportionally, making the light-curve peak a nickel-mass gauge. (For pure $^{56}\mathrm{Ni}$ , $Y_e=28/56=0.500$ ; a zone with $Y_e<0.5$ freezes out as $^{54}\mathrm{Fe}$ and $^{58}\mathrm{Ni}$ .)

Physics of the portrait

The idea that the elements are cooked inside stars and scattered by their explosions belongs, more than to anyone, to Fred Hoyle — who insisted the periodic table is a fossil record of stellar interiors and their catastrophes.

Origins · 1 minFred Hoyle and how this was discovered

Fred Hoyle

1915–2001 · stellar and explosive nucleosynthesis; architect of the B²FH synthesis

Fred Hoyle was the most audacious stellar theorist of the twentieth century. In 1946 he argued that the heavy elements could not be relics of the Big Bang but must be forged in stars, and in 1954 he predicted — from the mere fact that carbon is abundant — a specific excited state of $^{12}\mathrm{C}$ that lets the triple-alpha reaction proceed, a resonance experimenters then found where he said. The 1957 review by Margaret Burbidge, Geoffrey Burbidge, William Fowler and Hoyle (B²FH) laid out the full grammar of stellar nucleosynthesis, including the e-process of equilibrium iron-group synthesis and the explosive burning a shock ignites. With Fowler in 1964 he worked out the hydrodynamic timescale governing whether explosive reactions reach equilibrium or freeze — the $\tau_{\mathrm{HD}}\propto\rho^{-1/2}$ clock at the heart of this article. Hoyle was also famously wrong in public — a champion of steady-state cosmology who coined "Big Bang" as a term of derision — yet his conviction that nuclear physics written into starlight could be read back out defined the field. The Nobel that followed went to Fowler in 1983; Hoyle's omission remains one of the committee's most debated.

History

1946

— Hoyle argues the heavy elements are synthesized inside stars, not the early universe

1957

— Burbidge, Burbidge, Fowler and Hoyle (B²FH) map the stellar nucleosynthesis processes

1964

— Fowler and Hoyle set the hydrodynamic timescale governing explosive burning

1973

— Woosley, Arnett and Clayton compute explosive silicon, oxygen and carbon burning yields

1987

— SN 1987A shows the ⁵⁶Co gamma-ray lines, confirming explosive nickel production

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 15 and 18 · Clayton, Principles of Stellar Evolution and Nucleosynthesis (1983), Ch. 7 · Arnett, Supernovae and Nucleosynthesis (1996) · Woosley, Heger & Weaver, Rev. Mod. Phys. 74, 1015 (2002) · Thielemann, Nomoto & Hashimoto, ApJ 460, 408 (1996). The hydrodynamic-timescale coefficient is $\tau_{\mathrm{HD}}=(24\pi G\rho)^{-1/2}=445.8\,\rho^{-1/2}\,\mathrm{s}\approx446\,\rho^{-1/2}\,\mathrm{s}$ with $G=6.674\times10^{-8}\,\mathrm{cgs}$ ; note that the uniform-sphere free-fall time $\sqrt{3\pi/32G\rho}\approx2101\,\rho^{-1/2}\,\mathrm{s}$ is $4.71\times$ longer — a different definition, not to be conflated. Half-lives 6.08 d ( $^{56}\mathrm{Ni}$ ) and 77.2 d ( $^{56}\mathrm{Co}$ ) are standard; the $^{56}\mathrm{Co}$ tail slope is $2.5/(\ln 10)\times\ln2/77.2\,\mathrm{d}=0.0097\,\mathrm{mag\,day^{-1}}$ ; gamma lines 847/1238 keV are textbook. Burning-zone thresholds and $^{56}\mathrm{Ni}$ yields ( $\sim0.6\,M_\odot$ Type Ia; $0.05$ – $0.1\,M_\odot$ core-collapse), the NSE Saha law, $\alpha$ -rich freeze-out, and Arnett's rule are standard results in the literature.

Prerequisites. Advanced Burning Stages in Massive Stars · Helium Burning and the Triple-Alpha Process · Nuclear Reaction Rates and the Gamow Peak · White Dwarfs and Electron Degeneracy Pressure. Leads to. Asymptotic Giant Branch and Thermal Pulses.

Ask the tutor: Ask me to walk through the NSE Saha law step by step, or to show how the two clocks $\tau_{\mathrm{nuc}}$ and $\tau_{\mathrm{HD}}$ decide the fate of any given shell.

Why does the explosion make radioactive nickel instead of stable iron directly?

Because the weak interactions that would build stable iron are far too slow. On the seconds-long explosion timescale, electron captures and beta decays cannot change the electron fraction $Y_e$ from its inherited value near $0.5$ . With equal protons and neutrons, NSE settles on the most bound symmetric ( $Z=N$ ) nucleus, $^{56}\mathrm{Ni}$ — not $^{56}\mathrm{Fe}$ , which needs a neutron excess. Nature makes $^{56}\mathrm{Ni}$ promptly and lets it decay to iron over the following months, and that slow decay powers the visible light — the stable iron in your blood was, briefly, nickel.

What exactly is nuclear statistical equilibrium?

NSE is the nuclear analogue of chemical equilibrium. When the temperature is high enough (roughly $T_9\gtrsim5$ ), every strong and electromagnetic reaction runs so fast that each is balanced by its reverse. The abundance of any nucleus is then fixed not by reaction history but by statistical mechanics, through a Saha-like formula. The payoff is that the whole composition collapses onto a function of just three numbers — $T$ , $\rho$ and $Y_e$ — favouring the most tightly bound nuclei accessible at that $Y_e$ .

Why does density, not temperature, control the alpha-rich freeze-out?

Because the reassembly bottleneck is a three-body reaction. Turning free $\alpha$ particles back into heavier nuclei must pass through triple-alpha, $3\,{}^{4}\mathrm{He}\to{}^{12}\mathrm{C}$ , whose rate scales as $\rho^{2}$ — two independent encounters must overlap — while the expansion clock lengthens only as $\rho^{-1/2}$ . Lowering the density therefore collapses the reassembly rate far faster than it lengthens the available time, so the $\alpha$ particles are not swept up before the gas cools. Temperature sets whether everything dissolved into $\alpha$ particles; density sets whether they can be reassembled — low-density, high-entropy zones freeze out helium-rich, seeding $^{44}\mathrm{Ti}$ .

How does peak temperature stratify the ejecta into zones?

The shock deposits the most energy nearest the center, so peak temperature falls with radius. Each nucleus reacts only above a threshold temperature set by its Coulomb barrier, and because rates are exponentially temperature-sensitive, the "reacts"/"frozen" transition is sharp. Shells at $T_9\gtrsim5$ attain full NSE and make iron-group; $4\lesssim T_9\lesssim5$ gives incomplete silicon burning (Si, S, Ca plus some Ni); $3\lesssim T_9\lesssim4$ drives explosive O/Ne-burning (Si, S, Ar, Ca); below $T_9\sim2$ the composition is unchanged. One explosion writes concentric shells of distinct chemistry.

What is Arnett's rule and why is it so useful?

Arnett's rule says that at the peak of a supernova light curve the luminosity equals the instantaneous radioactive heating rate, because a homologously expanding, radiation-dominated ejecta brightens until the growing photon leakage balances the shrinking internal store. Its value is practical: since the decay heating is proportional to the $^{56}\mathrm{Ni}$ mass, the peak brightness measures that mass directly, independent of the diffusion history — letting astronomers read a nickel yield off a light curve and calibrate Type Ia supernovae as distance indicators.

Where do the intermediate-mass elements like silicon and calcium come from, if the hottest zones make iron?

From the cooler shells that never reach full equilibrium. In the incomplete-silicon-burning zone ( $4\lesssim T_9\lesssim5$ ) a quasi-equilibrium of a silicon group and an iron group leaves $^{28}\mathrm{Si}$ , $^{32}\mathrm{S}$ and $^{40}\mathrm{Ca}$ alongside some $^{56}\mathrm{Ni}$ . Cooler still, explosive O/Ne-burning ( $3\lesssim T_9\lesssim4$ ) makes the $\alpha$ -chain nuclei $^{28}\mathrm{Si}$ , $^{32}\mathrm{S}$ , $^{36}\mathrm{Ar}$ and $^{40}\mathrm{Ca}$ with almost no iron group — the $\alpha$ -elements whose ratio to iron tracks a galaxy's star-formation history.