Recently, it has received a great deal of attention to analyze the electromagnetic wave problems in dispersive media by using the finite-difference time-domain (FDTD) method. Accordingly, it is of great importance to employ a proper dispersion model which can fit the frequency-dependent permittivity of a medium considered. The reported dispersion models include Debye, Drude, Lorentz, modified Lorentz, quadratic complex rational function, complex-conjugate pole-residue (CCPR) models. The CCPR dispersion model has advantage over other dispersion models in the fact that accurate CCPR dispersion parameters can be simply extracted by using the powerful and robust vector fitting tool which has been widely used in the circuit theory. However, the arithmetic operation of CCPR-based FDTD implementation is involved with complex-valued numbers and thus its numerical computation is not efficient. In this work, we propose an accurate and efficient FDTD simulation for complex dispersive media. In specific, an accurate CCPR dispersion model is simply obtained using the vector fitting tool and then the CCPR dispersion model is converted to the modified Lorentz dispersion model which leads to the arithmetic operation of only realvalued numbers in its FDTD implementation. Numerical examples are used to illustrate the accuracy and efficiency of our dispersive FDTD simulation. INDEX TERMS Dispersion model, dispersive media, finite-difference time-domain (FDTD) method, human tissue, plasmonics.