计算资源有限的情况下（20核），先用gama点进行结构优化再用合理的K点跑是否可行

jbgande

如题，对80个以上的原子体系，直接用vasp官网手册推荐值331跑时间过长，经测试使用gama点跑的时间约为22小时，K点为221跑的时间为111小时，然后用331分别跑后的最低能量态有0.35ev的差别（先221后331所得能量更低为-459.72ev），结构也有些许不同。是否有必要先用221或直接用331进行结构优化。

nemolee

肯定是用vasp_gam版先跑啊，甚至如果资源这么有限的话就先往下做做，很多东西定性是问题不大的。

卡开发发

有必要看下能量和力收敛的情况再来考虑要不要降低精度。k点影响的不止是能量，还有力和电子结构不同，如果这样可能会导致构型优化的结果不同。如果极端情况下，gamma点计算相比最初没优化的结构偏离更高精度的结果更多，那么浪费的时间就更多了。

jbgande

卡开发发 发表于 2021-5-1 09:01
有必要看下能量和力收敛的情况再来考虑要不要降低精度。k点影响的不止是能量，还有力和电子结构不同，如果 ...

是的，我师兄也是这么说，不过用331对221或gama点优化后的结构进行再优化时，只用了10步左右。所以我尝试设置EDIFY =-0.02对先gama点后331的结果有一定优化，但也有0.25ev的差别。而且在gama点优化后的结构添加EDIFFG后跑了近400步才达到收敛。

jbgande

nemolee 发表于 2021-5-1 00:26
肯定是用vasp_gam版先跑啊，甚至如果资源这么有限的话就先往下做做，很多东西定性是问题不大的。

我也是这么想，就怕像下面那位仁兄说的对于电子结构影响太大。可能也需要进行这方面的测试再下结论。

卡开发发

jbgande 发表于 2021-5-1 16:32
是的，我师兄也是这么说，不过用331对221或gama点优化后的结构进行再优化时，只用了10步左右。所以我尝试 ...

这说法我没看明白，你gamma跑的时候不设置EDIFFG还是？
建议上最好k和动能截断都做好测试，体系稍大的话收敛性标准可以适当放宽一些，比如10～20meV/atom，但是要看看force和stress变化多大。实在不行可以试试一些LCAO程序，80原子应该能算得轻松些，平面波毕竟还是慢一些。

jbgande

卡开发发 发表于 2021-5-1 18:32
这说法我没看明白，你gamma跑的时候不设置EDIFFG还是？
建议上最好k和动能截断都做好测试，体系稍大的话 ...

好的，谢谢老师，我再测试一下，原来跑的时候ediffg用的就是默认值，后来尝试了较高的力收敛标准。

卡开发发

jbgande 发表于 2021-5-2 18:27
好的，谢谢老师，我再测试一下，原来跑的时候ediffg用的就是默认值，后来尝试了较高的力收敛标准。

默认EDIFFG是EDIFF*10，是以能量为收敛标准的，相对会比较宽松。