In recent work, electromagnetic propagation velocities for plane waves in dispersive metamaterials were calculated assuming frequency dispersion up to the second order. The three velocities were expressed in terms of dispersive coefficients under certain simplifying constraints. Frequency domains were found to exist around resonances where group and phase velocities are in opposition, implying possible negative index behavior.

In this paper, we incorporate in the derived equations physical models (including Debye, Lorentz and Condon) for material dispersion in permittivity, permeability and chirality in order to further examine the consequences of second-order dispersion leading to negative index for practical cases, and also evaluate the resulting phase and group indices.