计算化学公社

标题: 求助：那种+U方式对收敛性的影响最小 [打印本页]

作者
Author: guyingbing 时间: 2021-4-9 18:52
标题: 求助：那种+U方式对收敛性的影响最小
各位老师好，我是初学vasp的同学。我注意到VASP中有三种+U方式，由关键字：LDAUTYPE= 1 | 2 | 4控制。默认的+U方式是第二种。但是我在实践过程中发现，第2种+U方式，似乎很难以收敛。能请教一下老师用那种加+U的方式对体系收敛性的影响是最小的呢？

作者
Author: 卡开发发 时间: 2021-4-9 20:19
2的话主要是各向同性的，因此只要考虑U-J就行，这样会比较方便。先保证其他计算参数、结构、磁矩等合理，如果难收敛，试着用LMAXMIX=4（d电子）或6（f电子），也许会有帮助。如果还不行就试着修改迭代方面的参数例如AMIX和AMIX_mag，特定情况下BMIX和BMIX_mag会比较有用。

作者
Author: guyingbing 时间: 2021-4-9 20:26

卡开发发发表于 2021-4-9 20:19
2的话主要是各向同性的，因此只要考虑U-J就行，这样会比较方便。先保证其他计算参数、结构、磁矩等合理，如 ...

昨天和今天我花了好长时间扫描AMIX。但是它的收敛是震荡的。花了好长的时间+磁矩不+u可以收敛，+u不加磁矩可以收敛。但是，+u+磁矩就不收敛。昨天和今天我调整了好长时间的AMIX都不收敛。在AMIX=0.05时，+u+磁矩可以在某个时候收敛，但是考察收敛曲线，发现收敛的能量是跳跃的。

作者
Author: 卡开发发 时间: 2021-4-12 09:30

guyingbing 发表于 2021-4-9 20:26
昨天和今天我花了好长时间扫描AMIX。但是它的收敛是震荡的。花了好长的时间+磁矩不+u可以收敛，+u不加磁 ...

我不知道你怎么加磁矩，原则上要开自旋极化的选项且指定原子的局域磁矩。具体还是得看你的体系，牵扯到的因素太复杂暂时没办法做进一步讨论。

作者
Author: guyingbing 时间: 2021-4-13 13:31

卡开发发发表于 2021-4-12 09:30
我不知道你怎么加磁矩，原则上要开自旋极化的选项且指定原子的局域磁矩。具体还是得看你的体系，牵扯到的 ...

感谢老师的回复。我设定了NELMDL值，现在已经可以正常收敛了。我现在生出了一个新的问题。cp2k里面有没有类似于NELMDL的功能呢？

作者
Author: 卡开发发 时间: 2021-4-13 13:36

guyingbing 发表于 2021-4-13 13:31
感谢老师的回复。我设定了NELMDL值，现在已经可以正常收敛了。我现在生出了一个 ...

没有，因为算法不同。

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