research29 min studyLesson 66 in Physics

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

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 t>0t>0, and repulsion U>0U>0. 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 U=0U=0, the half-filled band can conduct. At t=0t=0, moving charge between two singly occupied sites creates one empty site and one double occupation at cost UU. 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.

H=−t∑⟨ij⟩,σ(ciσ†cjσ+cjσ†ciσ)+U∑ini↑ni↓H=-t\sum_{\langle ij\rangle,\sigma}(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma})+U\sum_i n_{i\uparrow}n_{i\downarrow}
Each unordered neighboring pair appears once. This one-orbital Hubbard Hamiltonian conserves total particle number; its parameters carry units of energy.↔ Scroll wide equations horizontally.

Virtual hopping leaves a magnetic interaction

At half filling with U≫tU\gg t, 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 t2/Ut^2/U. 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 t/Ut/U; solving a frustrated spin model at arbitrary size can still be difficult.

Heff(2)=J∑⟨ij⟩(Si⋅Sj−14),J=4t2UH_{\mathrm{eff}}^{(2)}=J\sum_{\langle ij\rangle}\left(\mathbf S_i\cdot\mathbf S_j-\frac14\right),\qquad J=\frac{4t^2}{U}
This is the second-order result. In the two-site or nearest-neighbor bipartite setting, omitted terms begin at order t⁴/U³ and can include additional interactions, not only a correction to J. Spins are dimensionless with Planck's constant set to one.↔ Scroll wide equations horizontally.

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.

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.

  1. Compute t/U=0.10/4.0=0.025t/U=0.10/4.0=0.025. Double-occupation admixture is perturbative in this small ratio.

  2. Use J=4(0.10)2/4.0 eV=0.010 eV=10 meVJ=4(0.10)^2/4.0\,\mathrm{eV}=0.010\,\mathrm{eV}=10\,\mathrm{meV}.

  3. For two spins, S1⋅S2\mathbf S_1\cdot\mathbf S_2 is −3/4-3/4 in the singlet and 1/41/4 in a triplet. The shifted effective Hamiltonian gives singlet energy −J-J and triplet energy zero.

The result

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.

Explore the model

Interactive lab

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.

Distinguish a controlled regime from a diverging oneGradient descent on f(x) = x²/2: xₙ = (1 − η)ⁿx₀. Playback advances two discrete iterations per second. Values outside |x| ≤ 4 are explicitly reported. Current readings: Position x0: 2; Objective x²/2: 2; Iteration behavior: Stable step. All parameters are adjustable with the labeled controls below.−4−202402.174.356.528.7Parameter xf(x) = x²/2n = 0
ObjectiveDiscrete iterate
Position x0
2
Objective x²/2
2
Iteration behavior
Stable step
n = 0 · t = 0 s

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.

Check your understanding

At half filling with U much larger than t, which low-energy scale does virtual hopping generate?

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.

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.