出现"ZBRENT"错误

卡开发发

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

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

lujinhui396

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

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

卡开发发

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-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 20:08
最近有点忙，这篇文献我只是简单看了一下，这篇文献的核心应该是采用了各向异性的U和J，也就是说在不同伸 ...

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

卡开发发

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