This paper is concerned with the electromagnetic theory of characteristic modes for general homogeneous material bodies using surface integral equations (SIEs), which are computationally more efficient than using volume integral equations (VIEs). However, the generalized eigenvalue equations (GEEs) based on SIEs are prone to generate spurious or extra modes, unlike their counterparts based on VIEs. In this work, several variant SIE-based GEEs are proposed to eliminate the spurious/extra modes and improve the numerical performances. These variant GEEs share a common point that the reactive power operator related to the material domain should be removed from the left-hand sides of GEEs by explicitly enforcing the boundary condition, because it seems sensitive to incur unwanted modes. These variant GEEs are diverse in choosing the active power operator related to the material domain on the right-hand sides to suit different medium parameters, which even can be dropped directly for lossless cases. Numerical results show that these variant GEEs do not generate spurious/extra resonant modes in observed frequency range, and at least one of the variants would have acceptable numerical accuracy.INDEX TERMS Characteristic modes, material bodies, power operators, spurious modes.