Bader电荷VS差分电荷VS Mulliken电荷

引言

Bader电荷、差分电荷和 Mulliken 电荷经常被放在同一段里，但它们解决的不是同一个问题。Bader 把电子密度按空间分区积分，差分电荷看吸附或成键前后的电子重排，Mulliken 则把电子按基组和轨道布居分配。三者可以互相旁证，不能把绝对数值混在一起硬比。本文重点讲清楚：什么时候看净得失，什么时候看空间重排，什么时候需要把电荷证据收住。

先把“电荷”这个词拆开看？

Bader电荷到底在积分什么？

Bader分析不是直接告诉某个原子“真实带几价”，而是按照电子密度的零通量面把空间划成不同区域，再对每个区域内的电子数积分。它适合比较同一方法、同一设置下某个原子或吸附物的净得失趋势。

在 VASP 里，Bader 通常基于 CHGCAR 或 AECCAR 组合后的电荷密度来做。真正需要记录的是使用了哪套电荷密度、赝势价电子数怎样设定、总电荷是否守恒，而不是只把一个 +0.30 e 或 -0.20 e 摘出来。

式(1) Δρ = ρ_total − ρ_A − ρ_B

差分电荷为什么不是一个原子电荷？

差分电荷看的不是某个原子的净电荷，而是电子密度在空间中哪里累积、哪里耗尽。图1中 Bader 与差分电荷并排出现时，前者更像定量旁证，后者更像空间路径：电子从哪一侧重排到哪一侧。

因此，差分电荷图不能直接读成“某原子得到 0.3 e”。黄色和蓝色区域说明的是相对于参考片段的密度变化，参考片段怎么拆、坐标是否一致、网格是否统一，都会影响结果。

比较稳的做法是把Bader 表格和差分电荷图配套使用：先用差分电荷说明重排方向，再用 Bader 给出净得失趋势，最后用吸附能、键长或 PDOS 判断这种转移是否真的和活性相关。

图1 Bader 电荷和差分电荷密度共同分析 CO₂RR 活性位。DOI: 10.1038/s41467-025-56534-1。

Mulliken电荷在看什么？

Mulliken电荷不是按实空间区域切分电子密度，而是把电子按基组函数和轨道布居分到不同原子上。它更像是在回答：在当前基组和布居规则下，哪些轨道把电子贡献到了哪个原子或片段。

放到 VASP 项目里看，Mulliken 通常来自局域基组、投影或额外后处理，而不是直接从 CHGCAR 做空间积分。因此它适合在同一程序、同一基组、同一投影设置下比较趋势，不适合拿绝对数值去和 Bader 电荷硬对。

为什么三种数值不能混在一张表里比较？

Mulliken 为什么换基组就会变？

Mulliken 电荷依赖基组重叠和轨道布居划分，换一套基组或换一种投影方式，数值就可能明显变化。它适合在同一计算框架内部看趋势，不适合拿来和 Bader 的绝对数值一一对应。

如果要比较不同吸附构型，先保证方法不变、基组不变、模型边界不变。只要分区规则变了，电荷数值的物理口径就变了。

式(2) q_Bader = Z_val − N_Bader

同一体系为什么会给出不同电荷数？

图2里的差分电荷等值面适合说明电子重排位置，而不是给出原子电荷。Bader 的分区边界、Mulliken 的基组布居、差分电荷的参考片段各不相同，所以同一体系得到不同数值并不奇怪。

真正该问的不是“哪个数值才是真的”，而是“这个方法回答的问题是否和结论匹配”。讨论净电荷转移，用 Bader 更直接；讨论空间极化和界面重排，用差分电荷更直观；讨论分子轨道布居，Mulliken 才有语境。

如果两个方法给出趋势相反，不要急着选一个支持自己结论。先检查参考态、网格、价电子数、基组和吸附构型，再判断这种差异是不是方法边界造成的。

图2 界面接触前后的空间电荷密度差分与线平均电荷变化。DOI: 10.1021/acsomega.3c05213。

后处理文件和参数该怎么留？

CHGCAR 相减最怕哪里出错？

差分电荷相减最怕三件事：FFT 网格不一致，原子坐标没有对齐，参考片段和总体系不是同一几何。只要其中一个条件不满足，Δρ 图里出现的“累积/耗尽”就可能是人为差值。

实操时建议把 total、slab、adsorbate 三个目录分开保存，并在相减前确认晶胞、坐标顺序、NGXF/NGYF/NGZF 或输出网格一致。相减脚本、等值面数值和颜色范围也要一起保留。

Mulliken 结果怎么留记录？

式(3) q_Mulliken = Z − P_μνS_μν

图3里的 Mulliken 电荷分布更适合放在轨道布居语境下读。它的重点不是和 Bader 数值逐项对表，而是说明在当前基组和布居规则下，电子被分配到哪些原子或片段上。

记录 Mulliken 结果时，至少要写清楚程序、基组或投影方式、布居分析方法、自旋通道以及被比较的原子/轨道范围。只给一个 q 值而不写这些条件，后面很难判断差异来自体系本身，还是来自布居规则。

如果同一组数据里同时出现 Bader 和 Mulliken，建议分成两列说明：Bader 用来讨论实空间净得失，Mulliken 用来讨论轨道布居趋势。两者可以互相参考，但不要把绝对电荷数值混成同一把尺子。

图3 Mulliken 电荷分布用于分子体系电子结构分析。DOI: 10.1021/acsomega.5c08549。

Bader 结果要检查哪些细节？

Bader 结果至少要检查总电荷守恒、价电子数设置、网格精度，以及是否需要使用 AECCAR 修正。对于含真空层、带电体系或强极化界面，还要留意电荷密度是否已经足够收敛。

不要把一张差分电荷图或一组 Bader 电荷当成最终机制证据。它们适合提出电子重排和净得失线索，但机制解释还需要吸附能、键长、PDOS、COHP 或反应自由能来支撑。

写机制时怎样把电荷证据放稳？

哪些话可以写，哪些话要收住？

比较稳的写法是：“Bader 分析显示吸附物净得电子增加”“差分电荷显示界面处电子累积”“Mulliken 布居提示某轨道电子占据变化”。这些说法都把方法边界交代清楚，没有把电荷数值写成真实价态。

需要收住的是绝对化表达。比如“该原子带 +1 价”“电子全部转移到吸附物”“电荷转移越多活性越高”，这些说法通常都需要更多证据，不能只靠电荷分析支撑。

式(4) Σq_i = Q_total

什么时候必须补吸附能、PDOS 或 COHP？

图4展示的是更完整的证据链：电荷极化可以解释位点附近电子环境变化，但吸附能和 PDOS 才能进一步说明这种变化是否有利于中间体吸附或轨道耦合。

如果结论只写“界面发生电荷重排”，差分电荷和 Bader 基本可以支撑；如果结论写到“活性提高”“吸附增强”“反应路径更优”，就必须补能量或轨道证据。

尤其要避免把“电荷转移多”直接等同于“催化更好”。有些体系中电子转移增加会增强关键吸附，也可能导致中间体过强吸附，反而拖慢后续步骤。

写论文时可以按这个顺序组织：先交代电荷方法和参考态，再描述净得失或空间重排，随后用结构、吸附能或电子态解释这种变化为什么与反应相关。

图4 电荷极化位点中吸附能、PDOS 与差分电荷的 DFT 分析。DOI: 10.1038/s41467-026-72465-x。

本章要点总结

要点1：Bader电荷是实空间分区积分，差分电荷是空间重排图，Mulliken 电荷是基组布居分析。

要点2：不同电荷方法的绝对数值不能混比，只能在同一方法、同一设置下比较趋势。

要点3：差分电荷适合看电子从哪里重排到哪里，Bader 适合给净得失的定量旁证。

要点4：电荷证据不能替代吸附能、键长、PDOS 或 COHP，机制解释要把证据链补完整。