计算化学公社

 找回密码 Forget password
 注册 Register

CP2K计算Gd-Ru-O体系SCF收敛难的问题

查看数: 317 | 评论数: 5 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-9-29 13:27

正文摘要:

各位老师好,我目前采用2024.1版本CP2K计算Gd16Ru48O128的(100)/(110)/(111)表面上的OER反应,计算级别为PBE+DZVP-MOLOPT-SR-GTH,采用Multiwfn得到CP2K输入文件,从Materials Project(MP)查到的结果是RuO2是无磁 ...

回复 Reply

Hedonism 发表于 Post on 2025-10-4 14:25:08
yxdd98 发表于 2025-10-2 22:19
使用vasp的话能收敛了吗?

可以收敛
yxdd98 发表于 Post on 2025-10-2 22:19:45
Hedonism 发表于 2025-10-2 20:35
后记:

该体系用CP2K在对角化条件下应该是收敛不了了。

使用vasp的话能收敛了吗?
Hedonism 发表于 Post on 2025-10-2 20:35:45
后记:

该体系用CP2K在对角化条件下应该是收敛不了了。

不自定义组态时,50步的SCF初猜过程中Gd只能收敛到E-2级别。

用&BS手动指定了组态后,20步以内即可到E-6级别,且SCF计算的第一步能量非常接近OT计算得到的能量,但在第二步开始就震荡了,本人无法解决,目前已改用VASP
Hedonism 发表于 Post on 2025-9-29 16:39:01
丁越 发表于 2025-9-29 16:20
1. 参考http://bbs.keinsci.com/thread-37196-1-1.html

2. 这种情况下只能通过对比设置磁矩(不同磁矩大 ...

感谢老师的回复,我继续按您回复的1和2尝试一下

我进一步准备用VASP先计算得到一个磁矩,后用&BS手动定义组态试一试。
丁越 发表于 Post on 2025-9-29 16:20:21
1. 参考http://bbs.keinsci.com/thread-37196-1-1.html

2. 这种情况下只能通过对比设置磁矩(不同磁矩大小)以及不设置磁矩情况下体系能量哪个低就将其作为基态,但是前提是你得把SCF弄收敛了。

3. 对于你开了smearing情况下,体系的电子是按照自旋轨道能量由低到高占据的,不会区分alpha或者beta轨道,这里自旋多重度的设置因此不会对电子占据造成影响,你设置一个最低的自旋多重度即可。
我建议要是多次尝试仍不收敛的话直接换程序,CP2K在算一些过渡金属氧化物体系就是存在难收敛的问题,况且你的还是挺复杂的氧化物。

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

GMT+8, 2026-1-25 13:51 , Processed in 0.179652 second(s), 25 queries , Gzip On.

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