请教ALGO的问题

Theor_Comp · 发表于 Post on 2019-9-3 14:48:25

请教大家：使用vasp在结构优化时，电子步不收敛，将ALGO由FAST改为ALL后计算了单点并保存了WAVECAR，之后读取WAVECAR进行结构优化，还需要把ALGO改为FAST吗？还是说ALGO=ALL也可以？？？

卡开发发 · 发表于 Post on 2019-9-3 17:12:26

可以改F或者N，已经收敛的情况A效率就有点低了。构型优化过后最好再做下SCF。

Theor_Comp · 发表于 Post on 2019-9-3 17:55:05

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

使用ALGO=A计算单点保存WAVECAR，之后读取WAVECAR进行结构优化，使用ALGO=A，结构优化也很快能够完成，这样的话，还有必要把ALGO改为FAST吗？（GGA=PE泛函下）

卡开发发 · 发表于 Post on 2019-9-3 18:00:32

Theor_Comp 发表于 2019-9-3 17:55
使用ALGO=A计算单点保存WAVECAR，之后读取WAVECAR进行结构优化，使用ALGO=A，结构优化也很快能够完成，这 ...

没。但还是建议做个SCF看下能量差别多大，以及自洽是否能正常收敛。

Theor_Comp · 发表于 Post on 2019-9-3 18:16:03

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

谢谢您。ALGO=A相比ALGO=FAST只是会减慢计算过程吗？
另外，您的意思是：在ALGO=A优化的结构上，再分别使用ALGO=A和ALGO=FAST来进行静态自恰??

卡开发发 · 发表于 Post on 2019-9-3 18:41:23

Theor_Comp 发表于 2019-9-3 18:16
谢谢您。ALGO=A相比ALGO=FAST只是会减慢计算过程吗？
另外，您的意思是：在ALGO=A优化的结构上，再分别 ...

恩，只是不同的算法，不过要注意：

The preconditioned conjugate gradient (IALGO = 58ALGO = 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，选一个就行。

Theor_Comp · 发表于 Post on 2019-9-3 20:14:57

卡开发发发表于 2019-9-3 18:41
恩，只是不同的算法，不过要注意：

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

非常感谢您！
您前面提到“建议做个SCF看下能量差别多大”，您的意思是：使用ALGO=A的优化的最后结构，来进行静态自恰计算，这时会调整体系电子的运动，看看SCF能量与上一步优化结构时的能量差别多少吗？

卡开发发 · 发表于 Post on 2019-9-3 20:47:00

Theor_Comp 发表于 2019-9-3 20:14
非常感谢您！
您前面提到“建议做个SCF看下能量差别多大”，您的意思是：使用ALGO=A的优化的最后结构， ...

对。有时候可能有点差别，差别比较大最好cp CONTCAR POSCAR继续再做。

Theor_Comp · 发表于 Post on 2019-9-3 21:08:00

卡开发发发表于 2019-9-3 20:47
对。有时候可能有点差别，差别比较大最好cp CONTCAR POSCAR继续再做。

1.ALGO = A) algorithm is recommended for insulators. 那ALGO=A能用在半导体计算中吗？
2.结构优化、静态自恰、非自恰计算中，需要保持ALGO都使用相同的吗？

卡开发发 · 发表于 Post on 2019-9-3 21:14:01

Theor_Comp 发表于 2019-9-3 21:08
1.ALGO = A) algorithm is recommended for insulators. 那ALGO=A能用在半导体计算中吗？
2.结构优化、 ...

1、自旋极化的情况最好用D。
2、一般不需要，除非是做特定的计算。

Theor_Comp · 发表于 Post on 2019-9-3 21:30:53

本帖最后由 Theor_Comp 于 2019-9-3 21:40 编辑

卡开发发发表于 2019-9-3 21:14
1、自旋极化的情况最好用D。
2、一般不需要，除非是做特定的计算。

您好，如果电子步收敛，而离子步一直不收敛怎么办呢？采用的是EDIFFG=-0.015。看OSZ文件时，能量基本不变了，但就是没有正常停止，一直在算。

卡开发发 · 发表于 Post on 2019-9-3 22:15:30

Theor_Comp 发表于 2019-9-3 21:30
您好，如果电子步收敛，而离子步一直不收敛怎么办呢？采用的是EDIFFG=-0.015。看OSZ文件时，能量基本不变 ...

那有可能是力没收敛，EDIFFG<0用的是力收敛标准。具体原因这个我不好说，也许是SCF收敛精度太低，也可能是收敛性测试没做好，或者是结构不太合理等等，看看从这几个方面进行排查。

Theor_Comp · 发表于 Post on 2019-9-4 09:55:54

卡开发发发表于 2019-9-3 22:15
那有可能是力没收敛，EDIFFG

好的，非常感谢。
另外，还有一个疑问，在一个任务中，读取上一步产生的WAVECAR（上一步的INCAR中LWAVE=TRUE）时，该任务INCAR中ISATRT=2,那LWAVE是为TRUE还是FALSE呢？如果为true，会不会覆盖掉前面的WAVECAR？？

柠檬鬼lemon · 发表于 Post on 2023-11-15 21:14:11

请问我一将ALGO=ALL或者Damped就停在了第六步不算了，如果用FAST或者Normal就算的很快是怎么回事呀，其他的参数我都没有变，这个是服务器的问题吗，我已经试过很多次了，就是只要一用ALL就卡住如下图这里，就一直停在这里了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 请教ALGO的问题

评分 Rate