计算化学公社

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

[VASP] 如何判断VASP的几何优化是否正常收敛

[复制链接 Copy URL]

183

帖子

2

威望

4933

eV
积分
5156

Level 6 (一方通行)

能科top50

跳转到指定楼层 Go to specific reply
楼主
各位老师好,本人用vasp算了一个几何优化的任务, INCAR中NELMIN和NELM分别为5和300,意思应该是最多优化跑300步。计算结束后,发现OUTCAR文件中出现了“aborting loop because EDIFF is reached”,能证明SCF已经收敛了吗?主要是在OSZICAR中发现SCF迭代也进行了300步,难道我的作业都是300步后正好收敛吗?还是说300步以后未收敛但是自动停止了呢?

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
顾顾顾 + 3 谢谢

查看全部评分 View all ratings

In defeat, malice. In victory, revenge.

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2020-7-26 08:38:33 | 只看该作者 Only view this author
你指的是几何优化中SCF有没收敛的情况?如果指的是这个问题,可以看“aborting loop because EDIFF is reached”。

关于消除VASP自带 IOp(5/13=1)的方案,我想了一种,可以一块讨论下:
vasp源码electron.F中585行:
  1.       IF (N>=INFO%NELM) THEN
  2.          IF (.NOT.INFO%LABORT) LABORT_WITHOUT_CONV=.TRUE.
  3.          INFO%LABORT=.TRUE.
  4.       ENDIF
复制代码

改成
  1.       IF (N>=INFO%NELM) THEN
  2.          IF (.NOT.INFO%LABORT) THEN
  3.            LABORT_WITHOUT_CONV=.TRUE.
  4.            IF (IO%IU0>=0) WRITE(IO%IU0,*) "SCF Failed!"
  5.            IF (IO%IU6>=0) WRITE(IO%IU6,*) "SCF Failed!"
  6.            STOP
  7.          ENDIF
  8.          INFO%LABORT=.TRUE.
  9.       ENDIF
复制代码

评分 Rate

参与人数
Participants 4
eV +15 收起 理由
Reason
顾顾顾 + 2 赞!
thanhtam + 3 赞!
zorow + 5 谢谢
physics_xw + 5 赞!

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

183

帖子

2

威望

4933

eV
积分
5156

Level 6 (一方通行)

能科top50

3#
 楼主 Author| 发表于 Post on 2020-7-26 14:35:58 | 只看该作者 Only view this author
卡开发发 发表于 2020-7-26 08:38
你指的是几何优化中SCF有没收敛的情况?如果指的是这个问题,可以看“aborting loop because EDIFF is reac ...

谢谢老师!看到“aborting loop because EDIFF is reached”就能证明收敛了吗?主要是好几个任务跑完都是300步,虽然出现了上述文字但我还是担心这是算够300步强制停止了...
看了您说的VASP自带 IOp(5/13=1),我搜到这个帖子:http://bbs.keinsci.com/thread-11217-1-3.html,看来我要仔细检查一下OUTCAR里面是不是有电子步没收敛就继续的情况...
In defeat, malice. In victory, revenge.

183

帖子

2

威望

4933

eV
积分
5156

Level 6 (一方通行)

能科top50

4#
 楼主 Author| 发表于 Post on 2020-7-26 19:00:55 | 只看该作者 Only view this author
卡开发发 发表于 2020-7-26 08:38
你指的是几何优化中SCF有没收敛的情况?如果指的是这个问题,可以看“aborting loop because EDIFF is reac ...

老师,我检查了INCAR文件,发现我NSW写了300...好像是离子步上限到了,我是不是该加大离子步继续算... NELMIN应该是每一步最多迭代300吧,我的问题是300步之后自动停止了,然而要怎么知道NSW上限300步结束以后是否收敛了呢。。。
In defeat, malice. In victory, revenge.

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

5#
发表于 Post on 2020-7-27 08:01:26 | 只看该作者 Only view this author
zorow 发表于 2020-7-26 14:35
谢谢老师!看到“aborting loop because EDIFF is reached”就能证明收敛了吗?主要是好几个任务跑完都是 ...

按照上面修改代码后,SCF超过最大步数还没达到收敛限就直接跳错了,避免了自带 IOp(5/13=1)。只是做的还不够好,但是要做得更好也不难,比如可以在合适的位置跳错,甚至可以在keywords里面加个logical来控制SCF失败是否继续进行计算。这里只是给个大概的思路而已。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2020-7-27 08:06:05 | 只看该作者 Only view this author
zorow 发表于 2020-7-26 19:00
老师,我检查了INCAR文件,发现我NSW写了300...好像是离子步上限到了,我是不是该加大离子步继续算... NE ...

NSW=300指的是离子步为300步,即构型优化最大步数为300;NELM指的是电子步最大步数,即自洽场迭代的最大步数;NELMIN指的是至少做多少个电子步,如果是负数指的是开始的非自洽至少做多少次。我不知道你到底要问哪个没收敛?
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

183

帖子

2

威望

4933

eV
积分
5156

Level 6 (一方通行)

能科top50

7#
 楼主 Author| 发表于 Post on 2020-7-27 21:58:44 | 只看该作者 Only view this author
卡开发发 发表于 2020-7-27 08:06
NSW=300指的是离子步为300步,即构型优化最大步数为300;NELM指的是电子步最大步数,即自洽场迭代的最大 ...

是我之前没搞懂... 每个离子步中SCF的步数很快就会收敛,我觉得是离子步没有收敛,到达300步后自动停止了。所以请问,如何判断离子步收敛与否呢...
In defeat, malice. In victory, revenge.

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2020-7-28 02:48:52 | 只看该作者 Only view this author
zorow 发表于 2020-7-27 21:58
是我之前没搞懂... 每个离子步中SCF的步数很快就会收敛,我觉得是离子步没有收敛,到达300步后自动停止了 ...

"reached required accuracy - stopping structural energy minimisation"
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

183

帖子

2

威望

4933

eV
积分
5156

Level 6 (一方通行)

能科top50

9#
 楼主 Author| 发表于 Post on 2020-7-28 12:08:14 | 只看该作者 Only view this author
卡开发发 发表于 2020-7-28 02:48
"reached required accuracy - stopping structural energy minimisation"

明白了,谢谢老师指教!
In defeat, malice. In victory, revenge.

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 05:58 , Processed in 0.168847 second(s), 24 queries , Gzip On.

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