Liquid Drop Model: Nuclear & Continuum Physics

• The liquid drop model is a macroscopic framework that treats nuclei as nearly incompressible fluid drops with short-range cohesive forces and long-range Coulomb repulsion.
• It underpins the semi-empirical mass formula and fission analyses by balancing surface energy with nonlocal repulsive interactions to explain nuclear size and stability.
• Extensions of the model to anisotropic, continuum, and exotic matter settings broaden its application, bridging nuclear physics with capillary and interface phenomenon studies.

The liquid drop model is a macroscopic description in which a finite body is treated as an incompressible or nearly incompressible drop with short-range cohesion, surface tension, and longer-range destabilizing interactions. In nuclear physics, it originated in the work of Gamow and of Bohr–Wheeler, and it remains the standard framework for describing smooth trends in nuclear binding, Coulomb-driven deformation, fission, and related variational shape problems. In modern usage, the same label also covers mathematically precise nonlocal shape functionals and several continuum free-surface models, but the nuclear case remains the canonical reference (Pino et al., 2024, Harsha, 2017, Frank et al., 2021).

1. Historical concept and physical content

In its classical nuclear form, the liquid drop model treats the nucleus as a tiny charged drop of quantum incompressible fluid. The underlying physical picture is that nuclear forces are short-ranged and saturating, so each nucleon interacts strongly only with near neighbors; as a result, the bulk binding grows roughly linearly with the nucleon number $A$ rather than quadratically. Approximate incompressibility implies an almost constant density and the familiar radius scaling $R \sim r_0 A^{1/3}$. Surface nucleons have fewer neighbors and are therefore less bound, while protons experience long-range Coulomb repulsion. These ingredients yield a macroscopic theory of size, stability, and deformation that does not resolve single-particle shell structure but captures bulk nuclear systematics (Harsha, 2017).

A second, closely related formulation treats the liquid drop model as a shape-optimization problem. There the nucleus is represented by a set $\Omega \subset \mathbb{R}^3$ of fixed volume, and the energy is written as a sum of surface area and Coulomb self-interaction. This variational form makes explicit the competition between a local compactness term and a nonlocal repulsive term, and it is the standard starting point for rigorous existence, nonexistence, bifurcation, and fission analyses (Frank et al., 2021, Pino et al., 2024).

2. Semi-empirical mass formula and nuclear binding systematics

The standard phenomenological realization of the model is the semi-empirical mass formula. In Bethe–Weizsäcker form, the binding energy is decomposed into volume, surface, Coulomb, symmetry, and pairing contributions,

$$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$

The volume term encodes saturation, the surface term penalizes small nuclei, the Coulomb term suppresses large proton-rich nuclei, the symmetry term penalizes neutron–proton imbalance, and the pairing term distinguishes even–even, odd-$A$, and odd–odd systems (Chaudhuri et al., 2020, Harsha, 2017).

A representative five-term fit used in one analysis gives

$$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$

with the Coulomb term written as $-a_c Z(Z-1)/A^{1/3}$ and the pairing term proportional to $A^{-1/2}$ (Harsha, 2017). Treating $A$ and $Z$ as continuous variables, analytic maximization of $R \sim r_0 A^{1/3}$0 yields local extrema at $R \sim r_0 A^{1/3}$1 and at $R \sim r_0 A^{1/3}$2, corresponding to the neighborhood of $R \sim r_0 A^{1/3}$3 and $R \sim r_0 A^{1/3}$4 (Harsha, 2017).

The model also clarifies a persistent misconception. In the smooth liquid-drop approximation, the binding-energy-per-nucleon curve has maxima near $R \sim r_0 A^{1/3}$5 and $R \sim r_0 A^{1/3}$6, but shell corrections are required to recover the experimental ordering $R \sim r_0 A^{1/3}$7. The high abundance of $R \sim r_0 A^{1/3}$8 is instead tied to stellar nucleosynthesis and nuclear statistical equilibrium, not to it being the global maximum of $R \sim r_0 A^{1/3}$9 (Harsha, 2017).

3. Variational Gamow model and equilibrium shapes

The mathematically precise liquid drop model writes the energy of a set $\Omega \subset \mathbb{R}^3$0 of fixed mass $\Omega \subset \mathbb{R}^3$1 as

$\Omega \subset \mathbb{R}^3$2

with the Coulombic case corresponding to $\Omega \subset \mathbb{R}^3$3 and $\Omega \subset \mathbb{R}^3$4. In this formulation, the perimeter favors compact shapes, while the Riesz term favors spreading and ultimately breakup. Under the scaling $\Omega \subset \mathbb{R}^3$5, the perimeter scales like $\Omega \subset \mathbb{R}^3$6 and the Riesz term like $\Omega \subset \mathbb{R}^3$7, so small masses are surface dominated and large masses are repulsion dominated (Frank et al., 2021, Choksi et al., 2018).

This competition produces a critical-mass structure. For the general Riesz model, there is a natural threshold $\Omega \subset \mathbb{R}^3$8 defined by equality of the energy of one ball and two infinitely separated half-mass balls. Existence of minimizers has been proved for all $\Omega \subset \mathbb{R}^3$9, and, conditionally on nonexistence for $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$0, balls are minimizers for $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$1 and unique minimizers for $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$2. For $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$3, there is nonexistence of minimizers for sufficiently large mass, reflecting fission into far-separated droplets (Frank et al., 2021).

Volume-constrained critical points satisfy the Euler–Lagrange equation

$$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$4

where $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$5 is the mean curvature and $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$6 is the Lagrange multiplier associated with the fixed-volume constraint (Pino et al., 2024). Balls always solve this equation. They are global minimizers for sufficiently small mass, but the set of critical points is much richer. In particular, there is a bifurcation at volume $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$7: a branch of axisymmetric but non-spherical equilibria emerges from the ball, with prolate-like shapes below the bifurcation and oblate-like shapes above it, together with an exchange of stability in the restricted axisymmetric class (Frank, 2019).

4. Fission, fragmentation, and radioactive decay

The model is central to macroscopic descriptions of fission. A recent phase-field treatment of the classical liquid drop functional computes full minimum-energy fission paths and numerically tracks the Bohr–Wheeler branch of non-spherical equilibria. In that formulation, the fissility parameter $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$8 controls the Coulomb strength, the ball loses local stability at $$B(A,Z)=a_v A-a_s A^{2/3}-a_c\frac{Z^2}{A^{1/3}}-a_{\text{sym}}\frac{(A-2Z)^2}{A}+\delta_{\text{pair}}.$$9, and a Businaro–Gallone point appears at $A$0, where the symmetric saddle changes index. The same study continues the Bohr–Wheeler branch over $A$1 with $A$2, and numerically confirms that as $A$3 the saddle approaches two touching equal spheres (Xu et al., 2023).

At finite temperature, liquid-drop input is used inside fragmentation models. In the canonical thermodynamical model, the fragment free energy includes a temperature-dependent surface coefficient $A$4, a Coulomb term with Wigner–Seitz correction, a symmetry term, and a Fermi-gas excitation contribution. Within that framework, the liquid–gas transition temperature is highly sensitive to the surface term: varying $A$5 from $A$6 to $A$7 MeV shifts the transition temperature by about $A$8 MeV, while turning off the temperature dependence of $A$9 shifts it by about $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$0 MeV. By contrast, varying $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$1 between $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$2 and $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$3 MeV or $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$4 over $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$5–$$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$6 MeV produces no significant shift in the transition temperature (Chaudhuri et al., 2020).

Generalized liquid-drop models also underpin barrier-penetration descriptions of charged-particle emission. In one extension to two-proton radioactivity, the macroscopic energy is written as

$$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$7

and WKB penetrability is combined with a spectroscopic factor to predict half-lives; the resulting model describes known ground-state $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$8 emitters well and has comparable accuracy to other approaches (Cui et al., 2020). For proton radioactivity, an effective liquid drop model with varying mass asymmetry shape and an effective inertial coefficient reproduces the available data with all deviations in $$a_v=15.7827\ \text{MeV},\quad a_s=17.9042\ \text{MeV},\quad a_c=0.724040\ \text{MeV},\quad a_a=23.7193\ \text{MeV},\quad a_p=11.0000\ \text{MeV},$$9 below $-a_c Z(Z-1)/A^{1/3}$0 and root-mean-square deviation $-a_c Z(Z-1)/A^{1/3}$1 (Sheng et al., 2014).

5. Three-flavor and exotic-matter extensions

The liquid-drop logic has also been extended beyond ordinary two-flavor nuclear matter. A three-flavor analogue, formulated for strangeon matter, assumes that strangeons are localized clusters containing $-a_c Z(Z-1)/A^{1/3}$2 quarks and that their interaction is short-ranged and saturating. Because three-flavor symmetry restoration makes asymmetry energy negligible and Coulomb effects are taken to be small, the strangeon-drop energy per baryon is modeled as

$-a_c Z(Z-1)/A^{1/3}$3

where $-a_c Z(Z-1)/A^{1/3}$4 is the strangeon mass per baryon in vacuum and the bulk and surface coefficients are derived phenomenologically from a Lennard–Jones-like interaction via a corresponding-states argument (Wang et al., 2017).

In that model, stability is governed by the competition between the bulk term and the positive surface correction. For $-a_c Z(Z-1)/A^{1/3}$5 GeV and inter-strangeon potential depth $-a_c Z(Z-1)/A^{1/3}$6 MeV, strangeon matter could be more stable than two-flavor nuclear matter, and the critical baryon number can be as low as $-a_c Z(Z-1)/A^{1/3}$7 (Wang et al., 2017). The authors connect this to possible strangeon stars, strangeon nuggets ejected in mergers, and relic strangeon matter synthesized during the QCD phase transition; these are presented as phenomenological implications of the assumed three-flavor liquid-drop stability window (Wang et al., 2017).

6. Mathematical generalizations and nonclassical equilibria

A major modern development is the extension of the model to anisotropic or externally confined settings. In anisotropic liquid-drop models, the isotropic perimeter is replaced by an anisotropic surface energy

$-a_c Z(Z-1)/A^{1/3}$8

while the repulsive term remains Riesz-type. For smooth elliptic surface tensions, Wulff shapes are critical points only in the isotropic case; for genuinely anisotropic smooth tensions they are not even critical. By contrast, for certain crystalline tensions the Wulff shape becomes the unique minimizer in the small-mass regime (Choksi et al., 2018). This sharply separates smooth elliptic anisotropy from crystalline anisotropy.

Another regularization introduces a background attraction,

$-a_c Z(Z-1)/A^{1/3}$9

The added long-range attractive potential restores existence of minimizers for arbitrary mass. In the small-$A^{-1/2}$0 regime, minimizers decompose into finitely many well-separated droplets, each a minimizer of the unconfined problem for its own mass, with separation scale $A^{-1/2}$1; this is formalized through the notion of generalized minimizers (Alama et al., 2017).

The spectrum of non-spherical compact equilibria is broader still. Besides the quadrupolar branch bifurcating from the ball, a recent construction proves the existence of large-volume embedded solutions whose geometry resembles a pearl necklace wrapped around a large circle and is close to a Delaunay unduloid. These are compact genus-one equilibria of the full liquid-drop Euler–Lagrange equation, and they have no analogue in the pure constant-mean-curvature problem, where Alexandrov’s theorem leaves only spheres (Pino et al., 2024).

7. Related continuum uses outside nuclear physics

Outside nuclear theory, the same label also denotes capillary drop models based on surface energy, effective interface potentials, and gradient dynamics. In one mesoscopic framework derived from density functional theory, the excess free energy of a film of thickness $A^{-1/2}$2 is

$A^{-1/2}$3

where $A^{-1/2}$4 is the binding potential. The equilibrium contact angle then satisfies

$A^{-1/2}$5

with $A^{-1/2}$6 the equilibrium film thickness. This formulation provides a direct link between microscopic interactions, wetting behavior, and mesoscopic drop profiles (Hughes et al., 2015).

Other continuum liquid-drop models study different free-boundary problems but retain the same energy-balance logic. A one-dimensional quasi-static model for a drop sliding down an inclined plane couples an elliptic capillary equilibrium problem to contact-line motion laws and admits existence, uniqueness, traveling waves, and homogenization results (Kim et al., 2012). An interface-formation model for coalescence removes the singularity of the conventional instantaneous-bridge description and predicts a non-monotone bridge-expansion speed for sufficiently small drops (Sprittles et al., 2014). A recent three-field gradient-dynamics model for volatile drops on thin porous substrates couples the drop height, vapor density, and substrate saturation through a single free-energy functional, thereby describing spreading, evaporation, imbibition, and the formation of a saturation halo (Hartmann et al., 2024).

Across these uses, the recurring structure is the same: a local interfacial term favors compactness, while additional nonlocal, entropic, or external contributions determine preferred size, shape, stability, and breakup pathways.