Partial charge density by VASP

VASP可以研究具体某个k点某条能带的局部电荷密度，并且把它投影到实空间。

具体的过程可以参考：

https://zhuanlan.zhihu.com/p/410514906

vasp官网上也给出了相关的参数解释：

https://www.vasp.at/wiki/index.php/Band_decomposed_charge_densities

但是如果打开了自旋，输出的电荷密度是上自旋+下自旋的，无法具体区分到某一个自旋轨道，在官网的wiki中也明确指出了：PARCHG中的第一部分是上自旋+下子旋，第二部分是上自旋-下子旋，所以只需要简单地处理一下就可以得到上下自旋分别的局域电荷密度。

写了一个简单的脚本，仅供参考，如果有问题敬请指正。

https://github.com/1241563482/Partial-Charge-Density

运行脚本，输入需要处理的文件名，就会得到up-PARCHG和down-PARCHG，分别对应两个自旋轨道。