哈特里-福克方程（），又称为HF方程，是一个应用变分法计算波函数的方程式，是量子物理、凝聚态物理学、量子化学中最重要的方程之一。HF方程形式上是单电子本征方程，求得的本征态是单电子波函数，即分子轨道。以HF方程为核心的数值计算方法称为“哈特里-福克方法”（-- 基于分子轨道理论的所有量子化学计算方法都是以HF方法为基础的。鉴于分子轨道理论在现代量子化学中的广泛应用，HF方程被视为现代量子化学的基石。 哈特里-福克近似. 哈特里-福克近似也称为分子轨道近似或单行列式近似 ，认为多电子体系波函数可以由体系分子轨道波函数构造的单个斯莱特行列式表示： formula_1 这里formula_2为电子数，formula_3是能量最低的formula_2个分子轨道波函数。而分子轨道波函数是由体系中所有原子轨道波函数经过线性组合构成的， formula_5 通常原子轨道波函数formula_6由计算时指定的基组给出。 在此近似下，由变分法，体系的多电子薛定谔方程可化为单电子方程 formula_7 其中formula_8为福克算符，formula_9为体系的分子轨道，formula_10是分子轨道formula_9的轨道能。这就是哈特里-福克方程。求解该方程，即求解使方程成立的一组分子轨道；这组分子轨道使得在HF近似下的体系总能量最低。 发展简史. 1927年，物理学家瓦尔特·海特勒和弗里茨·伦敦完成了氢气分子的量子力学计算之后，开启了量子化学的时代。从此，人们便开始尝试使用量子力学理论来解释化学物质结构和化学现象。 为了解决多电子体系薛定谔方程近似求解的问题，量子化学家道格拉斯·哈特里在1928年提出了哈特里假设，将每个电子看做是在其他所有电子构成的平均势场中运动的粒子，并且首先提出了迭代法的思路。哈特里根据他的假设，将体系电子哈密顿算子分解为若干个单电子哈密顿算子的简单代数和，每个单电子哈密顿算子中只包含一个电子的坐标，因而体系多电子波函数可以表示为单电子波函数的简单乘积，这就是哈特里方程。但是由于哈特里没有考虑电子波函数的反对称要求，他的哈特里方程实际上是非常不成功的。 1930年，哈特里的学生弗拉基米尔·福克和分别提出了考虑泡利原理的自洽场迭代方程和单行列式型多电子体系波函数，这就是今天的哈特里-福克方程。但是由于计算上的困难，HF方程诞生后整整沉寂了二十年，在1950年，量子化学家想到将分子轨道用原子轨道的线性组合来近似展开，而得到了结构的罗特汉方程。 1953年，美国的、和英国的约翰·波普尔花费两年时间使用手摇计算器分别独立地实现了对氮气分子的RHF自洽场计算，这是人类首次通过求解HF方程获得对化学结构的量子力学解释，也是量子化学计算方法第一次实际完成。在第一次成功之后，伴随着电脑技术的迅猛发展，HF方程与量子化学一道获得长足发展，在HF方程的基础上，人们发展出了高级量子化学计算方法，使得计算精度进一步提高，通过对HF方程电子积分的简化和参数化，人们大大缩减了量子化学的计算量，使得对超过1000个原子的中等大小分子的计算成为可能。 方程的导出. 哈特里-福克方程源于对多电子体系电子波函数的变分法处理。在玻恩-奥本海默近似条件下，一个多电子体系的电子运动与能量可以与原子核的运动和能量相互分离，这样利用电子哈密顿算符和多电子波函数便可以计算体系的电子能量。其能量的表达式为: formula_12 式中formula_13表示体系基态电子能量；formula_14表示体系的电子哈密顿算符， formula_15 根据作用方式，可以将formula_14分解为两部分 formula_17 其中formula_18为单电子算符 formula_19 描述单个电子的动能和原子核吸引势能；而formula_20为二电子算符 formula_21 描述电子间相互作用。 formula_22代表基态多电子波函数，是体系单电子分子轨道波函数为基函数组建的斯莱特行列式。构建formula_22的各个分子轨道相互之间是正交归一的，即约束条件为 formula_24 考虑此约束条件，应用拉格朗日乘数法对函数formula_25变分求极值。式中formula_26是拉格朗日待定因子，formula_27是formula_28的简写。 变分法的处理过程如下：令 formula_29 其中 formula_30 这里定义单电子积分记号 formula_31 以及二电子积分记号 formula_32 考虑到流动坐标的不可区分性，可以简化为： formula_33 同理，formula_34中的formula_35项有： formula_36 将两项相加，formula_34表示为： formula_38 若formula_39函数处于极值点，则变量formula_40向各个方向的微小变化都应该有formula_41。可以取formula_42沿虚轴变分，则在formula_34表达式中，第一项前会产生一个formula_44的系数，第二项复共轭会产生一个formula_45系数： formula_46 消去虚数单位，两式相加，可以消去表达式中的复共轭项： formula_47 在引入库仑算符和交换算符的概念之后，上述表达式可以改写为： formula_48 对任意formula_49上述等式均应成立，因而必须有: formula_50 整理得到： formula_51 定义福克算符 formula_52 方程可以表达为 formula_53 即哈特里-福克方程。为了求解，通过对分子轨道波函数进行酉变换处理，使得轨道能矩阵formula_54对角化，将一般的不可解的HF方程转化为正则哈特里-福克方程： formula_55 此方程形式上为本征方程，但是福克算符中的库仑算符和交换算符都与分子轨道有关，因此只能够通过自洽迭代的方法近似求解，即哈特里-福克自洽场（HF-SCF）方法。HF-SCF方法是组态相互作用方法、多体微扰理论、等现代量子化学计算方法的基础。 方程的解与自洽场方法. 正则哈特里-福克方程虽然具有简单的本征方程形式，但福克算子中的库仑算子和交换算子中含有所有formula_56的表达式，因而方程的实际形式非常复杂，无法求得精确的解析解，只能使用迭代法求解，即量子化学中的自洽场方法。 在实际操作中，人们会首先给定基组，将本征方程转化为矩阵方程。这一转变使分子轨道可以表示为formula_5，可以用一个变化的formula_58维向量formula_59来代表分子轨道，formula_58为基组中基函数的数量。同理，可以将福克算符转化为福克矩阵。然而由于基函数不一定正交，存在重叠积分 formula_61 最终HF方程的形式转化为 formula_62 此方程称为哈特里-福克-罗特汉方程，或罗特汉方程。这样，相对复杂的本征方程就转化为只需要进行简单代数计算就可以求解的矩阵本征方程，而原方程中复杂的积分则在上述转化过程中一次性完成了。 求解时首先需要给定一个分子轨道的初始猜测formula_63，用这套分子轨道计算福克矩阵formula_64，再求解哈特里-福克-罗特汉方程，得到一组新的分子轨道formula_63和相应的本征值矩阵formula_66，再使用新轨道formula_63计算新的福克矩阵formula_64……重复上述过程，直到formula_63的各个矩阵元数值不再有明显的变化，称作自洽迭代收敛，就得到了HF方程的解。 得到收敛的formula_63矩阵后，将这些系数代入formula_71与基函数结合，便获得了最终的分子轨道波函数以及体系电子总能量等各种性质。 方程的应用. HF方程在量子化学中有着广泛的应用，所有分子轨道理论的量子化学计算都是以HF方程为基础的。