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

Теория функционала плотности (Шаблон:Lang-en) — метод расчёта электронной структуры систем многих частиц в квантовой физике и квантовой химии. В частности, применяется для расчёта электронной структуры молекул и конденсированного вещества. Является одним из наиболее широко используемых и универсальных методов в вычислительной физике и вычислительной химии. Твёрдое тело рассматривается как система, состоящая из большого числа одинаково взаимодействующих между собой электронов, удерживаемых вместе решёткой из атомных ядер. Основная идея метода заключается в использовании понятия электронной плотности в основном состоянии, её распределение описывается одночастичным уравнением Шрёдингера.

Введение

Модель Томаса — Ферми

Шаблон:Main Методу теории функционала плотности предшествовала модель Томаса — Ферми, развитая Л. Томасом и Энрико Ферми в 1927 г. Они рассчитали энергию атома как сумму его кинетической энергии, представленной в виде функционала электронной плотности, и потенциальной энергии взаимодействия электронов с ядром и друг с другом; энергия взаимодействия также была выражена через электронную плотность.

Несмотря на заметную роль, которую модель Томаса — Ферми сыграла в развитии квантовой механики, её точность была недостаточной, поскольку не учитывалось обменное взаимодействие, в отличие, например, от метода Хартри — Фока. В 1928 г. Поль Дирак уточнил функционал энергии в модели Томаса — Ферми, добавив к нему слагаемое, описывающее обменное взаимодействие (это слагаемое также имело вид функционала электронной плотности).

Несмотря на это, для ряда применений модель Томаса — Ферми — Дирака не давала удовлетворительного результата. Основным источником погрешности являлось выражение кинетической энергии, приводящее к погрешности в вычислении обменной энергии. Кроме того, не учитывалась энергия электронной корреляции.

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

Хотя теория функционала плотности и базируется на ставшей классической модели Томаса — Ферми, надёжное теоретическое обоснование под неё было подведено только с формулировкой теорем Хоэнберга — Кона (названных так в честь Пьера Хоэнберга и Уолтера Кона).

В первой теореме доказано, что свойства основного состояния многоэлектронной системы определяются только электронной плотностью, зависящей от трех координат. Данная теорема сводит задачу об описании много-электронной системы из N электронов с 3N пространственными координатами к описанию функционала электронной плотности с тремя координатами.

Вторая теорема представляет собой вариационный принцип квантовой механики, сформулированный для функционала плотности, и утверждает, что энергия электронной подсистемы, записанная как функционал электронной плотности, имеет минимум, равный энергии основного состояния.

Первоначально теоремы Хоэнберга — Кона были сформулированы только для основного состояния электронной подсистемы в отсутствие магнитного поля. Они могут быть обобщены путём введения зависимости от времени, что позволяет использовать этот формализм для расчета состояний возбуждённых электронов.

Описание метода

Традиционные методы определения электронной структуры, в частности, метод Хартри — Фока и производные от него, описывают систему с помощью многоэлектронной волновой функции. Основная цель теории функционала плотности — при описании электронной подсистемы заменить многоэлектронную волновую функцию электронной плотностью. Это ведет к существенному упрощению задачи, поскольку многоэлектронная волновая функция зависит от $3N$ переменных — по 3 пространственных координаты на каждый из $N$ электронов, в то время как плотность — функция лишь трёх пространственных координат.

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

Описание двух последних взаимодействий и представляет собой основную сложность метода теории функционала плотности в формулировке Кона — Шэма. Простейшим приближением здесь является приближение локальной плотности, основанное на точном расчёте обменной энергии для пространственно однородного электронного газа, который может быть выполнен в рамках модели Томаса — Ферми, и из которого можно получить также и корреляционную энергию электронного газа.

Метод теории функционала плотности широко применяется для расчётов в физике твёрдого тела с 1970-х годов. В ряде случаев даже использование простого приближения локальной плотности дает удовлетворительные результаты, соответствующие экспериментальным данным, причём вычислительная сложность метода невысока относительно других подходов к проблеме многих частиц в квантовой механике. Тем не менее, долгое время метод был недостаточно точен для расчётов в области квантовой химии, пока в 1990-х годах не произошёл заметный сдвиг в описании обменного и корреляционного взаимодействий. В настоящее время метод теории функционала плотности является главным подходом в обеих областях. Впрочем, несмотря на прогресс в теории, все ещё имеются проблемы в приложении метода к описанию межмолекулярных сил, в особенности Ван-дер-Ваальсовых сил и дисперсионного взаимодействия, а также в расчётах ширины запрещённой зоны в полупроводниках.

Сложности с расчётом дисперсионного взаимодействия в рамках теории функционала плотности (которые возникают, как минимум, в том случае, когда этот метод не дополняется другими) делают метод теории функционала плотности малопригодным для систем, в которых дисперсионные силы являются преобладающими (например, при рассмотрении взаимодействия между атомами благородных газов) или систем, в которых дисперсионные силы имеют тот же порядок, что и другие взаимодействия (например, в органических молекулах). Решение этой проблемы является предметом современных исследований.

Формальное обоснование метода

Согласно приближению Борна — Оппенгеймера, которое применяется в большинстве расчётов электронной структуры, ядра, входящие в состав рассматриваемой системы, считаются неподвижными. Электростатический потенциал $V$, создаваемый этими «неподвижными» ядрами, является внешним для электронов. Стационарное состояние электронов описывается волновой функцией $\mathrm{\Psi }({\overrightarrow{r}}_1,\dots ,{\overrightarrow{r}}_N)$, которая является решением уравнения Шрёдингера

$H\mathrm{\Psi }=[T+V+U]\mathrm{\Psi }=[\sum _i^N-\frac{{\hslash }^2}{2m}\underset{i}{\overset{2}{\mathrm{\nabla }}}+\sum _i^{N}V({\overrightarrow{r}}_i)+\sum _{i<j}U({\overrightarrow{r}}_i,{\overrightarrow{r}}_j)]\mathrm{\Psi }=E\mathrm{\Psi },$

где $H$ — гамильтониан электронной подсистемы, $N$ — количество электронов, $U$ описывает электрон-электронное взаимодействие. Операторы $T$ и $U$ одинаковы для всех систем, в то время как вид $V$ зависит от конкретной системы. Как видно, основное отличие одночастичной задачи от задачи многих тел состоит в наличии слагаемого, описывающего электрон-электронное взаимодействие, $U$. Существует большое количество методов решения многочастичного уравнения Шрёдингера, основанных на разложении волновой функции с использованием определителя Слэтера. Простейший из них — метод Хартри — Фока, на основе которого развит ряд современных методов. Общей проблемой для них является значительная вычислительная трудоёмкость, из-за которой область применения метода Хартри — Фока и производных от него ограничена не слишком большими системами.

Метод теории функционала плотности в значительной степени решает проблему расчёта систем, включающих большое число частиц, путём сведения задачи о системе многих тел с потенциалом электрон-электронного взаимодействия $U$ к одночастичной задаче, в которой слагаемое $U$ отсутствует.

Плотность частиц, $n(\overrightarrow{r})$, с помощью которой и строится формализм теории функционала плотности, задается выражением:

$n(\overrightarrow{r})=N\int d^3{r}_2\int d^3{r}_3\dots \int d^3{r}_N{\mathrm{\Psi }}^{*}(\overrightarrow{r},{\overrightarrow{r}}_2,\dots ,{\overrightarrow{r}}_N)\mathrm{\Psi }(\overrightarrow{r},{\overrightarrow{r}}_2,\dots ,{\overrightarrow{r}}_N).$

Хоэнберг и Кон в 1964 показали, что это выражение может быть обращено: по заданной плотности частиц в основном состоянии, $n_0(\overrightarrow{r})$, можно найти соответствующую волновую функцию основного состояния ${\mathrm{\Psi }}_0({\overrightarrow{r}}_1,\dots ,{\overrightarrow{r}}_N)$. Иными словами, ${\mathrm{\Psi }}_0$ — единственный функционал от $n_0$, то есть

${\mathrm{\Psi }}_0={\mathrm{\Psi }}_0[n_0],$

а, следовательно, все остальные наблюдаемые физические величины $O$ также являются функционалами $n_0$:

$⟨O⟩[n_0]=⟨{\mathrm{\Psi }}_0[n_0]\mid O\mid {\mathrm{\Psi }}_0[n_0]⟩.$

В частности, для энергии основного состояния можно записать

$E_0=E[n_0]=⟨{\mathrm{\Psi }}_0[n_0]\mid T+V+U\mid {\mathrm{\Psi }}_0[n_0]⟩,$

где вклад внешнего потенциала $⟨{\mathrm{\Psi }}_0[n_0]\mid V\mid {\mathrm{\Psi }}_0[n_0]⟩$ может быть переписан через плотность частиц:

$V[n]=\int V(\overrightarrow{r})n(\overrightarrow{r})d^{3}r.$

Функционалы $T[n]$ и $U[n]$ одинаковы для всех систем, а $V[n]$, очевидно, зависит от вида рассматриваемой системы. Для заданной системы вид $V$ известен, и можно минимизировать функционал

$E[n]=T[n]+U[n]+\int V(\overrightarrow{r})n(\overrightarrow{r})d^{3}r$

относительно распределения плотности частиц $n(\overrightarrow{r})$, если, конечно, имеются выражения для $T[n]$ и $U[n]$. В результате минимизации получается плотность частиц в основном состоянии $n_0$, а вместе с ней и все наблюдаемые в основном состоянии величины.

Вариационная задача отыскания минимума функционала энергии $E[n]$ может быть решена с помощью метода множителей Лагранжа, как это и было сделано Коном и Шэмом в 1965 г. Таким образом, функционал энергии в приведённом выше выражении может быть записан как эффективный функционал плотности частиц в одночастичной системе:

$E_s[n]=⟨{\mathrm{\Psi }}_s[n]\mid T_s+V_s\mid {\mathrm{\Psi }}_s[n]⟩,$

где $T_s$ означает кинетическую энергию свободной частицы, а $V_s$ — эффективный внешний потенциал для электронной подсистемы. Ясно, что $n_s(\overrightarrow{r})\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}n(\overrightarrow{r})$ если $V_s$ взят в виде

$V_s=V+U+(T-T_s).$

Решение так называемых уравнений Кона — Шэма для вспомогательной системы, из которой исключено электрон-электронное взаимодействие,

$[-\frac{{\hslash }^2}{2m}{\mathrm{\nabla }}^2+V_s(\overrightarrow{r})]{\phi }_i(\overrightarrow{r})={\epsilon }_i{\phi }_i(\overrightarrow{r}),$

даёт орбитали ${\phi }_i$, по которым восстанавливается электронная плотность $n(\overrightarrow{r})$ исходной многочастичной системы:

$n(\overrightarrow{r})\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}{n}_s(\overrightarrow{r})=\sum _i^N|{\phi }_i(\overrightarrow{r}){|}^2.$

Эффективный одночастичный потенциал $V_s$ записывается как

$V_s=V+\int \frac{e^2{n}_s(\overrightarrow{r}{}^{\prime })}{|\overrightarrow{r}-\overrightarrow{r}{}^{\prime }|}{d}^3{r}^{\prime }+V_{\mathrm{X}\mathrm{C}}[n_s(\overrightarrow{r})],$

где второе слагаемое — слагаемое Хартри — описывает электрон-электронное кулоновское отталкивание, а последнее слагаемое $V_{\mathrm{X}\mathrm{C}}$ называется обменно-корреляционным потенциалом. Здесь $V_{\mathrm{X}\mathrm{C}}$ включает все многочастичные взаимодействия.

Поскольку слагаемое Хартри и член $V_{\mathrm{X}\mathrm{C}}$ зависят от плотности $n(\overrightarrow{r})$, которая зависит от ${\phi }_i$, которая, в свою очередь, зависит от $V_s$, решение самосогласованных уравнений Кона — Шэма может быть произведено с помощью итеративной процедуры последовательных приближений. Как правило, отталкиваясь от начального приближения для $n(\overrightarrow{r})$, рассчитывается соответствующее слагаемое $V_s$, для которого затем решаются уравнения Кона — Шэма, из которых получается ${\phi }_i$. Отсюда можно получить следующее приближение для плотности и т. д.

Приближения

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

$E_{\mathrm{X}\mathrm{C}}[n]=\int {\epsilon }_{\mathrm{X}\mathrm{C}}(n)n(r)d^{3}r.$

Приближение локальной спиновой плотности (LSDA) является непосредственным обобщением приближения локальной плотности, учитывающим спин электрона:

$E_{\mathrm{X}\mathrm{C}}[n_{\uparrow },n_{\downarrow }]=\int {\epsilon }_{\mathrm{X}\mathrm{C}}(n_{\uparrow },n_{\downarrow })n(r)d^{3}r.$

Достаточно точное выражение для плотности обменно-корреляционной энергии ${\epsilon }_{\mathrm{X}\mathrm{C}}(n_{\uparrow },n_{\downarrow })$ было получено с помощью квантового метода Монте-Карло при расчётах газа свободных электронов.

Метод обобщённого градиентного приближения (GGA) также является локальным, но, в отличие от метода локальной плотности, учитывает градиент плотности в точке рассмотрения:

$E_{\mathrm{X}\mathrm{C}}[n_{\uparrow },n_{\downarrow }]=\int {\epsilon }_{\mathrm{X}\mathrm{C}}(n_{\uparrow },n_{\downarrow },\overrightarrow{\mathrm{\nabla }}{n}_{\uparrow },\overrightarrow{\mathrm{\nabla }}{n}_{\downarrow })n(r)d^{3}r.$

Использование этого приближения дает хорошие результаты при расчете геометрии и энергии основного состояния молекул.

Существуют и более точные приближения, которые в значительной степени позволяют решить проблему вычисления функционала обменно-корреляционной энергии.

Обобщение на случай магнитного поля

Формализм метода теории функционала плотности нарушается в условиях наличия векторного потенциала, в частности, в присутствии магнитного поля. В этом случае не существует взаимно однозначного соответствия между электронной плотностью и внешним потенциалом (атомных ядер). Попытки обобщения формализма для учёта эффектов, связанных с магнитным полем, вылились в две разных теории: в теорию функционала плотности с учётом вектора плотности тока, и в теорию функционала плотности с учётом магнитного поля. В обоих случаях функционал обменно-корреляционной энергии обобщается и становится зависящим не только от электронной плотности. В первом подходе, развитом Vignale и Rasolt, помимо электронной плотности, аргументом является ещё и плотность тока. Во втором подходе (Salsbury, Grayce, Harris) дополнительным аргументом функционала служит магнитное поле, и вид функционала зависит от вида магнитного поля. Для обоих методов вычисление обменно-корреляционной энергии за рамками приближения локальной плотности (вернее, его обобщения на случай магнитного поля) оказалось весьма сложным.

Применения

На практике, метод Кона — Шэма может быть применён несколькими различными способами, в зависимости от цели исследования. В расчётах для физики твёрдого тела до сих пор широко используется приближение локальной плотности, вкупе с базисом плоских волн. Для расчётов электронной структуры молекул требуются более сложные выражения для функционалов. Так, большое число приближенных функционалов для расчёта обменно-корреляционного взаимодействия было развито для задач химии. Некоторые из них противоречат приближению пространственно однородного электронного газа, но, тем не менее, в пределе при переходе к электронному газу должны сводиться к приближению локальной плотности.

Для расчётов физических задач наиболее часто применяется, по-видимому, уточнённая обменная модель Perdew — Burke — Ernzerhof, однако известно, что она приводит к ошибкам в калориметрических параметрах, будучи приложенной к расчётам молекул в газовой фазе.

В расчётах квантовой химии одним из распространённых является вид обменно-корреляционного функционала, называемый BLYP (Becke, Lee, Yang, Parr). Еще более широко распространено приближение B3LYP, которое основано на гибридном функционале, в котором обменная энергия рассчитывается с привлечением точного результата, полученного методом Хартри — Фока.

В целом, текущее состояние метода теории функционала плотности таково, что невозможно оценить погрешность расчёта, не сравнивая его результаты с другими подходами или с результатами экспериментов.

Программное обеспечение, которое реализует метод теории функционала плотности

Abinit — использует метод псевдопотенциала ADF AIMPRO Atomistix Toolkit CADPAC CASTEP Шаблон:Iw CPMD CRYSTAL06 DACAPO DALTON deMon2K DFT++ DMol3 EXCITING: полный потенциал, метод линеаризованных присоединённых плоских волн (LAPW) Fireball GAMESS (UK) GAMESS (US) GAUSSIAN HyperChem JAGUAR TB-LMTO-47 — метод линеаризованных маффин-тин орбиталей

MOLCAS MOLPRO MPQC NRLMOL NWChem OCTOPUS OpenMX ORCA ParaGauss [1] Шаблон:Wayback PC GAMESS начиная с версии 6.4 PLATO — метод линейной комбинации атомных орбиталей ПРИРОДА PWscf (Quantum-ESPRESSO) — использует метод псевдопотенциала Q-Chem SIESTA[2] — линейно-масштабируемый метод теории функционала плотности, использует метод псевдопотенциала Spartan S/PHI/nX TURBOMOLE VASP — использует метод псевдопотенциала WIEN2k — полный потенциал, метод линеаризованных присоединённых плоских волн (LAPW)

Литература

• Марч Н., Кон В., Вашишта П. и др. Теория неоднородного электронного газа. — Шаблон:М: Мир, 1987.
• Dreizler R., Gross E. Density Functional Theory. — Plenum Press, New York, 1995.
• Koch W., Holthausen M. C. A Chemist’s Guide to Density Functional Theory. — ed. 2. — Weinheim: Wiley-VCH, 2002.
• Parr R. G., Yang W. Density-Functional Theory of Atoms and Molecules. — New York: Oxford University Press, 1989.
• Шаблон:Статья
• Шаблон:Статья

Ссылки

• Walter Kohn (Нобелевский лауреат) Видеозапись интервью с Вальтером о его работе по разработке теории функционала плотности, Vega Science Trust. Шаблон:Ref-en
• Klaus Capelle A bird’s-eye view of density-functional theory Шаблон:Ref-en
• Walter Kohn Nobel Lecture Шаблон:Ref-en
• Кон В. Нобелевская лекция.