Drude-Lorentz Model for Dispersion in Dielectrics

Drude and Lorentz (ca. 1900) developed a classical theory to account for the complex index of refraction and dielectric constant of materials, as well as their variations with the frequency of light. The model is based on treating electrons as damped harmonically bound particles subject to external electric fields. A highly simplified version of the model is given in this Demonstration, with results limited to a qualitative level. Still, the phenomena of normal and anomalous dispersion and their relation to the absorption of radiation can be quite reasonably accounted for. The classical parameters of the theory transform simply to their quantum analogs, so that the results remain valid in modern theories of materials science.
Usually the dielectric constant increases slowly with frequency—normal dispersion. However, in the neighborhood of an atomic transition the material exhibits anomalous dispersion, in which the dielectric constant decreases sharply with frequency, accompanied by absorption of light.
In this Demonstration, you can vary the parameters and . A checkbox highlights the region of anomalous dispersion.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The equation of motion for a harmonically bound classical electron interacting with an electric field is given by the Drude–Lorentz equation , where is the natural frequency of the oscillator and is the damping constant. For an incident electromagnetic field of frequency , at the point can be conveniently represented by a complex exponential . The steady-state solution, in complex form, is given by . The electric dipole moment of the electron, , corresponds to the macroscopic relation for the polarizability , where is the complex electric susceptibility , being the number of polarizable electrons per unit volume. It is sufficient to approximate , with transition frequencies lying in the optical region rad/s. Also, . The real part of the susceptibility gives the frequency dependence (dispersion) of the dielectric constant and index of refraction . The imaginary part represents the absorption coefficient. This function has the form of a Lorentzian.
In the more accurate quantum theory of dispersion, the frequency is replaced by a sum over several atomic transition frequencies and the damping parameters are determined by excited-state lifetimes.
The real and imaginary parts of the susceptibility are connected by the Kramers–Kronig relations: and , where signifies the Cauchy principal value of the integral.
[1] L. Rosenfeld, Theory of Electrons, New York: Dover Publications, 1965, pp. 68 ff.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+