Skip to content

From the Anderson Model to Kondo Exchange

In the local-moment regime, empty and doubly occupied Anderson states are energetically costly but still appear virtually. A Schrieffer–Wolff transformation block-diagonalizes the Hamiltonian so that these charge excursions are represented by operators acting entirely within the singly occupied sector. To second order in hybridization, the result is antiferromagnetic spin exchange plus spin-independent potential scattering.

The transformation is a matching calculation, not an RG flow. It determines a low-energy Hamiltonian, its cutoff, and the corresponding low-energy operators. The exchange runs only after this step, and the construction fails when either charge state is active rather than virtual.

Required background. The Anderson impurity model supplies the atomic sectors, charge gaps, and hybridization convention.

Helpful background. Integrating out heavy fields supplies the distinction between exact elimination and a low-energy expansion.

Write the Anderson Hamiltonian as

H=H0+HV,H=H_0+H_V,

with

H0=kσεkckσckσ+εdσndσ+Undnd,H_0 =\sum_{k\sigma}\varepsilon_kc_{k\sigma}^\dagger c_{k\sigma} +\varepsilon_d\sum_\sigma n_{d\sigma} +Un_{d\uparrow}n_{d\downarrow}, HV=kσ(Vkckσdσ+Vkdσckσ).H_V =\sum_{k\sigma} \left(V_kc_{k\sigma}^\dagger d_\sigma +V_k^\ast d_\sigma^\dagger c_{k\sigma}\right).

Let PP project onto the impurity states ,\lvert\uparrow\rangle,\lvert\downarrow\rangle and let Q=1PQ=1-P project onto the empty and doubly occupied sectors. Hybridization is off diagonal between these subspaces:

PHVP=0.PH_VP=0.

Choose an anti-Hermitian generator S=SS=-S^\dagger satisfying

[S,H0]=HV.[S,H_0]=-H_V.

For the one-orbital model, one explicit choice is

S=kσ[Vk(1ndσˉεkεd+ndσˉεkεdU)ckσdσh.c.].\begin{aligned} S=\sum_{k\sigma}\Bigg[ V_k\Bigg(& \frac{1-n_{d\bar\sigma}} {\varepsilon_k-\varepsilon_d} +\frac{n_{d\bar\sigma}} {\varepsilon_k-\varepsilon_d-U} \Bigg)c_{k\sigma}^\dagger d_\sigma -\mathrm{h.c.}\Bigg]. \end{aligned}

The occupation projectors select the intermediate charge. If the opposite spin is absent, cdc^\dagger d empties the orbital and its low-energy denominator approaches E=εdE_-=-\varepsilon_d. If the opposite spin is present, the Hermitian-conjugate process creates double occupancy and its cost is E+=εd+UE_+=\varepsilon_d+U.

Expanding the unitary transformation gives

H~=eSHeS=H0+HV+[S,H0]+[S,HV]+12[S,[S,H0]]+=H0+12[S,HV]+O ⁣(V3Ech2),\begin{aligned} \widetilde H &=e^SHe^{-S}\\ &=H_0+H_V+[S,H_0]+[S,H_V] +\frac12[S,[S,H_0]]+\cdots\\ &=H_0+\frac12[S,H_V] +O\!\left(\frac{V^3}{E_{\mathrm{ch}}^2}\right), \end{aligned}

where Ech=min(E,E+)E_{\mathrm{ch}}=\min(E_-,E_+) denotes the charge scale. Projecting the second-order term is equivalent, for a degenerate low-energy sector of energy EPE_P, to

Heff(2)=PHVQ1QH0QEPQHVP.H_{\mathrm{eff}}^{(2)} =-PH_VQ \frac{1}{QH_0Q-E_P} QH_VP.

This form makes the physics transparent: the state leaves PP, pays a positive excitation energy in QQ, and returns to PP. The full canonical transformation is systematic; the displayed approximation retains the leading low-energy block. Bravyi, DiVincenzo, and Loss 2011, §2, pp. 2795–2802 gives a modern general formulation of this block diagonalization.

Inside the singly occupied subspace, nd=1n_d=1 and the impurity bilinear decomposes as

PdαdβP=12δαβP+Sdσβα,P d_\alpha^\dagger d_\beta P =\frac12\delta_{\alpha\beta}P +\mathbf S_d\mathbin{\cdot} \boldsymbol\sigma_{\beta\alpha},

where

Sd=12dασαβdβ.\mathbf S_d =\frac12d_\alpha^\dagger \boldsymbol\sigma_{\alpha\beta}d_\beta.

Using this identity in P[S,HV]P/2P[S,H_V]P/2 separates a spin vector from a scalar. Up to a constant energy shift and a convention-dependent normal-ordering term,

Heff=Hbath+kk[JkkSdskk+Kkkσckσckσ],H_{\mathrm{eff}} =H_{\mathrm{bath}} +\sum_{kk'}\left[ J_{kk'}\, \mathbf S_d\mathbin{\cdot}\mathbf s_{kk'} +K_{kk'}\sum_\sigma c_{k\sigma}^\dagger c_{k'\sigma} \right],

with

skk=12ckασαβckβ.\mathbf s_{kk'} =\frac12c_{k\alpha}^\dagger \boldsymbol\sigma_{\alpha\beta}c_{k'\beta}.

Before the low-energy approximation, the exchange amplitude is

Jkk=VkVk(1εkεd+1εkεd1εkεdU1εkεdU),\begin{aligned} J_{kk'}=V_kV_{k'}^\ast\Bigg(& \frac1{\varepsilon_k-\varepsilon_d} +\frac1{\varepsilon_{k'}-\varepsilon_d}\\ &-\frac1{\varepsilon_k-\varepsilon_d-U} -\frac1{\varepsilon_{k'}-\varepsilon_d-U} \Bigg), \end{aligned}

and the combined scalar amplitude is

Kkk=VkVk4(1εkεd+1εkεd+1εkεdU+1εkεdU).\begin{aligned} K_{kk'}=\frac{V_kV_{k'}^\ast}{4}\Bigg(& \frac1{\varepsilon_k-\varepsilon_d} +\frac1{\varepsilon_{k'}-\varepsilon_d}\\ &+\frac1{\varepsilon_k-\varepsilon_d-U} +\frac1{\varepsilon_{k'}-\varepsilon_d-U} \Bigg). \end{aligned}

For a local orbital coupled uniformly to NN bath modes, choose the normalization

Vk=VN,c0σ=1Nkckσ.V_k=\frac{V}{\sqrt N}, \qquad c_{0\sigma}=\frac1{\sqrt N}\sum_k c_{k\sigma}.

The corresponding density of states is normalized per spin and per site,

ρ(ε)=1Nkδ(εεk),\rho(\varepsilon) =\frac1N\sum_k\delta(\varepsilon-\varepsilon_k),

so ρJ\rho J and ρK\rho K are dimensionless. Stating this normalization is essential when comparing exchange constants across continuum, lattice, and numerical conventions.

For external bath energies small compared with both charge gaps, JkkJ/NJ_{kk'}\simeq J/N and KkkK/NK_{kk'}\simeq K/N, where

J=2V2(1E+1E+)>0,J =2\lvert V\rvert^2 \left(\frac1{E_-}+\frac1{E_+}\right)>0, K=V22(1E1E+).K =\frac{\lvert V\rvert^2}{2} \left(\frac1{E_-}-\frac1{E_+}\right).

Equivalently,

Heff=Hbath+JSds(0)+Kσc0σc0σ,s(0)=12c0σc0.H_{\mathrm{eff}} =H_{\mathrm{bath}} +J\,\mathbf S_d\mathbin{\cdot}\mathbf s(0) +K\sum_\sigma c_{0\sigma}^\dagger c_{0\sigma}, \qquad \mathbf s(0)=\frac12c_0^\dagger\boldsymbol\sigma c_0.

These prefactors use the spin density with σ/2\boldsymbol\sigma/2. A source that omits the factor 1/21/2 from the spin-density definition prints an exchange coefficient smaller by two.

Schrieffer and Wolff 1966, Eqs. (4)–(13), pp. 491–492 performs this transformation in the original Anderson-model convention. The signs above have been translated to the common convention in which J>0J>0 in HK=JSdsH_K=J\mathbf S_d\cdot\mathbf s is antiferromagnetic.

Hewson 1993, ch. 1 gives a pedagogical Anderson-to-Kondo derivation and a useful translation among common Hamiltonian conventions.

The dedicated diagram can now be read as a derivation aid. Start with a singly occupied impurity and an incoming conduction electron of opposite spin. There are two time orderings:

  1. the impurity electron hops out first, leaving 0\lvert0\rangle at cost EE_-, and the bath electron hops in;
  2. the bath electron hops in first, creating \lvert\uparrow\downarrow\rangle at cost E+E_+, and one electron hops out.

Both return to charge one with the two spins exchanged. Fermionic signs and the spin decomposition make their exchange pieces add, whereas their scalar particle and hole pieces have opposite signs.

From a singly occupied Anderson orbital, an impurity electron can leave through an empty intermediate state of cost E minus or a bath electron can enter through a doubly occupied state of cost E plus; both paths return with exchanged spins and add in antiferromagnetic J.

The two second-order paths through the empty and doubly occupied charge sectors. Their exchange amplitudes add as E1+E+1E_-^{-1}+E_+^{-1}, while the scalar amplitudes subtract as E1E+1E_-^{-1}-E_+^{-1}. For a particle–hole-symmetric bath at εd=U/2\varepsilon_d=-U/2, the two costs are equal, so potential scattering vanishes but antiferromagnetic exchange remains. Original schematic for a spin-degenerate orbital and a normal bath; not to scale.

The same information has the following nonvisual form.

Virtual sectorEnergy costContribution to J/(2V2)J/(2\lvert V\rvert^2)Contribution to 2K/V22K/\lvert V\rvert^2
Empty impurity 0\lvert0\rangleEE_-+E1+E_-^{-1}+E1+E_-^{-1}
Double impurity \lvert\uparrow\downarrow\rangleE+E_++E+1+E_+^{-1}E+1-E_+^{-1}

For J>0J>0, the two-spin energies follow from

Sds=12(Stot23434).\mathbf S_d\mathbin{\cdot}\mathbf s =\frac12\left(S_{\mathrm{tot}}^2-\frac34-\frac34\right).

The singlet has eigenvalue 3/4-3/4 and the triplet has +1/4+1/4, so the singlet is lower by JJ. This is the operational meaning of an antiferromagnetic exchange sign.

For a particle–hole-symmetric bath at the symmetric impurity point,

εd=U2,E=E+=U2,\varepsilon_d=-\frac U2, \qquad E_-=E_+=\frac U2,

and therefore

J=8V2U,K=0.J=\frac{8\lvert V\rvert^2}{U}, \qquad K=0.

The cancellation of KK is an independent symmetry check. It is not a cancellation of charge fluctuations: both charge sectors still contribute equally to JJ.

Matching scale, errors, and transformed observables

Section titled “Matching scale, errors, and transformed observables”

The constant-coupling Kondo model is valid for external energies below both charge gaps. Its natural starting cutoff is of order

D0min(Dband,E,E+),D_0\lesssim \min(D_{\mathrm{band}},E_-,E_+),

with order-one scheme dependence absorbed into how JJ, KK, and the eventual Kondo scale are defined. Starting a constant-JJ RG flow at a much larger band cutoff double-counts energy ranges in which the Anderson charge states were still dynamical. The next page fixes a density-of-states convention and runs JJ below this matching scale.

The leading control conditions are

ΓE1,ΓE+1,ω,TE,E+.\frac{\Gamma_\ast}{E_-}\ll1, \qquad \frac{\Gamma_\ast}{E_+}\ll1, \qquad \lvert\omega\rvert,T\ll E_-,E_+.

If either ratio is order one, the projector PP has discarded an active state. A small numerical value of a formula for JJ does not rescue the mapping in mixed valence.

Block diagonalization also acts on observables. Every operator must be matched as

Oeff=P(O+[S,O]+12[S,[S,O]]+)P.O_{\mathrm{eff}} =P\left(O+[S,O] +\frac12[S,[S,O]]+\cdots\right)P.

For example, the impurity charge is fixed to one at zeroth order inside PP, but its transformed operator contains virtual bath–impurity pieces whose expectation value records O(V2/E±2)O(V^2/E_\pm^2) charge fluctuations. A tunnelling operator begins with an O(V/E±)O(V/E_\pm) composite operator in the Kondo theory. Matching only the Hamiltonian and then using the bare Anderson operator gives incorrect spectroscopy even if the energy spectrum is reproduced.

Nearby orbital multiplets, Hund coupling, spin–orbit coupling, or several conserved bath channels require a larger PP space. The result may be anisotropic Kondo exchange, Coqblin–Schrieffer exchange, orbital Kondo physics, or additional scalar and tensor operators. A superconducting or pseudogap bath changes the subsequent infrared flow, but it does not remove the requirement that charge elimination be controlled.

The chapter-wide map locates this derivation between Anderson charge dynamics and RG. The displayed one-orbital matching route has one active hybridization eigenchannel and therefore supplies k=1k=1. Genuine k>1k>1 screening instead enters through the separate microscopic lane with equivalent conserved baths and channel-resolved exchanges JaJ_a.

Separated charge gaps license the displayed one-orbital Schrieffer–Wolff path and its single active Kondo channel, while mixed valence blocks that reduction and genuine multichannel models require a separate channel-preserving matching path.

Matching and RG are distinct operations. For the one-orbital route, Schrieffer–Wolff eliminates separated charge sectors and supplies the k=1k=1 theory’s JJ, KK, transformed observables, and starting cutoff. A genuine multichannel route must instead supply and test the channel-resolved JaJ_a and any impurity-spin-assisted off-diagonal exchange. RG then evolves the applicable spin theory; its weak-coupling divergence is a handoff to a nonperturbative fixed-point solution. Original schematic, not to scale.

The impurity claim test matrix records mixed valence and an omitted nearby multiplet as decisive failures of this projection.

Using the right formula in the wrong regime. The denominators can be inserted algebraically even when EE_- or E+E_+ is comparable with Γ\Gamma. In that case the spin-only Hilbert space is not separated and the result is not controlled matching.

Hiding the spin-density convention. Whether σ/2\boldsymbol\sigma/2 is included in s\mathbf s changes the printed coefficient called JJ. State the Hamiltonian and compare the dimensionless product ρJ\rho J only after translating conventions.

Treating matching as the Kondo solution. Positive JJ identifies an antiferromagnetic low-energy interaction. It does not compute the universal strong-coupling phase shift, screening cloud, or Kondo-temperature prefactor.

Forgetting operator matching. A unitary transformation preserves matrix elements only when states and operators are transformed together. Bare impurity charge or tunnelling operators do not act purely inside the Kondo spin space.

1. Detune from particle–hole symmetry. Assume a particle–hole-symmetric bath and energy-independent local hybridization, then put εd=U/2+δ\varepsilon_d=-U/2+\delta with δ<U/2\lvert\delta\rvert<U/2. Express JJ and KK in terms of UU and δ\delta, and determine their parity under δδ\delta\to-\delta.

Solution

The charge gaps are E=U/2δE_-=U/2-\delta and E+=U/2+δE_+=U/2+\delta. Hence

J=2V2(1U/2δ+1U/2+δ)=8V2UU24δ2,J =2\lvert V\rvert^2 \left( \frac1{U/2-\delta} +\frac1{U/2+\delta} \right) =\frac{8\lvert V\rvert^2U} {U^2-4\delta^2}, K=V22(1U/2δ1U/2+δ)=4V2δU24δ2.K =\frac{\lvert V\rvert^2}{2} \left( \frac1{U/2-\delta} -\frac1{U/2+\delta} \right) =\frac{4\lvert V\rvert^2\delta} {U^2-4\delta^2}.

JJ is even and KK is odd in the detuning, exactly as particle–hole conjugation requires. Both expressions become unreliable as δU/2\lvert\delta\rvert\to U/2 because one charge gap closes.

2. Recover the antiferromagnetic level ordering. Diagonalize JSdsJ\mathbf S_d\cdot\mathbf s for two spin-1/21/2 degrees of freedom and show which state is lowered for J>0J>0.

Solution

Using Stot=Sd+s\mathbf S_{\mathrm{tot}}=\mathbf S_d+\mathbf s,

Sds=12[Stot(Stot+1)3434].\mathbf S_d\mathbin{\cdot}\mathbf s =\frac12 \left[ S_{\mathrm{tot}}(S_{\mathrm{tot}}+1) -\frac34-\frac34 \right].

For the singlet, Stot=0S_{\mathrm{tot}}=0 and the exchange energy is 3J/4-3J/4. For the triplet, Stot=1S_{\mathrm{tot}}=1 and it is +J/4+J/4. The singlet lies lower by JJ, so J>0J>0 is antiferromagnetic in the stated convention.

3. Test a proposed matching. Let U=8meVU=8\,\mathrm{meV}, εd=3meV\varepsilon_d=-3\,\mathrm{meV}, Γ=0.20meV\Gamma_\ast=0.20\,\mathrm{meV}, and a flat-band cutoff Dband=20meVD_{\mathrm{band}}=20\,\mathrm{meV}. Find E±E_\pm, D0maxD_0^{\max}, and the two charge-control ratios. Would the same mapping be controlled if εd=0.25meV\varepsilon_d=-0.25\,\mathrm{meV}?

Solution

For the first point,

E=3meV,E+=5meV,E_-=3\,\mathrm{meV}, \qquad E_+=5\,\mathrm{meV},

so a constant-coupling Kondo description should start no higher than an order-one multiple of D0max=3meVD_0^{\max}=3\,\mathrm{meV}, not at the full 20meV20\,\mathrm{meV} band cutoff. The ratios are Γ/E=0.067\Gamma_\ast/E_-=0.067 and Γ/E+=0.040\Gamma_\ast/E_+=0.040, giving a reasonably separated moment regime.

At εd=0.25meV\varepsilon_d=-0.25\,\mathrm{meV}, E=0.25meVE_-=0.25\,\mathrm{meV} and Γ/E=0.8\Gamma_\ast/E_-=0.8. Empty-state fluctuations are active, so the spin-only matching is uncontrolled even though UU remains much larger than Γ\Gamma_\ast.

  • Bravyi, S., DiVincenzo, D. P., and Loss, D. (2011). “Schrieffer–Wolff transformation for quantum many-body systems.” Annals of Physics 326, 2793–2826. doi:10.1016/j.aop.2011.06.004. Open PDF.
  • Hewson, A. C. (1993). The Kondo Problem to Heavy Fermions. Cambridge University Press, ch. 1. doi:10.1017/CBO9780511470752.
  • Schrieffer, J. R., and Wolff, P. A. (1966). “Relation between the Anderson and Kondo Hamiltonians.” Physical Review 149, 491–492. doi:10.1103/PhysRev.149.491.