Skip to content

Common Many-Body Hamiltonians

A many-body Hamiltonian is not fully specified by one line of algebra. The same formula can describe different physics when the Hilbert space, statistics, lattice, boundary conditions, filling, sign conventions, cutoff, or ensemble changes.

If the unresolved question is which family can predict a material observable, start with Choosing a Model for Quantum Matter. This page begins after that choice and supplies convention-aware formulas and canonical model links.

This page collects eleven standard Hamiltonian families:

  1. continuum bosons with contact interactions;
  2. continuum fermions with two-body interactions;
  3. tight-binding particles;
  4. the spinless tt–VV chain;
  5. the fermionic Hubbard model;
  6. the Bose–Hubbard model;
  7. the Heisenberg model;
  8. the transverse-field Ising model;
  9. the reduced BCS pairing model;
  10. the Anderson impurity model;
  11. the Kondo impurity model.

It is a convention-aware lookup sheet. The Model Encyclopedia owns convention-complete dossiers, comparisons, solution-status labels, and benchmark handoffs. The linked teaching pages own derivations, phase structure, methods, and applications. The Model-to-Volume Cross-Link Index records which page owns each layer and where planned applications will live.

Before using any row, state:

  • the Hilbert space and local degrees of freedom;
  • bosonic, fermionic, spin, or constrained statistics;
  • continuum or lattice geometry and spatial dimension;
  • open, periodic, twisted, or other boundary conditions;
  • fixed particle number or grand-canonical controls;
  • which bonds or ordered pairs are included in each sum;
  • whether spin operators or Pauli matrices are used;
  • the ultraviolet regulator or lattice spacing;
  • the approximation and energy window in which the model is valid.

The notation

K=H−μNK = H-\mu N

denotes a grand Hamiltonian. Many papers call KK simply “the Hamiltonian.” This sheet keeps HH and KK distinct whenever the distinction prevents a convention error.

ModelDegrees of freedomCompeting scalesCanonical continuation
continuum contact bosonsbosonic field ψ(x)\psi(\mathbf x)kinetic energy, gg, densityField Operators in Many-Body Models
continuum interacting fermionsspinful fermionic fieldskinetic energy, interaction, densityField Operators in Many-Body Models
tight bindingparticles in localized orbitalshopping, onsite energiesTight-Binding Chain
spinless tt–VV chainone spinless fermion mode per sitett, VV, fillingSpinless Fermion Chains
Hubbardspin-1/21/2 lattice fermionstt, UU, fillingHubbard Model
Bose–Hubbardlattice bosonstt, UU, μ\muBose–Hubbard Model
Heisenberglocalized spinsexchange JijJ_{ij}, fieldsHeisenberg Model
transverse-field Isingspin-1/21/2 sitesIsing exchange, transverse fieldTransverse-Field Ising Model
reduced BCStime-reversed fermion pairsattraction, cutoff, density of statesExact dossier; BCS Mean-Field Theory
Anderson impurityone correlated orbital and bath fermionsϵd\epsilon_d, UU, hybridization functionmodel dossier; teaching preview
Kondoimpurity spin and conduction fermionsJKJ_K, bandwidth, density of statesmodel dossier; screening and RG

Use bosonic Fock space over a continuum one-particle space, with equal-time algebra

[ψ(x),ψ†(y)]=δ(d)(x−y).[\psi(\mathbf x),\psi^\dagger(\mathbf y)] = \delta^{(d)}(\mathbf x-\mathbf y).

The density and total number are

n(x)=ψ†(x)ψ(x),N=∫ddx n(x).n(\mathbf x) = \psi^\dagger(\mathbf x)\psi(\mathbf x), \qquad N = \int d^d x\, n(\mathbf x).

A common effective Hamiltonian is

HB=∫ddx [ψ†(−ℏ2∇22m+Vext)ψ+g2ψ†ψ†ψψ](x).\begin{aligned} H_B = \int d^d x\, \bigg[ & \psi^\dagger \left( -\frac{\hbar^2\nabla^2}{2m} + V_{\mathrm{ext}} \right) \psi \\ & + \frac{g}{2} \psi^\dagger \psi^\dagger \psi \psi \bigg](\mathbf x). \end{aligned}

The corresponding grand Hamiltonian is

KB=HB−μN.K_B = H_B-\mu N.

The factor 1/21/2 avoids double counting pairs. The displayed contact term is normal ordered.

  • g>0g>0 denotes repulsion in the displayed convention.
  • Coincident-point products require a regulator and matching prescription.
  • In a dilute three-dimensional pseudopotential regime,
g=4πℏ2asm,g = \frac{4\pi\hbar^2a_s}{m},

where asa_s is the ss-wave scattering length.

  • This relation is not a cutoff-independent microscopic identity at arbitrary density, range, or dimension.
  • The model conserves NN when no source or conversion term is added.
LimitResult
g=0g=0ideal continuum Bose gas
weak depletion and low temperatureBogoliubov expansion around a condensate
strong short-range repulsion in one dimensionLieb–Liniger and Tonks–Girardeau regimes, with convention changes
periodic potential and localized Wannier basisBose–Hubbard model after projection

For the detailed operator form, momentum-space normalization, and regularization cautions, use Field Operators in Many-Body Models.

Continuum Fermions with Two-Body Interaction

Section titled “Continuum Fermions with Two-Body Interaction”

For components σ\sigma,

{ψσ(x),ψσ′†(y)}=δσσ′δ(d)(x−y).\{ \psi_\sigma(\mathbf x), \psi_{\sigma'}^\dagger(\mathbf y) \} = \delta_{\sigma\sigma'} \delta^{(d)}(\mathbf x-\mathbf y).

The total number is

N=∑σ∫ddx ψσ†(x)ψσ(x).N = \sum_\sigma \int d^d x\, \psi_\sigma^\dagger(\mathbf x) \psi_\sigma(\mathbf x). HF=∑σ∫ddx ψσ†(x)h0,σψσ(x)+12∑σ,σ′∫ddx ddy ψσ†(x)ψσ′†(y)×Vσσ′(x,y)ψσ′(y)ψσ(x).\begin{aligned} H_F ={}& \sum_\sigma \int d^d x\, \psi_\sigma^\dagger(\mathbf x) h_{0,\sigma} \psi_\sigma(\mathbf x) \\ &+ \frac12 \sum_{\sigma,\sigma'} \int d^d x\,d^d y\, \psi_\sigma^\dagger(\mathbf x) \psi_{\sigma'}^\dagger(\mathbf y) \\ &\qquad\qquad{}\times V_{\sigma\sigma'}(\mathbf x,\mathbf y) \psi_{\sigma'}(\mathbf y) \psi_\sigma(\mathbf x). \end{aligned}

A standard one-body operator is

h0,σ=−ℏ2∇22mσ+Vext,σ(x).h_{0,\sigma} = -\frac{\hbar^2\nabla^2}{2m_\sigma} + V_{\mathrm{ext},\sigma}(\mathbf x).

The grand Hamiltonian is

KF=HF−∑σμσNσ.K_F = H_F - \sum_\sigma \mu_\sigma N_\sigma.
  • Hermiticity requires the interaction kernel to obey the appropriate exchange-conjugation relation.
  • The factor 1/21/2 assumes the integration counts both ordered particle pairs.
  • For identical single-component fermions, a zero-range ss-wave contact term vanishes because two fields at the same point anticommute.
  • A two-component contact interaction can be written
Hint=g∫ddx n↑(x)n↓(x),H_{\mathrm{int}} = g \int d^d x\, n_\uparrow(\mathbf x) n_\downarrow(\mathbf x),

with the sign of gg setting attraction or repulsion in the displayed convention.

  • Coulomb, dipolar, finite-range, and effective interactions require their own kernels and regularization.
  • Random Phase Approximation owns the self-consistent density response and screened interaction generated from a specified direct kernel.
LimitResult
V=0V=0ideal Fermi gas
weak short-range attractionCooper instability and BCS reduction near the Fermi surface
periodic potential and one-band projectiontight-binding or Hubbard descriptions
low-energy quasiparticle regimeLandau Fermi-liquid description when applicable

The derivation, finite boundaries, graph form, and band connection are developed in Tight-Binding Model. This section is the compact convention reference.

Choose localized orbitals ∣i,σ⟩\lvert i,\sigma\rangle. The same one-body form can act on a one-particle Hilbert space or on a bosonic or fermionic Fock space. Filling and statistics must be specified before assigning many-body meaning.

Htb=−∑i,j,σtijciσ†cjσ+∑i,σviσniσ.H_{\mathrm{tb}} = - \sum_{i,j,\sigma} t_{ij} c_{i\sigma}^\dagger c_{j\sigma} + \sum_{i,\sigma} v_{i\sigma} n_{i\sigma}.

Hermiticity requires

tji=tij∗,viσ∈R.t_{ji} = t_{ij}^*, \qquad v_{i\sigma}\in\mathbb R.

For real nearest-neighbor hopping with each unordered bond counted once,

Htb=−t∑⟨i,j⟩,σ(ciσ†cjσ+cjσ†ciσ)+ϵ0∑i,σniσ.H_{\mathrm{tb}} = -t \sum_{\langle i,j\rangle,\sigma} \left( c_{i\sigma}^\dagger c_{j\sigma} + c_{j\sigma}^\dagger c_{i\sigma} \right) + \epsilon_0 \sum_{i,\sigma} n_{i\sigma}.

For a periodic chain with lattice spacing aa,

ϵ(k)=ϵ0−2tcos⁡(ka).\epsilon(k) = \epsilon_0 - 2t\cos(ka).

The sign of tt fixes where the band minimum occurs in this convention. It can be changed by a sublattice gauge transformation only for selected bipartite hopping networks.

The tight-binding term alone contains no interaction. Calling a partially filled tight-binding band a metal, insulator, magnet, or superconductor requires filling, dimensionality, symmetry, interactions, disorder, and boundary conditions.

Retain one spinless fermionic mode per site,

{ci,cj†}=δij,nj=cj†cj∈{0,1}.\{c_i,c_j^\dagger\} = \delta_{ij}, \qquad n_j = c_j^\dagger c_j \in \{0,1\}.

For LL sites, the fixed-NN sector has dimension (LN)\binom{L}{N}.

For an open chain,

KtV=−t∑j=1L−1(cj†cj+1+cj+1†cj)+V∑j=1L−1(nj−12)(nj+1−12)−μ∑j=1L(nj−12).\begin{aligned} K_{tV} = {}& -t \sum_{j=1}^{L-1} \left( c_j^\dagger c_{j+1} + c_{j+1}^\dagger c_j \right) \\ &+ V \sum_{j=1}^{L-1} \left(n_j-\frac12\right) \left(n_{j+1}-\frac12\right) \\ &- \mu \sum_{j=1}^{L} \left(n_j-\frac12\right). \end{aligned}

The centered convention makes particle-hole symmetry occur at μ=0\mu=0 on an even bipartite chain. Replacing the interaction by V∑jnjnj+1V\sum_jn_jn_{j+1} requires a corresponding chemical-potential and constant shift.

For a ring, include the LLth bond and specify

cL+1=eiϕc1.c_{L+1} = e^{i\phi}c_1.

The twist ϕ\phi, the filling N/LN/L, and whether the model is defined directly or obtained through the Jordan–Wigner Transformation are part of the Hamiltonian data.

At V=0V=0,

ϵ(k)=−2tcos⁡(ka)−μ.\epsilon(k) = -2t\cos(ka)-\mu.

At half filling for t>0t>0, the uniform thermodynamic chain is phase separated for V<−2tV<-2t, is a gapless Luttinger liquid for −2t<V<2t-2t<V<2t, and has gapped period-two charge order for V>2tV>2t. The endpoint at V=2tV=2t is a Berezinskii–Kosterlitz–Thouless transition.

After a staggered hopping-sign transformation, the XXZ parameters are

J=2t,Δ=V2t.J = 2t, \qquad \Delta = \frac{V}{2t}.

Spinless Fermion Chains owns the free solution, correlators, phase interpretation, transport twist, and finite-size parity caveats. XXZ Spin Chain owns the canonical spin-chain treatment.

Each site has four local states,

∣0⟩,∣↑⟩,∣↓⟩,∣↑↓⟩.\lvert0\rangle, \qquad \lvert\uparrow\rangle, \qquad \lvert\downarrow\rangle, \qquad \lvert\uparrow\downarrow\rangle.

For LL sites, the unrestricted Fock space has dimension 4L4^L.

HHub=−t∑⟨i,j⟩,σ(ciσ†cjσ+h.c.)+U∑ini↑ni↓+∑ivini.\begin{aligned} H_{\mathrm{Hub}} = & -t \sum_{\langle i,j\rangle,\sigma} \left( c_{i\sigma}^\dagger c_{j\sigma} + \mathrm{h.c.} \right) \\ & + U \sum_i n_{i\uparrow}n_{i\downarrow} + \sum_i v_i n_i. \end{aligned}

Here

ni=ni↑+ni↓.n_i = n_{i\uparrow} + n_{i\downarrow}.

The grand Hamiltonian is

KHub=HHub−μ∑ini.K_{\mathrm{Hub}} = H_{\mathrm{Hub}} - \mu \sum_i n_i.
  • t>0t>0 in the displayed nearest-neighbor convention.
  • U>0U>0 is repulsive and U<0U<0 attractive.
  • The filling is
n=NL.n = \frac{N}{L}.
  • Half filling is n=1n=1, because each site has two spin modes.
  • Each unordered bond appears once when the Hermitian conjugate is written explicitly.
LimitResult
U=0U=0spinful tight-binding fermions
t=0t=0independent atomic sites
repulsive U≫tU\gg t near half fillinglow-energy spin model with superexchange
U<0U<0onsite pairing tendency

At half filling in the large-repulsion regime, the leading exchange scale is

Jex=4t2UJ_{\mathrm{ex}} = \frac{4t^2}{U}

for the standard two-site hopping convention. Higher-order and lattice-dependent terms are omitted.

The full conventions, symmetries, limits, and strong-coupling derivation live in Hubbard Model.

Each lattice site is a bosonic mode with

ni=bi†bi∈{0,1,2,…}.n_i = b_i^\dagger b_i \in \{0,1,2,\ldots\}.

Numerical calculations often impose a maximum occupation as a truncation. That cutoff is not part of the exact model.

KBH=−t∑⟨i,j⟩(bi†bj+bj†bi)+U2∑ini(ni−1)−μ∑ini.\begin{aligned} K_{\mathrm{BH}} = & -t \sum_{\langle i,j\rangle} \left( b_i^\dagger b_j + b_j^\dagger b_i \right) \\ & + \frac{U}{2} \sum_i n_i(n_i-1) - \mu \sum_i n_i. \end{aligned}

For a fixed-NN calculation, omit the chemical-potential term and restrict the Hilbert space to

∑ini=N.\sum_i n_i = N.
  • U>0U>0 is onsite repulsion.
  • The factor ni(ni−1)/2n_i(n_i-1)/2 counts unordered onsite pairs.
  • t=0t=0 is the atomic limit.
  • The competition between t/Ut/U and μ/U\mu/U organizes the standard Mott-lobe diagram.
  • At integer filling, large U/tU/t suppresses number fluctuations; large t/Ut/U favors phase coherence.

The full conventions, limits, optical-lattice reduction, and mean-field boundary live in Bose–Hubbard Model. Bose–Hubbard Dimer owns the exact two-mode fixed-number specialization, while Bose–Hubbard Chain owns the one-dimensional soft-core model and its hard-core and Luttinger limits. The compact model card remains a lookup entry, while Quantum Phase Transitions owns the generic critical interpretation.

Place a spin-ss representation on each site:

H=⨂i=1LC2s+1.\mathcal H = \bigotimes_{i=1}^{L} \mathbb C^{2s+1}.

Let si\mathbf s_i be dimensionless spin operators satisfying

[sia,sjb]=iδijϵabcsic.[s_i^a,s_j^b] = i\delta_{ij} \epsilon^{abc} s_i^c. HHeis=∑⟨i,j⟩Jijsi⋅sj−∑ihi⋅si.H_{\mathrm{Heis}} = \sum_{\langle i,j\rangle} J_{ij} \mathbf s_i\cdot\mathbf s_j - \sum_i \mathbf h_i\cdot\mathbf s_i.

Each unordered bond is counted once. The field hi\mathbf h_i is measured in energy units.

For uniform nearest-neighbor exchange,

HHeis=J∑⟨i,j⟩si⋅sj.H_{\mathrm{Heis}} = J \sum_{\langle i,j\rangle} \mathbf s_i\cdot\mathbf s_j.

With the displayed plus sign:

  • J>0J>0 is antiferromagnetic;
  • J<0J<0 is ferromagnetic.

For two spin-1/21/2 sites,

s1⋅s2={−3/4,singlet,+1/4,triplet.\mathbf s_1\cdot\mathbf s_2 = \begin{cases} -3/4, & \text{singlet},\\ +1/4, & \text{triplet}. \end{cases} HXYZ=∑⟨i,j⟩(Jxsixsjx+Jysiysjy+Jzsizsjz)H_{XYZ} = \sum_{\langle i,j\rangle} \left( J_xs_i^xs_j^x + J_ys_i^ys_j^y + J_zs_i^zs_j^z \right)

is anisotropic and generally has less symmetry than the isotropic model. Dzyaloshinskii–Moriya exchange, longer-range bonds, single-ion anisotropy, and multi-spin interactions are distinct extensions.

For the uniform spin-1/21/2 chain with Jx=Jy=JJ_x=J_y=J and Jz=JΔJ_z=J\Delta, use XXZ Spin Chain for phase structure, Jordan–Wigner fermions, and integrability.

Use Heisenberg Model for exchange spectra, symmetries, limiting cases, and approximation methods.

For spin-1/21/2 sites with Pauli matrices,

HTFIM=−J∑⟨i,j⟩σizσjz−h∑iσix.H_{\mathrm{TFIM}} = -J \sum_{\langle i,j\rangle} \sigma_i^z\sigma_j^z - h \sum_i \sigma_i^x.

The Ising interaction and transverse-field terms do not commute.

The global spin flip

P=∏iσixP = \prod_i \sigma_i^x

satisfies

[HTFIM,P]=0.[H_{\mathrm{TFIM}},P] = 0.

It changes

σiz⟼−σiz.\sigma_i^z \longmapsto -\sigma_i^z.

If

siα=σiα2,s_i^\alpha = \frac{\sigma_i^\alpha}{2},

then

σizσjz=4sizsjz,σix=2six.\sigma_i^z\sigma_j^z = 4s_i^zs_j^z, \qquad \sigma_i^x = 2s_i^x.

Couplings must therefore be rescaled when changing notation.

For the infinite nearest-neighbor chain in the displayed convention, the critical point is

∣h∣=∣J∣.|h| = |J|.

Boundary conditions and parity sectors matter in the finite Jordan–Wigner solution. Use Ising Chain for the compact model card, Transverse-Field Ising Model for the model treatment and exact-solution preview, and Quantum Phase Transitions for general critical scaling.

Use spinful fermionic momentum modes. The reduced interaction scatters pairs of time-reversed states,

(k↑,−k↓)⟶(k′↑,−k′↓).(\mathbf k\uparrow,-\mathbf k\downarrow) \longrightarrow (\mathbf k'\uparrow,-\mathbf k'\downarrow).

Define

bk†=ck↑†c−k↓†.b_{\mathbf k}^\dagger = c_{\mathbf k\uparrow}^\dagger c_{-\mathbf k\downarrow}^\dagger. KBCS=∑k,σξkckσ†ckσ−V∑k,k′∈Sbk†bk′.\begin{aligned} K_{\mathrm{BCS}} = & \sum_{\mathbf k,\sigma} \xi_{\mathbf k} c_{\mathbf k\sigma}^\dagger c_{\mathbf k\sigma} \\ & - V \sum_{\mathbf k,\mathbf k'\in\mathcal S} b_{\mathbf k}^\dagger b_{\mathbf k'}. \end{aligned}

Here

ξk=ϵk−μ,\xi_{\mathbf k} = \epsilon_{\mathbf k}-\mu,

V>0V>0 denotes attraction in the displayed convention, and S\mathcal S is the pairing shell or model space.

The reduced interaction has two creation and two annihilation operators, so the exact reduced Hamiltonian conserves total particle number:

[KBCS,N]=0.[K_{\mathrm{BCS}},N] = 0.

The usual BCS mean-field approximation introduces

Δ=V∑k∈S⟨c−k↓ck↑⟩\Delta = V \sum_{\mathbf k\in\mathcal S} \left\langle c_{-\mathbf k\downarrow} c_{\mathbf k\uparrow} \right\rangle

and obtains quasiparticle energies

Ek=ξk2+∣Δ∣2.E_{\mathbf k} = \sqrt{ \xi_{\mathbf k}^2 + |\Delta|^2 }.

The unprojected mean-field state is not an eigenstate of NN. Do not confuse that approximation with number nonconservation in the exact reduced pairing Hamiltonian.

Use Reduced BCS Model for seniority sectors, Richardson equations, and the exact finite-level benchmark; BCS Mean-Field Theory for the Cooper instability, gap equation, and quasiparticle derivation; and BCS Model for the compact model card.

The model combines:

  • one spinful localized fermion orbital with states ∣0⟩\lvert0\rangle, ∣↑⟩\lvert\uparrow\rangle, ∣↓⟩\lvert\downarrow\rangle, and ∣2⟩\lvert2\rangle;
  • a Fock space of noninteracting bath fermions;
  • coherent bath–impurity tunneling.

With all one-particle energies measured relative to the equilibrium chemical potential,

HA=∑k,σϵkckσ†ckσ+ϵd∑σndσ+Und↑nd↓+∑k,σ(Vkckσ†dσ+Vk∗dσ†ckσ).\begin{aligned} H_{\mathrm A} &= \sum_{k,\sigma} \epsilon_k c_{k\sigma}^{\dagger}c_{k\sigma} + \epsilon_d \sum_\sigma n_{d\sigma} \\ &\quad+ U n_{d\uparrow}n_{d\downarrow} + \sum_{k,\sigma} \left( V_kc_{k\sigma}^{\dagger}d_\sigma + V_k^*d_\sigma^{\dagger}c_{k\sigma} \right). \end{aligned}

The atomic impurity energies are

E0=0,E↑=E↓=ϵd,E2=2ϵd+U.E_0=0, \qquad E_\uparrow=E_\downarrow=\epsilon_d, \qquad E_2=2\epsilon_d+U.

For repulsive UU, the atomic local-moment window is

ϵd<0<ϵd+U.\epsilon_d<0<\epsilon_d+U.

Hybridization means that impurity charge is not conserved even though total charge is.

The bath and tunneling amplitudes enter through

Δ(z)=∑k∣Vk∣2z−ϵk.\Delta(z) = \sum_k \frac{\lvert V_k\rvert^2}{z-\epsilon_k}.

For the retarded boundary value,

ΔR(ω)=Λ(ω)−iΓ(ω),\Delta^R(\omega) = \Lambda(\omega)-i\Gamma(\omega),

with

Γ(ω)=π∑k∣Vk∣2δ(ω−ϵk).\Gamma(\omega) = \pi \sum_k \lvert V_k\rvert^2 \delta(\omega-\epsilon_k).

Thus an unweighted bath density of states is not enough to specify the impurity problem. In the common wide-band convention, constant Γ\Gamma is the half-width of the noninteracting impurity resonance.

The Anderson Impurity Model dossier owns the compact specification, hybridization conventions, regime and exact-status map, observable–method dictionary, and finite benchmark. Anderson Impurity Model Preview owns the resonant-level derivation, impurity spectra, thermodynamics, transport, solution-method narrative, and DMFT connection. Effective Hamiltonians in Many-Body Systems owns the comparative Kondo reduction and its validity audit.

The model combines:

  • a localized impurity spin S\mathbf S;
  • a Fock space of conduction electrons;
  • a local exchange interaction at the impurity position.
HK=∑k,σϵkckσ†ckσ+JKS⋅sc(0),H_{\mathrm K} = \sum_{\mathbf k,\sigma} \epsilon_{\mathbf k} c_{\mathbf k\sigma}^\dagger c_{\mathbf k\sigma} + J_K \mathbf S\cdot \mathbf s_c(\mathbf 0),

where the local conduction-electron spin density is

sc(0)=12∑α,βψα†(0)σαβψβ(0).\mathbf s_c(\mathbf 0) = \frac12 \sum_{\alpha,\beta} \psi_\alpha^\dagger(\mathbf 0) \boldsymbol\sigma_{\alpha\beta} \psi_\beta(\mathbf 0).

The normalization of ψ(0)\psi(\mathbf0) determines the units assigned to JKJ_K. A bandwidth or ultraviolet cutoff and conduction density of states are part of the model.

With the displayed plus sign:

  • JK>0J_K>0 is antiferromagnetic and favors local singlet formation;
  • JK<0J_K<0 is ferromagnetic.

For an isolated pair of spin-1/21/2 objects,

S⋅s={−3/4,singlet,+1/4,triplet.\mathbf S\cdot\mathbf s = \begin{cases} -3/4, & \text{singlet},\\ +1/4, & \text{triplet}. \end{cases}

The many-electron Kondo problem is not merely this two-spin spectrum. Repeated low-energy scattering produces logarithmic scale dependence and, for antiferromagnetic exchange, a Kondo scale schematically of the form

TK∼Dexp⁡[−1ρtotJK],T_K \sim D \exp \left[ -\frac{1}{\rho_{\mathrm{tot}} J_K} \right],

up to convention-dependent prefactors and definitions. Here DD is a conduction-band cutoff and ρtot\rho_{\mathrm{tot}} is the spin-summed local density of states. With a per-spin density ρ0\rho_0, the same exponent is −1/(2ρ0JK)-1/(2\rho_0J_K). This formula is an asymptotic weak-coupling scale, not an exact universal equality.

Common extensions include potential scattering, anisotropic exchange, multiple channels, multiple impurities, and lattice versions. Each changes the symmetry and possible low-energy behavior.

The Kondo Model dossier owns the compact specification, normalization conventions, solvability map, controlled limits, observable–method dictionary, and finite benchmark. The Kondo Model Preview owns the one-loop flow, screening interpretation, thermodynamics, screening cloud, and low-energy local Fermi-liquid narrative.

Model as displayedTotal particle numberSpin symmetryTranslation symmetry
continuum bosonsconserveddepends on componentsif VextV_{\mathrm{ext}} and regulator permit
continuum fermionsconserveddepends on Vσσ′V_{\sigma\sigma'}if kernels are translation invariant
tight bindingconservedoften spin independentonly for a periodic uniform lattice
spinless tt–VV chainconservedno retained spin degree of freedomonly for a periodic uniform chain
HubbardconservedSU(2)SU(2) spin for the standard spin-independent formonly for a uniform lattice
Bose–Hubbardconservednot applicable for one componentonly for a uniform lattice
Heisenbergno particle numberglobal SU(2)SU(2) without anisotropy or fieldonly for a uniform exchange network
TFIMno particle numberdiscrete Z2\mathbb Z_2only for uniform couplings
reduced BCSconserved exactlyreduced by the pairing structuremomentum-space form assumes translation invariance
Anderson impuritytotal bath-plus-impurity number conservedSU(2)SU(2) for spin-independent parametersbroken by the impurity
Kondo impurityconduction number conservedSU(2)SU(2) for isotropic exchangebroken by the impurity

Additional conserved quantities can arise from integrability, point-group symmetry, bipartite structure, or special parameter choices.

Expanding a continuum field in localized orbitals,

ψ(x)≈∑iwi(x)ai,\psi(\mathbf x) \approx \sum_i w_i(\mathbf x)a_i,

turns kinetic and periodic-potential terms into hopping matrix elements and short-range interactions into onsite and extended lattice couplings. Truncating to one band is an approximation controlled only when higher bands remain sufficiently separated.

Adding an onsite interaction to spinful tight-binding fermions gives the standard single-band Hubbard structure:

HHub=Htb+U∑ini↑ni↓.H_{\mathrm{Hub}} = H_{\mathrm{tb}} + U \sum_i n_{i\uparrow}n_{i\downarrow}.

In one dimension, the Jordan–Wigner Transformation maps nearest-neighbor spinless hopping and density interaction to transverse and longitudinal spin exchange. With centered density and a staggered hopping-sign convention,

sjz=nj−12,J=2t,Δ=V2t.s_j^z = n_j-\frac12, \qquad J=2t, \qquad \Delta=\frac{V}{2t}.

For a periodic chain, the spin boundary induces a fermionic twist tied to total parity. The local bulk parameter map does not remove that finite-size boundary condition.

At exact half filling in the no-doublon sector and for U≫∣t∣U\gg |t|, virtual hopping through doubly occupied states produces antiferromagnetic exchange:

Heff=4t2U∑⟨i,j⟩(si⋅sj−14ninj)+⋯ .H_{\mathrm{eff}} = \frac{4t^2}{U} \sum_{\langle i,j\rangle} \left( \mathbf s_i\cdot\mathbf s_j - \frac14 n_i n_j \right) + \cdots.

Inside the singly occupied subspace, ninj=1n_i n_j=1 on occupied bonds, so the density term is an additive constant.

Effective Hamiltonians in Many-Body Systems derives the projection and explains why doped systems also retain projected hopping and same-order three-site terms.

Projecting an attractive interaction onto time-reversed states near a Fermi surface yields the reduced pairing model. The projection discards other scattering channels and introduces a cutoff or model shell.

When

−ϵd,ϵd+U≫Γ,-\epsilon_d, \quad \epsilon_d+U \gg \Gamma,

the empty and doubly occupied impurity states are virtual at low energy. A Schrieffer–Wolff transformation projects onto the singly occupied doublet and generates

JK≃2∣V∣2(1−ϵd+1ϵd+U)>0J_K \simeq 2\lvert V\rvert^2 \left( \frac1{-\epsilon_d} + \frac1{\epsilon_d+U} \right) >0

for the normalized local-orbital convention. The reduction fails near a charge degeneracy or in mixed valence.

The two charge denominators, potential-scattering term, and normalization conventions are derived in Effective Hamiltonians in Many-Body Systems.

Near a continuous transition, a spin Hamiltonian can produce a continuum order-parameter or quasiparticle field theory. That continuum theory is an effective long-wavelength description, not a literal replacement of Pauli matrices at every lattice scale.

For reproducible work, record:

  1. the exact Hamiltonian or grand Hamiltonian;
  2. the operator algebra and local Hilbert space;
  3. lattice graph, dimension, and boundary conditions;
  4. bond-counting convention;
  5. parameter signs and units;
  6. particle number, filling, or chemical potential;
  7. conserved symmetry sector;
  8. ultraviolet cutoff, basis truncation, or maximum occupation;
  9. energy, temperature, and momentum regime;
  10. approximation method and omitted terms;
  11. observable definitions;
  12. finite-size and thermodynamic-limit procedure.
  • Calling H−μNH-\mu N a fixed-NN Hamiltonian without stating the convention.
  • Counting each lattice bond twice while also writing an explicit Hermitian conjugate.
  • Switching between Pauli matrices and spin operators without rescaling couplings.
  • Using the same symbol JJ with opposite ferro- and antiferromagnetic sign conventions.
  • Treating g=4πℏ2as/mg=4\pi\hbar^2a_s/m as a universal bare coupling at every cutoff.
  • Writing a same-component fermion ss-wave contact term without addressing antisymmetry.
  • Calling a one-body tight-binding band a correlated many-body model.
  • Treating the spinless tt–VV chain as free when V≠0V\ne0.
  • Mixing centered and uncentered tt–VV conventions without shifting μ\mu.
  • Forgetting the parity-dependent boundary twist when the chain comes from periodic spins.
  • Forgetting that half filling in the spinful Hubbard model is N/L=1N/L=1.
  • Omitting the onsite-pair factor 1/21/2 in the Bose–Hubbard interaction.
  • Treating Jex=4t2/UJ_{\mathrm{ex}}=4t^2/U as exact outside the large-UU low-energy regime.
  • Assuming the reduced BCS Hamiltonian itself violates number conservation.
  • Treating the impurity charge as fixed in an Anderson model with nonzero hybridization.
  • Specifying an Anderson bath without the weighted hybridization function.
  • Replacing an Anderson impurity by a Kondo spin outside the local-moment regime.
  • Quoting a Kondo temperature without defining the bandwidth, density of states, and coupling convention.
  • Applying a model beyond the band, density, temperature, or interaction regime used to derive it.
  1. A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover (2003).
  2. A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press (2010).
  3. P. Coleman, Introduction to Many-Body Physics, Cambridge University Press (2015).
  4. J. Hubbard, “Electron Correlations in Narrow Energy Bands”, Proceedings of the Royal Society A 276, 238–257 (1963).
  5. M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, “Boson Localization and the Superfluid-Insulator Transition”, Physical Review B 40, 546–570 (1989).
  6. W. Heisenberg, “Zur Theorie des Ferromagnetismus”, Zeitschrift für Physik 49, 619–636 (1928).
  7. P. Pfeuty, “The One-Dimensional Ising Model with a Transverse Field”, Annals of Physics 57, 79–90 (1970).
  8. J. Bardeen, L. N. Cooper, and J. R. Schrieffer, “Theory of Superconductivity”, Physical Review 108, 1175–1204 (1957).
  9. P. W. Anderson, “Localized Magnetic States in Metals”, Physical Review 124, 41–53 (1961).
  10. J. Kondo, “Resistance Minimum in Dilute Magnetic Alloys”, Progress of Theoretical Physics 32, 37–49 (1964).
  11. J. R. Schrieffer and P. A. Wolff, “Relation between the Anderson and Kondo Hamiltonians”, Physical Review 149, 491–492 (1966).
  12. A. C. Hewson, The Kondo Problem to Heavy Fermions, Cambridge University Press (1993).
  13. P. Jordan and E. Wigner, “Über das Paulische Äquivalenzverbot”, Zeitschrift für Physik 47, 631–651 (1928).
  14. E. Lieb, T. Schultz, and D. Mattis, “Two Soluble Models of an Antiferromagnetic Chain”, Annals of Physics 16, 407–466 (1961).
  1. Hermiticity of hopping. Show that
Ht=−∑i,jtijci†cjH_t = - \sum_{i,j} t_{ij} c_i^\dagger c_j

is Hermitian when tji=tij∗t_{ji}=t_{ij}^*.

Solution

Take the adjoint:

Ht†=−∑i,jtij∗cj†ci.H_t^\dagger = - \sum_{i,j} t_{ij}^* c_j^\dagger c_i.

Relabel i↔ji\leftrightarrow j:

Ht†=−∑i,jtji∗ci†cj.H_t^\dagger = - \sum_{i,j} t_{ji}^* c_i^\dagger c_j.

If

tji∗=tij,t_{ji}^* = t_{ij},

then

Ht†=Ht.H_t^\dagger = H_t.

For a restricted bond sum, the same condition is implemented by writing the Hermitian-conjugate hopping explicitly.

  1. Onsite boson pair counting. Explain why the Bose–Hubbard interaction energy on a site with nn bosons is Un(n−1)/2Un(n-1)/2 rather than Un2/2Un^2/2.
Solution

The number of unordered pairs among nn particles is

(n2)=n(n−1)2.\binom n2 = \frac{n(n-1)}{2}.

If each onsite pair costs energy UU, then

Eint(n)=U(n2)=U2n(n−1).E_{\mathrm{int}}(n) = U \binom n2 = \frac{U}{2} n(n-1).

The expression vanishes for n=0n=0 and n=1n=1, as it should because no pair is present. The alternative Un2/2Un^2/2 contains an unphysical one-particle self-interaction unless compensated by a linear term.

  1. Number conservation in the reduced BCS model. Show by operator counting that the pair-scattering term commutes with total number.
Solution

For total fermion number NN,

[N,ci†]=ci†,[N,ci]=−ci.[N,c_i^\dagger] = c_i^\dagger, \qquad [N,c_i] = -c_i.

The pair-scattering monomial

ck↑†c−k↓†c−k′↓ck′↑c_{\mathbf k\uparrow}^\dagger c_{-\mathbf k\downarrow}^\dagger c_{-\mathbf k'\downarrow} c_{\mathbf k'\uparrow}

contains two creation and two annihilation operators. Using the commutator product rule, its net number change is

+1+1−1−1=0.+1+1-1-1 = 0.

Therefore the exact reduced BCS Hamiltonian commutes with NN. Number uncertainty enters the usual unprojected mean-field ansatz, not the exact pair-scattering operator.

  1. Exchange sign for two spins. For two spin-1/21/2 sites with
H=Js1⋅s2,H = J \mathbf s_1\cdot\mathbf s_2,

find the singlet and triplet energies and identify which is favored for each sign of JJ.

Solution

Use

s1⋅s2=12(stot2−s12−s22).\mathbf s_1\cdot\mathbf s_2 = \frac12 \left( \mathbf s_{\mathrm{tot}}^2 - \mathbf s_1^2 - \mathbf s_2^2 \right).

For two spin-1/21/2 objects,

s12=s22=34.\mathbf s_1^2 = \mathbf s_2^2 = \frac34.

The singlet has total spin 00, so

Es=−3J4.E_s = -\frac{3J}{4}.

The triplet has total spin 11, so

Et=J4.E_t = \frac{J}{4}.

Thus J>0J>0 favors the singlet and is antiferromagnetic, while J<0J<0 favors the triplet and is ferromagnetic in the displayed convention.

  1. Hubbard superexchange scale. Explain why the leading low-energy exchange near half filling scales as t2/Ut^2/U and state the standard coefficient.
Solution

In the large-UU singly occupied subspace, one hop creates a virtual doubly occupied site and an empty site. That intermediate state costs energy of order UU. A second hop returns the system to the low-energy subspace, possibly exchanging the two spins.

Second-order perturbation theory therefore produces an amplitude proportional to

t2U.\frac{t^2}{U}.

For the standard nearest-neighbor Hubbard hopping convention, adding the allowed virtual processes gives

Jex=4t2U.J_{\mathrm{ex}} = \frac{4t^2}{U}.

For U>0U>0, this is positive and therefore antiferromagnetic in the convention

HHeis=Jex∑⟨i,j⟩si⋅sj.H_{\mathrm{Heis}} = J_{\mathrm{ex}} \sum_{\langle i,j\rangle} \mathbf s_i\cdot\mathbf s_j.

The result assumes U≫tU\gg t, near-half filling, and projection to the low-energy subspace.

  1. Kondo exchange sign. Treat the impurity and one local conduction spin as two spin-1/21/2 objects. What energy splitting does JKS⋅sJ_K\mathbf S\cdot\mathbf s produce?
Solution

The same two-spin eigenvalues apply:

S⋅s={−3/4,singlet,+1/4,triplet.\mathbf S\cdot\mathbf s = \begin{cases} -3/4, & \text{singlet},\\ +1/4, & \text{triplet}. \end{cases}

Therefore

Es=−3JK4,Et=JK4,E_s = -\frac{3J_K}{4}, \qquad E_t = \frac{J_K}{4},

and the splitting is

Et−Es=JK.E_t-E_s = J_K.

Antiferromagnetic JK>0J_K>0 favors the singlet. The full Kondo effect additionally involves the continuum of conduction states and logarithmic renormalization; it is not exhausted by this two-spin calculation.

  1. Anderson impurity charge window. In the atomic limit, compare the energies of ∣0⟩\lvert0\rangle, ∣σ⟩\lvert\sigma\rangle, and ∣2⟩\lvert2\rangle. For U>0U>0, determine when the singly occupied states have the lowest energy.
Solution

The atomic energies are

E0=0,Eσ=ϵd,E2=2ϵd+U.E_0=0, \qquad E_\sigma=\epsilon_d, \qquad E_2=2\epsilon_d+U.

Single occupation lies below the empty state when

ϵd<0.\epsilon_d<0.

It lies below the doubly occupied state when

ϵd<2ϵd+U,\epsilon_d < 2\epsilon_d+U,

or

ϵd+U>0.\epsilon_d+U>0.

Therefore the atomic local-moment window is

−U<ϵd<0.-U<\epsilon_d<0.

Finite hybridization rounds the charge boundaries and makes the impurity occupancy fluctuate.