计算化学公社

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

[VASP] 求助:那种+U方式对收敛性的影响最小

[复制链接 Copy URL]

47

帖子

0

威望

552

eV
积分
599

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我是初学vasp的同学。我注意到VASP中有三种+U方式,由关键字:LDAUTYPE= 1 | 2 | 4控制。默认的+U方式是第二种。但是我在实践过程中发现,第2种+U方式,似乎很难以收敛。能请教一下老师用那种加+U的方式对体系收敛性的影响是最小的呢?

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2021-4-13 13:36:50 | 只看该作者 Only view this author
guyingbing 发表于 2021-4-13 13:31
感谢老师的回复。我设定了NELMDL值,现在已经可以正常收敛了。我现在生出了一个 ...

没有,因为算法不同。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

47

帖子

0

威望

552

eV
积分
599

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-4-13 13:31:03 | 只看该作者 Only view this author
卡开发发 发表于 2021-4-12 09:30
我不知道你怎么加磁矩,原则上要开自旋极化的选项且指定原子的局域磁矩。具体还是得看你的体系,牵扯到的 ...

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

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2021-4-12 09:30:13 | 只看该作者 Only view this author
guyingbing 发表于 2021-4-9 20:26
昨天和今天我花了好长时间扫描AMIX。但是它的收敛是震荡的。花了好长的时间+磁矩不+u可以收敛,+u不加磁 ...

我不知道你怎么加磁矩,原则上要开自旋极化的选项且指定原子的局域磁矩。具体还是得看你的体系,牵扯到的因素太复杂暂时没办法做进一步讨论。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

47

帖子

0

威望

552

eV
积分
599

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-4-9 20:26:21 | 只看该作者 Only view this author
卡开发发 发表于 2021-4-9 20:19
2的话主要是各向同性的,因此只要考虑U-J就行,这样会比较方便。先保证其他计算参数、结构、磁矩等合理,如 ...

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

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2021-4-9 20:19:12 | 只看该作者 Only view this author
2的话主要是各向同性的,因此只要考虑U-J就行,这样会比较方便。先保证其他计算参数、结构、磁矩等合理,如果难收敛,试着用LMAXMIX=4(d电子)或6(f电子),也许会有帮助。如果还不行就试着修改迭代方面的参数例如AMIX和AMIX_mag,特定情况下BMIX和BMIX_mag会比较有用。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

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

GMT+8, 2026-2-18 05:26 , Processed in 0.159723 second(s), 21 queries , Gzip On.

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