Numerical Simulation of Nonlinear Effectsin Multilayer Graphene Metasurfacesin Terahertz and Infrared Ranges

Using the developed numerical-analytical method the nonlinear effects in multilayer metasurfaces of graphene ribbons with finite length was investigated numerically at the resonant frequencies of the surface plasmon–polariton (SPP) modes in the terahertz (THz) and infrared (IR) frequency ranges. It is shown that the efficiency of nonlinear effects (the frequency multiplication) can be enhanced by increasing the number of array layers, its packing density, dielectric mirrors in substrates and tuned by changing the chemical potential of graphene (applying an external electric field).