4. 交换关联泛函与赝势工程选型

在DFT能带与态密度计算中，交换关联泛函和赝势是决定结果可信度的两个最底层选择。它们不像截断能或k点网格那样可以通过收敛测试自动逼近“正确答案”，而是直接改变所求解的物理问题本身。选错泛函，带隙可能偏差50%以上；选错赝势，同一泛函下态密度峰位也会漂移。本节从工程选型角度，梳理四类泛函的适用边界、两类赝势的取舍逻辑，以及如何评估它们对能带和态密度的影响。

一、四类泛函的适用场景与典型误差

LDA（局域密度近似）只依赖局域电子密度，计算量最小，对高密度均匀电子气体系表现尚可。但在半导体和绝缘体中，LDA系统性低估带隙，典型如Si的实验带隙1.17 eV，LDA仅给出约0.5 eV；对过渡金属氧化物如NiO，LDA甚至给出金属性基态，与实验的绝缘体行为定性矛盾。LDA适合作为快速预筛或与实验趋势做定性对比，不适合直接用于带隙工程。

GGA（广义梯度近似）引入密度梯度，其中PBE是最广泛使用的形式。PBE对晶格常数、结合能和弹性常数的描述通常优于LDA，但带隙仍系统性偏小，Si约0.6～0.7 eV。GGA适合结构优化、声子计算和大多数金属体系的态密度分析。对于需要可靠带隙的半导体，GGA只能作为起点，不能作为终点。

meta-GGA在密度和梯度之外引入动能密度，SCAN是代表性泛函。SCAN对半导体带隙的改进有限，但对强关联体系和表面能的描述更均衡。其工程价值在于：当GGA对某类体系的能量排序明显错误时，SCAN往往能给出更合理的趋势。代价是计算量约为PBE的2～3倍，且对数值设置更敏感。

杂化泛函通过混入精确交换提升带隙。HSE06是最常用的屏蔽杂化泛函，对Si带隙可给出约1.1～1.2 eV，接近实验值。但HSE06的计算量约为PBE的10～50倍，且对k点网格收敛更慢。工程上通常采用“PBE优化结构 + HSE06单点能带”的策略，而非全程使用杂化泛函。

二、PAW与超软赝势的工程对比

PAW（投影缀加波）和超软赝势（USPP）都通过引入补偿电荷来恢复芯区附近的电子密度，但实现路径不同。PAW在芯区使用全电子波函数与赝波函数的线性变换，理论上可以精确重建全电子密度；USPP则通过经验参数构造平滑赝波函数，芯区信息以拟合方式保留。

工程选型时需关注三个维度：精度——PAW对磁矩、超精细场和芯能级位移的描述通常更可靠，尤其对含d/f电子的体系；计算成本——同截断能下PAW比USPP慢约20%～40%，但PAW对截断能的收敛更快，实际总耗时差距可能缩小；可移植性——PAW数据集通常有明确的生成参数和验证流程，不同版本间一致性较好，而USPP的拟合质量更依赖生成者经验。

以VASP为例，PAW_PBE数据集是默认推荐；若使用USPP，需额外检查截断能是否足够（通常比PAW高20%～30%）。对于含氧体系，PAW对O 2p与金属d轨道的杂化描述更稳定，态密度中费米面附近的峰形更少出现伪影。

三、对能带和态密度的影响评估

泛函和赝势的选择会同时影响能带色散和态密度分布。评估时建议固定一个变量、观察另一个变量的响应。例如：固定PAW_PBE，切换LDA、PBE、SCAN、HSE06，记录带隙、导带底位置和价带顶有效质量的变化；再固定HSE06，切换PAW与USPP，观察态密度中过渡金属d峰与氧p峰的相对位置是否移动。

一个可操作的评估流程是：先做PBE结构优化，得到平衡晶格常数；再用HSE06读取该结构做单点能带计算；最后对比PBE和HSE06的态密度，重点看费米面附近3 eV范围内的峰位和峰高。若HSE06下带隙打开但态密度峰形与PBE基本一致，说明赝势选择对电子结构影响较小；若峰位明显移动，则需检查赝势的芯区处理是否合理。

四、工程选型检查清单

实际项目中，建议按以下顺序锁定泛函与赝势组合，避免反复返工：

• 明确目标量：若只关心态密度峰形和轨道成分，PBE+PAW通常足够；若需要带隙绝对值，必须上HSE06。
• 检查赝势版本：同一元素可能有多个PAW数据集（如标准、硬、软），硬版本精度高但截断能大，软版本适合大体系快速计算。
• 做一次交叉验证：用PBE和HSE06分别计算同一结构的带隙，若差异超过0.3 eV，说明该体系对泛函敏感，需谨慎报告绝对带隙。
• 记录截断能与k点：赝势切换后，原截断能可能不再适用，需重新做收敛测试，确保总能量变化小于1 meV/atom。

总结来说，泛函决定电子结构的物理框架，赝势决定数值实现的精度上限。工程选型没有“最好”，只有“匹配目标”。对于能带和态密度计算，推荐以PBE+PAW为基线，HSE06+PAW为带隙修正，SCAN或DFT+U处理强关联，并在报告中明确标注所用泛函和赝势版本，以保证结果可复现。