Papers
Topics
Authors
Recent
Search
2000 character limit reached

Fredkin Spin Chain

Updated 7 July 2026
  • Fredkin chain is a one-dimensional spin-½ quantum model with conditional three-site interactions and ground states organized by Dyck paths.
  • It employs Fredkin gate representations and controlled-SWAP logic to realize a frustration-free Hamiltonian with anomalous entanglement scaling.
  • Deformations and colored extensions reveal varied gap behaviors, phase transitions, and hidden Lie-algebra symmetries in the system.

The Fredkin chain is a one-dimensional spin-12\tfrac12 quantum chain with interaction of three nearest neighbors, introduced as a frustration-free model whose Hamiltonian can be expressed in terms of Fredkin gates and whose exact ground-state structure is organized by Dyck paths rather than Motzkin paths (Salberger et al., 2016). In its open-boundary form, the model has an exactly solvable Dyck-path ground state; in deformed, colored, stochastic, and periodic variants it exhibits area-weighted path ensembles, anomalous entanglement scaling, polynomially or exponentially small gaps depending on regime, slow constrained dynamics, and, for periodic boundary conditions, a hidden Lie-algebraic symmetry acting on a degenerate ground-state manifold (Pronko, 24 Jul 2025).

1. Definition and local interaction structure

In the basic spin-12\tfrac12 formulation, each site carries basis states |\uparrow\rangle and |\downarrow\rangle. The standard open-chain bulk Hamiltonian is written as

Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},

with

Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),

and

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),

the two-site singlet projector. The open chain is completed by boundary fields

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,

which pin the left edge up and the right edge down (Salberger et al., 2016, Adhikari et al., 2020).

The same interaction admits a Fredkin-gate representation. If Fi,j,kF_{i,j,k} is the controlled-SWAP gate acting with control on site ii, then one local term can be rewritten as

12\tfrac120

making explicit the relation between the chain and reversible controlled-swap logic (Salberger et al., 2016).

A closely related periodic formulation uses the three-site density

12\tfrac121

where

12\tfrac122

and

12\tfrac123

with site labels understood modulo 12\tfrac124 (Pronko, 24 Jul 2025). What distinguishes the Fredkin chain from standard nearest-neighbor chains is precisely this conditional three-site structure: a spin-0 component of two adjacent spins interacts with a third neighboring spin (Pronko, 24 Jul 2025).

2. Dyck paths, Fredkin moves, and exact ground states

The combinatorial description is obtained by identifying 12\tfrac125 with an up-step and 12\tfrac126 with a down-step. Equivalently, one may identify

12\tfrac127

so that spin strings become parenthesis words and lattice walks. A Dyck path is then a walk from 12\tfrac128 to 12\tfrac129 that never goes below the axis; the corresponding spin words are balanced parenthesis strings with nonnegative prefix height (Salberger et al., 2016).

The local equivalence relation is generated by the Fredkin moves

|\uparrow\rangle0

These moves allow an adjacent matched pair |\uparrow\rangle1 to move locally through neighboring unmatched parentheses. Basis states therefore decompose into equivalence classes |\uparrow\rangle2, and the open-chain boundary terms select the Dyck class |\uparrow\rangle3 (Salberger et al., 2016).

For even open chains, the exact ground state is the equal-weight superposition of all Dyck words,

|\uparrow\rangle4

with |\uparrow\rangle5 (Salberger et al., 2016). The model is frustration free because each local projector annihilates this state separately.

The same construction extends to colored higher-spin chains. With local Hilbert space |\uparrow\rangle6, one obtains properly colored Dyck paths, in which every matched up/down pair carries the same color, and the unique ground state becomes the equal-weight sum over all such properly colored Dyck paths. In this formulation the paper constructs an |\uparrow\rangle7-symmetric model, and the leading term of the entanglement entropy is proportional to the square root of the length of the lattice (Salberger et al., 2016).

The Dyck state also admits an exact MPS-like representation. Using infinite matrices

|\uparrow\rangle8

the amplitude of a word |\uparrow\rangle9 is |\downarrow\rangle0, which is nonzero exactly on Dyck words (Salberger et al., 2016).

3. Deformations, colors, and entanglement structure

A major generalization is the deformed Fredkin chain, in which the equal-weight Dyck superposition is replaced by an area-weighted one. In the uniform colored case the Hamiltonian is

|\downarrow\rangle1

and the exact ground state is

|\downarrow\rangle2

where |\downarrow\rangle3 is the area under the Dyck path |\downarrow\rangle4 (Salberger et al., 2016, Zhang et al., 2017). At |\downarrow\rangle5 this reduces to the undeformed Fredkin chain; for |\downarrow\rangle6, high-area paths dominate; for |\downarrow\rangle7, low-area paths dominate.

The deformation can be made site dependent while preserving frustration freeness. In the multi-parameter family, the local amplitudes |\downarrow\rangle8, |\downarrow\rangle9 must satisfy

Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},0

and the ground state is again an exact weighted superposition over colored Dyck paths (Zhang et al., 2017). A two-parameter reformulation writes the active two-site operator as

Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},1

so that the Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},2-deformed family lies on the unit circle

Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},3

with the original Fredkin point at Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},4 (Adhikari et al., 2020).

The entanglement structure depends sharply on Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},5 and on the number of colors. For the uniform model, the half-chain entropy is bounded for Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},6. At Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},7, the known undeformed results are logarithmic scaling for Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},8 and Hbulk=i=2N1Hi,Hi=Ui1Pi,i+1+Pi1,iDi+1,H_{\rm bulk}=\sum_{i=2}^{N-1} H_i, \qquad H_i = U_{i-1}P_{i,i+1}+P_{i-1,i}D_{i+1},9 scaling for Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),0. For Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),1, the colored chain has

Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),2

so the entropy is linear in the half-chain length, while the colorless chain remains bounded (Zhang et al., 2017). In the earlier deformation paper this same regime is described as a transition to a rainbow-like extensively entangled phase, with

Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),3

for Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),4 and bounded entropy for Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),5 (Salberger et al., 2016).

For Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),6, the ground-state magnetization is purely along Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),7,

Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),8

and Ui=12(1+σiz),Di=12(1σiz),U_i=\tfrac12(\mathbb 1+\sigma_i^z),\qquad D_i=\tfrac12(\mathbb 1-\sigma_i^z),9 exhibits a domain-wall structure (Udagawa et al., 2017). In parallel, bulk tensor networks for colorful area-deformed Fredkin states can be built exactly on an upside-down step pyramid, using Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),0 identical four-index tensors; in the Fredkin specialization the flat-step tiles of the Motzkin construction are removed, leaving only the Dyck-arch sector (Alexander et al., 2018).

4. Gap structure, correlations, and slow dynamics

The original open Fredkin chain is gapless with a polynomially small gap. By mapping the Hamiltonian in the balanced Dyck sector to a reversible Markov chain on colored Dyck paths, the energy of the first excited state above the ground state is bounded from below by Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),1, while a variational twisted-area state gives the upper bound

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),2

so the gap is polynomially small in Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),3 (Movassagh, 2016).

In the deformed regime Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),4, variational arguments yield much smaller upper bounds. For the colored chain,

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),5

while for the colorless chain,

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),6

The latter provides an example of bounded entanglement entropy together with a vanishing spectral gap (Zhang et al., 2017). A related analysis of the same regime proves that the finite-size gap decays at least exponentially with system size,

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),7

and develops general bounds for von Neumann and Rényi entropies and for the entanglement spectrum for arbitrary bipartitions (Udagawa et al., 2017).

Static correlations are strongly anisotropic in the colorless chain. The connected longitudinal correlator satisfies

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),8

whereas the transverse one remains finite,

Pi,i+1=14(1σiσi+1),P_{i,i+1}=\frac14(\mathbb 1-\vec\sigma_i\cdot \vec\sigma_{i+1}),9

and deep in the bulk

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,0

A local Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,1 quench therefore produces a normal light cone in Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,2-based observables but no light cone in the transverse channel (Dell'Anna et al., 2019).

The dynamical exponent is unusually large. Direct projector Monte Carlo for the lowest Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,3 excitation gives

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,4

with a consistent independent estimate

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,5

and interprets the excitation as an excited bond executing a constrained random walk through a Dyck background (Adhikari et al., 2020). A rigorous hydrodynamic bound identifies a pseudolocal conserved matched-parentheses charge and proves

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,6

establishing subdiffusion in the Fredkin dynamical universality class (McCarthy et al., 2024).

For the deformed Fredkin critical point at Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,7, the gap scales with critical exponents

Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,8

in the Hboundary=11+NN,H_{\rm boundary}=|\downarrow_1\rangle\langle\downarrow_1|+|\uparrow_N\rangle\langle\uparrow_N|,9 sector. Standard Kibble–Zurek scaling would give Fi,j,kF_{i,j,k}0, but because the overlaps Fi,j,kF_{i,j,k}1 scale with level-dependent exponents, the paper derives and numerically confirms the anomalous law

Fi,j,kF_{i,j,k}2

instead (Francica et al., 2022).

5. Periodic chain, degeneracy, and hidden Lie symmetry

With periodic boundary conditions, the boundary constraints are removed, the local Fredkin moves acquire periodic extensions, and the ground state becomes highly degenerate. The ground states are organized by equivalence classes Fi,j,kF_{i,j,k}3 of Dyck paths, where a path starts at Fi,j,kF_{i,j,k}4, ends at Fi,j,kF_{i,j,k}5, has Fi,j,kF_{i,j,k}6, Fi,j,kF_{i,j,k}7, and must hit the Fi,j,kF_{i,j,k}8-axis at least once. The corresponding class state is

Fi,j,kF_{i,j,k}9

and every vector in ii0 has

ii1

For each ii2 there is exactly one zero-energy ground state with cyclic eigenvalue ii3,

ii4

When ii5 is even and ii6, there is one additional zero-energy ground state with ii7,

ii8

Hence the ground-state degeneracy is ii9 for odd 12\tfrac1200 and 12\tfrac1201 for even 12\tfrac1202 (Pronko, 24 Jul 2025).

A central result in the periodic case is the construction of two nonlocal operators

12\tfrac1203

equivalently

12\tfrac1204

which commute with the Hamiltonian,

12\tfrac1205

and act as raising and lowering operators on the cyclic-invariant ground states,

12\tfrac1206

away from the extremal 12\tfrac1207 sectors (Pronko, 24 Jul 2025). More generally, the paper proves

12\tfrac1208

for a family of charge-12\tfrac1209 operators 12\tfrac1210.

From repeated commutators of 12\tfrac1211, the author conjectures that the symmetry algebra generated on the ground-state manifold is

12\tfrac1212

In this description,

12\tfrac1213

where 12\tfrac1214 is central. In the even-12\tfrac1215 case this formula exactly coincides with the analogous formula previously conjectured for the periodic Motzkin chain (Pronko, 24 Jul 2025).

For even 12\tfrac1216, the periodic chain also possesses

12\tfrac1217

which commutes with 12\tfrac1218, anticommutes with the cyclic shift, and exchanges the two 12\tfrac1219 parity sectors. This operator is excluded from the main Lie-algebra analysis because it does not commute with the cyclic shift, but it explains the spectral doubling in the 12\tfrac1220 sector (Pronko, 24 Jul 2025).

Several later developments place the Fredkin chain in broader settings. A two-parameter extension with Hamiltonian terms 12\tfrac1221 embeds the 12\tfrac1222-deformed model on the unit circle 12\tfrac1223, separating an inner ferromagnetic region from outer antiferromagnetic and dimerized regimes; throughout most of the non-ferromagnetic parts of the phase diagram, the ground state has Dyck word form, while two outer regions favor Dyck-word mismatches and exhibit exponentially small level spacing (Adhikari et al., 2020).

Generalized Fredkin Hamiltonians can also display Hilbert-space fragmentation and exact many-body scars. In one family,

12\tfrac1224

the minus sign introduces destructive interference preventing spin configurations from connected and therefore fractures the Hilbert space into invariant sub-spaces; within each fractured sub-sector the authors construct an exact zero-energy eigenstate with logarithmic or area-law entanglement entropy (Langlett et al., 2021).

A different deformation away from the stochastic point shows a sudden change from fast thermalization to slow metastable dynamics, connected to emergent kinetic constraints and fragmentation in an effective folded model. In that setting the authors identify non-thermal eigenstates akin to quantum many-body scars, with area-law entanglement at large system size (Causer et al., 2024). On the classical side, a stochastic Fredkin chain with generator

12\tfrac1225

inherits the same Dyck-path sector, equilibrium three-phase structure, slow dynamics, and active-inactive large-deviation transition (Causer et al., 2022).

A thin-cylinder realization in fractional quantum Hall physics produces a deformed Fredkin model for the 12\tfrac1226 Moore–Read state. After truncation and restriction to the root-state-connected Hilbert-space fragment, the effective qubit Hamiltonian becomes

12\tfrac1227

with exact ground state

12\tfrac1228

In the thin-cylinder regime this model reproduces the Moore–Read entanglement properties and geometric-quench response (Voinea et al., 2023).

The integrability question has also been addressed directly. A rigorous analysis of local conserved charges concludes that the Fredkin spin chain, under both periodic and open boundary conditions, has no nontrivial local conserved quantities beyond the obvious trivial ones; partially truncated Fredkin chains remain nonintegrable unless all three-site terms are removed, reducing the model to the pure XXX chain (Fan et al., 5 Sep 2025).

Finally, the name should be distinguished from a very different classical usage. A one-dimensional elementary cellular automaton with memory can realize Fredkin logic via glider collisions on a chain of binary cells, but this is not the quantum many-body Fredkin spin chain (Martinez et al., 2018).

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Fredkin chain.