计算化学公社

标题: 通过VASP的CINEB计算扩散能垒，不收敛的问题 [打印本页]

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

作者
Author: UW_0728. 时间: 前天 20:49
本帖最后由 UW_0728. 于 2026-1-21 21:14 编辑

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

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

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

作者
Author: liuzl 时间: 5 hour ago

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

谢谢回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)