Метод Хартри-Фока (HF) является фундаментальным приближением в квантовой химии для решения уравнения Шрёдингера для многолетных систем, где точное решение невозможно. Основная идея метода заключается в аппроксимации многоэлектронной волновой функции как антисимметричного произведения одноэлектронных орбиталей, что реализуется через детерминант Слэттера.
Приближение независимых частиц Электроны рассматриваются как движущиеся в эффективном поле, создаваемом ядрами и средним полем других электронов. Полная корреляция электронов игнорируется, что упрощает задачу до системы связанных одноэлектронных уравнений.
Детерминант Слэттера Многоэлектронная волновая функция Ψ строится из одноэлектронных орбиталей {ϕi} следующим образом:
$$ \Psi(\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_N) = \frac{1}{\sqrt{N!}} \begin{vmatrix} \phi_1(\mathbf{r}_1) & \phi_2(\mathbf{r}_1) & \dots & \phi_N(\mathbf{r}_1) \\ \phi_1(\mathbf{r}_2) & \phi_2(\mathbf{r}_2) & \dots & \phi_N(\mathbf{r}_2) \\ \vdots & \vdots & \ddots & \vdots \\ \phi_1(\mathbf{r}_N) & \phi_2(\mathbf{r}_N) & \dots & \phi_N(\mathbf{r}_N) \\ \end{vmatrix} $$
Такая форма обеспечивает антисимметричность волновой функции при обмене двух электронов, удовлетворяя принципу Паули.
Хартри-Фоковские уравнения Минимизация полной энергии системы по одноэлектронным орбиталям с использованием вариационного принципа приводит к системе интегродифференциальных уравнений:
F̂ϕi = εiϕi
Здесь F̂ — оператор Фока, который включает кинетическую энергию электрона, кулоновское взаимодействие с ядрами и среднее взаимодействие с другими электронами. εi — одноэлектронные энергии.
Оператор Фока F̂ определяется как сумма одноэлектронного гамильтониана ĥ и взаимодействий с другими электронами:
$$ \hat{F}(1) = \hat{h}(1) + \sum_{j=1}^{N} \left[ \hat{J}_j(1) - \hat{K}_j(1) \right] $$
где
Решение уравнений Хартри-Фока требует итеративного подхода:
Этот процесс называется самосогласованным полем (SCF, Self-Consistent Field).
Отсутствие корреляции электронов: HF учитывает только эффект обмена, игнорируя динамическую корреляцию, что приводит к завышенным значениям энергии.
Многообразие модификаций:
Энергия системы в приближении Хартри-Фока выражается через интегралы по одноэлектронным орбиталям:
$$ E_\text{HF} = \sum_i \langle \phi_i | \hat{h} | \phi_i \rangle + \frac{1}{2} \sum_{i,j} \left[ \langle \phi_i \phi_j | \hat{r}_{12}^{-1} | \phi_i \phi_j \rangle - \langle \phi_i \phi_j | \hat{r}_{12}^{-1} | \phi_j \phi_i \rangle \right] $$
Первый член отражает суммарную одноэлектронную энергию, второй — кулоновское отталкивание и обменное взаимодействие.
Метод Хартри-Фока является отправной точкой для большинства корреляционных методов квантовой химии, включая теорию возмущений (MP2, MPn), конфигурационное взаимодействие (CI) и методы плотностного функционала (DFT), где HF-орбитали часто используются как базис.
Метод HF обеспечивает баланс между вычислительной простотой и физической точностью, позволяя исследовать электронную структуру атомов, молекул и даже полимерных систем с разумными затратами ресурсов.