![]()
![]()
![]()
说明：本文华算科技主要介绍差分电荷密度在材料计算中的基本含义、图谱读法、物理化学意义和判断条件。
差分电荷密度通常写作 charge density difference，核心操作是在同一空间网格上比较组合体系电子密度与参考片段电子密度之和。以表面吸附为例，常见表达是 Δρ = ρsurface+adsorbate – ρsurface – ρadsorbate。这里的关键不是公式长什么样，而是三个条件：同一构型、同一网格、同一参考态。只有这三点明确，图中的正负区域才主要对应吸附或界面接触诱导的电子重新分布。
它首先回答的不是“某个原子得到了多少电子”，而是电子密度在空间中哪里增加、哪里减少。在金属表面吸附、单原子催化剂、二维异质结、缺陷掺杂和电池界面中，差分电荷密度可以帮助定位电荷重排发生在金属中心、配位原子、吸附键、缺陷位点还是层间界面附近。这个定位能力很重要，但它仍然是实空间图像，不是电荷布居数值。
结构模型是读差分电荷密度的起点。若没有先明确表面、吸附物、配位中心或异质界面的几何关系，后面的颜色等值面就很容易变成孤立图案。对于 Mo-Tren 这类配位体系，Mo 原子和 N 原子的空间位置决定了电子重排应该优先在哪里检查；对于 slab 吸附模型，吸附物高度、吸附位点、表面晶面和覆盖度都会影响差分图的解释范围。
差分电荷密度、Bader 电荷、Mulliken 电荷和 Hirshfeld 电荷必须分开。前者看空间区域，后者按分区规则给出数值估计。颜色区域更大、等值面更显眼，并非某个原子转移电荷更多；若要比较电荷转移量，需要另外说明采用哪一种电荷布居方法。
差分电荷密度最常被用来讨论电子积累和电子耗尽。一般来说，电子积累区域靠近吸附物和表面之间、金属中心与配体之间或层间界面附近，往往提示这些位置出现了电子云重排。电子耗尽区域则说明相对于参考片段，该空间范围内电子密度降低。这样的图像可以支持界面极化、吸附键形成、局域电子响应、缺陷诱导电荷重排等判断。
不过，电子积累不自动等于“成键更强”。如果积累区域位于两个原子之间，它可以作为相互作用位置的证据；若要讨论成键强弱，还要结合键长、吸附能、PDOS、COHP 或反应路径。对于催化体系，差分电荷密度更适合说明某个中间体在哪个方向被极化，而不是直接说明反应一定更快。
在吸附体系中，电子密度重排通常要和具体中间体绑定。N2O、CO、O2 这类小分子吸附后，差分图可以显示金属-吸附物键附近是否有积累区，也可以提示分子内部某个键是否被极化。这里的判断对象是吸附构型中的局部空间区域，而不是整块材料的全部电子性质。
对于含氧中间体，差分电荷密度还常被用来解释键长变化。若 O-O 键附近出现明显电子重新分布，同时优化结构中 O-O 键被拉长，这可以支持“吸附导致分子内部键被极化或活化”的判断。这里的证据适用条件是：键长变化说明几何响应，EDD 说明电子密度响应，二者互相支撑但不等同。若要进一步判断断键是否容易发生，还需要反应路径和能垒信息。
空间位置是差分电荷密度的优势。它能告诉读者重排发生在 Mo-C、Mo-O、N-Mo 或层间界面附近；但它不能单独区分成键态和反键态，也不能给出反应热力学或动力学结果。若一篇文章需要讨论催化活性，差分电荷密度只能承担“电子重排位置”这一类证据。
读图的第一步是确认颜色定义。不同论文可能用黄色表示电子积累、青色表示电子耗尽，也可能使用红蓝配色，甚至把正负颜色反过来。颜色必须服从原图图注或作者说明，不能凭经验直接解释。若图注没有说明正负颜色，正文里就不宜写得过满，只能把图作为局域电子重排的辅助线索。
第二步是确认等值面阈值。差分电荷密度图常设置 isovalue：阈值低时，等值面范围会变大，细碎区域增多；阈值高时，只剩下较强重排区域。两个体系若使用不同阈值，图面大小不能直接比较强弱。比较不同材料或不同吸附构型时，相同阈值、相同视角、相同参考态比颜色面积本身更关键。
第三步是确认参考态。吸附体系通常用整体吸附态减去表面片段和吸附物片段；异质结常用界面整体减去两个单独层；掺杂或缺陷体系则要说明是和完美晶体比较，还是和替换前后的片段比较。参考态一变，正负区域的含义就会跟着变化。
还要检查参考片段是否保持吸附态中的几何位置。若表面片段和吸附物片段重新优化后再相减，差值中会混入结构松弛带来的密度变化；若保持吸附态构型，图谱更接近“接触前后电子云响应”的比较。构型一致性会直接影响颜色区域能否被解释为界面相互作用，也会影响不同吸附构型之间的可比性。
第四步是检查它和其他图谱的分工。AIMD 轨迹回答有限温度和有限时间内结构是否保持，声子谱回答振动模式是否存在虚频，PDOS 回答能量轴上的轨道贡献，差分电荷密度回答实空间电子重排。同一篇文章里出现多类图谱，并非说明它们能互相替代，每种图只覆盖材料计算问题的一部分。
第一个误判是把差分电荷密度当成电荷转移量。图中的黄色或青色区域只能说明相对于参考态的电子密度变化位置，不能直接读出“某原子得到 0.2 e”。若要讨论转移数值，应使用 Bader、Hirshfeld 或 Mulliken 等布居分析，并说明不同方法的分区规则不同。Bader 电荷变化也不可当作形式价态变化，价态还要结合配位、磁矩和谱学证据。
这种限制在缺陷、掺杂和强关联体系中尤其重要。过渡金属 d 电子、氧空位局域态、表面自旋态和 U 值设置都会改变局域电子密度分布。若文章只给一张差分图，却没有说明自旋极化、U 值、泛函、真空层或覆盖度，读者很难判断颜色区域来自真实电子响应，还是来自模型设置差异。计算条件不是附属信息，它决定差分图能承担多强的解释。
第二个误判是把差分电荷密度直接写成反应更容易。反应是否容易发生，至少涉及初末态能量、过渡态或 NEB 能垒、零点能和熵修正、温度、电位、覆盖度和溶剂环境。差分电荷密度可以解释某个中间体为什么被极化，或某个键附近为什么出现电子重排，却不能替代反应路径计算。
第三个误判是把差分图替代 PDOS。PDOS 的横轴是能量，能够显示费米能级附近哪些原子轨道贡献较大、吸附物轨道和金属 d 态是否存在能量重叠；差分电荷密度的坐标是空间位置，用来显示电子密度在哪里增减。二者经常共同支持吸附相互作用分析，但回答的问题并不相同。
第四个误判是忽略模型条件。对于 slab 模型，真空层厚度、偶极修正、表面晶面、覆盖度和超胞尺寸会影响界面电荷分布；对于过渡金属体系，自旋极化、U 值和 vdW 修正可能改变吸附能、局域磁矩和态密度。写差分电荷密度时，判断范围应限定在具体模型和计算条件内，例如某个 O* 中间体、某个 Mo-N 配位环境、某个层间界面或某个缺陷位点的电子重排。
因此，差分电荷密度适合停在这些判断上：电子密度是否在吸附键附近积累，缺陷位点周围是否出现局域重排，界面两侧是否形成极化区域，某个吸附构型是否伴随分子内部键的电子密度改变。进一步讨论电荷转移量、轨道杂化、成键稳定化、反应能垒或材料稳定性时，应分别回到Bader/Hirshfeld 电荷、PDOS/COHP、吸附能、NEB、AIMD 或声子谱这些对应证据上。
