Prethermal Kinetically Constrained Dynamics

• Prethermal kinetically constrained dynamics are defined by long-lived, nearly stationary regimes emerging when local move restrictions inhibit rapid equilibration.
• This phenomenon is observed in open quantum systems and Rydberg gases where kinetic constraints yield metastable plateaus and hierarchical relaxation.
• Effective modeling using large-deviation theory and hydrodynamic descriptions provides quantitative insights into fragmentation and delayed thermalization.

Searching arXiv for recent and foundational papers on prethermal kinetically constrained dynamics, kinetic constraints, and related hydrodynamic/trajectory-space formulations. arXiv search: "prethermal kinetically constrained dynamics open quantum systems Rydberg facilitation" Prethermal kinetically constrained dynamics denotes long-lived, nearly-stationary regimes in which a many-body system evolves under strong dynamical constraints before eventually relaxing to its true steady state. In this class of phenomena, the dominant obstruction to equilibration is not a complicated equilibrium free-energy landscape but the restriction of allowed local moves, whether imposed directly, generated by facilitation rules, or emerging from strong interactions, noise, or conservation laws. Across open quantum spin systems, clean constrained quantum lattices, driven and stochastic kinetically constrained models, and experimentally realized Rydberg gases, the recurring structure is a simple or even trivial stationary measure together with highly nontrivial transient evolution featuring metastability, hierarchical relaxation, dynamical heterogeneity, and, in some cases, fragmentation or freezing (Everest et al., 2016, Lan et al., 2017, Valado et al., 2015).

1. Conceptual structure and defining features

Kinetically constrained dynamics is characterized by transition rules in which a local move is allowed only if a specified local condition is satisfied. In the canonical glass-theory setting, the static equilibrium measure is simple, yet the dynamics is highly nontrivial and slow. This separation between trivial thermodynamics and complex dynamics is central to both classical kinetically constrained models and their quantum or open-system analogues (Everest et al., 2016, Speck et al., 2010).

The “prethermal” qualifier refers to an intermediate regime in which observables relax within a restricted dynamical sector and appear stationary, or only very slowly evolving, long before the full system reaches its final steady state. In open constrained Rydberg systems this restricted sector can be produced by strong dephasing and resonance selection; in clean constrained quantum models it can arise from local facilitation rules, flippability constraints, or effective fragmentation of Hilbert space; in classical KCMs it appears as a long-lived inactive or low-activity regime in trajectory space (Valado et al., 2015, Lan et al., 2017, Bodineau et al., 2011).

A basic distinction is between hard and soft constraints. Hard constraints forbid certain moves exactly and can produce dynamical reducibility or true fragmentation of configuration space. Soft constraints suppress moves by large but finite factors, making the constrained regime effectively stable on experimentally or numerically accessible times while permitting eventual escape on much longer scales. This distinction is explicit in open quantum realizations where transition rates depend on energy mismatches, and it underlies the difference between transient prethermal plateaus and true frozen phases (Everest et al., 2016, Morningstar et al., 2020).

Another defining feature is the coexistence of active and inactive dynamical sectors. In trajectory-space formulations, activity or entropy production plays the role of an order parameter, and the relevant transition is a space-time transition rather than a thermodynamic one. This makes prethermal kinetically constrained dynamics fundamentally dynamical: what changes is the structure of trajectories, not the equilibrium ensemble of configurations (Speck et al., 2010, Speck et al., 2011).

2. Microscopic mechanisms that generate constraints

One general route to constrained prethermality is the emergence of effective classical jump dynamics in open quantum systems under strong dephasing. In the mechanism developed for quantum spin systems with Hamiltonian $H=H_C+H_Q$, rapid local dephasing $\gamma \gg \Omega$ suppresses coherences and projects the dynamics onto classical configurations $|\mathcal C_m\rangle$. The effective transition rates depend on the energy change $\delta E$ induced by $H_C$:

$$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$

Resonant moves with $\delta E=0$ occur at maximal rate, while off-resonant moves with $|\delta E| \gg \gamma$ are strongly suppressed. By engineering $H_C$, one therefore engineers the kinetic constraint itself (Everest et al., 2016).

A concrete experimental realization is the dissipative Rydberg gas in the strong-Zeno regime. There the effective local flip rate is

$$\Gamma_i(\Delta) = \frac{\Omega^2}{2\gamma}\; \frac{1}{1 + \left( \frac{\Delta - \frac{1}{\hbar}\sum_{j\neq i}V_{ij}n_j}{\gamma} \right)^2},$$

so the environment-dependent interaction shift $\gamma \gg \Omega$0 determines whether a site is blocked or facilitated. At $\gamma \gg \Omega$1 one obtains blockade constraints with blockade radius

$\gamma \gg \Omega$2

while for $\gamma \gg \Omega$3 one obtains facilitation at radius

$\gamma \gg \Omega$4

with resonant-shell width

$\gamma \gg \Omega$5

This yields controllable anti-correlated or correlated dynamics in a many-body system with an otherwise trivial long-time state (Valado et al., 2015).

In isolated quantum systems, constraints can be built directly into translationally invariant Hamiltonians. The constrained quantum lattice gas has bond operators $\gamma \gg \Omega$6 that allow hopping only when at least one common neighbor is empty, while the quantum dimer model permits dynamics only on flippable plaquettes. Around the Rokhsar–Kivelson points, these systems interpolate between fast, homogeneous thermalization and slow, metastable relaxation; on the slow side, most local moves are blocked and relaxation proceeds through rare cooperative processes (Lan et al., 2017).

A closely related constrained family arises in dual XXZ models. In the XPX model,

$\gamma \gg \Omega$7

the kinetic term acts only when the central spin is down. This produces fractonic-like mobility rules in which isolated $\gamma \gg \Omega$8 spins are nearly immobile while pairs of adjacent $\gamma \gg \Omega$9 spins can move, thereby generating slow and heterogeneous finite-energy relaxation (Zadnik et al., 2023).

Constraints can also coexist with additional slow degrees of freedom. In facilitated Rydberg chains coupled to trap vibrations, the constrained manifold of single and double neighboring excitations is dressed by phonons, producing heavy polaronic quasiparticles. The effective constrained Hamiltonian contains a two-level internal structure and a selective spin–phonon coupling that acts only for the two-excitation configuration, leading to band flattening and increased effective mass (Mazza et al., 2020).

3. Dynamical manifestations: metastability, plateaus, and heterogeneity

The hallmark phenomenology is a pronounced separation of timescales. In the open quantum reaction-diffusion model with bond-conserving resonant moves, excitations self-organize into mobile monomers, slow polymers, and immobile plaquettes. The plaquette density

$|\mathcal C_m\rangle$0

can overshoot its stationary value, while the imbalance

$|\mathcal C_m\rangle$1

shows a rapid initial decay followed by a much slower regime. The fast timescale is quoted as roughly

$|\mathcal C_m\rangle$2

after which assisted diffusion and rare cooperative rearrangements dominate. This is a prototypical prethermal plateau: local observables look nearly relaxed within a constrained sector, but the true infinite-temperature steady state is reached only after much longer evolution (Everest et al., 2016).

In clean constrained quantum systems, the same structure appears in two-step relaxation. For the constrained quantum lattice gas, the density autocorrelator

$|\mathcal C_m\rangle$3

rapidly decays from $|\mathcal C_m\rangle$4 to a nonzero plateau and only later relaxes to $|\mathcal C_m\rangle$5. The plateau is controlled by mobile vacancy dimers that reorganize the state on a timescale $|\mathcal C_m\rangle$6, followed by a metastable stage in which isolated vacancies block further motion and relaxation requires rarer many-body processes. In the quantum dimer model, large $|\mathcal C_m\rangle$7 similarly produces a plateau in the time-averaged dimer correlator and a delayed growth of the correlation length $|\mathcal C_m\rangle$8, indicating that extended correlations develop only after the system begins to escape the inactive manifold (Lan et al., 2017).

Dynamical heterogeneity is not merely a classical analogue but an explicit quantum signature. In the slow regime of the constrained quantum lattice gas, entanglement growth is spatially inhomogeneous: regions containing mobile vacancy dimers entangle rapidly, while regions containing isolated vacancies remain almost frozen for long times. The time-averaged state is well approximated by an effectively thermalized active region tensored with an almost unchanged inactive region during the metastable window. In the dimer model, active flippable regions spread while inert regions persist up to times of order $|\mathcal C_m\rangle$9 in the numerical data, and the correlation length grows only when the plateau in the two-time correlator begins to decay (Lan et al., 2017).

An atomistic glass-forming liquid shows the same trajectory-level structure. Defining a long-lived displacement indicator

$\delta E$0

the total excitation count

$\delta E$1

has a distribution with exponential low-mobility tails for observation times longer than the structural relaxation time. Trajectories biased toward low $\delta E$2 enter an inactive phase that remains at high overlap with the initial state for many structural relaxation times, while structural observables differ only weakly from those of the active phase. This demonstrates that long-lived constrained regimes can persist even in fully atomistic models with no nontrivial static phase structure (Speck et al., 2011).

4. Trajectory-space, hydrodynamic, and scaling descriptions

A unifying formalism is the large-deviation theory of trajectories. For activity $\delta E$3, the $\delta E$4-ensemble is defined by

$\delta E$5

with mean activity density $\delta E$6. In kinetically constrained models, $\delta E$7 exhibits a first-order singularity at $\delta E$8: $\delta E$9 selects active trajectories with finite activity density, whereas $H_C$0 selects inactive trajectories with vanishing activity density. This establishes active/inactive coexistence as a space-time phase transition and provides a direct dynamical underpinning for metastability and prethermal plateaus (Speck et al., 2010, Bodineau et al., 2011).

For driven KCMs, entropy production $H_C$1 plays an analogous role. The tilted generator

$H_C$2

yields a large-deviation function $H_C$3 for entropy production. In both driven FA and driven triangular lattice-gas models, $H_C$4 displays a first-order transition between finite- and zero-entropy-production trajectories, again showing that constrained systems naturally support long-lived low-activity regimes that are typical in biased ensembles and metastable in the unbiased dynamics (Speck et al., 2010).

The same space-time viewpoint survives in explicit activity-transition analyses of FA-1f and East models. With total activity $H_C$5, the scaled cumulant generating function

$H_C$6

exhibits a first-order dynamical transition. In finite-size scaling $H_C$7, the limiting function is linear, $H_C$8, on the active side and saturates at $H_C$9 on the inactive side, where $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$0 is an interface energy between active and inactive dynamical phases. This formalizes the blocking mechanism by which reducing activity is much cheaper than increasing it and explains why low-activity prethermal regimes are statistically prominent in constrained dynamics (Bodineau et al., 2011).

At coarse-grained scales, many constrained systems admit nonlinear hydrodynamic descriptions. Boundary-driven kinetically constrained lattice gases obey

$$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$1

with model-dependent nonlinear diffusion coefficient $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$2. For the Kob–Andersen model under the no-correlation approximation,

$$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$3

while for the spiral model

$$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$4

The paper shows that even infinitesimal driving generates non-negligible correlations that renormalize these diffusion coefficients, making hydrodynamic transport itself an emergent description of slow constrained relaxation rather than a trivial consequence of independent-site statistics (Teomy et al., 2016).

A recent rigorous extension constructs one-dimensional symmetric exclusion processes satisfying the gradient condition with diffusivities expressed in a Bernstein polynomial basis, yielding hydrodynamic limits of the form

$$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$5

For the continuously parametrised porous media construction, the macroscopic equation is a porous media equation with non-integer exponent $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$6, providing a controllable family of non-linear diffusion laws arising from microscopic kinetic constraints (Nahum, 16 Apr 2025).

A different scaling perspective comes from perturbative dynamics of the FA model on a Bethe lattice. There, a hierarchy of dynamical systems indexed by perturbation order $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$7 increasingly resolves local constrained structure and yields relaxation times $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$8 that become systematically closer to Monte Carlo results as $$\frac{1}{\Gamma_\nu} = 1 + \left(\frac{2\,\delta E}{s\,\gamma}\right)^2 .$$9 increases. The resulting scaling law relates time and size scales near the nonergodic transition through a minimal rearrangement scale $\delta E=0$0 and a dynamical exponent $\delta E=0$1, suggesting a universal relation $\delta E=0$2 within fixed facilitation class (Ohta, 2010).

5. Fragmentation, freezing, and localization-enhanced constrained dynamics

Kinetic constraints may produce merely delayed ergodicity, but they can also yield fragmentation or genuine freezing. In a one-dimensional dipole-conserving stochastic lattice gas with $\delta E=0$3 and four-site local updates, total charge

$\delta E=0$4

and total dipole moment

$\delta E=0$5

are conserved. As the mean density $\delta E=0$6 is tuned, the system passes from a weakly fragmented, subdiffusive thermalizing phase near half filling to a strongly fragmented frozen phase above

$\delta E=0$7

with strong evidence for $\delta E=0$8. In the thermalizing phase,

$\delta E=0$9

implying $|\delta E| \gg \gamma$0 and $|\delta E| \gg \gamma$1. In the frozen phase, the fraction of permanently frozen sites scales as

$|\delta E| \gg \gamma$2

and active bubbles cease to exchange particles with the background. This is not merely prethermal but truly non-ergodic; however, the weakly fragmented phase and the critical regime display prethermal-like slow relaxation with very large dynamical exponents (Morningstar et al., 2020).

In clean quantum constrained models, fragmentation is softer but still operationally significant. The dual XXZ family exhibits localized and delocalized ground-state phases, and on the localized side finite-energy relaxation is slow, autocorrelations display plateaus, and the relaxation time scales as

$|\delta E| \gg \gamma$3

This produces a genuine constrained prethermal regime even in nonintegrable deformations: the system first relaxes within kinetically restricted sectors and only later crosses over to full thermal behavior (Zadnik et al., 2023).

Disorder can stabilize these long-lived constrained sectors into many-body localized behavior. In the constrained hard-core boson model on a triangular ladder, density autocorrelations

$|\delta E| \gg \gamma$4

show plateaus already in the clean, strongly constrained regime, especially for initial states containing one vacancy dimer and one isolated hole. Adding uncorrelated disorder quickly produces non-decaying autocorrelations and localized eigenstates at disorder strengths much smaller than in the unconstrained model. A plausible implication is that kinetic constraints create a prethermal landscape of heavy, weakly hybridizing excitations, and weak disorder converts this landscape into a many-body localized phase more efficiently than in unconstrained systems (Royen et al., 2023).

Constraint-induced inactive phases can also arise in deterministic classical dynamics. In a periodically driven classical Heisenberg spin system with facilitation angle $|\delta E| \gg \gamma$5, the activity density

$|\delta E| \gg \gamma$6

undergoes a continuous transition from an active chaotic phase to a frozen absorbing phase. Extensive numerical scaling shows that the transition is in the directed percolation universality class in both one and two spatial dimensions, with critical exponents matching directed percolation. This indicates that near the critical constraint strength, long-lived quasi-inactive regimes are organized by absorbing-state criticality (Deger et al., 2022).

6. Experimental realizations, limitations, and open directions

The most direct experimental realization is the cold Rydberg gas with tunable blockade and facilitation. In the blockade regime, the excitation growth rate per atom collapses as the mean spacing $|\delta E| \gg \gamma$7 approaches the blockade radius $|\delta E| \gg \gamma$8, dropping by about four orders of magnitude between $|\delta E| \gg \gamma$9 and $H_C$0. In the facilitation regime at $H_C$1, the dynamics separates into a nucleation regime ($H_C$2), a facilitation regime ($H_C$3), and a saturation regime ($H_C$4), while the Mandel parameter

$H_C$5

becomes strongly positive during facilitation. Because the $H_C$6 Rydberg lifetime is $H_C$7, the experimentally accessible window is naturally an intermediate constrained regime rather than the true asymptotic steady state of the full open system (Valado et al., 2015).

Cold atoms in optical lattices, programmable Rydberg arrays, and superconducting qubit platforms are repeatedly identified as natural settings for realizing constrained hopping, density-dependent tunneling, dimer dynamics, or blockade/facilitation rules. In the vibrationally dressed Rydberg setting, the exaggerated length scales permit site-resolved monitoring of both spin and phonon degrees of freedom, making it possible in principle to track the real-space phonon cloud of a heavy constrained polaron and thereby directly observe a slow constrained quasiparticle rather than infer it indirectly (Lan et al., 2017, Mazza et al., 2020, Royen et al., 2023).

A persistent limitation is the distinction between transient glassiness and asymptotic non-ergodicity. In open systems with pure dephasing and coherent drive, the effective constrained dynamics often equilibrates to an infinite-temperature uniform distribution on each connected component, so any glassy behavior is transient unless extra dissipative channels stabilize low-density sectors. In clean isolated models such as the constrained quantum lattice gas and quantum dimer model, the evidence supports eventual thermalization consistent with ETH, even though the metastable regime can span many decades in time. By contrast, in fracton-like dipole-conserving stochastic models, the frozen phase is truly non-ergodic. It is therefore essential to distinguish delayed thermalization from exact fragmentation or absorbing-state arrest (Everest et al., 2016, Lan et al., 2017, Morningstar et al., 2020).

Open problems remain sharp even for the simplest classical one-facilitated models. Recent work on FA-1f and related variants emphasizes that in the interesting regime where the equilibrium density of facilitating vertices is small, many fundamental questions concerning the non-stationary evolution remain unsolved. Conjecturally, FA-1f on $H_C$8 should converge to equilibrium from any initial law containing at least one infection almost surely, and under mild exponential control on large healthy intervals the convergence should be exponentially fast; however, rigorous results at low infection density are still incomplete (Martinelli et al., 23 Oct 2025).

Taken together, the available results support a unified interpretation: prethermal kinetically constrained dynamics is the generic long-time but pre-asymptotic manifestation of facilitation, blocking, resonance selection, fragmentation, or conserved-mobility sectors in many-body evolution. Its observable content is a hierarchy of relaxation scales, plateaus in correlation functions, space-time coexistence of active and inactive regions, and transport laws renormalized by the constrained structure of microscopic trajectories (Speck et al., 2010, Nahum, 16 Apr 2025).