通过VASP的CINEB计算扩散能垒，不收敛的问题

liuzl · 发表于 Post on 3 day ago

计算的体系包括240个原子，96核计算，中间取点6个，一周多还没收敛，我通过脚本提取了其中一个点的能量，发现计算过程中有一个离子步的能量已经低于我设置的力收敛值0.05，于是我就停了，将所有点的CONTCAR复制为POSCAR，重新开始算，可是还是在震荡，这种情况合理吗？（其他五个点还没有达到力的收敛值，也存在力的震荡）

UW_0728. · 发表于 Post on 前天 20:49

本帖最后由 UW_0728. 于 2026-1-21 21:14 编辑

解决办法应该类似于解决几何优化不收敛的办法

比如检查SCF收敛精度是否足够（默认的1e-4我觉得太宽松了，怎么也得1e-5甚至更严）、积分格点精度是否足够（有没有用PREC = Accurate，这是几何优化强烈建议用的；以及必要时可以加ADDGRID=.TRUE.）、是否恰当考虑k点以及k点数目（这个不太重要，几何优化任务对k点敏感度不高）等等

另外，240个原子的体系用VASP算太慢、太费劲，推荐用CP2K替代

liuzl · 发表于 Post on 12 hour ago

UW_0728. 发表于 2026-1-21 20:49
解决办法应该类似于解决几何优化不收敛的办法

谢谢回复

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 通过VASP的CINEB计算扩散能垒，不收敛的问题

评分 Rate

浏览过的版块