Numerical modeling of active plasmonic metamaterials

The paper addresses numerical time-domain methods for modeling of active and passive dispersive media, needed for simulations of plasmonic metamaterials. The proposed algorithms differ from published results, as our models employ more general formalisms and are more computationally efficient. The frequency dispersion of the permittivity is considered as an arbitrary Pade approximant, its numerical implementation is more universal and effective for all known ADE and RC methods. The gain model is implemented for an arbitrary topology of transitions with the ADE method. The proposed dispersion models are in a good fit with spectroscopic data and are included into a database of optical materials at nanohub.org.