By Lorentz dispersion materials, we mean materials for which the frequency
dependence of the dielectric permittivity can be described by a sum of multiple
resonance Lorentzian functions: where

 ω0m are the resonant frequencies Gm is related to the oscillator strengths Γm is the damping coefficient ε∞ is the permittivity at infinite frequency X0 is the permittivity at ω = 0 .

In the lossless case Equation 20 is directly related to the Sellmeier equation which in
the three resonances can be presented as: In the lossy case, the Sellmeier equation can be written in a generalized form,
accounting for a non-zero damping coefficient Γm as well as for anisotropy in the
dispersion properties: There are different ways to implement Equation 20 into the FDTD formalism. Here we
consider the so-called polarization equation approach in the single resonance case.
It uses the dielectric susceptibility function: and the relation between the polarization and the electric field Py = ε0x (ω) Ey .
Taking the Fourier transform of the last equation leads to the following differential
equation: Then Equations 24 and 25 are solved numerically together with the modified Equation
26: The FDTD approach can also account for a large variety of materials such as Drude
dispersion materials, perfect metal, second-order, and third-order materials.

Lorentz model only supports 2D simulation. Lorentz_Drude material that covers
Lorentz model supports both 2D and 3D simulation.

## Drude Model

Drude material in OptiFDTD is marked as Where εr∞ is the permittivity for infinity frequency, ωp is the plasma frequency, and Γ
is the collision frequency or damping factor.

Drude model only supports 2D simulation, Lorentz_Drude model that covers Drude
model supports both 2D and 3D simulation.