WikiDer > Гибридный функционал

Гибридные функционалы являются классом приближений к обмен–корреляция энергия функциональный в теория функционала плотности (DFT), которые включают в себя часть точного обмена с Теория Хартри – Фока с остальной обменно-корреляционной энергией из других источников (ab initio или эмпирический). Точный функционал обменной энергии выражается через Орбитали Кона – Шама а не плотность, поэтому называется скрытый функционал плотности. Одна из наиболее часто используемых версий - B3LYP, что означает "Becke, 3-параметрический, Lee–Ян–Парр".

Источник

Гибридный подход к построению приближений функционала плотности был введен Аксель Беке в 1993 г.^[1] Гибридизация с обменом Хартри-Фока (HF) (также называемым точным обменом) обеспечивает простую схему для улучшения расчета многих молекулярных свойств, таких как энергии атомизации, длина облигаций и частоты колебаний, которые, как правило, плохо описываются простыми "ab initio" функционалами.^[2]

Метод

Гибридный обменно-корреляционный функционал обычно строится в виде линейная комбинация точного обменного функционала Хартри – Фока

${ displaystyle E _ { text {x}} ^ { text {HF}} = - { frac {1} {2}} sum _ {i, j} iint psi _ {i} ^ {* } ( mathbf {r} _ {1}) psi _ {j} ^ {*} ( mathbf {r} _ {2}) { frac {1} {r_ {12}}} psi _ { j} ( mathbf {r} _ {1}) psi _ {i} ( mathbf {r} _ {2}) , d mathbf {r} _ {1} , d mathbf {r} _ {2}}$

и любое количество явных функционалов плотности обмена и корреляции. Параметры, определяющие вес каждого отдельного функционала, обычно задаются путем подгонки предсказаний функционала к экспериментальным или точно рассчитанным термохимическим данным, хотя в случае «функционалов адиабатической связи» веса могут быть установлены априори.^[2]

B3LYP

Например, популярный B3LYP (Becke, 3-параметр, Lee – Yang – Parr)^[3][4] обменно-корреляционный функционал

${ displaystyle E _ { text {xc}} ^ { text {B3LYP}} = E _ { text {x}} ^ { text {LDA}} + a_ {0} (E _ { text {x}} ^ { text {HF}} - E _ { text {x}} ^ { text {LDA}}) + a _ { text {x}} (E _ { text {x}} ^ { text {GGA }} - E _ { text {x}} ^ { text {LDA}}) + E _ { text {c}} ^ { text {LDA}} + a _ { text {c}} (E _ { текст {c}} ^ { text {GGA}} - E _ { text {c}} ^ { text {LDA}}),}$

куда ${ displaystyle a_ {0} = 0.20}$, ${ displaystyle a _ { text {x}} = 0,72}$, и ${ displaystyle a _ { text {c}} = 0,81}$. ${ displaystyle E _ { text {x}} ^ { text {GGA}}}$ и ${ displaystyle E _ { text {c}} ^ { text {GGA}}}$ находятся обобщенные градиентные приближения: функционал обмена Becke 88^[5] и корреляционный функционал Ли, Янга и Парра^[6] для B3LYP и ${ displaystyle E _ { text {c}} ^ { text {LDA}}}$ это VWN приближение локальной плотности к функционалу корреляции.^[7]

Три параметра, определяющие B3LYP, были взяты без изменений из оригинальной подгонки Бекке аналогичного функционала B3PW91 к набору энергий атомизации, потенциалов ионизации, сродства протонов и полной атомной энергии.^[8]

PBE0

Функционал PBE0^[2][9]смешивает обменную энергию Пердью-Берк-Эрнцерхофа (PBE) и обменную энергию Хартри-Фока в установленном соотношении 3: 1 вместе с полной корреляционной энергией PBE:

${ displaystyle E _ { text {xc}} ^ { text {PBE0}} = { frac {1} {4}} E _ { text {x}} ^ { text {HF}} + { frac {3} {4}} E _ { text {x}} ^ { text {PBE}} + E _ { text {c}} ^ { text {PBE}},}$

куда ${ displaystyle E _ { text {x}} ^ { text {HF}}}$ - точный обменный функционал Хартри – Фока, ${ displaystyle E _ { text {x}} ^ { text {PBE}}}$ функционал обмена PBE, и ${ displaystyle E _ { text {c}} ^ { text {PBE}}}$ - функционал корреляции PBE.^[10]

HSE

НИУ ВШЭ (Heyd – Scuseria – Ernzerhof)^[11] обменно-корреляционный функционал использует функция-ошибка-экранированный Кулоновский потенциал для расчета обменной части энергии с целью повышения вычислительной эффективности, особенно для металлических систем:

${ displaystyle E _ { text {xc}} ^ { omega { text {PBEh}}} = aE _ { text {x}} ^ { text {HF, SR}} ( omega) + (1- а) E _ { text {x}} ^ { text {PBE, SR}} ( omega) + E _ { text {x}} ^ { text {PBE, LR}} ( omega) + E_ { text {c}} ^ { text {PBE}},}$

куда ${ displaystyle a}$ - параметр смешения, а ${ displaystyle omega}$ - настраиваемый параметр, контролирующий краткость взаимодействия. Стандартные значения ${ displaystyle a = 1/4}$ и ${ displaystyle omega = 0,2}$ (обычно обозначаемый как HSE06) показал хорошие результаты для большинства систем. Обменно-корреляционный функционал HSE вырождается в гибридный функционал PBE0 при ${ displaystyle omega = 0}$. ${ displaystyle E _ { text {x}} ^ { text {HF, SR}} ( omega)}$ - короткодействующий точный обменный функционал Хартри – Фока, ${ displaystyle E _ { text {x}} ^ { text {PBE, SR}} ( omega)}$ и ${ Displaystyle Е _ { текст {х}} ^ { текст {PBE, LR}} ( omega)}$ являются короткодействующими и дальнодействующими компонентами обменного функционала PBE, и ${ displaystyle E _ { text {c}} ^ { text {PBE}} ( omega)}$ это PBE^[10] корреляционный функционал.

Мета-гибрид GGA

Набор функционалов M06^[12][13] представляет собой набор из четырех мета-гибридных функционалов GGA и meta-GGA DFT. Эти функционалы строятся путем эмпирической подгонки их параметров, будучи ограниченными однородным электронным газом.

В семейство входят функционалы M06-L, M06, M06-2X и M06-HF, с разным количеством точного обмена для каждого из них. M06-L полностью локален без HF-обмена (поэтому его нельзя считать гибридным), M06 имеет 27% HF-обмен, M06-2X 54% и M06-HF 100%.

Преимущества и полезность каждого функционала

• M06-L: Быстрый, хороший для переходных металлов, неорганических и металлоорганических соединений.
• M06: Для основной группы, металлорганических соединений, кинетики и нековалентных связей.
• M06-2X: Основная группа, кинетика.
• M06-HF: TD-DFT с переносом заряда, системы, в которых самовзаимодействие является патологическим.

Набор дает хорошие результаты для систем, содержащих дисперсионные силы, что является одним из самых больших недостатков стандартных методов DFT. Коэффициенты масштабирования s6 для поправки на дисперсию на большие расстояния Гримма составляют 0,20, 0,25 и 0,06 для M06-L, M06 и M06-2X соответственно.

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