Dipole excitation of surface plasmon on a conducting sheet: Finite element approximation and validation

Abstract

We formulate and validate a finite element approach to the propagation of a slowly decaying electromagnetic wave, called surface plasmon-polariton, excited along a conducting sheet, e.g., a single-layer graphene sheet, by an electric Hertzian dipole. By using a suitably rescaled form of time-harmonic Maxwell's equations, we derive a variational formulation that enables a direct numerical treatment of the associated class of boundary value problems by appropriate curl-conforming finite elements. The conducting sheet is modeled as an idealized hypersurface with an effective electric conductivity. The requisite weak discontinuity for the tangential magnetic field across the hypersurface can be incorporated naturally into the variational formulation. We carry out numerical simulations for an infinite sheet with constant isotropic conductivity embedded in two spatial dimensions; and validate our numerics against the closed-form exact solution obtained by the Fourier transform in the tangential coordinate. Numerical aspects of our treatment such as an absorbing perfectly matched layer, as well as local refinement and a posteriori error control are discussed.Keywords: Time-harmonic Maxwell's equations, finite element method, surface plasmon-polariton, weak discontinuity on hypersurface 2010 MSC: 65N30, 78M10, 78M30, 78A45

Adaptive finite element simulations of waveguide configurations involving parallel 2D material sheets

“…For local refinement, the fit parameters we obtain are a local = 0.205333, b local = 1.25229×10 8 , c local = 1.57811, whereas for uniform refinement we get a unif = 0.206216, b unif = 85.8656, c unif = 0.612326. As expected [27], we obtain a much faster convergence rate (c ≈ 1.6) in the quantity of interest for local refinement as opposed to uniform refinement (c ≈ 0.6). We conclude that our computation of the energy transmission ratio with J (E) ≈ a local is reliable within 1%.…”

“…Next, we preform a parameter study of the energy transmition ratio for varying interlayer spacing d ranging from d min = 4 · 10 −12 to d max = 0.2, where d max corresponds to about 1/3 of the single-layer SPP wavelength, 2π/Re(k m,r ), or 1/30 of the free space wavelength [27]. The (interpolated) results are shown in Figure 5.…”

On the Wiener–Hopf Method for Surface Plasmons: Diffraction from Semiinfinite Metamaterial Sheet

“…Equations (16) then entail ± (ξ ) → 0 as ξ → ∞ in C ± . Thus, by (18) we assert that E(ξ ) → 0 as ξ → ∞ in C. By Liouville's theorem, E(ξ ) must be a constant everywhere, which has to be zero; E(ξ ) ≡ 0. Consequently, formula (3) results from Fourier-inversion of (19) by use of (13) and (16).…”

“…Next, we perform a finite element computation for the electric field in the presence of the semiinfinite sheet in 2D [18]. The computational domain is chosen to be large enough to account for distances, x, from the edge such that 0 ≤ k x ≤ 30.…”

Homogenization of plasmonic crystals: seeking the epsilon-near-zero effect

“…Building on computational methods that we developed for plasmonic problems [20,21], we propose the following computational approach for present purposes:…”