求助:结构优化时采用ALGO=F的优化算法导致任务失败（总能为正值）的原因

BADEGG · 发表于 Post on 2022-11-17 09:15:28

本帖最后由 BADEGG 于 2022-11-17 21:11 编辑

各位老师好，我在做一个O2吸附的结构优化时发现，采用ALGO=F的电子步优化算法时，能量会在DAV算法转换成RMM算法时由负值变为正值

202211170902059839..png (44.16 KB下载次数 Times of downloads: 23)

下载附件 Download

采用ALGO=F的output

2022-11-17 09:02 上传 Uploaded

而采用ALGO=N则不会出现这样的情况

202211170900349377..png (42 KB下载次数 Times of downloads: 29)

下载附件 Download

采用ALGO=N的output

2022-11-17 09:00 上传 Uploaded

我上传了我的POSCAR,KPOINTS和INCAR文件，泛函采用的是PBE。

INCAR (511 Bytes下载次数 Times of downloads: 0)

2022-11-17 08:58 上传 Uploaded

点击下载
Click to download

KPOINTS (25 Bytes下载次数 Times of downloads: 0)

2022-11-17 08:58 上传 Uploaded

点击下载
Click to download

POSCAR (11.52 KB下载次数 Times of downloads: 0)

2022-11-17 08:58 上传 Uploaded

点击下载
Click to download

问题一：为什么会出现这样的情况，我发现这样的错误似乎时随机的，因为我在同样的计算体系中，换一个吸附物种，比如CO2，我就能顺利的使用ALGO=F算法得到稳定的结构；
问题二：在计算中切换使用ALGO=F和ALGO=N算法对我计算结果的准确性是否有影响？
问题三：我通常将我的电子步优化步数确定为100，但是部分优化过程中出现至多6个离子步中电子步的不收敛（最终我依然能得到稳定的结构），那么在这6个离子步中没能找到稳定势能面是否对我最终得到的稳定结构产生影响？

jiangning198511 · 发表于 Post on 2022-11-17 15:19:18

对于问题二而言，不同的优化算法得到的结果应该是一样的或者差距很微小

BADEGG · 发表于 Post on 2022-11-17 21:08:11

jiangning198511 发表于 2022-11-17 15:19
对于问题二而言，不同的优化算法得到的结果应该是一样的或者差距很微小

好的，谢谢老师的解答

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

[VASP] 求助:结构优化时采用ALGO=F的优化算法导致任务失败（总能为正值）的原因