A frequency-dependent PSTD algorithm for general dispersive media

Based on the Fourier representation of spatial derivatives, the recursive c onvolution approaches, and the perfectly matched layer (PML), we develop a pseudospectral time-domain (PSTD) algorithm for electromagnetic problems in volving general dispersive media. As a result of the accurate representatio n of spatial derivatives, the PSTD algorithm for dispersive media requires only a small number of cells per minimum wavelength, significantly reducing the number of unknowns. Multidimensional numerical results confirm the eff icacy of the PSTD algorithm for large-scale inhomogeneous media.