VASP算bader电荷

**可以用来判断电子的得失情况**
需要用到chgsum.pl
这个脚本在vtst中有，是有一个算过渡态的脚本合集，这里我们不安装，直接解压缩安装包
[安装包下载点击](http://pan.liuyaoze.com/s/9Vfe "安装包下载点击")
解压缩 ` tar -xvzf vtstscripts.tgz`
在压缩包里可以找到chgsum.pl

![](/media/202312/2023-12-19_110810_0379130.2625695245877784.png)
这个时候这个脚本由于权限问题还不能执行
我们把这个文件夹的路径加到环境里
`vi ~/.bashrc`

![](/media/202312/2023-12-19_111809_5393430.9221513102284354.png)
激活
`source ~/.bashrc`
这个时候输入chg,再用tab键就可以补全出chgsum.pl了

看电荷算静态时需要加上参数
```
Bader Charge Analysis
LAECHG     = .TRUE.    (Write core charge into CHGCAR file)
LCHARG     = .TRUE.    (Write CHGCAR file)

```
得到AECCAR0 AECCAR1 AECCAR2文件
`chgsum.pl AECCAR0 AECCAR2 `
`bader CHGCAR -ref CHGCAR_sum`
得到ACF.dat文件，就有每个原子的bader电荷
这里需要跟原子本身最外层的电荷相比较来判断电子得失
`grep ZVAL POTCAR`

![](/media/202312/2023-12-19_111452_3455840.8495441675786674.png)

![](/media/202312/2023-12-19_111531_3481280.6982892454944095.png)
这里的Zr2O3就是两个Zr原子失去4个电子，变成Zr2+三个氧得到4个电子

自洽后，后处理Python代码:
```python

import numpy as np
import os

os.system('chgsum.pl AECCAR0 AECCAR2')
os.system('bader CHGCAR -ref CHGCAR_sum')

ff = open('ACF.dat', 'r')
a = ff.readlines()[2:-4]
#print(a[2:-4])

charge = []
for i in range(len(a)):
    charge.append(float(a[i].split()[4]))
#    print(a[i].split()[4])
print(charge)

f2 = open('POSCAR', 'r')
res = f2.readlines()
atoms = res[5].split()
num = list(map(int,res[6].split()))
print(atoms)
print(num)

ZVAL = []
with open('POTCAR', 'r') as f3:
    for line in f3:
        res = line.split()
        if 'ZVAL' in res:
            ZVAL.append(float(res[5]))
print(ZVAL)

f4 = open('result.txt', 'w')

#for i in range(len(num)):
#    for j in range(num[i]):
#        charge[]
atoms_new=[]
ZVAL_new = []
for i in range(len(num)):
    for j in range(num[i]):
        ZVAL_new.append(ZVAL[i])
        atoms_new.append(atoms[i])

diff = np.array(charge)-np.array(ZVAL_new)

for i in range(len(diff)):
    f4.write(atoms_new[i]+' '+str(diff[i])+'\n')
f4.close()
#print(ZVAL_new)

```

VASP算bader电荷

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱
work@liuyaoze.com

VASP算bader电荷

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱 work@liuyaoze.com

注册码获取邮箱
work@liuyaoze.com