求助计算态密度的切面和真空层设置

nanokey

计算二氧化钛衍生物的态密度，切了200晶面，表面加了20埃的真空层，然后计算态密度，加U和加了偶极矫正IDIPOL=3，但是结果还是除了问题，在费米能级附近存在新的峰，原胞计算不存在这种情况，具体参数和结果如下，请教各位到底哪里出了问题，谢谢！
SYSTEM = T200
ISTART =1
ICHARG = 11
LREAL = Auto
ENCUT = 400
PREC = Accurate
LWAVE = .FALSE.
LCHARG = .FALSE.
LORBIT=11
EMIN = -20
EMAX = 15
ISMEAR = -5
SIGMA = 0.05
NELM = 300
NELMIN = 4
EDIFF = 1.E-4
NSW = 0
IBRION = 2
ISIF = 2
EDIFFG = -0.05
ISYM = 0
POTIM = 0.3
IDIPOL=3
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 2 -1
LDAUU = 0 4 0
LDAUJ = 0 0.4 0
LDAUPRINT = 0
LMAXMIX = 4

喵星大佬

表面和块体的DOS不一样不是很正常嘛

nanokey

喵星大佬 发表于 2021-10-11 12:20
表面和块体的DOS不一样不是很正常嘛

我想消除这些小峰，请问应该怎么设置计算呢？