The Semi-Empirical Mass Formula: The Five-Term Bethe-Weizsäcker Binding Energy Equation

What the formula is: the nucleus as a liquid drop

The semi-empirical mass formula estimates the total binding energy B(A, Z) — the energy released when Z protons and N = A − Z neutrons assemble into a nucleus — using an analogy George Gamow first suggested: the nucleus behaves like a tiny drop of charged, nearly incompressible nuclear fluid. Carl Friedrich von Weizsäcker cast this into a mass formula in 1935, and Hans Bethe and Robert Bacher systematized it in 1936, so it is usually called the Bethe-Weizsäcker formula.

The model is semi-empirical because its structure is derived from physics but its numbers are fitted to data. Two experimental facts anchor it. First, nuclear density is roughly constant (~0.16 nucleons/fm³), so nuclear radius scales as R ≈ r₀·A^(1/3) with r₀ ≈ 1.2 fm — the drop is incompressible. Second, binding energy per nucleon, B/A, is nearly flat at about 8 MeV across most of the chart, evidence that the strong force saturates: each nucleon interacts only with its nearest neighbors, not with all others.

The mechanism: building the five terms

Each term corrects the picture of an idealized drop:

• Volume, +aV·A: if every nucleon bound the same number of neighbors, binding would be strictly proportional to volume, hence to A. This is the dominant, attractive term.
• Surface, −aS·A^(2/3): nucleons on the drop's skin lack neighbors on the outside. Surface area scales as R² ∝ A^(2/3), so this term subtracts binding, most importantly for light nuclei with high surface-to-volume ratio.
• Coulomb, −aC·Z(Z−1)/A^(1/3): the Z protons repel one another. A uniformly charged sphere of charge Ze and radius R has electrostatic energy ∝ Z²/R ∝ Z²/A^(1/3); using Z(Z−1) counts distinct pairs. Here aC = (3/5)·e²/(4πε₀r₀) ≈ 0.71 MeV.
• Asymmetry, −aA·(A−2Z)²/A: the Pauli exclusion principle makes it costly to have very unequal neutron and proton numbers; the cost grows quadratically in the imbalance (N − Z) = (A − 2Z).
• Pairing, δ: a quantum correction. Like nucleons prefer to couple into spin-zero pairs, so even-even nuclei get extra binding (+), odd-odd nuclei less (−), with δ ≈ ±aP·A^(−1/2).

Key quantities and a worked example: iron-56

Take iron-56: A = 56, Z = 26, N = 30. It is even-even, so δ = +aP·A^(−1/2). Using aV = 15.8, aS = 18.3, aC = 0.71, aA = 23.2, aP = 12 MeV:

• Volume: 15.8 × 56 = +884.8 MeV
• Surface: −18.3 × 56^(2/3) = −18.3 × 14.64 = −267.9 MeV
• Coulomb: −0.71 × 26×25 / 56^(1/3) = −0.71 × 650 / 3.826 = −120.6 MeV
• Asymmetry: −23.2 × (56 − 52)² / 56 = −23.2 × 16/56 = −6.6 MeV
• Pairing: +12 / √56 = +1.6 MeV

Sum: B ≈ 491.3 MeV, giving B/A ≈ 8.77 MeV/nucleon. The measured value for ⁵⁶Fe is 492.3 MeV (8.79 MeV/nucleon) — the formula lands within 0.2 percent. That flat maximum of B/A near A ≈ 56–62 is why iron and nickel sit at the bottom of the energy valley: nothing binds more tightly per particle.

How it is measured and applied

The SEMF's coefficients are extracted by least-squares fitting to thousands of precisely measured masses, tabulated in the Atomic Mass Evaluation (e.g. AME2020). Masses come from Penning-trap mass spectrometry, which compares cyclotron frequencies of trapped ions to reach relative precision of 10⁻¹⁰–10⁻¹¹, and from reaction/decay Q-values. Because B(A, Z) fixes the nuclear mass, the formula predicts each nucleus's mass to a few hundred keV.

Its applications are broad:

• Valley of stability: minimizing M(A, Z) over Z at fixed A gives Z_optimal ≈ A / (1.98 + 0.015·A^(2/3)), predicting the most stable isobar and which nuclei β-decay toward it.
• Fission and fusion energetics: the balance of surface versus Coulomb terms determines whether a nucleus gains energy by splitting. The formula yields the fissility parameter x = (Z²/A)/(2aS/aC) ≈ (Z²/A)/50; nuclei above Z²/A ≈ 47–49 face vanishing fission barriers.
• Astrophysics and reactor physics use it to estimate Q-values, drip lines, and energy release.

Comparison to related models and regimes

The SEMF is a macroscopic model and should be contrasted with its cousins:

• Pure liquid-drop model: the SEMF is the liquid-drop model plus the asymmetry and pairing terms, which are quantum (Pauli/shell) in origin. Drop the last two and you cannot even explain why nuclei favor N ≈ Z.
• Nuclear shell model: a microscopic theory built on individual nucleon orbitals and magic numbers (2, 8, 20, 28, 50, 82, 126). The SEMF is smooth and misses the sharp extra binding at magic numbers — e.g. it under-binds doubly-magic ⁴⁰Ca, ²⁰⁸Pb by several MeV. Modern mass models (Weizsäcker-Skyrme, finite-range droplet model, HFB) add shell and deformation corrections on top of the drop.
• Regime of validity: the SEMF works best for medium-to-heavy nuclei (A ≳ 20–30). For very light nuclei (deuteron, helium-4) the surface term dominates and the smooth picture breaks down; ⁴He is far more bound than the formula predicts.

Significance, famous cases, and open questions

In five fitted numbers the Bethe-Weizsäcker formula captures the gross energetics of nearly all nuclei — a rare feat of physical compression. Its landmark consequences are the binding-energy curve peaking near iron, the explanation that fusing light nuclei or fissioning heavy ones both climb toward that peak and release energy, and the quantitative valley of stability. Bohr and Wheeler used exactly this surface-versus-Coulomb competition in their 1939 liquid-drop theory of fission, published just months after fission's discovery.

Famous limitations remain instructive. The formula's smooth ⁵⁶Fe-versus-⁶²Ni comparison actually gets subtle: ⁶²Ni has slightly higher B/A (8.795 MeV) than ⁵⁶Fe (8.790 MeV), a shell effect the SEMF cannot resolve. Open frontiers include the neutron drip line and superheavy elements, where extrapolating the asymmetry and Coulomb terms is uncertain and shell stabilization (the predicted "island of stability" near Z ≈ 114–120) dominates — regimes the plain SEMF cannot reach without microscopic corrections.

The five terms of the semi-empirical mass formula: physical origin, A/Z dependence, sign, and approximate coefficient.

Term	Physical origin	Scaling	Sign	Coefficient (MeV)
Volume	Short-range saturating strong force; each nucleon binds its neighbors	∝ A	+ (binds)	aV ≈ 15.8
Surface	Nucleons at the drop's surface have fewer neighbors	∝ A^(2/3)	− (reduces)	aS ≈ 18.3
Coulomb	Mutual electrostatic repulsion of Z protons	∝ Z(Z−1)/A^(1/3)	− (reduces)	aC ≈ 0.71
Asymmetry	Pauli/isospin cost of unequal N and Z	∝ (A−2Z)²/A	− (reduces)	aA ≈ 23.2
Pairing	Extra binding when protons/neutrons pair up	∝ ±A^(−1/2)	± (varies)	aP ≈ 12