Теория функционала плотности

Теория функционала плотности (англ. Density Functional Theory, DFT) — это метод квантово-механического моделирования, используемый для исследования электронной структуры многоэлектронных систем (атомов, молекул, кристаллов). В отличие от традиционных методов квантовой химии, основанных на волновой функции, DFT оперирует электронной плотностью как основной переменной, что позволяет существенно сократить вычислительные затраты при сохранении приемлемой точности для широкого класса задач. Теория базируется на теоремах Хоэнберга — Кона и уравнениях Кона — Шэма, сформулированных в 1964—1965 годах.

История

Предпосылки и ранние работы

Квантово-механическое описание многоэлектронных систем традиционно основывалось на уравнении Шрёдингера, решение которого для волновой функции Ψ(r₁, r₂, …, rₙ) становится экспоненциально сложным с ростом числа электронов N. Уже в 1927 году Ллевеллин Томас и Энрико Ферми независимо предложили модель, в которой кинетическая энергия электронов выражалась через их плотность, однако эта модель не учитывала обменные и корреляционные эффекты и была непригодна для точных расчётов.

Теоремы Хоэнберга — Кона (1964)

В 1964 году Пьер Хоэнберг и Уолтер Кон опубликовали работу, заложившую математические основы DFT. Они доказали две ключевые теоремы:

1. Первая теорема Хоэнберга — Кона: внешний потенциал v(r) (создаваемый ядрами) однозначно (с точностью до аддитивной константы) определяется электронной плотностью n(r) основного состояния. Следовательно, все свойства системы являются функционалами плотности.
2. Вторая теорема Хоэнберга — Кона: существует универсальный функционал энергии E[n], который достигает минимума при точной плотности основного состояния. Это позволило сформулировать вариационный принцип для плотности.

Эти теоремы доказали, что плотность n(r) содержит всю информацию о системе, но не давали практического способа вычисления функционала.

Уравнения Кона — Шэма (1965)

В 1965 году Уолтер Кон и Лю Шэм предложили метод, сделавший DFT вычислительно реализуемым. Они ввели вспомогательную систему невзаимодействующих электронов, движущихся в эффективном потенциале, который воспроизводит ту же плотность, что и реальная система. Кинетическая энергия невзаимодействующих электронов вычисляется точно, а остальные вклады (обмен и корреляция) объединяются в обменно-корреляционный функционал Eₓc[n]. Уравнения Кона — Шэма имеют вид одноэлектронных уравнений Шрёдингера, что позволяет решать их итерационно.

Развитие и признание

В 1970—1980-е годы были разработаны первые приближения для обменно-корреляционного функционала: приближение локальной плотности (LDA) и обобщённое градиентное приближение (GGA). В 1990-е годы появились гибридные функционалы (например, B3LYP), сочетающие DFT с точным обменом из метода Хартри — Фока. В 1998 году Уолтер Кон получил Нобелевскую премию по химии за развитие теории функционала плотности. К началу XXI века DFT стала доминирующим методом в квантово-химических расчётах твёрдого тела и молекулярных систем.

Математическая формулировка

Основные переменные

В DFT состояние системы описывается не волновой функцией, а электронной плотностью n(r), которая для системы из N электронов определяется как:

n(r) = N ∫ |Ψ(r, r₂, …, rₙ)|² dr₂ … drₙ

Плотность интегрируется до полного числа электронов: ∫ n(r) dr = N.

Функционал энергии

Полная энергия системы в приближении Кона — Шэма записывается как:

E[n] = Tₛ[n] + Eₑₓₜ[n] + Eₕ[n] + Eₓc[n]

где:

• Tₛ[n] — кинетическая энергия системы невзаимодействующих электронов;
• Eₑₓₜ[n] — энергия взаимодействия электронов с внешним потенциалом (ядрами);
• Eₕ[n] — энергия кулоновского взаимодействия электронов (прямая часть, энергия Хартри);
• Eₓc[n] — обменно-корреляционная энергия, содержащая все квантовые эффекты, не учтённые в предыдущих членах.

Уравнения Кона — Шэма

Для минимизации функционала E[n] решаются одноэлектронные уравнения:

[ -½ ∇² + vₑff(r) ] ψᵢ(r) = εᵢ ψᵢ(r)

где эффективный потенциал vₑff(r) = vₑₓₜ(r) + vₕ(r) + vₓc(r) включает в себя внешний потенциал, потенциал Хартри и обменно-корреляционный потенциал. Плотность восстанавливается из орбиталей Кона — Шэма: n(r) = Σᵢ |ψᵢ(r)|² (суммирование по занятым орбиталям). Уравнения решаются самосогласованно до достижения сходимости.

Обменно-корреляционные функционалы

Точность DFT определяется качеством приближения для Eₓc[n]. Различают несколько уровней приближения, образующих «лестницу Иакова» (по классификации Джона Пердью):

Приближение локальной плотности (LDA)

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

Eₓc[LDA] = ∫ εₓc(n(r)) n(r) dr

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

Обобщённое градиентное приближение (GGA)

GGA добавляет зависимость от градиента плотности ∇n(r):

Eₓc[GGA] = ∫ εₓc(n(r), |∇n(r)|) n(r) dr

Популярные функционалы GGA: PBE (Perdew-Burke-Ernzerhof), BLYP, PW91. GGA обычно точнее LDA для молекул и поверхностей.

Мета-GGA

Мета-GGA функционалы дополнительно учитывают кинетическую энергию орбиталей (или лапласиан плотности). Примеры: TPSS, SCAN.

Гибридные функционалы

Гибридные функционалы смешивают обмен из DFT с точным обменом из метода Хартри — Фока:

Eₓc[hybrid] = a Eₓ[HF] + (1-a) Eₓ[DFT] + Eₓc[DFT]

Наиболее известный гибридный функционал — B3LYP (использует три параметра), широко применяемый в органической химии.

Двойные гибриды и RPA

Более точные методы включают двойные гибриды (например, B2PLYP) и приближение случайных фаз (RPA), которые учитывают корреляцию на более высоком уровне, но требуют больших вычислительных ресурсов.

Применение

Физика твёрдого тела

DFT является основным методом расчёта зонной структуры кристаллов, плотности состояний, упругих констант, фононных спектров и магнитных свойств. С её помощью предсказывают свойства новых материалов, включая полупроводники, сверхпроводники и топологические изоляторы.

Квантовая химия

В молекулярных расчётах DFT используется для оптимизации геометрии, расчёта энергий реакций, спектров (ИК, УФ, ЯМР), потенциалов ионизации и сродства к электрону. Гибридные функционалы, такие как B3LYP, обеспечивают точность, сопоставимую с методами пост-Хартри-Фока, при значительно меньших затратах.

Нанотехнологии и материаловедение

DFT применяется для моделирования наночастиц, кластеров, графена, углеродных нанотрубок, двумерных материалов и интерфейсов. Она позволяет изучать адсорбцию молекул на поверхностях, каталитические процессы и диффузию.

Биохимия и фармакология

DFT используется для расчёта реакционных механизмов ферментов, моделирования активных центров металлоферментов и виртуального скрининга лекарственных кандидатов. Однако для больших биомолекул (белков) DFT часто комбинируют с методами молекулярной механики (QM/MM).

Ограничения и критика

Точность обменно-корреляционных функционалов

Ни один универсальный функционал не является точным для всех систем. LDA и GGA могут давать систематические ошибки: завышение стабильности некоторых структур, неправильное описание дисперсионных взаимодействий (вандерваальсовых сил) и сильнокоррелированных систем (например, оксидов переходных металлов). Для коррекции дисперсии разработаны поправки (DFT-D, DFT-D3, DFT-D4).

Самодействие

DFT в стандартных приближениях страдает от ошибки самодействия: электрон взаимодействует сам с собой через потенциал Хартри. Это приводит к неверному описанию локализованных состояний, например, в дефектах или радикалах. Разработаны методы частичного исправления (например, SIC-DFT).

Вычислительная сложность

Хотя DFT масштабируется как O(N³) (N — число базисных функций), для систем с тысячами атомов требуются суперкомпьютеры. Для больших систем (более 10⁴ атомов) применяют линеаризованные методы (линейное масштабирование, O(N)).

Альтернативные методы

Для систем с сильными корреляциями (например, высокотемпературные сверхпроводники) DFT может давать неверные результаты. В таких случаях применяют методы DFT+U, динамическую теорию среднего поля (DMFT) или квантовые методы Монте-Карло.

Программное обеспечение

DFT реализована в большинстве современных квантово-химических пакетов: VASP, Quantum ESPRESSO, Gaussian, ORCA, CP2K, ABINIT, SIESTA, FHI-aims, NWChem. Выбор пакета зависит от типа задачи (молекулы, кристаллы, поверхности) и доступных вычислительных ресурсов.

Источники

1. Hohenberg P., Kohn W. Inhomogeneous Electron Gas // Physical Review. — 1964. — Vol. 136, № 3B. — P. B864–B871.
2. Kohn W., Sham L. J. Self-Consistent Equations Including Exchange and Correlation Effects // Physical Review. — 1965. — Vol. 140, № 4A. — P. A1133–A1138.
3. Parr R. G., Yang W. Density-Functional Theory of Atoms and Molecules. — Oxford University Press, 1989.
4. Perdew J. P., Burke K., Ernzerhof M. Generalized Gradient Approximation Made Simple // Physical Review Letters. — 1996. — Vol. 77, № 18. — P. 3865–3868.
5. Becke A. D. Density-functional thermochemistry. III. The role of exact exchange // Journal of Chemical Physics. — 1993. — Vol. 98, № 7. — P. 5648–5652.
6. Martin R. M. Electronic Structure: Basic Theory and Practical Methods. — Cambridge University Press, 2004.