Case Study 4.1 — Why Iron Is the End of the Line: The Peak of the Binding Energy Curve

"The iron peak is the graveyard of nuclear burning in stars. Everything lighter fuses toward iron; everything heavier was built by neutron capture."

The Question

Every introductory astronomy textbook states that "fusion stops at iron." The binding energy per nucleon curve peaks near $A \approx 56$, and therefore nuclei lighter than iron release energy by fusing together, while nuclei heavier than iron release energy by splitting apart (fission). The peak of the $B/A$ curve is the reason iron and nickel are the most abundant heavy elements in the universe. But the story is more nuanced than it first appears. Let us examine it carefully.

The Peak Region in Detail

The experimental binding energies per nucleon for nuclei near the peak (AME2020 data):

Notice how tightly clustered these values are: the entire range from $^{56}$Ni to $^{62}$Ni spans only 0.15 MeV/nucleon. The top of the $B/A$ curve is remarkably flat — a "plateau" rather than a sharp peak. This flatness is itself significant: it means that nuclei in the range $50 \lesssim A \lesssim 65$ are all nearly equally stable, and the detailed ordering depends on subtle nuclear structure effects far below the SEMF's resolution.

The first surprise: $^{62}$Ni, not $^{56}$Fe, has the highest binding energy per nucleon. The value $B/A = 8.7945$ MeV for $^{62}$Ni exceeds $B/A = 8.7903$ MeV for $^{56}$Fe by about 4 keV/nucleon. This is a small but real difference, well within the precision of modern mass measurements (Penning trap measurements achieve uncertainties of eV or better).

The second surprise: $^{58}$Fe also has a slightly higher $B/A$ than $^{56}$Fe.

So why does every textbook say "iron"?

There is also a third candidate that is sometimes mentioned: $^{58}$Fe ($B/A = 8.7922$ MeV), which exceeds $^{56}$Fe but falls below $^{62}$Ni. And $^{60}$Ni ($B/A = 8.7807$ MeV) is also in the running. The entire "peak region" spans a range of only about 0.015 MeV/nucleon — the flattest part of the binding energy curve, where the macroscopic SEMF prediction and the microscopic shell corrections are competing at the limit of our ability to distinguish them.

The resolution of the "$^{56}$Fe vs. $^{62}$Ni" question turns out to involve stellar physics, not just nuclear physics. The answer to "which nucleus is most tightly bound?" is $^{62}$Ni. The answer to "what do stars produce at the end of their lives?" is $^{56}$Fe, but for reasons that have nothing to do with which nucleus sits at the absolute peak of the $B/A$ curve.

What the SEMF Tells Us

Using the SEMF with our standard parameters ($a_V = 15.75$, $a_S = 17.80$, $a_C = 0.711$, $a_{\text{sym}} = 23.7$, $a_P = 11.2$ MeV), we can find the mass number $A_{\text{peak}}$ that maximizes $B/A$ along the valley of stability. Taking the derivative $d(B/A)/dA = 0$ along the stability line $Z = Z_{\text{stable}}(A)$ and solving (numerically or approximately) gives $A_{\text{peak}} \approx 58$–$63$, depending on the exact fitting parameters.

The SEMF correctly predicts the existence of the peak and its approximate location, but it cannot distinguish between $^{56}$Fe, $^{58}$Fe, $^{62}$Ni, and their neighbors — the differences in $B/A$ are 4–5 keV/nucleon, while the SEMF accuracy is $\sim 30$ keV/nucleon. The detailed structure of the peak requires shell model calculations.

The SEMF also reveals why the peak exists:

• For $A \lesssim 60$: The surface term dominates the Coulomb term. $B/A$ increases with $A$ because the surface-to-volume ratio decreases.
• For $A \gtrsim 60$: The Coulomb term dominates the surface term. $B/A$ decreases with $A$ because the growing proton number makes the Coulomb penalty per nucleon ever larger.

The peak occurs where $d(B/A)/dA = 0$: the point at which the marginal decrease in surface energy per nucleon is exactly balanced by the marginal increase in Coulomb energy per nucleon.

Why Stars Care About Iron, Not Nickel

If $^{62}$Ni has the highest $B/A$, why does stellar nucleosynthesis produce iron, not nickel? The answer involves several layers of physics:

1. Silicon burning produces $^{56}$Ni, not $^{56}$Fe. In the final stage of nuclear burning in massive stars (Chapter 22), silicon undergoes quasi-equilibrium nuclear statistical processes at temperatures $T \sim 3$–$5 \times 10^9$ K. The equilibrium strongly favors nuclei with $N = Z$ (because the timescale is too short for weak interactions to establish the $N > Z$ preference of the valley of stability). The most bound $N = Z$ nucleus in this mass range is $^{56}$Ni ($Z = N = 28$, doubly magic). So silicon burning produces $^{56}$Ni.

2. $^{56}$Ni decays to $^{56}$Fe after the explosion. $^{56}$Ni is unstable — it undergoes electron capture to $^{56}$Co ($t_{1/2} = 6.075$ days), which in turn decays to $^{56}$Fe ($t_{1/2} = 77.236$ days). These decays occur after the core-collapse supernova explosion has already ejected the material. The gamma rays from $^{56}$Co decay are, in fact, what powers the light curve of Type Ia and core-collapse supernovae during the first few months.

3. The relevant thermodynamic quantity is not $B/A$ but rather the Helmholtz free energy at the relevant temperature and density. At the extreme conditions of stellar cores ($T > 10^9$ K, $\rho > 10^7$ g/cm$^3$), the equilibrium composition depends on temperature, density, and the electron fraction $Y_e = Z/A$. The equilibrium favors $^{56}$Ni (and nearby nuclei) at the conditions where silicon burning occurs, even though $^{62}$Ni has a higher $B/A$ at zero temperature. The electron fraction $Y_e$ is crucial: in the core of a massive star, electron capture on protons ($p + e^- \to n + \nu_e$) reduces $Y_e$ below 0.5, favoring neutron-rich nuclei. The freeze-out value of $Y_e$ determines the precise mix of iron-peak isotopes ejected in a supernova.

4. Fusion beyond $^{56}$Ni requires photodisintegration. To build nuclei heavier than $^{56}$Ni through charged-particle fusion, you need temperatures so high that photodisintegration (the reverse process, where a gamma ray breaks apart a nucleus) becomes faster than fusion. This creates a bottleneck: the core is in nuclear statistical equilibrium (NSE), and the equilibrium composition at the relevant temperatures clusters around the iron peak.

The Iron Peak in the Solar Abundance Pattern

The solar system abundances show a prominent peak at $A = 56$ (iron), with significant abundances of $^{54}$Fe, $^{56}$Fe, $^{57}$Fe, $^{58}$Ni, and $^{60}$Ni. This "iron peak" is the direct signature of nuclear statistical equilibrium in stellar cores. The detailed abundances within the peak are set by the freeze-out conditions (temperature, density, $Y_e$) when the explosive material cools and NSE breaks down.

The iron peak is flanked by: - A deep trough at $A \approx 45$ (the "titanium-vanadium valley"), where abundances are low because these nuclei are not produced efficiently by any burning stage. - A slow decline from $A \approx 62$ to $A \approx 90$, populated primarily by the weak $s$-process in massive stars and the $\alpha$-rich freeze-out.

The detailed shape of the iron peak abundance pattern encodes information about the conditions in the supernova core at the moment of explosion: the temperature, the density, the electron fraction $Y_e$, and the entropy. Nuclear astrophysicists use the observed iron-peak abundances as constraints on supernova models — the nuclear physics (binding energies, decay rates, electron capture rates) and the astrophysics (explosion mechanism, neutrino transport) must work together to reproduce the observed pattern.

Energy Release in Fusion and Fission

The $B/A$ curve has direct, quantitative implications for energy release in nuclear reactions. The energy released per nucleon when fusing two nuclei to form a product near the iron peak, or when fissioning a heavy nucleus into fragments near the iron peak, is simply the difference in $B/A$ between the reactants and products.

Fusion example: Consider the fusion of four protons into $^4$He (the net effect of the $pp$-chain in the Sun). The binding energy per nucleon of $^4$He is $B/A = 7.074$ MeV, while free protons have $B/A = 0$. The energy release per nucleon is therefore about 7 MeV — or about 26.7 MeV total for each helium nucleus formed. This is approximately 0.7% of the rest-mass energy of the four protons ($4 \times 938.3$ MeV), converted to radiation and neutrinos. This 0.7% mass-to-energy conversion efficiency is what powers the Sun.

Fission example: In the fission of $^{235}$U into two fragments near $A \approx 118$ (a symmetric split; in practice fission is asymmetric), $B/A$ increases from about 7.59 MeV to about 8.50 MeV. The energy release per nucleon is $\sim 0.9$ MeV, or about $0.9 \times 235 \approx 210$ MeV total per fission. This energy — released as kinetic energy of the fission fragments, prompt neutrons, and gamma rays — is what powers nuclear reactors and nuclear weapons.

The fact that both fusion and fission release energy, and that the energy release in both cases can be read directly from the $B/A$ curve, is one of the most consequential predictions of the SEMF.

Quantitative Analysis

Let us verify the SEMF's prediction for the peak location analytically. Along the stability line, with $Z_{\text{stable}} \approx A / (2 + \eta A^{2/3})$ where $\eta = a_C / (4 a_{\text{sym}})$, the binding energy per nucleon is approximately:

$$\frac{B}{A} \approx a_V - a_S A^{-1/3} - \frac{a_C}{4} \frac{A^{2/3}}{(1 + \eta A^{2/3}/2)^2}$$

The first two terms (volume minus surface) give a monotonically increasing function. The third term (Coulomb, evaluated along the stability line) gives a contribution that increases with $A$. The peak occurs where:

$$\frac{d(B/A)}{dA} = \frac{a_S}{3} A^{-4/3} - \frac{a_C}{6} \frac{A^{-1/3}}{(1 + \eta A^{2/3}/2)^2} + \ldots = 0$$

Solving numerically with our standard parameters gives $A_{\text{peak}} \approx 58$, consistent with the experimental observation.

Lessons

1. The SEMF correctly predicts the existence and approximate location of the $B/A$ peak. This is a macroscopic result — it follows from the competition between surface and Coulomb energies and does not require shell structure.

2. The detailed structure of the peak ($^{62}$Ni vs. $^{56}$Fe) is a shell effect. The SEMF cannot resolve differences of $\sim 5$ keV/nucleon. $^{62}$Ni benefits from the $Z = 28$ proton shell closure.

3. Stellar nucleosynthesis produces $^{56}$Ni (doubly magic, $N = Z$), which decays to $^{56}$Fe. The dominance of iron is a nuclear structure effect (magic numbers) combined with weak interaction physics (electron capture decay chain).

4. The "iron is the end of the line" statement is correct in the following precise sense: Exothermic charged-particle nuclear reactions cannot proceed beyond the iron peak under conditions of nuclear statistical equilibrium. Nuclei heavier than the iron peak are built by neutron capture ($s$-process and $r$-process), which circumvents the Coulomb barrier entirely. This is the subject of Chapter 22 (nucleosynthesis) and Chapter 35 (nuclear astrophysics).

5. The $B/A$ curve is the most important single plot in nuclear physics. It explains why fusion powers stars (light nuclei fuse to release energy), why fission powers reactors (heavy nuclei split to release energy), and why the universe is made mostly of hydrogen and helium (building heavier nuclei from lighter ones releases energy only up to the iron peak, and the universe has not had enough time or the right conditions to convert all light elements to iron).

The Iron Core and Stellar Death

When a massive star ($M \gtrsim 8 M_\odot$) exhausts silicon burning in its core and builds up an iron-nickel core, it has reached the end of exothermic nuclear burning. The iron core cannot release energy by either fusion or fission (both would require energy input). As the core grows beyond the Chandrasekhar mass ($\sim 1.4 M_\odot$), electron degeneracy pressure can no longer support it against gravity. The core collapses in about 0.1 seconds to nuclear density, rebounds, and launches a shock wave that (in most cases) tears the star apart in a core-collapse supernova.

The iron peak in the $B/A$ curve is therefore not merely an academic curiosity — it is the reason massive stars die. Without the maximum in $B/A$, stellar evolution would proceed differently: either all matter would fuse to the heaviest possible nucleus (if $B/A$ increased without bound) or fusion would stop at some lighter element (if the peak occurred elsewhere). The specific location of the peak near $A = 56$ determines the mass of the iron core at the moment of collapse, which in turn determines the energetics of the supernova explosion, the mass of the neutron star remnant, and the nucleosynthetic yields of elements heavier than iron.

The SEMF, through its prediction of the $B/A$ peak, provides the nuclear physics foundation for understanding stellar death — one of the most dramatic events in the universe.

The iron peak is where nuclear physics meets astrophysics. The shape of the binding energy curve — predicted by five numbers in the SEMF — determines the ultimate fate of matter in the universe.