Pd(111)ORR计算遇到问题

goosecar

各位老师们，我在计算微弱晶格拉伸的Pd111面在进行ORR计算的时候遇到了slab模型和吸附态模型的能量差距非常大的问题。下面详述具体计算步骤：
      首先我先从网站上下载了Pd的初始晶胞，并进行优化和切面扩胞后，在VESTA中进行了约为0.86%的晶格拉伸(实验得到的数据)，然后ISIF=2，固定底部两层原子进行优化。起初我在结构网站上看到Pd是有微弱的磁性，所以ISPIN=2进行优化（但是后来我百度和看帖子发现Pd是没有磁性的）。目前能量是正常的（-315eV左右），但是在吸附一个氧原子达到收敛标准的时候能量是-340eV左右，差距非常大。然后思考了一下会不会是添加了磁矩的影响，但是把ISPIN=1试了一试，发现还是老问题。下面是我的POSCAR和INCAR。注：前后改变了ispin，结构并能没有跑散。
slab模型INCAR
ISTART=1
ICHARG=1
ISPIN=2
PREC=N
LORBIT=10
ALGO=F
NELM=80
EDIFF=1E-5
ENCUT=520
LREAL=A

IBRION=2
NSW=400
ISIF=2
EDIFFG=-0.05

ISMEAR=0
SIGMA=0.05
吸附模型INCAR
ISTART=1
ICHARG=1
ISPIN=2
PREC=N
LORBIT=10
MAGMOM= 64*0.39 0
ALGO=N
NELM=80
EDIFF=1E-5
ENCUT=520
IVDW=11
LREAL=A

IBRION=2
NSW=200
ISIF=2
EDIFFG=-0.05

ISMEAR=0
SIGMA=0.05
求助各位老师，谢谢

goosecar

这是吸附态的O的POSCAR