3. 密度泛函理论工程化理解

密度泛函理论（DFT）是量子材料电子结构计算的主力工具，但真正决定一次计算能否收敛、结果是否可信的，往往不是“会不会用软件”，而是对三个核心环节的工程化理解：Hohenberg-Kohn 定理给出了理论合法性，Kohn-Sham 方程给出了可执行的计算框架，交换关联泛函则决定了精度与成本的平衡点。本节从工程视角把这三者串成一条可操作的判断链。

一、Hohenberg-Kohn 定理：为什么电子密度可以替代波函数

Hohenberg-Kohn 第一定理指出：对于处于外势 v_ext(r) 中的相互作用电子体系，基态电子密度 n(r) 与 v_ext(r) 之间存在一一对应关系（相差一个常数除外）。这意味着基态密度唯一确定了体系的所有基态性质，包括总能量、受力、能带结构。第二定理则给出变分原理：存在一个普适泛函 F[n]，使得总能量泛函 E[n]=F[n]+∫v_ext(r)n(r)dr 在真实基态密度处取最小值。

工程上要抓住两点：其一，F[n] 是普适的，与具体材料无关，因此可以“一次构造、处处使用”；其二，定理只保证存在性，并未给出 F[n] 的显式形式。这正是 DFT 所有近似（尤其是交换关联泛函）的根源。实际计算中我们并不直接对 n(r) 做变分，而是借助 Kohn-Sham 方案把问题转化为可解的独立粒子方程。

二、Kohn-Sham 方程：把多体问题变成可执行的自洽循环

Kohn-Sham 的核心思想是引入一个无相互作用的辅助体系，使其基态密度与真实相互作用体系相同。此时电子密度由 Kohn-Sham 轨道 ψ_i 构造：n(r)=Σ_i f_i |ψ_i(r)|²，其中 f_i 为占据数。对应的单粒子方程为：

[-ħ²/2m ∇² + v_eff(r)] ψ_i(r) = ε_i ψ_i(r)
v_eff(r) = v_ext(r) + v_H(r) + v_xc(r)
v_H(r) = e² ∫ n(r')/|r-r'| dr'
v_xc(r) = δE_xc[n]/δn(r)

其中 v_H 是 Hartree 势，v_xc 是交换关联势。由于 v_eff 依赖 n(r)，而 n(r) 又依赖 ψ_i，方程必须自洽求解。工程实现上通常采用平面波基组加赝势，流程为：给定初始密度 → 构造 v_eff → 对角化 Kohn-Sham 矩阵得到 ψ_i 和 ε_i → 更新密度 → 判断能量与密度收敛。收敛判据一般取总能量变化小于 1E-5 eV 至 1E-6 eV，密度残差小于 1E-4 至 1E-5。截断能 E_cut 和 k 点网格是影响成本的两个最直接参数：E_cut 翻倍，平面波数量约增 8 倍；k 点从 4×4×4 加密到 8×8×8，计算量约增 8 倍。

三、交换关联泛函层级：精度与成本的阶梯

交换关联能 E_xc 是 DFT 中唯一需要近似的部分。按 Jacob's Ladder 的层级，从低到高依次为：

• LDA（局域密度近似）：E_xc 只依赖局域密度 n(r)。对均匀电子气精确，计算最便宜，但会高估结合能、低估晶格常数和带隙。典型误差：晶格常数偏小约 1%–3%，半导体带隙低估 30%–50%。
• GGA（广义梯度近似）：引入密度梯度 ∇n(r)，如 PBE、PW91。显著改善结合能和晶格常数，是当前材料计算的主力泛函。PBE 对大多数金属和半导体的晶格常数误差可控制在 1% 以内，但带隙仍严重低估。
• meta-GGA：进一步引入动能密度 τ(r)，如 SCAN、R2SCAN。对分子和固体的结合能、晶格常数有更好平衡，成本比 GGA 高约 1.5–2 倍。
• 杂化泛函：混合部分精确交换，如 PBE0、HSE06。HSE06 对半导体带隙的预测误差可降到 0.1–0.3 eV，但计算成本比 PBE 高 10–100 倍，且对 k 点收敛更敏感。
• 双杂化及 RPA：引入未占据轨道或响应函数，精度最高，成本也最高，通常只用于小体系或基准验证。

选择泛函时，工程判断应围绕目标量展开：若只关心结构弛豫和声子，PBE 通常足够；若关心带隙和缺陷能级，HSE06 或 meta-GGA 更合适；若关心强关联体系（如过渡金属氧化物），可能需要 DFT+U 或杂化泛函。注意 DFT+U 的 U 值不是第一性原理参数，需通过拟合实验带隙或约束 DFT 确定，典型 U 值在 3–8 eV 之间。

四、工程选择与验证：从参数到结果判断

实际计算中，泛函选择只是第一步，参数收敛同样决定结果可信度。以硅的带隙计算为例：PBE 给出约 0.6 eV（实验值 1.17 eV），HSE06 给出约 1.1–1.2 eV，接近实验。若用 PBE 做缺陷形成能计算，带隙低估会导致电荷转移能级位置错误，进而误判缺陷的稳定性。此时应至少用 HSE06 对关键缺陷做单点能校正，或采用 DFT+U 处理局域 d/f 电子。

验证策略上，建议做三件事：第一，对同一体系用两种泛函（如 PBE 和 HSE06）计算目标量，观察差异是否在可接受范围；第二，检查 E_cut 和 k 点收敛曲线，确保总能量变化小于 1 meV/atom；第三，对关键结果与实验或其他高精度方法（如 GW）交叉比对。若 PBE 与 HSE06 的带隙差异超过 0.5 eV，说明体系对交换关联近似敏感，需谨慎使用 PBE 结论。

总结来看，Hohenberg-Kohn 定理提供了理论合法性，Kohn-Sham 方程提供了可执行的自洽框架，交换关联泛函则决定了精度与成本的平衡。工程化理解的关键在于：明确目标量对泛函的敏感度，用收敛测试和交叉验证控制数值误差，并在成本允许时优先选择更高层级泛函对关键结论做校正。只有这样，DFT 计算才能从“跑得动”走向“信得过”。