Fredrickson–Andersen j-Spin Model

• The Fredrickson–Andersen j-spin model is a kinetically constrained system where a spin flips only if a threshold number of neighboring sites are active, leading to nontrivial dynamics despite a trivial equilibrium measure.
• It can be formulated using either Ising spins or occupation variables, with its dynamics mapping onto bootstrap percolation and mode-coupling theories to yield precise critical exponents and scaling laws.
• The model spans diverse settings—from Bethe lattices and finite-dimensional systems to quantum spin chains—demonstrating how local facilitation influences glassy behavior and phase transitions.

The Fredrickson–Andersen $j$-spin facilitated model is a family of kinetically constrained spin models in which the equilibrium measure is trivial or non-interacting, while the dynamics is restricted by a facilitation rule: a site can change state only if at least $j$ neighboring sites are in a facilitating state. In the literature the same parameter is often denoted $f$, so the “$j$-spin” and “$f$-spin” formulations are equivalent at the level of the constraint $j\leftrightarrow f$. Depending on convention, the microscopic variables are either Ising spins $\sigma_i\in\{-1,+1\}$ or occupation variables $\eta_x\in\{0,1\}$; in both representations, glassy behavior is generated by dynamic facilitation rather than by nontrivial static interactions (Sellitto et al., 2010).

1. Microscopic definition and notational conventions

In the occupancy formulation on $\mathbb Z^d$, the configuration space is $\Omega=\{0,1\}^{\mathbb Z^d}$, with $j$0 interpreted as an empty, infected, or excited site and $j$1 as a filled, healthy, or inactive site. A standard FA-$j$2f dynamics is a continuous-time Markov process in which each site attempts to resample from a Bernoulli$j$3 law when its Poisson clock rings, but the update is accepted only if at least $j$4 nearest neighbors are empty. The corresponding constraint indicator is

$j$5

The Dirichlet form can be written as

$j$6

and the generator as

$j$7

in the two-dimensional presentation (Hartarsky et al., 2022).

In the Ising-spin formulation often used on Bethe lattices and random regular graphs, one takes $j$8 with a non-interacting Hamiltonian

$j$9

so the equilibrium up-spin probability is

$f$0

The single-spin flip rate is Glauber- or Metropolis-like,

$f$1

but only if the facilitation constraint is satisfied. In the FA model on a Bethe lattice studied with random pinning, the rule is that the number of nearest neighbors in state $f$2 must be greater than or equal to $f$3, with the concrete choice $f$4 and $f$5 corresponding to a two-spin facilitated model in the usual FA classification (Ikeda et al., 2015).

These two descriptions are equivalent up to convention. A common source of confusion is the identity of the facilitating state: some papers use empty sites or zeros, others down spins or defects. The invariant content is the same. The Hamiltonian is deliberately simple, and the slow dynamics arises because mobility can propagate only through neighborhoods that already contain enough facilitating sites. This separation between trivial thermodynamics and nontrivial kinetics is one of the defining features of the FA family (Hartarsky et al., 2022).

2. Bethe-lattice formulation and bootstrap-percolation structure

On locally tree-like graphs, especially random regular graphs and Bethe lattices, the FA $f$6-spin model admits an exact description in terms of bootstrap percolation or $f$7-core–type recursion relations. This is the canonical mean-field setting for the model. On a 4-regular random graph, for example, the glass order parameter can be taken as the fraction $f$8 of frozen spins or, dynamically, as the long-time limit of the persistence $f$9, where $j$0 is the probability that a spin has never flipped up to time $j$1 (Fennell et al., 2014).

The bootstrap-percolation interpretation is central. For uniform facilitation on a Bethe lattice, dynamical arrest corresponds to the appearance of an infinite cluster of frozen spins. For $j$2, the transition is of bootstrap-percolation type and is hybrid: the frozen fraction jumps discontinuously, but critical fluctuations diverge. For $j$3 or $j$4, the transition is continuous and equivalent to standard percolation. In facilitated mixtures on a $j$5 Bethe lattice, the distribution

$j$6

interpolates between these regimes. The resulting frozen fraction $j$7 is obtained exactly from a recursive solution involving rooted-tree probabilities $j$8 and $j$9, and the critical line is

$f$0

For $f$1, the transition is discontinuous with $f$2; for $f$3, it is continuous with $f$4 when $f$5 and $f$6 when $f$7 (Sellitto et al., 2010).

A closely related dynamical viewpoint is obtained from exact equations for occupation correlations on the Bethe lattice. In the oriented FA model on a $f$8-ary rooted tree with facilitation $f$9, the long-time blocked fraction $j\leftrightarrow f$0 satisfies

$j\leftrightarrow f$1

This is simultaneously the bootstrap-percolation fixed-point equation and the equation that determines the plateau of the dynamical occupation correlation. The same structure governs the asymptotic memory kernel of the model, which makes the bootstrap-percolation mapping not merely geometric but dynamically exact in the scaling regime near the transition (Önder et al., 2019).

The Bethe-lattice formulation therefore supplies more than a mean-field caricature. It identifies the precise blocked structures, yields exact recursion relations for the nonergodicity parameter, and separates two distinct geometries of frozen clusters: compact bootstrap-percolation clusters in the discontinuous regime and fractal standard-percolation clusters in the continuous regime. This distinction reappears in the critical exponents and in the character of relaxation.

3. Mean-field criticality, memory kernels, and higher-order singularities

Near the Bethe-lattice dynamical transition, the FA $j\leftrightarrow f$2-spin model displays the full mode-coupling-theory scaling scenario. For the occupation correlation

$j\leftrightarrow f$3

the Mori–Zwanzig equation takes the form

$j\leftrightarrow f$4

In the scaling regime near the critical point, the memory kernel becomes asymptotically exact, and the deviation $j\leftrightarrow f$5 from the plateau satisfies the standard MCT scaling equation. The exponents $j\leftrightarrow f$6 and $j\leftrightarrow f$7 obey

$j\leftrightarrow f$8

and the characteristic times diverge as

$j\leftrightarrow f$9

For the case $\sigma_i\in\{-1,+1\}$0, $\sigma_i\in\{-1,+1\}$1, the exponent parameter is $\sigma_i\in\{-1,+1\}$2, with $\sigma_i\in\{-1,+1\}$3, $\sigma_i\in\{-1,+1\}$4, and $\sigma_i\in\{-1,+1\}$5 (Önder et al., 2019).

Random pinning reveals a higher-order version of the same structure. On a 4-regular Bethe lattice with facilitation $\sigma_i\in\{-1,+1\}$6, pinning a fraction $\sigma_i\in\{-1,+1\}$7 of spins from an equilibrium configuration generates a line of dynamical transitions $\sigma_i\in\{-1,+1\}$8 in the $\sigma_i\in\{-1,+1\}$9–$\eta_x\in\{0,1\}$0 plane. For $\eta_x\in\{0,1\}$1, the transition is discontinuous in the frozen fraction and obeys the square-root law

$\eta_x\in\{0,1\}$2

the hallmark of an $\eta_x\in\{0,1\}$3 singularity. The line terminates at

$\eta_x\in\{0,1\}$4

where the order parameter becomes continuous but singular,

$\eta_x\in\{0,1\}$5

and the relaxation becomes logarithmic,

$\eta_x\in\{0,1\}$6

The susceptibility

$\eta_x\in\{0,1\}$7

shows the corresponding $\eta_x\in\{0,1\}$8 and $\eta_x\in\{0,1\}$9 divergences, with $\mathbb Z^d$0 on the $\mathbb Z^d$1 line and $\mathbb Z^d$2 at the terminal point (Ikeda et al., 2016).

The terminal point also admits a geometrical interpretation in terms of avalanches and self-induced disorder. In the pinned Bethe-lattice model, changing $\mathbb Z^d$3 or unblocking a single site can trigger cascades of unblocking events with size distribution

$\mathbb Z^d$4

identical in scaling form to the Random Field Ising Model at its continuous hysteresis transition. This gives a microscopic explanation for the $\mathbb Z^d$5 order-parameter exponent and the $\mathbb Z^d$6 divergence of fluctuations (Ikeda et al., 2016).

A related Bethe-lattice analysis of random pinning showed that the same terminal point has correlation-length exponents matching the inhomogeneous mode-coupling-theory prediction for an $\mathbb Z^d$7 singularity. In the $\mathbb Z^d$8, $\mathbb Z^d$9 model, the terminal point is at

$\Omega=\{0,1\}^{\mathbb Z^d}$0

and the point-to-set correlation length scales as $\Omega=\{0,1\}^{\mathbb Z^d}$1 along the $\Omega=\{0,1\}^{\mathbb Z^d}$2 line, $\Omega=\{0,1\}^{\mathbb Z^d}$3 in the generic approach to the $\Omega=\{0,1\}^{\mathbb Z^d}$4 point, and $\Omega=\{0,1\}^{\mathbb Z^d}$5 along the tangent direction to the transition line (Ikeda et al., 2015).

Taken together, these results show that the mean-field FA $\Omega=\{0,1\}^{\mathbb Z^d}$6-spin model is not merely qualitatively glassy. It realizes the full $\Omega=\{0,1\}^{\mathbb Z^d}$7 and $\Omega=\{0,1\}^{\mathbb Z^d}$8 phenomenology, with exact exponent relations, logarithmic higher-order relaxation, and RFIM-type avalanche criticality.

4. Finite-dimensional behavior and rigorous time-scale asymptotics

In finite dimensions, the behavior depends sharply on the facilitation threshold. In two dimensions, the FA-1f model is non-cooperative, the FA-2f model is cooperative, and FA-3f and FA-4f have $\Omega=\{0,1\}^{\mathbb Z^d}$9 and $j$00 with positive probability for any $j$01. The two-dimensional review emphasizes that FA-1f and FA-2f are therefore the nontrivial ergodic cases, with radically different low-density asymptotics (Hartarsky et al., 2022).

For FA-1f in stationarity, the spectral gap has dimension-dependent power-law behavior. In $j$02,

$j$03

in $j$04,

$j$05

and in $j$06,

$j$07

The persistence function $j$08 has exponential tails with the corresponding characteristic scales: $j$09 while $j$10 satisfies the matching lower bounds $j$11, $j$12, and $j$13 in dimensions $j$14, $j$15, and $j$16, respectively (Shapira, 2020).

For FA-2f, the asymptotics are much more singular. In $j$17, the infection time of the origin in the stationary process obeys

$j$18

and the same leading asymptotics hold for $j$19. In $j$20, since $j$21, this becomes

$j$22

with corresponding high-probability bounds for $j$23. This was the first sharp threshold result for a critical KCM and explicitly tied the FA-2f timescale to the sharp threshold constant for 2-neighbor bootstrap percolation (Hartarsky et al., 2020).

The two-dimensional overview reformulates this result as

$j$24

for stationary FA-2f on $j$25, and interprets the mechanism in terms of rare mobile droplets whose probability is $j$26. This settled several long-standing controversies in the physics literature over the correct FA-2f scaling regime (Hartarsky et al., 2022).

A persistent misconception is that all FA models should share a single low-temperature universality class. The rigorous results show otherwise. FA-1f has polynomial scales and mobile single-defect dynamics; FA-2f has stretched-exponential scales governed by critical droplets and bootstrap-percolation constants; higher facilitation in low dimension can produce genuine frozen configurations with positive probability. The facilitation threshold is therefore not a perturbative parameter but a structural classifier of the dynamics.

5. Front propagation, transport, and trajectory-space phase behavior

A particularly detailed picture is available for the nonequilibrium front in one-dimensional FA-1f. Starting from a configuration that is fully occupied on the left half-line and has a single zero at the origin, one tracks the leftmost zero $j$27. For $j$28, where

$j$29

the process seen from the front has a unique invariant measure $j$30, and there exist $j$31 and $j$32 such that

$j$33

with

$j$34

Thus the front moves ballistically and has Gaussian fluctuations at scale $j$35 (Blondel et al., 2018).

This front motion controls finite-volume relaxation. On an interval $j$36 with empty boundary conditions at 0 and $j$37, FA-1f exhibits cutoff at time

$j$38

with window $j$39 for a positive constant $j$40. More precisely, for $j$41, the worst-case total variation distance satisfies

$j$42

for suitable $j$43 and all large $j$44 (Ertul, 2021).

Transport by a tagged particle in FA environments has a similarly sharp characterization. For noncooperative KCM in the $j$45-zeros class, the tracer diffusion coefficient satisfies

$j$46

For FA-1f, which corresponds to $j$47 in this class, this gives

$j$48

This rigorously confirmed the $j$49 prediction made earlier for the one-spin facilitated FA model (Blondel, 2013).

The FA family also supports trajectory-space phase transitions. In the one-dimensional FA model, the large deviations of the dynamical activity $j$50 define a scaled cumulant generating function $j$51 with a first-order singularity at $j$52, corresponding to coexistence between active and inactive dynamical phases. In a two-dimensional spin-model representation of trajectories, the susceptibility at coexistence diverges as

$j$53

in the long-time dynamical limit $j$54, while for isotropic aspect ratio $j$55 it follows standard first-order finite-size scaling $j$56 (Jack et al., 2019).

These results indicate that FA dynamics supports several distinct macroscopic structures: moving fronts in real space, tracer subdynamics in a fluctuating constrained environment, and active–inactive coexistence in trajectory space. They are different manifestations of the same facilitation principle.

6. Random environments, pinned variants, and broader connections

Quenched heterogeneity modifies FA $j$57-spin dynamics in a way that sharply separates local and global timescales. In the mixed-threshold FA model on $j$58, each site is assigned a random threshold $j$59 with $j$60 and $j$61, so the system is a quenched mixture of FA-1f and FA-2f constraints. For the associated bootstrap percolation, the emptying time of the origin has deterministic scaling

$j$62

in two dimensions, and more generally $j$63 in $j$64. For the stochastic KCM, by contrast, the local emptying time has random exponents: $j$65 with $j$66 depending on the disorder realization. At the same time, the global relaxation time is dominated by rare FA-2f-like islands and satisfies

$j$67

almost surely (Shapira, 2018).

A different cross-disciplinary extension arises in noisy quantum many-body dynamics. In the strong-noise limit of an exchange-coupled qubit chain without conserved spin component, the second moments of operators map onto an FA-1f process with binary variables $j$68. The effective local rates are

$j$69

and the stationary excitation probability is $j$70. In this representation, operator spreading becomes front propagation in FA-1f. In one dimension, the asymptotic front fluctuations are Gaussian, and numerical evidence in two dimensions supports the conjecture that the long-time fluctuations belong to the Kardar–Parisi–Zhang universality class (Rowlands et al., 2018).

Random pinning, mixed facilitation, and quantum mappings all reinforce the same conclusion: the FA $j$71-spin model is not a narrowly defined toy process but a flexible dynamical framework. It supports mean-field glass transitions with $j$72 and $j$73 structure, sharp finite-dimensional threshold laws, trajectory-space first-order transitions, and mappings to problems outside classical glass theory. A plausible implication is that dynamic facilitation organizes a broader universality structure than any single microscopic realization, provided the local rule still enforces mobility through neighboring mobile sites.