Kinetically-Constrained Hopping Process

• Kinetically-Constrained Hopping Process (KCHP) is a stochastic particle model where hopping is permitted only when a prescribed local pattern of vacancies is present.
• It serves as a lattice-based transport theory that illustrates how local kinetic constraints generate nontrivial, density-dependent diffusion and cooperative dynamics.
• Recent analyses using classical-to-quantum mapping and hydrodynamic limits reveal cases with exact transport coefficients despite non-gradient current structures.

Searching arXiv for recent and foundational papers on kinetically-constrained hopping processes and related kinetically constrained lattice gases. A kinetically-constrained hopping process (KCHP) is a stochastic interacting-particle system in which particle motion is governed not only by hard-core exclusion but also by local kinetic constraints that depend on the surrounding configuration. In the explicit usage of the term, a representative KCHP is the triangular-chain or triangular-lattice model in which a particle can move only if the other two sites of a triangle are unoccupied; more broadly, the literature treats closely related kinetically constrained lattice gases, especially the Kob–Andersen model, as canonical realizations of the same constrained-transport mechanism (McRoberts et al., 23 Jul 2025, Arita et al., 2017). Within the wider theory of kinetically constrained models, KCHP is best understood as the conservative, hopping-based analogue of vacancy-facilitated dynamics, rather than as a synonym for all kinetically constrained models (Hartarsky et al., 2024).

1. Conceptual definition and placement within kinetically constrained dynamics

The defining feature of KCHP is that hopping is conditional. A move is legal only when a prescribed local pattern of vacancies is present, so the admissibility of a hop depends on more than the occupancy of the target site. This distinguishes KCHP from unconstrained exclusion processes, in which the only local requirement is that the destination be empty.

This local-facilitation structure places KCHP within the broader family of kinetically constrained dynamics studied in glassy systems. The general framework of kinetically constrained models uses binary occupation variables on a lattice and imposes update rules that require a suitable neighborhood to be empty before a local event can occur. The 2024 monograph on kinetically constrained models states explicitly that the acronym KCHP is not used there, and identifies kinetically constrained lattice gases and Kob–Andersen-type conservative models as the closest direct analogues of hopping processes (Hartarsky et al., 2024). A common misconception is therefore to identify KCHP with the more widely studied spin-flip KCMs; the relation is closer and more specific: KCHP corresponds to the conservative, motion-based side of the same facilitation principle.

Another recurrent misconception is that simple equilibrium statistics should imply simple transport. The KCHP literature shows the opposite: trivial Bernoulli equilibrium can coexist with strongly nontrivial density-dependent diffusion, non-gradient microscopic currents, and severe transport suppression at high density (Arita et al., 2017, Teomy et al., 2016).

2. Microscopic realizations and local constraint rules

In the triangular-chain and triangular-lattice KCHPs, the microscopic degrees of freedom are hard-core occupation variables, with $\sigma_i=1$ for an occupied site and $\sigma_i=0$ for a hole. A triangle is selected at random, and if that triangle contains exactly one particle, that particle hops to one of the other two sites with equal probability. Equivalently, a particle can move only if the other two sites of the triangle are unoccupied. The constraint is therefore a three-site gate: hopping between two sites is conditioned on the occupancy of a third site (McRoberts et al., 23 Jul 2025).

A closely related and foundational realization is the two-dimensional Kob–Andersen kinetically constrained lattice gas on the square lattice. Here particles obey hard-core exclusion and attempt symmetric nearest-neighbor hops, but a move is permitted only if the departure site has at least two empty neighbors and the target site also has at least two empty neighbors. On the square lattice, the relevant nontrivial choice is $m=2$. The instantaneous rightward current can be written as

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$

where $H_{i,j}(\tau)$ is a binary kinetic constraint factor encoding whether the move is blocked or allowed (Arita et al., 2017).

The same constrained-hopping idea appears in the hydrodynamic study of kinetically constrained lattice gases driven by reservoirs. There the occupancy variable is

$n_\alpha(\vec r,t)\in\{0,1\},$

and the averaged density is

$\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$

For the two-dimensional $m=2$ Kob–Andersen model, the allowed-move indicator $K^{KA}_{\alpha,\hat d(\vec r,t)}$ depends on the three neighbors of the origin and the three neighbors of the target site, again expressing the requirement that sufficient local vacancies exist both before and after the hop (Teomy et al., 2016).

These models are representative rather than exhaustive. They illustrate the two central microscopic ingredients of KCHP: conservative particle motion and configuration-dependent mobility.

3. Equilibrium structure, conserved density, and non-gradient character

A notable structural property of the canonical lattice realizations is that equilibrium is simple despite constrained dynamics. In the square-lattice Kob–Andersen model, detailed balance implies that the stationary equilibrium measure is a product Bernoulli measure at density $\rho$. For distinct sites,

$\sigma_i=0$0

and the compressibility is

$\sigma_i=0$1

Thus the equal-time static measure has no nontrivial spatial correlations even though the dynamics is highly cooperative (Arita et al., 2017).

The same theme appears in the broader hydrodynamic treatment of kinetically constrained lattice gases. Their equilibrium state is trivial in the sense that there are no equal-time correlations between occupancies at different sites, but transport under weak driving still develops nontrivial structure. The literature distinguishes the bulk diffusion coefficient $\sigma_i=0$2, defined from transport under a density gradient, from the self-diffusion coefficient

$\sigma_i=0$3

defined from the long-time mean-squared displacement of a tagged particle in equilibrium; these are generally very different objects in kinetically constrained systems (Teomy et al., 2016).

The principal analytical difficulty is that these models are typically non-gradient. In a gradient lattice gas, the microscopic current is the discrete gradient of a local function, and the diffusion coefficient can often be extracted directly. In the Kob–Andersen model, the current contains nontrivial constraint factors, so it is not expressible as a simple lattice gradient of a local observable. This prevents a direct microscopic identification of the macroscopic diffusion coefficient and is the central reason KCHP transport is mathematically delicate (Arita et al., 2017).

4. Hydrodynamic limit and diffusion coefficients

Under coarse-graining, KCHP-like models exhibit hydrodynamic diffusion with density-dependent mobility. For the constrained lattice gases studied with boundary reservoirs, weakly driven large systems are described by the nonlinear diffusion equation

$\sigma_i=0$4

while in isotropic two-dimensional form one writes

$\sigma_i=0$5

The diffusion coefficient is model-dependent and encapsulates the effect of the local kinetic rule on macroscopic transport (Teomy et al., 2016, Arita et al., 2017).

A simple and widely used approximation is the no-correlation estimate

$\sigma_i=0$6

obtained by replacing local occupancies by the local density. For the two-dimensional Kob–Andersen model this gives

$\sigma_i=0$7

and for the spiral model

$\sigma_i=0$8

In the Kob–Andersen case, the same expression also arises as the smallest-subspace evaluation of the Varadhan–Spohn variational formula, where it is an upper bound rather than an exact result (Teomy et al., 2016, Arita et al., 2017).

For non-gradient models, the more systematic route is the Varadhan–Spohn representation

$\sigma_i=0$9

with the infimum taken over local test functions. Restricting the minimization to a finite-dimensional local subspace supported on a finite set $m=2$0 yields

$m=2$1

with $m=2$2, and if $m=2$3 then $m=2$4. This provides a controlled sequence of improving upper bounds on the diffusion coefficient (Arita et al., 2017).

The representative diffusion expressions established in the literature are summarized below.

Model	Diffusion coefficient	Status
Kob–Andersen, 2D $m=2$5	$m=2$6	no-correlation estimate; also smallest-subspace upper bound
Spiral model	$m=2$7	no-correlation estimate
Triangular chain	$m=2$8	exact
Triangular lattice	$m=2$9	exact

Numerical and analytical comparisons show that the no-correlation approximation is often good but not exact in general. The true density profiles in boundary-driven simulations lie below the no-correlation prediction, which means the true $$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$0 is smaller than $$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$1. The discrepancy does not vanish in the thermodynamic limit and is attributed to out-of-equilibrium spatial correlations that appear even under infinitesimal driving (Teomy et al., 2016). At high density, the Kob–Andersen upper bounds vanish algebraically, while simulations suggest stronger suppression and possibly non-analytic behavior, with

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$2

which is characteristic of transport slowing near jamming (Arita et al., 2017).

5. Exact hydrodynamics from classical-to-quantum mapping

A major recent development is the exact hydrodynamic analysis of triangular KCHPs using a classical-to-quantum correspondence. The stochastic evolution of the configuration distribution $$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$3 is written as

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$4

with $$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$5 interpreted as a quantum spin Hamiltonian. Under the mapping,

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$6

and the equilibrium distribution corresponds to the exact product-state ground state

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$7

Near-equilibrium hydrodynamics is then extracted from the low-energy sector of the spin Hamiltonian (McRoberts et al., 23 Jul 2025).

The analysis rotates the spins so that $$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$8 becomes the Holstein–Primakoff vacuum and expands the Hamiltonian as

$$P^{(1,0)}_{i,j}(\tau)=\tau_{i,j}(1-\tau_{i+1,j})H_{i,j}(\tau),$$9

The quadratic part gives a one-magnon dispersion $H_{i,j}(\tau)$0, from which the bare diffusion constant is read off. For the triangular chain,

$H_{i,j}(\tau)$1

so

$H_{i,j}(\tau)$2

For the triangular lattice,

$H_{i,j}(\tau)$3

so

$H_{i,j}(\tau)$4

The central result is that these bare spin-wave values are exact: the one-magnon self-energy begins at too high an order in momentum to renormalize the $H_{i,j}(\tau)$5 diffusion pole, so the non-interacting spin-wave theory predicts the exact diffusion constant for these models (McRoberts et al., 23 Jul 2025).

The same work conjectures that exactness extends to all KCHPs with hard-core occupancy, parity symmetry, and three-site gates. By contrast, four-site-gate models admit an additional $H_{i,j}(\tau)$6 magnon vertex at order $H_{i,j}(\tau)$7, which produces a self-energy correction of order $H_{i,j}(\tau)$8 and hence renormalizes the diffusion constant downward. Numerical results in the same study support this distinction: the L1R1 chain is consistent with the three-site-gate exactness conjecture, whereas the OR chain shows a visible reduction below the bare spin-wave prediction (McRoberts et al., 23 Jul 2025).

6. Related continuum formulations and broader mathematical context

The KCHP mechanism also has a continuum analogue. A continuum individual-based model of hopping and coalescing particles considers microscopic states $H_{i,j}(\tau)$9 in the configuration space

$n_\alpha(\vec r,t)\in\{0,1\},$0

with a generator

$n_\alpha(\vec r,t)\in\{0,1\},$1

where $n_\alpha(\vec r,t)\in\{0,1\},$2 describes coalescence of a pair $n_\alpha(\vec r,t)\in\{0,1\},$3 into $n_\alpha(\vec r,t)\in\{0,1\},$4, and $n_\alpha(\vec r,t)\in\{0,1\},$5 describes hopping of a particle from $n_\alpha(\vec r,t)\in\{0,1\},$6 to $n_\alpha(\vec r,t)\in\{0,1\},$7. The jump and coalescence rates depend on the entire surrounding configuration through multiplicative factors such as

$n_\alpha(\vec r,t)\in\{0,1\},$8

with $n_\alpha(\vec r,t)\in\{0,1\},$9. If $\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$0, these factors are $\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$1, so surrounding particles suppress the event rate. This is a continuum form of repulsive kinetic constraint or crowding-induced inhibition (Pilorz, 2015).

Using a Vlasov-type scaling, that continuum model passes from a hierarchy for correlation functions to a mesoscopic nonlinear kinetic equation for the density $\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$2. The resulting evolution contains gain and loss terms for both coalescence and jumps, each multiplied by exponential suppression factors of the form

$\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$3

and the corresponding initial-value problem has a unique local classical solution for bounded nonnegative initial data (Pilorz, 2015). This suggests that the KCHP idea is not confined to lattice gases: the essential structure is configuration-dependent inhibition of motion, with or without additional reaction channels.

From the broader mathematical viewpoint, kinetically constrained models are defined by update families $\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$4 and constraints

$\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$5

with generator

$\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$6

The general theory shows that ergodicity, mixing, and spectral properties are controlled by bootstrap-percolation-type closure mechanisms. In that setting, conservative kinetically constrained lattice gases and Kob–Andersen models are the natural nearest relatives of KCHP, while East and FA-$\rho(\vec r,t)=\langle n_\alpha(\vec r,t)\rangle.$7f models provide nonconservative analogues of the same facilitation principle (Hartarsky et al., 2024).

Taken together, these developments characterize KCHP as a transport theory of vacancy-mediated mobility. Its distinctive combination of simple equilibrium measure, constrained microscopic hopping, non-gradient current structure, and density-dependent hydrodynamics makes it a central model class for studying how local blocking and facilitation generate emergent diffusion, correlation corrections, and, in some cases, exact transport coefficients.