A general algorithm for computing the bidimensional spectral Green's dyad in multilayered complex bianisotropic media: the equivalent boundary method

A systematic method to obtain the bidimensional spectral dyadic Green's function of stratified planar structures with arbitrary complex bianisotropic layers is developed. The method is based on the uniqueness and equivalence electromagnetic theorems. A first-order partial differential formulation for the electromagnetic field inside each layer is used. An explicit algorithm makes it possible to go from the single-layer formulas to the general multilayer matrix formulation. The perturbative nature of the method provides good numerical efficiency and straightforward determination of asymptotic behavior.>