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

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

Для розрахунку енергетичної зонної структури кристалів останнім часом набув поширення метод апріорних атомних псевдопотенціалів (ПП). У загальних рисах цей підхід ґрунтується на самоузгодженому пошуку ПП у наближенні функціонала локальної спінової густини. Стартова точка цієї процедури базується на релятивіському рівнянні Дірака для хвильової функції Gl (r) і Fl (r):

dF_l (r)/dr - (g/r)^. F_l (r) + a^. [E_l - V(r)]G_l (r) = 0 (1a)

dG_l (r)/dr + (g/r)^. G_l (r) - a[(2/a² ) + E_l - V(r)]^. F_l (r) = 0, (1b)

де l = a^-1 =137.07 – обернене значення константи надтонкої структури; g - ненульове ціле число. Розв’язки рівняння (1) визначають густину заряду:

r(r) = å [|G_l (r)|² + | F_l (r)|² ]. (2)

E_l < E_F

Використовуючи процедуру нелінійної інтерполяції, визначимо псевдопотенціал:

V_ps ^{(l )} (r) = V_ост (r) + V_ioн ^{(l )} (r)+ V_со ^{(l )} (r), (3)

де V_ост (r) – остовний потенціал, V^(l) _ioн (r) визначає іонну корекцію псевдопотенціалу і V^(l) _со (r) – спін-орбітальна корекція. Кожний з вищенаведених доданків може бути виражений в аналітичній формі і тому відповідні матричні елементи можна точно обчислити (замість числового інтегрування):

V_îńň (r) = (-Z_v ^. e² /r)^. C_i ^. erf[a_i ^1/2. r]; (4)

V_łîí (r) = (A_i +r^2. A_i+3 )^. exp(-g_i ^. r² ). (5)

Інтерполяційні коефіціенти A_i , a_i , C_I ,g_i визначаються як розв`язки самоузгодженого рівняння Дірака – Хартрі – Фока – Слетера для конкретних атомів з відповідними орбітальними числами з подальшою нелінійною інтерполяцією (4, 5).

Повний псевдопотенціал є сумою нелокальних псевдопотенціалів, які перекриваються і розміщені в точках t _p,q . Взаємодія електрона з остовами, які описуються періодичними іонними ПП, визначається оператором:

V_ps (r , r ') = å V_q,s (r - R _p - t _q,s , t - R _p - t _q,s ), (6)

p,q,s

де R _p – вектор прямої ґратки, який визначає розташування елементарної комірки; t _q,s – вектор, який визначає розташування s -го іона сорту q в елементарній комірці. При обчисленні матричних елементів секулярного рівняння в базисі плоских хвиль необхідно перейти від r -простору до оберненого G -простору за допомогою Фур’є перетворення:

<k +G _i |V_ps (r , r^’ )|k +G _j >=N^-1 exp[-iq (R _p + t _q,s )]^. V^-1.

ò dr ^. dr^{’ .} exp[-i(k +G _j ) ^. r ]^. V_q (r , r ')^. exp[i(k +G _j ) ^. r '], (7)

де q = G _i - G _j і V – об’єм першої зони Бріллюена. Враховуючи, що q _j ^. t _i = 2^. p^. d_ij (t _i , q _j – основні вектори прямої й оберненої ґратки) знаходимо, що при довільних значеннях p exp(-iqR p) = 1, а тому сума по p просто дає множник N i остання рівність у локальному наближенні набуває вигляду:

<k +G _i |V_ps (r , r^’ )|k +G _j >= å exp[-iq ^. t _q,s )]^. V^-1. ò d³ r ^. exp[-i(q ^. r ]^. V_q (r )(8)

q, s

Зауважимо, що тут інтегрування по r здійснюється в основній сфері кристала.

Розглянемо локальну частину псевдопотенціала. Форм-фактор потенціалу іона (4) дорівнює:

V_q = V^-1. ò V_îńň (r )^. exp(-iqr )^. d³ r.

Для обчислення цього виразу використаємо загальну процедуру. Запишемо розклад плоских хвиль за Релеєм:

exp(iqr ) = (2l +1)^. i^{l .} j_l (|q ^. r |^. P

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать