View source for Drude model

[[File:Electrona in crystallo fluentia.svg|thumb|300 px|right|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<ref>{{cite journal |last= Drude |first= Paul|year= 1900 |title= Zur Elektronentheorie der Metalle |url= http://www3.interscience.wiley.com/cgi-bin/fulltext/112485959/PDFSTART |journal= [[Annalen der Physik]] | volume= 306 | issue=3 | pages=566–613 |doi= 10.1002/andp.19003060312|bibcode = 1900AnP...306..566D }}{{dead link|date=February 2019|bot=medic}}{{cbignore|bot=medic}}</ref><ref>{{cite journal |last= Drude |first= Paul |year= 1900 |title= Zur Elektronentheorie der Metalle; II. Teil. Galvanomagnetische und thermomagnetische Effecte |url= http://www3.interscience.wiley.com/cgi-bin/fulltext/112485893/PDFSTART |journal= [[Annalen der Physik]] |volume= 308 |issue=11 |pages=369–402  |doi= 10.1002/andp.19003081102|bibcode = 1900AnP...308..369D }}{{dead link|date=February 2019|bot=medic}}{{cbignore|bot=medic}}</ref> by [[Paul Karl Ludwig Drude|Paul Drude]] to explain the transport properties of [[electron]]s in materials (especially metals). The model, which is an application of [[kinetic theory of gases|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,

:<math>\frac{d}{dt}\langle\mathbf{p}(t)\rangle = q\left(\mathbf{E}+\frac{\langle\mathbf{p}(t)\rangle\times\mathbf{B}}{m}\right) - \frac{\langle\mathbf{p}(t)\rangle}{\tau},</math>

and a linear relationship between [[current density]] {{math|'''J'''}} and electric field {{math|'''E'''}},

:<math>\mathbf{J} = \left( \frac{n q^2 \tau}{m} \right) \mathbf{E}.</math>

Here {{mvar|t}} is the time, ⟨'''p'''⟩ is the average momentum per electron and {{mvar|q, n, m}}, and {{mvar|τ}} 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 be true.<ref group="note" name=":0">{{Harvnb|Ashcroft|Mermin|1976|pp=6–7}}</ref><ref>{{cite book | author = Edward M. Purcell | year = 1965 | title = Electricity and Magnetism | url = https://archive.org/details/electricitymagne00purc | url-access = registration | publisher = McGraw-Hill | pages = [https://archive.org/details/electricitymagne00purc/page/117 117–122] | isbn = 978-0-07-004908-6}}</ref><ref>{{cite book | author = David J. Griffiths | year = 1999 | title = Introduction to Electrodynamics | publisher = Prentice-Hall | pages = 289 | isbn = 978-0-13-805326-0 | url = https://archive.org/details/introductiontoel00grif_0 }}</ref>

The model was extended in 1905 by [[Hendrik Antoon Lorentz]] (and hence is also known as the '''Drude–Lorentz model'''){{Citation needed|date=August 2018}} and is a [[Classical physics|classical]] model. Later it was supplemented with the results of quantum theory in 1933 by [[Arnold Sommerfeld]] and [[Hans Bethe]], leading to the [[Free electron model|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.<ref group="note" name=":1">{{Harvnb|Ashcroft|Mermin|1976|pp=2-6}}</ref>

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 {{mvar|τ}}, and the nature of the collision partner of the electron does not matter for the calculations and conclusions of the Drude model.<ref name=":1" group="note" />

After a collision event, the velocity (and direction) of the electron only depends on the local temperature distribution and is completely independent of the velocity of the electron before the collision event.<ref name=":1" group="note" />

== Explanations ==

=== DC field ===
The simplest analysis of the Drude model assumes that electric field {{math|'''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 {{math|''d'''''p'''}} between collisions, which occur on average every {{mvar|τ}} seconds.<ref name=":0" group="note" />

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

:<math>\Delta\langle\mathbf{p}\rangle= q \mathbf{E} \tau.</math>

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

:<math>\langle\mathbf{p}\rangle = q \mathbf{E} \tau.</math>

Substituting the relations

:<math>\langle\mathbf{p}\rangle = m \langle\mathbf{v}\rangle,</math>
:<math>\mathbf{J} = n q \langle\mathbf{v}\rangle,</math>
results in the formulation of Ohm's law mentioned above:

:<math>\mathbf{J} = \left( \frac{n q^2 \tau}{m} \right) \mathbf{E}.</math>

=== Time-varying analysis ===
[[File:DrudeResponse.gif|thumb|300 px|right|Drude response of current density to an AC electric field.]]

The dynamics may also be described by introducing an effective drag force. At time {{math|''t'' {{=}} ''t''<sub>0</sub> + ''dt''}} the average electron's momentum will be

:<math>\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),</math>

because, on average, a fraction of {{math|1 − {{sfrac|''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.<ref group="note" name=":2">{{Harvnb|Ashcroft|Mermin|1976|pp=|p=11}}</ref>

With a bit of algebra and dropping terms of order {{math|''dt''<sup>2</sup>}}, this results in the differential equation

:<math>\frac{d}{dt}\langle\mathbf{p}(t)\rangle = q\mathbf{E} - \frac{\langle\mathbf{p}(t)\rangle}{\tau},</math>

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

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

for {{math|''p''(''t'')}}. The [[steady state]] solution, {{math|{{sfrac|''d'' ⟨'''p'''⟩|''dt''}} {{=}} 0}}, is then

:<math>\langle\mathbf{p}\rangle = q \tau \mathbf{E}.</math>

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

:<math>\langle\mathbf{p}\rangle = m \langle\mathbf{v}\rangle,</math>
:<math>\mathbf{J} = n q \langle\mathbf{v}\rangle,</math>

and the material can be shown to satisfy Ohm's law <math>\mathbf{J} = \sigma_0 \mathbf{E}</math> with a [[Direct current|DC]]-conductivity {{math|''σ''<sub>0</sub>}}:

:<math>\sigma_0 = \frac{n q^2 \tau}{m}</math>

=== AC field ===
[[File:DrudeModelComplexConductivity.png|thumb|300 px|right|Complex conductivity for different frequencies assuming that {{math|''τ'' {{=}} 10<sup>−5</sup>}} and that {{math|''σ''<sub>0</sub> {{=}} 1}}.]]
The Drude model can also predict the current as a response to a time-dependent electric field with an angular frequency {{mvar|ω}}. The complex conductivity is
:<math>\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}.</math>

Here it is assumed that:

:<math>E(t) = \Re\left(E_0 e^{i\omega t}\right);</math>
:<math>J(t) = \Re\left(\sigma(\omega) E_0 e^{i\omega t}\right).</math>
In engineering, {{mvar|i}} is generally replaced by {{mvar|−i}} (or {{mvar|−j}} ) in all equations, which reflects the phase difference with respect to origin, rather than delay at the observation point traveling in time. 
The imaginary part indicates that the current lags behind the electrical field. This happens because the electrons need roughly a time {{mvar|τ}} 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 {{math|''σ''(''ω'')}} are shown in the graph.

If a sinusoidally varying electric field with frequency <math>\omega</math> 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
: <math>\epsilon = \frac {D}{\epsilon_0 E} = 1 + \frac {P}{\epsilon_0 E} </math>
where <math>D</math> is the [[Electric displacement field|electric displacement]] and <math>P</math> is the [[polarization density]].

The polarization density is written as
:<math>P(t)=\Re\left(P_0e^{i\omega t}\right) </math>
and the polarization density with ''n'' electron density is
: <math>P = - n e x</math>
After a little algebra the relation between polarization density and electric field can be expressed as
: <math>P = - \frac{ne^2}{m\omega^2} E</math>
The frequency dependent dielectric function of the solid is
: <math>\epsilon(\omega) = 1 - \frac {n e^2}{\epsilon_0m \omega^2}</math>
At a resonance frequency <math>\omega_{\rm p}</math>, called the '''plasma frequency''', the dielectric function changes sign from negative to positive and real part of the dielectric function drops to zero.
: <math>\omega_{\rm p} = \sqrt{\frac{n e^2}{\epsilon_0 m}} </math>
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.<ref name="Kittel2">{{cite book|title=Introduction to Solid State Physics|author=C. Kittel|publisher=Wiley & Sons|year=1953–1976|isbn=978-0-471-49024-1}}</ref> 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.

=== 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 {{mvar|τ}} or the frequency dependence for {{math|''σ''(''ω'')}} 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 {{math|''τ''<sup>−1</sup>}} is in the infrared frequency range, where other features that are not considered in the Drude model (such as [[Electronic band structure|band structure]]) play an important role.<ref name="Dressel">{{cite journal|author1=M. Dressel |author2=M. Scheffler |year=2006|title=Verifying the Drude response|journal=[[Annalen der Physik]]|volume=15|issue=7–8|pages=535–544|doi=10.1002/andp.200510198|bibcode = 2006AnP...518..535D }}</ref> But for certain other materials with metallic properties, frequency-dependent conductivity was found that closely follows the simple Drude prediction for {{math|''σ''(''ω'')}}. These are materials where the relaxation rate {{math|''τ''<sup>−1</sup>}} is at much lower frequencies.<ref name="Dressel"/> This is the case for certain [[Doping (semiconductor)|doped semiconductor]] single crystals,<ref>{{cite journal|author1=M. van Exter |author2=D. Grischkowsky | year=1990 |title=Carrier dynamics of electrons and holes in moderately doped silicon|journal=[[Physical Review B]]|volume=41|issue=17|pages=12140–12149|doi=10.1103/PhysRevB.41.12140|bibcode = 1990PhRvB..4112140V |hdl=11244/19898 |url=https://shareok.org/bitstream/11244/19898/1/okds_Grischkowsky_PRB_1990-06-15.pdf}}</ref> high-mobility [[2DEG|two-dimensional electron gases]],<ref>{{cite journal|author1=P. J. Burke |author2=I. B. Spielman |author3=J. P. Eisenstein |author4=L. N. Pfeiffer |author5=K. W. West |year=2000|title=High frequency conductivity of the high-mobility two-dimensional electron gas|journal=[[Applied Physics Letters]]| volume=76| issue=6| pages=745–747 |doi=10.1063/1.125881|bibcode = 2000ApPhL..76..745B |url=https://authors.library.caltech.edu/2073/1/BURapl00b.pdf }}</ref> and [[Heavy fermion|heavy-fermion metals]].<ref>{{cite journal|author1=M. Scheffler |author2=M. Dressel |author3=M. Jourdan |author4=H. Adrian |year=2005|title=Extremely slow Drude relaxation of correlated electrons|journal=[[Nature (journal)|Nature]]|volume=438|pages=1135–1137|doi=10.1038/nature04232|issue=7071|bibcode = 2005Natur.438.1135S |pmid=16372004}}</ref>

== Accuracy of the model ==
Historically, the Drude formula was first derived in a limited way, namely by assuming that the charge carriers form a [[Classical physics|classical]] [[ideal gas]]. [[Arnold Sommerfeld]] considered quantum theory and extended the theory to the [[free electron model]], where the carriers follow [[Fermi–Dirac distribution]]. Amazingly, the conductivity predicted turns out to be the same as in the Drude model, as 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]]<ref name=":2" group="note" /> 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, as demonstrated by the Hall effect.

One note of trivia surrounding the theory is that in his original paper, Drude made an error, estimating the Lorenz number of Wiedemann–Franz law to be in fact twice of what it classically should have been, thus making it seem in agreement with the experimental measurement. Another surprise is that 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 actually about 100 times bigger than Drude's calculation.<ref group="note">{{Harvnb|Ashcroft|Mermin|1976|pp=|p=23}}</ref>

==See also==
*[[Free electron model]]
*[[Arnold Sommerfeld]]
*[[Electrical conductivity]]

== Notes ==
<references group="note" />

==References==
{{Reflist}}

=== General ===

* {{cite book|last=Ashcroft|first=Neil|ref=harv|author-link=Neil Ashcroft|last2=Mermin|first2=N. David|author-link2=N. David Mermin|title=Solid State Physics|publisher=Holt, Rinehart and Winston|location=New York|year=1976|isbn=978-0-03-083993-1|pages=|url=https://archive.org/details/solidstatephysic00ashc}}

==External links==
*{{Cite book|title=Electrical Measurement, Signal Processing, and Displays|last=Heaney|first=Michael B|work=|publisher=CRC Press|year=2003|isbn=9780203009406|editor-last=Webster|editor-first=John G.|location=|pages=|chapter=Electrical Conductivity and Resistivity|script-chapter=|chapter-url=https://www.researchgate.net/publication/309188334}}

{{Atomic models}}

[[Category:Condensed matter physics]]