Устранение полюсов детерминанта ККР в расчетах зонной структуры

Детерминант секулярного уравнения метода Кориги, Кона и Ростокера (ККР) для расчета зонной структуры имеет полюсы при $E= |k+G_n|^2$, где $G_n$ - векторы обратной решетки кристалла. Показано, что можно получить уравнение с матрицей, элементы и детерминант которой являются непрерывными функциями энергии на выбранном интервале. Это позволяет устранить трудности вычисления спектра вблизи полюсов и автоматизировать вычисление спектра во всем интервале энергий для данного волнового вектора $k$.