Manna Sandpile Model: Dynamics & Universality

• The Manna sandpile model is a stochastic process where active sites redistribute particles, leading to absorbing-state transitions and avalanche phenomena.
• Its conserved and Abelian formulations highlight critical dynamics, with density serving as a control parameter that underpins hyperuniformity and anomalous transport.
• Hydrodynamic analysis links activity with diffusion and conductivity, establishing scaling relations that bridge self-organized and absorbing-state criticality.

The Manna sandpile model is a family of stochastic sandpile models in which unstable or active sites redistribute two particles randomly to nearest neighbours. Two formulations are central in the literature: the conserved Manna sandpile, also called a fixed-energy or conserved stochastic sandpile, where the total mass is strictly conserved and the control parameter is the density $\rho$; and the Abelian Manna Model, in which slow driving and boundary dissipation generate self-organized criticality. Across these formulations, the model is a canonical setting for absorbing-state phase transitions, avalanche scaling, hydrodynamic coarse graining, hyperuniformity, and universality questions involving directed percolation, conserved directed percolation, and interface depinning (Chatterjee et al., 2017, Huynh et al., 2011).

1. Microscopic definitions and principal variants

In the one-dimensional conserved Manna sandpile, one considers a periodic chain of $L$ sites, each carrying an unbounded integer mass $m_i \ge 0$ or $n_X \ge 0$. A site is active iff $m_i>1$ or equivalently $n_X \ge 2$. The continuous-time dynamics is random sequential: an active site topples at rate $1$, loses two particles, and each particle is sent independently to one of the two nearest neighbours with probability $1/2$. The total mass $N=\sum_i m_i$ is strictly conserved, and the global density is $\rho=N/L$. In steady state for $L$0, the order parameter is the activity,

$L$1

the mean density of active sites (Chatterjee et al., 2017, Mukherjee et al., 2022).

The fixed-energy formulation exhibits an active-to-absorbing transition. For $L$2 the density of active sites decays to zero and the system reaches an absorbing configuration; for $L$3 activity persists in the infinite-size limit. Any configuration with $L$4 for all $L$5 is absorbing, so the model has infinitely many absorbing states in the thermodynamic limit (Basu et al., 2012).

The Abelian Manna Model is the slowly driven counterpart. On a lattice with integer heights $L$6 or $L$7, a site is unstable if $L$8 or $L$9. When the configuration is quiescent, one particle is added at a uniformly random site; if this creates an unstable site, an avalanche begins. Each unstable site topples by removing two particles and sending them independently to randomly chosen nearest neighbours. In open geometries, particles sent to virtual neighbours are lost at the boundary. Because the toppling operators commute, $m_i \ge 0$0, the final stable configuration and the total number of topplings do not depend on the order of updates (Chen et al., 2023, Huynh et al., 2011).

This distinction between conserved and slowly driven dynamics is structurally important. In the conserved model, density is the tuning parameter of an absorbing-state transition. In the Abelian model, slow driving and dissipation organize the system into a critical state with avalanche observables such as size $m_i \ge 0$1, duration $m_i \ge 0$2 or $m_i \ge 0$3, and area $m_i \ge 0$4 (Huynh et al., 2011).

2. Hydrodynamic structure of the conserved model

A central result for conserved stochastic sandpiles is that their coarse-grained density dynamics closes through the activity. Under a local-steady-state or local equilibrium ansatz, the one-dimensional conserved Manna sandpile obeys

$m_i \ge 0$5

or, in the presence of a weak biasing field $m_i \ge 0$6,

$m_i \ge 0$7

Reading off the diffusive and drift currents yields

$m_i \ge 0$8

and

$m_i \ge 0$9

Thus the bulk diffusion coefficient is the derivative of the activity, whereas the conductivity equals the activity itself (Chatterjee et al., 2017, Tapader et al., 2020).

Macroscopic Fluctuation Theory and an additivity principle imply an equilibrium-like Einstein relation between transport and density fluctuations,

$n_X \ge 0$0

For the conserved Manna sandpile this becomes

$n_X \ge 0$1

Here $n_X \ge 0$2 is the scaled variance of subsystem mass, defined by $n_X \ge 0$3 for a macroscopic subsystem of volume $n_X \ge 0$4 (Chatterjee et al., 2017).

The same hydrodynamic structure defines an equilibrium-like chemical potential and large-deviation principle. Writing $n_X \ge 0$5 for the free-energy density, one has

$n_X \ge 0$6

so that

$n_X \ge 0$7

The coarse-grained density $n_X \ge 0$8 in a subvolume $n_X \ge 0$9 then satisfies

$m_i>1$0

This gives the conserved Manna sandpile an equilibrium-like thermodynamic description at the level of subsystem mass fluctuations, despite its nonequilibrium microscopic dynamics (Chatterjee et al., 2017).

Near the critical density $m_i>1$1, with $m_i>1$2,

$m_i>1$3

Since $m_i>1$4, one obtains

$m_i>1$5

and the dynamic exponent

$m_i>1$6

or equivalently $m_i>1$7. For the one-dimensional unrestricted-height Manna sandpile, the quoted estimates are

$m_i>1$8

while inserting $m_i>1$9 and $n_X \ge 2$0 into the scaling relation yields $n_X \ge 2$1 (Chatterjee et al., 2017).

3. Dynamic correlations, fluctuation suppression, and hyperuniformity

The conserved Manna sandpile also displays anomalous temporal correlations. Let $n_X \ge 2$2 be the cumulative bond current across bond $n_X \ge 2$3, with instantaneous current $n_X \ge 2$4. In the thermodynamic limit, for short times $n_X \ge 2$5,

$n_X \ge 2$6

with

$n_X \ge 2$7

Far from criticality, $n_X \ge 2$8 and $n_X \ge 2$9; near criticality, $1$0, so fluctuation growth is further suppressed. At very late times $1$1, ordinary diffusive growth $1$2 is recovered (Mukherjee et al., 2022).

This temporal suppression is reflected in long-time correlations and low-frequency power spectra. The current autocorrelation decays as

$1$3

Consequently, the bond-current power spectrum and subsystem-mass power spectrum behave as

$1$4

and their exponents satisfy the continuity-equation identity

$1$5

The literature characterizes the case $1$6 as dynamic hyperuniformity (Mukherjee et al., 2022).

The fluctuation relations are unusually explicit. In the long-time, large-subsystem limit,

$1$7

so the integrated current variance per unit space-time equals $1$8. Together with $1$9, this gives the exact Green-Kubo form

$1/2$0

hence again

$1/2$1

As $1/2$2, this vanishes linearly in $1/2$3, consistent with hyperuniform scaling of static density fluctuations (Mukherjee et al., 2022).

Tagged-particle transport is likewise tied directly to activity. If $1/2$4 is the displacement of a tagged particle, then

$1/2$5

and for $1/2$6 the self-diffusion coefficient defined by $1/2$7 obeys

$1/2$8

Thus $1/2$9 vanishes as $N=\sum_i m_i$0 near criticality, while the bulk diffusion coefficient $N=\sum_i m_i$1 diverges as $N=\sum_i m_i$2 (Mukherjee et al., 2022). A plausible implication is that bulk and tagged-particle transport probe different aspects of the same activity field.

4. Relaxation, anomalous transport, and front propagation

Long-wavelength density relaxation is diffusive only away from criticality. For a uniformly supercritical background with finite correlation length $N=\sum_i m_i$3, linearizing the nonlinear diffusion equation gives

$N=\sum_i m_i$4

In a finite system of size $N=\sum_i m_i$5, the slowest mode decays on a time scale of order $N=\sum_i m_i$6. Simulations of step, wedge, and Gaussian initial profiles collapse under the diffusive scaling $N=\sum_i m_i$7 (Tapader et al., 2020).

Near criticality, simple diffusion breaks down. Since

$N=\sum_i m_i$8

the bulk diffusivity diverges: $N=\sum_i m_i$9 Once the perturbation wavelength $\rho=N/L$0 becomes comparable to the correlation length $\rho=N/L$1, transport is anomalous. Finite-size or finite-wavelength scaling then yields

$\rho=N/L$2

which is the same relation as $\rho=N/L$3 (Tapader et al., 2020).

On a critical background, a localized excess-density perturbation spreads self-similarly. With $\rho=N/L$4 and $\rho=N/L$5, the asymptotic equation is

$\rho=N/L$6

A self-similar solution

$\rho=N/L$7

exists with

$\rho=N/L$8

and

$\rho=N/L$9

Hence the perturbation width grows anomalously,

$L$00

Monte Carlo data for Gaussian and delta-like initial pulses collapse onto the predicted scaling function over several decades in $L$01 (Tapader et al., 2020).

The same hydrodynamic equation describes invasion fronts between subcritical and supercritical regions. Numerically integrating $L$02 with $L$03 reproduces a moving front, the boundary-layer law

$L$04

and the front position

$L$05

For global density $L$06 on a finite ring, invaded regions freeze at $L$07; for $L$08, the whole system becomes active (Tapader et al., 2020).

5. Avalanches, finite-size scaling, and substrate correlations in the Abelian Manna Model

The slowly driven Abelian Manna Model is a standard self-organized critical sandpile. Its primary avalanche observables are the size $L$09, the duration $L$10 or $L$11, and the area $L$12. The standard finite-size-scaling ansatz is

$L$13

with cutoffs

$L$14

For $L$15,

$L$16

Moment analysis across four one-dimensional and eight two-dimensional lattices yields lattice-independent exponents within errors, supporting universality across different geometries (Huynh et al., 2011).

The quoted exponent estimates are, in $L$17,

$L$18

with

$L$19

In $L$20,

$L$21

and

$L$22

Particle conservation in the bulk implies the exact scaling relation

$L$23

equivalently $L$24, and the numerical results satisfy the narrow-joint-distribution identity

$L$25

within errors (Huynh et al., 2011).

Corrections to scaling are notable in the Manna universality class. In the moment fits, leading corrections were empirically found to be $L$26 and $L$27, and with these terms the fits reached high goodness-of-fit $L$28 for all non-trivial moments $L$29. Moment ratios

$L$30

are also lattice-independent, especially in $L$31 (Huynh et al., 2011).

Two-dimensional studies of the Abelian Manna Model further show that weak substrate correlations matter specifically for the largest avalanches. After each non-zero avalanche, one may destroy correlations by shuffling configurations across an ensemble of replicas in a way that preserves the one-point density profile but destroys two-point correlations between different $L$32-layers. In the unshuffled model, the best estimates are

$L$33

In the fully shuffled model, moment analysis yields instead

$L$34

for the large-$L$35 tail. However, the distributions show two scaling regimes: for small avalanches up to

$L$36

the exponent remains $L$37, while for large avalanches the tail crosses to $L$38 with cutoff $L$39. No single power law plus cutoff describes the entire range (Chen et al., 2023).

This establishes that destroying local substrate correlations does not alter small and intermediate avalanches, but changes the scaling of the largest, system-wide events. The result is a crossover and a breakdown of simple finite-size scaling in the decorrelated model (Chen et al., 2023).

6. Effective theories, critical densities, and universality issues

Several complementary descriptions of the Manna model have been developed. For critical absorbing states in one dimension, inverse Ising inference reconstructs an effective equilibrium Hamiltonian

$L$40

Below $L$41, the inferred interaction has a short-range repulsive core and a long-range attractive tail. At criticality, the interaction becomes purely repulsive,

$L$42

with the conjectured exact exponent

$L$43

Standard Fourier analysis then gives the window-variance exponent

$L$44

in $L$45, and the associated Gaussian field theory is

$L$46

implying

$L$47

This description links critical absorbing states to hyperuniformity and to quenched-EW-type depinning models (Machon, 2018).

At the continuum level, the stochastic field theory of the Manna or C-DP class admits an exact mapping to an overdamped elastic interface in short-range correlated quenched disorder. After a change of variables, the interface height $L$48 obeys

$L$49

where $L$50 is a retarded memory term. The mapping shows that this memory term is RG-irrelevant, so the universal large-scale dynamics is that of the quenched Edwards-Wilkinson depinning model. In this formulation, the sandpile order parameter maps to the mean interface velocity, and avalanche size is the total interface advance (Doussal et al., 2014).

The model is also central to the distinction among several density notions. In a quenched version of the Manna sandpile on directed environments, the stationary density and transition density are not equal: $L$51 while $L$52 is the unique positive solution of

$L$53

In the same setting,

$L$54

so $L$55. For the annealed Manna model, the rigorous bound $L$56 is quoted, simulations suggest $L$57, $L$58, and it is conjectured that $L$59 (Fey et al., 2012).

Universality remains contested in part of the literature. One numerical study of selected $L$60-dimensional fixed-energy sandpiles, using natural homogeneous initial states, reported

$L$61

and argued that the critical behavior converges to directed percolation after very long times (Basu et al., 2012). By contrast, the hydrodynamic treatment of conserved stochastic sandpiles derives the scaling relation

$L$62

which ordinary directed percolation in $L$63 does not satisfy: with its quoted exponents $L$64, $L$65, the relation would predict $L$66, not the measured $L$67 (Chatterjee et al., 2017). Additional numerical evidence for a distinct Manna class comes from the random rotational sandpile model, where a continuous crossover culminates at the maximally disordered point $L$68 with exponents

$L$69

matching the best-known Manna values within errors (Bhaumik et al., 2015).

Taken together, these results show that the Manna sandpile model is not only a stochastic redistribution process but also a nexus connecting absorbing-state criticality, self-organized criticality, hydrodynamic fluctuation theory, hyperuniformity, and depinning-like field theories. The principal unresolved issue is not whether the model is critical, but how its conserved dynamics should be situated relative to directed percolation, conserved directed percolation, and quenched Edwards-Wilkinson universality descriptions.