Bader 电荷分析

Northword2021年1月21日
  • VASP
  • 电子结构
大约 2 分钟

Bader 电荷分析

Todo

理查德·贝德(Richard Bader)开发了一种将分子分解为原子的直观方法。他对原子的定义纯粹基于电子电荷密度。Bader 使用所谓的零磁通表面来划分原子。零通量表面是 2D 表面,其上电荷密度垂直于表面。通常在分子系统中,电荷密度在原子之间达到最小值,这是将原子彼此分开的自然位置。除了作为分子中原子可视化的直观方案外,Bader 的定义通常也可用于电荷分析。例如,Bader 体积内的电荷与原子的总电子电荷很接近。电荷分布可用于确定相互作用的原子或分子的多极矩。Bader 的分析也被用来定义原子的硬度,可以用来量化从原子中去除电荷的成本。

步骤

结构优化

自洽

LAECHG =.TRUE.
LCHARG = .TRUE.
NSW    = 0
IBRION = -1 (前面有了NSW = 0, 这个也可以不设置)

提示

AECCAR0AECCAR2 是由 LAECHG=.TRUE. 控制输出的。

Bader

chgsum.pl AECCAR0 AECCAR2   #output: CHGCAR_sum
bader CHGCAR -ref CHGCAR_sum

生成以下输出文件:ACF.datBCF.datAtomVolumes.dat

ACF.dat 包含每个原子的坐标,根据 Bader 分区与其相关的电荷,根据 Bader 分区的整体百分比以及到表面的最小距离。如果使用了伪电位,则应将该距离与核心区域的最大截止半径进行比较。

BCF.dat 包含每个 Bader 最大值的坐标,该体积内的电荷,最近的原子以及到该原子的距离。

AtomVolumes.dat 包含已分配给每个原子的每个卷的编号。

数据分析

Todo

参考资料