费米能级计算：DFT框架下的电子结构精度与参数策略

费米能级计算决定了材料电子行为的基准线。这个项目在研究n型氧化锌半导体的掺杂特性时，费米能级相对于导带底的位置直接关系到载流子浓度和导电类型的判断。项目组需要通过VASP软件完成电子结构计算，精确定位费米能级，而其中的参数选择远比教科书描述的复杂。

费米能级计算的首要前提，是布里渊区积分的充分收敛。项目组在ZnO原胞的计算中发现，使用4×4×4的Monkhorst-Pack K点网格时，费米能级位于-3.27 eV（以真空能级为参考），当网格加密到12×12×12时，费米能级稳定在-3.41 eV，变化幅度达到0.14 eV。对于金属体系，这种敏感性更加剧烈——项目组在铜的费米能级测试中，K点从8×8×8增至20×20×20，费米能级漂移了0.08 eV后才趋于收敛。

差距不会说谎。项目组同时测试了Tetrahedron方法和Gaussian展宽方法对费米能级的影响。在半导体ZnO中，Tetrahedron方法给出的费米能级比Gaussian展宽（σ=0.05 eV）低0.03 eV，差异不大。但在金属铜中，两种方法的差异扩大到0.06 eV。项目组认定，对于金属体系的费米能级计算，Tetrahedron with Blöchl corrections是更可靠的选择，而半导体体系使用Methfessel-Paxton展宽（σ=0.1 eV）即可满足精度需求。

费米能级计算中最根本的挑战，来自DFT对带隙的系统性低估。项目组在计算ZnO时，PBE泛函给出的带隙仅为0.73 eV，而实验值为3.37 eV。这意味着PBE计算得到的费米能级相对于导带底的位置严重偏离实际。对于n型掺杂的ZnO，PBE预测费米能级进入导带，呈现出简并半导体特征，但实际上掺杂浓度不足以使费米能级穿越带隙。

项目组转向了HSE06杂化泛函。HSE06将ZnO带隙修正至2.85 eV，虽然仍低于实验值，但费米能级的定位已经合理——在Al掺杂浓度为1×10¹⁸ cm⁻³时，费米能级位于导带底下方0.15 eV，与Hall效应实验测量的载流子浓度吻合。这一结果与Ceder课题组在Physical Review B上发表的系统性研究一致，即HSE06对宽禁带半导体的费米能级定位具有明显改善。

对于计算资源不足以支撑HSE06的体系，项目组验证了DFT+U方法作为替代方案。对Zn的3d轨道施加U=10 eV的Hubbard参数后，带隙提升至1.60 eV，虽然仍不及实验值，但费米能级相对于价带顶的位置已经修正到合理范围。值得注意的是，U值的选择必须基于物理论证而非简单拟合——项目组采用线性响应方法自洽确定U值，避免了人为调参引入的非物理偏移。

费米能级计算不仅是定位一个能量值，更是理解费米面附近电子态的分布。项目组在分析铜的费米面时，发现PBE计算的态密度在费米能级处存在明显的尖峰，对应于sp带与d带的交叉区域。通过投影态密度分析，确认d带中心位于费米能级下方2.1 eV处，与Hammer-Nørskov的d带中心模型给出的2.0 eV高度一致。

这个参数在后续的催化活性预测中发挥了关键作用。项目组根据d带中心位置预测了铜表面CO吸附强度，理论值与实验值的偏差控制在0.1 eV以内。需要强调的是，投影态密度的精度同样依赖于K点密度，项目组在费米面附近将K点加密至24×24×24后，d带中心的收敛精度达到了0.02 eV量级，这一精度足以支撑后续的催化剂筛选工作。

费米能级计算必须坦承自身的局限。标准DFT无法描述强关联体系中的Mott绝缘体行为——项目组在尝试计算NiO的费米能级时，PBE预测其为金属，而实验上NiO是带隙为4.0 eV的绝缘体。即使采用DFT+U（U=5.3 eV），费米能级的绝对位置仍存在0.3 eV的不确定性。对于这类体系，DFT+DMFT或GW近似才是更恰当的选择。

此外，费米能级的绝对值参考点选择也值得警醒。项目组在计算表面体系的功函数时，需要将费米能级与真空能级对齐。在有限平板模型中，偶极校正的引入会使费米能级偏移0.1-0.2 eV。被证明有效的做法是，在平板计算中始终开启偶极校正（`IDIPOL=3, LDIPOL=.TRUE.`），并在功函数报告中明确标注校正前后的差异。这套规范在项目组后续的六个表面体系中均得到了验证。