第一性原理计算说明：本文主要介绍差分电荷密度图在第一性原理计算中的定义来源、参考态选择、等值面读法、剖面图读法、积分曲线和适用范围。
差分电荷密度比较的是哪两个电子密度
差分电荷密度常写成 Δρ(r)，它不是某个原子自带的电荷数，而是两个电子密度函数相减后的空间分布。吸附模型里常用 ρsurface+ads(r) – ρsurface(r) – ρads(r)，界面模型里常用 ρA/B(r) – ρA(r) – ρB(r)。参考态一换，Δρ(r) 的物理含义也会改变：固定几何的片段相减，强调成键后的电子重排；用弛豫后的独立片段相减，结构松弛贡献已经进入差分结果。
读差分图时，第一件事不是盯着颜色，而是确认参与相减的电子密度来自同一套超胞、同一网格、同一泛函和同一赝势设置。周期性体系里的 ρ(r) 是连续场，黄色、蓝色或红色等值面只是把某个 Δρ 数值切出来显示。颜色代表电子密度增减的位置，不自动等同于整数电荷从一个原子跑到另一个原子。若论文没有写正负号颜色定义，正文里就不能把某种颜色直接写成电子积累或耗尽。
Mg 体相案例把参考态问题放得很直观：无 H 的晶格里，间隙区域已经有由金属电子形成的电荷积累结构；加入 H 后，高对称间隙周围出现新的局域响应。这类图适合说明 Δρ(r) 的起点：差分图读的是体系变化后的密度响应，不是孤立原子的价电子清单。如果把无 H 和含 H 的两幅图当成普通结构图，容易忽略等值面水平、片段参考和网格对齐。
差分电荷密度还要和总电荷密度分开。总电荷密度告诉我们电子云主要在哪里，通常在原子核附近数值很大；差分电荷密度把相互作用前后的差别放大，常在键区、界面、缺陷邻域或吸附位点附近更醒目。所以它适合讨论成键诱导的极化、界面偶极、缺陷屏蔽和吸附物活化，不适合单独承担价态、反应速率或催化活性的完整判断。后续结论应回到同一结构模型里的其他计算量。
等值面/剖面先读正负号，再读空间位置
三维等值面读起来很像分子轨道云图，但它的数值含义不同。等值面阈值设得低，远处微弱振荡会被放大；阈值设得高，只剩最强的局域积累或耗尽区域。同一组结构若使用不同阈值，等值面体积就不能直接比较。可比较的是同一色标、同一等值面水平、同一参考态下的相对空间分布。换成不同色标时，图形大小会跟着可视化设定改变。
Mn-H 案例适合练习读图顺序。先定中心原子和缺陷位点，再读红色积累、蓝色耗尽与周围 Mn 原子的相对位置；接着看灰色截面上的正负等值线是否穿过 H-Mn 键区、间隙中心或远离成键方向的区域。右侧 RDF 曲线把空间扰动压成以 H 为中心的径向信号，峰高反映局域积累强弱，峰宽反映扰动延伸范围。这些曲线信息需要和左侧结构位置互相参照。
剖面图的优点是能显示正负值在同一平面内的连续变化。三维图上看似相连的等值面，在剖面中可能被分成几个局域峰；三维图上不醒目的耗尽区，剖面色标可能给出清晰负值区域。差分图的空间位置必须和原子坐标一起读，否则很容易把表面真空区的微弱振荡、周期镜像附近的响应或可视化阈值造成的形状变化写成成键证据。
Si 的补充案例提醒我们，差分电荷密度不只服务于金属表面吸附。H 位于 HCP Si 超胞时，等值面沿 Si 原子方向出现拉长，剖面图里正值峰集中在局部区域。这类信号可以支持“缺陷周围电子密度发生局域重排”的说法；若要继续写成 Si-H 成键增强，需要键长、PDOS 或 COHP 等电子结构证据同步约束。
Hg 体系里的负值区域更显眼，电子积累和耗尽在剖面上同时出现。这时若只截取正等值面，读者会误以为 H 周围主要是积累；若只看二维色块，又可能忽略三维周期结构中的重复响应。正负号、色标范围和剖面位置需要同时写进读图说明。否则同一张图会被解释得过窄。
积分曲线把形状变成数值，但会丢掉方向信息
很多文章会把差分电荷密度沿 z 方向积分，得到平面平均 Δρ(z)，再对 z 积分得到累积电荷 ΔQ(z)。界面、电极和二维材料里，这种曲线能定位偶极层、屏蔽距离和电荷重排的主区域。体相缺陷里，若扰动近似围绕缺陷中心展开，也可以转成径向分布函数 RDF，用峰位、峰高和峰宽比较不同金属基体对同一 H 缺陷的电子响应。三者分别对应扰动位置、强弱和延伸范围。
RDF 曲线的好处是把一大组三维电荷扰动压成可比较的峰。Cr、V、Li、Ti、Mn 等曲线在峰高和峰宽上不同，说明同一个 H 缺陷在不同 FCC 金属里诱导的屏蔽响应并不相同。但曲线化的代价也很明确：径向平均会弱化方向差异，某些沿特定键向延伸的电子积累，在 RDF 中可能只表现为峰宽变化。
把每条 RDF 除以自身最大值，再减去平均曲线，得到的是相对平均响应的偏离。这种处理不再追求“哪个点颜色更深”，而是比较不同元素基体的峰位、极值数量和振荡范围。贵金属样品在曲线振荡上更复杂，早期过渡金属则更接近另一类响应。差分电荷密度图可以变成描述符，前提是归一化方式、中心点、网格尺寸和构型类别写得足够明确。
平面平均或 RDF 积分不能替代原图。曲线告诉我们“沿某个坐标或半径的净变化”，原图告诉我们“变化发生在结构里的哪个位置”。界面体系若只给 ΔQ(z)，读者仍需要知道界面法向、真空层、哪一侧是 A 层、哪一侧是 B 层；缺陷体系若只给 RDF，读者仍需要知道 H 坐在八面体间隙还是四面体间隙。
结论要回到模型、分区电荷和能量响应
差分电荷密度最常见的误读，是把可视化图像直接写成电荷转移量。等值面体积不是电子数，颜色深浅也不是原子净电荷。要给出数值，通常需要 Bader、Hirshfeld、DDEC 或片段积分；要讨论成键，还要连接键长、PDOS、COHP、ELF 或局域态分布。差分图负责定位电子密度在哪里重排，分区电荷负责给出选定原子盆或片段中的积分数值。两类结果的单位和积分区域不同。
线性关系图给了一个很好的限制条件：RDF 峰值与 H 的 Bader 电荷正相关，与金属电荷半径呈负相关倾向。这说明差分电荷密度构造出的曲线与分区电荷之间存在可比较的数值联系。但这种关系来自同一批 FCC 金属、同一 H 间隙缺陷和同一套计算设置，不能直接搬到氧化物表面、二维异质结或带电缺陷模型里当成通用数值关系。适用体系改变时，相关性需要重新验证。
相关矩阵把 Bader 电荷、Bader 半径、RDF 峰宽、RDF 最大值与晶格常数、体模量、原子密度和 H 结合能排在同一张图里。这里的重点不是把某个相关系数写成因果，而是看哪些电子密度特征与结构尺度、压缩响应或结合能变化同步。差分电荷密度进入材料筛选时，最有价值的部分正在这里：它把局域电子响应转成能与能量、结构和元素性质比较的数值特征。这种转化适合筛选线索，机制判断仍需能量路径和电子态证据继续约束。
Al 和 Mn 的 NEB 能垒差异说明，电子密度特征可以提示缺陷迁移环境的软硬差别，真正的动力学结论仍然来自能量路径。Al 中 H 的 Bader 电荷体积更大，对应的迁移能垒较低；Mn 中 H 的电荷体积更小，能垒更高。这组结果适合写成“电荷空间延展与迁移能垒存在对应”，不适合写成“差分电荷密度图已经给出扩散速率”。
实际读图时，可以按一个固定顺序检查：参考态是否明确，正负号和色标是否明确，等值面阈值是否统一，剖面是否穿过关键原子或键区，平面平均或 RDF 是否说明积分方向，Bader 或片段积分是否来自同一结构，能量和态密度证据是否支持相同物理图像。差分电荷密度图真正能回答的问题，是电子密度在相互作用后如何重新分布，以及这种重排在空间上靠近哪个结构单元。
