1， 在HSE-dos计算的文件夹，新建两个文件夹，分别vb, cb;

2，进入vb, 将上一级文件夹的WAVECAR，CHGCAR，KPOINTS，INCAR，POSCAR，POTCAR都拷贝过来，INCAR修改如下：

ISTART = 1

ICHARG = 1

LPARD=.TRUE.

NBMOD=-3

EINT = -1 #表示从Ef-1到Ef

LSEPB=.FALSE.

LSEPK=.FALSE.

其它参数保留。

3，计算，获得band-decomposed CHGCAR.  拖到vesta作图；

同理，进入cb

4. 将HSE-dos计算的文件夹的WAVECAR，CHGCAR，KPOINTS，INCAR，POSCAR，POTCAR都拷贝过来，INCAR修改如下：

ISTART = 1

ICHARG = 1

LPARD=.TRUE.

NBMOD=-3

EINT = 5 #表示从Ef到Ef+5

LSEPB=.FALSE.

LSEPK=.FALSE.

其它参数保留。

5，计算，获得band-decomposed CHGCAR.  拖到vesta作图；

---

还有其它两种方式，供参考

1. 指定K点及能带所对应的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

LSEPB=.TRUE.

LSEPK=.TRUE.

2. 指定能量区间的Partial Charge。这种模式适合于分析某个能量区间内的波函数的性质。需要注意的是，需要仔细看一下价带顶和导带底位置的具体能量值，价带顶可能不在0处。

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.