Lorentz oscillator model explained

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.

Derivation of electron motion

The model is derived by modeling an electron orbiting a massive, stationary nucleus as a spring-mass-damper system. 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.

Equation of motion for electron oscillator:

where

r is the displacement of charge from the rest position,t is time,\tau is the relaxation time/scattering time,k is a constant factor characteristic of the spring,m is the effective mass of the electron,

• $\backslash omega\_0\; =\; \backslash sqrt$ is the resonance frequency of the oscillator,

e is the elementary charge,E(t) is the electric field.

For time-harmonic fields:$$\backslash mathbf\; E(t)\; =\; \backslash mathbf\; E\_0\; e^$$$$\backslash mathbf\; r(t)\; =\; \backslash mathbf\; r\_0\; e^$$

The stationary solution of this equation of motion is:$$\backslash mathbf\; r(\backslash omega)\; =\; \backslash frac\; \backslash mathbf\; E(\backslash omega)$$

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.

Dipole moment

The displacement,

r , induces a dipole moment, p , given by$$\backslash mathbf\; p(\backslash omega)\; =\; -e\; \backslash mathbf\; r(\backslash omega)\; =\; \backslash hat\backslash alpha(\backslash omega)\; \backslash mathbf\; E(\backslash omega)\; .$$\hat \alpha(\omega) is the polarizability of single oscillator, given by$$\backslash hat\; \backslash alpha(\backslash omega)\; =\; \backslash frac\; \backslash frac\; .$$

Three distinct scattering regimes can be interpreted corresponding to the dominant denominator term in the dipole moment:

Regime	Condition	Dispersion Scaling	Phase Shift
Thomson scattering	\omega ^{2} \gg (\omega)/(\tau) , \omega ^{2} _{0}	1	0°
Shneider-Miles scattering	(\omega)/(\tau) \gg	\omega_0^2-\omega^2		\omega ^{2}	90°
Rayleigh scattering	\omega ^{2 \gg} _{0} \omega ^{2} , (\omega)/(\tau)	\omega ^{4}	180°

Polarization

The polarization

P 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$$\backslash mathbf\; P\; =\; N\; \backslash mathbf\; p\; =\; N\; \backslash hat\; \backslash alpha(\backslash omega)\; \backslash mathbf\; E(\backslash omega)\; .$$

Electric displacement

The electric displacement

D is related to the polarization density P by$$\backslash mathbf\; D\; =\; \backslash hat\backslash varepsilon\; \backslash mathbf\; E\; =\; \backslash mathbf\; E\; +\; 4\backslash pi\; \backslash mathbf\; P\; =\; (1\; +\; 4\backslash pi\; N\; \backslash hat\; \backslash alpha)\; \backslash mathbf\; E$$

Dielectric function

The complex dielectric function is given the following (in Gaussian units):$$\backslash hat\; \backslash varepsilon(\backslash omega)\; =\; 1\; +\; \backslash frac\; \backslash frac$$where

4 \pi N e ^{2/m} = \omega ^{2} _{p} and \omega _{p} 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$$\backslash hat\; \backslash varepsilon(\backslash omega)\; =\; \backslash varepsilon\_\; +\; \backslash sum\_\; \backslash chi\_^(\backslash omega;\; \backslash omega\_)$$where$$\backslash chi\_^(\backslash omega;\; \backslash omega\_)\; =\; \backslash frac$$and

\varepsilon _{infty} 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, s _{j} = \omega ^{2} _{p} f _{j} and f _{j} is related to the strength of the j th absorption mechanism, \Gamma _{j} = 1 /\tau .

Separating the real and imaginary components,$$\backslash hat\; \backslash varepsilon(\backslash omega)=\; \backslash varepsilon\_1(\backslash omega)\; +\; i\; \backslash varepsilon\_2(\backslash omega)=\; \backslash left[\backslash varepsilon\_\{\backslash infty\}\; +\; \backslash sum\_\{j\}\; \backslash frac\{s\_\{j\}\; (\backslash omega\_\{0,j\}^2\; -\; \backslash omega^2)\}\{\backslash left(\backslash omega\_\{0,j\}^\{2\}\; -\; \backslash omega^\{2\}\backslash right)^\{2\}\; +\; \backslash left(\backslash Gamma\_\{j\}\; \backslash omega\backslash right)^2\}\; \backslash right]\; +\; i\; \backslash left[\backslash sum\_\{j\}\; \backslash frac\{s\_\{j\}\; (\backslash Gamma\_\{j\}\; \backslash omega)\}\{\backslash left(\backslash omega\_\{0,j\}^\{2\}\; -\; \backslash omega^\{2\}\backslash right)^\{2\}\; +\; \backslash left(\backslash Gamma\_\{j\}\; \backslash omega\backslash right)^\{2\}\}\; \backslash right]$$

Complex conductivity

The complex optical conductivity in general is related to the complex dielectric function$$\backslash hat\; \backslash sigma(\backslash omega)\; =\; \backslash frac\; \backslash left(\backslash hat\backslash varepsilon(\backslash omega)\; -\; 1\backslash right)$$

Substituting the formula of

\hat\varepsilon(\omega) in the equation above we obtain$$\backslash hat\; \backslash sigma(\backslash omega)\; =\; \backslash frac\; \backslash frac$$

Separating the real and imaginary components,$$\backslash hat\; \backslash sigma(\backslash omega)\; =\; \backslash sigma\_1(\backslash omega)\; +\; i\; \backslash sigma\_2(\backslash omega)\; =\backslash frac\; \backslash frac\; -\; i\; \backslash frac\; \backslash frac$$

See also

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