The Hubbard model and Mott localization
Separate interaction-driven charge localization from a filled-band insulator and derive superexchange in a controlled limit.
On this page
The question for today
What you’ll understand
- Explain why one electron per site need not produce a metal.
- Identify the small parameter behind a spin-only description.
- Calculate the leading antiferromagnetic exchange scale.
A partially filled noninteracting band suggests mobile charge, but electrons also repel one another. The Hubbard model isolates this competition in an intentionally small Hamiltonian: particles hop between sites and pay an energy when opposite spins share a site. Its simplicity makes every assumption inspectable. It does not make the general many-electron problem easy. We will use an exactly understood atomic limit and a controlled expansion, rather than infer a complete phase diagram from a sketch.
A minimal competition with a large Hilbert space
Take one orbital per site, two spin states, nearest-neighbor hopping , and repulsion . The local occupations are empty, up, down, and doubly occupied. Hopping favors states spread over the lattice, whereas the interaction penalizes double occupancy. The model omits other orbitals, long-range Coulomb terms, lattice vibrations, and disorder. Those omissions can matter for a material, so a successful qualitative mechanism is not automatically a quantitative fit.
Half filling means one electron per site on average, not a prohibition against empty or doubly occupied sites. At , the half-filled band can conduct. At , moving charge between two singly occupied sites creates one empty site and one double occupation at cost . Spin orientations remain degenerate in this atomic limit. Charge localization therefore does not require the noninteracting band to be full, and it does not by itself choose a magnetic pattern.
Virtual hopping leaves a magnetic interaction
At half filling with , low-energy configurations have almost one particle per site. A single hop leaves that subspace, but a second hop can return to it. The intermediate double occupation is virtual: it contributes to an amplitude without becoming a long-lived low-energy charge excitation. Second-order perturbation theory produces an exchange proportional to . This expansion requires the charge-excitation denominator to remain large compared with the hopping.
For the two-site problem, the spin singlet can lower its energy by admixing doubly occupied configurations; the triplet does not gain the same energy. The resulting exchange is antiferromagnetic. After an additive constant, the effective Hamiltonian is a Heisenberg spin model. Its spin dynamics occur on a scale much smaller than the charge penalty. Deriving that scale is controlled at small ; solving a frustrated spin model at arbitrary size can still be difficult.
Localization, symmetry, and unresolved regimes
A band insulator can arise without interactions when all available low-energy bands are filled. A Mott insulator instead involves suppressed charge motion caused by interactions in a setting where band counting alone would allow a metal. In practice, magnetic order can also reconstruct bands. Separating these mechanisms requires attention to symmetry, temperature, spectral gaps, and how the state evolves as parameters change. The label should summarize evidence rather than replace it.
Away from half filling, holes can move without always creating the same expensive double occupation, and the spin-only description no longer contains all low-energy degrees of freedom. Intermediate coupling in two dimensions remains a demanding many-body problem. Mean-field, finite-cluster, Monte Carlo, and tensor-network calculations have distinct errors and accessible regimes. Agreement in a controlled limit is a useful calibration; it cannot settle every parameter region that the approximation does not control.
Work it through
Two separated energy scales
In a two-site half-filled model take t = 0.10 eV and U = 4.0 eV. Estimate the leading singlet-triplet splitting and assess the expansion parameter.
Compute . Double-occupation admixture is perturbative in this small ratio.
Use .
For two spins, is in the singlet and in a triplet. The shifted effective Hamiltonian gives singlet energy and triplet energy zero.
The leading splitting is 10 meV, far below the 4 eV atomic charge penalty. It is an asymptotic estimate for this Hamiltonian, not a measured material gap.
The living notebook
Explore the model
Interactive lab · optimization
Distinguish a controlled regime from a diverging one
Adjust step size η in gradient descent on x²/2 and inspect xₙ = (1−η)ⁿx₀. This separate mathematical example makes the domain of an iterative approximation visible. It does not diagonalize the Hubbard model, calculate double occupancy, or identify a Mott transition.
- Position x0
- 2
- Objective x²/2
- 2
- Iteration behavior
- Stable step
Gradient descent on f(x) = x²/2: xₙ = (1 − η)ⁿx₀. Playback advances two discrete iterations per second. Values outside |x| ≤ 4 are explicitly reported. 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 forEvery approximation needs a validity condition. Here stability is explicit; in superexchange the relevant small parameter is t/U.
Pause & reason
Check your understanding
Your turn
Try it yourself
Keeping U fixed, halve t in the worked example. Find the new leading exchange and explain why the spin-only approximation improves.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Count powers of t separately in J and in the expansion parameter.
Compare with the solution
With t = 0.050 eV, J = 4 × 0.050²/4.0 eV = 0.0025 eV = 2.5 meV, one quarter of the original exchange. Meanwhile t/U falls from 0.025 to 0.0125, suppressing omitted virtual-charge corrections. Lower spin energy does not mean the atomic charge penalty has fallen.
Keep the idea
Take it with you
- Interactions can localize charge in a partially filled band.
- Virtual charge motion produces spin exchange even when real charge motion is costly.
- A controlled strong-coupling expansion does not solve the entire Hubbard phase diagram.
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.