求助计算Eu掺杂二氧化钛能带时出问题

夜凯

第一步结构优化：NSW=100  IBRION=2  ISIF=3  ICHARGE=2
输出的CONTCAR调整对称性并改成第二部的POSCAR
第二步：NSW=100  IBRION=2  ISIF=2  ICHARGE=2 输出的CONTCAR改后续计算的POSCAR
第三步：ICHARGE=2  NSW=0  IBRION=-1  KPOINTS加密  输出CHGCAR和WAVECAR
前三步：KPOINTS是Mnkhorst-Pack (第二张图)
第四步：ICHARGE=11  NSW=0  IBRION=-1  ISTART=1 并复制第三步输出的CHGCAR和WAVECAR
第四步的KPOINTS是(第三张图)




我画的能带图
文献画的能带图
请求各位大神帮忙看看，指出我的问题所在。