📚 课程目录第3章多粒子体系与近似方法

3. 多粒子体系与近似方法

真实材料中的一块晶体，哪怕只有微米尺度，也包含约 \(10^{18}\) 个原子核与电子。若直接求解这个多粒子体系的薛定谔方程，变量维度会随粒子数线性增长，计算量呈指数爆炸。因此，第一性原理计算的核心并不是"硬解"多体方程，而是通过一系列物理上合理、数学上可控的近似，把不可解的问题逐步转化为可计算的问题。本节要理清的就是这条近似链：玻恩-奥本海默近似、单电子近似、Hartree-Fock 方法，以及密度泛函理论（DFT）各自解决了什么、又留下了什么。

一、玻恩-奥本海默近似：先冻结原子核

多粒子体系的总哈密顿量包含电子动能、原子核动能、电子-电子相互作用、核-核相互作用以及电子-核相互作用。其中原子核质量至少是电子质量的 1836 倍，运动速度远慢于电子。玻恩-奥本海默近似（Born-Oppenheimer approximation）据此把电子运动和核运动分离：在求解电子状态时，把原子核视为固定在瞬时位置上；原子核的坐标只作为参数进入电子哈密顿量。

这一步的物理含义是：电子几乎瞬时地"跟随"原子核构型调整，而原子核感受到的是电子基态形成的势能面。于是总问题被拆成两层——给定核构型求电子基态能量，再在势能面上处理核的运动（振动、声子、分子动力学）。这一步没有引入电子-电子相互作用的简化，只是把核自由度暂时冻结，因此它是后续所有电子结构方法共同的前提。

二、单电子近似：把多电子问题拆成单体问题

冻结原子核后，剩下的仍是多电子问题。电子之间存在库仑排斥，每个电子的运动都与其他电子耦合，波函数同时依赖所有电子坐标。单电子近似的核心思想是：把每个电子看作在"原子核 + 其他所有电子平均场"中独立运动，用一组单电子轨道（自旋轨道）的乘积或行列式来近似多电子波函数。

这里必须区分两个层次：一是"平均场"思想，即用其他电子的平均密度产生的势替代瞬时相互作用；二是"波函数形式"，即如何构造满足泡利不相容原理的多电子波函数。前者是物理近似，后者是数学表示。Hartree 最早提出简单乘积形式，但没有考虑电子交换反对称性；Fock 和 Slater 用 Slater 行列式替代乘积，自动满足反对称要求，这才形成 Hartree-Fock 方法。

三、Hartree-Fock 方法：交换能精确，关联能缺失

Hartree-Fock（HF）方程可以写成单电子形式：

F̂ ψ_i(r) = ε_i ψ_i(r)

F̂ = ĥ + Σ_j (Ĵ_j − K̂_j)

ĥ：单电子动能 + 核吸引势
Ĵ_j：第 j 个电子产生的库仑势（经典排斥）
K̂_j：交换算符（来自反对称波函数，无经典对应）
ε_i：第 i 个单电子轨道能量

HF 方法通过自洽场（SCF）迭代求解：先猜一组初始轨道，构造平均场，解方程得到新轨道，再重新构造平均场，直到轨道和能量收敛。它的优点是交换能（exchange energy）被精确处理，因为交换项直接来自 Slater 行列式的反对称性。但 HF 只包含平均场，电子之间的瞬时关联（动态关联）被完全忽略，导致总能量偏高，键能、带隙等性质误差明显。这个缺失的部分称为关联能 \(E_c = E_{exact} - E_{HF}\)。

四、密度泛函理论：从波函数转向电子密度

HF 仍然以多电子波函数为基本变量，计算量随电子数增长很快。密度泛函理论（DFT）换了一个基本变量：电子密度 \(n(r)\)。Hohenberg-Kohn 定理指出，基态电子密度唯一决定外势，从而决定体系所有基态性质；Kohn-Sham 定理则构造了一个与真实体系具有相同密度的非相互作用参考体系，把问题转化为一组单电子方程：

[ −½∇² + v_ext(r) + v_H(r) + v_xc(r) ] φ_i(r) = ε_i φ_i(r)

v_ext：原子核外势
v_H：Hartree 势（电子密度产生的经典库仑势）
v_xc：交换-关联势，包含所有多体效应
n(r) = Σ_i |φ_i(r)|²  （占据轨道求和）

DFT 与 HF 的关系可以这样理解：两者都采用单电子方程和 SCF 迭代，但 HF 的交换项是精确的非局域算符，而 DFT 把所有交换和关联效应都塞进一个待近似的 \(v_{xc}\) 中。若 \(v_{xc}\) 精确，DFT 原则上给出精确基态密度和能量；实际计算中 \(v_{xc}\) 必须近似，常用 LDA、GGA（如 PBE）、meta-GGA、杂化泛函等。杂化泛函正是把 HF 的精确交换按一定比例混入 DFT，以改善带隙和局域化描述。

五、四者关系与选择依据

把这条链串起来看：玻恩-奥本海默近似分离核与电子，单电子近似把多电子问题化为单体问题，HF 和 DFT 是单电子近似下的两条不同技术路线。HF 的交换是精确的，但缺少关联；DFT 把交换和关联都放进 \(v_{xc}\)，理论上更完整，但 \(v_{xc}\) 必须近似。实际材料计算中，DFT 因计算效率高、对大多数基态性质描述合理而成为主流；HF 及其后处理（MP2、CCSD(T)）更多用于小分子高精度基准；杂化泛函则试图在两者之间取得平衡。

理解要点：近似不是"错误"，而是有明确适用边界的模型。判断一个计算结果是否可信，首先要问：用了哪个近似？交换-关联泛函是什么？体系是否强关联？带隙、磁性、缺陷能级等性质对泛函选择是否敏感？

常见误区：把 DFT 的 Kohn-Sham 轨道能量直接当作真实准粒子能级。严格来说，Kohn-Sham 带隙与实验带隙之间存在"导数不连续"问题，需要 GW 等准粒子方法修正。HF 轨道能量也不等于电离能（Koopmans 定理只是近似）。

从工程实践角度，选择方法时可以按以下顺序判断：先确认体系是否含重元素（需考虑相对论效应）、是否强关联（如过渡金属氧化物、f 电子体系）；再确认目标性质是基态能量、结构、带隙还是激发态；最后根据精度和成本权衡选择泛函或后 HF 方法。这条近似链不是历史陈列，而是每一次计算设置背后的物理依据。