Bose-Hubbard model 2026-10-25
A lattice boson Hamiltonian operator combining hopping and on-site repulsion. The on-site interaction counts distinct particle pairs through . It is the bosonic counterpart of the Hubbard model; projecting a continuum contact interaction onto localized orbitals gives when intersite overlaps are negligible.
A direct intersite Coulomb exchange integral acts without an intermediate doubly occupied state. In a repulsive one-band Hubbard model, by contrast, virtual hopping into a doubly occupied state lowers the singlet energy by relative to the triplet. The latter is antiferromagnetic kinetic exchange, even when hopping is directly between neighbouring magnetic orbitals. Ligand-mediated superexchange introduces an intermediate ligand orbital as well.
The first term assigns energy to occupying the localized orbital at site . The hopping term moves an electron between neighbouring sites, producing kinetic delocalization and band dispersion; Hermitian conjugate hops make the Hamiltonian Hermitian. The final term costs only when both spin states occupy the same local orbital. This is a Hubbard model with an on-site Coulomb penalty, not a classical interaction between magnetic dipoles.
At one electron per site and large , a real hop is energetically costly, but a virtual hop into a doubly occupied intermediate state is possible. Opposite-spin states can exploit these virtual processes and lower their energy, while parallel spins cannot occupy the same one-orbital site. The resulting singlet–triplet splitting produces antiferromagnetic Hubbard kinetic exchange.
There is a terminology qualification. The paper calls this overlap-mediated process “direct exchange” because the hopping is directly between neighbouring magnetic orbitals. In the stricter distinction between direct Coulomb exchange versus kinetic exchange, the displayed Hamiltonian omits an explicit intersite Coulomb exchange integral. The antiferromagnetic mechanism derived below is kinetic exchange; a genuine intersite direct exchange term is separate and can compete with it. This distinction also separates the model from ligand-mediated superexchange in part f.
The microscopic electric interaction is spin independent to leading order, but the Pauli exclusion principle ties the spatial symmetry of an electronic state to its spin symmetry. Different spin combinations therefore have different Coulomb and kinetic energies. Projecting those energies onto localized low-energy moments often gives an approximately rotationally invariant bilinear exchange Hamiltonian operator . For spin one half, the spin singlet state and triplet state energies are distinguished by and , so any isotropic two-level spin singlet state–triplet state splitting can be written in this form plus a constant. More general spins or orbital degeneracy can require extra interactions.
An antiferromagnetic oxide commonly has superexchange through a ligand such as oxygen. Virtual hopping between half-filled metal orbitals lowers the energy of antiparallel spins more strongly because the corresponding intermediate doubly occupied configurations are allowed by the Pauli exclusion principle. The simple one-orbital Hubbard model limit gives , equivalent to negative in the question's minus-sign convention. Bond angles and orbital occupancy can change the sign in more complicated oxides, so superexchange is not universally antiferromagnetic irrespective of geometry.
Let . On a site, empty and doubly occupied states have spin zero, while either singly occupied state has spin . Thus
The on-site spin-square normal-ordering identity rewrites the interaction as . Its contact coefficient is in this parameterization, rather than the convention often used for the Hubbard model.
Using a fermionic coherent-state trace and time slicing, replace each normal-ordered monomial by its Grassmann field symbol. With ,
The independent fields have antiperiodic thermal boundary conditions. Normal ordering is crucial: the explicit subtraction removes the spin-square contraction, so there is no extra chemical-potential shift from it.
The Hubbard model combines band motion and local repulsion:
Hopping favors itinerant states, while positive penalizes simultaneous opposite-spin occupation on one site. In a uniform mean-field treatment, a spin imbalance reduces the interaction cost but raises the band-filling energy. The Stoner criterion for instability of the unpolarized metal is when is the density of states per site, per energy, for one spin. Using a total two-spin density of states changes the displayed numerical convention.
The ordered mean-field state has exchange-split spin bands and a spontaneous magnetization, with the kinetic and interaction costs balancing at equilibrium. This is a model of itinerant ferromagnetism, not a guarantee that every repulsive Hubbard system becomes ferromagnetic. Filling, lattice structure and competing order matter: near half filling on a bipartite lattice, antiferromagnetism and a Mott insulator can be favored instead. The approximation also omits important spin fluctuations and strong-correlation effects.
The Hubbard model combines hopping with a local repulsion,
Hopping favours a broad itinerant band; the interaction penalizes opposite-spin double occupancy. For uniform Stoner decoupling of the Hubbard interaction, let be the spin imbalance per site. The interaction is , so polarization lowers it by , at the expense of kinetic energy.
If is the noninteracting density of states per spin per site, small polarization costs . Therefore
is the mean-field ferromagnetic instability criterion. The exchange splitting is , and the paramagnetic response is enhanced as . A density of states including both spins changes the written factor in the criterion.
The model explains how itinerant bands can polarize without imposing pre-existing localized atomic spins. It is not a claim that every large- metal becomes ferromagnetic: near half filling the Hubbard model can instead form a Mott insulator with antiferromagnetic superexchange of scale . Band filling, bandwidth and fluctuations matter.[ref] The interaction-energy and kinetic-cost expansion above is the stated mean-field derivation of the enhancement.[ref]
Specify the spin counting first. Let denote the unpolarized single-spin density of states per atom per energy, and put . At fixed total occupation , the two Hubbard model mean-field energies can be written
The common shift can be absorbed into the chemical potential. Lowering the up-spin band by adds carriers, while raising the down-spin band removes the same number. Thus
If is the atom number density, the magnetization is . The Stoner-enhanced magnetic susceptibility in SI volume units is therefore
Here is the noninteracting Pauli paramagnetism result. If the paper's counts both spins, and the same result is
A susceptibility per atom omits . Stating the spin and volume conventions avoids an otherwise spurious factor of two or a dimensional mismatch.
The Stoner criterion is , equivalently with spin-summed . To see the instability, transfer electrons per atom from one spin band to the other. The band-filling cost is , whereas the interaction energy changes by . Their sum becomes negative at the threshold. The unpolarized state then ceases to be a local minimum, and nonlinear band filling determines the spontaneous ferromagnetism. The negative continuation of the displayed linear susceptibility beyond the threshold is not a stable paramagnetic response.
Use and dimensionless spin operators. At one site, the canonical anticommutation relations imply the on-site spin-square normal-ordering identity
Indeed the empty and doubly occupied states have spin zero, while either singly occupied state has . Thus the interaction in normal-ordered form is . There is no residual chemical-potential shift: that contraction has cancelled.
Insert fermionic coherent state resolutions between imaginary-time evolution factors. With and , the grand-canonical coherent-state action is
The independent barred and unbarred Grassmann variables satisfy antiperiodic boundary conditions at . Adjacent time-slice labels in the normal-ordered symbol fix the continuum prescription and thermal normalization. The interaction coefficient is , not the of a differently normalized Hubbard model.
The t-J model approximates the low-energy sector of a repulsive Hubbard model with near one electron per site. Large on-site repulsion suppresses double occupancy. Virtual hopping out of and back into that low-energy sector produces antiferromagnetic superexchange with . One conventional effective Hamiltonian is
with projected electron operator . At half filling, nearest-neighbor hopping cannot move an electron without entering a forbidden doubly occupied state, leaving the leading spin-exchange model. Doped holes can move and couple the charge dynamics to these spin correlations. Additional terms, including three-site processes at order , are omitted in this simplified Hamiltonian.
The exchange term energetically favors a nearest-neighbor spin singlet. With in the no-double-occupancy subspace, . A square-lattice pairing ansatz can therefore use bond singlets with opposite signs along the two axes. This d-wave bond order has
It changes sign under a quarter-turn and vanishes along the zone diagonals. It is an even-parity, spin-singlet pairing pattern, not a statement of triplet pairing or necessarily broken time reversal. Phase-sensitive interference experiments distinguish this sign structure from a merely anisotropic nonnegative gap in cuprate superconductors.[ref]
The model captures strong local repulsion, antiferromagnetism and mobile doped carriers. Its interacting phase diagram is not solved by writing down the Hamiltonian: a mean-field d-wave ansatz is a proposed approximation, not a proof that every parameter choice has that exact ground state or a complete microscopic explanation of every high-temperature superconductor.
Let be the single-spin density of states per site per energy, and let the on-site Hubbard model repulsion have energy . For a weak uniform magnetic field , the half-splitting of the two spin bands is , where . Constant density of states gives . Therefore
This is the response per site; multiplying by the site number density gives the dimensionless SI volume magnetic susceptibility. With spin-summed density of states , the denominator is . Its vanishing is the mean-field Stoner criterion. The unpolarized linear-response expression is stable only for .
t-J model 2026-10-20
An effective low-energy lattice model obtained from a strongly repulsive Hubbard model near one electron per site. The retained Hilbert space excludes double occupancy. Its usual nearest-neighbor Hamiltonian has projected hopping and antiferromagnetic superexchange . Virtual intermediate states with one double occupancy cost energy , giving to leading order. Three-site hopping terms of the same expansion order are omitted in this simplified version.
At half filling the projected hopping is blocked. Doped holes restore charge motion coupled to spin correlations. Nearest-neighbor singlet pairing can support a d-wave bond order ansatz on a square lattice, but the Hamiltonian alone is not a proof of a unique superconducting ground state for all parameters.