Comparative Analysis of Analytic and Numerical Versions of the Laplacian Representation Approach for the Asymptotic Part of the Layered-Medium Green’s Function in the Mixed Potential Formulation

We report on developments in the evaluation of matrix elements of the electric and magnetic field operators involving the asymptotic (large transverse wave-number or small transverse distances) components of the mixedpotential Green’s function of a layered medium. Subtracting these asymptotic terms significantly accelerates numerical computation of the Sommerfeld-type integrals required in constructing Green’s function [1]. The described method is applicable to planar conducting structures embedded in a multi-layered medium. Ability of a fast and accurate evaluation of matrix elements for such structures is critical in many applications, especially in solving forward problems in design and optimization of electromagnetic metasurfaces.