advanced32 min studyLesson 530 in Physics

Kondo screening and impurity physics

Follow the weak-coupling flow of a localized spin in a metal and interpret the Kondo scale as a many-body crossover.

On this page

What you’ll understand

  • Specify exchange and density-of-states normalizations in a Kondo Hamiltonian.
  • Integrate the leading logarithmic flow and estimate its breakdown scale.
  • Distinguish screening, a crossover, and extensions with different channel counts.

A localized spin in a metal is not an isolated Curie moment at arbitrarily low temperature. Antiferromagnetic exchange with a sea of conduction electrons generates collective screening and a new, often exponentially small energy scale. The impurity is coupled to many occupied and unoccupied electronic states, so a picture of a permanently bound isolated electron pair is incomplete. The Kondo problem illustrates how a weak microscopic interaction can grow as progressively lower-energy states are examined.

Declare the factor of two before defining a coupling

Consider one spin-half impurity and one spin-degenerate metallic screening channel, with a smooth approximately constant local density of states near the Fermi energy. Use dimensionless impurity spin S and conduction spin density s. For this lesson choose an interaction 2JKS⋅s(0)2J_K\mathbf S\cdot\mathbf s(0), explicitly including a factor of two. Positive JK is antiferromagnetic. This normalization differs from the coefficient of an unordered intersite spin bond earlier in the module; comparing bare numbers without the factor would be incorrect.

Let rho be the local density of electronic states per spin for that channel, in inverse energy units, and define g=ρJKg=\rho J_K. The band energies are measured relative to the Fermi energy, with initial cutoff D0 representing the half-bandwidth in the flat-band idealization. The local conduction spin uses Pauli matrices divided by two. Changing either the exchange convention or from a per-spin to a two-spin density of states changes the numerical coefficient in the flow equation and the exponent unless g is translated consistently.

H=∑kσϵkckσ†ckσ+2JKS⋅s(0),s(0)=12c0†σc0,g=ρJK,ℓ=ln⁡D0DH=\sum_{k\sigma}\epsilon_kc_{k\sigma}^\dagger c_{k\sigma}+2J_K\mathbf S\cdot\mathbf s(0),\qquad \mathbf s(0)=\frac12c_0^\dagger\boldsymbol\sigma c_0,\qquad g=\rho J_K,\quad \ell=\ln\frac{D_0}{D}
c0 is the normalized local channel orbital at the impurity, S and s are dimensionless, and JK and D have energy units. rho is the corresponding local per-spin density of states, so g and the infrared scale parameter ell are dimensionless.↔ Scroll wide equations horizontally.

Logarithms make weak antiferromagnetic exchange grow

Virtual electron and hole processes near the edges of the remaining band renormalize the low-energy spin-flip amplitude. Their denominators produce logarithms when shells of high-energy states are eliminated. To leading logarithmic order the isotropic coupling grows at a rate proportional to its square as the cutoff decreases. An initially positive g therefore increases, whereas a negative isotropic ferromagnetic coupling approaches zero from below. These are weak-coupling statements and cannot be used quantitatively after g becomes order one.

Integrating the leading flow yields a denominator that vanishes at an exponentially reduced cutoff. That formal divergence is the failure of perturbation theory, not an infinite observable or a physical thermodynamic singularity. Its location defines a leading-logarithmic estimate of the Kondo scale. Higher-order scaling and operational definitions introduce algebraic and numerical prefactors, so the exponential estimate should not be treated as a precise experimental temperature without specifying the convention and additional corrections.

dgdℓ=2g2+O(g3),g(ℓ)=g01−2g0ℓ,kBTK(1)∼D0exp⁡ ⁣(−12g0),0<g0≪1\frac{dg}{d\ell}=2g^2+O(g^3),\qquad g(\ell)=\frac{g_0}{1-2g_0\ell},\qquad k_BT_K^{(1)}\sim D_0\exp\!\left(-\frac1{2g_0}\right),\quad 0<g_0\ll1
The integrated expression retains only the 2g² term and applies while the running coupling stays small. TK(1) denotes a leading-log scale, with D0 and kBTK expressed in the same energy units. The sign is for ell increasing toward the infrared.↔ Scroll wide equations horizontally.

Screening is a many-body low-temperature crossover

For the ordinary single-channel spin-half antiferromagnetic model with a metallic bath, the low-temperature state is screened and has Fermi-liquid behavior. The impurity contribution to entropy crosses from approximately kBln⁡2k_B\ln2 for a free moment to zero at low temperature. Its susceptibility no longer follows an unlimited Curie divergence. Screening involves correlated conduction states over many scales and positions; the full conduction band is not replaced by one isolated singlet bond.

This evolution is a crossover for one impurity, not a bulk finite-temperature ordering transition. A perturbative logarithmic increase of impurity scattering must therefore not be extrapolated into a literal divergent zero-temperature resistance. The resistance minimum also depends on phonons and other backgrounds and need not equal a uniquely defined Kondo temperature. Larger impurity spins or additional independent screening channels can produce under- or overscreening, and a vanishing or structured bath density of states can change the problem qualitatively. Those extensions are outside the assumptions of the simple flow used here.

An exponentially generated energy scale

Use D0=1.0 eV and g0=0.10 in the leading-log estimate, ignoring its unknown prefactor. Find the scale kBTK(1), convert it to kelvin using kB=0.08617 meV/Kk_B=0.08617\,\mathrm{meV/K}, and identify where the one-loop running coupling formally fails.

  1. The exponent is −1/(2g0)=−5-1/(2g_0)=-5. Thus D0e−5=0.006738 eV=6.738 meVD_0e^{-5}=0.006738\,\mathrm{eV}=6.738\,\mathrm{meV}. A modest change in the initial coupling would change this scale exponentially.

  2. Dividing by Boltzmann’s constant gives TK(1)≃6.738/0.08617=78.19 KT_K^{(1)}\simeq6.738/0.08617=78.19\,\mathrm K. This is the chosen leading-log scale, not a high-precision temperature prediction for a material.

  3. The denominator of the integrated coupling vanishes at ℓ∗=1/(2g0)=5\ell_*=1/(2g_0)=5, corresponding to D=D0e−5D=D_0e^{-5}. Perturbation theory already loses quantitative control before that point as g grows to order one; the formal pole marks no actual divergent thermodynamic quantity.

The result

The leading exponential scale is 6.738 meV or about 78.19 K, and the one-loop pole occurs at ell=5 where weak-coupling extrapolation fails.

Explore the model

Interactive lab

Interactive lab · function

Inspect the quadratic weak-coupling beta function

For each fresh comparison, Reset and set BOTH Coefficient a=2 and Position x. Compare x=0.10 and x=0.20, predicting y=0.020 and 0.080. Interpreting x as g and y as its leading infrared flow rate shows a fourfold increase when g doubles. This y=ax² graph does not integrate the flow, calculate a Kondo temperature, screen an impurity, or represent an interacting many-spin state.

Inspect the quadratic weak-coupling beta functionThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. 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 forA growing positive quadratic beta function generates strong scale dependence, but the graph is only the local weak-coupling rate.

Check your understanding

What does the pole in the leading-log running coupling mean for the ordinary single-channel spin-half Kondo model?

Try it yourself

With g0=0.10 and D0=1.0 eV, compute the one-loop running coupling after lowering the cutoff to D=D0e−2D=D_0e^{-2}. Compare it with the bare coupling and calculate the cutoff numerically using e−2≈0.1353. Explain why the same formula cannot be trusted all the way to ell=5.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Read ell from the logarithm of the cutoff ratio. Use the integrated expression, not a single Euler step with the initial beta function.

Compare with the solution

Here ℓ=2\ell=2, so g(2)=0.10/[1−2(0.10)(2)]=0.10/0.60=1/6≃0.1667g(2)=0.10/[1-2(0.10)(2)]=0.10/0.60=1/6\simeq0.1667. The cutoff is D=0.1353 eVD=0.1353\,\mathrm{eV}, and the coupling has grown by a factor about 1.667. A single initial-rate step would give 0.14 and miss the rate’s growth. The one-loop expression is still an approximation even here; as ell approaches five, the predicted g becomes large and neglected higher terms cease to be small. Its formal divergence cannot be used as a physical infinite interaction or an exact resistance prediction.

Take it with you

  • A Kondo flow requires explicit exchange, channel, and density-of-states conventions.
  • Weak antiferromagnetic coupling generates an exponentially small strong-coupling scale.
  • Single-channel spin-half screening is a many-body crossover, with different possibilities in other impurity models.

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.