🔢 mathematics

The Liquid Drop Model with a Yukawa Potential: Existence of Minimizers and Sharp Stability of the Ball

This paper establishes the existence of minimizers and the sharp stability of the ball for a liquid drop model perturbed by a Yukawa potential, demonstrating that the competition between surface tension and repulsion is governed by two independent parameters (screening rate and volume) and overcoming technical challenges posed by the kernel's lack of homogeneity.

Original authors: Lia Bronsard, Kenneth DeMason, Ihsan Topaloglu

Published 2026-08-21

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Lia Bronsard, Kenneth DeMason, Ihsan Topaloglu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world of atomic nuclei and the complex fluids of chemistry, matter often faces a fundamental tug-of-war. On one side, there is a force that tries to pull a collection of particles into the tightest, most compact shape possible, much like a drop of water pulling itself into a sphere to minimize its surface area. This is the pull of surface tension. On the other side, there is a repulsive force. If the particles carry an electric charge, they push away from one another, trying to spread the group out as much as possible. In the vacuum of space, this repulsion is strong and long-range, often preventing the formation of a single, stable drop if the group gets too large. However, in many real-world environments, such as the dense matter inside a neutron star or a salty solution in a laboratory, the repulsion is not infinite. It is screened, or dampened, by the surrounding medium. Imagine a charged particle trying to push away its neighbor, but a crowd of other particles stands between them, absorbing and weakening the push. This screening effect changes the rules of the game entirely, turning a simple problem of shape into a complex puzzle of stability.

A team of mathematicians has recently solved a critical piece of this puzzle, determining exactly when a cluster of matter will hold its shape as a perfect sphere and when it will break apart or deform. Their work focuses on a specific model of how these forces interact, using a mathematical description of the screened repulsion that decays exponentially. In the absence of this screening, it was long believed that if a cluster of charged matter grew beyond a certain size, no stable shape could exist at all; the repulsion would simply tear it apart. However, the researchers found that the presence of screening changes this outcome dramatically. They proved that if the screening is strong enough, stable spherical shapes can exist at any size, no matter how large the cluster becomes. This overturns the old assumption that large clusters must inevitably disintegrate.

The team also established precise boundaries for when these spherical shapes are the only possible solution. They showed that for very small clusters, the sphere is always the unique, most stable shape, regardless of how the screening is tuned. But as the cluster grows, a tipping point is reached. The researchers calculated a sharp threshold, a specific volume limit that depends on the strength of the screening. Below this limit, the sphere is not just a stable shape; it is the only shape that minimizes the energy of the system. Above this limit, the sphere becomes unstable, and the system will naturally seek out a different, non-spherical form to lower its energy. This threshold is not a vague estimate but a precise, computable value that the authors derived in a closed form, meaning it can be calculated exactly for any given level of screening.

One of the most significant aspects of this work is how it handles the mathematical difficulty of the problem. Previous attempts to understand these systems often relied on the assumption that the forces involved behaved in a uniform, scalable way. The researchers in this study had to overcome the fact that the screened repulsion does not behave uniformly; it changes character depending on the distance between particles. By developing new techniques to handle this lack of uniformity, they were able to prove the existence of stable shapes in situations where earlier methods failed. They demonstrated that for a wide range of conditions, a stable solution exists, and they identified exactly when that solution is a perfect sphere. Their findings provide a rigorous foundation for understanding the shapes of nuclear clusters in neutron stars and the behavior of charged particles in electrolytes, confirming that the competition between surface tension and screened repulsion leads to a predictable and mathematically precise set of outcomes.

Technical Summary: The Liquid Drop Model with a Yukawa Potential

Problem Statement
This paper investigates the variational problem associated with the liquid drop model of the atomic nucleus, modified by a nonlocal repulsive term defined through the Yukawa kernel. The energy functional is given by:
$F_\alpha(E) := P(E) + \int_E \int_E \frac{e^{-\alpha|x-y|}}{|x-y|} \, dx \, dy$
where $E \subset \mathbb{R}^3$ is a set of finite perimeter with fixed volume $|E| = m$ , $P(E)$ denotes the perimeter, and $\alpha > 0$ is the screening parameter. The kernel $Y_\alpha(z) = e^{-\alpha|z|}/|z|$ arises in physical models involving screened Coulomb interactions, such as charged phases in electrolytes or the inner crust of neutron stars (nuclear pasta).

The central questions addressed are the existence of global minimizers, the characterization of these minimizers (specifically whether they are balls), and the stability of the ball as a critical point under volume constraints. Unlike the classical unscreened case ( $\alpha=0$ , Gamow's model), where minimizers are conjectured to exist only below a critical mass and fail to exist above it, the presence of the exponential decay in the Yukawa kernel introduces a new length scale ( $\ell = 1/\alpha$ ) that qualitatively alters the competition between surface tension and repulsion.

Methodology
The authors employ a combination of geometric measure theory, spectral analysis, and variational techniques. Key methodological components include:

Boundary Representation: Following the approach of Muratov, Novaga, and Simon, the nonlocal term is rewritten as a boundary integral. This isolates an effective surface tension term and a non-negative defect term, allowing the authors to analyze the energy scaling and subadditivity properties.
Scaling and Homogeneity Analysis: The paper addresses the lack of homogeneity in the Yukawa kernel (unlike the Riesz interaction in the $\alpha=0$ case). By utilizing specific scaling arguments and comparing the energy of a set to its scaled versions, the authors derive strict binding inequalities necessary for proving the existence of minimizers.
Spectral Analysis of the Second Variation: To determine the stability of the ball, the authors expand the Yukawa potential in spherical harmonics. They compute the eigenvalues of the associated integral operator on the sphere, relating them to modified Bessel functions. This allows for the explicit calculation of the second variation of the energy functional.
Selection Principle and Regularity: To establish local minimality, the authors use a selection principle argument. They construct a sequence of nearly spherical sets with better regularity properties, allowing the application of quantitative stability estimates derived from the second variation.

Key Contributions and Results

Existence of Minimizers for All Volumes (High Screening):
Theorem 1.1 establishes that if the screening rate $\alpha$ exceeds an explicit critical threshold, $\alpha > (2\pi)^{1/3}$ , minimizers exist for every volume $m > 0$ . This contrasts with previous results requiring large volumes and resolves the non-existence issue found in the unscreened model for large masses. Furthermore, the map $m \mapsto m^{-1}e_\alpha(m)$ is shown to be strictly decreasing.
Existence of Minimizers for Small Volumes (Any Screening):
Theorem 1.2 proves that for any screening rate $\alpha \geq 0$ , minimizers exist provided the volume $m$ is sufficiently small (specifically $m < \bar{m}$ , where $\bar{m}$ is an explicit constant derived from the unscreened case). This extends the existence results of Frank and Nam to the screened case.
Uniqueness of the Ball for Small Volumes:
Theorem 1.3 demonstrates that for sufficiently small volumes (uniformly in $\alpha$ ), the unique global minimizer (up to translation) is the ball. This is proven by showing that for small perturbations of the ball, the gain in perimeter energy outweighs the loss in nonlocal energy, utilizing the positive definiteness of the Yukawa kernel to simplify the deficit estimates.
Sharp Stability Threshold:
Theorem 1.4 provides a sharp, computable volume threshold $m^*(\alpha)$ for the stability of the ball as a volume-constrained critical point.
- If $\alpha \geq (2\pi)^{1/3}$ , the ball is stable for all volumes ( $m^*(\alpha) = +\infty$ ).
- If $\alpha < (2\pi)^{1/3}$ , the ball is stable if and only if $m \leq m^*(\alpha)$ .
  The threshold is given in closed form via the ratio of eigenvalues of the second variation operator.
Local Minimality and Rigidity:
Theorem 1.5 establishes the $L^1$ -local minimality of the ball below the stability threshold. It proves that if $m < m^*(\alpha)$ , the ball is a local minimizer, and if $m > m^*(\alpha)$ , the ball is unstable (there exists a sequence of sets converging to the ball in $L^1$ with strictly lower energy).

Significance and Claims
The paper claims to resolve the existence and stability questions for the Yukawa liquid drop model across the full range of parameters, overcoming technical difficulties posed by the lack of homogeneity in the kernel.

Qualitative Shift: The results highlight a qualitative shift from the unscreened model: screening prevents the "blow-up" or non-existence of minimizers at large masses when the screening is strong enough.
Sharpness: The stability threshold $m^*(\alpha)$ is derived in closed form and is shown to be sharp, recovering the known unscreened value as $\alpha \to 0$ .
Generality: The proofs are designed to be robust, relying on the specific properties of the Yukawa kernel (such as positive definiteness and exponential decay) rather than homogeneity, suggesting the potential applicability of these techniques to other nonlocal functionals with similar decay properties.

The authors note that while the results are stated for dimension $n=3$ , similar results are expected to hold in higher dimensions, though the kernel analysis is most established in the three-dimensional case. The work connects the mathematical analysis of these functionals to physical phenomena like nuclear pasta and colloidal stability, providing rigorous justification for the stability of spherical clusters under screened repulsion.

More like this