Computing diffraction anomalies as nonlinear eigenvalue problems

TL;DR

This paper introduces an efficient numerical method using nonlinear eigenvalue problems and contour-integral techniques to accurately compute diffraction anomalies in periodic electromagnetic structures, improving over traditional methods.

Contribution

The paper develops a novel contour-integral-based approach to solve nonlinear eigenvalue problems for diffraction anomalies, enabling comprehensive analysis over frequency ranges.

Findings

Successfully applied to periodic cylinder arrays

Accurately identifies diffraction anomalies

Reduces computational effort compared to iterative methods

Abstract

When a plane electromagnetic wave impinges upon a diffraction grating or other periodic structures, reflected and transmitted waves propagate away from the structure in different radiation channels. A diffraction anomaly occurs when the outgoing waves in one or more radiation channels vanish. Zero reflection, zero transmission and perfect absorption are important examples of diffraction anomalies, and they are useful for manipulating electromagnetic waves and light. Since diffraction anomalies appear only at specific frequencies and/or wavevectors, and may require the tuning of structural or material parameters, they are relatively difficult to find by standard numerical methods. Iterative methods may be used, but good initial guesses are required. To determine all diffraction anomalies in a given frequency interval, it is necessary to repeatedly solve the diffraction problem for many frequencies. In this paper, an efficient numerical method is developed for computing diffraction anomalies. The method relies on nonlinear eigenvalue formulations for scattering anomalies and solves the nonlinear eigenvalue problems by a contour-integral method. Numerical examples involving periodic arrays of cylinders are presented to illustrate the new method.