Nuclear Binding Energy Calculator — Mass Defect, the Iron Peak and Reaction Q-Values
A nucleus weighs less than the parts it is made of. That is not a measurement error; it is the whole of nuclear energy in one sentence. This nuclear binding energy calculator takes the difference — the mass defect — and turns it into energy through E = Δm·c², reporting the total, the per-nucleon figure that decides whether a nuclide would rather fuse or split, and the mass excess and packing fraction that nuclear data tables quote alongside them.
How the mass defect becomes energy
Work with atomic masses, as every published table does, and the constituents are Z hydrogen atoms plus N neutrons, so the Z electrons cancel between the two sides:
Δm = Z·m_H + N·m_n − M_atom, then E_B = Δm × 931.49410242 MeV/u
For iron-56, 26 hydrogen atoms and 30 neutrons come to 56.463398 u, against a measured 55.934936 u. The missing 0.528462 u is worth 492.259 MeV, which across 56 nucleons is 8.7903 MeV per nucleon. Less than one per cent of the mass has gone, and it is worth roughly fifty million times the energy of a chemical bond.
Why the curve peaks near iron
Plot binding energy per nucleon against mass number and you get the most important graph in nuclear physics. It climbs steeply through the light elements, flattens into a broad plateau around A = 50–68, and declines slowly through the heavy ones. Two effects are competing: the strong force is short-ranged and saturates, so attraction grows roughly with A, while electrostatic repulsion between protons is long-ranged and grows with Z². Nickel-62 at 8.7945 MeV is the champion, with iron-56 a hair behind — iron wins the astrophysical argument only because supernova nucleosynthesis produces more of it.
Everything else can move downhill toward that plateau. Light nuclei get there by fusion; heavy nuclei get there by fission. Uranium-235 sits at 7.591 MeV per nucleon, about 1.2 MeV below the peak, and multiplying that shortfall across 235 nucleons is where the roughly 173 MeV of a fission event comes from.
The semi-empirical mass formula
When no measured mass exists, the Weizsäcker formula models the nucleus as a charged liquid drop and estimates the binding energy from Z and A alone:
E_B = a_V·A − a_S·A^(2/3) − a_C·Z(Z−1)/A^(1/3) − a_A·(A−2Z)²/A + δ
With the standard coefficients, iron-56 breaks down as +882.000 volume, −260.543 surface, −120.796 Coulomb, −6.771 asymmetry and +1.494 pairing, totalling 495.384 MeV. That is +3.124 MeV above the measured value, a deviation of +0.635 % — the model slightly over-binds nuclei near the peak.
Q-values, and the units that trip people up
A reaction's Q-value is the mass that disappears across it: Q = [Σm(reactants) − Σm(products)]·c². Positive means exothermic. For U-235 + n → Ba-141 + Kr-92 + 3n the answer is +173.28 MeV, which over 235 g/mol works out at 7.11 × 10¹³ J/kg, or about 17.0 kilotons of TNT per kilogram of uranium. Deuterium–tritium fusion yields only 17.59 MeV per event, but spread over five nucleons instead of 236 it is several times more energy per kilogram of fuel.
9.648533 × 10⁷ kJ/mol. Iron-56's 492.26 MeV is therefore 4.7496 × 10¹⁰ kJ/mol — equivalently 4.7496 × 10¹³ J/mol. The same digits appear in print under both labels, and only one of them is kilojoules.Atomic mass or nuclear mass?
Almost always atomic. Tables list atomic masses, and the atomic convention cancels the electrons to within a few hundred eV. If you feed a genuine bare nuclear mass in while the atomic basis is selected, the missing electrons inflate the binding energy by about 0.511 MeV per proton — 13 MeV for iron. Because nothing in nature exceeds 8.7945 MeV per nucleon, the calculator flags any result above 8.9 as an input mistake rather than a discovery.
Mass excess and packing fraction
Two older quantities describe the same data. The mass excess Δ = (M − A·u)·c² measures how far a mass sits from a whole number of mass units; for iron-56 it is about −60.61 MeV. Aston's packing fraction (M − A)/A is the same idea normalised per nucleon, −1.162 × 10⁻³ for iron, and plotting it was how the shape of the binding curve was first discovered in the 1920s — before anyone knew what held a nucleus together at all.