计算化学公社

标题: vasp非自洽计算成功收敛的结构，进行自洽计算会出现非常大的原子受力 [打印本页]

作者
Author: Shana 时间: 2023-4-12 15:42
标题: vasp非自洽计算成功收敛的结构，进行自洽计算会出现非常大的原子受力
本帖最后由 Shana 于 2023-4-16 21:51 编辑

各位老师好，小弟使用vasp计算也已经快两年了，还是第一次碰到这种情况。最近在计算双钙钛矿Cs2AgBiCl6的声子谱，需要高精度优化结构，优化精度是1E-6 eV/A，截断能是ENCUT=700 eV, PREC = acc, ADDGRID=T，k点设置是gamma center 9*9*9，通过离子步迭代后，成功达到力的收敛精度，即每个原子受力小于1E-6 eV/A。但是我复制CONTCAR 为POSCAR，继续运行结构优化，第一个离子步的力如图片一所示，可以看到有几个特殊位点的受力远大于之前的收敛标准，想请问一下各位老师，这是怎么回事？我试过提高截断能，提高k点密度，仍然没有效果。
感谢各位老师的指导，谢谢！

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