Fredkin Spin Chain

• The Fredkin spin chain is a one-dimensional spin-½ quantum lattice model defined via three-site interactions and Dyck-path combinatorics, leading to an exactly solvable, frustration-free ground state.
• Colored and area-deformed variants introduce rich entanglement behaviors—from logarithmic to volume-law scaling—and feature exact tensor-network representations that reveal unconventional quantum dynamics.
• Studies highlight anomalous spectral-gap scaling, nonergodic transport properties, and Hilbert space fragmentation with nonlocal symmetries, distinguishing it from conventional integrable models.

The Fredkin spin chain is a one-dimensional spin-$\tfrac12$ quantum lattice model with genuine three-site interactions, most commonly defined with open boundaries that pin the first spin up and the last spin down, and whose frustration-free ground state is exactly expressible in terms of Dyck-path combinatorics. In its original form the Hamiltonian can be written through controlled-SWAP (Fredkin) gates; in colored and area-deformed variants, the same combinatorial structure yields logarithmic, square-root, or volume-law entanglement, polynomially or faster closing gaps, anomalous transport, and a variety of nonergodic phenomena (Salberger et al., 2016, Zhang et al., 2017).

1. Hamiltonian and local constraint structure

A standard open-chain formulation on $N$ sites uses the Fredkin gate

$$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$

with boundary projectors fixing the left edge to $\uparrow$ and the right edge to $\downarrow$. One convenient expression is

$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$

$$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$

where $P_j^\pm=(1\pm \sigma_j^z)/2$ (Salberger et al., 2016). Equivalent projector formulations use singlet projectors on nearest-neighbor bonds conditioned by an adjacent spin, for example

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$

with $U_i=\tfrac12(1+\sigma_i^z)$, $N$0, and $N$1 (Adhikari et al., 2018).

The local moves implemented by these projectors are the Fredkin moves: adjacent $N$2 pairs can “bubble” past a neighboring spin only when the relevant control spin is in the correct state. This local kinetic rule is the quantum counterpart of a constrained nearest-neighbor exchange process. The Hamiltonian is frustration-free: each local projector annihilates the exact ground state (Dell'Anna et al., 2019).

The colored generalization enlarges the on-site basis to

$N$3

or equivalently $N$4, and adds nearest-neighbor color-matching terms that force each matched up-step and down-step to carry the same color (Zhang et al., 2017). In the original SU$N$5 formulation, the higher-spin half-integer extensions are realized as properly colored Dyck-path models with $N$6 colors (Salberger et al., 2016).

2. Dyck-path ground states and exact representations

The defining combinatorial map identifies $N$7 with an up-step and $N$8 with a down-step. A basis state is therefore a lattice walk, and the open-boundary ground space selects Dyck paths: walks from $N$9 to $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$0 that never cross below the horizontal axis. For the uncolored chain of length $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$1, the unique ground state is the equal-weight Catalan superposition

$$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$2

while the colored version is

$$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$3

(Salberger et al., 2016, Movassagh, 2016).

Area-deformed Fredkin chains reweight the Dyck configurations by the area $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$4 under the path. In the uniform deformation, one exact ground-state convention is

$$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$5

while another uses $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$6; both are explicit weighted Dyck-path superpositions differing by parameterization conventions (Zhang et al., 2017, Alexander et al., 2018). More generally, the multiparameter deformation assigns position-dependent weights while preserving frustration-freeness through the constraint $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$7 (Zhang et al., 2017).

These exact states admit explicit tensor-network descriptions. The original uncolored model has an exact matrix-product-state representation built from shift matrices $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$8 and $$F_{i,j,k} =|\downarrow\rangle\langle\downarrow|_i\otimes I_{j,k} +|\uparrow\rangle\langle\uparrow|_i\otimes \mathrm{SWAP}_{j,k},$$9 satisfying

$\uparrow$0

with an exact truncation of the bond dimension to $\uparrow$1 for a chain of length $\uparrow$2 (Salberger et al., 2016). For the area-deformed “rainbow” states, an exact two-dimensional bulk tensor network on an inverted step pyramid of identical four-leg tensors $\uparrow$3 reproduces the weighted Dyck-path amplitudes, with each horizontal bulk arrow contributing a factor $\uparrow$4; the construction uses only $\uparrow$5 copies of a single tensor, but the resulting network is not compressible into an MPS of finite bond dimension (Alexander et al., 2018).

3. Entanglement structure, correlations, and phase organization

The Fredkin chain is notable for entanglement scaling that goes well beyond the standard one-dimensional area-law paradigm. In the undeformed uncolored model, the half-chain entanglement entropy grows logarithmically: $\uparrow$6 or more explicitly

$\uparrow$7

(Salberger et al., 2016, Sugino et al., 2018). In the undeformed colored model, the von Neumann entropy scales as

$\uparrow$8

giving a square-root violation of the area law (Sugino et al., 2018).

For the Rényi entropies of the colored undeformed chain, the asymptotics are non-analytic at $\uparrow$9: for $\downarrow$0 they obey a volume law, at $\downarrow$1 they scale as $\downarrow$2, and for $\downarrow$3 they scale as $\downarrow$4 (Sugino et al., 2018). This establishes an entanglement transition in Rényi index rather than in a microscopic coupling.

In the uniform area-deformed chain, the phase structure is controlled by the deformation parameter $\downarrow$5 and the number of colors $\downarrow$6 (Zhang et al., 2017).

Regime	Half-chain entanglement	Gap characterization
$\downarrow$7, any $\downarrow$8	$\downarrow$9	gapped
$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$0, $$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$1	$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$2	gapless
$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$3, $$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$4	$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$5	gapless
$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$6, $$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$7	bounded $$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$8	gapless
$$H = H_{\partial}+\sum_{j=1}^{N-2}H_j,\qquad H_{\partial}=P_1^-+P_N^+,$$9, $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$0	$$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$1	gapless

For $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$2 and $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$3, one proves

$$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$4

and, in particular, for the half-chain length $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$5,

$$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$6

which is a full volume law (Zhang et al., 2017). The tensor-network description of this regime is the “rainbow phase,” in which dominant configurations are concentric arcs crossing the bipartition extensively (Alexander et al., 2018).

The static correlation structure is also unconventional. In the undeformed spin-$$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$7 chain, the connected longitudinal correlator $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$8 decays to zero at large separation, but the connected transverse correlator $$H_j=(1-F_{j,j+1,j+2})+(1-\sigma^x_{j+2}F_{j+2,j+1,j}\sigma^x_{j+2}),$$9 does not: at the edge,

$P_j^\pm=(1\pm \sigma_j^z)/2$0

and deep in the bulk the corresponding limit is $P_j^\pm=(1\pm \sigma_j^z)/2$1 (Dell'Anna et al., 2019). Accordingly, the cluster decomposition property holds along the $P_j^\pm=(1\pm \sigma_j^z)/2$2 axis and fails in the $P_j^\pm=(1\pm \sigma_j^z)/2$3 plane (Dell'Anna et al., 2019).

For colorful Fredkin chains, long-distance entanglement persists between disjoint subsystems. For two edge spins, the negativity is nonzero for $P_j^\pm=(1\pm \sigma_j^z)/2$4 and has the asymptotic limit

$P_j^\pm=(1\pm \sigma_j^z)/2$5

while for $P_j^\pm=(1\pm \sigma_j^z)/2$6 the reduced state is PPT-entangled. The edge mutual information remains finite for all $P_j^\pm=(1\pm \sigma_j^z)/2$7, whereas in the colorless case $P_j^\pm=(1\pm \sigma_j^z)/2$8 it vanishes (Dell'Anna, 2019). This is consistent with the failure of cluster decomposition in transverse observables.

4. Spectral gap, transport, and anomalous dynamics

The original Fredkin chain is gapless in the thermodynamic limit, but its finite-size gap closes unusually slowly. For the balanced sector of the spin chain on $P_j^\pm=(1\pm \sigma_j^z)/2$9 sites, rigorous bounds give

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$0

hence

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$1

with the lower bound derived via a mapping to Glauber dynamics on Dyck paths and the Diaconis–Saloff-Coste comparison theorem, and the upper bound obtained from a Brownian-excursion variational state (Movassagh, 2016).

In the deformed chain, the gap can close much faster. For the uniform $F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$2, $F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$3 phase,

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$4

a super-exponential upper bound in $F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$5 (Zhang et al., 2017). In the colorless $F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$6, $F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$7 case one has

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$8

which is exponentially small even though the entanglement remains bounded (Zhang et al., 2017). A related one-parameter deformation proved

$F_i=U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},$9

and established a domain-wall magnetization profile with $U_i=\tfrac12(1+\sigma_i^z)$0 and $U_i=\tfrac12(1+\sigma_i^z)$1 on the left half, $U_i=\tfrac12(1+\sigma_i^z)$2 on the right half for spin-$U_i=\tfrac12(1+\sigma_i^z)$3 (Udagawa et al., 2017).

Direct dynamical simulations reveal slow relaxation with a large dynamical exponent. Quantum Monte Carlo and finite-size scaling studies reported

$U_i=\tfrac12(1+\sigma_i^z)$4

interpreting the slow mode as a “canted bond” executing a constrained random walk through a Dyck-word background (Adhikari et al., 2020). More recently, a pseudolocal matched-parentheses conserved charge was identified in the periodic Fredkin chain, yielding an exact variational upper bound

$U_i=\tfrac12(1+\sigma_i^z)$5

and therefore a rigorous hydrodynamic lower bound

$U_i=\tfrac12(1+\sigma_i^z)$6

for Fredkin-constrained transport (McCarthy et al., 2024).

The response to local perturbations is also anisotropic. After a local $U_i=\tfrac12(1+\sigma_i^z)$7 quench, the change in $U_i=\tfrac12(1+\sigma_i^z)$8 connected correlators shows a conventional light cone, whereas the $U_i=\tfrac12(1+\sigma_i^z)$9 connected correlators respond almost immediately at all distances, reflecting pre-existing long-range transverse correlations (Dell'Anna et al., 2019).

Adiabatic driving across the deformed critical point at $N$00 gives a further anomaly. In that setting, finite-size scaling of the first gap yields

$N$01

with

$N$02

and the excitation density obeys

$N$03

rather than the standard Kibble–Zurek form (Francica et al., 2022).

5. Deformations, generalized descendants, and nonergodic sectors

Several families interpolate away from the original frustration-free point. A two-parameter generalization with projector

$N$04

has the exactly solvable $N$05-deformed line on the unit circle

$N$06

where the ground state remains an area-weighted Dyck-word superposition (Adhikari et al., 2020). Outside that line, the same study identified ferromagnetic, antiferromagnetic, dimerized, and non-Dyck “cat’s-ear” regions, with exponentially small gaps in the non-Dyck domes according to exact diagonalization for small $N$07 (Adhikari et al., 2020).

A different interpolation

$N$08

connects the Fredkin point $N$09 to the ferromagnetic and antiferromagnetic Heisenberg chains at $N$10 (Adhikari et al., 2018). Its zero-temperature phase diagram contains an $N$11-FM region for $N$12, a dimer plus Ising-AF region for $N$13, and a Heisenberg-AF regime for $N$14, with a crossover between two dimer patterns at $N$15 (Adhikari et al., 2018). The Fredkin point is “brittle” in the sense that an infinitesimal $N$16 opens a gap and destroys the gapless frustration-free behavior (Adhikari et al., 2018).

Generalized Fredkin Hamiltonians also realize Hilbert-space fragmentation and exact scars. In one spin-$N$17 family,

$N$18

the dynamics preserve the order of an emergent sequence of red and blue particles, so the Hilbert space fractures into disconnected Krylov subsectors labeled by that fixed color string and the hole number (Langlett et al., 2021). Each subsector contains one exact closed-form zero-energy scar state for arbitrary $N$19, and the half-chain Rényi-2 entropy satisfies

$N$20

(Langlett et al., 2021).

A deformed stochastic Fredkin chain with parameters $N$21 exhibits a sharp dynamical crossover between fast thermalization and slow metastable dynamics near its stochastic point $N$22 (Causer et al., 2024). In the large-$N$23 limit, an emergent folded Hamiltonian fragments the Hilbert space into exponentially many disconnected pieces; exact diagonalization finds a set of $N$24 non-thermal eigenstates with low entanglement and anomalous local observables, while variational MPS shows that their maximum bipartite entanglement saturates with system size (Causer et al., 2024). This suggests that Fredkin-type kinetic constraints provide a direct route from frustration-free combinatorics to weak ergodicity breaking.

6. Periodic-chain symmetries, conservation laws, and integrability status

Under periodic boundary conditions, the Fredkin chain changes qualitatively. The periodic Hamiltonian can be written as

$N$25

with indices understood modulo $N$26 (Pronko, 24 Jul 2025). The ground state is no longer unique: the zero-energy subspace is spanned by vectors built from equivalence classes $N$27 of Dyck paths related by Fredkin moves, and its dimension is $N$28 for odd $N$29 and $N$30 for even $N$31 (Pronko, 24 Jul 2025).

A notable recent development is the construction of two nonlocal operators $N$32 commuting with the periodic Hamiltonian. Acting on the ground-state multiplet, they serve as raising and lowering operators in $N$33, and the Lie algebra generated by nested commutators of $N$34 closes on

$N$35

with rank $N$36 (Pronko, 24 Jul 2025). In even length, the decomposition of total $N$37 into Cartan generators plus a central element coincides with a corresponding formula proposed for the periodic Motzkin chain (Pronko, 24 Jul 2025).

These nonlocal symmetries do not imply Yang–Baxter integrability. A separate analysis showed that, under both periodic and open boundary conditions, the Fredkin spin chain has no nontrivial local conserved charges beyond the Hamiltonian itself; the same conclusion holds for truncated variants obtained by deleting subsets of the three-site interaction terms (Fan et al., 5 Sep 2025). The proof proceeds by an exhaustive Pauli-string commutator analysis using a column method to exclude $N$38-local charges for all $N$39 under periodic boundaries and boundary-supported local charges under open boundaries (Fan et al., 5 Sep 2025). The present status is therefore that the Fredkin chain is partially solvable in its ground-state sector, but nonintegrable in the local-conservation-law sense.

The coexistence of exact frustration-free structure, combinatorial solvability, anomalous entanglement, subdiffusive transport bounds, fragmentation phenomena, and nonintegrability makes the Fredkin spin chain a distinctive class of constrained quantum many-body systems rather than a conventional integrable spin chain.