vasp理论计算：从INCAR参数收敛到电子结构解读的一套交付标准

vasp理论计算交付给合作方之前，我必做的三件事是：确认截断能收敛、确认k点收敛、确认力与应力同时收敛。这套标准不是凭经验拍脑袋，是早年一个MoS2项目因为k点取稀了导致价带顶劈裂看错，被合作方打回重算后定下的规矩。从那以后，参数收敛检查成了交付清单的第一栏。

我的INCAR为什么先定ENCUT再定k点

参数收敛有顺序：先能量截断，再k点密度，最后smearing和收敛阈值，因为ENCUT错了后面全白搭。MoS2体系里Mo的PAW赝势ENMAX≈400eV，S的ENMAX≈270eV。我按惯例取ENCUT=1.3×max(ENMAX)=520eV，然后在400/450/520/600eV四个值上做单点能扫描。

结果显示：400→450eV总能下降0.08eV/原子，450→520eV只降0.012eV/原子，520→600eV降0.003eV/原子。按”相邻两个点能量差<0.01eV/原子”的收敛判据，520eV已经稳。这里有个可复用要点：ENCUT不需要盲目取高，1.3倍ENMAX是经验下界，真正决定取值的是能量差收敛曲线，不是赝势手册上的推荐值。

ENCUT还有个隐藏雷区：它对晶胞体积敏感。做ISIF=3完全弛豫时，若ENCUT偏低，程序会用”软赝势”把原子推到不真实的小体积来降总能，优化出的晶格常数能偏小1-2%。我交付钙钛矿BaTiO3时专门对比过：ENCUT=400eV时a轴3.94Å，520eV时3.97Å，实验3.99Å——差的那0.02Å就是截断误差。所以收敛扫描不能只在固定体积下单点能扫，最好连体积一起验证。可复用要点升级成一句：ENCUT的收敛判据要在弛豫后的体积上复算一遍，否则结构参数会偷偷跑偏。

k点密度是把双刃剑

定好ENCUT=520eV后做k点收敛。MoS2单层在xy平面周期、z方向真空15Å（ISIF=3优化时真空层不参与应力收敛），k点用Γ-centered：3×3×1、5×5×1、7×7×1、9×9×1四档。总能随k点变化：3×3→5×5降0.04eV，5×5→7×7降0.015eV，7×7→9×9降0.004eV。7×7×1达标。

但电子结构计算（能带）不能只看总能收敛。MoS2的价带顶在K点、导带底也在K点（直接带隙），K点附近能带色散对k点密度敏感。我把高对称路径K-Γ-M-K上的能带对照7×7和9×9的结果，带隙差0.03eV，可接受。这里的高下立判：结构优化用7×7×1够，但能带计算我仍用9×9×1配合line-mode KPOINTS，因为色散精度比总能更吃k点。关于k点与截断能的协同收敛，我也会参考vasp理论计算里对二维材料收敛曲线的讨论来交叉验证。

line-mode能带的具体画法有个容易翻车的点：高对称路径必须连续且覆盖价带顶/导带底所在的k点。MoS2的带边都在K点，路径若只画Γ-M-K而漏掉K-Γ段，就会看不到直接带隙。我现在固定用VASPKIT的能带路径自动生成，它按晶体对称性挑高对称点和路径，避免人为漏点。另外k点密度在line-mode里也要够：我习惯在每段路径上取至少40个点，否则色散曲线被拉成折线，有效质量算出来能差一倍。这条经验在狄拉克材料、拓扑绝缘体这类带边尖锐的体系更致命。

磁性体系里自旋与U不能省

MoS2本身是非磁的，但我并行算了一个1T’相MoS2（拓扑相变后含Mo的局域d电子），这个体系必须开ISPIN=2，且Mo的d电子有自相互作用误差，要加LDA+U。U(Mo)=3.0eV是我扫描2.0/3.0/4.0eV后定的：2.0eV时带隙0.9eV偏窄，4.0eV时1.4eV偏宽，实验ARPES给1.1eV，3.0eV吻合。

自旋初始化用MAGMOM显式给定（Mo 2.0，S 0.0），否则程序从顺磁初态出发容易收敛到非物理的低磁矩态。我踩过的坑是某次没写MAGMOM，Mo的磁矩收敛到0.3μB，带隙比应有值宽了0.2eV——整套电子结构全偏。磁性体系的U值和初始磁矩，必须写死在INCAR里，这不是可选项。

电子结构解读从DOS开始而非能带

交付电子结构时，我习惯先画总态密度（DOS）再画能带，因为DOS能立刻暴露问题是来自价带还是导带、是金属性还是带隙异常。MoS2的DOS显示价带顶由S-3p主导、导带底由Mo-4d(z²)主导，这与能带的反转一致——直接带隙且带边轨道成分明确，是后续载流子有效质量计算的前提。

PDOS进一步显示Mo的d轨道在导带底有强峰，说明电子有效质量小、迁移率高，与实验场效应迁移率约200 cm²/V·s定性相符。这一步的解读要点是：能带只给能量位置，PDOS才给轨道来源，二者结合才能把”计算带隙”翻译成”可理解的物理”。

把DOS能带翻译成有效质量时，我通常用能带曲率μ⁻¹∝∂²E/∂k²在带边取值。MoS2 K点直接带隙处的电子有效质量约0.45m₀，和ARPES实验的0.4-0.5m₀对得上——这是DFT能给到的”可实验验证”的硬指标。但要提醒边界：PBE的带隙系统性偏低（MoS2 PBE给1.6eV，实验2.0eV左右），做输运讨论时若直接用PBE带边曲率，迁移率会偏乐观。我的做法是带隙用HSE06或GW修正（scissor修正也行），但有效质量的曲率形状仍取自PBE，因为PBE的带形通常是对的，只是整体下移了带边。

三个交付前必查的收敛指标

第一，EDIFF与EDIFFG分开看。电子步EDIFF=1e-5保证总能稳定，离子步EDIFFG=-0.02 eV/Å保证力收敛。只查能量不查力，优化后的结构可能还差0.05Å没到位。

第二，应力收敛。二维材料ISIF=3优化时，真空层方向的应力必须接近0（<1 kB），否则真空层被压缩、层间耦合被引入假象。

第三，赝势版本一致。同项目里Mo和S的PAW必须用同一套（如vasp.5.4的Mo_sv、S默认），混用不同版本PAW会导致总能不可比。我每次交付都附一张赝势版本清单。

回过头看，理论计算交付的底线

回过头看，vasp理论计算交付的底线不是”跑完了”，而是”每一步都能量化重复”。截断能、k点的收敛曲线，U值的扫描结果，磁矩的初始化值——这些不是过程记录，是可重复性的担保。被证明有效的是把收敛检查做成交付前的强制关卡，而不是算完才补。

可复用要点总结：ENCUT看能量差收敛曲线定，k点在总能和电子结构两个层面分别验证，磁性体系U和MAGMOM写死，DOS先于能带解读。这套标准在二维材料、钙钛矿、氧化物三种体系上都站得住，它的适用边界只在强关联体系（如Mott绝缘体）需要扩展——那时HSE06或GW才是下一层，PBE+U已经触到天花板。这套vasp理论计算的收敛关卡，后来我移植到了所有第一性原理交付里。

交付时我还会附一张”计算成本表”：每步的核时、内存、耗时，让合作方知道这个带隙值花了多少算力。有次一个钙钛矿项目客户嫌慢，我拿出成本表说明HSE06带隙虽准但贵20倍，客户权衡后选PBE+scissor，省下的机时做了三组掺杂对比。把成本透明化，反而让合作更顺。

理论计算交付到我手里，早就不是”跑完发图”。收敛曲线、U扫描、磁矩、成本表，四件套齐了客户才放心。这套标准看着繁琐，但它让”可重复”从口号变成动作——这正是第一性原理计算最该交付的东西。