分析Bader电荷是否必须考虑core charge?

archer

http://theory.cm.utexas.edu/henkelman/code/bader/
用Henkelman的脚本计算VASP的Bader电荷：
可以直接bader chargefile，或者chgsum.pl AECCAR0 AECCAR2；bader CHGCAR -ref CHGCAR_sum
最终结果输出到ACF.dat
后者考虑了Core charge和valance charge，网上的教程都写的是后者，并且电子得失用POTCAR的ZVAL减去ACF.dat的电荷。
分析Bader电荷是否必须考虑core charge? 为什么后者考虑了Core charge，只减去价电子ZVAL？

卡开发发

注意-ref的参数，你实际上可以发现计算过程还是加了AECCAR0和AECCAR2，这其实还是意味着使用了冻芯近似下的全电子密度。也就是说，实际上，处理bader电荷（AIM电荷）分为两个过程，1、划分原子盆，2、原子盆积分，如果原子盆划分的情况我们使用赝密度就会导致虚假的核临界点，因此划分应当采用冻芯近似下的全电子密度，也就是你的AECCAR0和AECCAR2合成出来的那个，而得到各个盆的区域后，实际在各个盆的积分采用的是CHARGCAR，也就是赝密度所对应的网格。因此这样积分得到总电子数是和赝势对应的ZVAL求和是一致的。

archer

卡开发发 发表于 2021-6-16 22:19
注意-ref的参数，你实际上可以发现计算过程还是加了AECCAR0和AECCAR2，这其实还是意味着使用了冻芯近似下的 ...

谢谢您的回复！我应该明白了，-ref 就是指定原子盆的划分以全电子密度为依据。