Computing Library › Quantum Simulation
Quantum Simulation

Simulating the Heisenberg Model

The quantum Heisenberg spin model, its role as a simulation benchmark, and how its dynamics are implemented on qubits.

The model

The Heisenberg Hamiltonian couples neighboring spins isotropically: H = J sum over neighboring pairs of (X_i X_j + Y_i Y_j + Z_i Z_j), optionally with anisotropy or an external field. Positive J favors antiferromagnetic order; negative J favors ferromagnetic. It captures the physics of localized magnetic moments in insulators.

Why it is a benchmark

Kronos motion — pid vs model

The Heisenberg model is exactly solvable in one dimension by the Bethe ansatz, giving a rigorous reference for validating quantum simulators. In two dimensions and with frustration it becomes hard classically, so it marks the boundary where simulation earns its keep.

Digital implementation

The pair term X_i X_j + Y_i Y_j + Z_i Z_j exponentiates to a well-known two-qubit gate. On a 1D chain the terms split into even and odd bonds, each set internally commuting, so a second-order Trotter step is a layer of even-bond gates, a layer of odd-bond gates, and back. This regular brickwork is efficient and shallow per step.

python
import numpy as np
from scipy.linalg import expm

X = np.array([[0,1],[1,0]]); Y = np.array([[0,-1j],[1j,0]]); Z = np.diag([1,-1])

def pair_term(J):
    XX = np.kron(X, X); YY = np.kron(Y, Y); ZZ = np.kron(Z, Z)
    return J * (XX + YY + ZZ)

def bond_gate(J, dt):
    return expm(-1j * pair_term(J) * dt)  # two-qubit Trotter gate

What is measured

Simulations track spin correlation functions , the spreading of entanglement after a quench, and dynamical structure factors comparable to neutron-scattering experiments. Ground-state energies come from VQE or phase estimation; dynamics come from Trotterized evolution.

Beyond one dimension

Two-dimensional and frustrated Heisenberg models suffer the classical sign problem in quantum Monte Carlo, making them genuine candidates for quantum advantage. Accurate simulation would clarify exotic phases such as spin liquids. The Heisenberg model thus serves both as a trusted calibration target and as a frontier problem, which is why it appears in nearly every quantum-simulation hardware demonstration.