24.12.2014 Views

Set of supplementary notes.

Set of supplementary notes.

Set of supplementary notes.

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

28 CHAPTER 2. SOMMERFELD THEORY<br />

Response to an external potential<br />

The aim <strong>of</strong> this calculation is to estimate the response <strong>of</strong> a free electron gas to a perturbing<br />

potential. The perturbing potential could be caused by charges outside the metal, but it could<br />

also be due to extra charges placed inside the metal.<br />

We begin with the free electron gas in a metal, without an externally applied perturbing<br />

potential. The electrostatic potential in the metal, V 0 (r) is connected to the charge distribution<br />

ρ 0 (r) via<br />

∇ 2 V 0 (r) = − ρ 0(r)<br />

ɛ 0<br />

(2.23)<br />

In the simplest model <strong>of</strong> a metal, we consider the positive background charge to be smeared<br />

out homogeneously throughout the metal. The electron gas moves on top <strong>of</strong> this positive<br />

background. This is the plasma or ’Jellium’ model for a metal. ρ 0 = 0 everywhere in this case. 1<br />

In the presence <strong>of</strong> a perturbing potential V ext (r), the electron charge density ρ(r) will redistribute,<br />

ρ(r) = ρ 0 (r) + δρ(r), causing a correction to the potential V (r) = V 0 (r) + δV (r):<br />

∇ 2 δV (r) = − δρ(r)<br />

ɛ 0<br />

(2.24)<br />

In order to make progress, we need to link the charge density redistribution δρ to the applied<br />

potential V ext . For long-wavelength perturbations, it is plausible that in a region surrounding<br />

the position r the perturbing potential effectively just shifts the free electron energy levels,<br />

which is equivalent to assuming a spatially varying Fermi energy. This is the essence <strong>of</strong> the<br />

Thomas Fermi approximation.<br />

Thomas-Fermi approximation<br />

The Thomas-Fermi theory <strong>of</strong> screening starts with the Hartree approximation to the Schrödinger<br />

equation. The Hartree approximation is to replace the many-body pairwise interaction between<br />

the electrons by a set <strong>of</strong> interactions between a single electron and the charge density made up<br />

from all the other electrons, i.e. by a one-body potential for the i th electron<br />

U coul (r) = −<br />

e ∫<br />

4πɛ 0<br />

dr ′ ρ(r ′ )<br />

|r − r ′ | = e2 ∑<br />

∫<br />

4πɛ 0<br />

j≠i<br />

dr ′ |ψ j(r ′ )| 2<br />

|r − r ′ |<br />

, (2.25)<br />

where the summation is over all the occupied states ψ j . This clearly takes into account<br />

the averaged effect <strong>of</strong> the Coulomb repulsion due to all the other electrons. This introduces<br />

enormous simplicity, because instead <strong>of</strong> needing to solve an N-body problem, we have a (selfconsistent)<br />

one-body problem. It contains a lot <strong>of</strong> important physics, and turns out to be an<br />

approximation that is good provided the electron density is high enough. We will discuss better<br />

theories later, in the special topic <strong>of</strong> the electron gas.<br />

We shall treat the case <strong>of</strong> “jellium”, where the ionic potential is spread out uniformly to<br />

neutralise the electron liquid. Note: the average charge density is therefore always zero! The<br />

1 The correction to the charge density, δρ, does not include those charges (ρ ext ) which may have been placed<br />

inside the metal to set up the perturbing potential. They obey ∇ 2 V ext = −ρ ext /ɛ 0 .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!