Lorentz oscillator model

Electrons are bound to the atomic nucleus analogously to springs of different strengths, AKA springs that are not isotropic, AKA anisotropic.

The Lorentz oscillator model describes the optical response of bound charges. The model is named after the Dutch physicist Hendrik Antoon Lorentz. It is a classical, phenomenological model for materials with characteristic resonance frequencies (or other characteristic energy scales) for optical absorption, e.g. ionic and molecular vibrations, interband transitions (semiconductors), phonons, and collective excitations.<ref name = "Lorentz1909"/><ref name = "Dressel2002"/>

Derivation of electron motion

The model is derived by modeling an electron orbiting a massive, stationary nucleus as a spring-mass-damper system.<ref name = "Dressel2002"/><ref name = "Almog2011"/><ref name = "Colton2020"/> The electron is modeled to be connected to the nucleus via a hypothetical spring and its motion is damped by via a hypothetical damper. The damping force ensures that the oscillator's response is finite at its resonance frequency. For a time-harmonic driving force which originates from the electric field, Newton's second law can be applied to the electron to obtain the motion of the electron and expressions for the dipole moment, polarization, susceptibility, and dielectric function.<ref name = "Colton2020"/>

Equation of motion for electron oscillator: <math display="block">\begin{align}

\mathbf F_\text{net} =\mathbf F_\text{damping} + \mathbf F_\text{spring} + \mathbf F_\text{driving} &= m\frac{\mathrm d^2 \mathbf r}{\mathrm dt^2} \\[1ex]
\frac{ -m}{ \tau} \frac{\mathrm d\mathbf r}{\mathrm dt} - k \mathbf r - {e} \mathbf E(t) &= m\frac{\mathrm d^2 \mathbf r}{\mathrm dt^2} \\[1ex]
\frac{\mathrm d^2 \mathbf r}{\mathrm dt^2} + \frac{ 1}{ \tau} \frac{\mathrm d\mathbf r}{\mathrm dt} + \omega_0^2 \mathbf r\; &= \; \frac{-e}{m}  \mathbf E(t)

\end{align}</math>

where

• <math>\mathbf r</math> is the displacement of charge from the rest position,
• <math>t</math> is time,
• <math> \tau</math> is the relaxation time/scattering time,
• <math>k</math> is a constant factor characteristic of the spring,
• <math>m</math> is the effective mass of the electron,
• <math display="inline">\omega_0 = \sqrt{k / m }</math> is the resonance frequency of the oscillator,
• <math>e</math> is the elementary charge,
• <math>\mathbf E(t)</math> is the electric field.

For time-harmonic fields: <math display="block">\mathbf E(t) = \mathbf E_0 e^{- i \omega t}</math> <math display="block">\mathbf r(t) = \mathbf r_0 e^{- i \omega t}</math>

The stationary solution of this equation of motion is: <math display="block">\mathbf r(\omega) = \frac{\frac{-e}{m}} {\omega_0^2 - \omega^2 - i \omega/\tau} \mathbf E(\omega) </math>

The fact that the above solution is complex means there is a time delay (phase shift) between the driving electric field and the response of the electron's motion.<ref name = "Colton2020"/>

Dipole moment

The displacement, <math>\mathbf r </math>, induces a dipole moment, <math>\mathbf p</math>, given by <math display="block">\mathbf p(\omega) = -e \mathbf r(\omega) = \hat\alpha(\omega) \mathbf E(\omega) . </math>

<math>\hat \alpha(\omega)</math> is the polarizability of single oscillator, given by <math display="block">\hat \alpha(\omega) = \frac{e^2}{m} \frac{1}{(\omega_0^2 - \omega^2) - i \omega/\tau} .</math>

Polarization

The polarization <math>\mathbf P </math> is the dipole moment per unit volume. For macroscopic material properties N is the density of charges (electrons) per unit volume. Considering that each electron is acting with the same dipole moment we have the polarization as below <math display="block">\mathbf P = N \mathbf p = N \hat \alpha(\omega) \mathbf E(\omega) . </math>

Electric displacement

The electric displacement <math>\mathbf D </math> is related to the polarization density <math>\mathbf P </math> by <math display="block">\mathbf D = \hat\varepsilon \mathbf E = \mathbf E + 4\pi \mathbf P = (1 + 4\pi N \hat \alpha) \mathbf E </math>

Dielectric function

Lorentz oscillator model. The real (blue solid line) and imaginary (orange dashed line) components of relative permittivity are plotted for a single oscillator model with parameters <math>\omega_0 = \mathrm{23.8\,THz}</math> (12.6 μm), <math>s/\omega_0^2 = \mathrm{3.305}</math>, <math>\Gamma/\omega_0 = \mathrm{0.006}</math>, and <math>\varepsilon_{\infty} = \mathrm{6.7}</math>. These parameters approximate hexagonal silicon carbide.<ref name = "Spitzer1959"/>

The complex dielectric function is given by <math display="block">\hat \varepsilon(\omega) = 1 + \frac{4\pi N e^2}{m} \frac{1}{(\omega_0^2 - \omega^2) - i \omega/\tau} </math> where <math>4\pi N e^2/m = \omega_p^2 </math> and <math> \omega_p </math> is the so-called plasma frequency.

In practice, the model is commonly modified to account for multiple absorption mechanisms present in a medium. The modified version is given by<ref name = "Zhang1998"/> <math display="block">\hat \varepsilon(\omega) = \varepsilon_{\infty} + \sum_{j} \chi_{j}^{L}(\omega; \omega_{0,j}) </math> where <math display="block">\chi_{j}^{L}(\omega; \omega_{0,j}) = \frac{s_j}{\omega_{0,j}^2 - \omega^2 - i \Gamma_j \omega} </math> and

• <math>\varepsilon_{\infty}</math> is the value of the dielectric function at infinite frequency, which can be used as an adjustable parameter to account for high frequency absorption mechanisms,
• <math>s_{j} = \omega_p^{2} f_{j}</math> and <math>f_{j}</math> is related to the strength of the <math>j</math>th absorption mechanism,
• <math>\Gamma_{j} = 1/\tau</math>.

Separating the real and imaginary components, <math display="block">\hat \varepsilon(\omega) = \varepsilon_1(\omega) + i \varepsilon_2(\omega) = \left[ \varepsilon_{\infty} + \sum_{j} \frac{s_{j} (\omega_{0,j}^2 - \omega^2)}{\left(\omega_{0,j}^{2} - \omega^{2}\right)^{2} + \left(\Gamma_{j} \omega\right)^2} \right]

 + i \left[ \sum_{j} \frac{s_{j} (\Gamma_{j} \omega)}{\left(\omega_{0,j}^{2} - \omega^{2}\right)^{2} + \left(\Gamma_{j} \omega\right)^{2}} \right]</math>

Complex conductivity

The complex optical conductivity in general is related to the complex dielectric function <math display="block"> \hat \sigma(\omega) = \frac{\omega}{4\pi i} \left(\hat\varepsilon(\omega) - 1\right) </math>

Substituting the formula of <math> \hat\varepsilon(\omega)</math> in the equation above we obtain <math display="block">\hat \sigma(\omega) = \frac{N e^2}{m} \frac{\omega}{\omega/\tau + i \left(\omega_0^2 - \omega^2 \right)} </math>

Separating the real and imaginary components, <math display="block">\hat \sigma(\omega) = \sigma_1(\omega) + i \sigma_2(\omega) = \frac{N e^2}{m} \frac{\frac{\omega^2}{\tau}}{\left(\omega_0^2 - \omega^2\right)^2 + \omega^2 / \tau^2 } - i \frac{N e^2}{m} \frac{\left(\omega_0^2 - \omega^2\right) \omega}{\left(\omega_0^2 - \omega^2\right)^2 + \omega^2/\tau^2}</math>

See also

• Cauchy equation
• Sellmeier equation
• Forouhi–Bloomer model
• Tauc–Lorentz model
• Brendel–Bormann oscillator model