VASP软件对电荷相关的计算和处理有多个方面，这篇文章介绍基本的电荷密度、Bader电荷以及电子局域函数的计算和数据处理方法。

1. 文件介绍：

• CHG、CHGCAR: charge density

• AECCAR0, AECCAR1, AECCAR2： core charge and valence chage

• ELFCAR：electron localization function

2. 计算方法：

通过一步静态计算即可获取以上所有文件，相关INCAR关键参数设置：

• LCHARG=.TRUE. （VASP默认.TRUE.），将电荷密度写入CHG, CHGCAR文件

• LAECHG = .TRUE. (VASP默认.FLASE.)， 得到包含charge信息的AECCAR0（core charge）, AECCAR1, AECCAR2(valance charge)

• LELF = .TRUE. （VASP默认.FALSE.），得到写有电子局域函数的文件ELFCAR。

3. 结果处理软件（使用方法参考软件附带链接）：

•  处理CHGCAR文件推荐软件: LEV00 、VESTA

• VESTA直接根据系统版本选择安装使用，不具体介绍

• LEV00:通过链接下载，根据机器安装的数学库和编译器修改Makefile中的“LIBS”和“FCOMPL”两项（一般用ifort即可），并进行安装（执行make），得到lev00可执行程序（根据需求设置lev00的环境变量，这里不赘述）。
  

lev00使用：

直接运行可执行程序lev00,根据提示设置选项即可。运行完成后，可以得到.dat或.gnu文件，可以通过Origin、gnuplot等软件绘图。

注意事项：（1）除了用到CHGCAR文件，还需要读取OUTCAR，POSCAR(或CONTCAR)文件。（2）如果OUTCAR文件中包含用虚线包围的警示信息，需要删除这些警示信息后再使用，否则会遇到读取文件的错误。（3）lev00除了可以处理CHGCAR以外，还可以用于DOS、coulomb potential等VASP计算的数据处理，将会在后面具体介绍。  

• 处理AECCAR文件推荐软件: Bader 

• 安装：Bader软件可以通过链接直接下载到ifort编译的可执行程序（bader），也可以下载源代码。（设置bader的环境变量)

• 软件执行命令(两步即可完成)： 需要借用VTST中的处理脚本chgsum.pl 

• chgsum.pl AECCAR0 AECCAR2  (产生CHGCAR_sum文件)

• bader CHGCAR -ref CHGCAR_sum （得到ACF.dat, BCF.dat以及AVF.dat三个文件）

ACF.dat文件可以获取每个原子的bader电荷。        

• 处理ELFCAR文件推荐软件: VESTA  根据系统版本选择安装使用。

软件相关链接：

LEV00:  https://nms.kcl.ac.uk/lev.kantorovitch/codes/lev00/index.html

VESTA：https://jp-minerals.org/vesta/en/

Bader ：http://theory.cm.utexas.edu/henkelman/code/bader/

VTST：http://theory.cm.utexas.edu/vtsttools/scripts.html