Bader电荷为什么通常不是整数？

说明：本文华算科技主要介绍 Bader 原子盆电子数产生小数的原因、它与整数电子计数和氧化态的关系，以及分区电荷的数值收敛与比较条件。

Bader电荷中的电子数怎样定义？

一份标为“Bader 电荷”的结果，可能列出原子区域内的电子数 N_A，也可能列出扣除电子后的原子净电荷 q_A。两者的符号含义相反：区域内电子增加时，N_A 增大，q_A 减小。电子数密度 ρ(r) 描述单位体积内的电子数，积分区域和中性原子的参考电子数决定了输出数值。

Bader 分区沿电子密度的梯度轨线划分空间，通常把终止于同一原子核附近密度极大值的轨线归入一个原子盆 Ω_A。相邻盆之间满足∇ρ·n=0 的零通量条件，n 是分界面的法向单位矢量。对盆内密度积分，得到 N_A=∫_ΩAρ(r)d³r；全电子密度下，净电荷写为 q_A/e=Z_A−N_A，Z_A 是核电荷数，e 是正的元电荷。

零通量条件约束的是密度梯度的法向分量，没有要求盆内积分只能取整数。这里也没有电子运动的时间信息；分界面由静态密度的几何形状定义。电子密度在键区可以保持非零，成键轨道也可以跨越多个原子盆。将完整密度限制在一个盆内积分，所得电子数由区域大小和区域内的密度共同确定。

例如，对水分子作 B3LYP/cc-pVDZ 计算时，一组 QTAIM 结果给出氧的净电荷 −1.173 e，即氧盆内平均有9.173 个电子。两个对称等价的氢盆各分得约 0.4135 个电子，三盆电子数相加仍为 10。图中的 HF+J 电子构型只表示联合概率的一个极大位置，图内小球的数量与这组 B3LYP 平均布居分别描述不同统计对象。

整数电子怎样形成小数原子布居？

对固定空间区域作理想电子数计数，每次得到的是 0、1、2……这样的整数；量子态可以对多个计数结果赋予非零概率。若 P_A(n) 是原子盆内有 n 个电子的概率，密度积分给出的平均电子数满足 N_A=Σ_nnP_A(n)。各个 n 都是整数，加权平均却可以取小数。整个分子有确定的总电子数，仍允许各个原子盆的电子数存在量子涨落。

在单粒子轨道表示中，ρ(r)=Σ_if_i|ψ_i(r)|²。一个归一化轨道在完整空间的积分为 1，但在某个原子盆内的积分可以介于 0 和 1 之间。即使所有轨道占据数 f_i 都取整数，跨原子的轨道密度仍会给出小数布居。这种区域平均描述没有把一颗电子分割成几块，也没有赋予单个电子小于 e 的电荷。

CaH₂ 的全局极大构型中，每个氢附近有两个电子；另一些局部极大构型中，一个或两个氢附近只有一个电子。这些构型在量子态中具有各自的权重，平均布居包含整个分布的贡献，仅统计最可能构型的小球会丢失其余构型的信息。

电子展宽改变轨道占据和电子密度，在金属计算中会影响分区电荷；它只是一项计算条件。分子、绝缘体和整数占据的电子态也能产生小数 Bader 电荷，原子盆内的轨道积分已经允许这种结果。有限温度或分子动力学还可以引入结构随时间变化的平均，那属于另一种统计；静态结构中的小数布居本身不要求原子随时间振动。

整数氧化态与总电荷守恒分别约束什么？

NaCl 中的 Na 和 Cl 通常分别标为+1、−1 氧化态。这个计数采用化学键的离子近似，把异核键电子按归属规则分配给键合原子。Bader 方法则保留连续电子密度，再按零通量面积分。两套规则使用不同的电子归属对象：氧化态记录形式电子计数，原子盆净电荷记录指定密度分区后的电子增减。

NaCl 的一项网格研究中，Na 的 Bader 净电荷随网格加密逐渐收敛到约 +0.88 e。横轴使用网格点数的倒数，越靠左代表网格越密；两种积分方法都趋于一个小数值。晶体中的 Na⁺ 标签与这个小数结果可以同时使用，各自保留离子近似与实空间积分的定义。

当全部原子盆覆盖整个体系时，Σ_AN_A=N_total，相应有 Σ_Aq_A=Q_total。中性 NaCl 晶胞内，Na 的正电荷与 Cl 的负电荷相互抵消；守恒约束的是总和，每个原子的数值仍可保留小数。若存在非核吸引盆，电子积分的求和也须包含这些区域，遗漏区域会造成统计缺项。

把每个 Bader 数值四舍五入成整数，会删去配位、键区极化和结构变化引起的连续布居差，并可能破坏净电荷总和。数值接近某个整数，也没有自动给出氧化态归属。对称等价原子之间的均分有时恰好产生整数电子数，例如中性 H₂ 的两个等价氢盆各分得一个电子。

怎样判断Bader电荷是否已经收敛？

离散网格把连续密度写成有限个采样值。原子盆分界面穿过网格单元时，把整个单元归给一侧会产生积分误差；加权积分允许按单元体积分数分配贡献。这里的分数权重用于近似积分区域，与原子盆电子数的量子平均来自不同的数学操作。

判断收敛时，比较的是网格加密后同一个电荷值的变化，目标由待研究的电荷差决定。NaCl 研究对近网格法和加权法分别外推得到 +0.878 e、+0.881 e；两者相差 0.003 e，原文将差异联系到近网格法的区域误分配。数字保留到小数点后三位，没有自动赋予三位小数同样的计算精度；积分算法也参与数值误差。

赝势和 PAW 计算中，用于寻找原子盆的参考密度与实际积分的密度可以分别指定。用重构全电子密度确定分区，再积分价电子密度时，净电荷须用相应的价电子参考数 Z_val,A 减去盆内价电子数。若积分的是核电子与价电子的总密度，则使用核电荷数 Z_A；将全电子积分与价电子参考数相减，会引入整组核电子的偏差。

比较掺杂前后、不同吸附位点或不同晶面时，密度类型、价电子参照和积分协议应保持一致；泛函、U 值、结构和自旋态的变化则可能改变真实密度。对磁性体系，原子电子数来自两种自旋密度之和；自旋密度之差在同一区域内的积分统计局域自旋布居。