4. 密度泛函理论核心
密度泛函理论（DFT）之所以能成为材料计算的主力工具，核心在于它把“多电子波函数”这个高维对象，替换成了“电子密度”这个三维函数。理解这条替换链上的四个关键环节——Hohenberg-Kohn 定理、Kohn-Sham 方程、交换关联泛函、自洽场迭代——就理解了几乎所有 DFT 计算程序在做什么。
一、Hohenberg-Kohn 定理：密度足以决定一切
1964 年，Hohenberg 和 Kohn 给出两条定理，为 DFT 奠定了严格基础。
第一定理（存在性）：体系的外势 \(V_{ext}(\mathbf{r})\) 由基态电子密度 \(n(\mathbf{r})\) 唯一确定（相差一个常数）。外势确定后，哈密顿量随之确定，于是基态波函数、能量、所有可观测量的期望值都被 \(n(\mathbf{r})\) 唯一决定。换句话说，电子密度不是“一个可观测的附属量”，而是携带了体系全部基态信息的基本变量。
第二定理（变分性）：存在一个普适的能量泛函 \(E[n]\)，对任意满足粒子数守恒的试探密度 \(n(\mathbf{r})\)，都有 \(E[n] \ge E_0\)，等号仅在真实基态密度处成立。这给了我们一条可操作的路线：不去解 3N 维波函数，而是在三维密度空间里做变分极小化。
但定理只保证了泛函存在，没有给出它的显式形式。真正的困难在于电子—电子相互作用中的交换与关联部分，无法写成密度的简单解析式。
二、Kohn-Sham 方程：把难题“打包”进交换关联项
1965 年，Kohn 和 Sham 提出一个巧妙方案：构造一个无相互作用的辅助体系，使其基态密度与真实相互作用体系完全相同。这样，动能的主要部分可以用单粒子轨道精确写出，剩下的复杂多体效应全部塞进一个未知项——交换关联能 \(E_{xc}[n]\)。
总能量泛函可写为：
E[n] = T_s[n] + E_H[n] + E_ext[n] + E_xc[n]
T_s[n] : 无相互作用体系的动能
E_H[n] : Hartree 能（经典电子间库仑排斥）
E_ext[n]: 外势（离子实）贡献
E_xc[n] : 交换关联能（含所有多体修正）对密度变分后，得到一组形式上像单电子薛定谔方程的本征方程，即 Kohn-Sham 方程：
[ -ħ²/2m ∇² + V_eff(r) ] ψ_i(r) = ε_i ψ_i(r)
V_eff(r) = V_ext(r) + V_H(r) + V_xc(r)
V_xc(r) = δE_xc[n] / δn(r)
n(r) = Σ_i f_i |ψ_i(r)|²这里 \(f_i\) 是轨道占据数。注意一个关键点：KS 轨道 \(\psi_i\) 和本征值 \(\varepsilon_i\) 本身没有严格的物理意义（除最高占据轨道能量在精确泛函下等于电离能），它们只是构造密度和能量的工具。真正有物理意义的是总能量和电子密度。
三、交换关联泛函：DFT 精度的真正来源
KS 方程把多体问题化成了单粒子问题，代价是 \(E_{xc}[n]\) 形式未知。所有 DFT 计算的精度差异，几乎都来自对 \(E_{xc}\) 的近似。常用近似按“Jacob 阶梯”从低到高排列：
| 层级 | 代表泛函 | 依赖量 | 典型适用 |
|---|---|---|---|
| LDA | PZ、VWN | 局域密度 n(r) | 金属、简单固体，键长偏短 |
| GGA | PBE、PW91 | n 与 ∇n | 通用首选，表面、缺陷 |
| meta-GGA | SCAN、TPSS | n、∇n、动能密度 τ | 兼顾分子与固体 |
| 杂化 | PBE0、HSE06 | 混入精确交换 | 带隙、半导体、氧化物 |
实践中的选择原则：金属和结构优化用 PBE；半导体带隙、过渡金属氧化物用 HSE06 等杂化泛函；强关联体系（如 NiO、铜氧化物）需要 DFT+U 或更高级方法。泛函选择不是“越高级越好”，而是精度与计算成本的权衡——杂化泛函成本通常是 PBE 的几十倍。
四、自洽场迭代：从猜测密度到收敛
KS 方程是非线性的：\(V_{eff}\) 依赖 \(n\)，而 \(n\) 又来自方程的解。因此必须迭代求解，这就是自洽场（SCF）循环。
SCF 循环的每一步都有实践要点：
- 初始密度：通常由孤立原子密度叠加得到，也可用赝势提供的原子密度。
- 密度混合：直接把新密度代入下一步极易振荡甚至发散。实际程序用 Broyden 或 Pulay 混合，保留若干步历史密度信息来加速收敛。混合参数 α 一般取 0.1～0.4。
- 收敛判据：常用总能量变化小于 1×10⁻⁵ eV/atom，或密度变化小于 1×10⁻⁶。金属体系因费米面附近态密度高，需要更小的混合参数和更多迭代步。
- 收敛困难：若迭代 100 步仍不收敛，可尝试降低 α、增加展宽（smearing）、检查初始磁矩设置或结构是否合理。
五、四个环节如何串成一条链
把本节内容串起来看：Hohenberg-Kohn 定理告诉我们“用密度就够了”；Kohn-Sham 方程给出“怎么用密度算”；交换关联泛函决定“算得准不准”；自洽场迭代解决“方程怎么解”。四者缺一不可，共同构成 DFT 计算的理论骨架。
一个典型的 PBE 计算流程示例（以硅为例）：
1. 读入晶体结构（Si 金刚石结构，a = 5.43 Å）
2. 选择 PBE 泛函 + 平面波赝势
3. 设置截断能 500 eV，k 点网格 8×8×8
4. 初始密度由原子密度叠加
5. SCF 迭代：
- 构造 V_eff
- 解 KS 方程得 ψᵢ, εᵢ
- 计算 n_new
- 混合后判断收敛
6. 收敛后输出：总能量、带隙、电荷密度对硅，PBE 给出的带隙约 0.6 eV，而实验值约 1.1 eV——这正是 LDA/GGA 的“带隙低估”问题，也是杂化泛函或 GW 方法存在的理由。理解这一点，就能明白为什么泛函选择是 DFT 计算中最需要经验判断的一步。
掌握这四个环节后，面对任何 DFT 计算输入文件，你都能看懂：泛函选了什么、截断能和 k 点如何影响精度、SCF 为什么可能不收敛、结果偏差可能来自哪里。这是从“会点软件”走向“理解计算”的关键一步。