VASP能级对齐问题

北大-陶豫

请教一下各位大佬，我计算了几个不同体系的 DOS 和 PDOS，有没有什么办法把它们的能量对齐呢？我想看一下是不是不同体系中的不同结构导致碘的 5p 轨道能量改变，改变情况是上升还是下降，这样。应该怎么操作比较好呢？VASP 里面给出的能级，是不是把 Fock 矩阵对角化得到的能级，减去 Fermi 能级，得到的？我们能不能跳过“减去 Fermi 能级” 这步，直接得到 Fock 矩阵对角化得到的能级，然后将这一数值在不同体系中比对？
我在网上看到一种说法说，可以根据真空能级或者深能级芯电子对齐。但我还有问题：算真空能级的话要重新修改 POSCAR，计算表面而不是体相结构，打开LOCPOT得到真空区的 KS potential，是这样吗，为什么需要这么麻烦而不能像我之前说的直接用 Fock 矩阵对角化得到的能级呢？芯电子的话好像直接被赝掉了没有能级，最低的价电子能级应该还是会轻微受到化学环境的影响的吧，可以用来比较吗？
谢谢！

北大-陶豫

自己写了一个脚本来处理这个问题。见附件。

脚本会找到体系中能量最低的几条能带（默认是最低的宽度大于2eV的禁带以下的所有能带，这一标准可以通过修改代码中的变量ForbiddenBand_MinWidth来更改），将这些能带在Γ点处本征值的平均值，作为体系的能量零点（便于不同体系之间比较）， 然后重新计算所有能带的能量，输出到输出文件（默认是BAND_0.0.dat）。


详见代码开头的注释。