计算化学公社

标题: 用不同的ALGO参数优化出来的结构能量差别很大。 [打印本页]

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

作者
Author: 卡开发发 时间: 2019-11-30 03:25
1、

结构很难收敛

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

作者
Author: sunhao 时间: 2019-11-30 19:31
谢谢老师的解答。我先按照您说的试试。非常感谢。

作者
Author: sunhao 时间: 2019-12-2 09:06

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

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

作者
Author: sunhao 时间: 2019-12-2 09:22

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

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

作者
Author: 卡开发发 时间: 2019-12-2 12:25
本帖最后由卡开发发于 2019-12-2 12:26 编辑

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

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

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

作者
Author: ghifi37 时间: 2019-12-2 15:23

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

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

作者
Author: sunhao 时间: 2019-12-2 16:08

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

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

作者
Author: sunhao 时间: 2019-12-2 16:21

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

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

作者
Author: sunhao 时间: 2019-12-2 16:22

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

作者
Author: waitingseven 时间: 2019-12-2 17:29
simga=0.2是不是有点儿大了，设置成0.05试试，相应的也可以调整一下kpoints，还有既然有磁性，电子密度mix参数是不是应该设置一下

作者
Author: 卡开发发 时间: 2019-12-2 18:21

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

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

作者
Author: sunhao 时间: 2019-12-2 20:06

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

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

作者
Author: sunhao 时间: 2019-12-4 07:59

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

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

作者
Author: 卡开发发 时间: 2019-12-4 08:50

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

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

作者
Author: sunhao 时间: 2019-12-4 09:30

卡开发发发表于 2019-12-4 08:50
1、后续如果吸附的体系比较大的话可以适当放宽一点。
2、做收敛性测试的目的之一就是保证能量和力计算可 ...

如图，我在大师兄教程上看到的。他说的意思是只要不是变胞优化，就可以统一使用正常的相对较低的ENCUT值。所以我在想是不是可以用大的ENCUT值优化好的晶胞参数固定住，再用小的ENCUT值优化一下，后面的吸附统一用小的ENCUT值去计算，这样能量有可比性。但不知道这种想法对不对。。主要是ENCUT设置成800计算起来还是很吃力的。。。后续吸附计算都用800比较浪费时间。希望老师能够指导指导。

作者
Author: 卡开发发 时间: 2019-12-4 09:54

sunhao 发表于 2019-12-4 09:30
如图，我在大师兄教程上看到的。他说的意思是只要不是变胞优化，就可以统一使用正常的相对较低的ENCUT值 ...

无论是晶格常数计算还是一般的构型优化，stress和force对动能截断都是有依赖的，测试的目的就是为了获得可靠合理的参数，高或低应当有个定量的标准。

作者
Author: sunhao 时间: 2019-12-4 10:11

卡开发发发表于 2019-12-4 09:54
无论是晶格常数计算还是一般的构型优化，stress和force对动能截断都是有依赖的，测试的目的就是为了获得 ...

好的，我明白了。非常感谢老师的指导！

作者
Author: sunhao 时间: 2019-12-11 10:44

卡开发发发表于 2019-12-4 09:54
无论是晶格常数计算还是一般的构型优化，stress和force对动能截断都是有依赖的，测试的目的就是为了获得 ...

老师您好，我先用ISIF=3优化好了晶格，然后用ISIF=2固定住晶格之后在优化结构。用ISIF=2优化好的CONTCAR复制成POSCAR，想检验一下是不是优化好了。但是发现不是一个离子步结束，而且能量差别依然很大，磁矩也差别很大。请问老师这可能是什么原因呢？我的体系比较复杂，我设置的力的收敛标准是-0.05.

作者
Author: 卡开发发 时间: 2019-12-11 11:19

sunhao 发表于 2019-12-11 10:44
老师您好，我先用ISIF=3优化好了晶格，然后用ISIF=2固定住晶格之后在优化结构。用ISIF=2优化好的CONTCAR ...

1、一直用ISIF=3就行，没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出来，然后按照那个趋势去做初猜，看看能量是否有较大差异。

作者
Author: sunhao 时间: 2019-12-11 11:31

卡开发发发表于 2019-12-11 11:19
1、一直用ISIF=3就行，没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出 ...

好的老师，我试试。谢谢您

作者
Author: sunhao 时间: 2019-12-11 11:43

卡开发发发表于 2019-12-11 11:19
1、一直用ISIF=3就行，没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出 ...

老师您能加我一下好友吗？我想单独请教你一下。我的微信号：wxid_cjt40nkp3h9712

作者
Author: Adair 时间: 2022-1-4 17:26

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

老师，结构优化的时候是不是一般选择ALGO=Normal比较好？

作者
Author: 卡开发发 时间: 2022-1-4 23:36

Adair 发表于 2022-1-4 17:26
老师，结构优化的时候是不是一般选择ALGO=Normal比较好？

是，一般建议ALGO=N，体系很大的话，如果数值还比较稳定，ALGO=V或F也行。

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