Fe(110)晶面吸附计算电子步难收敛问题如何解决

玄武一甩 · 发表于 Post on 2022-1-6 11:20:24

最近在做Fe(110)晶面的结构优化、吸附计算工作，发现Fe（110）电子步很难收敛，通常200步以内无法收敛，目前的晶面结构是4层，4*4，共256个原子，INCAR文件如下如所示，想问下大家是否有参数需要添加或修改，来加快电子步收敛

Aristotler · 发表于 Post on 2022-1-6 11:53:59

SIGMA=0.1 或者等于0.2试一下或者ALGO=48；你的问题在于Fe你用了自旋，并且截断能太大了（400或者450足以）

Aristotler · 发表于 Post on 2022-1-6 11:55:47

本帖最后由 Aristotler 于 2022-1-6 11:57 编辑

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

玄武一甩 · 发表于 Post on 2022-1-6 12:16:10

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

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

Para1lel · 发表于 Post on 2022-1-6 16:14:24

Fe有磁性，设个初始磁矩比较好，MAGMOM这个参数设一下。

wsz · 发表于 Post on 2022-1-7 02:38:05

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

玄武一甩 · 发表于 Post on 2022-1-7 20:33:48

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

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

Para1lel · 发表于 Post on 2022-1-8 08:42:09

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

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

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

[VASP] Fe(110)晶面吸附计算电子步难收敛问题如何解决

浏览过的版块