Finite difference time domain (FDTD) simulations were performed with a home made workstation with Dual CPU (AMD EPYC 7302 16 core × 2), 512 GB of DDR4 memories. For the spectrally refined modeling, a high accuracy simulations were performed to reveal detailed spectral features. Commercial software FDTD solution (ANSYS/Lumerical) was used for numerical modeling. For a wide wavelength region from visible to MIR, the Palik's optical permittivity data were used for Au, Cr, SiO 2 , Al 2 O 3 . For the ITO properties which were difficult to measure experimentally over entire vis-MIR spectral range, we setup a model based on the 3D conductive and plasma based on Drude models described next. For the Drude's free electron model in two equivalent presentations: