计算化学公社

标题: FeP211 slab cp2k scf 不收敛，麻烦指点一下谢谢 [打印本页]

作者
Author: Janson_Shi 时间: 2023-7-2 16:00
标题: FeP211 slab cp2k scf 不收敛，麻烦指点一下谢谢
FeP211 slab 加了磁性，但是不收敛。此前FeP单胞cell优化成功了。麻烦指点一下，因为后面还要加Ag cluster 做分子动力学，希望有经验的高手传授下，万分感激。

作者
Author: sobereva 时间: 2023-7-3 01:19
不要用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）

作者
Author: Janson_Shi 时间: 2023-7-3 14:23

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

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

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

作者
Author: PoorChaos 时间: 2023-7-5 15:46
本帖最后由 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

作者
Author: Janson_Shi 时间: 2023-7-5 20:20

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

好的谢谢经验哈

我试试先

作者
Author: Janson_Shi 时间: 2023-7-5 23:26

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

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

作者
Author: PoorChaos 时间: 2023-7-5 23:51
本帖最后由 PoorChaos 于 2023-7-5 23:57 编辑

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

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

试试看

作者
Author: Janson_Shi 时间: 2023-7-18 11:19

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

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

作者
Author: PoorChaos 时间: 2023-7-18 19:16

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

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

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