计算化学公社

标题: 出现"ZBRENT"错误 [打印本页]

作者
Author: lujinhui396 时间: 2014-11-30 10:33
标题: 出现"ZBRENT"错误
采用vasp优化时，IBRION=2，出现" ZBRENT: increasing intervall"
" ZBRENT: bracketing found"
" ZBRENT: interpolating"
" ZBRENT: can't locate minimum, use default step"
在这种情况下，收敛得到的结果是否可信？如何避免这种情况？

作者
Author: 卡开发发 时间: 2014-11-30 13:22
（1）你可以看看是否达到计算的预设精度，如果达到了原则上应该结果是可信的。
（2）可以试试把CONTCAR拷贝成POSCAR继续优化，也可以尝试减小POTIM，或者试试修改IBRION即换一种离子步的代数。

作者
Author: lujinhui396 时间: 2014-12-2 15:38

卡开发发发表于 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
请问这两个该如何取舍呢？

作者
Author: 卡开发发 时间: 2014-12-2 17:01

lujinhui396 发表于 2014-12-2 15:38
我采用了CP CONTCAR POSCAR 继续优化，得到收敛的结构a,但是中间的离子步依然出现ZBRENT的警告;换了一种 ...

有警告的应该不怎么可靠，后面的应该可以吧，不过我不确定电子步的设置是否合理。也看到有在INCAR增加ADDGRID=.True的情况。上面写错了，应该是适当增加POTIM，不过这种方法有时候可能不奏效，我自己以前试过效果并不好。

作者
Author: lujinhui396 时间: 2014-12-2 19:58

卡开发发发表于 2014-12-2 17:01
有警告的应该不怎么可靠，后面的应该可以吧，不过我不确定电子步的设置是否合理。也看到有在INCAR增加ADD ...

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

作者
Author: 卡开发发 时间: 2014-12-2 22:03

lujinhui396 发表于 2014-12-2 19:58
增加ADDGRID=.TRUE.的方法我试过了，也不行。电子步的设置，指的是什么？收敛判据为EDIFF=1E-5，NELMIN=5 ...

我指的是截断能和以及精度相关的，包括KPOINTS。

作者
Author: lujinhui396 时间: 2014-12-3 09:23

卡开发发发表于 2014-12-2 22:03
我指的是截断能和以及精度相关的，包括KPOINTS。

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

作者
Author: 卡开发发 时间: 2014-12-3 14:58

lujinhui396 发表于 2014-12-3 09:23
截断能做过测试，没问题。就是KPOINTS,参考文献。精度是EDIFFG=-0.02
现在在优化好的结构上进行自洽计算 ...

电子步收敛都这么难，还是先看看结构本身是否有问题。

作者
Author: 冰山 时间: 2014-12-4 09:30
貌似EDIFFG=-0.02有点小，还有KPOINT也是要测试一下。
建议加强收敛标准，再优化看看结构有变化没？

作者
Author: lujinhui396 时间: 2014-12-10 16:34

卡开发发发表于 2014-11-30 13:22
（1）你可以看看是否达到计算的预设精度，如果达到了原则上应该结果是可信的。
（2）可以试试把CONTCAR拷 ...

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

作者
Author: 卡开发发 时间: 2014-12-10 16:49

lujinhui396 发表于 2014-12-10 16:34
可以试试把CONTCAR拷贝成POSCAR继续优化,
如果还出现这个警告，是不是要继续cp CONTCAR POSCAR

如果已经做过这个我觉得可能用处不大，试试其他方法吧。我总觉得U的设置好像有些问题，这几天抽空我做个测试。

作者
Author: lujinhui396 时间: 2014-12-10 17:15

卡开发发发表于 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

作者
Author: 卡开发发 时间: 2014-12-10 18:39

lujinhui396 发表于 2014-12-10 17:15
U的设置那里出现问题了呢？这是我新采用的U的设置：LDAU= .TRUE.
LDAUTYPE= 2
LDAUL=2 -1

我的意思是按说LDAUTYPE=2只取决于Ueff这一项，我先找找原因。

作者
Author: lujinhui396 时间: 2014-12-10 20:24
@卡开发发

作者
Author: lujinhui396 时间: 2014-12-11 09:58

卡开发发发表于 2014-12-10 18:39
我的意思是按说LDAUTYPE=2只取决于Ueff这一项，我先找找原因。

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

作者
Author: 卡开发发 时间: 2014-12-11 13:17

lujinhui396 发表于 2014-12-11 09:58
我在下面发了个U的文献，我有点看不懂，你看看

Ok，我来看一下，然后把你那个测试跑一下。

作者
Author: lujinhui396 时间: 2014-12-11 14:58

卡开发发发表于 2014-12-2 17:01
有警告的应该不怎么可靠，后面的应该可以吧，不过我不确定电子步的设置是否合理。也看到有在INCAR增加ADD ...

请问这个POTIM对于CG算法，是不是只决定第一步的步长，后面的步长是通过trialstep决定的？

作者
Author: 卡开发发 时间: 2014-12-11 19:14
本帖最后由卡开发发于 2014-12-11 19:28 编辑

lujinhui396 发表于 2014-12-11 14:58
请问这个POTIM对于CG算法，是不是只决定第一步的步长，后面的步长是通过trialstep决定的？

不是，你看的可能是准牛顿法的，CG的方法应该是第一步是寻找下降方向，然后后面的Hessian都是*POTIM进行阻尼，手册原话：The conjugate gradient method requires a line minimization, which is performed in several steps: i) first a trial step into the search direction (scaled gradients) is done, with
the length of the trial step controlled by the POTIM parameter (section 6.23).
6.23的公式显然STEP是按照阻尼的结果算的。

作者
Author: 卡开发发 时间: 2014-12-11 20:08

lujinhui396 发表于 2014-12-11 09:58
我在下面发了个U的文献，我有点看不懂，你看看

最近有点忙，这篇文献我只是简单看了一下，这篇文献的核心应该是采用了各向异性的U和J，也就是说在不同伸展方向(磁量子数不同)的相互作用的两个d或f的U和J采用的值不相同，从而U和J就变成了张量。话说回来，对于VASP的Type=2的LDAU方法，还是只有U-J有意义。
Note: in Dudarev’s approach the parameters U and J do not enter seperately, only the difference (U− J) is meaningfull.

作者
Author: lujinhui396 时间: 2014-12-11 20:31

卡开发发发表于 2014-12-11 20:08
最近有点忙，这篇文献我只是简单看了一下，这篇文献的核心应该是采用了各向异性的U和J，也就是说在不同伸 ...

1 那他采用的LDAUTYPE应该是多少呢？为什么他保持U-J一样，仅改变了U和J的数值，能量就发生变化了呢？
2 还有一个问题，就是POTIM的设置。当初始结构严重偏离平衡结构时，是应该增大还是减小POTIM呢？为什么？

作者
Author: 卡开发发 时间: 2014-12-11 20:55

lujinhui396 发表于 2014-12-11 20:31
1 那他采用的LDAUTYPE应该是多少呢？为什么他保持U-J一样，仅改变了U和J的数值，能量就发生变化了呢？
2 ...

1、这种方法我在手册中没找到，不知是不是自己编程来实现的；虽然文献U-J一样，但这里不存在Ueff的概念，这里的U和J指的都是U和J张量的模值之差是恒定的。这种不同于我们采用的Type=2的形式，采用Type=2的时候按手册的说法只有Ueff有意义；
2、初始结构严重偏离的时候建议还是采用稍大些的POTIM这个时候应该能够较快接近稳定结构，太大的话就可能造成振荡。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)