计算化学公社

标题: VASP对Fe3C的Slab模型优化使用111的K点原子乱飞，使用221的K点则计算正常 [打印本页]

作者
Author: KangBao 时间: 2024-7-4 10:59
标题: VASP对Fe3C的Slab模型优化使用111的K点原子乱飞，使用221的K点则计算正常
我目前的工作流程：
1网上获取Fe3C的晶胞进行ISIF=3的结构优化，晶格常数与实验值对比及文献误差均在3%以内；
2.将优化后的Fe3C晶胞用MS切面，扩建超胞，添加15A真空层构建Slab模型，模型大小为13 11 21；
3.仅放开表面3层原子，其余固定，对Slab模型进行结构优化用于表面吸附计算。
我在Slab模型优化时发现了K点的问题，当我使用221的K点得到正常优化结果，但消耗机时很长，约8h，之后的吸附计算耗时就更长了，约15h，明显不合理，我咨询了导师及师兄，导师说111K点足够了，但当我尝试111 计算时，无论调整INCAR中的ENCUT,IBRION,POTIM等或是增加真空层，均会出现C原子向上大幅乱移的现象，检查了赝势等也没有问题。
下面展示：优化前后结构图，优化的INCAR，POSCAR，KPOINTS。
我查询了论坛，没有找到有关K点对结构影响的讨论，期待老师们的回复

第一次发帖经验不足，如果有做的不好的地方请老师们指出，期待老师们的回复！

SYSTEM=Fe3C
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.05
ALGO=N
NELM=200
ENCUT=450
EDIFF=1E-5

IDIPOL=.TRUE.
IDIPOL=3

LWAVE=.F.
LCHARG=.F.

IBRION=2
#POTIM=0.1
EDIFFG=-0.02
NSW=200
ISIF=2

NCORE=12

KPOINTS
0
Gamma
2 2 1
0 0 0

POSCAR
1.0
   13.2568998337       0.0000000000       0.0000000000
      0.0000000000       10.7523002625       0.0000000000
      0.0000000000       0.0000000000       20.9848995209
C Fe
24 72
S
Cartesian
1.657112479 7.448225726 0                F F F
8.285562396 7.448225726 0                F F F
4.971337438 0.896526783 0.523153551       F F F
11.59978735 0.896526783 0.523153551       F F F

作者
Author: KangBao 时间: 2024-7-4 11:02
我在使用CINEB计算表面负载原子迁移能垒时发现，当我使用K点111的时候，第一步的力有8点多，但是我用221的时候力只有0.0几，K点的影响有这么大吗？

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