Abstract

Accurate numerical solutions for electromagnetic fields near sharp corners and edges are important for nanophotonics applications that rely on strong near fields to enhance light-matter interactions. For cylindrical structures, the singularity exponents of electromagnetic fields near sharp edges can be solved analytically, but in general the actual fields can only be calculated numerically. In this paper, we use a boundary integral equation method to compute electromagnetic fields near sharp edges, and construct the leading terms in asymptotic expansions based on numerical solutions. Our integral equations are formulated for rescaled unknown functions to avoid unbounded field components, and are discretized with a graded mesh and properly chosen quadrature schemes. The numerically found singularity exponents agree well with the exact values in all the test cases presented here, indicating that the numerical solutions are accurate.

© 2017 Optical Society of America

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

Yumao Wu and Ya Yan Lu
J. Opt. Soc. Am. A 28(6) 1191-1196 (2011)

High order integral equation method for diffraction gratings

Wangtao Lu and Ya Yan Lu
J. Opt. Soc. Am. A 29(5) 734-740 (2012)

New boundary integral equations for computer-aided design of dielectric waveguide circuits

Kazuo Tanaka and Masaaki Kojima
J. Opt. Soc. Am. A 6(5) 667-674 (1989)

Accurate simulation of two-dimensional optical microcavities with uniquely solvable boundary integral equations and trigonometric Galerkin discretization

Svetlana V. Boriskina, Phillip Sewell, Trevor M. Benson, and Alexander I. Nosich
J. Opt. Soc. Am. A 21(3) 393-402 (2004)

Analyzing bowtie structures with sharp tips by a vertical mode expansion method

Hualiang Shi, Ya Yan Lu, and Qiang Du
Opt. Express 26(24) 32346-32352 (2018)