Heisenberg XXZ spin-chain operator

The XXZ Heisenberg Hamiltonian on a chain of $n$ sites is

$H = J\sum_{i=1}^{n-1}\big(X_iX_{i+1} + Y_iY_{i+1} + \Delta\, Z_iZ_{i+1}\big).$

Special points. $\Delta = 1$ recovers the isotropic Heisenberg (XXX) model, exactly solvable by Bethe's 1931 ansatz; $\Delta = 0$ gives the XY model, mappable to free fermions via a Jordan–Wigner transformation; $|\Delta| \to \infty$ recovers the classical Ising limit (the XX terms become negligible relative to ZZ).

Symmetry. Writing $S^z_{\text{tot}} = \tfrac12\sum_i Z_i$, the XX+YY term can be re-expressed with raising/lowering operators $S^+_iS^-_{i+1}+S^-_iS^+_{i+1}$, which conserve the number of up-spins — a $U(1)$ symmetry, so $[H, S^z_{\text{tot}}] = 0$. This block-diagonalizes $H$ by total magnetization sector, which is why the model is tractable analytically (Bethe ansatz solves each magnetization sector separately) and why it is a natural target for excitation-number-conserving VQE ansätze.

Physical content. In the gapless regime $-1 \le \Delta \le 1$, the XXZ chain is a Luttinger liquid with power-law spin correlations; for $\Delta > 1$ the ground state is gapped and Néel (antiferromagnetically) ordered, and for $\Delta < -1$ it is ferromagnetically ordered. The Lieb–Schultz–Mattis theorem constrains this gap structure for half-integer spin chains with translational and $U(1)$ symmetry — a spin-1/2 chain at generic filling cannot have both a unique gapped ground state and these symmetries simultaneously.

Use as a testbed. Because exact Bethe-ansatz energies are known at every $\Delta$ for finite chains, the XXZ operator is a standard benchmark for VQE and quantum-simulation error analysis: any claimed ground-state energy can be checked against the exact solution rather than only against exact diagonalization of the same finite instance.