▎ 摘　　要

The accurate time-domain numerical modeling of graphene surface conductivity due to the interband contribution is presented in this article. Initially, the lower frequency limit for the interband conductivity inclusion is studied, highlighting that even in the far-infrared regime, it should not be ignored. Then, a precise vector-fitting technique is utilized to decompose the analytical conductivity function into complex conjugate pairs. Every pair is connected to a complex surface current, while the frequency dispersion is evaluated via a proper recursive convolution scheme. Finally, a straightforward algebraic manipulation is conducted to eliminate the complex terms from the updated equations of the electromagnetic components, concerning the finite-difference time-domain (FDTD) algorithm. The precision and efficiency of the proposed methodology are thoroughly validated by comparing the numerically extracted surface wave propagation characteristics to those obtained in terms of analytical expressions.