---
title: Anyonic Tight-Binding Lattice Model
url: https://www.emergentmind.com/topics/anyonic-tight-binding-model
type: topic
---

# Anyonic Tight-Binding Lattice Model

Searching arXiv for recent and foundational papers on anyonic tight-binding lattice models.
An anyonic tight-binding model is a lattice model in which particles propagate by nearest-neighbor hopping on a discrete lattice while the microscopic fields obey fractional exchange algebra, so that the kinetic term is intrinsically statistical rather than a free bosonic or fermionic hopping problem. In one-dimensional realizations, this structure appears either directly in number-conserving Fock parafermion chains or, after a fractional Jordan–Wigner transformation, as occupation-dependent Peierls phases and correlated hopping in anyon-Hubbard systems. Two paradigmatic formulations are the Fock-parafermion chain, where a formally quadratic Hamiltonian is nevertheless nonintegrable and supports statistics-induced bound states, and the extended anyon-Hubbard chain at unit filling, where the statistical angle reorganizes the bosonized low-energy theory and the commensurate phase diagram [1811.09422] [2410.00089].

## 1. Operator algebra and the meaning of “anyonic” hopping

The defining feature of an anyonic tight-binding model is not merely lattice discreteness but the fact that exchange statistics are encoded directly in the operator algebra. In the parafermionic formulation, one starts from parafermion operators $\{\hat\gamma_j\}$ of order $p$ obeying
$$
\hat\gamma_j \hat\gamma_l=\omega^{\mathrm{sgn}(j-l)}\hat\gamma_l\hat\gamma_j,\qquad \omega=e^{2\pi i/p},
$$
together with
$$
\hat\gamma_j^p=1,\qquad \hat\gamma_j^\dagger=\hat\gamma_j^{p-1}.
$$
For $p=2$ these reduce to Majorana operators. The particle-like degrees of freedom are Fock parafermions $\hat F_j,\hat F_j^\dagger$, for which each site has a local $p$-dimensional Fock space
$$
|m_j\rangle=\hat F_j^{\dagger m}|0\rangle,\qquad 0\le m\le p-1,
$$
and onsite number operator
$$
\hat N_j=\sum_{l=1}^{p-1}\hat F_j^{\dagger l}\hat F_j^l.
$$
The total number $\hat N=\sum_j\hat N_j$ is conserved, giving the model a genuine $U(1)$ symmetry. On different sites, Fock parafermions satisfy
$$
\hat F_j\hat F_l=\omega^{\mathrm{sgn}(l-j)}\hat F_l\hat F_j,\qquad
\hat F_j^\dagger \hat F_l=\omega^{-\mathrm{sgn}(l-j)}\hat F_l\hat F_j^\dagger,
$$
or equivalently
$$
\hat F_j\hat F_l=e^{i\pi\kappa\,\mathrm{sgn}(l-j)}\hat F_l\hat F_j,
\qquad \kappa=\frac{2}{p}.
$$

In the anyon-Hubbard formulation, the lattice anyon operators satisfy the deformed exchange relation
$$
\hat a_l \hat a_{l'}^\dagger - e^{i\theta\,\mathrm{sgn}(l-l')} \hat a_{l'}^\dagger \hat a_l =0
\qquad (l\neq l').
$$
Here $\theta=0$ gives the bosonic limit, while $\theta=\pi$ gives “pseudo-fermions”: fermionic exchange for $l\neq l'$, but not ordinary fermions onsite. The exact fractional Jordan–Wigner transformation
$$
\hat a_l = \hat b_l \exp\!\left(i\theta \sum_{j<l}\hat n_j\right)
$$
maps the anyons to bosons, and the kinetic term becomes
$$
\hat H_{\rm kin} = -J\sum_l \left[\hat b_l^\dagger \hat b_{l+1} e^{i\theta \hat n_l} +\text{H.c.} \right].
$$
This occupation-dependent phase is the core tight-binding manifestation of anyonic statistics: the exchange algebra is traded for a dynamical gauge phase attached to hopping rather than a removable global phase [1811.09422] [2410.00089].

## 2. The number-conserving Fock-parafermion chain

The simplest number-conserving anyonic tight-binding model studied for Fock parafermions is the one-dimensional nearest-neighbor chain
$$
\hat H=-t\sum_j\left(\hat F_j^\dagger \hat F_{j+1}+\hat F_{j+1}^\dagger \hat F_j\right),
\qquad t>0,
$$
usually with $t=1$. Although formally quadratic in $\hat F_j$, this Hamiltonian is not a free-particle problem for $p>2$. It conserves $\hat N$, is not inversion invariant because the FP algebra is oriented, is not time-reversal invariant, and is not particle-hole symmetric; however, it is invariant under the combined transformation $\hat F_j\to \hat F^\dagger_{-j}$. Because of this combined symmetry, densities above $(p-1)/2$ can be related to lower ones.

The single-particle sector has the usual tight-binding dispersion
$$
\varepsilon(k)=-2t\cos k,\qquad k=\frac{2\pi m}{L},
$$
with Fourier modes
$$
\hat F_k=\frac{1}{\sqrt L}\sum_j e^{ikj}\hat F_j.
$$
For $N>1$, however, the momentum-space operators are not again Fock parafermions and do not satisfy a simple canonical algebra. The formal expression $\hat H=\sum_k \varepsilon(k)\hat F_k^\dagger\hat F_k$ is therefore misleading, because
$$
[\hat F_k^\dagger \hat F_k,\hat F_q^\dagger \hat F_q]\neq C_k\delta_{kq}.
$$
From the viewpoint of ordinary bosons or fermions, the exchange strings generate effective interactions even though the Hamiltonian is bilinear in $\hat F$.

The generalized Jordan–Wigner or Fradkin–Kadanoff transformation makes this locality structure explicit. Introducing onsite Weyl hard-core boson operators $\hat B_j$ and diagonal phase operators
$$
\hat U_j=\mathrm{diag}(1,\omega,\omega^2,\dots,\omega^{p-1}),
$$
one writes
$$
\hat F_j=\left(\prod_{k=1}^{j-1}\hat U_k\right)\hat B_j,
$$
so that
$$
\hat F_j^\dagger \hat F_{j+1}=\hat B_j^\dagger \hat U_j \hat B_{j+1},
$$
and hence
$$
\hat H=-t\sum_j\left[\hat B_j^\dagger \hat U_j \hat B_{j+1} +\hat U_j^\dagger \hat B_j \hat B_{j+1}^\dagger\right].
$$
The mapped representation is local and suitable for DMRG, but the hopping amplitude is now occupation dependent through $\hat U_j$: quadratic anyonic hopping becomes a correlated-hopping problem of conventional commuting onsite operators.

A further structural distinction is the relation between Fock parafermions and fractionalized fermions. For even $p=2m$, the clustered operator
$$
\hat f_j=\hat F_j^m
$$
obeys canonical fermion relations,
$$
\{\hat f_j,\hat f_l\}=0,\qquad \{\hat f_j,\hat f_l^\dagger\}=\delta_{jl},\qquad \hat f_j^2=0.
$$
Thus one ordinary fermion can be viewed as a bound cluster of $m$ parafermions; for $p=6$, $\hat f_j=\hat F_j^3$. This is why the $p=6$ chain is interpreted as a model of fractionalized fermions [1811.09422].

## 3. Statistics-induced spectral structure and evidence for nonintegrability

The two-body problem already displays nontrivial consequences of anyonic statistics. With basis states
$$
|n_1,n_2\rangle=\hat F_{n_1}^\dagger \hat F_{n_2}^\dagger|0\rangle,
$$
the wavefunction can be written as
$$
|\Psi\rangle=\sum_{1\le n_1\le n_2\le L} a(n_1,n_2)|n_1,n_2\rangle,
$$
with scattering form
$$
a(n_1,n_2)=
\begin{cases}
A\left(e^{i(kn_1+qn_2)}+e^{i\theta}e^{i(qn_1+kn_2)}\right), & n_1<n_2,\\
B\,e^{i(k+q)n_1}, & n_1=n_2.
\end{cases}
$$
The energy remains additive,
$$
E(k,q)=-2t\cos k-2t\cos q,
$$
but periodic boundary conditions now impose statistically twisted quantization,
$$
k=\frac{2\pi}{L}\left(\lambda_k-\frac1p\right)-\frac{\theta}{L},\qquad
q=\frac{2\pi}{L}\left(\lambda_q-\frac1p\right)+\frac{\theta}{L}.
$$
The momenta are thus shifted by a statistical twist of order $\kappa\pi/L$, analogous to twisted boundary conditions in continuum anyon models.

Exact diagonalization reveals a stronger effect: branches of two-body bound states in addition to the scattering continuum. For $p=3$ they appear as two “mustaches” outside the continuum in $E$ versus center-of-mass momentum $K=k+q$. Their hallmark is that the wavefunction weight $a(n_1,n_2)$ is exponentially localized in the relative coordinate $d=n_2-n_1$, strongly peaked at $d=0$, which requires complex $k,q$. The binding is strongest near $K=\pi$ and weakens as the branch approaches the continuum lobes. The physical origin is purely statistical: there is no explicit interaction term, but anyonic exchange phases effectively generate attraction in certain channels. As $p$ increases and $\kappa\to 0$, the bound-state branches become less visible and merge into the continuum.

The nonintegrability claim for the quadratic chain is numerical rather than a theorem. The evidence comes from level-spacing statistics at fixed particle number with open boundary conditions. Integrable systems should show Poisson statistics,
$$
P_{\rm P}(s)=e^{-s},
$$
whereas nonintegrable systems show Wigner–Dyson level repulsion,
$$
P_{\rm WD}(s)\sim A s^\beta e^{-Bs^2},
\qquad
P_{\rm WD}(s\to0)\sim s^\beta,\quad \beta>0.
$$
For $p=3,6$, increasing system size yields clear level repulsion; for $p=6$ the statistics converge rapidly toward a GOE-like form, while $p=3$ also shows repulsion although finite-size convergence is slower. By contrast, $p=2$ gives Poisson-like statistics, as expected for an integrable free-fermion chain. The paper also notes that a small hopping inhomogeneity can drive a GOE-to-GUE crossover, consistent with the complex nature of the Hamiltonian and symmetry sensitivity of random-matrix classes. The resulting conceptual point is precise: a quadratic anyonic model need not be integrable because the exchange algebra generates effective many-body interactions [1811.09422].

## 4. Gapless regimes, anyonic correlations, and fractionalized fermions

The many-body problem of the Fock-parafermion chain was analyzed by DMRG on systems up to $L=288$ with open boundary conditions. In the cases studied, the neutral gap $\Delta_0$ at fixed $N$ scales as $\Delta_0\sim L^{-1}$ except for $p=3$ at unit filling $N/L=1$, where it remains finite, indicating a gapped phase. Entanglement entropy in the gapless regimes fits
$$
S(\rho_\ell)=a+\frac{c}{6}\log\!\left[\frac{2L}{\pi}\sin\!\left(\frac{\pi\ell}{L}\right)\right]
$$
with $c=1$, establishing a $c=1$ Luttinger-liquid universality class. The gapped $p=3$, $N/L=1$ state instead obeys an area law and is interpreted as an “anyonic Mott-like” phase.

The basic one-body anyonic correlator is
$$
G_1(j,l)=\langle \Psi_{\rm GS}|\hat F_j^\dagger \hat F_l|\Psi_{\rm GS}\rangle.
$$
In the local $\hat B_j,\hat U_j$ representation it becomes a string correlator,
$$
G_1(j,l)=\left\langle \hat B_j^\dagger \big(\hat U_j\cdots \hat U_{l-1}\big)\hat B_l\right\rangle
\qquad (j<l),
$$
which is the direct lattice manifestation of anyonic statistics. In the gapless phases,
$$
|G_1(x,x+r)|\sim r^{-\alpha_1}.
$$
For $p=3$, $\alpha_1$ stays around $0.67$-$0.74$ over a broad density range; for $p=6$, $\alpha_1$ decreases with density from about $0.46$ to $0.34$. In the gapped $p=3$, $N/L=1$ case, $G_1$ decays exponentially as $e^{-r/\xi}$.

These exponents were compared with the anyonic Luttinger-liquid theory of Calabrese, Mintchev, and Bellazzini, according to which
$$
|G_1(x,x+r)|\sim r^{-(\kappa^2 K+1/K)/2}.
$$
For noninteracting anyons, $K=\kappa^{-1}$ and therefore $\alpha_1=\kappa=2/p$. The agreement improves at larger densities. The same framework predicts for
$$
G_2(j,l)=\langle \hat F_j^{\dagger 2}\hat F_l^2\rangle
$$
the decay
$$
|G_2(x,x+r)|\sim r^{-\alpha_2},\qquad \alpha_2=4\alpha_1=\frac{8}{p}.
$$
For $p=6$ the agreement is good; for $p=3$ there are visible deviations, heuristically ascribed to truncation effects from the very small onsite Hilbert space.

The $p=6$ case is distinguished by the fermionic cluster operator $\hat f_j=\hat F_j^3$. Its correlator
$$
G_3(j,l)=\langle \Psi_{\rm GS}|\hat f_j^\dagger \hat f_l|\Psi_{\rm GS}\rangle
$$
also decays algebraically,
$$
|G_3(x,x+r)|\sim r^{-\alpha_3},
$$
with Luttinger-liquid-inspired prediction
$$
\alpha_3=9\alpha_1=\frac{18}{p},
$$
so that $\alpha_3=3$ for $p=6$. Numerically the fitted exponents approach this value at higher densities. This means that although $\hat f_j$ are genuine fermion operators, their two-point function decays roughly as $r^{-3}$ rather than as $r^{-1}$ for free fermions.

The momentum-space signatures are likewise anyonic. The anyonic momentum distribution
$$
n_0(k)=\frac1L\sum_{j,l} e^{ik(j-l)}\,G_1(j,l)
$$
is not symmetric under $k\to-k$ and has a peak at positive momentum. The expected peak position is
$$
k_{\max}=\kappa \pi \frac{N}{L}.
$$
This relation is motivated from twisted momentum quantization, from the exchange-phase accumulation in $G_1$, and from Luttinger-liquid theory, though significant deviations appear for $N/L\gtrsim 1/4$. The fermionic momentum distribution
$$
n_1(k)=\frac1L\sum_{j,l} e^{ik(j-l)}G_3(j,l)
$$
is again asymmetric under $k\to-k$ and peaked at nonzero momentum, unlike the standard free-fermion case.

For interpretation, the $p=6$ chain was compared phenomenologically with two counter-propagating Laughlin edge modes at filling $\nu=1/3$. The low-energy Hamiltonian was taken as
$$
\hat H=\frac{v}{2}\int \left[(\partial_x \hat\theta)^2+(\partial_x \hat\phi)^2\right]\,dx,
$$
with
$$
[\hat\phi(x_1),\hat\theta(x_2)]= i\frac{2\pi}{p}\,\Theta_{\rm H}(x_2-x_1).
$$
The correspondence is not exact: in the Laughlin-edge theory right and left movers carry opposite statistical parameter and the full theory is time-reversal invariant, whereas the lattice chain has both chiralities but only one statistical parameter and is not time-reversal invariant. Still, the scaling of the fermionic correlator for $p=6$ at larger density agrees quantitatively with Wen’s hydrodynamic prediction for electron correlations at the edge of a $\nu=1/3$ Laughlin state [1811.09422].

## 5. Extended anyon-Hubbard chains and the interacting phase diagram

A second major realization of the anyonic tight-binding idea is the one-dimensional extended Hubbard model for lattice anyons,
$$
\hat H = \sum_l \left[ - J\, \hat a_l^\dagger \hat a_{l+1}^{\phantom\dagger} +\text{H.c.} +\frac{U}{2}\hat n_l(\hat n_l-1) +V\,\hat n_l\hat n_{l+1} \right],
$$
analyzed at average unit filling $\langle \hat n_l\rangle=1$ and, for the main bosonization treatment, under the local constraint $n_l\le 2$. After the fractional Jordan–Wigner map, the kinetic term acquires the density-dependent gauge phase
$$
\hat H_{\rm kin} = -J\sum_l \left[ \hat b_l^\dagger \hat b_{l+1} e^{i\theta \hat n_l} +\text{H.c.} \right].
$$
At unit filling with maximal occupancy two, the local Hilbert space is mapped to spin-1 through
$$
\hat S_l^z = -\hat n_l +1,
$$
and the kinetic term takes the occupation-dependent form
$$
\hat H_{\rm kin} = -\frac12\sum_l \left[ t(\hat S_l^z,\hat S_{l+1}^z)\, \hat S_l^- \hat S_{l+1}^+ +\text{H.c.} \right].
$$
The nontrivial amplitudes are
$$
t(-1,0)=2J e^{i\theta},\qquad t(0,1)=J,\qquad t(-1,1)=J\sqrt2\,e^{i\theta},\qquad t(0,0)=J\sqrt2.
$$
This is a concrete microscopic realization of hopping in a dynamical lattice gauge field.

Bosonization yields a two-channel sine-Gordon theory in terms of $\Phi_\pm$ and $\Theta_\pm$, with low-energy Hamiltonian containing
$$
\Delta\left(\partial_x\Theta_+\partial_x\Phi_+ +\partial_x\Theta_-\partial_x\Phi_-\right),
\qquad
\Delta=\frac{2Ja}{\pi}\sin\theta.
$$
The angle $\theta$ therefore enters not only through the cosine perturbations and Luttinger parameters $K_\pm$, but also through a current-density coupling that produces left-right asymmetry in dynamics. A naively allowed perturbation $\sin(\sqrt{4\pi}\Phi_+)$ is absent because a modified inversion symmetry $\tilde{\mathcal I}= \hat U(\theta)\mathcal K\mathcal I$ forbids it, protecting a full $c=1$ Gaussian transition line.

At unit filling, the phase diagram contains four gapped insulating phases and two gapless regimes. The insulating phases are characterized by distinct pinning patterns and order parameters:

| Phase | Pinned fields | Diagnostic |
|---|---|---|
| Mott insulator (MI) | $\Phi_+=0,\ \Theta_-=0$ | $\mathcal O_{\rm MOTT}$ |
| Charge-density wave (CDW) | $\Phi_+=\frac{\sqrt\pi}{2},\ \Phi_-=0$ | $\mathcal O_{\rm CDW}$ |
| Haldane insulator (HI) | $\Phi_+=\frac{\sqrt\pi}{2},\ \Theta_-=0$ | $\mathcal O_{\rm HI}$ |
| Dimerized state | $\Phi_+=0,\ \Phi_-=0$ | $\mathcal O_{\rm DIMER}$ |

The corresponding diagnostics are
$$
\mathcal O_{\rm MOTT} = \left\langle e^{i\pi\sum_{i<k<j}\delta n_k} \right\rangle,
\qquad
\mathcal O_{\rm CDW} = (-1)^{|i-j|} \langle \delta\hat n_i \delta\hat n_j\rangle,
$$
$$
\mathcal O_{\rm HI} = \left\langle \delta\hat n_i e^{i\pi\sum_{i\le k<j}\delta\hat n_k} \delta\hat n_j \right\rangle,
\qquad
\mathcal O_{\rm DIMER} = \left\langle \hat H^{\rm bond}_{\rm even} - \hat H^{\rm bond}_{\rm odd} \right\rangle,
$$
with $\delta \hat n_l=\hat n_l-1$. For small $V$ there are also superfluid (SF) and pair-superfluid (PSF) regions; the paper identifies them numerically and does not give a full bosonization description.

The bosonization predicts two principal critical lines. A Gaussian $c=1$ line occurs when $g_1=0$, yielding the weak-coupling estimate
$$
U_{c=1}^c(\theta) \sim 2V+\frac{4\,\mathrm{Re}[J_{zz}(\theta)]}{\pi}.
$$
An Ising $c=\frac12$ line arises from competition between $g_2\cos(\sqrt{4\pi}\Phi_-)$ and $g_3\cos(\sqrt{4\pi}\Theta_-)$ near $K_-=1$, with estimate
$$
U_{\mathrm{Ising}}^c(\theta) \sim 2V -\frac{5\,\mathrm{Re}[J_{zz}(\theta)]}{2\pi} -3\pi\,\mathrm{Re}[\tilde J(\theta)].
$$
Numerically the transition lines are fit as
$$
U_{c=1}^c = -0.34 + 2.02\cos\theta,
\qquad
U_{\mathrm{Ising}}^c = -3.51 -1.73\cos\theta.
$$
Their opposite trends explain why the Haldane insulator shrinks and the dimerized phase grows as $\theta\to\pi$. For each fixed $V$, the MI, HI, CDW, and dimer phases meet at a common crossing point, interpreted as a multicritical line in $(U,V,\theta)$ space with bosonization suggestion $c=\frac32$.

The numerical analysis combines finite DMRG, infinite DMRG, and exact diagonalization for level spectroscopy. Phase boundaries and universality classes are identified through order parameters, fidelity susceptibility, entanglement entropy, correlation length, and level crossings. The reported critical behavior is: HI–CDW Ising with $c=1/2$ and exponents $\nu\approx 1.03$, $\beta\approx 0.247$; dimer–MI also Ising at large $\theta$; HI–MI Gaussian with $c=1$; CDW–dimer continuous with $c=1$ and exponents close to $\nu\approx 1.5$, $\beta\approx 0.5$; SF/PSF boundaries of BKT type; and MI–CDW first order at large $V$ [2410.00089].

## 6. Scope, misconceptions, and relation to topological lattice modeling

Several distinctions are central to the subject. First, an anyonic tight-binding model is not defined by nontrivial band topology alone. The parafermion chain and the extended anyon-Hubbard model are anyonic because fractional statistics are present in the microscopic algebra or, equivalently after mapping, in the density-dependent gauge structure of the hopping. This is why a formally quadratic chain can be nontrivial, and why in the interacting anyon-Hubbard case the statistical angle reshapes the competition among commensurate instabilities. A useful concise formulation is that the characteristic lattice signature of one-dimensional anyonic tight binding is statistical transmutation encoded as local, occupation-dependent phases in nearest-neighbor hopping [1811.09422] [2410.00089].

Second, not every tight-binding model relevant to fractionalized many-body physics is itself an anyonic tight-binding model. The three-orbital moiré model for twisted transition metal dichalcogenide homobilayers is explicitly presented as a strictly local tight-binding basis with nontrivial band Chern numbers, short-range interactions, and complex symmetry-fixed phases, and it is discussed as a promising parent lattice model for fractional Chern insulators. At the same time, it explicitly states that anyons are not directly studied, and that it does not derive topological order, quasiparticle braiding statistics, emergent gauge fields, or anyonic excitations. Its relevance is therefore indirect: it provides a local topological lattice framework from which anyonic phases might emerge after interactions, rather than realizing anyonic exchange algebra at the microscopic level [2403.15546].

Third, the evidentiary status of different claims is not uniform. In the parafermion chain, the operator algebra, Fradkin–Kadanoff map, one-particle dispersion, and existence of canonical clustered fermions for even $p$ are exact. The two-body scattering ansatz and twisted quantization are analytical, while the existence and structure of bound states are supported by exact diagonalization and wavefunction analysis. Nonintegrability is strongly supported by level-spacing statistics rather than rigorously proven, the $c=1$ low-energy description is numerically established through gap scaling and entanglement, and the comparison to Laughlin-edge physics is heuristic or phenomenological rather than a microscopic derivation. In the extended anyon-Hubbard chain, the bosonized field theory and symmetry analysis organize the phase structure, while the phase boundaries and universality classes are determined numerically.

Within condensed-matter theory, the significance of the anyonic tight-binding model is twofold. It isolates the effect of anyonic statistics without requiring explicit quartic interactions, showing that exchange and occupancy constraints alone can generate bound states, nonintegrability, asymmetric momentum distributions, and unconventional correlation exponents. It also furnishes a lattice laboratory for fractionalized fermions and for parafermion platforms related to fractional quantum Hall edges proximitized by superconductors and magnets. The general lesson is that in anyonic lattice systems, “quadratic” does not imply free behavior, and “tight binding” does not reduce the problem to ordinary single-particle band theory.

Source: https://www.emergentmind.com/topics/anyonic-tight-binding-model