电子结构计算晶体：VASP能带与态密度计算的参数收敛实战

一篇文献报道了二维MoS₂单层的带隙为1.89eV，另一篇说是1.78eV，还有一篇给了1.92eV。同一个材料，三组数据，差距超过0.1eV。对于带隙这个量级的物理量来说，0.1eV的偏差已经足以改变对材料光电性能的判断。项目组在复现过程中发现，问题的根源不在VASP本身，而在于电子结构计算晶体时的参数收敛是否做到了位。

ENCUT：不是越大越好，而是要找到”够了”的那个值

电子结构计算的第一步是确定平面波截断能ENCUT。VASP默认的ENCUT值通常是POTCAR中ENMAX的1.3倍，对于大多数元素来说在300到400eV之间。但”默认”不代表”够了”。

项目组对MoS₂单层做了ENCUT从300eV到700eV的收敛测试，步长50eV。以700eV的结果为参考基准，总能量的收敛情况如下：300eV时偏差达0.23eV/atom，这在表面能计算中是不可接受的；400eV时偏差降到0.05eV/atom，已经进入化学精度范围（约1kcal/mol）；500eV时偏差进一步降到0.008eV/atom。

但如果只看总能量收敛就确定ENCUT，会错过一个关键问题：带隙对ENCUT的收敛速度比总能量慢得多。在400eV下，PBE泛函给出的带隙是1.72eV；在500eV下是1.78eV；到600eV才稳定在1.80eV。也就是说，总能量在400eV就已经收敛了，但带隙需要到600eV才趋于稳定。

这个现象的原因是带隙涉及导带底和价带顶的能级差，这两个能级对基组的完备性比总能量更敏感。导带底附近的态通常更加弥散，需要更高的平面波截断才能准确描述。项目组的经验是：电子结构计算晶体的ENCUT不能只看总能量收敛，必须以带隙（或你所关心的具体电子结构特征）的收敛作为判断标准。对于含3d过渡金属的体系，ENCUT建议不低于520eV；对于只含轻元素的体系，450eV通常是足够的。

KPOINTS：Γ中心还是Monkhorst-Pack？

k点网格的选择是电子结构计算晶体时的第二个关键参数。对于MoS₂这种六方晶系，很多人习惯用Γ中心的网格，因为高对称点Γ在倒空间的原点，直觉上觉得这样更”对称”。但项目组的测试结果推翻了这个直觉。

在相同的k点密度下（以倒空间分辨率0.02Å⁻¹为基准，对应约12×12×1的网格），Γ中心网格给出的带隙是1.76eV，而Monkhorst-Pack网格给出的带隙是1.80eV。后者更接近实验报道的1.88eV（考虑PBE通常低估带隙约0.1eV）。

差异的来源是六方晶系的布里渊区采样。Γ中心网格在六方晶系中会在布里渊区边界附近产生不对称的采样点分布，导致某些高对称方向上的能带计算出现系统偏差。Monkhorst-Pack网格的采样点在倒空间中是均匀分布的，对于六方晶系更加合适。一个实用的判断规则是：如果晶格不是完全的立方对称，优先使用Monkhorst-Pack网格，除非你要计算的物理量对Γ点有特殊敏感性（比如光学性质中的直接跃迁）。

能带计算中的”幽灵”能级

完成ENCUT=600eV、12×12×1 Monkhorst-Pack网格的静态计算后，项目组提取了能带结构。在Γ-K方向的高对称线上，发现了一个意外的平带，位于费米能级以上约0.5eV处。这个平带在文献中从未被报道过。

排查过程花了两天。先怀疑是真空层不够厚——MoS₂单层的周期性计算需要在z方向加足够的真空层来隔离相邻层间的相互作用。当前用的是15Å，增加到20Å后平带依然存在。然后怀疑是赝势问题——换了GW势（PAW_PBE_GW）重新计算，平带消失了。

结论是：标准PBE赝势在处理Mo的4d电子时，在半芯态的描述上存在不足，导致出现了一个虚假的能级。对于含4d/5d过渡金属的体系，如果遇到无法解释的能带特征，换用GW势或包含半芯态的更硬赝势重新验证是必要的。这个问题在轻元素体系中几乎不会遇到，但在Mo、W、Nb等元素的电子结构计算中并不罕见。

态密度的”美化”陷阱

态密度图是最容易被”美化”的电子结构分析结果。项目组注意到很多文献中的DOS曲线极其光滑，而原始数据其实有大量噪声。平滑处理本身没有问题——通常用高斯展宽，展宽参数SIGMA取0.05到0.10eV是合理的——但展宽过大会抹掉真实的物理特征。

以MoS₂为例，在SIGMA=0.05eV时，价带顶附近可以清晰看到两个分立的峰，分别对应Mo的4d轨道与S的3p轨道的杂化态。当SIGMA增大到0.20eV时，这两个峰合并成了一个宽峰，完全丢失了轨道杂化的细节信息。如果在论文中使用SIGMA=0.20eV的DOS图，不仅无法分辨轨道成分，甚至可能得出错误的成键分析结论。

项目组的建议是：DOS计算的SIGMA不应超过0.10eV，如果噪声仍然太大，应该增加k点密度而不是增大展宽。k点从12×12增加到18×18后，SIGMA=0.05eV下的DOS曲线就已经足够平滑，不需要牺牲分辨率来换取光滑的曲线。

VASP电子结构计算晶体本质上是一系列参数收敛测试的系统工程。ENCUT要验证到带隙收敛而不只是总能量收敛，k点要关注倒空间采样均匀性而非单纯的密度数值，赝势的选择在过渡金属体系中不能图省事。这些环节中的任何一步偷工减料，最终的能带和态密度都可能是一个看似漂亮但偏离物理真实的计算结果。

更多VASP计算技术案例可参考科研学术网VASP计算栏目，欢迎访问科研学术网首页获取DFT计算服务支持。