出现"ZBRENT"错误

lujinhui396

采用vasp优化时，IBRION=2，出现" ZBRENT: increasing intervall"
" ZBRENT: bracketing found"
" ZBRENT: interpolating"
" ZBRENT: can't locate minimum, use default step"
在这种情况下，收敛得到的结果是否可信？如何避免这种情况？

卡开发发

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这个时候应该能够较快接近稳定结构，太大的话就可能造成振荡。

lujinhui396

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

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

卡开发发

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.

卡开发发

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是按照阻尼的结果算的。

lujinhui396

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

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

卡开发发

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

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

lujinhui396

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

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

lujinhui396

@卡开发发  

卡开发发

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

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

lujinhui396

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

卡开发发

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

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

lujinhui396

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

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

冰山

貌似EDIFFG=-0.02有点小，还有KPOINT也是要测试一下。
建议加强收敛标准，再优化看看结构有变化没？

卡开发发

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

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