Приближение локальной плотности - Local-density approximation

Приближения локальной плотности (LDA) представляют собой класс приближений к обмен –корреляция (XC) энергия функциональный в теория функционала плотности (DFT), которые зависят исключительно от значения электронная плотность в каждой точке пространства (а не, например, производные плотности или Орбитали Кона – Шама ). Многие подходы могут дать локальные приближения к энергии XC. Однако наиболее успешными локальными приближениями являются те, которые были получены из однородный электронный газ (HEG) модель. В этом отношении LDA обычно является синонимом функционалов, основанных на приближении HEG, которые затем применяются к реалистичным системам (молекулам и твердым телам).

В общем, для спин-неполяризованной системы приближение локальной плотности для обменно-корреляционной энергии записывается как

${ Displaystyle E_ {xc} ^ { mathrm {LDA}} [ rho] = int rho ( mathbf {r}) epsilon _ {xc} ( rho ( mathbf {r})) mathrm {d} mathbf {r} ,}$

куда ρ это электронная плотность и ε_xc - обменно-корреляционная энергия, приходящаяся на одну частицу однородный электронный газ плотности заряда ρ. Обменно-корреляционная энергия линейно разлагается на обменные и корреляционные члены:

${ Displaystyle E_ {xc} = E_ {x} + E_ {c} ,}$

так что отдельные выражения для E_Икс и E_c ищутся. Член обмена принимает простую аналитическую форму для HEG. Точно известны только предельные выражения для корреляционной плотности, что приводит к множеству различных приближений для ε_c.

Приближения локальной плотности важны при построении более сложных приближений к обменно-корреляционной энергии, таких как обобщенные градиентные приближения (GGA) или гибридные функционалы, поскольку желательным свойством любого приближенного обменно-корреляционного функционала является то, что он воспроизводит точные результаты HEG для неизменяющихся плотностей. По сути, LDA часто являются явным компонентом таких функционалов.

Приложения

Приближения локальной плотности, как и в случае с GGA, широко используются физики твердого тела в ab-initio исследованиях DFT для интерпретации электронных и магнитных взаимодействий в полупроводниковых материалах, включая полупроводниковые оксиды и спинтроника. Важность этих вычислительных исследований проистекает из сложности системы, которая приводит к высокой чувствительности к параметрам синтеза, что требует анализа на основе первых принципов. Предсказание Уровень Ферми и зонная структура в легированных полупроводниковых оксидах часто выполняется с использованием LDA, включенного в пакеты моделирования, такие как CASTEP и DMol3.^[1] Однако недооценка Ширина запрещенной зоны ценности, часто связанные с LDA и GGA приближения могут привести к ложным предсказаниям проводимости, опосредованной примесями, и / или опосредованного носителями магнетизма в таких системах.^[2] Начиная с 1998 г., применение Теорема Рэлея для собственных значений привело к в основном точным расчетам запрещенной зоны материалов с использованием потенциалов LDA.^[3][4] Непонимание второй теоремы DFT, по-видимому, объясняет большую часть недооценки ширины запрещенной зоны расчетами LDA и GGA, как объясняется в описании теория функционала плотности, в связи с утверждениями двух теорем ДПФ.

Однородный электронный газ

Приближение для ε_xc в зависимости только от плотности может развиваться множеством способов. Наиболее успешный подход основан на однородном электронном газе. Это построено путем размещения N взаимодействующие электроны в объеме, V, с положительным фоновым зарядом, сохраняющим нейтральность системы. N и V затем уносятся на бесконечность таким образом, чтобы сохранялась плотность (ρ = N / V) конечно. Это полезное приближение, поскольку полная энергия состоит только из вкладов кинетической энергии и обменно-корреляционной энергии, и что волновая функция выражается в терминах плоских волн. В частности, при постоянной плотности ρ, плотность обменной энергии пропорциональна ρ^⅓.

Функционал обмена

Плотность обменной энергии ГЭГ известна аналитически. LDA для обмена использует это выражение в приближении, что обменная энергия в системе с неоднородной плотностью получается путем поточечного применения результатов HEG, в результате чего получается выражение^[5][6]

${ displaystyle E_ {x} ^ { mathrm {LDA}} [ rho] = - { frac {3} {4}} left ({ frac {3} { pi}} right) ^ { 1/3} int rho ( mathbf {r}) ^ {4/3} mathrm {d} mathbf {r} .}$

Функционал корреляции

Аналитические выражения для корреляционной энергии ГЭГ доступны в пределах высокой и низкой плотности, соответствующих бесконечно-слабой и бесконечно сильной корреляции. Для ГЭГ с плотностью ρ, предел высокой плотности корреляционной плотности энергии равен^[5]

${ displaystyle epsilon _ {c} = A ln (r_ {s}) + B + r_ {s} (C ln (r_ {s}) + D) ,}$

и нижний предел

${ displaystyle epsilon _ {c} = { frac {1} {2}} left ({ frac {g_ {0}} {r_ {s}}} + { frac {g_ {1}} { r_ {s} ^ {3/2}}} + точки right) ,}$

где параметр Вигнера-Зейтца ${ displaystyle r_ {s}}$ безразмерен.^[7] Он определяется как радиус сферы, охватывающей ровно один электрон, деленный на радиус Бора. Параметр Вигнера-Зейтца ${ displaystyle r_ {s}}$ связана с плотностью как

${ displaystyle { frac {4} {3}} pi r_ {s} ^ {3} = { frac {1} { rho}} .}$

На основе теории возмущений многих тел предложено аналитическое выражение для всего диапазона плотностей. Рассчитанные корреляционные энергии согласуются с результатами квантовый Монте-Карло симуляция с точностью до 2 миллихартри.

• Функционал корреляции Чачиё

${ displaystyle epsilon _ {c} = a ln left (1 + { frac {b} {r_ {s}}} + { frac {b} {r_ {s} ^ {2}}} верно).}$ ^[8]

Параметры ${ displaystyle a}$ и ${ displaystyle b}$ не от эмпирической подгонки к данным Монте-Карло, но из теоретического ограничения, что функционал приближается к пределу высокой плотности. Формула Чачиё более точна, чем стандартная функция соответствия VWN.^[9] в атомная единица, ${ Displaystyle а = { гидроразрыва { ln (2) -1} {2 pi ^ {2}}}}$. Выражение в закрытой форме для ${ displaystyle b}$ существует; но удобнее использовать числовое значение: ${ displaystyle b = 20.4562557 = exp ({ text {C}} / 2a)}$. Здесь, ${ displaystyle { text {C}}}$ был точно вычислен с использованием интеграла в замкнутой форме и дзета-функции (уравнение 21, Г. Хоффман, 1992).^[10] ${ Displaystyle { текст {C}} = { tfrac { ln (2)} {3}} - { tfrac {3} {2 pi ^ {2}}} left [ zeta (3) + { tfrac {22} {9}} - { tfrac { pi ^ {2}} {3}} + { tfrac {32 ln (2)} {9}} - { tfrac {8 ln ^ {2} (2)} {3}} right] + { tfrac {2 (1- ln 2)} { pi ^ {2}}} left [ ln ({ tfrac {4 } { alpha pi}}) + left langle ln R_ {0} right rangle _ { text {av}} - { tfrac {1} {2}} right].}$ Сохраняя прежнюю функциональную форму,^[11] параметр ${ displaystyle b}$ также был адаптирован к моделированию Монте-Карло, обеспечивая лучшее согласование. Также в этом случае ${ displaystyle r_ {s}}$ должно быть либо в атомных единицах, либо делиться на радиус Бора, что делает его безразмерным параметром.^[7]

Таким образом, формула Чачиё представляет собой простой (также точный) функционал корреляции из первого принципа для DFT (однородной электронной плотности). Тесты на кривых дисперсии фононов ^[12] дают достаточную точность по сравнению с экспериментальными данными. Его простота также подходит для вводных курсов теории функционала плотности.^[13][14]

Сравнение нескольких функционалов корреляционной энергии LDA и квантового моделирования Монте-Карло

Точный квантовый Монте-Карло моделирование энергии ГЭГ было выполнено для нескольких промежуточных значений плотности, что, в свою очередь, позволило получить точные значения плотности корреляционной энергии.^[15] Наиболее популярные LDA для корреляционной плотности энергии интерполируют эти точные значения, полученные в результате моделирования, воспроизводя точно известное предельное поведение. Различные подходы, использующие разные аналитические формы для ε_c, создали несколько LDA для корреляционного функционала, в том числе

• Воско-Вилк-Нусаир (VWN) ^[16]
• Perdew-Zunger (PZ81) ^[17]
• Коул-Пердью (CP) ^[18]
• Пердью-Ван (PW92) ^[19]

Предваряя их и даже формальные основы самого ДПФ, был получен корреляционный функционал Вигнера. пертурбативно от модели HEG.^[20]

Спиновая поляризация

Распространение функционалов плотности на спин-поляризованный системы просты для обмена, где известен точный спин-скейлинг, но для корреляции необходимо использовать дополнительные приближения. Спин-поляризованная система в DFT использует две спиновые плотности: ρ_α и ρ_β с ρ = ρ_α + ρ_β, а форма приближения локальной спиновой плотности (LSDA) имеет вид

${ displaystyle E_ {xc} ^ { mathrm {LSDA}} [ rho _ { alpha}, rho _ { beta}] = int mathrm {d} mathbf {r} rho ( mathbf {r}) epsilon _ {xc} ( rho _ { alpha}, rho _ { beta}) .}$

Для обменной энергии точный результат (не только для приближений локальной плотности) известен в терминах спин-неполяризованного функционала:^[21]

${ displaystyle E_ {x} [ rho _ { alpha}, rho _ { beta}] = { frac {1} {2}} { bigg (} E_ {x} [2 rho _ { alpha}] + E_ {x} [2 rho _ { beta}] { bigg)} .}$

Спин-зависимость плотности корреляционной энергии приближается путем введения относительной спиновой поляризации:

${ displaystyle zeta ( mathbf {r}) = { frac { rho _ { alpha} ( mathbf {r}) - rho _ { beta} ( mathbf {r})} { rho _ { alpha} ( mathbf {r}) + rho _ { beta} ( mathbf {r})}} .}$

${ Displaystyle zeta = 0 ,}$ соответствует диамагнитной спин-неполяризованной ситуации с равными${ Displaystyle альфа ,}$ и ${ Displaystyle бета ,}$ спиновые плотности, тогда как ${ displaystyle zeta = pm 1}$ соответствует ферромагнитной ситуации, когда одна спиновая плотность обращается в нуль. Плотность спиновой корреляционной энергии для заданных значений полной плотности и относительной поляризации, ε_c(ρ,ς), построен так, чтобы интерполировать крайние значения. Несколько форм были разработаны совместно с корреляционными функционалами LDA.^[16][22]

Наглядные расчеты

Расчеты LDA разумно согласуются с экспериментальными значениями.

Обменно-корреляционный потенциал

Обменно-корреляционный потенциал, соответствующий обменно-корреляционной энергии для приближения локальной плотности, определяется выражением^[5]

${ displaystyle v_ {xc} ^ { mathrm {LDA}} ( mathbf {r}) = { frac { delta E ^ { mathrm {LDA}}} { delta rho ( mathbf {r}) )}} = epsilon _ {xc} ( rho ( mathbf {r})) + rho ( mathbf {r}) { frac { partial epsilon _ {xc} ( rho ( mathbf { r}))} { partial rho ( mathbf {r})}} .}$

В конечных системах потенциал LDA убывает асимптотически с экспоненциальной формой. Это ошибка; истинный обменно-корреляционный потенциал кулоновским образом спадает гораздо медленнее. Искусственно быстрый распад проявляется в количестве орбиталей Кона – Шэма, которые потенциал может связать (то есть, сколько орбиталей имеют энергию меньше нуля). Потенциал LDA не может поддерживать ряд Ридберга, и те состояния, которые он связывает, имеют слишком высокую энергию. Это приводит к HOMO энергия слишком высока, так что любые прогнозы для потенциал ионизации на основе Теорема Купманса бедны. Кроме того, LDA дает плохое описание богатых электронами разновидностей, таких как анионы где он часто не может связать дополнительный электрон, ошибочно предсказывая, что виды будут нестабильными.^[17][23]

Рекомендации