ab-initio - VASP 中的形成能间隙 He-W

Question

我想计算包含间隙 He 原子（八面体位置）的 W 超晶胞在 VASP 中的形成能 Ef。我运行程序，它返回非常高的值（15/25 eV）。这是我的印加：

System = W
PREC = accurate #cutoff + wrap around errors.
IBRION = 2
ISIF= 3 #whether stress tensor is calculated, what is allowed to relax.
NSW = 151
SIGMA = 0.19
ISMEAR = 1
ENCUT = 400 #eV
ISPIN = 1
EDIFF = 1E-5
EDIFFG = -0.025

这是我的 2x2x2 超级单体 POSCAR：

    critic2 | <derived>
1.0
   6.3696000000    0.0000000000    0.0000000000
   0.0000000000    6.3696000000    0.0000000000
   0.0000000000    0.0000000000    6.3696000000
He W
1  16
Direct
0.2500000000 0.2500000000  0.000000
   0.00000000    0.00000000    0.00000000
   0.25000000    0.25000000    0.25000000
   0.00000000    0.00000000    0.50000000
   0.25000000    0.25000000    0.75000000
   0.00000000    0.50000000    0.00000000
   0.25000000    0.75000000    0.25000000
   0.00000000    0.50000000    0.50000000
   0.25000000    0.75000000    0.75000000
   0.50000000    0.00000000    0.00000000
   0.75000000    0.25000000    0.25000000
   0.50000000    0.00000000    0.50000000
   0.75000000    0.25000000    0.75000000
   0.50000000    0.50000000    0.00000000
   0.75000000    0.75000000    0.25000000
   0.50000000    0.50000000    0.50000000
   0.75000000    0.75000000    0.75000000

与其他作者相比，我期望值接近 5/6 eV。