Bader电荷分析全解：密度分区与电子数积分、配位基团统计及带电界面参照态

说明：本文华算科技主要介绍 Bader 电荷分析的密度分区与电子数积分、配位基团的电荷统计、带电界面的参照态，以及密度表示和网格对结果的影响。

Bader分析的基本定义是什么？

分子或晶体成键后，电子密度连续分布在整个空间，原子之间没有预设的球面分隔。DFT 求得的电子数密度 ρ(r)记录单位体积内的电子数；Bader 分析沿密度梯度划分空间，再统计每个区域含有多少电子。它处理的是已获得的电子密度，原子电荷并非额外一次自洽计算求出的独立电子结构。

从空间一点沿∇ρ 的上升方向延伸，终止于同一密度极大值的点组成一个吸引盆，普通分子中的极大值通常位于原子核附近。相邻盆之间的分隔面满足∇ρ·n=0，n 为面法向。零通量指密度梯度的法向分量为零，面上的 ρ 可以非零；梯度线也没有描述电子随时间运动的轨迹。

对原子 A 的盆 Ω_A 积分，得到N_A=∫_ΩAρ(r)d³r。在全电子表示中，净电荷写为 q_A=e(Z_A−N_A)，Z_A 是核电荷数，e 是正的元电荷。N_A 与 q_A 是两个量：前者为盆内电子数，后者包含原子核的正电贡献。计算结果表中的“charge”可能使用其中任一种约定，正负号必须由定义确定。

碳的全电子基准为 Z=6，盆内电子数少于 6 时净电荷为正，多于 6 时为负。非整数电子数来自对连续密度作空间积分；一个电子的量子态可以分布在多个盆内。常规原子盆覆盖全部积分空间时，ΣN_A 恢复模型电子总数，Σq_A 恢复相应子系统的总电荷；若另有真空区域或非核吸引盆，其积分也应参与总量核对。

配位结构中的原子净电荷怎样比较？

稀土分离材料 NCU-1 的计算模型中，Pr³⁺、Dy³⁺ 和 Al³⁺ 与孔内含氧基团配位，金属原子的 Bader 净电荷分别为 +2.31e、+2.17e 和 +2.09e。形式价态 +3来自电子归属的化学计数规则；三个非整数值来自各自密度分区后的电子数积分。不同元素的核电荷和价电子基准不同，比较时应使用已换算的净电荷。

金属附近的差分密度等值面描述相对选定组分参照的积累与耗尽，Bader 数值则来自整个金属盆的积分。等值面覆盖的体积还受所选等值大小影响，不能按色块大小换算电子数。若要讨论配位后金属获得多少电子，比较对象应为吸附态与明确定义的参照态；直接用 3 减去上述净电荷，会把形式价态计数与空间积分混用。

三种离子的吸附能依次为 −3.29、−2.58 和 −1.66 eV，分别包含对应模型的成键、形变及参照能贡献。净电荷记录电子分配，吸附能记录能量差，二者没有通用的单调换算式。Pr 与 Al 的净电荷差还涉及不同元素和配位环境，不能将差值写成同一金属位点的得失电子量。

对于含多个原子的吸附物或配体，可将选定原子的盆作为一个片段。硫酸根基团的Q(SO₄)=q_S+Σ_i=1…4q_O,i统计硫和四个氧，未计入相邻金属。钾镁矾中，氧盆与硫、镁、钾的盆相接，形状由晶体电子密度决定；将硫和氧的盆合并后，所得基团区域也没有球形要求。

SO₄²⁻ 的形式电荷为 −2，晶体中这一组原子盆的积分净电荷可以偏离 −2e。片段名单固定后，不同配位状态的 Q 差值才具有相同统计对象。原子盆随密度重新划分，电荷差包含盆内密度和盆形的共同变化，不能由它追溯某个电子经过哪个原子。

带电界面中的电子来自哪里？

固定电子总数的中性吸附模型中，吸附物电子数增加时，其余区域必须补偿相应变化。电化学界面若与电子库交换电子，则金属和水层可能同时增电子。恒电位计算通过调整电子数维持所选电极电位，模型中的电荷变化还包括外电路供电子，不能全部归为水层与金属之间的内部转移。

溶剂化 jellium 模型以连续补偿电荷平衡显式体系的附加电荷，并通过溶剂介电响应屏蔽电场。整个计算单元电中性时，显式原子盆的净电荷和仍可非零：补偿背景属于模型中的另一项电荷。若将背景遗漏，却要求所有原子净电荷之和为零，就会误判积分结果。

Au(111)/水界面可分别求和金属片段与水层片段的 Bader 电子数，并将过量质子归入水层。以功函数 4.4 eV 的状态为参照，ΔN_e 为所选片段相对该状态的电子数变化；在原子名单和核电荷固定时，ΔQ=−eΔN_e。正的 ΔN_e 表示参照态以上的电子积累，对应净电荷向负方向变化。

左侧红色实心点由 DFT 的逐原子 Bader 结果求和；黑色空心点来自机器学习预测的原子电荷。水层与金属的增量随电位改变而不同，右侧散点比较两种模型的片段增量偏差。这里比较的是指定 Volmer 初态的静态电荷响应；电流或反应速率则涉及电子转移、质子运动与活化自由能，静态 ΔN_e 没有时间分母。

数值误差会怎样进入原子电荷？

赝势或 PAW 计算常将冻结芯电子与价电子分别处理。平滑价电子密度在核附近的形状可能缺少全电子密度的峰。含芯电子贡献的重构密度可用于确定原子盆，再对指定价电子密度积分；此时净电荷采用 q_A=e(Z_A,val−N_A,val)，Z_A,val 为对应势的价电子基准。若积分包含全部电子，基准才改用核电荷数。

用于划分空间的参考密度与实际求积分的密度可以不同，二者承担的计算任务须匹配。全电子基准配价电子积分会产生包含芯电子数的假电荷。比较不同 PAW 数据集时，即使元素相同，价电子基准也可能改变；相同的输出积分值未必对应相同净电荷。

连续积分在计算中写为 N_A≈Σ_j∈ΩAρ_jw_j，w_j 是体积或数值积分权重。网格点的归属与积分权重共同影响结果。苯的分区图在分子平面附近保持六重对称，包含远离原子核的稀疏网格后，白圈中的部分点发生归属错误；颜色标识不同原子盆，颜色本身没有电荷正负含义。

该粗网格例子中，误分约占网格点的 1%，但主要发生在低密度区域，对盆内电子数的误差约为 10⁻³ 个电子。点数占比与电子数误差采用不同权重，不能把 1% 的点误分直接记成 1% 电荷误差。积分收敛还受梯度追踪算法影响，限制上升方向的网格偏差可能在加密后仍然保留。

收敛曲线的纵轴接近 5.95，按原论文定义是碳盆的平均电子数；全电子净电荷应再按 6−N_C 换算。橙线对应均匀网格，蓝线对应适应原子附近密度变化的积分网格，黑线是另一积分实现的参照。两类网格在粗采样时的偏差不同，不能将某一图中的总点数当成所有材料通用的收敛阈值。

若两个吸附位点的净电荷只差几百分之一个 e，网格加密造成的电荷变化应小于待比较差值。比较前后的密度泛函、芯价处理、自旋态和结构参照须一致；研究电位响应时还需指定同一电位标尺与片段名单。否则，净电荷差可能来自密度表示或模型电子数的改变，而非同一配位环境下的电子重排。