Bader 电荷分析
大约 2 分钟
Bader 电荷分析
Todo
理查德·贝德（Richard Bader）开发了一种将分子分解为原子的直观方法。他对原子的定义纯粹基于电子电荷密度。Bader 使用所谓的零磁通表面来划分原子。零通量表面是 2D 表面，其上电荷密度垂直于表面。通常在分子系统中，电荷密度在原子之间达到最小值，这是将原子彼此分开的自然位置。除了作为分子中原子可视化的直观方案外，Bader 的定义通常也可用于电荷分析。例如，Bader 体积内的电荷与原子的总电子电荷很接近。电荷分布可用于确定相互作用的原子或分子的多极矩。Bader 的分析也被用来定义原子的硬度，可以用来量化从原子中去除电荷的成本。
步骤
结构优化
自洽
LAECHG =.TRUE.
LCHARG = .TRUE.
NSW = 0
IBRION = -1 (前面有了NSW = 0, 这个也可以不设置)
提示
AECCAR0 和 AECCAR2 是由 LAECHG=.TRUE. 控制输出的。
Bader
chgsum.pl AECCAR0 AECCAR2 #output: CHGCAR_sum
bader CHGCAR -ref CHGCAR_sum
生成以下输出文件：ACF.dat，BCF.dat，AtomVolumes.dat。
ACF.dat 包含每个原子的坐标，根据 Bader 分区与其相关的电荷，根据 Bader 分区的整体百分比以及到表面的最小距离。如果使用了伪电位，则应将该距离与核心区域的最大截止半径进行比较。
BCF.dat 包含每个 Bader 最大值的坐标，该体积内的电荷，最近的原子以及到该原子的距离。
AtomVolumes.dat 包含已分配给每个原子的每个卷的编号。
数据分析
Todo