K点密度收敛测试标准方法：VASP总能收敛与能带精度的权衡

总能收敛到1 meV/atom以内了，K点测试就算通过了——这是VASP用户中最普遍的误解之一。能量对K点密度的收敛速度比电子结构快得多，当总能还在1 meV量级漂移的时候，能带可能还在0.1 eV量级变化。对于需要报告带隙到0.01 eV精度的项目来说，这个差距不容忽视。

能量收敛和电子结构收敛是两回事

一个简单的测试就能说明问题。用Si的PBE计算，从5×5×5的K点网格逐步加密到15×15×15：

总能从7×7×7开始就基本收敛了（变化<1 meV），但带隙从7×7×7到13×13×13还在持续变化，直到11×11×11以上才趋于稳定。如果只看总能收敛就停下，你的带隙可能差了0.04-0.05 eV——对Si来说这是约10%的相对误差。

金属体系更敏感。体心立方Fe的费米面附近态密度（DOS）对K点密度的要求在金属体系中出名的严格——用15×15×15可能还不足以精确描述费米面的细节。有些文献用了21×21×21甚至更高密度来保证磁矩和电子比热的收敛。

KSPACING还是手动K点网格？

VASP 5.4之后引入了KSPACING参数，自动根据倒格矢长度生成K点网格。它的优势是：不管你的晶胞是立方还是斜方，KSPACING都能给出密度均匀的网格。对比之下，手动设置KPOINTS文件时如果只写一个K×K×K，对于a≠b≠c的晶胞实际上采样密度在三个方向上是不均等的。

KSPACING的推荐取值：

但KSPACING也有局限。对于需要精确控制K点在倒空间分布的场合（比如计算能带时沿特定高对称路径），手动指定KPOINTS仍然不可替代。同时，KSPACING生成的网格在某些奇数网格（如1×3×5）可能不对称，需要在OUTCAR里检查最终生成的K点分布。

折中方案是：SCF计算用KSPACING做收敛测试确定最佳密度，能带/态密度计算手动写KPOINTS沿特定路径。

正确的收敛测试流程

不是把K点从3×3×3一路加到21×21×21，那是暴力计算。有效率的做法是：

第一步：从KSPACING=0.5开始，每步除以√2。 序列是0.5 → 0.35 → 0.25 → 0.18 → 0.13 Å⁻¹。这个等比递进保证每次变化量递减，不会浪费计算资源。

第二步：监测总能和带隙两条收敛曲线。 单独看总能收敛不够。带隙（或费米能级处的DOS值）是更严格的判据。

第三步：确定收敛判据。 对于常规能带计算，带隙变化<0.02 eV即可接受；对于金属，费米能级处DOS变化<5%即可接受；对于需要报告带隙精度到0.01 eV的工作，带隙判据收紧到<0.005 eV，对应K点密度可能需要翻倍。

第四步：对收敛的网格做一次+1验证。 比如你认为11×11×11已经收敛，就用13×13×13验证：如果总能差<1 meV/atom且带隙变化<0.01 eV，确认收敛。

特殊体系的K点陷阱

低维材料（2D/1D）： 真空层方向的K点只需设为1——没有必要在真空方向采样。KPOINTS写成N×N×1（2D）或N×1×1（1D）。这个看似简单的规则被忽略的频率超乎想象——很多算石墨烯和MoS₂的文章，真空方向K点写成了和面内一样密，白白浪费了几倍的计算量。

超胞计算： 做了2×2×2超胞后，K点可以从原胞的8×8×8降到4×4×4——超胞的倒格矢缩小了对应倍数，K点密度的绝对分辨率保持一致。但如果超胞后仍然用原来的K点密度，问题也不大（精度更高），只是计算量不必要地增加了。

六方晶系： 对于hcp、wurtzite等六方结构，Gamma-centered网格（KPOINTS首字母为G）比Monkhorst-Pack（首字母为M）更适合——因为六方布里渊区的高对称点和Gamma-centered网格的分布更一致。选错的话能带收敛会慢一个数量级。

K点密度的选择在很大程度上决定了计算结果的可信度。科研学术网（https://www.keyanxueshu.com）VASP计算的标准化参数库和收敛测试模板，覆盖从晶体到低维材料的常见计算场景。