Drude model

Drude model electrons (shown here in blue) constantly bounce between heavier, stationary crystal ions (shown in red).

The Drude model of electrical conduction was proposed in 1900^[1][2] by Paul Drude to explain the transport properties of electrons in materials (especially metals). The model, which is an application of kinetic theory, assumes that the microscopic behavior of electrons in a solid may be treated classically and looks much like a pinball machine, with a sea of constantly jittering electrons bouncing and re-bouncing off heavier, relatively immobile positive ions.

The two most significant results of the Drude model are an electronic equation of motion,

$\frac{d}{dt}\mathbf{p}(t) = q\left(\mathbf{E}+\frac{\mathbf{p}(t)\times\mathbf{B}}{m}\right) - \frac{\mathbf{p}(t)}{\tau},$

and a linear relationship between current density J and electric field E,

$\mathbf{J} = \left( \frac{n q^2 \tau}{m} \right) \mathbf{E}.$

Here t is the time and p, q, n, m, and τ are respectively an electron's momentum, charge, number density, mass, and mean free time between ionic collisions (that is, the mean time an electron has traveled since the last collision, not the average time between collisions). The latter expression is particularly important because it explains in semi-quantitative terms why Ohm's law, one of the most ubiquitous relationships in all of electromagnetism, should be true.^[3][4][5]

The model was extended in 1905 by Hendrik Antoon Lorentz (and hence is also known as the Drude–Lorentz model) and is a classical model. Later it was supplemented with the results of quantum theory in 1933 by Arnold Sommerfeld and Hans Bethe, leading to the Drude–Sommerfeld model.

Assumptions

The Drude model considers the metal to be formed of a mass of positively charged ions from which a number of "free electrons" were detached. These may be thought to have become delocalized when the valence levels of the atom came in contact with the potential of the other atoms.^[6]

The Drude model neglects any long-range interaction between the electron and the ions or between the electrons. The only possible interaction of a free electron with its environment is via instantaneous collisions. The average time between subsequent collisions of such an electron is τ, and the nature of the collision partner of the electron does not matter for the calculations and conclusions of the Drude model.^[6]

Explanations

DC field

The simplest analysis of the Drude model assumes that electric field E is both uniform and constant, and that the thermal velocity of electrons is sufficiently high such that they accumulate only an infinitesimal amount of momentum dp between collisions, which occur on average every τ seconds.^[3]

Then an electron isolated at time t will on average have been traveling for time τ since its last collision, and consequently will have accumulated momentum

$\Delta\langle\mathbf{p}\rangle= q \mathbf{E} \tau.$

During its last collision, this electron will have been just as likely to have bounced forward as backward, so all prior contributions to the electron's momentum may be ignored, resulting in the expression

$\langle\mathbf{p}\rangle = q \mathbf{E} \tau.$

Substituting the relations

$\langle\mathbf{p}\rangle = m \langle\mathbf{v}\rangle,$

$\mathbf{J} = n q \langle\mathbf{v}\rangle,$

results in the formulation of Ohm's law mentioned above:

$\mathbf{J} = \left( \frac{n q^2 \tau}{m} \right) \mathbf{E}.$

Time-varying analysis

The dynamics may also be described by introducing an effective drag force. At time t = t₀ + dt the average electron's momentum will be

$\langle\mathbf{p}(t_0+dt)\rangle=\left( 1 - \frac{dt}{\tau} \right) \left(\langle\mathbf{p}(t_0)\rangle + q\mathbf{E}dt\right),$

because, on average, a fraction of 1 − dt/τ of the electrons will not have experienced another collision, and the ones that have will contribute to the total momentum to only a negligible order.^[7]

With a bit of algebra and dropping terms of order dt^2, this results in the differential equation

$\frac{d}{dt}\langle\mathbf{p}(t)\rangle = q\mathbf{E} - \frac{\langle\mathbf{p}(t)\rangle}{\tau},$

where ⟨p⟩ denotes average momentum and q the charge of the electrons. This, which is an inhomogeneous differential equation, may be solved to obtain the general solution of

$\langle\mathbf{p}(t)\rangle = q \tau \mathbf{E}(1-e^{-t/\tau}) + \langle\mathbf{p(0)}\rangle e^{-t/\tau}$

for p(t). The steady state solution (d ⟨p⟩/dt = 0) is then

$\langle\mathbf{p}\rangle = q \tau \mathbf{E}.$

As above, average momentum may be related to average velocity and this in turn may be related to current density,

$\langle\mathbf{p}\rangle = m \langle\mathbf{v}\rangle,$

$\mathbf{J} = n q \langle\mathbf{v}\rangle,$

and the material can be shown to satisfy Ohm's law with a DC-conductivity σ₀:

$\mathbf{J} = \left( \frac{n q^2 \tau}{m} \right) \mathbf{E}.$

Complex conductivity for different frequencies assuming that τ = 10⁻⁵ and that σ₀ = 1.

The Drude model can also predict the current as a response to a time-dependent electric field with an angular frequency ω, in which case

$\sigma(\omega) = \frac{\sigma_0}{1 + i\omega\tau}= \frac{\sigma_0}{1 + \omega^2\tau^2}- i\omega\tau\frac{\sigma_0}{1 + \omega^2\tau^2}.$

Here it is assumed that

$E(t) = \Re(E_0 e^{i\omega t});$

$J(t) = \Re(\sigma(\omega) E_0 e^{i\omega t}).$

In other conventions, i is replaced by −i in all equations. The imaginary part indicates that the current lags behind the electrical field, which happens because the electrons need roughly a time τ to accelerate in response to a change in the electrical field. Here the Drude model is applied to electrons; it can be applied both to electrons and holes; i.e., positive charge carriers in semiconductors. The curves for σ(ω) are shown in the graph.

Drude response in real materials

The characteristic behavior of a Drude metal in the time or frequency domain, i.e. exponential relaxation with time constant τ or the frequency dependence for σ(ω) stated above, is called Drude response. In a conventional, simple, real metal (e.g. sodium, silver, or gold at room temperature) such behavior is not found experimentally, because the characteristic frequency τ⁻¹ is in the infrared frequency range, where other features that are not considered in the Drude model (such as band structure) play an important role.^[8] But for certain other materials with metallic properties, frequency-dependent conductivity was found that closely follows the simple Drude prediction for σ(ω). These are materials where the relaxation rate τ⁻¹ is at much lower frequencies.^[8] This is the case for certain doped semiconductor single crystals,^[9] high-mobility two-dimensional electron gases,^[10] and heavy-fermion metals.^[11]

Accuracy of the model

Historically, the Drude formula was first derived in an incorrect way, namely by assuming that the charge carriers form an ideal gas. It is now known that they follow Fermi–Dirac distribution and have appreciable interactions, but amazingly, the result turns out to be the same as the Drude model because, as Lev Landau derived in 1957, a gas of interacting particles can be described by a system of almost non-interacting 'quasiparticles' that, in the case of electrons in a metal, can be well modelled by the Drude equation.

This simple classical Drude model provides a very good explanation of DC and AC conductivity in metals, the Hall effect, and thermal conductivity (due to electrons) in metals near room temperature. The model also explains the Wiedemann–Franz law of 1853. However, it greatly overestimates the electronic heat capacities of metals. In reality, metals and insulators have roughly the same heat capacity at room temperature. The model can be applied to positive (hole) charge carriers, as demonstrated by the Hall effect.

One note of trivia surrounding the theory is that in his original paper Drude made a conceptual error, estimating electrical conductivity to in fact be only half of what it classically should have been.^[12]

See also

• Free electron model
• Arnold Sommerfeld
• Classical and quantum conductivity
• Electrical conductivity

References

1. Drude, Paul (1900). "Zur Elektronentheorie der metalle". Annalen der Physik 306 (3): 566. Bibcode:1900AnP...306..566D. doi:10.1002/andp.19003060312.
2. Drude, Paul (1900). "Zur Elektronentheorie der Metalle; II. Teil. Galvanomagnetische und thermomagnetische Effecte". Annalen der Physik 308 (11): 369. Bibcode:1900AnP...308..369D. doi:10.1002/andp.19003081102.
3. Neil Ashcroft; David Mermin (1976). Solid State Physics. Saunders College. pp. 6–7. ISBN 0-03-083993-9.
4. Edward M. Purcell (1965). Electricity and Magnetism. McGraw-Hill. pp. 117–122. ISBN 978-0-07-004908-6.
5. David J. Griffiths (1999). Introduction to Electrodynamics. Prentice-Hall. p. 289. ISBN 978-0-13-805326-0.
6. Neil Ashcroft; N. David Mermin (1976). Solid State Physics. Saunders College. pp. 2–6. ISBN 0-03-083993-9.
7. Neil Ashcroft; N. David Mermin (1976). Solid State Physics. Saunders College. p. 11. ISBN 0-03-083993-9.
8. M. Dressel, M. Scheffler (2006). "Verifying the Drude response". Ann. Phys. 15 (7–8): 535–544. Bibcode:2006AnP...518..535D. doi:10.1002/andp.200510198.
9. M. van Exter, D. Grischkowsky (1990). "Carrier dynamics of electrons and holes in moderately doped silicon". Phys. Rev. B 41 (17): 12140–12149. Bibcode:1990PhRvB..4112140V. doi:10.1103/PhysRevB.41.12140.
10. P. J. Burke, I. B. Spielman, J. P. Eisenstein, L. N. Pfeiffer, K. W. West (2000). "High frequency conductivity of the high-mobility two-dimensional electron gas". Appl. Phys. Lett. 76 (6): 745–747. Bibcode:2000ApPhL..76..745B. doi:10.1063/1.125881.
11. M. Scheffler, M. Dressel, M. Jourdan, H. Adrian (2005). "Extremely slow Drude relaxation of correlated electrons". Nature 438 (7071): 1135–1137. Bibcode:2005Natur.438.1135S. doi:10.1038/nature04232.
12. Neil W. Ashcroft; N. David Mermin (1976). Solid State Physics. Saunders College. p. 23. ISBN 0-03-083993-9.

External links

• DrudeLorentz.com Database featuring online plotting and parameterisation of the Drude-Lorentz models of common metals