Title

Authors

Document Type

Article

Publication Date

9-14-2009

Department

Chemistry and Biochemistry

Abstract

We present a new approach to the numerical solution of Maxwell’s equations in the case of spatially‐varying electric permittivity and/or magnetic permeability, based on Krylov subspace spectral (KSS) methods. KSS methods for scalar equations compute each Fourier coefficient of the solution using techniques developed by Gene Golub and Gérard Meurant for approximating elements of functions of matrices by Gaussian quadrature in the spectral, rather than physical, domain. We show how they can be generalized to coupled systems of equations, such as Maxwell’s equations, by choosing appropriate basis functions that, while induced by this coupling, still allow efficient and robust computation of the Fourier coefficients of each spatial component of the electric and magnetic fields. We also discuss the implementation of appropriate boundary conditions for simulation on infinite computational domains, and how discontinuous coefficients can be handled.