Sandpile solitons via smoothing of superharmonic functions
Nikita Kalinin, Mikhail Shkolnikov
TL;DR
The paper develops a smoothing theory for integer-valued superharmonic functions on $Z^2$ to construct canonical, minimal refinements $\theta_F$ of a piecewise-linear $F$ and uses these refinements to model sandpile solitons. It proves the existence (and uniqueness up to translation) of solitons for all rational slopes, and shows how solitons, triads, and nodes arise from the interaction of waves under a Least Action Principle on infinite domains. The analysis leverages reductions to cylinders, monotonicity properties, and discrete potential-theoretic tools to establish stabilization of the smoothing sequences for base functions $\Psi_{edge}$, $\Psi_{vertex}$, and $\Psi_{node}$, yielding explicit soliton profiles of the form $3+\Delta\theta_{\min(p'x+q'y,0)}$. These results connect discrete sandpile dynamics with tropical geometry and fractal pattern formation, providing a rigorous bridge between combinatorial smoothing and observed self-reproducing sandpile patterns.
Abstract
Let $F:\mathbb Z^2\to \mathbb Z$ be the pointwise minimum of several linear functions. The theory of smoothing of integer-valued superharmonic function allows us to prove that under certain conditions there exists the pointwise minimal superharmonic function which coincides with $F$ "at infinity". We develop such a theory to prove the existence of so-called solitons (or strings) in a certain sandpile model, studied by S. Caracciolo, G. Paoletti, and A. Sportiello. Thus we made a step towards understanding the phenomena of the identity in the sandpile group for a square where solitons appear according to experiments. We prove that sandpile states, defined using our smoothing procedure, move changeless when we send waves (that is why we call them solitons), and can interact, forming triads and nodes.
