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

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

История

Ранние работы Дугласа Хартри

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

Развитие Владимиром Фоком

В 1930 году советский физик Владимир Фок математически строго обобщил метод Хартри, введя в рассмотрение детерминант Слейтера. Это позволило автоматически учесть принцип Паули и антисимметричность волновой функции по отношению к перестановке двух электронов. В результате появились обменные члены, отсутствовавшие в исходной модели Хартри. Полученный подход получил название метода Хартри-Фока.

Дальнейшее совершенствование

В 1950-х годах, с развитием вычислительной техники, метод стал применяться для расчёта молекул. В 1951 году Рутан и Холл независимо предложили представить молекулярные орбитали в виде линейной комбинации атомных орбиталей (ЛКАО), что привело к созданию метода Хартри-Фока-Рутана — наиболее распространённой вычислительной реализации. В 1998 году Джон Попл и Уолтер Кон получили Нобелевскую премию по химии за развитие вычислительных методов в квантовой химии, включая метод Хартри-Фока.

Основные принципы

Волновая функция и детерминант Слейтера

В методе Хартри-Фока волновая функция системы из N электронов записывается в виде детерминанта Слейтера:

`` Ψ(1,2,...,N) = (1/√N!) * det[φ₁(1) φ₂(2) ... φₙ(N)] ``

где φᵢ — одноэлектронные спин-орбитали, а det — детерминант. Такая форма автоматически удовлетворяет принципу антисимметрии: при перестановке двух электронов знак волновой функции меняется на противоположный.

Уравнения Хартри-Фока

Уравнения Хартри-Фока представляют собой систему одночастичных уравнений:

`` F̂ φᵢ = εᵢ φᵢ ``

где F̂ — оператор Фока, εᵢ — орбитальная энергия i-й орбитали. Оператор Фока включает:

• Кинетическую энергию электрона (T̂)
• Потенциальную энергию притяжения к ядрам (V̂ₙₑ)
• Кулоновский оператор (Ĵ), описывающий среднее электростатическое отталкивание от всех остальных электронов
• Обменный оператор (K̂), возникающий из-за антисимметричности волновой функции и не имеющий классического аналога

Самосогласованное поле

Уравнения решаются итерационно: начальное приближение для орбиталей используется для вычисления оператора Фока, затем решается система уравнений, находятся новые орбитали, и процесс повторяется до достижения сходимости (когда разница между последовательными решениями становится меньше заданного порога). Этот метод называется методом самосогласованного поля (SCF).

Классификация и варианты

Ограниченный и неограниченный методы

• Ограниченный метод Хартри-Фока (RHF) — применяется для систем с чётным числом электронов, где все электроны спарены. Электроны с разными спинами занимают одну и ту же пространственную орбиталь.
• Неограниченный метод Хартри-Фока (UHF) — допускает разные пространственные орбитали для электронов с разными спинами. Используется для систем с неспаренными электронами (радикалы, возбуждённые состояния).
• Метод Хартри-Фока с ограничением по спину (ROHF) — компромиссный вариант, сохраняющий правильную спиновую симметрию при частично заполненных оболочках.

Метод Хартри-Фока-Рутана

В этом подходе молекулярные орбитали раскладываются по базисному набору атомных орбиталей (обычно гауссовых функций). Уравнения Хартри-Фока преобразуются в матричную форму:

`` FC = SCε ``

где F — матрица Фока, C — матрица коэффициентов разложения, S — матрица перекрывания, ε — диагональная матрица орбитальных энергий. Решение сводится к диагонализации матрицы Фока.

Применение

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

Метод Хартри-Фока является отправной точкой для большинства расчётов электронной структуры молекул. Он позволяет:

• Рассчитывать равновесные геометрии молекул
• Оценивать энергии связей и потенциалы ионизации
• Определять дипольные моменты и колебательные спектры
• Анализировать распределение электронной плотности

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

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

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

Атомная физика

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

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

Отсутствие электронной корреляции

Основной недостаток метода Хартри-Фока — пренебрежение мгновенными корреляциями движения электронов. Каждый электрон движется в усреднённом поле остальных, что не учитывает эффекта «избегания» электронов друг друга. Это приводит к завышению энергии системы (на 0.5–2% от полной энергии, что может составлять десятки ккал/моль для химических реакций).

Проблемы с возбуждёнными состояниями

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

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

Вычислительная сложность метода Хартри-Фока-Рутана формально составляет O(N⁴), где N — число базисных функций, что ограничивает его применение для больших систем (более нескольких сотен атомов). Для уменьшения сложности разработаны методы с использованием разреженных матриц и интегральных приближений.

Интересные факты

• Владимир Фок впервые опубликовал свою работу по методу самосогласованного поля в 1930 году в журнале Zeitschrift für Physik, и она была написана на немецком языке.
• Метод Хартри-Фока часто называют «методом среднего поля» в квантовой механике, хотя строго это неверно, так как обменный член не имеет классического аналога.
• В 1950-х годах метод был впервые реализован на электронных вычислительных машинах для расчёта атомов. Первая программа для молекул была написана в 1960-х годах.
• Метод Хартри-Фока используется как эталон для тестирования более приближённых методов, таких как теория функционала плотности (DFT).
• В 1970-х годах была показана связь метода Хартри-Фока с теорией функционала плотности через обменно-корреляционный функционал.

Источники

1. Хартри Д. Расчёты атомных структур. — М.: Издательство иностранной литературы, 1960.
2. Фок В. А. Начала квантовой механики. — Л.: Кубуч, 1932.
3. Слэтер Дж. Электронная структура молекул. — М.: Мир, 1965.
4. Попл Дж. Приближённые методы расчёта молекул. — М.: Мир, 1976.
5. Левин И. Квантовая химия. — М.: Мир, 1982.
6. Szabo A., Ostlund N. S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. — Dover Publications, 1996.