静态自洽电子结构计算：从入门到绕坑的实战笔记

做第一性原理计算的人，大概没有谁绕得过静态自洽这一关。它看起来只是”VASP 输入文件里 ISTART=0.NSW=0.然后跑就是了”，但真正上手才发现——收敛判据设多少合适?k点密度怎么选?ENCUT 太小结果漂移，太大又跑不动。这些问题没人手把手教，只能在反复调参里慢慢摸出门道。

什么是静态自洽

所谓静态自洽(Self-Consistent Field，SCF)，本质是求解 Kohn-Sham 方程的迭代过程。体系总能量和电子密度相互依赖：给定初始电荷密度，解方程得到波函数，计算新电荷密度，再代回去重新求解——直到新旧密度之差低于设定阈值，SCF 收敛，体系达到基态。

这个过程是所有后续计算(能带、态密度、声子谱、光学性质)的起点。SCF 没跑稳，后面全是空中楼阁。

几个最容易卡住的地方

收敛困难是最常见的坑。 EDIFF 和 EDIFFG 两个参数控制收敛精度，默认值 1E-4 通常够用，但某些磁性体系或含过渡金属的体系需要更严格的标准。这里有个经验：如果能量在 1E-3 量级来回跳动，别急着调阈值，先检查 POTCAR 是否一致、INCAR 里有没有漏掉 ISMEAR 和 SIGMA(治金体系和半导体体系对 smearing 的敏感度完全不同)。

k点网格也是隐形的坑。 金属体系通常需要 0.03 Å⁻¹ 以上的密度，半导体和绝缘体可以放宽到 0.05.gamma-centered 和 Monkhorst-Pack 两种方法对立方晶系差别不大，但对六方或斜方体系，选错可能导致布里渊区采样不足，能带看起来歪歪扭扭。

关于 CHGCAR 和 WAVECAR 的理解也很关键。 静态自洽完成后，会生成 CHGCAR(电荷密度)和 WAVECAR(波函数)。前者用于后续非自洽计算，后者可以用于加速 AIMD 或其他需要良好初猜的静态计算。ISTART=1 时 VASP 会读取 WAVECAR 作为初猜——但这里有个陷阱：如果你的 INCAR 改了 ENMAX 或 ISIF，WAVECAR 的兼容性会被破坏，收敛行为完全不可预期。

实操建议

对于大多数体系，我的建议是：先跑一个低精度版本(ENCUT=400.k点稀疏)，确认结构没问题、体系没有对称性报错，再逐步加密参数。分两步走比一次性跑高精度失败后找不到问题根源要高效得多。

如果做的是磁性体系，别忘了设置 ISPIN=2 和 MAGMOM。初始磁矩的设置直接影响自旋极化方向的收敛方向——我见过不少人在 ISPIN=2 但 MAGMOM 全设成 0 的情况下跑出非磁性基态，完全偏离预期。

结语

静态自洽计算是看似简单、实则细节满满的环节。收敛判据、k点采样、smearing 参数，每一个选择背后都有物理依据。建议在正式计算前先找同类型体系的文献参数做参考，有个基准再微调，比闭着眼睛试错要节省大量机时。

如果你在静态自洽阶段遇到了奇怪的能量震荡或电荷密度不收敛，欢迎到科研学术网交流具体参数和体系背景，同样的问题往往有不同的解决路径。