Superexchange from a strong-coupling expansion
Derive antiferromagnetic exchange from virtual double occupancy and benchmark its coefficient against an exact dimer.
On this page
The question for today
What you’ll understand
- Project a half-filled repulsive Hubbard model into its singly occupied sector.
- Derive the sign and coefficient of the leading exchange.
- State the limits of a spin-only effective Hamiltonian.
Suppressing real charge motion does not make an interaction-driven insulator magnetically inert. A pair of electrons can virtually visit a doubly occupied configuration and return, lowering some spin states more than others. This superexchange process survives at energies far below the atomic charge gap. Its derivation is a controlled strong-coupling expansion for a specified one-orbital repulsive model; its positive exchange is not a universal sign rule for every magnetic material.
Separate low and high charge sectors
At exactly one electron per site and U much larger than the hopping magnitude, define a projector onto the subspace with no empty or doubly occupied sites. A single hop leaves that subspace by creating a neighboring hole and doublon. Consequently the projected hopping vanishes to first order. The intermediate configuration costs U in the simple uniform model, so a second hop can return to the low-energy subspace with an energy shift of order hopping squared over U.
Eliminating the higher-energy sector gives a negative second-order operator because its energy denominators are positive relative to the low-energy reference. That negative sign represents energy lowering through virtual admixture. The operators still retain fermionic signs and spin exchange; treating the process as classical random jumps loses the relative signs that distinguish the singlet from the triplet. The expansion applies to low energies and temperatures compared with the charge gap, with corrections organized in powers of hopping over interaction.
Read the coefficient from spin symmetry
For two sites and two electrons, the spin triplet cannot couple to an on-site pair because one orbital with two electrons supports only a spin singlet. The low-energy singlet does mix with doubly occupied states and gains energy. Rotational invariance therefore fixes the projected bond interaction to an antiferromagnetic Heisenberg form plus a constant. Use dimensionless spin operators with eigenvalues one half, and count each undirected bond once; these choices fix the numerical coefficient.
The resulting exchange is four times squared hopping magnitude over U. Keeping the density offset makes the triplet bond energy zero and the singlet bond energy minus the exchange, directly matching the Hubbard reference energies. Dropping that constant changes absolute energies but not their separation. For complex spin-independent hopping the leading magnitude remains squared modulus, while spin-orbit-dependent hopping can produce anisotropic exchange. Multiple orbitals, Hund coupling and direct exchange can change the sign or structure and are absent here.
Benchmark the approximation before using a lattice
The Hubbard dimer can be diagonalized exactly in its singlet sector. Its lower singlet energy contains a square root that expands to the second-order result plus higher even powers of hopping. This exact benchmark reveals the leading error without invoking a thermodynamic phase transition. Agreement for a small hopping ratio establishes the effective bond scale, but does not solve a frustrated quantum spin lattice or establish its ordered state.
At finite doping the low-energy subspace also contains mobile holes, so a spin-only Hamiltonian is insufficient and projected hopping must remain. Even at half filling, longer virtual loops generate additional interactions at higher order. Exchange sets a useful temperature scale after division by Boltzmann constant, but an ordering temperature also depends on dimension, anisotropy and fluctuations. The two-level lab can illustrate singlet-state superposition algebra only after a basis is declared; it cannot perform the strong-coupling projection or represent an interacting electron lattice.
Work it through
Check the exchange against the dimer
A Hubbard dimer has and . Compute the leading exchange, approximate singlet energy and exact singlet energy. Estimate the relative error in the predicted splitting.
The small ratio is . The exchange is , so the projected triplet is at zero and the approximate singlet is at .
The exact energy is . It is less negative than the second-order estimate, consistent with the positive fourth-order correction.
The exact splitting is . The excess predicted splitting is approximately , or about of the exact value. This checks a bond-scale expansion, not a bulk ordering temperature.
The leading exchange is 20 meV; the exact singlet-triplet splitting is 19.8 meV, giving about a 1% second-order error.
The living notebook
Explore the model
Interactive lab · quantum
Two spin configurations as a limited basis analogy
Reset and set BOTH Mixing angle θ = 90° and Relative phase φ = 0°, then Reset and set θ = 90° and φ = 180°. Interpret the basis states as |↑↓⟩ and |↓↑⟩ only for comparing symmetric and antisymmetric superpositions; predict equal basis probabilities in both. This pure two-state lab does not include double occupancy, compute exchange energies or solve a many-electron Hubbard model.
- Probability |0⟩
- 71.1%
- Probability |1⟩
- 28.9%
- Bloch (x, y, z)
- 0.78, 0.45, 0.42
A normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forRelative phase distinguishes spin superpositions even when their basis probabilities agree.
Pause & reason
Check your understanding
Your turn
Try it yourself
Take and in the same half-filled strong-coupling convention. Compute J and the scale J/kB using . Does the result determine a Néel temperature?
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The exchange is an energy. Converting it to kelvin gives a microscopic scale, not an automatic transition temperature.
Compare with the solution
The exchange is . Its temperature scale is . This does not determine a Néel temperature: ordering depends on lattice dimension, frustration, anisotropy and fluctuations, and an ideal two-dimensional isotropic short-range Heisenberg model does not develop conventional finite-temperature continuous-spin order. The calculation only identifies the leading low-energy exchange scale under the stipulated single-orbital assumptions.
Keep the idea
Take it with you
- Virtual charge excitations produce low-energy spin interactions at second order.
- The positive coefficient 4|t|²/U requires the stated orbital, filling and bond conventions.
- An exchange scale does not solve the spin lattice or determine its ordering temperature.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.