Liquid Drop Model

Definition

The liquid drop model treats the atomic nucleus as a charged drop of incompressible nuclear fluid. The protons and neutrons inside it behave like molecules in a tiny droplet of water: each one feels only its immediate neighbors through the short-range strong force, the surface costs energy the way surface tension does, and the positive charge of the protons tries to blow the drop apart.

From those few analogies you can write down a formula for how tightly any nucleus is bound. That formula — the semi-empirical mass formula, also called the Bethe–Weizsäcker formula — is the centerpiece of the model:

B = a_V·A − a_S·A^(2/3) − a_C·Z²/A^(1/3) − a_A·(A−2Z)²/A ± δ

Here A is the mass number (total nucleons), Z is the proton number, and B is the total binding energy — the energy you would have to supply to pull the nucleus completely apart into free nucleons. It is "semi-empirical" because the shape of each term comes from physics, but the five coefficients are fitted to measured masses.

The five energy terms

Each term encodes a different piece of physics. Read them as a tug-of-war: the first term wants to bind the drop together, the rest chip away at that binding.

Term	Expression	Coefficient	Physical meaning
Volume	+ a_V·A	≈ 15.8 MeV	Each nucleon binds to its neighbors; total binding scales with the number of nucleons (the bulk).
Surface	− a_S·A^(2/3)	≈ 18.3 MeV	Surface nucleons have fewer neighbors — like surface tension, this costs energy. Scales with surface area ∝ A^(2/3).
Coulomb	− a_C·Z²/A^(1/3)	≈ 0.71 MeV	Protons mutually repel. The electrostatic energy of a uniformly charged sphere of radius ∝ A^(1/3).
Asymmetry	− a_A·(A−2Z)²/A	≈ 23.2 MeV	Nuclei prefer N ≈ Z. A Pauli-exclusion penalty for filling proton and neutron levels unevenly.
Pairing	± δ	≈ 12/√A MeV	Like-nucleon pairs gain binding. +δ for even-even, 0 for odd-A, −δ for odd-odd nuclei.

The volume term is the only one with a plus sign — it is the engine of binding. Everything else is a correction that explains why the curve has the shape it does.

How it works: the binding-energy curve

Divide the binding energy by the number of nucleons and you get B/A, the binding energy per nucleon. This is the single most important plot in nuclear physics, because it tells you which way energy flows.

For light nuclei, B/A climbs steeply. Adding a nucleon adds full volume binding (a_V per nucleon) while the surface penalty, which scales only as A^(2/3), gets diluted as the drop grows. So fusing light nuclei releases energy.

For heavy nuclei, the Coulomb term takes over. Because it grows as Z²/A^(1/3) — faster than the linear volume term — every extra proton pays a steeper and steeper electrostatic price. So B/A bends back down. Splitting heavy nuclei releases energy.

The two trends cross in a broad maximum near A ≈ 56–62. The peak of the curve sits at about 8.8 MeV per nucleon around iron-56 (nickel-62 is fractionally higher by total binding, but Fe-56 is the famous landmark). Iron-56 is the bottom of the energy valley: you cannot release energy by fusing it or by splitting it. That is precisely why the cores of massive stars stall at iron and then collapse.

Worked example: binding energy of iron-56

Let's plug iron-56 into the formula. For Fe-56, A = 56 and Z = 26 (so N = 30). Using the standard coefficients a_V = 15.8, a_S = 18.3, a_C = 0.71, a_A = 23.2 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²/56^(1/3) = −0.71×676/3.83 = −125.3 MeV
Asymmetry:  −23.2 × (56−52)²/56 = −23.2×16/56   =  −6.6 MeV
Pairing:    even-even → +12/√56 = +12/7.48      =  +1.6 MeV
                                                 ───────────
Total binding B ≈ 884.8 − 267.9 − 125.3 − 6.6 + 1.6 ≈ 486.6 MeV
Per nucleon:  B/A = 486.6 / 56 ≈ 8.69 MeV/nucleon

The measured value is 492.3 MeV total, 8.79 MeV/nucleon. Our five-term estimate lands within about 1.2 percent of reality — from a formula you can evaluate on the back of an envelope. That is the headline achievement of the liquid drop model: thousands of nuclei, five constants, ~1 percent error.

Predicting fission: the deformation contest

The model's most dramatic prediction is nuclear fission. Imagine slowly squashing a spherical nucleus into an ellipsoid. Two terms respond:

• Surface energy rises. Deforming a sphere increases its surface area, so the surface term (which penalizes area) gets larger. This resists deformation.
• Coulomb energy falls. Stretching the drop spreads the protons farther apart on average, lowering electrostatic energy. This encourages deformation.

Expand both to second order in the deformation and you find a fission barrier. Whether the barrier exists at all depends on the ratio of Coulomb to surface energy — captured by the fissility parameter Z²/A. Bohr and Wheeler (1939) showed that when

Z²/A  ≳  2·a_S / a_C  ≈  47   (critical, spontaneous limit)

the barrier vanishes entirely and the nucleus is unstable to instantaneous fission. Real heavy nuclei have lower values (uranium-235 has Z²/A ≈ 36) so they need a small kick — a neutron — to climb over the barrier. As a practical rule of thumb, fission becomes energetically favorable (positive Q) once Z²/A exceeds roughly 17, which is why everything heavier than about niobium is, in principle, unstable to splitting given enough time.

Variants and refinements

• Bethe–Weizsäcker (1935). The classic five-term form above. Still the version taught and used for quick estimates.
• Bohr–Wheeler fission theory (1939). Adds the deformation-dependent surface and Coulomb energies, yielding the fission barrier and the fissility parameter.
• Liquid-drop-plus-shell (Strutinsky, 1967). The macroscopic-microscopic method: take the smooth liquid-drop energy and add a "shell correction" so the magic numbers reappear. This is the workhorse for predicting superheavy elements and fission barriers today.
• Droplet and finite-range models (Myers–Swiatecki, Möller–Nix). Refine the surface (a diffuse skin, not a sharp edge), add a curvature term, and treat the neutron skin. These push global mass-fit errors below ~0.6 MeV.
• Collective vibrations and rotations. Treating the drop's surface as oscillating gives the collective model of Aage Bohr and Ben Mottelson — quantized surface "phonons" and giant dipole resonances.

Liquid drop model vs the shell model

The liquid drop model and the shell model are complementary. One captures the smooth bulk trend; the other captures the quantum bumps.

Aspect	Liquid drop model	Shell model
Picture of the nucleus	Structureless charged fluid	Nucleons in quantized orbitals
Best at	Smooth bulk binding, fission, masses away from closed shells	Magic numbers, spins, parities, excited states
Magic numbers (2, 8, 20, 28, 50, 82, 126)	Misses them entirely	Predicts them from shell closures
Number of free parameters	~5 global constants	Many (interaction matrix elements)
Fission and deformation	Natural — its signature success	Hard without large-scale computation
Year / pioneers	1935 · Gamow, Weizsäcker, Bethe	1949 · Goeppert Mayer, Jensen

Common pitfalls and misconceptions

• "It's literally a drop of liquid." No — it's an analogy. Nucleons are quantum particles, not classical molecules. The model works because the strong force is short-ranged and saturating, which makes nuclear matter behave like an incompressible fluid in bulk.
• Forgetting the minus signs. Four of the five terms reduce binding. Students who add everything get nonsense. Only the volume term binds; the rest are penalties.
• Expecting magic numbers. The liquid drop curve is smooth. If you need the extra stability at Z or N = 50, 82, 126, you must add shell corrections — the drop alone is blind to them.
• Confusing total B with B/A. Total binding energy keeps rising with A (more nucleons, more binding); it is binding per nucleon that peaks at iron-56 and then declines.
• Misreading the Coulomb scaling. It is Z²/A^(1/3), not Z²/A. The A^(1/3) is the nuclear radius (R ∝ A^(1/3)); the Z² is every proton pair repelling every other.
• Thinking fission needs Z²/A above 47. That's the spontaneous, barrier-free limit. Induced fission (uranium plus a neutron) happens far below it, around Z²/A ≈ 36, because the neutron supplies the energy to surmount the barrier.

Where the model shows up

• Nuclear energy. The fission Q-value and barrier estimates that underpin reactor and weapon physics come straight from the surface–Coulomb competition.
• Stellar nucleosynthesis. The iron-56 peak explains why fusion in stars releases energy up to iron and no further, setting the stage for core-collapse supernovae.
• The valley of stability. Minimizing the mass formula over Z at fixed A predicts the most stable isotope of each element — the line of beta stability.
• Superheavy element hunting. Liquid-drop-plus-shell models predict the fission barriers and the rumored "island of stability" near Z ≈ 114–126.
• Beta-decay energetics. The parabola of mass vs Z at fixed A (from the asymmetry and Coulomb terms) tells you which way a nucleus will beta-decay and how fast.
• Mass tables and astrophysical r-process. Global mass formulas seed the network calculations that model how heavy elements form in neutron-star mergers.

Performance and derivation analysis

Why does such a crude picture work so well? The deep reason is the saturation of nuclear forces. The strong force is short-ranged, so a nucleon deep inside the nucleus interacts only with a fixed number of neighbors regardless of how big the nucleus is. That fixed binding-per-nucleon is exactly the volume term. The constant nuclear density (about 0.16 nucleons/fm³ everywhere) is the same statement: nuclear matter is nearly incompressible, just like a liquid.

Each correction then has a clean geometric origin. The surface term subtracts the binding the surface nucleons are missing, and surface area scales as R² ∝ A^(2/3). The Coulomb term is the textbook electrostatic self-energy of a uniformly charged sphere, (3/5)·(Z²e²)/(4πε₀R), with R ∝ A^(1/3) — that is where Z²/A^(1/3) comes from, and the 3/5 sphere factor is baked into a_C ≈ 0.71 MeV. The asymmetry term is a Fermi-gas result: filling proton and neutron levels unequally forces nucleons into higher kinetic-energy states, costing energy that grows as (N−Z)²/A.

The numbers earn their keep. With a_V ≈ 15.8, a_S ≈ 18.3, a_C ≈ 0.71, a_A ≈ 23.2 MeV and δ ≈ 12/√A, the formula reproduces measured total binding energies — which range from a few MeV for light nuclei to about 1800 MeV for the heaviest — with a root-mean-square error of just a few MeV, i.e. well under 1 percent across most of the chart of nuclides. The largest residuals cluster exactly at the magic numbers, which is the model honestly telling you where its fluid analogy breaks and quantum shell structure begins. Five constants, one elegant idea, the whole periodic table of nuclei — that is the liquid drop model's enduring appeal.