Static Equilibria of Perturbed Spheres: A Single-Harmonic Class Map, a Parity Obstruction, and a Certified Counter for the Mono-Monostatic Regime

Technical Summary: Static Equilibria of Perturbed Spheres

Problem Statement
The paper addresses the classical problem of determining the number and stability types (stable $S$, unstable $U$, and saddle $H$) of static equilibria for a homogeneous convex body resting on a horizontal plane under gravity. While Várkonyi and Domokos [1, 2] proved the existence of bodies for any class $(S, U)$ with $S, U \ge 1$ (including the mono-monostatic Gömböc where $S=U=1$), their result does not provide an elementary closed-form description of which class a specific shape realizes. The paper focuses on the simplest nondegenerate perturbations of a sphere: homogeneous bodies whose boundaries are defined by a unit sphere radially perturbed by a single real spherical harmonic $Y^m_\ell$ of degree $\ell \ge 2$ and order $1 \le m \le \ell$. The central question is to determine the exact equilibrium class $(S, U)$ for these single-harmonic shapes.

Methodology
The approach combines analytical geometry, spherical harmonic theory, and certified numerical computation:

1. Reduction to Critical Points: The author utilizes the established reduction that static equilibria correspond to the critical points of the squared centroid-to-surface distance function, $\rho(u) = |x(u) - c|^2$.
2. Centroid Analysis (Lemma 1): A key analytical step involves determining the location of the centroid $c$. The paper proves that for a single tesseral harmonic $Y^m_\ell$ with $m \ge 2$, the symmetry of the perturbation forces the centroid to remain exactly at the origin ($c=0$) for all amplitudes within the convex range. For $m=1$, the centroid shift is $O(\epsilon^3)$ and vanishes only under specific parity conditions.
3. Exact Reduction (Lemma 2): Because $c=0$ for $m \ge 2$, the distance function $\rho(u)$ becomes exactly $(1 + \epsilon g(u))^2$. Since this is a strictly increasing function of the harmonic $g$, the critical points of the physical body are in exact bijection with the critical points of the harmonic $Y^m_\ell$ itself, preserving stability types and indices. This holds for the entire convex range, not just asymptotically.
4. Critical Point Enumeration (Proposition 1): The author counts the critical points of $Y^m_\ell$ by exploiting its separable structure ($P^m_\ell(\cos \theta) \cos(m\phi)$). They analyze the zeros of the associated Legendre functions and their derivatives to classify local maxima, minima, and saddles, including the behavior at the poles (which act as degenerate "monkey saddles" for $m \ge 3$).
5. Certified Numerical Counter (Section 7): For the multi-harmonic regime (where the single-harmonic reduction fails and the centroid shifts), the author develops a certified equilibrium counter. This method uses interval arithmetic to enclose the centroid-to-surface distance, applies the Krawczyk operator for existence and uniqueness certification, and utilizes interval Hessians for stability classification. Stereographic charts are used to handle polar singularities.

Key Results

• The Equilibrium-Class Theorem (Theorem 1): For a body perturbed by a single harmonic $Y^m_\ell$ ($\ell \ge 2, 1 \le m \le \ell$), the number of stable and unstable equilibria is exactly:
  $$S = U = m(\ell - m + 1)$$
  The remaining critical points are $2m(\ell-m)$ ordinary saddles and two poles of index $1-m$. The Poincaré–Hopf balance $S - H + U = 2$ is satisfied exactly.
• Diagonal Reachability (Corollary 1): Single-harmonic perturbations can only realize equilibrium classes on the diagonal $S=U$. The reachable values are $\{2, 3, 4, \dots\}$.
• No Single-Harmonic Mono-monostaticity (Corollary 2): Since $m(\ell - m + 1) = 1$ implies $\ell=1$ (a rigid translation, not a shape change), no single spherical harmonic of degree $\ell \ge 2$ can produce a mono-monostatic body ($S=U=1$).
• Parity Obstruction (Corollary 3): Any perturbation composed solely of even-degree harmonics results in a centrally symmetric body. Such bodies have equilibria in antipodal pairs, forcing $S$ and $U$ to be even numbers. Consequently, centrally symmetric bodies cannot be mono-monostatic.
• Certified Verification: The author validated their analytical formula against numerical simulations for eight specific single-harmonic bodies, achieving exact agreement. Furthermore, the certified counter was applied to a known analytic parameterization of a near-spherical body (previously ambiguous in literature), rigorously confirming it as mono-monostatic with $(S, U, H) = (1, 1, 0)$.

Significance and Claims
The paper claims to provide a "complete constructive answer" to the equilibrium classification problem for the simplest nondegenerate shapes, contrasting with the non-constructive existence proofs of Várkonyi and Domokos.

1. Structural Insight: It establishes that mono-monostaticity is an intrinsically multi-harmonic phenomenon requiring odd-parity content; it cannot be achieved by a single mode or by centrally symmetric shapes.
2. Exactness vs. Approximation: The paper highlights that for $m \ge 2$, the reduction from the physical body to the harmonic is an exact identity, not a leading-order approximation. This allows for a closed-form count without asymptotic error terms.
3. Resolution of Ambiguity: The paper addresses the practical difficulty of counting equilibria in near-spherical bodies, where seed-based methods often undercount and grid-based methods overcount. The proposed certified counter removes this ambiguity, providing a provably correct count that neither under- nor over-estimates, thereby settling specific questions regarding analytic parameterizations that previous methods left ambiguous.

The work is presented as a rigorous mapping of the "ends" of the problem: an exact closed-form solution where a single mode dominates, and a certified decision procedure for the complex multi-harmonic regime where the Gömböc class resides.