用不同的ALGO参数优化出来的结构能量差别很大。

sunhao

各位老师好，我想请教两个问题。
（1）我用vasp计算的体系有磁性，并考虑了U值，由于结构很难收敛，调整ALGO=Fast，ALGO=Normal，ALGO=Very_Fat，但是发现他们都能收敛时输出的能量差别非常大，不知道用哪个能量。
（2）优化好的结构做静态计算不收敛，或者是收敛了但和优化结构输出的能量不一致，差别较大。

卡开发发

1、

结构很难收敛

调整ALGO没有用处啊，ALGO是控制SCF中对角化的方法的。很难收敛的情况RMM-DIIS确实可能不稳定，相对而言Davidson对角化更好一些。DFT+U的情况LMAXMIX尽量设置一下。
2、检查下结构弛豫是否每一步SCF都收敛了，如果是这样，也尽量多cp CONTCAR POSCAR算几次。

sunhao

谢谢老师的解答。我先按照您说的试试。非常感谢。

sunhao

卡开发发 发表于 2019-11-30 03:25
1、调整ALGO没有用处啊，ALGO是控制SCF中对角化的方法的。很难收敛的情况RMM-DIIS确实可能不稳定，相对而言 ...

老师你好，我调整了一下初始结构，现在收敛起来比较容易。但是优化好的结构做静态计算，输出的能量差别还是很大。而且用优化好的结构，把CONTCAR改成POSCAR，其它什么都不变，提交任务后并不是计算一步就达到力的收敛标准，任务结束了，而是跑了很多步，且能量与原来也相差很大。请问这个问题是不是还是初始结构没有优化好造成的？应该怎么解决啊？

sunhao

卡开发发 发表于 2019-11-30 03:25
1、调整ALGO没有用处啊，ALGO是控制SCF中对角化的方法的。很难收敛的情况RMM-DIIS确实可能不稳定，相对而言 ...

这个是我的INCAR，您看看有什么问题没

卡开发发

sunhao 发表于 2019-12-2 09:06
老师你好，我调整了一下初始结构，现在收敛起来比较容易。但是优化好的结构做静态计算，输出的能量差别还 ...

而是跑了很多步，且能量与原来也相差很大

1、做好动能截断的收敛性测试。
2、在1基础上多拷贝几次CONTCAR为POSCAR。
3、必须保证每一步SCF都收敛了。

ghifi37

sunhao 发表于 2019-12-2 09:06
老师你好，我调整了一下初始结构，现在收敛起来比较容易。但是优化好的结构做静态计算，输出的能量差别还 ...

VASP在静态计算中第一步基本上得不到以前收敛的能量，哪怕一点东西都不改。不知道这算不算个bug

sunhao

卡开发发 发表于 2019-12-2 12:25
1、做好动能截断的收敛性测试。
2、在1基础上多拷贝几次CONTCAR为POSCAR。
3、必须保证每一步SCF都收 ...

老师您说的动能截断的收敛测试就是用不同的ENCUT值去优化结构吧？随着ENCUT增大，相邻两个ENCUT值优化出来的能量差ΔE小于多少时，算是收敛了呢？

sunhao

ghifi37 发表于 2019-12-2 15:23
VASP在静态计算中第一步基本上得不到以前收敛的能量，哪怕一点东西都不改。不知道这算不算个bug

静态自洽不就只有一步吗？结构优化的能量跟静态自洽的能量有一点点差别我觉得还可以接受吧，但我的这个差别在0.5eV左右，所以算出来的能量根本没法用。我觉得应该是结构优化没有优化好。我正在调整ENCUT参数测试。

sunhao

ghifi37 发表于 2019-12-2 15:23
VASP在静态计算中第一步基本上得不到以前收敛的能量，哪怕一点东西都不改。不知道这算不算个bug

静态自洽不就只有一步吗？结构优化的能量跟静态自洽的能量有一点点差别我觉得还可以接受吧，但我的这个差别在0.5eV左右，所以算出来的能量根本没法用。我觉得应该是结构优化没有优化好。我正在调整ENCUT参数测试。

waitingseven

simga=0.2是不是有点儿大了，设置成0.05试试，相应的也可以调整一下kpoints，还有既然有磁性，电子密度mix参数是不是应该设置一下

卡开发发

sunhao 发表于 2019-12-2 16:21
静态自洽不就只有一步吗？结构优化的能量跟静态自洽的能量有一点点差别我觉得还可以接受吧，但我的这个差 ...

1、收敛性测试手册上提到过标准（1meV/atom，50eV间隔），但也看你的体系尺寸大小和计算能力。
2、变胞优化的时候，FFT网格尺寸应该是不调整的，所以这应该是导致晶格优化和自洽结果不同的一个原因，这也是上面提到最好cp CONTCAR POSCAR的原因，动能截断比较充分的时候问题会不太显著。不管怎样，收敛性测试应该还是要做好，因为会影响到force、stress计算的质量，进一步影响到构型优化的稳定性和效率。
3、ISPIN=2的时候尽量做好MAGMOM的初猜。

sunhao

waitingseven 发表于 2019-12-2 17:29
simga=0.2是不是有点儿大了，设置成0.05试试，相应的也可以调整一下kpoints，还有既然有磁性，电子密度mix ...

好的，谢谢老师。0.2是有点大，但是我的体系比较大，查看了一下OUTCAR里面的entropy T*S这一项，除以原子个数，也小于0.001。后面我在减小sigma计算一下对比对比。mix参数是什么啊？之前没有用过，我看看vasp手册再跟您讨论吧，谢了。

sunhao

卡开发发 发表于 2019-12-2 18:21
1、收敛性测试手册上提到过标准（1meV/atom，50eV间隔），但也看你的体系尺寸大小和计算能力。
2、变胞 ...

1.我用您说的方法做了收敛性测试，发现ENCUT需要在800时才基本满足1meV/atom，大多数人变胞优化ENCUT用到700就够了，我这种情况正常吗？
2.如果我在后面做吸附的时候，是都按ENCUT=800固定晶格参数计算，还是说可以把吸附前的结构固定晶格参数用ENCUT=500再优化一下，在此基础上，后面的吸附计算都用ENCUT=500呢？
多谢老师指导！

卡开发发

sunhao 发表于 2019-12-4 07:59
1.我用您说的方法做了收敛性测试，发现ENCUT需要在800时才基本满足1meV/atom，大多数人变胞优化ENCUT用到 ...

1、后续如果吸附的体系比较大的话可以适当放宽一点。
2、做收敛性测试的目的之一就是保证能量和力计算可靠。从而构型优化更加顺利，结构优化应当使用充分一些的动能截断。