VASP对Fe3C的Slab模型优化使用111的K点原子乱飞，使用221的K点则计算正常

KangBao

我目前的工作流程：
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

KangBao

我在使用CINEB计算表面负载原子迁移能垒时发现，当我使用K点111的时候，第一步的力有8点多，但是我用221的时候力只有0.0几，K点的影响有这么大吗？