计算化学公社

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

[VASP] 出现"ZBRENT"错误

[复制链接 Copy URL]

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
采用vasp优化时,IBRION=2,出现" ZBRENT: increasing intervall"
" ZBRENT: bracketing found"
" ZBRENT: interpolating"
" ZBRENT: can't locate minimum, use default step"

在这种情况下,收敛得到的结果是否可信?如何避免这种情况?


3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2014-11-30 13:22:53 | 只看该作者 Only view this author
(1)你可以看看是否达到计算的预设精度,如果达到了原则上应该结果是可信的。
(2)可以试试把CONTCAR拷贝成POSCAR继续优化,也可以尝试减小POTIM,或者试试修改IBRION即换一种离子步的代数。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2014-12-2 15:38:54 | 只看该作者 Only view this author
卡开发发 发表于 2014-11-30 13:22
(1)你可以看看是否达到计算的预设精度,如果达到了原则上应该结果是可信的。
(2)可以试试把CONTCAR拷 ...

我采用了CP CONTCAR POSCAR 继续优化,得到收敛的结构a,但是中间的离子步依然出现ZBRENT的警告;换了一种方法,即采用IBRION=1,也收敛了,得到结构b.
结构a:键角为144.4 能量为-323.9
结构b:键角为147.8,能量为-323.68
请问这两个该如何取舍呢?

3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

4#
发表于 Post on 2014-12-2 17:01:23 | 只看该作者 Only view this author
lujinhui396 发表于 2014-12-2 15:38
我采用了CP CONTCAR POSCAR 继续优化,得到收敛的结构a,但是中间的离子步依然出现ZBRENT的警告;换了一种 ...

有警告的应该不怎么可靠,后面的应该可以吧,不过我不确定电子步的设置是否合理。也看到有在INCAR增加ADDGRID=.True的情况。上面写错了,应该是适当增加POTIM,不过这种方法有时候可能不奏效,我自己以前试过效果并不好。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2014-12-2 19:58:11 | 只看该作者 Only view this author
卡开发发 发表于 2014-12-2 17:01
有警告的应该不怎么可靠,后面的应该可以吧,不过我不确定电子步的设置是否合理。也看到有在INCAR增加ADD ...

增加ADDGRID=.TRUE.的方法我试过了,也不行。电子步的设置,指的是什么?收敛判据为EDIFF=1E-5,NELMIN=5,NELM为默认。

3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

6#
发表于 Post on 2014-12-2 22:03:26 | 只看该作者 Only view this author
lujinhui396 发表于 2014-12-2 19:58
增加ADDGRID=.TRUE.的方法我试过了,也不行。电子步的设置,指的是什么?收敛判据为EDIFF=1E-5,NELMIN=5 ...

我指的是截断能和以及精度相关的,包括KPOINTS。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2014-12-3 09:23:34 | 只看该作者 Only view this author
卡开发发 发表于 2014-12-2 22:03
我指的是截断能和以及精度相关的,包括KPOINTS。

截断能做过测试,没问题。就是KPOINTS,参考文献。精度是EDIFFG=-0.02
现在在优化好的结构上进行自洽计算,不收敛。精度没变。

3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

8#
发表于 Post on 2014-12-3 14:58:48 | 只看该作者 Only view this author
lujinhui396 发表于 2014-12-3 09:23
截断能做过测试,没问题。就是KPOINTS,参考文献。精度是EDIFFG=-0.02
现在在优化好的结构上进行自洽计算 ...

电子步收敛都这么难,还是先看看结构本身是否有问题。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

29

帖子

0

威望

132

eV
积分
161

Level 3 能力者

9#
发表于 Post on 2014-12-4 09:30:55 | 只看该作者 Only view this author
貌似EDIFFG=-0.02有点小,还有KPOINT也是要测试一下。
建议加强收敛标准,再优化看看结构有变化没?

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2014-12-10 16:34:52 | 只看该作者 Only view this author
卡开发发 发表于 2014-11-30 13:22
(1)你可以看看是否达到计算的预设精度,如果达到了原则上应该结果是可信的。
(2)可以试试把CONTCAR拷 ...

可以试试把CONTCAR拷贝成POSCAR继续优化,
如果还出现这个警告,是不是要继续cp CONTCAR POSCAR

3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

11#
发表于 Post on 2014-12-10 16:49:58 | 只看该作者 Only view this author
lujinhui396 发表于 2014-12-10 16:34
可以试试把CONTCAR拷贝成POSCAR继续优化,
如果还出现这个警告,是不是要继续cp CONTCAR POSCAR

如果已经做过这个我觉得可能用处不大,试试其他方法吧。我总觉得U的设置好像有些问题,这几天抽空我做个测试。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2014-12-10 17:15:57 | 只看该作者 Only view this author
卡开发发 发表于 2014-12-10 16:49
如果已经做过这个我觉得可能用处不大,试试其他方法吧。我总觉得U的设置好像有些问题,这几天抽空我做个 ...

U的设置那里出现问题了呢?这是我新采用的U的设置:LDAU= .TRUE.
LDAUTYPE= 2
LDAUL=2 -1
LDAUU=4.5 0.0
LDAUJ=1.0 0.0
LDAUPRINT=2

3622

帖子

3

威望

1万

eV
积分
18442

Level 6 (一方通行)

第一原理惨品小作坊

13#
发表于 Post on 2014-12-10 18:39:50 | 只看该作者 Only view this author
lujinhui396 发表于 2014-12-10 17:15
U的设置那里出现问题了呢?这是我新采用的U的设置:LDAU= .TRUE.
LDAUTYPE= 2
LDAUL=2 -1

我的意思是 按说LDAUTYPE=2只取决于Ueff这一项,我先找找原因。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2014-12-10 20:24:22 | 只看该作者 Only view this author
@卡开发发  

PhysRevB.86.2051261.pdf

640.61 KB, 下载次数 Times of downloads: 145

58

帖子

0

威望

88

eV
积分
146

Level 2 能力者

15#
 楼主 Author| 发表于 Post on 2014-12-11 09:58:29 | 只看该作者 Only view this author
卡开发发 发表于 2014-12-10 18:39
我的意思是 按说LDAUTYPE=2只取决于Ueff这一项,我先找找原因。

我在下面发了个U的文献,我有点看不懂,你看看

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 06:00 , Processed in 0.175538 second(s), 24 queries , Gzip On.

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