计算化学公社

标题: 请教ALGO的问题 [打印本页]

作者
Author:
Theor_Comp    时间: 2019-9-3 14:48
标题: 请教ALGO的问题
请教大家:使用vasp在结构优化时,电子步不收敛,将ALGO由FAST改为ALL后计算了单点并保存了WAVECAR,之后读取WAVECAR进行结构优化,还需要把ALGO改为FAST吗?还是说ALGO=ALL也可以???


作者
Author:
卡开发发    时间: 2019-9-3 17:12
可以改F或者N,已经收敛的情况A效率就有点低了。构型优化过后最好再做下SCF。
作者
Author:
Theor_Comp    时间: 2019-9-3 17:55
卡开发发 发表于 2019-9-3 17:12
可以改F或者N,已经收敛的情况A效率就有点低了。构型优化过后最好再做下SCF。

使用ALGO=A计算单点保存WAVECAR,之后读取WAVECAR进行结构优化,使用ALGO=A,结构优化也很快能够完成,这样的话,还有必要把ALGO改为FAST吗? (GGA=PE泛函下)
作者
Author:
卡开发发    时间: 2019-9-3 18:00
Theor_Comp 发表于 2019-9-3 17:55
使用ALGO=A计算单点保存WAVECAR,之后读取WAVECAR进行结构优化,使用ALGO=A,结构优化也很快能够完成,这 ...

没。但还是建议做个SCF看下能量差别多大,以及自洽是否能正常收敛。
作者
Author:
Theor_Comp    时间: 2019-9-3 18:16
卡开发发 发表于 2019-9-3 18:00
没。但还是建议做个SCF看下能量差别多大,以及自洽是否能正常收敛。

谢谢您。ALGO=A相比ALGO=FAST只是会减慢计算过程吗?
另外,您的意思是:在ALGO=A优化的结构上,再分别使用ALGO=A和ALGO=FAST来进行静态自恰??
作者
Author:
卡开发发    时间: 2019-9-3 18:41
Theor_Comp 发表于 2019-9-3 18:16
谢谢您。ALGO=A相比ALGO=FAST只是会减慢计算过程吗?
另外,您的意思是:在ALGO=A优化的结构上,再分别 ...

恩,只是不同的算法,不过要注意:
The preconditioned conjugate gradient (IALGO = 58, ALGO = A) algorithm is recommended for insulators. The best stability is usually obtained if the number of bands equals half the number of electrons (non spin polarized case).

不用分别用ALGO=A和F做SCF,选一个就行。
作者
Author:
Theor_Comp    时间: 2019-9-3 20:14
卡开发发 发表于 2019-9-3 18:41
恩,只是不同的算法,不过要注意:

不用分别用ALGO=A和F做SCF,选一个就行。

非常感谢您!
您前面提到“建议做个SCF看下能量差别多大”,您的意思是:使用ALGO=A的优化的最后结构,来进行静态自恰计算,这时会调整体系电子的运动,看看SCF能量与上一步优化结构时的能量差别多少吗?
作者
Author:
卡开发发    时间: 2019-9-3 20:47
Theor_Comp 发表于 2019-9-3 20:14
非常感谢您!
您前面提到“建议做个SCF看下能量差别多大”,您的意思是:使用ALGO=A的优化的最后结构, ...

对。有时候可能有点差别,差别比较大最好cp CONTCAR POSCAR继续再做。
作者
Author:
Theor_Comp    时间: 2019-9-3 21:08
卡开发发 发表于 2019-9-3 20:47
对。有时候可能有点差别,差别比较大最好cp CONTCAR POSCAR继续再做。

1.ALGO = A) algorithm is recommended for insulators. 那ALGO=A能用在半导体计算中吗?
2.结构优化、静态自恰、非自恰计算中,需要保持ALGO都使用相同的吗?
作者
Author:
卡开发发    时间: 2019-9-3 21:14
Theor_Comp 发表于 2019-9-3 21:08
1.ALGO = A) algorithm is recommended for insulators. 那ALGO=A能用在半导体计算中吗?
2.结构优化、 ...

1、自旋极化的情况最好用D。
2、一般不需要,除非是做特定的计算。
作者
Author:
Theor_Comp    时间: 2019-9-3 21:30
本帖最后由 Theor_Comp 于 2019-9-3 21:40 编辑
卡开发发 发表于 2019-9-3 21:14
1、自旋极化的情况最好用D。
2、一般不需要,除非是做特定的计算。

您好,如果电子步收敛,而离子步一直不收敛怎么办呢?采用的是EDIFFG=-0.015。看OSZ文件时,能量基本不变了,但就是没有正常停止,一直在算。
作者
Author:
卡开发发    时间: 2019-9-3 22:15
Theor_Comp 发表于 2019-9-3 21:30
您好,如果电子步收敛,而离子步一直不收敛怎么办呢?采用的是EDIFFG=-0.015。看OSZ文件时,能量基本不变 ...

那有可能是力没收敛,EDIFFG<0用的是力收敛标准。具体原因这个我不好说,也许是SCF收敛精度太低,也可能是收敛性测试没做好,或者是结构不太合理等等,看看从这几个方面进行排查。
作者
Author:
Theor_Comp    时间: 2019-9-4 09:55
卡开发发 发表于 2019-9-3 22:15
那有可能是力没收敛,EDIFFG

好的,非常感谢。
另外,还有一个疑问,在一个任务中,读取上一步产生的WAVECAR(上一步的INCAR中LWAVE=TRUE)时,该任务INCAR中ISATRT=2,那LWAVE是为TRUE还是FALSE呢?如果为true,会不会覆盖掉前面的WAVECAR??
作者
Author:
柠檬鬼lemon    时间: 2023-11-15 21:14
请问我一将ALGO=ALL或者Damped就停在了第六步不算了,如果用FAST或者Normal就算的很快是怎么回事呀,其他的参数我都没有变,这个是服务器的问题吗,我已经试过很多次了,就是只要一用ALL就卡住如下图这里,就一直停在这里了




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3