在计算化学和材料科学中，Hirshfeld 电荷、Mulliken 电荷和 Bader 电荷是三种常用的电荷分析方法，用于量化分子或固体中原子的电荷分布。它们的理论基础、计算方法和适用场景各有差异。

一、Mulliken 电荷（Mulliken Charge）

理论基础

        基于分子轨道理论（MO 理论），通过原子轨道在分子轨道中的贡献来划分电荷。
        核心假设：分子轨道由原子轨道线性组合而成（LCAO-MO），电荷按原子轨道在分子轨道中的系数平方分配。

优缺点

• 优点：
  •         计算简单，广泛应用于早期量子化学计算（如 Gaussian、MS等软件默认支持）。
  •         适用于定性分析电荷分布趋势（如判断原子的正负电性）。
• 缺点：
  •         严格依赖基组选择，结果可能因基组大小或极化函数的引入而显著变化。
  •         假设成键电子完全平均分配，忽略电负性差异，对极性键或离子键体系误差较大（如 NaCl 中 Na 的 Mulliken 电荷可能远小于 +1）。
  •         对离域体系（如共轭分子、金属固体）的电荷划分不合理。

典型应用

•         小分子体系的电荷定性分析（如判断有机分子中原子的亲电 / 亲核活性）。
•         对比不同结构或取代基的电荷变化趋势。




二、Hirshfeld 电荷（Hirshfeld Charge）




理论基础

        基于独立原子模型（Independent Atom Model, IAM），通过将分子中的电子密度与自由原子的电子密度对比，划分原子间的电荷。
        核心思想：每个原子在分子中的电荷等于其自由原子的电子数减去分子中该原子区域内的电子数。

优缺点

• 优点：
  •         无需人为假设成键电子分配，基于实际电子密度划分，比 Mulliken 更合理。
  •         对极性体系（如离子化合物）和共价体系均适用，结果更符合化学直觉（如 NaCl 中 Na 的 Hirshfeld 电荷接近 +1）。
  •         对基组敏感性低于 Mulliken 电荷。
• 缺点：
  •         依赖自由原子的电子密度模型，可能高估原子的局域性，对强共价键或离域体系（如金属、共轭分子）的电荷划分仍存在争议。
  •         计算复杂度高于 Mulliken 电荷，需先计算自由原子的电子密度。

  典型应用

  •         有机分子、主族化合物的电荷定量分析（如氨基酸、金属配合物）。
  •         对比不同力场或计算方法的电荷一致性。
  • 
    

  三、Bader 电荷（Bader Charge）

• 
  

• 理论基础

          基于量子化学中的原子在分子（Atoms in Molecules, AIM）理论，由 Richard Bader 提出，通过电子密度的拓扑分析划分原子区域。
          核心概念：以电子密度梯度的临界点（如键临界点、核临界点）定义原子间的边界，每个原子区域为Voronoi 型分区，满足电荷守恒。

• 优缺点

  • 优点：
    •         严格基于电子密度的物理拓扑，无需人为假设，是目前最严谨的电荷划分方法之一。
    •         适用于强相互作用体系（如离子键、金属键、氢键）和离域体系（如半导体、金属固体）。
    •         结果不依赖基组（仅依赖电子密度的收敛性），具有较好的可转移性。
  • 缺点：
    •         计算复杂度高，需额外进行 AIM 分析（常用软件如 AIMAll、TopMod）。
    •         对弱相互作用体系（如范德华复合物）的边界划分可能存在歧义。
    • 
      

  典型应用

  •         固体材料（如氧化物、金属合金）的电荷分析。
  •         强极性分子（如酸碱复合物）、过渡金属配合物的电荷定量研究。
  •         化学反应机理中的电荷转移分析（如催化反应、电子转移过程）。
  • 
    

  四、三种电荷的对比总结

  特征	Mulliken 电荷	Hirshfeld 电荷	Bader 电荷
  理论基础	分子轨道理论（均摊法）	独立原子模型（电子密度权重）	AIM 理论（电子密度拓扑）
  电子密度依赖性	间接（通过轨道系数）	直接（对比自由原子密度）	直接（基于实际密度分布）
  基组敏感性	高	中	低
  适用体系	小分子、定性分析	有机分子、极性化合物	强相互作用体系、固体
  结果合理性	定性趋势合理，定量误差大	定量更合理（尤其极性键）	最接近实际电荷分布
  计算复杂度	低（软件默认支持）	中（需自由原子密度）	高（需拓扑分析）

  
  

  • 五、选择建议

  •         若只需定性对比电荷趋势或计算资源有限，选择 Mulliken 电荷。
  •         若需定量分析有机 / 主族化合物的电荷，且避免均摊法的缺陷，选择 Hirshfeld 电荷。
  •         若研究强相互作用体系、固体材料或精确电荷转移，选择 Bader 电荷（需配合 AIM 分析工具）。
  • 
    

  实际应用中，常结合多种方法交叉验证，例如同时计算 Mulliken 和 Bader 电荷，对比分析体系的电荷离域程度或计算方法的可靠性。