（二）Bader电荷分析

https://blog.sciencenet.cn/blog-3469498-1280199.html

具体方法如下：

一、INCAR中加入参数LAECHG = .TRUE.进行全电子计算。（可同时加入LELF=.TRUE.进行ELF计算）

LELF=.TRUE.
LAECHG=.TRUE.

LVTOT=.TRUE (Write total electrostatic potential into LOCPOT or not)

二、基于脚本chgsum.pl获得全电子电荷密度文件CHGCAR_sum。

chmod +x chgsum.pl

./chgsum.pl AECCAR0 AECCAR2

三、在CHGCAR和CHGCAR_sum文件中加入ISQ的坐标（用H原子表示）。

FCC-Al                                                        
   1.00000000000000    
     3.457999    0.000000    0.000000
     0.000000    3.457999    0.000000
     0.000000    0.000000    3.457999
   Al  H
     4  8
Direct
  0.000000  0.000000  0.000000
  0.000000  0.500000  0.500000
  0.500000  0.000000  0.500000
  0.500000  0.500000  0.000000
  0.250000  0.250000  0.250000
  0.250000  0.750000  0.250000
  0.250000  0.250000  0.750000
  0.750000  0.250000  0.250000
  0.750000  0.750000  0.250000
  0.750000  0.250000  0.750000
  0.250000  0.750000  0.750000
  0.750000  0.750000  0.750000

四、基于Bader分析小软件得到各个原子/ISQs的电荷量。

chmod +x bader./bader CHGCAR -ref CHGCAR_sum

 其中CHARGE就是Al原子和间隙ISQ的电荷量，ATOMIC VOL是Al原子和间隙ISQ的体积

注意事项：要保证各个原子/ISQs坐标的格式的统一。

全电子计算开关LAECHG = .TRUE.的说明：http://blog.sciencenet.cn/blog-2827406-1198494.html

Bader电荷分析的软件及其命令：http://theory.cm.utexas.edu/henkelman/code/bader/

chgsum.pl

bader