A high accuracy algorithm for 3D periodic electromagnetic scattering
We develop a highly accurate numerical method for scattering of 3D electromagnetic waves by doubly periodic structures. We approximate scattered fields using the Mller boundary integral formulation of Maxwells equations. The accuracy is achieved as singularities are isolated through the use of partitions of unity, leaving smooth, periodic integrands that can be evaluated with high accuracy using trapezoid sums. The removed singularities are resolved through a transformation to polar coordinates. The method relies on the ideas used in the free space scattering algorithm of Bruno and Kunyansky.
periodic scattering periodic helmholtz greens function ewald method mller integral equation