Gamow's Liquid Drop Functional

• Gamow's Liquid Drop Functional is a model that represents atomic nuclei as fluid drops, balancing short-range surface tension against long-range Coulomb repulsion.
• It employs rigorous variational methods and extensions such as anisotropic, nonlocal, and fractional perimeter formulations to capture complex nuclear behaviors.
• The model is refined with phenomenological corrections that improve predictions of nuclear fission thresholds and astrophysical phase separations.

Gamow’s Liquid Drop Functional, first proposed to model the bulk energetics of atomic nuclei, encapsulates the competition between short-range surface tension forces and long-range Coulombic repulsion. This functional and its modern generalizations play a central role in nuclear theory, variational analysis, and the mathematical physics of many-body systems. The form and consequences of Gamow’s functional have been investigated from nuclear phenomenology to rigorous variational calculus, as well as in extensions incorporating nonlocal and anisotropic effects, background potentials, and high-dimensional analogues.

1. Mathematical Formulation and Physical Context

The original Gamow liquid drop energy functional models a nucleus as an incompressible fluid drop and can be written in geometric measure-theoretic terms as

$E(\Omega) = P(\Omega) + D(\Omega),$

where $\Omega \subset \mathbb{R}^3$ represents the domain occupied by nuclear matter, $P(\Omega)$ denotes the perimeter (surface area), and

$$D(\Omega) = \frac{1}{2}\iint_{\Omega \times \Omega} \frac{1}{|x-y|} \, dx\,dy$$

is the Coulomb self-energy of a uniformly charged distribution. The minimization is performed with a fixed volume constraint, $|\Omega| = V$, which is proportional to the total number of nucleons.

The functional generalizes in higher dimensions by replacing the classical perimeter and Coulomb kernel with

$$E_{\lambda,N}(\Omega) = \operatorname{Per}(\Omega) + \frac{1}{2} \iint_{\Omega \times \Omega} \frac{1}{|x-y|^\lambda} \, dx\,dy,$$

where $\lambda \in (0,N)$ is the Riesz potential exponent, and $N$ is the ambient dimension.

2. Existence and Uniqueness of Minimizers: Spherical Symmetry, Critical Mass, and Fission

For small volumes, the perimeter term dominates and the competition results in minimizers that are Euclidean balls, as established via the isoperimetric inequality and the Riesz rearrangement principle (Frank et al., 2021, Chodosh et al., 9 Jan 2024). Recent progress (Chodosh et al., 9 Jan 2024) rigorously confirms that for $|\Omega|\leq 1$ (normalized units), the unique minimizer of $E(\Omega)$ is the round ball. The proof employs new isoperimetric-type inequalities linking perimeter, Coulomb energy, and volume:

$$|\Omega|^3 < \frac{3}{16\pi} P(\Omega)^2 D(\Omega),$$

and

$$|\Omega|^2 \leq \frac{1}{12\pi} P(\Omega) D^\partial(\Omega),$$

where $D^\partial(\Omega)$ incorporates boundary-interior interactions. Equality holds only for balls, establishing rigidity.

A critical volume (or mass) threshold $V^*$ arises from energetic consideration: Above $V^*$, the repulsion outweighs surface tension and single-ball configurations energetically favor splitting, corresponding physically to nuclear fission. For the classical model with the Coulomb kernel, this threshold is computed as $m^* \approx 3.512$ (Frank et al., 2021). For $V > V^*$, minimizers do not exist in the space of connected drops; minimizing sequences break up into disjoint pieces (dichotomy), modeling the onset of fission.

3. Extension to Anisotropic and Nonlocal Energies

Recent work extends Gamow’s functional by substituting the isotropic perimeter with an anisotropic surface energy,

$$F(\Omega) = \int_{\partial^* \Omega} f(\nu_\Omega) d\mathcal{H}^{n-1}$$

where $f$ is a one-homogeneous convex function describing surface tension with preferred orientation (Choksi et al., 2018, Misiats et al., 2019). The nonlocal functional becomes

$E_f(\Omega) = F(\Omega) + V(\Omega)$

where $V(\Omega)$ is the (possibly anisotropic) Riesz potential. For smooth, elliptic $f$, the minimizer is the Wulff shape only if $f$ is isotropic; otherwise, no set other than the Euclidean ball satisfies the Euler–Lagrange equation for the Coulombic nonlocal term (Misiats et al., 2019). In the small nonlocality regime, minimizers approach the Wulff shape, with the deviation controlled in $C^1$ norm at a rate proportional to the strength of the repulsion. For crystalline surface tensions (e.g., in 2D, squares), Wulff shapes can be minimizers for small mass.

The fractional perimeter generalizes the surface term further:

$$P_s(\Omega) = \iint_{\mathbb{R}^N \times \mathbb{R}^N} \frac{|\chi_\Omega(x) - \chi_\Omega(y)|^2}{|x - y|^{N+s}} dx\,dy,$$

with $s \in (0,1)$. Nonlocal variants such as

$E_{s,g}(\Omega) = P_s(\Omega) + V_g(\Omega)$

enable modeling of long-range correlations in the interface energy (Novaga et al., 2021, Novaga et al., 2023). Existence and regularity of minimizers are preserved for kernels decaying sufficiently rapidly, with spherical symmetry recovered in the large volume limit if the kernel decay overcomes the fractional perimeter.

4. Phenomenological and Microscopic Refinements

Empirical nuclear mass formulas incorporate corrections to Gamow’s original functional, reflecting shell effects, deformation, pairing, and valence-nucleon interactions. The semi-empirical mass formula (Bethe–Weizsäcker) supplements the functional with macroscopic volume, surface, Coulomb, symmetry, and pairing energies:

$$B_{LDM} = a_v A - a_s A^{2/3} - a_c \frac{Z(Z-1)}{A^{1/3}} - a_{sym} \frac{(N-Z)^2}{A} \pm \delta,$$

where $A = N + Z$, and $\delta$ is the pairing correction.

Modern phenomenological formulas (Gangopadhyay, 2010) add shell corrections via a one-body Hamiltonian-like term,

$$B_{bunc} = \sum_{i=1,2}\sum_\alpha \epsilon_\alpha^i \mathcal{N}^i n_\alpha^i,$$

and neutron–proton valence interaction terms,

$$B_{np} = a_n N_p + a_p N_n + a_{np}^2 (N_p + N_n)^2 + a_{np}^3 (N_p + N_n)^3,$$

using 50 adjustable parameters. These corrections enable sub-MeV RMS deviation across thousands of nuclei, demonstrating that explicit modeling of shell closures and n–p correlations is essential for predictive power beyond the macroscopic LDM.

Deformation dependence is critical: Global mass formula fits show markedly lower RMS deviations for prolate deformed nuclei compared to spherical or semi-magic nuclei (Hirsch et al., 2011). The macroscopic symmetry energy surface term is highly sensitive to deformation, necessitating group-specific parameterizations and inclusion of valence nucleon and shell effects for genuine global accuracy.

5. Fission, Critical Points, and Energy Landscape Topology

Beyond absolute minimizers, the landscape of Gamow’s functional features non-minimizing volume-constrained critical points, saddles, and bifurcating branches central to modeling nuclear fission. For intermediate volumes $V \in (\alpha_0,10)$ where $\alpha_0 \approx 3.512$, mountain pass critical points exist between a ball and configurations of two balls infinitely far apart (Chambers et al., 11 Sep 2025). These represent the energy barrier to fission, with the min–max value constructed using paths in the space of Caccioppoli sets satisfying the volume constraint:

$$L = \inf_{\sigma \in \Lambda_{V_0}} \sup_{t \in [0,1]} \mathcal{E}(\sigma(t)).$$

Geometric measure theory ensures non-degenerate (multiplicity one) limits for critical sequences. The bifurcation analysis (Frank, 2019) reveals that for volumes above a threshold (e.g., $A=10$), the ball loses stability, and a smooth family of non-spherical, cylindrically symmetric critical points emerges, matching physical expectations for fission intermediates. The transition is rigorously shown to be transcritical, with stability exchanged between spherical and non-spherical branches—quantified using second variation and spherical harmonic decompositions.

6. Low-Density Limit, Microphase Separation, and Astrophysical Applications

In environments with low nuclear density and a uniform electron background, such as the outer crust of neutron stars, the functional describes microphase separation into numerous isolated droplets ("gnocchi phase") (Frank et al., 18 Jul 2025). The ground state energy per unit volume admits a precise two-term asymptotic expansion:

$$e(\rho)= \mu_*\,\rho + m_*^{2/3}\,e_{\rm Jel}\,\rho^{4/3} + o(\rho^{4/3}),$$

with $\mu_*$ the isolated droplet energy, $m_*$ the optimal droplet mass, and $e_{\rm Jel}$ the Jellium ground state energy. Here, the leading term is set by individual droplet energetics, while the first correction captures the optimal spatial arrangement and Coulomb interactions (Jellium).

7. Broader Generalizations, Open Problems, and Impact

The functional framework supports broad generalizations, including background attractive potentials restoring minimizer existence for arbitrary mass (Alama et al., 2017), fractional perimeter models for collective effects, and anisotropic or crystalline surface tensions (Choksi et al., 2018, Misiats et al., 2019). The theory informs applications in phase separation, pattern formation, and the modeling of "pasta phases" in astrophysics.

Open problems include characterizing minimizer uniqueness and existence in higher dimensions, with non-Coulombic kernels, and in regimes where the kernel's moment determines critical mass (Merlet et al., 2021). Stability analyses, energy expansions, and geometric rigidity results continue to broaden both the mathematical foundation and physical applicability of Gamow-type liquid drop models.

---

This survey reflects the modern understanding and refinement of Gamow’s Liquid Drop Functional, integrating phenomenological corrections, rigorous variational analysis, multidimensional extensions, and astrophysical relevance—a testament to the central role of combined surface tension and long-range repulsion functionals in the mathematical physics of nuclear matter.