Метод Хартри — Фока

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск

Метод Хартри — Фока — в квантовой механике приближённый метод решения уравнения Шрёдингера путём сведения многочастичной задачи к одночастичной в предположении, что каждая частица двигается в некотором усреднённом самосогласованном поле, создаваемом всеми остальными частицами системы. Известно, что решение уравнения Шрёдингера позволяет получить целый ряд сведений о свойствах системы, в том числе и об её энергетическом спектре.

Метод был впервые предложен английским физиком Дугласом Хартри в 1927 году, однако содержал существенные недостатки и был впоследствии улучшен советским физиком В. А. Фоком. В отличие от Хартри, использовавшего метод самосогласованного поля с пробной волновой функцией в виде произведения одноэлектронных функций, В. А. Фок предложил в качестве пробной функции брать слэтеровский детерминант, что позволило автоматически учитывать антисимметрию полной волновой функции квантовомеханической системы по электронным переменным.[1]

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

Метод Хартри — Фока также применяется для исследования физических свойств смешанных кристаллов (например, для построения моделей распределения ионов замещения по узлам кристаллической решётки и расчета тензоров градиента электрических полей).

Введение[править | править вики-текст]

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

Построение самосогласованного поля может осуществляться либо методом последовательных приближений (изначально предложенным Хартри) или прямым вариационным методом.

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

Метод Хартри — Фока[править | править вики-текст]

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

Волновая функция многоэлектронной системы выбирается в виде детерминанта Слэтера. Уравнения Хартри — Фока являют собой одноэлектронные уравнения типа уравнения Шрёдингера, которым соответствуют орбитали \varphi_j, отвечающие минимальным значениям энергии молекулярной системы. В простейшем случае уравнения Хартри — Фока имеют вид:

\hat F[\{\varphi_j\}](1)= \hat H^{\mathrm{core}}(1)+\sum_{j=1}^{n/2}[2\hat J_j(1)-\hat K_j(1)],

где фокиан \hat F[\{\varphi_j\}](1) является оператором Гамильтона для одного электрона, находящегося в самосогласованном поле. Фокиан состоит из: суммы одноэлектронного оператора \hat H^{\mathrm{core}}(1), равного сумме оператора кинетической энергии электрона (1) и оператора потенциальной энергии его взаимодействия со всеми ядрами:

\hat H^{\mathrm{core}}(1)=-\frac{1}{2}\nabla^2_1-\sum_\alpha\frac{Z_\alpha}{r_{1\alpha}};

и суммы операторов (2\hat J_j(1)-\hat K_j(1)), определяющих взаимодействие рассматриваемого электрона (1) с усреднённым полем остальных электронов. Действие двух последних операторов на орбиталь \varphi_j определяется следующими соотношениями:

\hat J_i(1)\varphi_j(1)=\varphi_j(1)\int\frac{|\varphi_i(2)|^2}{r_{12}}\,dv_2 — оператор Кулона, учитывающий взаимодействие с орбиталью j-го электрона,
\hat K_i(1)\varphi_j(1)=\varphi_i(1)\int\frac{\varphi_i^*(2)\varphi_j(2)}{r_{12}}\,dv_2 — обменный оператор.

Основным недостатком метода является то, что он не учитывает корреляционную энергию для электронов.

Метод Хартри — Фока — Боголюбова[править | править вики-текст]

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

Метод Хартри — Фока — Дирака[править | править вики-текст]

Метод Хартри — Фока — Дирака, или Дирака — Хартри — Фока — это релятивистское обобщение метода Хартри — Фока, в основе которого лежит уравнение Дирака.

Метод Хартри — Фока — Слейтера[править | править вики-текст]

Решение уравнений Хартри — Фока значительно упрощается, если заменить обменные слагаемые (то есть члены, обязанные своим существованием антисимметрии волновой функции) некоторым усреднённым значением. Тогда они сведутся к добавлению некоторого эффективного потенциала в одноэлектронное уравнение Шрёдингера. Для вычисления этого эффективного потенциала можно воспользоваться приближением свободных электронов. Такое приближение, предложенное Джоном Слейтером,[2] а позже им же обобщённое для случая взимодействий между произвольным числом состояний, представленных слэтеровскими детерминантами,[3] носит название метода Хартри — Фока — Слейтера.

Аналогичное приближение для метода Дирака — Хартри — Фока носит название метода Дирака — Фока — Слейтера.

Метод Хартри — Фока — Рутана[править | править вики-текст]

Метод Хартри — Фока — Рутана (ХФР) — это алгебраический подход к решению уравнений Хартри — Фока, в котором неизвестные одноэлектронные функции-орбитали ищутся в виде линейных комбинаций функций заданного вида — атомных орбиталей (приближение ЛКАО).

Литература[править | править вики-текст]

  • Хартри Д. Расчёты атомных структур. — М.: ИИЛ, 1960. — 256 с.
  • Фок В. А. Начала квантовой механики. — М.: Наука, 1976. — 376 с. — Часть IV, § 3. — стр. 273—279.
  • Слэтер Дж. Методы самосогласованного поля для молекул и твердых тел / Пер. с англ. — М.: Мир, 1978. — 664 с.
  • Майер И. Избранные главы квантовой химии: доказательства теорем и вывод формул. — М.: БИНОМ. Лаборатория знаний, 2006. — 384 с. — Глава 6. Метод Хартри — Фока. — стр. 197—267.
  • Давыдов А. С. Квантовая механика. — 2-ое изд. — М.: Наука, 1973. — 704 с. — § 75. — стр. 347—353.
  • Мессиа А. Квантовая механика / Пер с франц. — М.: Наука, 1979. — Т. 2. — Глава XVIII. — стр. 254—290.

Примечания[править | править вики-текст]

  1. Давыдов А.С. Квантовая механика. — М.: Государственное издательство физико-математической литературы, 1963. — С. 391—397. — 748 с. — 35 000 экз.
  2. Slater J.C. A Simplification of the Hartree-Fock Method (англ.). — 1951. — Т. 51. — № 3. — С. 385—390. — DOI:10.1103/PhysRev.81.385
  3. Slater J.C. A Generalized Self-Consistent Field Method (англ.). — 1953. — Т. 91. — № 3. — С. 528—530. — DOI:10.1103/PhysRev.91.528

См. также[править | править вики-текст]