---
title: 'Zero-Range Process: Dynamics and Condensation'
url: https://www.emergentmind.com/topics/zero-range-process
type: topic
---

# Zero-Range Process: Dynamics and Condensation

A zero-range process (ZRP) is a conservative interacting particle system in which the jump rate out of a site depends only on the number of particles at that departure site. In its standard form, particles move on a lattice or graph, the local rate function \(g\) encodes the interaction, and the resulting dynamics admit a remarkably broad range of behaviors: exact stationary product measures, diffusive or hyperbolic hydrodynamic limits, condensation transitions, metastable condensate motion, non-equilibrium stationary states, and integrable multispecies generalizations [2008.06386] [2002.09214] [1507.03797].

## 1. Definition and core structure

In a general ZRP, a configuration is an occupation field \(\eta=(\eta_x)_{x\in\Lambda}\) with \(\eta_x\in\mathbb N_0\), and one particle moves from \(x\) to \(y\) at a rate determined by \(g(\eta_x)\) and a jump kernel \(p(x,y)\). A standard generator, appearing in several formulations, is
\[
(Lf)(\eta)=\sum_{x}\sum_{y} p(x,y)\,g(\eta_x)\big[f(\eta^{x,y})-f(\eta)\big],
\]
where \(\eta^{x,y}\) is the configuration after moving one particle from \(x\) to \(y\) [2002.09214]. This framework includes finite lattices, infinite lattices such as \(\mathbb Z^d\), complete graphs, open systems with reservoirs, long-jump kernels, and quenched random environments [1002.2994] [2201.07573].

A central structural assumption is often \(g(0)=0\), together with monotonicity of \(g\). When \(g\) is nondecreasing, the process is attractive, so monotone couplings and order-preservation become available; this is a key ingredient in hydrodynamic limits, strong local equilibrium, and coupling-based mixing arguments [2008.06386] [2412.16714]. The rate function can be linear, as in \(g(k)=k\); constant on occupied sites, as in \(g(k)=\mathbf 1_{\{k>0\}}\); decreasing, as in \(g(n)=(n/(n-1))^b\); or even system-size dependent, such as
\[
g_L(n)=
\begin{cases}
c,& n\le aL,\\
1,& n>aL,
\end{cases}
\qquad g_L(0)=0,
\]
which produces a discontinuous condensation transition [1507.03797] [1408.6864]. There are also constructions with rapidly growing, even superlinear, rates on \(\mathbb Z^d\), extending the standard bounded-increment regime [2004.12410]. This range of choices is one reason the ZRP functions less as a single model than as a flexible class of exactly analyzable interacting particle systems.

## 2. Stationary product measures, fugacity, and ensembles

A defining feature of many ZRPs is the existence of stationary product measures. Writing
\[
g(k)!:=\prod_{j=1}^k g(j),\qquad g(0)!:=1,
\]
the one-site grand-canonical law is typically
\[
\nu_\phi(\eta_x=k)=\frac{1}{Z(\phi)}\frac{\phi^k}{g(k)!},
\qquad
Z(\phi)=\sum_{k=0}^\infty \frac{\phi^k}{g(k)!},
\]
with product measure \(\bigotimes_x \nu_\phi\) [2002.09214]. The mean density is
\[
R(\phi)=\phi \frac{Z'(\phi)}{Z(\phi)},
\]
and when \(R\) is invertible, its inverse \(\Phi\) maps density to activity [2201.07573]. A recurrent identity is
\[
\mathbb E_{\nu_\rho}[g(\eta_x)]=\Phi(\rho),
\]
which is the microscopic origin of the nonlinear flux \(\Phi(\rho)\) in hydrodynamic equations [2002.09214].

This structure persists in several nonstandard settings. In right-biased disordered ZRPs on \(\mathbb Z\), the invariant product measure has site-dependent fugacity \(\beta/\alpha(x)\), yielding
\[
\mu_\beta^\alpha=\bigotimes_{x\in\mathbb Z}\theta_{\beta/\alpha(x)},
\]
with a maximal invariant measure at the critical fugacity \(\beta=c=\inf_x\alpha(x)\) [2008.06386]. On randomly oriented bistochastic graphs, the usual product distribution remains stationary because the divergence-free orientation balances incoming and outgoing mass in the quenched environment [2002.09214]. In boundary-driven long-jump ZRPs, the non-equilibrium stationary state is again product, but with site-dependent fugacities \(\varphi_N(x)\) solving a linear traffic equation [2201.07573].

The canonical ensemble is obtained by conditioning the product measure on fixed total mass, and ensemble equivalence or non-equivalence becomes a major question in condensation theory. In finite-capacity models, the single-site partition function is truncated,
\[
\mathcal Z_i(z)=\sum_{n=0}^{C_i} f_i(n) z^n,
\]
and the paper on finite compartments emphasizes that grand-canonical analysis remains exact even in condensed phases, unlike the infinite-capacity case where grand-canonical descriptions may fail above the critical density [1312.7251]. This suggests that the product-measure formalism is robust, but the thermodynamic interpretation of those products depends sharply on geometry, capacity constraints, and the asymptotics of \(g\).

## 3. Hydrodynamic and hydrostatic limits

The macroscopic evolution of ZRPs is controlled by the same activity function \(\Phi\) or \(\sigma\) that appears in equilibrium. For the symmetric ZRP on the discrete torus \( \mathbb T_N^d \) in dimensions \(d=1,2\), under diffusive scaling the empirical density converges to the nonlinear diffusion
\[
\partial_t f=\mathtt A:\nabla^2 \sigma(f),
\]
with \(\mathtt A\) determined by the second moment of the jump kernel [2412.16714]. A notable recent development is a quantitative consistency–stability method that yields convergence rates, uniform in time, in a modulated Monge–Kantorovich distance and in relative entropy, while avoiding block estimates [2412.16714].

Microscopic asymmetry does not by itself force hyperbolic behavior. On a two-lane torus with randomly oriented, divergence-free directed edges, a totally asymmetric and non-reversible ZRP still has a diffusive quenched hydrodynamic limit,
\[
\partial_t \rho(t,u)=\kappa\,\partial_u^2 \Phi(\rho(t,u)),
\]
because the local zero-flux constraint removes macroscopic drift [2002.09214]. The main technical novelty there is an environment-adapted tile decomposition and a second-order local equilibrium that neutralizes the usual non-gradient obstruction in entropy methods [2002.09214].

By contrast, in right-biased disordered one-dimensional ZRPs the natural scale is hyperbolic, and the macroscopic law is a scalar conservation equation
\[
\partial_t \rho(t,x)+\partial_x f(\rho(t,x))=0,
\]
with truncated flux \(f(\rho)=(p-q)c\) above the critical density \(\rho_c\) [2008.06386]. That theory extends beyond subcritical profiles: it includes hydrodynamics for supercritical data, pseudo-equilibria above \(\rho_c\), and strong local equilibrium results showing convergence to the critical measure in supercritical regions [2008.06386].

Open and nonlocal settings produce further variants. For the long-jump boundary-driven ZRP on \(\Lambda_N=\{1,\dots,N-1\}\), the hydrostatic profile is governed by regional fractional operators, and the stationary current satisfies a nonlocal Fick law in terms of the fugacity profile \(\phi_*(u)\) [2201.07573]. For the symmetric finite-segment ZRP with slow boundary creation at site \(1\) and slow annihilation at site \(N-1\), the hydrostatic limit is explicit, and the proposed diffusive hydrodynamic equation is
\[
\partial_t \rho=\Delta \Phi(\rho),
\]
with boundary conditions that depend on \(\theta\): Robin-type when \(\theta=1\), and Neumann when \(\theta>1\) [2006.13479].

## 4. Condensation and its mechanisms

Condensation is the phenomenon whereby, above a critical density, a macroscopic fraction of the mass concentrates on a vanishing fraction of sites. In the classical homogeneous setting with decreasing rates
\[
g(0)=0,\qquad g(1)=1,\qquad g(n)=\left(\frac{n}{n-1}\right)^b\quad(n\ge2),
\]
the one-site weights satisfy \(w(n)=n^{-b}\), the critical density is finite for \(b>2\), and in the supercritical regime a single site carries the excess mass \(N-\rho_c L\) while the background remains distributed as the critical product measure [1507.03797]. This is the canonical condensing ZRP on which metastability theory has been built.

The condensation mechanism, however, is not confined to decreasing rates. In a generalized conserved interacting ZRP with pairwise updates, condensation appears for a range of transfer bias \(p\), but the stationary measure does not factorize because the rates depend on both occupancies of the interacting pair [1603.05403]. In short-range dynamics the condensate-phase bump has universal scaling exponent \(\gamma^{\text{SR}}\approx1.00\), whereas in long-range mean-field dynamics the exponent is non-universal and depends on \(p\) [1603.05403]. This shows that condensation can survive even after the standard product-measure structure is broken.

Disorder can either create condensation or suppress it. In randomly perturbed ZRPs with
\[
g_x(n)=\exp\big(\sigma \xi_x(n)+b\,n^{-\gamma}\big),
\]
the homogeneous condensation criterion is altered drastically: with any \(\sigma>0\), condensation in the thermodynamic limit survives only for \(0<\gamma<1/2\), while for \(\gamma\in[1/2,1]\) the quenched critical density becomes infinite and condensation disappears [1201.1264]. By contrast, in a ZRP whose one-particle stationary distribution is multifractal, condensation is induced by quenched inhomogeneity even though the microscopic hopping is the “pure chipping” case \(u(n)=1\); the non-condensed mass scales algebraically as
\[
N_{\text{out}}(L)\sim L^{\gamma(a)},
\]
with exponent determined by the disorder strength [1609.03237].

A particularly sharp modification is the fast-rate model in which one occupation level \(A\) carries a diverging rate \(\theta_L\), while the rest of the rates are constant. In that model the critical density is finite, local marginals converge to the critical product law above \(\rho_c\), and the condensed phase depends on the scaling \(\theta_L\sim L^\gamma\): for \(\gamma\in(0,2)\), the excess mass splits into many clusters of typical size
\[
C_L=\sqrt{(\rho-\rho_c)z(1)\theta_L},
\]
whereas for \(\gamma>2\) it collapses to a single condensate site [2404.02590]. Finite capacities provide yet another mechanism: with site capacities \(C_i<\infty\), condensation can coexist with exact grand-canonical analysis and exhibit self-blocking and glassy relaxation [1312.7251]. Taken together, these results show that condensation is not equivalent to one specific tail condition on \(g\); it can arise, disappear, or change structure under disorder, capacity constraints, or even a single diverging fast rate.

## 5. Metastability, coarsening, and mixing

Once condensation occurs, the dominant dynamical questions concern motion of the condensate, coarsening toward the condensed state, and the time to equilibrium. For the one-dimensional symmetric condensing ZRP on a ring at fixed supercritical density, the condensate location is metastable on the time scale
\[
\theta_L=L^{1+b},
\]
and the rescaled condensate location converges to a symmetric pure-jump Lévy process on the unit torus, with jump kernel determined by capacities of a single random walk on the lattice [1507.03797]. The same analysis yields equilibration estimates within metastable wells, capacity asymptotics, and a uniform bound on well-exit rates [1507.03797].

In a size-dependent ZRP with
\[
g_L(n)=
\begin{cases}
c,& n\le aL,\\
1,& n>aL,
\end{cases}
\]
the free-energy landscape is explicitly metastable. The canonical free energy has competing fluid and condensed minima, and within the condensed region there is a dynamical transition line
\[
\rho_{\mathrm{dyn}}=\rho_{\mathrm{trans}}+a
\]
separating two distinct condensate-relocation mechanisms: evaporation–reformation through a fluid saddle, and splitting–coalescence through a two-condensate saddle [1408.6864]. The associated transition times satisfy Arrhenius-type asymptotics governed by the relevant large-deviation barriers [1408.6864].

Before equilibrium is reached, condensing ZRPs exhibit coarsening. On the complete graph, for
\[
u_k=1+\frac{b}{k^\sigma},\qquad 0<\sigma<1,
\]
the coarsening scale is \(t^{1/(1+\sigma)}\), and both at criticality and in the condensed phase finite-time corrections are essential to the correct scaling description [1608.01625]. At criticality the scaling function develops a front at
\[
x_0=\big[b(1+\sigma)\big]^{1/(1+\sigma)},
\]
with Gaussian smoothing at finite times; in the condensed phase the excess mass sits in a scaling bump whose amplitude approaches a selected value only slowly [1608.01625].

Mixing theory reveals a parallel picture in mean-field non-condensing models. For the unit-rate ZRP on the complete graph, with density \(m/n\to\rho>0\), the worst-case total-variation distance exhibits cutoff at
\[
n\left(\rho+\frac{1}{2}\rho^2\right),
\]
and more generally the mixing time from a given initial state depends explicitly on the largest initial occupancies [1804.04608]. For mean-field ZRPs with unbounded monotone rates \(r\), bounded density, and
\[
R(k)=\sum_{i=1}^k \frac{1}{r(i)},
\]
the mixing time is asymptotic to \(R(\hat x)\), where \(\hat x\) is the largest initial height, and the chain has cutoff at \(R(m)\) when \(R(m)\gg\log n\); the Poincaré constant is \(\Theta(1)\) in mean field and compares to the underlying graph spectral gap on general geometries [2104.10478].

## 6. Boundary-driven, multispecies, and related extensions

Several extensions show how far the ZRP framework can be pushed without losing analytic control. In open one-dimensional systems with injection at the left boundary, hopping to the right, decay in the bulk, and absorption at the target site, the arrival process can be computed exactly in the bulk-dynamics case. If each particle hops with rate \(p\) and decays with rate \(\mu\), then the asymptotic arrival rate at site \(N+1\) is
\[
\lambda_{\mathrm{arr}}=\alpha\left(\frac{p}{p+\mu}\right)^N,
\]
and the arrival counting process is Poisson in the independent-particle regime [1002.2994]. This connects ZRPs to first-passage and queueing problems as directly as to equilibrium statistical mechanics.

Multispecies and integrable variants replace scalar occupation fields by species-labeled multiplicities. In the \(n\)-species totally asymmetric ZRP on a periodic chain, smaller species have priority to hop, the process is realized as the image of an \(n\)-line process with uniform steady state, and the stationary weights admit a matrix product formula
\[
P(\sigma_1,\dots,\sigma_L)=\mathrm{Tr}\bigl(X_{\sigma_1}\cdots X_{\sigma_L}\bigr),
\]
where the \(X_{\sigma_i}\) are built from a \(q=0\) oscillator algebra and combinatorial \(R\)-matrices of \(U_q(\widehat{\mathfrak{sl}_L})\) at \(q=0\) [1511.09168]. In a different totally asymmetric special limit with \(g(n)=1\) for \(n>0\), algebraic Bethe ansatz yields explicit transition probabilities that also serve as exponential generating functions for directed random walks on multi-dimensional simplicial lattices [1703.07731].

The ZRP also appears as a dynamic environment. In a nearest-neighbor symmetric ZRP on \(\mathbb Z\) started from equilibrium, if particles on the negative axis are initially infected and healthy particles become infected upon sharing a site with an infected one, then the infection front
\[
R_t=\sup\{x\in\mathbb Z:\xi_t(x)>0\}
\]
has positive and finite asymptotic speed bounds [1808.07907]. A central ingredient is a space-time decoupling estimate derived by sprinkling, which controls correlations between trajectory functionals supported in sufficiently separated boxes [1808.07907]. This suggests that ZRPs are not only models of transport and aggregation, but also natural interacting backgrounds for reaction, spread, and propagation phenomena.

Source: https://www.emergentmind.com/topics/zero-range-process