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

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

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

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

卡开发发 发表于 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.’ 矛盾了呢？

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

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

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

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

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

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

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

对

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

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

这个就是该认真看说明书了。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 |，显然，力应当比较绝对值。这些洋码都是来自说明书。