FeP211 slab cp2k scf 不收敛，麻烦指点一下谢谢

Janson_Shi · 发表于 Post on 2023-7-2 16:00:02

FeP211 slab 加了磁性，但是不收敛。此前FeP单胞cell优化成功了。麻烦指点一下，因为后面还要加Ag cluster 做分子动力学，希望有经验的高手传授下，万分感激。

sobereva · 发表于 Post on 2023-7-3 01:19:39

不要用PSOLVER ANALYTIC，PSOLVER MT好得多。XY周期性完全没必要用偶极校正

如果不是无gap或gap极小的体系，没必要smearing

盒子偏小，考虑用2 2 1 k点

注意恰当设置MAGNETIZATION

仔细看下文
解决CP2K中SCF不收敛的方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

Janson_Shi · 发表于 Post on 2023-7-3 14:23:08

sobereva 发表于 2023-7-3 01:19
不要用PSOLVER ANALYTIC，PSOLVER MT好得多。XY周期性完全没必要用偶极校正

如果不是无gap或gap极小的体 ...

非常感谢sob的建议，但是还是很难收敛，即使我按照VASP的magnetization数值来设置Fe的磁矩和自旋多重度，但是还是收敛不了。

PoorChaos · 发表于 Post on 2023-7-5 15:46:06

本帖最后由 PoorChaos 于 2023-7-5 15:50 编辑

Janson_Shi 发表于 2023-7-3 14:23
非常感谢sob的建议，但是还是很难收敛，即使我按照VASP的magnetization数值来设置Fe的磁矩和自旋多重度， ...

对于MAGETIZATION设置，我的经验
先将print level 改为medium 跑单点能
看看输出文件mulliken电荷那一段相关原子SPIN是多少，比如2.123，回到inp相应的原子MAGETIZATION改为2.5试试
泛函试试revPBE 开启OT，我做氧化铁或者Ni掺杂氧化铁，同一个结构，用revPBE 很容易收敛
强烈建议用Multiwfn生成inp

Janson_Shi · 发表于 Post on 2023-7-5 20:20:39

PoorChaos 发表于 2023-7-5 15:46
对于MAGETIZATION设置，我的经验
先将print level 改为medium 跑单点能
看看输出文件mulliken电荷那一 ...

好的谢谢经验哈

我试试先

Janson_Shi · 发表于 Post on 2023-7-5 23:26:53

PoorChaos 发表于 2023-7-5 15:46
对于MAGETIZATION设置，我的经验
先将print level 改为medium 跑单点能
看看输出文件mulliken电荷那一 ...

不好意思哈，按照你给的建议，还是不能收敛。OT 加 REVPRB 再加上 multiwfn 生成的 input 文件，scf还是收敛不了

PoorChaos · 发表于 Post on 2023-7-5 23:51:43

本帖最后由 PoorChaos 于 2023-7-5 23:57 编辑

Janson_Shi 发表于 2023-7-5 23:26
不好意思哈，按照你给的建议，还是不能收敛。OT 加 REVPRB 再加上 multiwfn 生成的 input 文件，scf还是 ...

好吧。我用multiwfn生成的inp
ALGORITHM 这里改成了IRAC

试试看

Janson_Shi · 发表于 Post on 2023-7-18 11:19:40

PoorChaos 发表于 2023-7-5 23:51
好吧。我用multiwfn生成的inp
ALGORITHM 这里改成了IRAC

牛人呀果然收敛了，谢谢帮忙哈

PoorChaos · 发表于 Post on 2023-7-18 19:16:39

Janson_Shi 发表于 2023-7-18 11:19
牛人呀果然收敛了，谢谢帮忙哈

CP2K官网/论坛有一些相关的可以看看

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] FeP211 slab cp2k scf 不收敛，麻烦指点一下谢谢