计算化学公社

 找回密码 Forget password
 注册 Register
Views: 21490|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] VASP几何优化中电子步数的设置问题

[复制链接 Copy URL]

689

帖子

2

威望

4103

eV
积分
4832

Level 6 (一方通行)

各位好,请教二个问题
VASP程序中电子步数默认为60步,并且手册中说明这个参数一般不建议修改,在几何优化(EDIFF=1e-5,EDIFFG=-0.01)过程中经常出现SCF在60步内不收敛,VASP程序会继续优化不停止,就像高斯中的IOP(5/13=1)。既然高斯中的IOP(5/13=1)应该禁止使用,那么在VASP中是否也应该在SCF收敛的条件(加大步数,加大SIGMA值)下再优化结构?还有,在最初的几步几何优化中SCF在60步内不收敛,优化若干步后SCF就能收敛了,这样条件下得到的几何结构是否合理可信?

26

帖子

0

威望

388

eV
积分
414

Level 3 能力者

2#
发表于 Post on 2017-3-9 15:23:25 | 只看该作者 Only view this author
如果在60步不收敛,将IBRION=2改为1,cp CONTCAR POSCAR继续优化,如此反复,最终至结构收敛为止!

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
obaica + 2

查看全部评分 View all ratings

689

帖子

2

威望

4103

eV
积分
4832

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2017-3-9 15:48:30 | 只看该作者 Only view this author
mylifecg 发表于 2017-3-9 15:23
如果在60步不收敛,将IBRION=2改为1,cp CONTCAR POSCAR继续优化,如此反复,最终至结构收敛为止!

SCF收敛性不用管了吗?

3621

帖子

3

威望

1万

eV
积分
18426

Level 6 (一方通行)

第一原理惨品小作坊

4#
发表于 Post on 2017-3-9 16:36:28 | 只看该作者 Only view this author
jiangning198511 发表于 2017-3-9 15:48
SCF收敛性不用管了吗?

楼上可能误认为是离子步难收敛的情形了。

一般还是先检查结构和磁性以及和截断能、k点等参数,若不收敛可以修改ALGO或者MIX的方案,尽量还是不要将SIGMA修改过大。

若60步当中出现SCF收敛,结果很可能是不可靠的,实际再进行自洽很有可能能量与离子弛豫有差别。有时候通过cp CONTCAR再计算(注意不要读取波函数和电荷密度)能起作用,但有时候即便这样做,新的结构做出的SCF还是会出现不收敛,不太提倡这样的做法。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

689

帖子

2

威望

4103

eV
积分
4832

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2017-3-9 16:46:48 | 只看该作者 Only view this author
卡开发发 发表于 2017-3-9 16:36
楼上可能误认为是离子步难收敛的情形了。

一般还是先检查结构和磁性以及和截断能、k点等参数,若不收 ...

我碰到的是优化初始几步SCF不收敛,等优化若干步,SCF在60步能就收敛了,这样的结果是否可信

3621

帖子

3

威望

1万

eV
积分
18426

Level 6 (一方通行)

第一原理惨品小作坊

6#
发表于 Post on 2017-3-9 19:29:23 | 只看该作者 Only view this author
jiangning198511 发表于 2017-3-9 16:46
我碰到的是优化初始几步SCF不收敛,等优化若干步,SCF在60步能就收敛了,这样的结果是否可信

cp CONTCAR POSCAR之后再算一次自洽确保SCF收敛以及得到的energy、force、stress仍然落在收敛范围内。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 17:10 , Processed in 0.172726 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list