Translation-Invariant Free-Fermion Hamiltonians

Updated 5 December 2025

Translation-Invariant Free-Fermion Hamiltonians are operators on lattice fermion Fock spaces that remain invariant under lattice shifts, enabling explicit spectral resolution via Fourier transform.
They enforce frustration-free and locality constraints, leading to quadratic band touchings and excluding linear Dirac or Weyl dispersions.
Their algebraic and topological structures underpin robust bulk-boundary correspondence and have applications in quantum codes and entanglement scaling.

A translation-invariant free-fermion Hamiltonian is an operator on a lattice-fermion Fock space whose quadratic form and coefficients are translation-invariant under lattice shifts. These Hamiltonians play a central role in condensed matter, statistical mechanics, quantum information, and topological phases. Their structure allows both explicit spectral solutions by Fourier methods and powerful algebraic/topological classification. Recent work rigorously delineates their exact solvability, constraints deriving from symmetries and frustration-freeness, emergent edge behavior, and applications to codes and many-body entanglement scaling.

1. Algebraic Structure and Diagonalization

The generic translation-invariant quadratic (free) fermionic Hamiltonian on a Bravais lattice $\Lambda\subset\mathbb{Z}^d$ with $N$ -orbital unit cells is

$H = \sum_{R,R'\in\Lambda} c_{R}^\dagger h_{R-R'} c_{R'} = \sum_{k\in \mathrm{BZ}} c_k^\dagger h(k)\,c_k,$

where $c_{R}$ is a vector of fermionic annihilation operators, $h_{r}$ is a finite-range Hermitian hopping matrix, and the Bloch Hamiltonian is $h(k) = \sum_{r} e^{-ik\cdot r} h_r$ (Masaoka et al., 17 Mar 2025). In 1D, pairing terms can be incorporated, leading to

$H = \sum_{m,n}\frac{1}{2} \Psi_m^\dagger \begin{pmatrix} h_{m-n} & \Delta_{m-n} \ \Delta_{n-m}^* & -h_{n-m} \end{pmatrix} \Psi_n,$

with Nambu spinor $\Psi_n$ (Jones et al., 2022). Translation invariance always allows full spectral resolution via discrete Fourier transform; each momentum sector diagonalizes independently, yielding an explicit band structure.

2. Frustration-Free and Locality Constraints

A translation-invariant free-fermion Hamiltonian $H = \sum_{R} H_R$ is called frustration-free if there exists a (Slater-determinant) ground state $|\Phi\rangle$ such that $H_R|\Phi\rangle = 0$ for all $R$ (Masaoka et al., 17 Mar 2025, Ono et al., 18 Mar 2025). The necessary and sufficient algebraic condition is that, for all local positive (unoccupied) modes $\psi_{R,\alpha}$ and local negative (occupied) modes $\phi_{R,\beta}$ arising from diagonalizing $h_R$ , their anticommutator vanishes:

$\{\psi_{R,\alpha},\phi^{\dagger}_{R',\beta}\} = 0 \quad \forall R, R', \alpha, \beta.$

In momentum space, let $P(k)$ be the projector onto occupied (valence) bands. The frustration-free condition is then

$P(k)\,H(k)\,P(k) = 0 \quad \forall k.$

This algebraic constraint has strong implications for the allowed low-energy structure. Specifically, all band touchings in translation-invariant, frustration-free, local free-fermion systems are at least quadratic in momentum; linear (Dirac or Weyl) cones are excluded (Ono et al., 18 Mar 2025). This result excludes the realization of Dirac/Weyl semimetals within this framework.

3. Symmetry, Topology, and Critical Phenomena

Inversion symmetry plays a crucial role: translation-invariant quadratic Hamiltonians with inversion-symmetric dispersion relations $\varepsilon(k) = \varepsilon(-k)$ yield ground states with symmetric two-point correlators. If inversion is spontaneously broken in the ground state—for instance,

$\operatorname{Im} \langle a_j a^\dagger_k \rangle \neq 0,$

then gaplessness and algebraically decaying correlations are enforced (Kadar, 2016). Explicitly, the spectral gap $\Delta = \min_p |\varepsilon(p)|$ vanishes and the decay $\langle a_j a^\dagger_k \rangle \sim |j-k|^{-\alpha}$ follows.

A central topological invariant in 1D is the winding number, given by the argument of the complex symbol $f(z)=\sum_r t_r z^r$ , $z=e^{ik}$ , which counts the phase winding of $f(e^{ik})$ over the Brillouin zone,

$\omega = \frac{1}{2\pi}\oint_{0}^{2\pi} d k\,\partial_k \arg f(e^{ik})\in\mathbb{Z},$

with bulk-boundary correspondence: $\omega$ localized zero-energy edge modes (Jones et al., 2022).

4. Exact Solvability, Transfer Matrices, and Nonlinear Constructions

Not all exactly solvable translation-invariant free-fermion models are manifestly quadratic. The quantum chain of local four-Majorana terms

$H_{\mathrm{u}} = \sum_{m=1}^{2M} \psi_m \psi_{m+1} \psi_{m+3} \psi_{m+4}$

admits nonlocal, highly nonlinear raising and lowering operators obeying canonical Clifford algebra relations. Upon transformation, $H$ becomes diagonal in these nonlocal Bogoliubov modes and is isospectral to a free-fermion quadratic Hamiltonian (Fendley, 2019). The commuting family of classical transfer matrices

$T_M(u) = \sum_{s=0}^{\lfloor(M+2)/3\rfloor}(-u)^s Q^{(s)}$

encodes all conserved quantities, and the spectrum is obtained as all possible choices of sign for sums of single-particle energies. In the uniform case, the chain is gapless with a multicritical point at dynamical exponent $z=3/2$ and admits a supersymmetry generated by trilinear fermionic charges, resulting in exponential degeneracy of free-fermion multiplets.

5. Band-Touching Phenomena and Finite-Size Scaling

Quadratic (or softer) band touchings at points $k^*$ are universal in frustration-free, translation-invariant free-fermion systems with locality. The Bloch Hamiltonian structure enforces that, in a finite-size $L\times L$ Brillouin zone, the single-particle gap scales as $O(1/L^2)$ (Masaoka et al., 17 Mar 2025, Ono et al., 18 Mar 2025). The honeycomb-lattice frustration-free model realizes a quadratic node at $K,K'$ ; the many-body charge-neutral gap scales as $O(1/L^2)$ even in the presence of frustration-free interactions, while variational estimates using the Single-Mode Approximation produce an anomalously small $O(1/\ln L)$ pseudo-gap due to quantum-metric divergence at the node. This quadratic touchings phenomenon precludes realizing linear-dispersion (Dirac/Weyl) points in such models.

6. Bulk-Boundary Correspondence, Edge Modes, and Topological Features

For translation-invariant free-fermion Hamiltonians whose symbol function $f(z)$ has nontrivial winding, bulk topological order is connected to robust edge modes. The explicit Wiener–Hopf solution for edge-mode wavefunctions demonstrates that edge-localized zero-modes exist in number $|\omega|$ , with their spatial decay governed by analytic properties (roots or singularities) of $f(z)$ (Jones et al., 2022). In systems with longer-range hopping, edge-mode wavefunctions display algebraic, not exponential, localization, and finite-size energy splitting scales as power laws set by the order of the leading singularity of $1/f(1/z)$. Coupling decay exponent $\alpha$ determines whether the topological classification survives; for $\alpha>1$ , bulk-boundary correspondence and local invariants remain intact even in the absence of a strict gap.

For quantum Hall systems, Kitaev's distance-modulated construction extracts a purely edge-local Hamiltonian from a gapped, translation-invariant free-fermion bulk model. The resulting edge Hamiltonian hosts edge modes whose spectral flow captures the sum of bulk Chern numbers, yielding a strictly chiral branch per bulk occupied band (Du et al., 4 Dec 2025).

7. Applications: Subsystem Codes, Entanglement, and Many-Body Physics

Translation-invariant free-fermion Hamiltonians underpin exactly solvable quantum many-body and coding models. The solvability of a free-fermion spin Hamiltonian depends on the topology of its frustration graph: a translation-invariant model is solvable iff its frustration graph is a line graph $L(R)$ of some root graph $R$ . This is efficiently testable via graph-theoretic algorithms (Chapman et al., 2022). The spectral gaps of such codes are governed by the skew energy and median eigenvalue of the oriented root graph; optimal thermal robustness arises in low-dimensional, odd-coordination-number lattices.

Entanglement renormalization via explicit wavelet-based MERA circuits allows for provably accurate, analytic construction of the ground states of translation-invariant free-fermion models, both in 1D and 2D (with Fermi surfaces), maintaining rigorous control of local correlation errors (Haegeman et al., 2017). These methods extend to any translation-invariant model diagonalizable by block Fourier transform plus a half-shift, including systems with complex Fermi surfaces.

References: