Case Study 2: Onion Shells — The Nuclear Structure of a Dying Star

The Star: A Massive Object at the End of Its Life

Consider a star born with $25 \, M_\odot$ — roughly 25 times the mass of the Sun. It formed from a cloud of gas and dust about 7 million years ago, at the same time as dinosaurs still walked the Earth (in a human-timescale analogy). For 99.99% of its life, it burned hydrogen placidly on the main sequence. Now, in the last few hours before its death, its core is a layered structure of extraordinary complexity — a series of concentric shells, each undergoing a different nuclear burning stage, each at a different temperature and density, and each producing a different set of elements.

This is the onion-shell structure of a pre-supernova star, and it is the most complete physical realization of stellar nucleosynthesis theory.

The Layers in Detail

Modern stellar evolution codes (MESA, KEPLER, FRANEC, GENEC) can follow the evolution of a massive star from birth to core collapse, tracking the detailed composition of every mass element. The results for a $25 \, M_\odot$ star at the pre-supernova stage (Woosley, Heger, and Weaver, 2002; Sukhbold et al., 2016) give the following structure:

Layer 1: The Iron Core ($M_r < 1.5 \, M_\odot$)

• Temperature: $T \approx 5 \times 10^9$ K ($k_BT \approx 430$ keV)
• Density: $\rho \approx 3 \times 10^9$ g/cm$^3$
• Radius: $\sim 3{,}000$ km (smaller than the Earth)
• Composition: ${}^{56}$Fe ($\sim 70\%$ by mass), ${}^{54}$Fe ($\sim 10\%$), ${}^{58}$Ni ($\sim 5\%$), other iron-peak species, free protons and neutrons
• Electron fraction: $Y_e \approx 0.42$ (neutron-rich, due to electron captures during silicon burning)

The iron core is inert — no exothermic nuclear reactions are possible. It is supported against gravity by electron degeneracy pressure, growing in mass as the silicon-burning shell feeds it from above. When the iron core exceeds the Chandrasekhar mass ($\sim 1.4 \, M_\odot$ for $Y_e \approx 0.42$), collapse is inevitable.

Nuclear physics at work: The iron-peak composition is the product of nuclear statistical equilibrium (NSE). At the extreme temperatures of the core, all reactions and their reverses are in balance, and the composition is determined by the binding energies alone. The dominance of iron-peak nuclei reflects the maximum of $B/A$ near $A \sim 56$.

Layer 2: The Silicon-Burning Shell ($1.5 < M_r < 2.0 \, M_\odot$)

• Temperature: $T \approx 3.5 \times 10^9$ K
• Density: $\rho \approx 10^8$ g/cm$^3$
• Composition: ${}^{28}$Si ($\sim 35\%$), ${}^{32}$S ($\sim 20\%$), ${}^{36}$Ar, ${}^{40}$Ca, and increasing amounts of iron-peak species as the shell burns
• Duration of current burning: $\sim 1$ day

The silicon-burning shell is where photodisintegration rearrangement is actively converting silicon-group nuclei to iron-peak nuclei. The alpha particles, protons, and neutrons released by photodisintegration are captured faster than they can be produced, building up the nuclear mass number step by step. This shell is feeding iron-peak material onto the core, increasing its mass toward the Chandrasekhar limit.

Layer 3: The Oxygen-Burning Shell ($2.0 < M_r < 3.5 \, M_\odot$)

• Temperature: $T \approx 2 \times 10^9$ K
• Density: $\rho \approx 10^7$ g/cm$^3$
• Composition: ${}^{28}$Si ($\sim 40\%$), ${}^{32}$S ($\sim 25\%$), ${}^{16}$O (residual, $\sim 10\%$), ${}^{34}$S, ${}^{38}$Ar, ${}^{40}$Ca
• Duration: $\sim 6$ months

Active ${}^{16}$O $+ {}^{16}$O fusion reactions are producing silicon-group nuclei. The released protons, alpha particles, and neutrons drive secondary reactions that diversify the composition beyond the primary channels.

Layer 4: The Neon/Carbon-Burning Shell ($3.5 < M_r < 5.5 \, M_\odot$)

• Temperature: $T \approx (0.8\text{--}1.5) \times 10^9$ K
• Density: $\rho \approx 10^6$ g/cm$^3$
• Composition: ${}^{16}$O ($\sim 50\%$), ${}^{20}$Ne ($\sim 20\%$), ${}^{24}$Mg ($\sim 10\%$), ${}^{12}$C ($\sim 5\%$), ${}^{23}$Na
• Duration: C burning $\sim 600$ years; Ne burning $\sim 1$ year

This zone shows the products of carbon burning (${}^{20}$Ne, ${}^{23}$Na, ${}^{24}$Mg) mixed with residual oxygen from helium burning. The neon photodisintegration and alpha-capture cycle ($2 \, {}^{20}$Ne $\to {}^{16}$O $+ {}^{24}$Mg) operates in the hotter inner region of this shell.

Layer 5: The Helium-Burning Shell ($5.5 < M_r < 8 \, M_\odot$)

• Temperature: $T \approx 2 \times 10^8$ K
• Density: $\rho \approx 10^4$ g/cm$^3$
• Composition: ${}^{4}$He ($\sim 60\%$), ${}^{12}$C ($\sim 15\%$), ${}^{16}$O ($\sim 10\%$), ${}^{14}$N ($\sim 5\%$), ${}^{22}$Ne
• Duration: $\sim 5 \times 10^5$ years

The helium shell is actively producing ${}^{12}$C via the triple-alpha process and converting some to ${}^{16}$O via ${}^{12}$C$(\alpha,\gamma){}^{16}$O. The presence of ${}^{14}$N is a fossil of CNO-cycle hydrogen burning in the previous layer — ${}^{14}$N is slowly being converted to ${}^{22}$Ne via ${}^{14}$N$(\alpha,\gamma){}^{18}$F$(e^+\nu){}^{18}$O$(\alpha,\gamma){}^{22}$Ne, a reaction chain that is also the main source of free neutrons for the s-process (slow neutron capture, Chapter 23).

Layer 6: The Hydrogen-Burning Shell ($8 < M_r < 10 \, M_\odot$)

• Temperature: $T \approx 5 \times 10^7$ K
• Density: $\rho \approx 10$ g/cm$^3$
• Composition: ${}^{1}$H ($\sim 30\%$), ${}^{4}$He ($\sim 65\%$), ${}^{14}$N ($\sim 1\%$), CNO trace species
• Duration: Active throughout post-main-sequence evolution

The hydrogen-burning shell operates via the CNO cycle (the star is massive enough that CNO dominates over pp). The conversion of C and O to ${}^{14}$N is evident in the composition.

Layer 7: The Hydrogen Envelope ($10 < M_r < 15\text{--}20 \, M_\odot$)

• Temperature: $T \sim 10^4$ K (photosphere)
• Density: $\rho \sim 10^{-6}$–$10^{-9}$ g/cm$^3$
• Composition: ${}^{1}$H ($\sim 70\%$), ${}^{4}$He ($\sim 28\%$), metals ($\sim 2\%$) — essentially primordial composition
• Radius: $\sim 500$–$1000 \, R_\odot$ (the star is a red supergiant)

The envelope has never been processed by nuclear reactions. Its composition reflects the material from which the star formed — modified slightly by convective dredge-up of CNO-processed material from the hydrogen-burning shell.

Note: Some massive stars ($M > 25\text{--}30 \, M_\odot$) lose their hydrogen envelopes to stellar winds before core collapse, becoming Wolf-Rayet stars. In such cases, the outermost layer is the helium shell, and the resulting supernova is classified as a Type Ib or Ic (hydrogen-deficient).

The Physics Written in the Layers

Temperature and Density Gradients

The onion-shell structure encodes a remarkable set of thermodynamic conditions. From the iron core to the hydrogen envelope:

• Temperature spans five orders of magnitude: $5 \times 10^9$ K to $10^4$ K
• Density spans fifteen orders of magnitude: $3 \times 10^9$ g/cm$^3$ to $10^{-6}$ g/cm$^3$
• Pressure spans twenty orders of magnitude

Each burning shell maintains approximate hydrostatic equilibrium: the pressure gradient exactly balances gravity. But the stability is precarious — the entire structure is evolving on the nuclear timescale of the innermost burning shell, which at the silicon-burning stage is measured in hours.

The Mass Budget

A $25 \, M_\odot$ star yields approximately:

Of the total initial mass, about $5$–$10 \, M_\odot$ may be lost to stellar winds before collapse, and the supernova ejects most of the remainder. Only $\sim 1.4 \, M_\odot$ is retained as the neutron star remnant.

Observational Evidence

The onion-shell structure is not merely theoretical — it is confirmed by observations:

1. Supernova spectra. As the ejecta expand and thin out after the explosion, successively deeper layers become visible. Type II supernova spectra show hydrogen lines (from the envelope) fading over weeks, followed by helium, then oxygen and calcium, and eventually iron-peak lines — a spectral "peeling" of the onion.

2. Supernova remnants. Young supernova remnants like Cassiopeia A show distinct chemical layers in X-ray observations. Iron-rich knots (from the innermost layers) are spatially separated from silicon-rich and oxygen-rich regions, confirming the layered nucleosynthetic origin.

3. Nucleosynthesis yields. Integrated yields from stellar evolution models reproduce the observed cosmic abundance pattern when averaged over a realistic initial mass function (the distribution of stellar birth masses). The agreement is remarkable and constitutes strong evidence that the stellar nucleosynthesis framework is fundamentally correct.

4. SN 1987A. The detection of gamma-ray lines from ${}^{56}$Co decay in SN 1987A (by the Solar Maximum Mission and the Gamma-Ray Observatory) provided direct evidence that $\sim 0.07 \, M_\odot$ of ${}^{56}$Ni was synthesized in the explosion — confirming that the iron-peak products of silicon burning were present and were ejected.

Cassiopeia A: The Onion Shell in X-Rays

The most vivid observational confirmation of the onion-shell structure comes from the young supernova remnant Cassiopeia A (Cas A), located $\sim 3.4$ kpc away in the constellation Cassiopeia. The supernova that produced Cas A occurred around 1680 CE (it was not reliably observed, possibly due to interstellar dust absorption).

The Chandra X-ray Observatory has imaged Cas A with extraordinary spatial resolution, revealing:

• Iron-rich ejecta (from the iron core and silicon-burning shell) concentrated in the interior of the remnant, appearing as bright knots in the Fe K-alpha line at 6.6 keV.
• Silicon-rich ejecta (from oxygen and silicon burning) forming a thick shell surrounding the iron.
• Oxygen and neon-rich filaments (from carbon and neon burning) at larger radii.
• Hydrogen and helium largely absent, consistent with the progenitor having lost most of its hydrogen envelope before exploding (likely a Type IIb supernova from a $\sim 15$--$25 \, M_\odot$ progenitor).

Remarkably, some iron-rich knots are found outside the silicon-rich shell — they have been ejected at higher velocities than the overlying layers, a phenomenon called "jet-like" or "finger-like" overturn. This radial mixing, confirmed by three-dimensional simulations of the supernova explosion, demonstrates that while the onion-shell structure accurately describes the pre-supernova star, the explosion itself is highly asymmetric and partially inverts the layering. Understanding this mixing remains an active area of research.

Numerical Exercise

The Chandra image of Cas A shows iron-rich knots moving at $\sim 5{,}000$ km/s, while the outer blast wave expands at $\sim 6{,}000$ km/s. At the age of Cas A ($\sim 340$ years), the blast wave radius should be approximately:

$$R_{\text{blast}} = v_{\text{blast}} \times t \approx 6{,}000 \text{ km/s} \times 340 \text{ yr} \approx 6.4 \times 10^{13} \text{ km} \approx 2.1 \text{ pc}$$

The observed angular radius of $\sim 2.5'$ at a distance of 3.4 kpc gives $R \approx 2.5$ pc — in reasonable agreement, showing mild deceleration as the blast wave sweeps up interstellar material.

The Final Seconds

In the last day of the star's life, the silicon-burning shell is actively converting its fuel to iron-peak material, growing the iron core at a rate of $\sim 0.1 \, M_\odot$ per day. The neutrino luminosity from thermal processes ($\sim 10^{45}$ erg/s) far exceeds the photon luminosity ($\sim 10^{39}$ erg/s) — the star is "shining" primarily in neutrinos, invisible to optical telescopes.

When the iron core reaches $\sim 1.4$–$1.8 \, M_\odot$ (depending on $Y_e$ and temperature), the electron degeneracy pressure can no longer support it. Collapse begins, and within $\sim 0.3$ seconds, the core contracts from $\sim 3{,}000$ km to $\sim 15$ km — from Earth-sized to city-sized. The density rises to $\rho \sim 3 \times 10^{14}$ g/cm$^3$ (nuclear density), the temperature exceeds $10^{11}$ K, and a neutron star is born.

The story of what happens next — the bounce, the shock wave, the neutrino-driven explosion, and the explosive nucleosynthesis that occurs in the ejected material — is the subject of Chapter 23.

Discussion Questions

1. The iron core of the pre-supernova star contains $\sim 1.5 \, M_\odot$ compressed to a radius of $\sim 3{,}000$ km. Calculate the average density and compare it to (a) the density of nuclear matter ($\rho_0 \approx 2.3 \times 10^{14}$ g/cm$^3$), (b) the density of a white dwarf ($\sim 10^6$ g/cm$^3$), and (c) the density of water.

2. Why are the burning timescales so different from the Kelvin-Helmholtz timescale $\tau_{\text{KH}} = GM^2/(RL) \sim 3 \times 10^7$ years? Which timescale controls the star's evolution at each stage?

3. If you could somehow observe a massive star during its final day of silicon burning, what would it look like from the outside? Would any surface property reveal the drama occurring in the core? (Hint: the photon luminosity is essentially unchanged, but the neutrino luminosity is not.)

4. The onion-shell model assumes spherical symmetry. Modern three-dimensional simulations show that convective shells can develop large-scale asymmetries, including SASI (standing accretion shock instability) and convective overshooting. How might these asymmetries affect the nucleosynthetic yields?

Significance

The onion-shell structure of a pre-supernova star is perhaps the most beautiful illustration of nuclear physics in action. Each layer is a laboratory, each composition a record of specific nuclear reactions at specific temperatures. The entire structure — from the inert iron core to the pristine hydrogen envelope — can be understood quantitatively from the nuclear physics developed in this textbook: the binding energy curve (Chapter 1), the Gamow peak and reaction rate formalism (Chapter 21), the Breit-Wigner resonance formula (Chapter 18), and the nuclear reaction channels catalogued in this chapter.

The star is, in a very real sense, a nuclear physics experiment — conducted by gravity, over millions of years, at temperatures and densities we cannot achieve in the laboratory. And it writes its results in the periodic table of elements.

A $25 \, M_\odot$ star is born, burns hydrogen for seven million years, burns helium for half a million years, burns carbon for six centuries, burns oxygen for half a year, burns silicon for a day — and collapses in a third of a second. The acceleration is not poetic; it is quantitative, and every timescale can be calculated from nuclear cross sections and neutrino emission rates. The universe keeps good books.