如何使用vasp计算Pz-band-center

zravenclaw

求助一下老师们，我现在在重复一篇文献，别的问题都解决了，但是d图中的Pz-band-center不知道怎么算，请老师们帮忙看一下

卡开发发

没有详细的文献信息？我猜测应该是直接用体系Pz的分波直接算带中心，当然前提是如果体系本身可能朝向符合通常化学上那样的定义，即Pz轨道的原子都在一个xy平面上。

但是问题又来了：
1、进行带中心的计算时，积分上下限如何确定？比如更大的空带和PDOS计算范围这个Pz带中心会收敛吗？还是会无止境上移或下移？
2、作者给的PDOS是展宽过后的，因为使用了展宽导致没办法判断Pz到底某些地方有没有gap，假定gap就在CBM和VBM之间，这时候积分上下限到底是只取到VBM/E_fermi？还是更高的能级都有算进去？
这些可能都会导致复现不了文献。