Суперпозиція ЛКАО і псевдопотенціалу для розрахунку структури монокристалів CoW2

після m-го ітераційного кроку. Досягнута точність, менша ніж = 0,07% між вхідними і вихідними параметрами ітераційного кроку, служила головним критерієм самоузгодження. Власні значення енергії були стабільними потужністю до 0,003 атомних одиниць енергії. Процедура діагоналізації гамільтоніана здійснювалася QL методом [7].

Для підвищення точності опису густини електронних станів діагоналізація гамільтоніану здійснювалася для 64 рівновіддалених точок в 1/16-ій частині ЗБ. Числові обчислення проводилися з використанням методу теадраедрів. Вираз для інтегралів перекриття, які використовувалися при обчисленні:

= /4p ijDBy (32)

= /4p ijDBz (33)

= /4p ijDAx (34)

= /4p ijDAxy (35)

= /4p ijDAz [DAxDBx + 1/2(ai + bj )] (36)

= 3/4p ijDAz [DAxDBx + 1/2(ai + bj )] (37)

= 3/4p ijDAxDBy (38)

= 3/4p ijDAxDBz (39)

= 3/4p ijDAyDBx (40)

= 3/4p ij [DAyDBy + 1/2(ai + bj )] (41)

= 3/4p ijDAyDBz (42)

= 3/4p ijDAzDBx (43)

= 3/4p ijDAzDBy (44)

= 3/4p ij [DAzDBz + 1/2(ai + bj )], (45)

де DBx = Dx - Bx ł DAx = Dx - Ax.

Після подібної процедури як і у випадку плоских хвиль, можна отримати кінцеві вирази для матричних елементів секулярного рівняння:

= [ Ea - E(k) ] [ Sqq''+ exp(-iktq ).Sqq' +

Bqq' + exp(-iktq )Aqq', (47)

де

Sqq' = < n | n', r + l - l' > ; (48a)

Bqq' = < n | V(r ) - Va(r) | n' > ; (48b)