# Drude model

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 behaviour 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,

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

Here t is the time, ⟨**p**⟩ is the average momentum per electron and q, n, m, and τ are respectively the electron charge, number density, mass, and mean free time between ionic collisions (that is, the mean time an electron has traveled since the last collision). 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 hold.^{[note 1]}^{[3]}^{[4]}

The model was extended in 1905 by Hendrik Antoon Lorentz (and hence is also known as the **Drude–Lorentz model**)^{[citation needed]} to give the relation between the thermal conductivity and the electric conductivity of metals (see Lorenz number), 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.

## History[edit]

This section may be confusing or unclear to readers. In particular, history section should concern mainly historical aspects. It should avoid comparing back and forth the historical, Drude and modern perspectives. It can be merged with assumptions. (November 2020) (Learn how and when to remove this template message) |

German physicist Paul Drude proposed his model in 1900 when it was not clear whether atoms existed, and it was not clear what atoms were on a microscopic scale.^{[5]} The first direct proof of atoms through the computation of the Avogadro number from a microscopic model is due to Albert Einstein, the first modern model of atom structure dates to 1904 and the Rutherford model to 1909.
Drude starts from the discovery of electrons in 1897 by J.J. Thomson and assumes as a simplistic model of solids that the bulk of the solid is composed of positively charged scattering centers, and a sea of electrons submerge those scattering centers to make the total solid neutral from a charge perspective.^{[note 2]}

In modern terms this is reflected in the valence electron model where the sea of electrons is composed of the valence electrons only,^{[6]} and not the full set of electrons available in the solid, and the scattering centers are the inner shells of tightly bound electrons to the nucleus. The scattering centers had a positive charge equivalent to the valence number of the atoms.^{[note 3]}
This similarity added to some computation errors in the Drude paper, ended up providing a reasonable qualitative theory of solids capable of making good predictions in certain cases and giving completely wrong results in others.
Whenever people tried to give more substance and detail to the nature of the scattering centers, and the mechanics of scattering, and the meaning of the length of scattering, all these attempts ended in failures.^{[note 4]}

The scattering lengths computed in the Drude model, are of the order of 10 to 100 inter-atomic distances, and also these could not be given proper microscopic explanations. In modern terms, there are experiments in which electrons can travel for meters in a solid in the same manner as they would travel in free space, and this shows how a purely classical model cannot work.^{[7]}

Drude scattering is not electron-electron scattering which is only a secondary phenomenon in the modern theory, neither nuclear scattering given electrons can be at most be absorbed by nuclei. The model remains a bit mute on the microscopic mechanisms, in modern terms this is what is now called the "primary scattering mechanism" where the underlying phenomenon can be different case per case.^{[note 5]}

The model gives better predictions for metals, especially in regards to conductivity,^{[note 6]} and sometimes is called Drude theory of metals. This is because metals have essentially a better approximation to the free electron model, i.e. metals do not have complex band structures, electrons behave essentially as free particles and where, in the case of metals, the effective number of de-localized electrons is essentially the same as the valence number.^{[note 7]}

The same Drude theory, despite inconsistencies which baffled most physicists of the period, was the major one accepted to explain solids until the introduction in 1927 of the Drude-Sommerfeld model.

A few more hints of the correct ingredients of a modern theory of solids was given by the following:

- The Einstein solid model and the Debye model, suggesting that the quantum behaviour of exchanging energy in integral units or quantas was an essential component in the full theory especially with regard to specific heats, where the Drude theory failed.
- In some cases, namely in the Hall effect, the theory was making correct predictions if instead of using a negative charge for the electrons a positive one was used. This is now interpreted as holes (i.e. quasi-particles that behave as positive charge carriers) but at the time of Drude it was rather obscure why this was the case.
^{[note 8]}

Drude used Maxwell–Boltzmann statistics for the gas of electrons and for deriving the model, which was the only one available at that time. By replacing the statistics with the correct Fermi Dirac statistics, Sommerfeld significantly improved the predictions of the model, although still having a semi-classical theory that could not predict all results of the modern quantum theory of solids.^{[note 9]}

Nowadays Drude and Sommerfeld models are still significant to understanding the qualitative behaviour of solids and to get a first qualitative understanding of a specific experimental setup.^{[note 10]} This is a generic method in solid state physics, where it is typical to incrementally increase the complexity of the models to give more and more accurate predictions. It is less common to use a full-blown quantum field theory from first principles, given the complexities due to the huge numbers of particles and interactions and the little added value of the extra mathematics involved (considering the incremental gain in numerical precision of the predictions).^{[8]}

## Assumptions[edit]

- Drude used the kinetic theory of gases applied to the gas of electrons moving on a fixed background of "ions"; this is in contrast with the usual way of applying the theory of gases as a neutral diluted gas with no background. The number density of the electron gas was assumed to be

- where
*Z*is the effective number of de-localized electrons per ion, for which Drude used the valence number,*A*is the atomic mass number, is the amount of substance concentration of the "ions", and*N*_{A}is the Avogadro constant. - Considering the average volume available per electron as a sphere:
- The quantity is a parameter that describes the electron density and is often of the order of 2 or 3 times the Bohr radius, for alkali metals it ranges from 3 to 6 and some metal compounds it can go up to 10.
- The densities are of the order of 100 times of a typical classical gas. Despite that, Drude applied the kinetic theory of a dilute gas, ignoring electron–electron and electron–ion interactions aside from collisions.
^{[note 11]}

- The Drude model considers the metal to be formed of a collection of positively charged ions from which a number of "free electrons" were detached. These may be thought to be the valence electrons of the atoms that have become delocalized due to the electric field of the other atoms.
^{[note 12]} - The Drude model neglects long-range interaction between the electron and the ions or between the electrons; this is called the independent electron approximation.
^{[note 12]} - The electrons move in straight lines between one collision and another; this is called free electron approximation.
^{[note 12]} - The only interaction of a free electron with its environment was treated as being collisions with the impenetrable ions core.
^{[note 12]} - The average time between subsequent collisions of such an electron is τ, with a memoryless Poisson distribution. The nature of the collision partner of the electron does not matter for the calculations and conclusions of the Drude model.
^{[note 12]} - After a collision event, the distribution of the velocity and direction of an electron is determined by only the local temperature and is independent of the velocity of the electron before the collision event.
^{[note 12]}The electron is considered to be immediately at equilibrium with the local temperature after a collision.

Removing or improving upon each of these assumptions gives more refined models, that can more accurately describe different solids:

- Improving the hypothesis of the Maxwell–Boltzmann statistics with the Fermi–Dirac statistics leads to the Drude–Sommerfeld model.
- Improving the hypothesis of the Maxwell–Boltzmann statistics with the Bose–Einstein statistics leads to considerations about the specific heat of integer spin atoms
^{[9]}and to the Bose–Einstein condensate. - A valence band electron in a semiconductor is still essentially a free electron in a delimited energy range (i.e. only a "rare" high energy collision that implies a change of band would behave differently); the independent electron approximation is essentially still valid (i.e. no electron–electron scattering), where instead the hypothesis about the localization of the scattering events is dropped (in layman terms the electron is and scatters all over the place).
^{[10]}

## Mathematical treatment[edit]

### DC field[edit]

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 *d***p** between collisions, which occur on average every τ seconds.^{[note 1]}

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

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

Substituting the relations

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

### Time-varying analysis[edit]

The dynamics may also be described by introducing an effective drag force. At time *t* = *t*_{0} + *dt* the electron's momentum will be:

where can be interpreted as generic force (e.g. Lorentz Force) on the carrier or more specifically on the electron. is the momentum of the carrier with random direction after the collision (i.e. with a momentum ) and with absolute kinetic energy

- .

On average, a fraction of of the electrons will not have experienced another collision, the other fraction that had the collision on average will come out in a random direction and will contribute to the total momentum to only a factor which is of second order.^{[note 13]}

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

The second term is actually an extra drag force or damping term due to the Drude effects.

### Constant electric field[edit]

At time *t* = *t*_{0} + *dt* the average electron's momentum will be

and then

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

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

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

and the material can be shown to satisfy Ohm's law with a DC-conductivity *σ*_{0}:

### AC field[edit]

The Drude model can also predict the current as a response to a time-dependent electric field with an angular frequency ω. The complex conductivity is

Here it is assumed that:

In engineering, i is generally replaced by −i (or −j ) in all equations, which reflects the phase difference with respect to origin, rather than delay at the observation point traveling in time.

^{[note 14]}

Given

And the equation of motion above

substituting

Given

defining the complex conductivity from:

We have:

The imaginary part indicates that the current lags behind the electrical field. This 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.

If a sinusoidally varying electric field with frequency is applied to the solid, the negatively charged electrons behave as a plasma that tends to move a distance *x* apart from the positively charged background. As a result, the sample is polarized and there will be an excess charge at the opposite surfaces of the sample.

The dielectric constant of the sample is expressed as

where is the electric displacement and is the polarization density.

The polarization density is written as

and the polarization density with *n* electron density is

After a little algebra the relation between polarization density and electric field can be expressed as

The frequency dependent dielectric function of the solid is

^{[note 15]}

Given the approximations for the included above

- we assumed no electromagnetic field: this is always smaller by a factor v/c given the additional Lorentz term in the equation of motion
- we assumed spatially uniform field: this is true if the field does not oscillate considerably across a few mean free paths of electrons. This is typically not the case: the mean free path is of the order of Armstrongs corresponding to wavelengths typical of X rays.

Given the Maxwell equations without sources (which are treated separately in the scope of plasma oscillations)

then

or

which is an electromagnetic wave equation for a continuous homogeneous medium with dielectric constant in the helmoltz form

where the refractive index is and the phase velocity is therefore the complex dielectric constant is

which in the case can be approximated to:

At a resonance frequency , called the **plasma frequency**, the dielectric function changes sign from negative to positive and real part of the dielectric function drops to zero.

The plasma frequency represents a plasma oscillation resonance or plasmon. The plasma frequency can be employed as a direct measure of the square root of the density of valence electrons in a solid. Observed values are in reasonable agreement with this theoretical prediction for a large number of materials.^{[11]} Below the plasma frequency, the dielectric function is negative and the field cannot penetrate the sample. Light with angular frequency below the plasma frequency will be totally reflected. Above the plasma frequency the light waves can penetrate the sample, a typical example are alkaline metals that becomes transparent in the range of ultraviolet radiation. ^{[note 16]}

### Thermal conductivity of metals[edit]

One of the most spectacular successes of the Drude model is due to the prediction of the Wiedemann-Franz law this was due to a set of circumstances and errors in the Drude paper. Namely Drude predicted the value of the Lorenz number:

where the real values are mostly in the range between 2 and 3 for room temperatures between 0 and 100 degrees Celsius.^{[note 17]}

^{[note 15]}

First of all solids can conduct heat given the movement of the charge carrier and given the movement of atoms or ions as per Drude model. Conductors have free charge carriers namely electrons where insulators essentially don't, ions are present in both. Given good conductivity of both electricity and heat from metals and not from semiconductors, the conductivity shall be given by the charge carriers.

We define thermal current density as the flux per unit time of thermal energy across a unit area perpendicular to the flow

where k is the thermal conductivity Considering a one dimensional model the Energy of electrons depends on the temperature at the location of the collision If we imagine a gradient of temperature where the temperature drops in the positive x direction, the average electron speed is null, but the electrons coming from the higher energy size will have had the last collision on average in and come with energies the ones coming from the lower energy size will come with energies .

The total flux (half from left and half from right) is given by

and therefore in the limit of a mean free path that is small the quantity reduces to the derivative over x.

And therefore

Expanding to 3 degrees of freedom and given

or with the thermal conductivity

Here we will not take into account that the velocity v is also dependent on the temperature and therefore on the position, this will not contribute significantly. We will also not define precisely what is the energy transported by a specific group of electrons. diving by the conductivity and therefore getting rid of

Now drude introduced here two errors namely he used the classical statistical mechanics formula for which is an overestimate by a factor of 100 and the average energy still from classical mechanics which is an underestimate of a factor 100, namely the charge carriers move much faster than atoms and than what Drude could ever have imagined.

In total it remains:

This is half of the Drude result above, given that Drude also underestimated the conductivity by a factor of two due to underestimating by a factor of 2.

This was due to the fact that Drude estimated the time between last and next collision (which is actually ) as the mean time between two collisions (which is instead ) for a Poisson distribution.^{[note 18]}

### Thermopower[edit]

A generic temperature gradient when switched on in a thin bar will trigger a current of electrons towards the lower temperature side, given the experiments are done in an open circuit manner this current will accumulate on that side generating an electric field countering the electric current. This field is called thermoelectric field:

and Q is called thermopower. The estimates by Drude are a factor of 100 low given the direct dependency with the specific heat.

where the typical thermopowers at room temperature are 100 times smaller of the order of micro-volts.^{[note 19]}

^{[note 20]}

From the simple one dimensional model

Expanding to 3 degrees of freedom

The mean velocity due to the Electric field (given the equation of motion above at equilibrium)

To have a total current null we have

And as usual in the Drude case

where the typical thermopowers at room temperature are 100 times smaller of the order of micro-Volts.^{[note 19]}

## Drude response in real materials[edit]

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 *τ*^{−1} 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.^{[12]} 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 *τ*^{−1} is at much lower frequencies.^{[12]} This is the case for certain doped semiconductor single crystals,^{[13]} high-mobility two-dimensional electron gases,^{[14]} and heavy-fermion metals.^{[15]}

## Accuracy of the model[edit]

Historically, the Drude formula was first derived in a limited way, namely by assuming that the charge carriers form a classical ideal gas. Arnold Sommerfeld considered quantum theory and extended the theory to the free electron model, where the carriers follow Fermi–Dirac distribution. The conductivity predicted is the same as in the Drude model because it does not depend on the form of the electronic speed distribution.

The Drude model provides a very good explanation of DC and AC conductivity in metals, the Hall effect, and the magnetoresistance^{[note 13]} in metals near room temperature. The model also explains partly 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 also be applied to positive (hole) charge carriers.

In his original paper, Drude made an error, estimating the Lorenz number of Wiedemann–Franz law to be twice what it classically should have been, thus making it seem in agreement with the experimental value for the specific heat is about 100 times smaller than the classical prediction but this factor cancels out with the mean electronic speed that is about 100 times bigger than Drude's calculation.^{[note 21]}

## See also[edit]

## Citations[edit]

- ^
^{a}^{b}Ashcroft & Mermin 1976, pp. 6–7 **^**Ashcroft & Mermin 1976, pp. 2–3**^**Ashcroft & Mermin 1976, pp. 3 page note 4 and fig. 1.1**^**Ashcroft & Mermin 1976, pp. 3 page note 7 and fig. 1.2**^**Ashcroft & Mermin 1976, pp. 3 page note 6**^**Ashcroft & Mermin 1976, pp. 8 table 1.2**^**Ashcroft & Mermin 1976, pp. 5 table 1.1**^**Ashcroft & Mermin 1976, pp. 15 table 1.4**^**Ashcroft & Mermin 1976, pp. 4**^**Ashcroft & Mermin 1976, pp. 2**^**Ashcroft & Mermin 1976, pp. 4- ^
^{a}^{b}^{c}^{d}^{e}^{f}Ashcroft & Mermin 1976, pp. 2–6 - ^
^{a}^{b}Ashcroft & Mermin 1976, p. 11 **^**Ashcroft & Mermin 1976, pp. 16- ^
^{a}^{b}Ashcroft & Mermin 1976, pp. 17 **^**Ashcroft & Mermin 1976, pp. 18 table 1.5**^**Ashcroft & Mermin 1976, pp. 18 table 1.6**^**Ashcroft & Mermin 1976, pp. 25 prob 1- ^
^{a}^{b}Ashcroft & Mermin 1976, pp. 25 **^**Ashcroft & Mermin 1976, pp. 24**^**Ashcroft & Mermin 1976, p. 23

## References[edit]

**^**Drude, Paul (1900). "Zur Elektronentheorie der Metalle".*Annalen der Physik*.**306**(3): 566–613. Bibcode:1900AnP...306..566D. doi:10.1002/andp.19003060312.^{[dead link]}**^**Drude, Paul (1900). "Zur Elektronentheorie der Metalle; II. Teil. Galvanomagnetische und thermomagnetische Effecte".*Annalen der Physik*.**308**(11): 369–402. Bibcode:1900AnP...308..369D. doi:10.1002/andp.19003081102.^{[dead link]}**^**Edward M. Purcell (1965).*Electricity and Magnetism*. McGraw-Hill. pp. 117–122. ISBN 978-0-07-004908-6.**^**David J. Griffiths (1999).*Introduction to Electrodynamics*. Prentice-Hall. pp. 289. ISBN 978-0-13-805326-0.**^**"Niels bohr Nobel Lecture" (PDF).**^**springer, ed. (2009). ""Free" Electrons in Solids".*Free electrons in solid*. pp. 135–158. doi:10.1007/978-3-540-93804-0_6. ISBN 978-3-540-93803-3.**^**"Detail of the inelastic mean free paths in the kinetic energy range below 250 eV, for Ag, Al, GaAs, NA, PMMA, and Si". pp. fig 3.5.**^**"Solid State Physics, Lec ture 3: Drude Theory and Sommerfeld Free Electron".**^**Einstein (1924). Sitzungsberichte der Preussischen Akademie der Wissenschaften, Physikalisch-mathematische Klasse (ed.). "Quantum Theory of the Monatomic Ideal Gas": 261–267. Cite journal requires`|journal=`

(help)**^**"Solid State Physics, Lecture17: Dynamics of Electrons in Bands".**^**C. Kittel (1953–1976).*Introduction to Solid State Physics*. Wiley & Sons. ISBN 978-0-471-49024-1.- ^
^{a}^{b}M. Dressel; M. Scheffler (2006). "Verifying the Drude response".*Annalen der Physik*.**15**(7–8): 535–544. Bibcode:2006AnP...518..535D. doi:10.1002/andp.200510198. **^**M. van Exter; D. Grischkowsky (1990). "Carrier dynamics of electrons and holes in moderately doped silicon" (PDF).*Physical Review B*.**41**(17): 12140–12149. Bibcode:1990PhRvB..4112140V. doi:10.1103/PhysRevB.41.12140. hdl:11244/19898. PMID 9993669.**^**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" (PDF).*Applied Physics Letters*.**76**(6): 745–747. Bibcode:2000ApPhL..76..745B. doi:10.1063/1.125881.**^**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. PMID 16372004. S2CID 4391917.

### General[edit]

- Ashcroft, Neil; Mermin, N. David (1976).
*Solid State Physics*. New York: Holt, Rinehart and Winston. ISBN 978-0-03-083993-1.

## External links[edit]

- Heaney, Michael B (2003). "Electrical Conductivity and Resistivity". In Webster, John G. (ed.).
*Electrical Measurement, Signal Processing, and Displays*. CRC Press. ISBN 9780203009406.