Integral method for echelles covered with lossless or absorbing thin dielectric layers

High order integral equation method for diffraction gratings

Conventional integral equation methods for diffraction gratings require lattice sum techniques to evaluate quasi-periodic Green's functions. The boundary integral equation Neumann-to-Dirichlet map (BIE-NtD) method in Wu and Lu [J. Opt. Soc. Am. A 26, 2444 (2009)], [J. Opt. Soc. Am. A 28, 1191 (2011)] is a recently developed integral equation method that avoids the quasi-periodic Green's functions and is relatively easy to implement. In this paper, we present a number of improvements for this method, including a revised formulation that is more stable numerically, and more accurate methods for computing tangential derivatives along material interfaces and for matching boundary conditions with the homogeneous top and bottom regions. Numerical examples indicate that the improved BIE-NtD map method achieves a high order of accuracy for in-plane and conical diffractions of dielectric gratings.

Boundary integral equation Neumann-to-Dirichlet map method for gratings in conical diffraction

Boundary integral equation methods for diffraction gratings are particularly suitable for gratings with complicated material interfaces but are difficult to implement due to the quasi-periodic Green's function and the singular integrals at the corners. In this paper, the boundary integral equation Neumann-to-Dirichlet map method for in-plane diffraction problems of gratings [Y. Wu and Y. Y. Lu, J. Opt. Soc. Am. A26, 2444 (2009)] is extended to conical diffraction problems. The method uses boundary integral equations to calculate the so-called Neumann-to-Dirichlet maps for homogeneous subdomains of the grating, so that the quasi-periodic Green's functions can be avoided. Since wave field components are coupled on material interfaces with the involvement of tangential derivatives, a least squares polynomial approximation technique is developed to evaluate tangential derivatives along these interfaces for conical diffraction problems. Numerical examples indicate that the method performs equally well for dielectric or metallic gratings.

Solving conical diffraction grating problems with integral equations

Off-plane scattering of time-harmonic plane waves by a plane diffraction grating with arbitrary conductivity and general surface profile is considered in a rigorous electromagnetic formulation. Integral equations for conical diffraction are obtained involving, besides the boundary integrals of the single and double layer potentials, singular integrals, the tangential derivative of single-layer potentials. We derive an explicit formula for the calculation of the absorption in conical diffraction. Some rules that are expedient for the numerical implementation of the theory are presented. The efficiencies and polarization angles compared with those obtained by Lifeng Li for transmission and reflection gratings are in a good agreement. The code developed and tested is found to be accurate and efficient for solving off-plane diffraction problems including high-conductive gratings, surfaces with edges, real profiles, and gratings working at short wavelengths.

Analyzing diffraction gratings by a boundary integral equation Neumann-to-Dirichlet map method

For analyzing diffraction gratings, a new method is developed based on dividing one period of the grating into homogeneous subdomains and computing the Neumann-to-Dirichlet (NtD) maps for these subdomains by boundary integral equations. For a subdomain, the NtD operator maps the normal derivative of the wave field to the wave field on its boundary. The integral operators used in this method are simple to approximate, since they involve only the standard Green's function of the Helmholtz equation in homogeneous media. The method retains the advantages of existing boundary integral equation methods for diffraction gratings but avoids the quasi-periodic Green's functions that are expensive to evaluate.

Chebyshev collocation Dirichlet-to-Neumann map method for diffraction gratings

For diffraction gratings with layered refractive index profiles, the Fourier modal method is widely used. However, it is quite expensive to calculate the eigenmodes for each layer, especially when the structure involves absorptive media. We develop an efficient method that avoids the eigenvalue problems based on the so-called Dirichlet-to-Neumann (DtN) map. For each layer, the DtN map is an operator that maps the wave field to its normal derivative on one period of the boundaries of the layer, and it is approximated by a matrix. An efficient procedure for computing the DtN map is developed based on a Chebyshev collocation method and a fourth-order finite difference method for discretizing the uniform and the periodic directions, respectively. The efficiency and accuracy of our method are illustrated by numerical examples.

Adaptive finite-element method for diffraction gratings

A second-order finite-element adaptive strategy with error control for one-dimensional grating problems is developed. The unbounded computational domain is truncated to a bounded one by a perfectly-matched-layer (PML) technique. The PML parameters, such as the thickness of the layer and the medium properties, are determined through sharp a posteriori error estimates. The adaptive finite-element method is expected to increase significantly the accuracy and efficiency of the discretization as well as reduce the computation cost. Numerical experiments are included to illustrate the competitiveness of the proposed adaptive method.

Differential method applied for photonic crystals

The classic differential method is applied for modeling the diffraction of light from two-dimensional photonic crystals that consist of dielectric cylindrical objects. Special attention is paid to mutual interpenetration of consecutive layers. Two algorithms for dealing with a stack of repetitive layers are discussed, namely, the eigenvalue technique and the S-matrix algorithm. Their advantages and limitations are analyzed, and times required for their implementation are compared.

Numerical performance of finite-difference modal methods for the electromagnetic analysis of one-dimensional lamellar gratings

The numerical performance of a finite-difference modal method for the analysis of one-dimensional lamellar gratings in a classical mounting is studied. The method is simple and relies on first-order finite difference in the grating to solve the Maxwell differential equations. The finite-difference scheme incorporates three features that accelerate the convergence performance of the method: (1) The discrete permittivity is interpolated at the lamellar boundaries, (2) mesh points are located on the permittivity discontinuities, and (3) a nonuniform sampling with increased resolution is performed near the discontinuities. Although the performance achieved with the present method remains inferior to that achieved with up-to-date grating theories such as rigorous coupled-wave analysis with adaptive spatial resolution, it is found that the present method offers rather good performance for metallic gratings operating in the visible and near-infrared regions of the spectrum, especially for TM polarization.