搜索
1.计算。在静态计算中，加入以下参数，其中NGX参数与POSCAR中的晶格参数成比例，太小bader计算不会收敛，所以，要收敛之后，结果才可靠。#baderchargeLAECHG=.TRUE.NGX=100NGY=100NGZ=200NGXF=150NGYF=150NGZF=3002.算完之后，生成文件AECCAR0AECCAR1AECCAR2将chgsum脚本上传到计算文件夹将其变成可执行脚本：chmod777chgsum.pl然后使用命令：./chgsum.plAECCAR0AECCAR2生成CHARGCAR_sum文件3.将chgsum脚本上传到计算文件夹然后将bader变成可执行文件：chmod777bader输入命令：./baderCHGCAR-refCHGCAR_sum生成文件：ACF.datAVF.datBCF.datACF.dat就包含价电子的电荷信息：4.判断计算的可靠性把ACF.dat文件拖到Origin中，CHARGE这一列的数字就是每个原子对应的电荷，而最下面那一行，NUMBEROFELECTRONS：126，代表整个体系的总电量有126e。我们可以通过对比总电量和所有原子加和电量之间的差值来判断是否收敛：我们把CHGCAR这一列数据复制到Excel中，选中所有数据，按求和符号，或者在另一个格子中输入=sum()，Enter，得到总电荷是125.9585126-125.9585=0.0415(e)0.0415/126=0.03%我们发现：CHGCAR这一列数据这和+VACUUMCHARGE这个数值=NUMBEROFELECTRONS而VACUUMCHARGE这个值是真空的电荷，如果这个值不为零，就意味这个有电子转移到真空中。标准的情况应该是：所以，判断bader计算是否收敛，也可以通过判断VACUUMCHARGE这个值是否为0来判断。5.应用：怎么用bader分析研究体系的电荷转移总量。首先，我们可以从ACF.dat文件的最后一行NUMBEROFELECTRONS知道这个体系的总电量。我们计算一个体系在电荷转移前后的bader电荷，然后在ACF文件中读取其电量。如果两个文件中的原子数相等，就可以直接作差，得到电荷转移的量，如果原子数不同，就拿对应原子的电量作差。