paper

Robust fast direct integral equation solver for quasi-periodic scattering problems with a large number of layers

arXiv:1410.5003 · doi:10.1364/OE.23.001775

Abstract

We present a new boundary integral formulation for time-harmonic wave diffraction from two-dimensional structures with many layers of arbitrary periodic shape, such as multilayer dielectric gratings in TM polarization. Our scheme is robust at all scattering parameters, unlike the conventional quasi-periodic Green's function method which fails whenever any of the layers approaches a Wood anomaly. We achieve this by a decomposition into near- and far-field contributions. The former uses the free-space Green's function in a second-kind integral equation on one period of the material interfaces and their immediate left and right neighbors; the latter uses proxy point sources and small least-squares solves (Schur complements) to represent the remaining contribution from distant copies. By using high-order discretization on interfaces (including those with corners), the number of unknowns per layer is kept small. We achieve overall linear complexity in the number of layers, by direct solution of the resulting block tridiagonal system. For device characterization we present an efficient method to sweep over multiple incident angles, and show a speedup over solving each angle independently. We solve the scattering from a 1000-layer structure with unknowns to 9-digit accuracy in 2.5 minutes on a desktop workstation.

27 pages, 7 figures, submitted to Opt. Expr. Note that Figs. 4-6 are high resolution but appear blurry on certain PDF viewers (in this case, switch to a different PDF viewer)

References in corpus (3)

Cited by in corpus (9)