Skip to main content

Physics Library

Free Physics Study Material for CSIR NET, IIT JAM, GATE, JEST & TIFR

by Pravegaa Education

Chapter 6

Free Electron Theory

The Drude Model

The assumptions

Drude imported the kinetic theory of gases wholesale. The conduction electrons of a metal are treated as a classical gas moving through a fixed array of positive ion cores:

  1. The valence electrons detach from their parent atoms and move freely through the crystal; the crystal is held together by the electrostatic attraction between the ion cores and this electron gas.
  2. Between collisions the electrons do not interact with each other (the independent electron approximation) or with the ions (the free electron approximation). The potential from the ions is taken to be perfectly uniform.
  3. Collisions are instantaneous events that randomise the electron's velocity. They occur with probability \(\mathrm{d}t/\tau\) in any interval \(\mathrm{d}t\), where \(\tau\) is the relaxation time, independent of the electron's position and velocity.
  4. After a collision the electron emerges with a velocity drawn from the Maxwell–Boltzmann distribution at the local temperature.

Assumption 4 is the one Sommerfeld will overturn in the section on failures of Drude, and the Sommerfeld Correction; assumptions 1–3 survive intact.

What Drude thought the collisions were
Drude imagined electrons bouncing off the ion cores, which gives a mean free path of order the interatomic spacing. That picture is wrong: a perfectly periodic lattice does not scatter electrons at all — an insight that has to wait for Bloch's theorem. Real scattering comes from departures from periodicity: phonons, impurities, defects and boundaries, exactly as for phonons in the section in Chapter 4. This is why measured mean free paths in pure metals at low temperature can reach centimetres, thousands of times the interatomic spacing.

Drift velocity and conductivity

Between collisions an electron in a field \(\vec{E}\) accelerates at \(-e\vec{E}/m\). If it has been travelling freely for a time \(\tau\) on average since its last collision — which randomised its velocity, so contributing nothing to the average — then the mean velocity of the electron gas is

\begin{equation} \boxed{\;\vec{v}_{d} = -\frac{e\tau}{m}\vec{E} \equiv -\mu\vec{E}, \qquad \mu = \frac{e\tau}{m} . \;} \end{equation}
The quantity \(\mu\) is the mobility, with units m²/V/s. Note that \(\vec{v}_{d}\) is a small net drift superposed on the fast random motion; it is not the speed at which electrons actually move.

Conduction in the Drude model. Between collisions (red dots) the electron travels in a straight line; each collision randomises its direction. Without a field the displacement averages to zero. With a field the trajectory is almost unchanged — the drift is enormously exaggerated here — but it acquires a small systematic component along \(-\vec{E}\), and that component is the current.

The current density follows immediately. If there are \(n\) electrons per unit volume each drifting at \(v_{d}\), then in time \(t\) every electron within a distance \(v_{d}t\) of a cross-section \(A\) passes through it, so the charge transported is \(Q = neAv_{d}t\) and

\begin{equation} \vec{J} = \frac{I}{A} = -ne\vec{v}_{d} = \frac{ne^{2}\tau}{m}\vec{E} . \end{equation}
This is Ohm's law \(\vec{J} = \sigma\vec{E}\) in local form, with
The Drude conductivity
\begin{equation} \sigma = \frac{ne^{2}\tau}{m} = ne\mu, \qquad \rho = \frac{1}{\sigma} = \frac{m}{ne^{2}\tau}, \qquad \mu = \frac{e\tau}{m} = \frac{\sigma}{ne} . \end{equation}
Everything in the section on the Drude Model, the section on response at Finite Frequency: the Plasma Frequency and the section on thermal Conductivity and the Wiedemann–Franz Law follows from this one result. The macroscopic form is recovered by writing \(R = \rho L/A\) and \(E = V/L\).
Mobility is \(\sigma/ne\), not \(\sigma ne\)
The relation is \(\mu = \sigma/ne\), so that \(\sigma = ne\mu\). Writing \(\mu = \sigma ne\) is dimensionally impossible — check it: \(\sigma\) has units /Ω/m and \(ne\) has units C/m³, so \(\sigma/ne\) gives m²/V/s while \(\sigma ne\) does not. Dimensional analysis settles this in five seconds and is worth doing every time.

How fast do electrons actually move?

The two velocities in this problem differ by ten orders of magnitude, and keeping them apart is essential.

Worked example 6.1. Drift velocity in a copper wire
A copper wire of diameter 0.16 cm carries 10 A. Copper has density 8920 kg/m³ and atomic weight 63.5, with one conduction electron per atom. Find the current density and the drift velocity.

The electron density is

\[ n = \frac{\rho N_{A}}{M} = \frac{8920\times6.023\times10^{26}}{63.5} = \text{8.46×10²⁸ m⁻³}, \]
using \(N_{A} = \text{6.023×10²⁶ kmol⁻¹}\) to match the mass in kg. With \(r = \text{8×10⁻⁴ m}\),
\[ J = \frac{I}{\pi r^{2}} = \frac{10}{\pi\left(8\times10^{-4}\right)^{2}} = \text{4.97×10⁶ A/m²}, \]
and therefore
\[ v_{d} = \frac{J}{ne} = \frac{4.97\times10^{6}}{8.46\times10^{28}\times1.6\times10^{-19}} = \text{3.67×10⁻⁴ m/s} . \]
Three very different speeds
That drift velocity is 0.37 mm/s — an electron takes about three quarters of an hour to travel a metre of wire. Compare:
SpeedOrder of magnitudeWhat it is
drift velocity \(v_{d}\)10⁻⁴ m/snet motion producing the current
Fermi velocity \(v_{F}\)10⁶ m/sactual electron speed at \(T=0\)
signal speed10⁸ m/spropagation of the field, \(\sim c\)
A light switches on immediately not because electrons rush to it, but because the electric field establishes itself throughout the circuit at nearly the speed of light and every electron in the wire starts drifting almost at once. Any answer giving \(v_{d}\) larger than about 10⁻² m/s for an ordinary wire should be suspected of a power-of-ten error.
Worked example 6.2. Mobility and relaxation time in silver
A silver wire has resistivity \(\rho = \text{1.54×10⁻⁸ Ωm}\) and \(n = \text{5.8×10²⁸ m⁻³}\). Find the mobility and the relaxation time, and the drift velocity in a field of 1 V/cm.

From the relation above,

\[ \mu = \frac{\sigma}{ne} = \frac{1}{\rho n e} = \frac{1}{1.54\times10^{-8}\times5.8\times10^{28}\times1.6\times10^{-19}} = \text{6.99×10⁻³ m²/V/s}, \]
\[ \tau = \frac{\mu m}{e} = \frac{6.99\times10^{-3}\times9.1\times10^{-31}}{1.6\times10^{-19}} = \text{3.97×10⁻¹⁴ s} . \]
With \(E = \text{100 V/m}\), \(v_{d} = \mu E = \text{0.70 m/s}\) — but note that this field is enormous for a metal: it would drive \(J = \sigma E = \text{6.5×10⁹ A/m²}\) and vaporise the wire. Realistic fields in a metal are millivolts per metre, which is why realistic drift velocities are the fractions of a millimetre per second found above.
Worked example 6.3. Resistivity of sodium from its structure
Sodium is body-centred cubic with two atoms per conventional cell and atomic radius 1.85 Å, contributing one electron each. Taking \(\tau = \text{3×10⁻¹⁴ s}\), find the resistivity.

For bcc the atoms touch along the body diagonal, \(\sqrt{3}a = 4R\), so

\[ a = \frac{4R}{\sqrt{3}} = \frac{4\times1.85}{1.732} = \text{4.27 Å}, \qquad n = \frac{2}{a^{3}} = \frac{2}{\left(4.27\times10^{-10}\right)^{3}} = \text{2.57×10²⁸ m⁻³} . \]
Then
\[ \rho = \frac{m}{ne^{2}\tau} = \frac{9.1\times10^{-31}} {2.57\times10^{28}\times\left(1.6\times10^{-19}\right)^{2} \times3\times10^{-14}} = \text{4.6×10⁻⁸ Ωm}, \]
against a measured 4.7×10⁻⁸ Ωm at room temperature. Good agreement — but see the section on failures of Drude, and the Sommerfeld Correction for why the value of \(\tau\) that makes this work is not one Drude could have justified.
Worked example 6.4. Current through an aluminium wire
Aluminium is trivalent, atomic weight 27, density 2.7 g/cm³, with mean collision time \(\tau = \text{4×10⁻¹⁴ s}\). Find the current when 3 V is applied across 20 m of wire of cross-section 2 mm², and find the lattice constant.

Each atom contributes three electrons, so

\[ n = \frac{3N_{A}\rho}{M} = \frac{3\times6.02\times10^{23}\times2.7}{27} = \text{1.81×10²³ cm⁻³} = \text{1.81×10²⁹ m⁻³}. \]
The resistance is
\[ R = \frac{m\,L}{ne^{2}\tau A} = \frac{9.1\times10^{-31}\times20} {1.81\times10^{29}\times\left(1.6\times10^{-19}\right)^{2} \times4\times10^{-14}\times2\times10^{-6}} = \text{4.9×10⁻² Ω}, \]
so \(I = V/R = 3/0.049 = \text{61 A}\).

For the lattice constant, aluminium is fcc with four atoms per cube:

\[ \frac{4}{a^{3}} = \frac{N_{A}\rho}{M} = \text{6.02×10²² cm⁻³} \quad\Longrightarrow\quad a = \left(\frac{4}{6.02\times10^{22}}\right)^{1/3} = \text{4.05×10⁻⁸ cm} = \text{4.05 Å}, \]
the accepted value. A useful check: lattice constants of metals are always a few Å. An answer below 1 Å is smaller than a single atom and must be wrong.

Response at Finite Frequency: the Plasma Frequency

The equation of motion

Drude's model handles an oscillating field just as easily as a static one. Write the average momentum \(\vec{p} = m\vec{v}_{d}\) and note that collisions damp it on the timescale \(\tau\):

\begin{equation} m\frac{\mathrm{d}\vec{v}}{\mathrm{d}t} = -e\vec{E} - \frac{m\vec{v}}{\tau} . \end{equation}
With the field switched off this gives \(v(t) = v(0)e^{-t/\tau}\): the drift decays exponentially, falling to \(1/e\) of its initial value in one relaxation time. That is what \(\tau\) means operationally.

For a driving field \(\vec{E} = \vec{E}_{0}e^{-i\omega t}\), look for a steady response \(\vec{v} = \vec{v}_{0}e^{-i\omega t}\). Substituting,

\begin{equation} -i\omega m\vec{v}_{0} = -e\vec{E}_{0} - \frac{m\vec{v}_{0}}{\tau} \quad\Longrightarrow\quad \vec{v}_{0} = \frac{-e\tau\vec{E}_{0}}{m\left(1-i\omega\tau\right)}, \end{equation}
and since \(\vec{J} = -ne\vec{v}\),
\begin{equation} \boxed{\;\sigma(\omega) = \frac{\sigma_{0}}{1-i\omega\tau} = \frac{\sigma_{0}}{1+\omega^{2}\tau^{2}} + i\,\frac{\omega\tau\,\sigma_{0}}{1+\omega^{2}\tau^{2}}, \qquad \sigma_{0} = \frac{ne^{2}\tau}{m} . \;} \end{equation}
Sign conventions in the AC response
The relation above goes with the convention \(\vec{E}\propto e^{-i\omega t}\). If you write \(\vec{E}\propto e^{+i\omega t}\) instead, the result is \(\sigma(\omega) = \sigma_{0}/(1+i\omega\tau)\). The two are complex conjugates and describe identical physics, but mixing them — quoting \(\sigma_{0}/(1-i\omega\tau)\) alongside \(e^{+i\omega t}\) — puts the sign of the reactive part backwards, turning inductive behaviour into capacitive. State the convention, then keep it.

Two limits matter.

Low frequency, \(\omega\tau\ll1\). Then \(\sigma\to\sigma_{0}\), real: the metal responds in phase with the field and dissipates energy. Since \(\tau\sim\text{10⁻¹⁴ s}\), this covers everything up to the far infrared.

High frequency, \(\omega\tau\gg1\). Then

\begin{equation} \sigma(\omega) \to \frac{\sigma_{0}}{-i\omega\tau} = \frac{ine^{2}}{m\omega}, \end{equation}
purely imaginary. The current is 90° out of phase with the field and no energy is absorbed: the electrons simply oscillate, too rapidly for collisions to intervene. The metal behaves as a lossless, inertial medium.

The plasma frequency

That inertial response has a dramatic optical consequence. The dielectric function of a medium with conductivity \(\sigma\) is \(\epsilon(\omega) = 1 + i\sigma(\omega)/\epsilon_{0}\omega\), so in the collisionless limit

\begin{equation} \epsilon(\omega) = 1 - \frac{ne^{2}}{\epsilon_{0}m\omega^{2}} = \boxed{\;1 - \frac{\omega_{p}^{2}}{\omega^{2}}, \qquad \omega_{p} = \sqrt{\frac{ne^{2}}{\epsilon_{0}m}} \;} \end{equation}
where \(\omega_{p}\) is the plasma frequency — the natural frequency at which the whole electron gas oscillates against the ion background. Its quantum is the plasmon of the table in Chapter 4.
Why metals are shiny, and where they stop being shiny
For \(\omega < \omega_{p}\), \(\epsilon < 0\), the refractive index is imaginary, and an incident wave cannot propagate: it is totally reflected. For \(\omega > \omega_{p}\), \(\epsilon > 0\) and the metal becomes transparent.

Since \(\omega_{p}\) for a typical metal corresponds to a photon energy of \(5\)–15 eV, the crossover lies in the ultraviolet. Visible light is below it and is reflected — hence metallic lustre — while the alkali metals are genuinely transparent in the ultraviolet, one of the early confirmations of the free-electron picture.

Worked example 6.5. Plasma frequency of copper
Copper has \(n = \text{8.5×10²⁸ m⁻³}\). Find \(\omega_{p}\), the corresponding photon energy, and the wavelength at which copper becomes transparent.

\[ \omega_{p} = \sqrt{\frac{ne^{2}}{\epsilon_{0}m}} = \sqrt{\frac{8.5\times10^{28}\times\left(1.602\times10^{-19}\right)^{2}} {8.854\times10^{-12}\times9.11\times10^{-31}}} = \text{1.64×10¹⁶ rad/s}, \]
so \(\hbar\omega_{p} = \text{10.8 eV}\) and
\[ \lambda_{p} = \frac{2\pi c}{\omega_{p}} = \text{115 nm}, \]
deep in the ultraviolet. Everything in the visible range (\(400\)–700 nm) is therefore reflected, which is what makes copper a mirror. Its reddish colour comes not from this mechanism but from interband transitions involving the filled \(d\) states, which the free-electron model ignores entirely.

Two consequences worth knowing

How fast a metal screens a field

Put a metal between capacitor plates and switch on a field. Charge flows to cancel it, and the cancellation is not instantaneous. Combining charge conservation \(\partial\rho/\partial t + \nabla\cdot\vec{J} = 0\) with \(\vec{J} = \sigma\vec{E}\) and Gauss's law \(\nabla\cdot\vec{E} = \rho/\epsilon_{0}\),

\begin{equation} \frac{\partial\rho}{\partial t} = -\frac{\sigma_{0}}{\epsilon_{0}}\rho \quad\Longrightarrow\quad \rho(t) = \rho(0)e^{-t/\tau_{c}}, \qquad \tau_{c} = \frac{\epsilon_{0}}{\sigma_{0}} = \frac{\epsilon_{0}m}{ne^{2}\tau} = \boxed{\;\frac{1}{\omega_{p}^{2}\tau}\;} \end{equation}
Note that this is shorter than \(\tau\) itself, by the large factor \((\omega_{p}\tau)^{2}\). A metal screens an applied field far faster than an individual electron survives between collisions.

Kinetic inductance

The drifting electrons carry kinetic energy, and a system that stores energy in proportion to \(I^{2}\) has an inductance. For a wire of length \(l\) and cross-section \(A\) carrying \(I = neAv_{d}\),

\begin{equation} \frac{1}{2}L_{k}I^{2} = \left(nAl\right)\cdot\frac{1}{2}mv_{d}^{2} \quad\Longrightarrow\quad \boxed{\;L_{k} = \frac{m\,l}{ne^{2}A} \;} \end{equation}
This kinetic inductance is negligible in ordinary wires, where the magnetic inductance dominates, but it becomes important in thin films, in superconducting transmission lines (where \(\tau\to\infty\) and there is no resistive term to mask it), and in the design of superconducting detectors. Notice that it is exactly \(\sigma_{0}^{-1}\tau\,l/A\): the inertial counterpart of the resistance \(R = l/\sigma_{0}A\).

Thermal Conductivity and the Wiedemann–Franz Law

Metals conduct heat as well as charge, and the same electrons do both. That suggests the two conductivities should be related, and the empirical Wiedemann–Franz law (1853) says they are: the ratio \(K/\sigma\) is proportional to \(T\), with a constant of proportionality that is very nearly the same for all metals.

\begin{equation} \frac{K}{\sigma T} = L, \qquad L = \text{the Lorenz number}. \end{equation}
Explaining \(L\) was Drude's greatest success, and — as it turned out — his luckiest accident.

Drude's calculation

Treat the electrons as a classical gas and use the kinetic-theory result of the relation in Chapter 4, \(K = \tfrac{1}{3}Cv\lambda\), with \(\lambda = v\tau\). Classically the electronic heat capacity per unit volume is \(C = \tfrac{3}{2}nk_{B}\) and the mean square speed satisfies \(m\langle v^{2}\rangle = 3k_{B}T\), so

\begin{equation} K = \frac{1}{3}\left(\frac{3}{2}nk_{B}\right)\langle v^{2}\rangle\tau = \frac{1}{2}nk_{B}\cdot\frac{3k_{B}T}{m}\cdot\tau = \frac{3nk_{B}^{2}T\tau}{2m} . \end{equation}
Dividing by \(\sigma T = ne^{2}\tau T/m\), the density, the relaxation time and the mass all cancel:
\begin{equation} L_{\text{classical}} = \frac{K}{\sigma T} = \frac{3}{2}\left(\frac{k_{B}}{e}\right)^{2} = \text{1.11×10⁻⁸ WΩ/K²} . \end{equation}

That every material-dependent quantity cancels is the substance of the law, and it is a genuine triumph of the model. The measured values cluster around 2.2×10⁻⁸  to 2.7×10⁻⁸ WΩ/K² — a factor of two above the relation above.

\(\langle v^{2}\rangle\), not the one-dimensional value
The equipartition result is \(\tfrac{1}{2}m\langle v^{2}\rangle = \tfrac{3}{2}k_{B}T\), i.e. \(v_{\text{rms}} = \sqrt{3k_{B}T/m}\). Writing \(v_{\text{rms}} = \sqrt{k_{B}T/m}\) — the value for a single component — gives \(L = \tfrac{1}{2}(k_{B}/e)^{2}\) and destroys the agreement. Some circulated notes quote the one-component form and then use the three-dimensional one in the Lorenz-number step, which is internally inconsistent even though the final number comes out right.
Remark. Drude originally obtained \(L \approx \text{2.2×10⁻⁸ }\), in near-perfect agreement with experiment, and it made his reputation. He had made two compensating errors: he used a heat capacity a hundred times too large (the classical \(\tfrac{3}{2}nk_{B}\) instead of the tiny Fermi-suppressed value of the relation in Chapter 5) and a mean square velocity a hundred times too small (the classical thermal value instead of \(v_{F}^{2}\)). The two errors cancel in the ratio almost exactly. It is one of the most instructive accidents in the history of physics: a right answer is not by itself evidence that a model is right.

Sommerfeld's correction

Redo the calculation with the correct quantum ingredients. From the relation in Chapter 5 the electronic heat capacity per unit volume is \(C = \tfrac{\pi^{2}}{2}nk_{B}T/T_{F}\), and the electrons that carry the heat are those at the Fermi surface, moving at \(v_{F}\) with \(v_{F}^{2} = 2k_{B}T_{F}/m\). Then

\begin{equation} K = \frac{1}{3}\left(\frac{\pi^{2}}{2}\frac{nk_{B}T}{T_{F}}\right) \left(\frac{2k_{B}T_{F}}{m}\right)\tau = \frac{\pi^{2}nk_{B}^{2}T\tau}{3m}, \end{equation}
and dividing by \(\sigma T\) as before:
The Lorenz number
\begin{equation} \boxed{\;L = \frac{K}{\sigma T} = \frac{\pi^{2}}{3}\left(\frac{k_{B}}{e}\right)^{2} = \text{2.45×10⁻⁸ WΩ/K²} \;} \end{equation}
independent of \(n\), \(\tau\), \(m\) and the material. This is in excellent agreement with experiment for good metals at room temperature and below the Debye temperature. Numerically, \(k_{B}/e = \text{8.62×10⁻⁵ V/K}\), so \(L = 3.290\times\left(8.62\times10^{-5}\right)^{2}\).
Getting \(L\) right numerically
\(\left(k_{B}/e\right)^{2} = \left(1.381\times10^{-23}/1.602\times10^{-19}\right)^{2} = \left(8.62\times10^{-5}\right)^{2} = 7.43\times10^{-9}\), and \(\pi^{2}/3 = 3.29\), giving \(2.45\times10^{-8}\). A value near \(2.8\times10^{-8}\) — which appears in several circulated solutions — is an arithmetic slip, and it matters because it changes the apparent theory-to-experiment ratio from about \(1.08\) to \(1.26\) and so invites the wrong conclusion about how well the model works.
Worked example 6.6. Lorenz number from measurements
A copper wire of length 0.5 m and diameter 0.3 mm has resistance 0.12 Ω at 293 K, and its thermal conductivity is 390 W/m/K. Find the Lorenz number and compare with theory.

From \(R = L/\sigma\pi r^{2}\) with \(r = \text{0.15×10⁻³ m}\),

\[ \sigma = \frac{L}{R\pi r^{2}} = \frac{0.5}{0.12\times\pi\times\left(0.15\times10^{-3}\right)^{2}} = \text{5.89×10⁷ /Ω/m} . \]
Hence
\[ L_{\text{expt}} = \frac{K}{\sigma T} = \frac{390}{5.89\times10^{7}\times293} = \text{2.26×10⁻⁸ WΩ/K²} . \]
Against the Sommerfeld value 2.45×10⁻⁸ , the theoretical figure is higher by a factor \(2.45/2.26 = 1.08\) — agreement to 8%, which is about as well as a free-electron model can be expected to do.

Continue reading — it's free

Register with your mobile number to unlock the full chapter and all 100+ chapters of the Physics Library.

Register / Login with OTP

Already registered? Use the same button — just verify your number.