“…In the context of FDTD studies, several algorithms for wave propagation in media that exhibit both anisotropy and dispersive tensor elements have been proposed and implemented, indicatively, via recursive convolution (RC) [31], piecewise linear recursive convolution (PLRC) [32], z-transform [33,34], and auxiliary differential equation (ADE) [35] formulations. In these works, the permittivity tensors of the simulated materials, typically magnetized plasma and ferrites, exhibit dispersion types incapable of describing the frequency dependence of LC tensor elements and, hence, they are not suitable for the numerical investigation of LC-based photonic or plasmonic devices.…”