能带分解电荷密度

Band decomposed charge density：（VASP）

首先，所谓的Partial Charge是针对与Total Charge来说的，指的是某个能量范围、某个K点或者某个特定的态所对应的电荷密度。在文献中最常见的是价带顶部，导带底部，表面态或者局域态所对应的 Partial Charge。通过分析这些态所对应的Partial Charge，可以得到体系的一些性质，比如局域态具体的是局域在哪个原子上等。进行Partial Charge Analysis的第一步是进行自洽的计算，得到体系的电子结构。这一步的计算采用通常的INCAR和KPOINTS文件。在自洽计算结束后，保存WAVECAR文件。

第二步是画出能带结构，以决定你需要画哪条能带的那个K点的态所对应的Partial Charge。决定画哪条能带，或者那些感兴趣的K点之后，有如下几种方法计算不同的Partial Charge。如果计算价带顶端的Partial Charge，则需要首先通过能带结构图确定价带的能带标号。需要注意，进行Partial Charge分析必须要保留有自洽计算的WAVECAR才可以。

第一种Partial Charge分析的INCAR，这样的INCAR给出的是指定能带，指定K点所对应的Partial Charge。分析导带、价带等的Partial Charge特性，通常采用的都是这种模式。

ISTART = 1 （job : 0-new 1-cont 2-samecut
ICHARG = 1 charge: 1-file 2-atom 10-const）
LPARD=.TRUE.
IBAND= 20 21 22 23 （能带的序号，即第几条）
KPUSE= 1 2 3 4 （某段能带（两个高对称点间）之间的第几个K点，在KPOINTS中需要输入这两个高对称点的坐标）
LSEPB=.TRUE.
LSEPK=.TRUE.

第二种Partial Charge分析的INCAR，这样的INCAR给出的是在[-10.3 -5.1]能量之间的Partial Charge。这种模式适合于分析某个能量区间内的波函数的性质。
ISTART = 1 （job : 0-new 1-cont 2-samecut
ICHARG = 1 charge: 1-file 2-atom 10-const）
LPARD=.TRUE.
EINT = -10.3 -5.1
LSEPB=.FALSE.
LSEPK=.FALSE.

第三种Partial Charge分析的INCAR，这样的INCAR给出的是从[Ef-1.0 Ef]能量之间的Partial Charge。这种模式最利于分析费米面附近的波函数的性质。
ISTART = 1 （job : 0-new 1-cont 2-samecut
ICHARG = 1 charge: 1-file 2-atom 10-const）
LPARD=.TRUE.
NBMOD=-3
EINT = -1
LSEPB=.FALSE.
LSEPK=.FALSE.

计算能带的部分电荷密度时，除了所需的四个文件（INCAR、KPOINTS、POSCAR、POTCAR），还需将WAVECAR、CHGCAR拷入文件夹。在INCAR中加入上述第一种情况的那些参数，在KPOINTS中修改成该条能带中某两个高对称点的坐标。