Drude dispersion in ADE FDTD revisited

Studies of-double negative materials and of metals at optical frequencies have lead to renewed interest in modelling of materials characterised by Drude or multi-pole Lorentz-Drude frequency dispersion. The auxiliary differential equation method of treatment of Drude dispersion in the finite difference time domain method is revisited, and a substantially more accurate formulation is derived.