3. 密度泛函理论(DFT)核心：Hohenberg-Kohn定理、Kohn-Sham方程、交换关联泛函(LDA/GGA)
好，咱们今天来啃DFT这块硬骨头。说实话，我刚入行那会儿，觉得DFT就是个黑盒子——扔进去一个结构，出来一堆能量和能带。但做材料筛选，不懂核心原理，你连算出来的结果对不对都不敢确认。
密度泛函理论，说白了就是解决一个核心问题：多电子体系的薛定谔方程太难解了。电子之间互相排斥，你推我我推你，直接求解的计算量随电子数指数增长。DFT的高明之处在于——它把问题简化了。
核心思想：用电子密度 ρ(r) 代替多电子波函数 Ψ(r₁, r₂, ..., rₙ)，把3N维问题降到3维。
3.1 Hohenberg-Kohn定理：DFT的基石
1964年，Hohenberg和Kohn发了篇论文，奠定了整个DFT大厦的地基。两个定理，简洁有力。
第一定理：电子密度决定一切
系统的基态电子密度 ρ₀(r) 唯一决定了外势场 v(r)，进而决定了系统的所有性质。换句话说：知道密度，就知道一切。
这个定理的证明其实很巧妙——反证法。假设两个不同的外势场能产生相同的基态密度，最后会推出矛盾。嗯，这里不展开数学细节，你只要记住结论就行。
我在做高熵合金筛选时，经常需要快速判断不同构型的稳定性。有了这个定理，我只需要关注电子密度的变化，而不必去纠结每个电子的具体位置。这在实际计算中省了不少事。
第二定理：变分原理
对于任何尝试密度 ρ̃(r)，对应的能量 E[ρ̃] ≥ E₀（基态能量）。等号成立当且仅当 ρ̃ = ρ₀。
这给了我们一个实用的方法：不断调整电子密度，让能量降到最低，找到的就是基态。
我的经验：在实际计算中，我们很少直接对密度做变分。Kohn-Sham方程的引入，让这件事变得可操作得多。但理解这个变分思想很重要——它告诉你，DFT本质上是一个优化问题。
3.2 Kohn-Sham方程：把多体问题变成单体问题
Hohenberg-Kohn定理告诉我们理论可行，但没告诉具体怎么算。1965年，Kohn和Sham给出了一个漂亮的方案。
他们的想法是这样的：引入一个假想的无相互作用体系，这个体系中的电子在一个有效势场中运动，产生的电子密度与真实相互作用体系完全相同。
这样一来，复杂的多体问题就变成了一个简单的单体问题。Kohn-Sham方程长这样：
[-½∇² + v_eff(r)] ψ_i(r) = ε_i ψ_i(r)
其中 v_eff(r) = v_ext(r) + v_Hartree(r) + v_xc(r)
- v_ext(r): 外势（原子核势场）
- v_Hartree(r): 哈特里势（电子间的经典库仑相互作用）
- v_xc(r): 交换关联势（包含所有非经典相互作用）这个方程需要自洽求解：
- 猜一个初始电子密度 ρ_in
- 构造有效势 v_eff
- 解Kohn-Sham方程，得到轨道 ψ_i
- 由轨道构造新的电子密度 ρ_out
- 混合 ρ_in 和 ρ_out，得到新的输入密度
- 重复直到收敛
这个过程，我们叫它自洽场（SCF）迭代。我刚开始做计算时，经常遇到SCF不收敛的情况——要么是初始猜测太离谱，要么是混合参数没调好。后来我养成了一个习惯：先用粗网格算几步，再切到细网格，收敛会稳定很多。
注意：Kohn-Sham方程中的轨道 ψ_i 并不是真实电子的波函数，而是准粒子的波函数。它们的物理意义是辅助量，不要过度解读。但它们的能量 ε_i 在Koopmans定理的意义下，可以近似对应电离能和电子亲和能。
3.3 交换关联泛函：LDA与GGA
Kohn-Sham方程把问题都推到了 v_xc 上。这个交换关联泛函包含了所有我们不知道的东西——交换能、关联能、自相互作用修正等等。
问题是：交换关联泛函的精确形式是未知的。我们需要做近似。这就好比做菜，你知道最终要什么味道，但调料配方得自己摸索。
LDA：局域密度近似
LDA假设：在r处的交换关联能密度，只取决于该点的电子密度 ρ(r)。
E_xc^LDA[ρ] = ∫ ρ(r) ε_xc(ρ(r)) dr其中 ε_xc(ρ) 是从均匀电子气模型精确计算出来的。对于均匀电子气，这个泛函是精确的。
LDA的优点：简单、稳定、对很多体系效果不错。缺点：高估结合能、低估晶格常数。我在算金属的弹性常数时，LDA给出的结果往往偏硬。
GGA：广义梯度近似
GGA在LDA的基础上，引入了密度梯度 ∇ρ(r) 的信息：
E_xc^GGA[ρ] = ∫ ρ(r) ε_xc(ρ(r), |∇ρ(r)|) dr常用的GGA泛函有PBE、PW91等。PBE是我用得最多的，它在固体和表面体系上表现均衡。
| 泛函 | 特点 | 适用场景 | 我的评价 |
|---|---|---|---|
| LDA | 简单、计算快 | 金属、简单固体 | 入门级，但别太依赖 |
| PBE | 均衡、通用 | 固体、表面、分子 | 我的首选，90%场景够用 |
| RPBE | 改进吸附能 | 表面催化 | 算吸附时用这个 |
| PBEsol | 改进晶格常数 | 固体结构优化 | 算晶格常数时推荐 |
我曾经在筛选锂离子电池正极材料时，用PBE算出来的电压总是偏低0.3V左右。后来换成PBEsol，结果跟实验对上了。所以，选泛函不是一劳永逸的事，要根据你的具体问题来定。
避坑指南：我曾经在算二维材料时，直接用PBE算带隙，结果跟实验差了1eV以上。后来才知道，PBE系统性地低估带隙。对于这类问题，要么用杂化泛函（如HSE06），要么用GW方法。记住：没有万能泛函，每个近似都有它的适用范围。
3.4 知识体系总览
下面这张图，是我自己总结的DFT核心逻辑。每次做计算前，我都会在脑子里过一遍这个流程：
这张图把DFT的脉络理得很清楚：从问题出发，经过HK定理的理论保证，到KS方程的具体实现，最后落到交换关联泛函的近似选择。每一步都有它的物理意义和实际考量。
做材料筛选时，我一般这样操作：先用PBE快速扫描一批候选材料，筛选出有潜力的几个，再用更精确的方法（如杂化泛函或考虑自旋极化）做精细计算。这样既保证了效率，又不牺牲精度。
实用建议：如果你刚开始做DFT计算，我建议你先用PBE跑通整个流程。等你对结果有了感觉，再尝试不同的泛函。记住：计算只是工具，物理直觉才是核心。多跟实验数据对比，慢慢你就能判断哪些结果是可靠的，哪些需要进一步验证。