Calculate ionization energies and analyze binding properties of multi-level impurity centers (double donors, double acceptors) with higher Coulomb charge. Use when working with substitutional impurities that donate or accept more than one carrier, such as Group VI elements in silicon or germanium, where hydrogen-like approximation becomes less accurate due to tight binding.