计算化学公社

标题: Fe(110)晶面吸附计算电子步难收敛问题如何解决 [打印本页]

作者
Author: 玄武一甩 时间: 2022-1-6 11:20
标题: Fe(110)晶面吸附计算电子步难收敛问题如何解决
最近在做Fe(110)晶面的结构优化、吸附计算工作，发现Fe（110）电子步很难收敛，通常200步以内无法收敛，目前的晶面结构是4层，4*4，共256个原子，INCAR文件如下如所示，想问下大家是否有参数需要添加或修改，来加快电子步收敛 (, 下载次数 Times of downloads: 14)

作者
Author: Aristotler 时间: 2022-1-6 11:53
SIGMA=0.1 或者等于0.2试一下或者ALGO=48；你的问题在于Fe你用了自旋，并且截断能太大了（400或者450足以）

作者
Author: Aristotler 时间: 2022-1-6 11:55
本帖最后由 Aristotler 于 2022-1-6 11:57 编辑

KPOINTS 也可以改小点，模型原子数太多，一般吸附质吸附后，扩胞，相邻晶胞中吸附质距离大于10埃就行。没必要把基底ab面设置过大

作者
Author: 玄武一甩 时间: 2022-1-6 12:16

Aristotler 发表于 2022-1-6 11:53
SIGMA=0.1 或者等于0.2试一下或者ALGO=48；你的问题在于Fe你用了自旋，并且截断能太大了（400或者450足以 ...

哇，谢谢，我来尝试下，因为我的吸附质分子比较大，所以扩胞就大了点，原子数有点多

作者
Author: Para1lel 时间: 2022-1-6 16:14
Fe有磁性，设个初始磁矩比较好，MAGMOM这个参数设一下。

作者
Author: wsz 时间: 2022-1-7 02:38
建议首先尝试以下方法，如果还不能收敛，可以参考此贴http://bbs.keinsci.com/thread-13065-1-1.html的其他方法。
(1)设置MAGMOM、AMIX;
(2)增加积分格点PREC=Accurate;
(3)使用ALGO=Normal;
(4)金属体系可以使用ISMEAR=1 or 2。

作者
Author: 玄武一甩 时间: 2022-1-7 20:33

Para1lel 发表于 2022-1-6 16:14
Fe有磁性，设个初始磁矩比较好，MAGMOM这个参数设一下。

请问这个参数设置多少比较好

作者
Author: Para1lel 时间: 2022-1-8 08:42

玄武一甩发表于 2022-1-7 20:33
请问这个参数设置多少比较好

未成对电子数*1.5，具体可以去看看大师兄的教程。

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