“…The conformal FDTD method on non-orthogonal grids [15][16][17], counter path FDTD method [18,19], and subgridding method [20][21][22] can represent curved interfaces suitably, but they are relatively difficult to implement and increase the memory and computation time. A different approach using effective permittivities (EPs), which derives from interface interpolations based on Ampere's and Faraday's integration laws, can reduce the error of the permittivity model on coarse grids in a simple implementation and at low computational cost [23][24][25][26][27][28][29][30][31][32][33][34][35]. The EP method, however, is not effective for plasmonic materials because the SP resonance condition is changed on interfaces due to interpolated values of EPs.…”