[Paper Review] Classifying Minimum Energy States for Interacting Particles (I) -Spherical Shells

Cameron Davies, Tongseok Lim+1 more|arXiv (Cornell University)|Jan 1, 2021

Mathematical Biology Tumor Growth46 references2 citations

TL;DR

This paper classifies minimum energy states for particles on spherical shells with long-range attraction and short-range repulsion governed by power-law potentials. It proves that for attraction exponent α ≥ α_Δⁿ(β) and repulsion β ≥ 2, equidistribution on a regular n-simplex minimizes energy, with unique minimizers up to rigid motions when strict inequality holds; at (α,β) = (2,4), minimizers are characterized by matching first and second moments with the spherical shell, and nonlinear stability is established via d_α-Lyapunov stability without linearization.

ABSTRACT

Particles interacting through long-range attraction and short-range repulsion given by power-laws have been widely used to model physical and biological systems, and to predict or explain many of the patterns they display. Apart from rare values of the attractive and repulsive exponents $(\al,\bt)$, the energy minimizing configurations of particles are not explicitly known, although simulations and local stability considerations have led to conjectures with strong evidence over a much wider region of parameters. For a segment $\bt=2 \bt\ge2$, a unimodal threshold $2<\al_{\Delta^n}(\bt) \le \max\{\bt,4\}$ exists such that equidistribution of particles over a unit diameter regular $n$-simplex minimizes the energy if and only if $\al \ge \al_{\Delta^n}(\bt)$ (and minimizes uniquely up to rigid motions if strict inequality holds). At the point $(\al,\bt)=(2,4)$ separating these regimes, we show the minimizers all lie on a sphere and are precisely characterized by sharing all first and second moments with the spherical shell. Although the minimizers need not be asymptotically stable, our approach establishes $d_\al$-Lyapunov nonlinear stability of the associated ($d_2$-gradient) aggregation dynamics near the minimizer in both of these adjacent regimes -- without reference to linearization. The $L^\al$-Kantorovich-Rubinstein distance $d_\al$ which quantifies stability is chosen to match the attraction exponent.

Motivation & Objective

• To identify energy-minimizing configurations of particles on spherical shells under long-range attraction and short-range repulsion.
• To determine the critical threshold α_Δⁿ(β) that separates regimes of equidistribution versus non-equidistributed minimizers.
• To establish nonlinear Lyapunov stability for aggregation dynamics near minimizers without relying on linearization.
• To characterize minimizers at the critical point (α,β) = (2,4) via moment matching with the spherical shell.
• To unify the analysis of energy minimization and dynamical stability across adjacent parameter regimes.

Proposed method

• The analysis uses the L^α-Kantorovich-Rubinstein distance d_α to quantify stability, matching the attraction exponent α.
• A unimodal threshold α_Δⁿ(β) is derived for β ≥ 2, determining when equidistribution on a regular n-simplex minimizes energy.
• The paper employs variational methods and moment constraints to characterize minimizers at (α,β) = (2,4), showing they match the spherical shell in first and second moments.
• Nonlinear d_α-Lyapunov stability is proven for the d_2-gradient aggregation dynamics near minimizers, using energy-based arguments without linearization.
• The method leverages symmetry and geometric constraints to establish uniqueness and stability across parameter regimes.

Experimental results

Research questions

• RQ1 For which values of the attraction exponent α and repulsion exponent β is equidistribution on a regular n-simplex the unique energy minimizer on a spherical shell?
• RQ2 What characterizes the energy minimizers at the critical point (α,β) = (2,4), and how do they relate to the spherical shell?
• RQ3 How can nonlinear stability of aggregation dynamics be established near minimizers without linearization?
• RQ4 What is the role of the L^α-Kantorovich-Rubinstein distance in quantifying stability for power-law interactions?
• RQ5 How does the threshold α_Δⁿ(β) govern the transition between equidistributed and non-equidistributed minimizers?

Key findings

• For β ≥ 2 and α ≥ α_Δⁿ(β), the energy is minimized by equidistribution on a regular n-simplex, uniquely up to rigid motions when α > α_Δⁿ(β).
• At the critical point (α,β) = (2,4), all minimizers share the same first and second moments as the spherical shell, providing a precise geometric characterization.
• The d_α-Lyapunov stability of the d_2-gradient aggregation dynamics is established near minimizers in both regimes adjacent to (2,4), without linearization.
• The threshold α_Δⁿ(β) satisfies 2 < α_Δⁿ(β) ≤ max{β, 4}, defining the boundary between equidistribution and non-equidistribution in energy minimization.
• The stability analysis is intrinsic to the interaction potential, using the L^α-Kantorovich-Rubinstein distance that matches the attraction exponent α.