Search

In this paper, we consider numerical simulations of the nonlocal optical response of metallic nanostructure arrays inside a dielectric host, which is of particular interest to the nanoplasmonics community due to many unusual properties and potential applications. Mathematically, it is described by Maxwell's equations with discontinuous coefficients coupled with a set of Helmholtz-type equations defined only on the domains of metallic nanostructures. To solve this challenging problem, we develop an efficient multiscale method consisting of three steps. First, we extend the system into the domain occupied by the dielectric medium in a novel way and result in a coupled system with rapidly oscillating coefficients. A rigorous analysis of the error between the solutions of the original system and the extended system is given. Second, we derive the homogenized system and define the multiscale approximate solution for the extended system by using the multiscale asymptotic method. Third, to fix the inaccuracy of the multiscale asymptotic method inside the metallic nanostructures, we solve the original system in each metallic nanostructure separately with boundary conditions given by the multiscale approximate solution. A fast algorithm based on the $LU$ decomposition is proposed for solving the resulting linear systems. By applying the multiscale method, we obtain the results that are in good agreement with those obtained by solving the original system directly at a much lower computational cost. Numerical examples are provided to validate the efficiency and accuracy of the proposed method.

In this paper, we consider the frequency-domain Maxwell's equations coupled to a nonlocal Drude model which describes the nonlocal optical response in metallic nanostructures. We prove the existence and uniqueness of weak solutions to the coupled equations. A Galerkin finite element method based on the Raviart--Thomas and N\'{e}d\'{e}lec elements is proposed to solve the equations and the associated error estimates are given. This is the first work on the mathematical and numerical analysis of this model. Numerical examples are presented to verify our theoretical analysis.