Метод Хартри-Фока: SCF и вариационный принцип

Метод Хартри-Фока - фундамент современной вычислительной квантовой химии. Он позволяет рассчитать электронную структуру атомов и молекул, не решая точное многочастичное уравнение Шрёдингера, а заменяя взаимодействие каждого электрона со всеми остальными усреднённым эффективным полем. Этот подход оказался настолько удачным, что даже сегодня - в эпоху метода функционала плотности и post-HF методов - Хартри-Фок остаётся стартовой точкой большинства квантово-химических расчётов. Покрутите калькулятор ниже: он показывает, как SCF-итерации сходятся к минимуму вариационной энергии для атома гелия, - а затем разберём каждый шаг.

Проблема многоэлектронного уравнения Шрёдингера

Для атома с $N$ электронами точное уравнение Шрёдингера:

$\hat{H}\Psi(\mathbf{r}_1,\dots,\mathbf{r}_N) = E\,\Psi(\mathbf{r}_1,\dots,\mathbf{r}_N)$

содержит полный гамильтониан, в котором кроме кинетики и притяжения к ядру стоят члены парного кулоновского отталкивания:

$\hat{H} = \sum_{i=1}^N \left(-\frac{1}{2}\nabla_i^2 - \frac{Z}{r_i}\right) + \sum_{i<j}\frac{1}{r_{ij}}.$

Именно двойная сумма $\sum_{i<j} 1/r_{ij}$ делает задачу неразделимой: координаты электронов перепутаны, и нельзя записать $\Psi$ как произведение одночастичных функций. Для атома гелия аналитическое решение ещё можно приблизить пертурбативно, но уже для лития точная волновая функция - трёхчастичный объект, который не вычислить в закрытой форме.

Выход, предложенный Дугласом Хартри в 1928 году и усовершенствованный Фоком в 1930-м, - заменить точное взаимодействие на задачу в среднем поле.

Вариационный принцип и волновая функция Хартри-Фока

Основа метода - вариационный принцип: для любой пробной волновой функции $\Phi$ выполняется неравенство

$E[\Phi] = \frac{\langle\Phi|\hat{H}|\Phi\rangle}{\langle\Phi|\Phi\rangle} \geq E_0,$

где $E_0$ - точная энергия основного состояния. Это значит: чем лучше $\Phi$ описывает реальную функцию, тем ниже $E[\Phi]$. Задача - минимизировать $E[\Phi]$ по параметрам пробной функции.

Рассчитать для своей задачи

В методе Хартри-Фока пробная волновая функция строится как слэтеровский детерминант одноэлектронных орбиталей $\{\phi_i(\mathbf{r})\}$:

$$\Psi_{\text{HF}} = \frac{1}{\sqrt{N!}} \begin{vmatrix} \phi_1(\mathbf{x}_1) & \phi_2(\mathbf{x}_1) & \cdots & \phi_N(\mathbf{x}_1) \\ \phi_1(\mathbf{x}_2) & \phi_2(\mathbf{x}_2) & \cdots & \phi_N(\mathbf{x}_2) \\ \vdots & & \ddots & \vdots \\ \phi_1(\mathbf{x}_N) & \cdots & & \phi_N(\mathbf{x}_N) \end{vmatrix}$$

Детерминантная форма автоматически обеспечивает антисимметрию волновой функции по перестановке электронов - принцип Паули выполняется по построению.

Уравнения Хартри-Фока и самосогласованное поле

Минимизация $E[\Psi_{\text{HF}}]$ по орбиталям $\phi_i$ приводит к системе псевдоодночастичных уравнений Хартри-Фока:

$\hat{f}\,\phi_i = \varepsilon_i\,\phi_i,$

где $\hat{f}$ - оператор Фока:

$\hat{f} = \hat{h} + \sum_{j=1}^N (\hat{J}_j - \hat{K}_j).$

Здесь $\hat{h} = -\frac{1}{2}\nabla^2 - Z/r$ - одноэлектронный гамильтониан (кинетика плюс притяжение к ядру), $\hat{J}_j$ - кулоновский оператор (среднее отталкивание от $j$-го электрона), $\hat{K}_j$ - обменный оператор, возникающий из антисимметрии слэтеровского детерминанта.

Ловушка: оператор Фока $\hat{f}$ сам зависит от орбиталей $\{\phi_i\}$, которые мы ищем. Поэтому уравнения Хартри-Фока нелинейны - их надо решать итеративно методом самосогласованного поля (SCF):

1. Задать начальное приближение для орбиталей (например, водородоподобные функции).
2. Построить оператор Фока $\hat{f}$ в текущем приближении.
3. Решить уравнения $\hat{f}\phi_i = \varepsilon_i\phi_i$, получить новые орбитали.
4. Проверить сходимость (изменение энергии или плотности меньше порога).
5. Если не сошлось - вернуться на шаг 2 с новыми орбиталями.

Рассчитать для своей задачи

Сходимость обычно достигается за 10-30 итераций для небольших молекул. Критерий: изменение полной энергии между итерациями $|\Delta E| < 10^{-6}$ а.е. (для геометрии) или $|\Delta E| < 10^{-8}$ а.е. (для свойств).

Эффективный заряд и экранирование: пример атома гелия

Для иллюстрации рассмотрим атом гелия ($Z=2$, два 1s-электрона). Пробная орбиталь: $\phi(r) = (\zeta^3/\pi)^{1/2} e^{-\zeta r}$, где $\zeta = Z_{\text{eff}}$ - эффективный заряд, вариационный параметр.

Вариационная энергия:

$E(Z_{\text{eff}}) = Z_{\text{eff}}^2 - 2Z\cdot Z_{\text{eff}} + \frac{5}{8}Z_{\text{eff}}.$

Первый член - кинетическая энергия двух 1s-электронов, второй - их притяжение к ядру с зарядом $Z=2$, третий - среднее кулоновское отталкивание между электронами (аналитическое значение $5Z_{\text{eff}}/8$).

Минимизация $dE/dZ_{\text{eff}} = 0$ даёт:

$Z_{\text{eff}}^* = Z - \frac{5}{16} = 2 - 0{,}3125 = 1{,}6875.$

Физический смысл: каждый электрон экранирует ядро от другого на $5/16 \approx 0{,}31$ единицы заряда. HF-энергия при $Z_{\text{eff}}^*$:

$E_{\text{HF}} = -\left(\frac{27}{16}\right)^2 \approx -2{,}8477 \text{ а.е.}$

Экспериментальное значение $E_0 = -2{,}9037$ а.е. Разница $\approx 0{,}056$ а.е. - это корреляционная энергия, которую метод Хартри-Фока принципиально не учитывает: он описывает движение каждого электрона в усреднённом поле, но не мгновенное уклонение электронов друг от друга.

Базисные наборы: от STO до cc-pVTZ

На практике орбитали $\phi_i$ разлагают по базисным функциям:

$\phi_i = \sum_\mu c_{\mu i}\,\chi_\mu,$

что превращает уравнения Хартри-Фока в алгебраическую задачу на собственные значения (матричные уравнения Рутаана-Холла). Качество результата определяется выбором базиса:

• STO-3G - минимальный базис, функции Слэтера аппроксимированы тремя гауссианами. Быстро, но грубо; годится для обзорных расчётов.
• 6-31G - расщеплённый базис: внутренние орбитали описываются 6 гауссианами, валентные - двумя наборами (3+1). Стандарт для органических молекул.
• 6-31G(d,p) - плюс поляризационные функции d на тяжёлых атомах и p на водороде. Учитывает деформацию орбиталей при образовании связей.
• cc-pVTZ - корреляционно-согласованный базис Данинга тройной дзета-качества. Для точных расчётов термохимии и спектроскопии.

Вычислительная стоимость метода Хартри-Фока масштабируется формально как $O(N^4)$ по числу базисных функций $N$ (четыре индекса в двухэлектронных интегралах), на практике линейно-масштабирующиеся алгоритмы снижают это до $O(N^{2..3})$ для крупных систем.

Предел Хартри-Фока и пост-HF методы

Метод Хартри-Фока даёт $\sim 99{,}5\%$ полной энергии, но оставшиеся $0{,}5\%$ - корреляционная энергия - принципиально важны для большинства химических задач: разрывы связей, барьеры реакций, межмолекулярные взаимодействия.

Для учёта электронной корреляции разработаны post-HF методы:

• MP2 (теория возмущений Мёллера-Плессе второго порядка) - самый дешёвый post-HF, $O(N^5)$. Восстанавливает $\sim 80\%$ корреляционной энергии.
• CCSD(T) - метод связанных кластеров с одиночными, двойными и пертурбативными тройными возбуждениями. «Золотой стандарт» квантовой химии, $O(N^7)$.
• CI и CASSCF - методы многоконфигурационного взаимодействия конфигураций. Необходимы при разрыве связей и для вырожденных состояний.
• DFT - метод функционала плотности: учитывает часть корреляции за счёт обменно-корреляционного функционала при формальной стоимости $O(N^3)$; основан на теореме Хоэнберга-Кона, а не на слэтеровском детерминанте.

Частые ошибки

• Путать HF-орбитали с реальными орбиталями. Орбитали Хартри-Фока - математические объекты, вспомогательные функции, а не физически наблюдаемые величины. Только плотность $\rho(\mathbf{r}) = \sum_i |\phi_i|^2$ наблюдаема.
• Забывать о корреляционной энергии. Когда метод Хартри-Фока даёт «неправильный» результат (например, диссоциация H2 в два атома He вместо H+H), причина - отсутствие корреляции, а не ошибка расчёта.
• Неверно конвертировать единицы. В атомных единицах 1 а.е. энергии (хартри) = 27,211 эВ = 627,5 ккал/моль. Типичная ошибка - смешивать кДж/моль и а.е.
• Игнорировать симметрию молекулы. Использование точечной группы симметрии ускоряет расчёт в разы и предотвращает конвергенцию к возбуждённому состоянию вместо основного.
• Интерпретировать орбитальные энергии $\varepsilon_i$ как энергии ионизации без теоремы Купманса. Теорема Купманса ($IE_i \approx -\varepsilon_i$) - приближение: она игнорирует релаксацию орбиталей после удаления электрона.

FAQ

Чем метод Хартри отличается от метода Хартри-Фока?

Метод Хартри (1928) записывает волновую функцию как простое произведение одноэлектронных орбиталей $\Psi = \phi_1\phi_2\cdots\phi_N$. Это нарушает принцип Паули: функция не антисимметрична. Фок в 1930 году заменил произведение слэтеровским детерминантом, что автоматически вводит обменное взаимодействие и гарантирует антисимметрию. Именно обменный интеграл $K_{ij} = \langle\phi_i\phi_j|r_{12}^{-1}|\phi_j\phi_i\rangle$ отличает метод Хартри-Фока от метода Хартри: он снижает энергию параллельных спинов (ферми-яма) и лежит в основе правил Хунда.

Что такое ограниченный (RHF) и неограниченный (UHF) Хартри-Фок?

В ограниченном HF (RHF) электроны с противоположными спинами занимают одни и те же пространственные орбитали - по две штуки на орбиталь. Это корректно для молекул с чётным числом электронов в замкнутой оболочке. Неограниченный HF (UHF) использует разные пространственные орбитали для $\alpha$- и $\beta$-электронов: необходим для радикалов, открытых оболочек и состояний при разрыве связи. UHF точнее описывает спиновую поляризацию, но страдает от загрязнения спиновым квантовым числом (spin contamination).

Почему метод Хартри-Фока не учитывает дисперсионные силы Ван-дер-Ваальса?

Дисперсионные взаимодействия (силы Лондона) обусловлены мгновенными флуктуациями электронной плотности и их корреляцией между молекулами. В методе HF каждый электрон движется в усреднённом поле, мгновенных флуктуаций нет по построению - корреляционная энергия равна нулю. Поэтому HF не воспроизводит связывание в комплексах, удерживаемых Ван-дер-Ваальсом (благородные газы, стэкинг ароматических колец). Для этих систем необходимы либо post-HF методы (MP2, CCSD), либо DFT с дисперсионными поправками (DFT-D3).

Коротко

Метод Хартри-Фока решает задачу многоэлектронной структуры, заменяя точное межэлектронное отталкивание на усреднённое эффективное поле: каждый электрон движется независимо в поле ядра и среднего заряда остальных. Волновая функция строится как слэтеровский детерминант - это вводит обменное взаимодействие и принцип Паули. Уравнения Хартри-Фока нелинейны и решаются итеративным методом самосогласованного поля (SCF): орбитали уточняются до самосогласования. Метод даёт $\sim 99{,}5\%$ полной энергии, но пропускает электронную корреляцию, что исправляется post-HF методами (MP2, CCSD(T)) или DFT.