Abstract:There is a growing interest to use plasmonic nanostructures for many applications including clean energy (e.g., solar cells), biological and chemical sensing. The optical behavior of metals is completely different from that observed at low frequencies. While at RF, metal behaves like a perfect electric conductor (with a negligible surface impedance) which is modeled by a surface electric current, at optical frequencies, the penetration of the waves cannot be ignored. Capturing this effect requires volumetric discretization of the structure. A plasmonic material has a complex permittivity with a negative real and positive imaginary parts (assuming e <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">‒iwt</sup> time convention). This results in a penetration depth of 1/k <inf xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">0</inf> n <inf xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">i</inf> , where k <inf xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">0</inf> is the wavevector of the background material and n <inf xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">i</inf> is the imaginary part of the refractive index of the plasmonic object. In optical wavelengths where the wavevector is very large, penetration depth becomes very small. Consequently, to model and consider the evanescent waves inside the plasmonic structures, a numerical method with a very fine meshing strategy near the plasmon-dielectric interface is required. Existing modeling numerical methods can be classified into differential equation methods (e.g., Finite Element Method (FEM) and Finite Difference Time Domain (FDTD)) and integral equation (IE) methods (e.g., method of moments or MoM). In differential equation methods, the whole space should be discretized which will result in a large number of unknowns. In contrast, the integral equation methods only discretize the desired objects. Hence, if there is a large white space between objects, integral equation methods are preferred over differential equation methods. However, as soon as the size of the problem increases, and consequently the number of unknowns rises, traditional IE methods such as MoM will consume a huge memory and calculation time increases drastically. Also, these methods have serious problem dealing with multiscale multiphysics problems. Equivalence principle algorithm (EPA), as an efficient domain decomposition method, invokes equivalence principle to decompose the solution domain into subdomains and then solve each of them independently. Later, it stitches them together to calculate the solution for the original problem. The entire procedure can be accelerated by application of fast techniques such as the multilevel fast multipole algorithm (MLFMA). It provides a level of parallelization and also alleviates the need to build and save large matrices. Therefore, EPA (relatively) needs lesser amount of memory and it provides more robust solution for multiscale multiphysics problems. Although EPA has been used for several applications, but its application for modeling and analyzing nanoplasmonic structures has not been reported before and it will be presented in here.

Simplified numerical modeling for Fano interference-induced asymmetric light reflectance effect using equivalent medium theory

Finite-Element-Based Generalized Impedance Boundary Condition for Modeling Plasmonic Nanostructures

Study on Spontaneous Emission in Complex Multilayered Plasmonic System Via Surface Integral Equation Approach with Layered Medium Green'S Function

Modeling and Simulation of Active Plasmonics with the FDTD Method by Using Solid State and Lorentz–Drude Dispersive Model

Multiphysics modeling of plasmonic nanodevices

Fano Response Induced by the Interference Between Localized Plasmons and Interface Reflections in Metal-Insulator-metal Waveguide Structure

An Efficient Numerical Method To Analyze Nanoplasmonic Structures

Investigation on Plasmonic Responses in Multilayered Nanospheres Including Asymmetry and Spatial Nonlocal Effects

Time-domain Simulation for Plasmonic Devices

Simulation of Plasmonics Nanodevices with Coupled Maxwell and Schrödinger Equations Using the FDTD Method

Surface Plasmons Based on Nonlocal and Size-Dependent Effects of Metallic Nanoparticles

Asymmetric light reflectance by Fano resonance between Fresnel reflection and localized surface plasmons

A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects

A Boundary Integral Equation Formalism for Modeling Multiple Scattering of Light from 3D Nanoparticles Incorporating Nonlocal Effects

Theory, Figures of Merit, and Design Recipe of the Plasmonic Structure Composed of a Nano-Slit Aperture Surrounded by Surface Corrugations

Modeling photomolecular effect using generalized boundary conditions for Maxwell equations

Fano Resonances Based on Nanoscale Plasmonic Structure

Tailoring the Coupling Between Localized and Propagating Surface Plasmons: Realizing Fano-Like Interference and High-Performance Sensor

Multiple Fano Interferences in A Plasmonic Metamolecule Consisting of Asymmetric Metallic Nanodimers

An Improved FDTD Method to Calculate Nonlocal Response in Plasmonics

Reconfigurable Lateral Optical Force Achieved by Selectively Exciting Plasmonic Dark Modes Near Fano Resonance