材料能带DFT计算：带隙预测与缺陷态分析的工程实践

材料能带DFT计算做半导体方向的人基本都跑过，但跑到能直接拿去指导实验的程度，中间的弯弯绕绕比大多数人预想的多。一个二硫化钼场效应管沟道材料项目，目标是用DFT预测单层和双层二硫化钼的带隙大小，以及应变调控对带隙的调制幅度。项目启动时给的截止日期是三周，实际用了七周，主要时间都耗在带隙算不准这个问题上。

PBE带隙低估：修正策略的选择

PBE对二硫化钼单层的带隙预测是1.4 eV，实验结果约1.8到1.9 eV，低估了0.4到0.5 eV。这个偏差对于应变能不能把带隙调到目标区间这个工程判断有直接影响。项目在这里的抉择是：不上杂化泛函，改用剪刀算符修正——把PBE的能带结构沿能量轴剪开一段距离，让带隙匹配实验值，然后在修正后的能带基础上做应变响应计算。

这个方法在文献里争议很大，但在工程实践中，如果你的核心问题不是绝对带隙是多少而是带隙随应变的变化趋势，剪刀算符修正的误差通常是可接受的。这个项目的结论是双轴拉伸应变到5%时带隙缩小约0.15 eV，直接和间接带隙之间的竞争导致应变超过3%后间接到直接转变消失——这个结论在后续实验中得到了验证。

缺陷能级：电荷修正不能省

二硫化钼在实际器件中总有硫空位，这个空位在能带隙中引入缺陷态，影响载流子迁移率。计算缺陷能级的标准方法是计算带缺陷超胞的总能，然后用公式推导缺陷形成能和能级位置。但这里有个经典难题：周期边界条件下点缺陷的带电超胞计算会遇到静电相互作用和参考势能漂移。

这个项目用的是Freysoldt电荷修正方法，对带电缺陷超胞的总能进行长程静电修正。修正量在0.1到0.3 eV之间，对于缺陷能级是否在带隙上半部这个判断，修正前后的结论可能完全相反。材料能带DFT计算做到缺陷态这一步，精度控制的标准要比纯能带计算高一个档次。

回过头梳理，这个项目最值得记录的经验是：能带计算不是跑一个标准的SCF加NSCF流程就结束了。你要问自己这个带隙数值要用来做什么判断，然后再决定精度投资应该投在哪个环节。