FDTD calculation of scattering from frequency-dependent materials

An efficient method to include frequency-dependent materials in finite difference time domain calculations based on the recursive evaluation of the convolution of the electric field and the susceptibility function has previously been presented. The method has been applied to various materials, including those with the Debye, Drude, and Lorentz forms of complex permittivity, and to anisotropic magnetized plasmas. Previous demonstrations of this approach have been confined to total field calculations in one dimension. In this paper the recursive convolution method is extended to three-dimensional scattered field calculations. The accuracy of the method is demonstrated by calculating scattering from spheres of various sizes composed of three different types of frequency-dependent materials