abacus+wannier90拟合64个水分子的能带拟合结果很差

liubing

这个是我的wannier.win文件
num_bands         =   512
num_wann          =   384

dis_win_min       =  -10.4
dis_win_max       =  23.2
dis_froz_min      =  -10.4
dis_froz_max      =  16.67

num_iter      =  0

exclude_bands = 1-64

begin projections
O:sp3
H:s
end projections


后续我也会给出pdos和eig来佐证我的投影子和窗口 pdos这个图片上传不了... 但是大概是这样的由pdos可知64条底带是由o的s态提供因此我在O投影子的选择上p和sp3（H的投影子始终为s）均进行了尝试但是效果都不理想  同样对于迭代次数也进行了计算 0次和100次的wannier迭代 效果也是不行 现在有点束手无策了 遂求助各位 不知道有没有熟悉wannier计算的大佬能指点一二 感谢各位大佬
  63    1  -19.599974724586
   64    1  -19.480380894019
   65    1  -10.266680876607
这个是eig文件我是在第65条能带开始计算的 算到64+384=448能带
447    1   16.647049410482
  448    1   16.812133391141
  449    1   16.886232567784

大概是这里 然后我的num_bands这个部分的设置我也有点小疑惑 就是网上的教程说就在目标能带往上拓展两个能带就可以 还有的说这个部分要包含所有投影子有贡献的能带 但是我选择的是原子轨道基来进行计算 要是把所有投影子都有贡献的能带都选取的话 总能带太多了 不知道这个部分的设置要怎么能弄
ps：在.wout文件中也可以看出来整个体系的离散化不好 而且很奇怪的是 当我选择sp3投影的时候离散化不好的是H原子轨道 选择p3投影子时候离散化不好的是O原子轨道 也非常的奇怪 在某些轨道中心我的spread都可以高达10+