计算化学公社

标题: 求助 VASP 优化结构是否收敛 [打印本页]

作者
Author: yaoyaoy 时间: 2016-9-29 15:17
标题: 求助 VASP 优化结构是否收敛
本帖最后由 yaoyaoy 于 2016-9-29 15:24 编辑

大家好，有一个问题想请大家帮忙。
我现在用VASP 优化结构，虽然显示reached required accuracy - stopping structural energy minimisation，但是从total force 看唯一优化的原子的受力是负值，如图，
VASP默认的是每个原子上的受力小于0.02eV/A, 在OUTCAR中显示的total force 会有负值，这个负值代表什么意思呢，要如何与0.02eV/A比较呢？那这个结构算是优化完成了吗？

作者
Author: 卡开发发 时间: 2016-9-29 18:26
这个就是该认真看说明书了。VASP默认是EDIFFG = EDIFF*10，其中EDIFFG defines the break condition for the ionic relaxation loop，默认的收敛标准不是力， If EDIFFG is negative it has a different meaning: In this case the relaxation will stop if all forces are smaller than | EDIFFG |，显然，力应当比较绝对值。这些洋码都是来自说明书。

作者
Author: yaoyaoy 时间: 2016-9-29 19:08

卡开发发发表于 2016-9-29 18:26
这个就是该认真看说明书了。VASP默认是EDIFFG = EDIFF*10，其中EDIFFG defines the break condition for th ...

谢谢

那像x方向的total force = -0.026678中的负号是代表受x轴负方向的力吗？

作者
Author: 卡开发发 时间: 2016-9-29 22:57

yaoyaoy 发表于 2016-9-29 19:08
谢谢
那像x方向的total force = -0.026678中的负号是代表受x轴负方向的力吗？

对

作者
Author: yaoyaoy 时间: 2016-9-30 10:30

卡开发发发表于 2016-9-29 22:57
对

多谢您的回答，谢谢

还有一个问题想问您，就是现在优化的结构的结果显示reached required accuracy - stopping structural energy minimisation，但是能量的变化和受力明明没有收敛，这是为什么呢？有什么方法可以调整使得这两个结果收敛吗？

作者
Author: 卡开发发 时间: 2016-9-30 17:45

yaoyaoy 发表于 2016-9-30 10:30
多谢您的回答，谢谢
还有一个问题想问您，就是现在优化的结构的结果显示reached required accuracy - ...

应该不会吧，如果EDIFFG设置的是正值的话应该per atom的能量差要达到手里标准才会结束计算吧。

作者
Author: yaoyaoy 时间: 2016-9-30 20:21

卡开发发发表于 2016-9-30 17:45
应该不会吧，如果EDIFFG设置的是正值的话应该per atom的能量差要达到手里标准才会结束计算吧。

优化的时候只有一个吸附的原子的z坐标在优化，但是得到的z方向的total force = 0.052，而且能量变化也没有达到标准，这是为什么呢？麻烦您给解答一下，谢谢。

作者
Author: 卡开发发 时间: 2016-9-30 23:35

yaoyaoy 发表于 2016-9-30 20:21
优化的时候只有一个吸附的原子的z坐标在优化，但是得到的z方向的total force = 0.052，而且能量变化也没 ...

EDIFFG设置正值的话就和Force的收敛判据没关系了，能量则按照per atom的变化来计算啊，不知道你的体系几个原子？不管是否固定，总能的贡献都来自于电子能和核排斥。

作者
Author: yaoyaoy 时间: 2016-10-1 09:14
本帖最后由 yaoyaoy 于 2016-10-1 10:30 编辑

卡开发发发表于 2016-9-30 23:35
EDIFFG设置正值的话就和Force的收敛判据没关系了，能量则按照per atom的变化来计算啊，不知道你的体系几 ...

原来能量要看每个原子的能量的变化，，之前我一直在关注total free energy的值的变化，那上面粘贴的结构优化结束得到的 free energy TOTEN 值的变化没有小于EDIFFG=1E - 04，这个是不是和手册里写的 ‘If the change in the total (free) energy is smaller than EDIFFG between two ionic steps relaxation will be stopped.’ 矛盾了呢？

作者
Author: 卡开发发 时间: 2016-10-1 14:06

yaoyaoy 发表于 2016-10-1 09:14
原来能量要看每个原子的能量的变化，，之前我一直在关注total free energy的值的变化，那上面粘贴的结构 ...

我不知道你的系统有多少原子数目啊，其实你也可以想到肯定是要除掉atom number的，否则不同大小的系统计算精度也就无法比较了。

作者
Author: yaoyaoy 时间: 2016-10-9 14:49

卡开发发发表于 2016-10-1 14:06
我不知道你的系统有多少原子数目啊，其实你也可以想到肯定是要除掉atom number的，否则不同大小的系统计 ...

嗯，除去原子数目以后达到了精度。谢谢您的帮助。

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