求助VASP能带分解计算的方法

Kmetsch

新人初次接触VASP计算，尝试复现文献工作。
查阅的文献在如图a的体系中进行了能带的分解计算，如文献所述：

(b)Decomposed atomic band structure of Nbm-MoS2. Blue and red dots represent the bands decomposed on the orbitals of MoS2 (leftpanel) and the N-benzylmaleimide molecules (right panel), respectively.

是不是说作者对修饰后的复合体系的能带进行了计算，然后在后处理步骤中将MoS2的部分和修饰分子的部分进行了分解，再以不同颜色绘制成图b的样子？
但是参考vasp手册和vaspkit的例子，能带计算需要在体系的高对称性K点下，读取静态计算的CHGCAR进行非自洽计算。对于这样的体系，六方MoS2自身的对称性已经被修饰分子破坏了，是如何获取其高对称K点进行计算的呢？
另外，参考vaspkit的例子（MoS2的S空位掺杂），对于超晶胞的能带计算，需要进行反折叠处理来获得干净的能带图，在如图a这样的体系中，文献采用了3*3*1的超晶胞，是不是也需要进行反折叠计算呢？
文献doi: 10.1021/acsnano.2c07670

北大-陶豫

在能带曲线的每个点，可以计算该点的波函数里各个原子的贡献权重，这些信息保存在PROCAR文件里。通过vaspkit等脚本可以提取，绘制你这种“能带投影图”。vaspkit官方QQ群（498683602）的群文件里有画图模版。