差分电荷密度图怎么看？从空间形貌、界面偶极到轨道成键与实验物性的完整对应

说明：本文华算科技主要介绍差分电荷密度的计算对象、空间形貌、界面偶极及其与轨道成键和实验物性的对应范围。

差分电荷密度是什么？

差分电荷密度通常写成 Δρ(r) = ρ_AB(r) − ρ_A(r) − ρ_B(r)。ρ_AB来自组合体系，ρ_A与ρ_B来自两个片段。三个密度应采用同一晶胞、同一原子坐标、同一赝势与网格，片段计算仅移去另一部分原子。若片段另行弛豫，Δρ中会叠加几何改变造成的密度差，界面成键与构型变化便难以区分。

Δρ大于零的位置表示组合后电子密度增加，小于零的位置表示电子密度减少。红、黄、蓝、青都只是绘图配色，正负号必须依据图例或原文定义。单位常见 e Å⁻³、e Bohr⁻³，二维平面平均后变为 e Å⁻¹或按面积归一的写法。跨论文比较颜色面积之前，等值面阈值和单位必须一致。

原子球与键线给出空间坐标，彩色等值面标出 |Δρ| 达到指定阈值的区域。阈值调低会出现更多稀疏区域，调高后只剩峰值附近的密度重排。等值面体积受阈值控制，颜色块变大并未自动对应更多转移电子。切片保留平面内的连续数值，更适合辨认键轴、孤对电子方向和界面法向。

PEDOT/Te 模型把三种空间信息排在同一幅图中：界面附近红蓝区域反映跨界面重排，PEDOT 链内部的红蓝交替对应 π 轨道向 σ 键区域的密度调整，右侧态密度用于检查组合前后的电子态变化。原文采用 0.005 e Å⁻³等值面，Te 与第一层 PEDOT 间距约 3.6–4.0 Å，跨界面转移较弱，却存在明显的链内 π/σ 密度再分配。大面积彩色区域由此可能主要来自极化或泡利排斥。

组合体系与中性片段含有相同电子总数时，全空间积分 ∫Δρ(r)dr理论上接近零。局部积累与局部耗尽仍可形成化学键、界面偶极或分子内部极化。零总积分约束的是整个晶胞，某个片段获得多少电子则依赖空间分区、积分面或原子盆划分。

积累区/耗尽区的形貌对应哪些电子重排

相邻原子之间出现沿键轴连续的积累区，并伴随原子外侧耗尽，常与成键电子共享相符；积累区偏向电负性较高的一端，则带有极化键特征。键间积累的位置与方向比单个颜色块的峰值更有物理意义。若红蓝区域主要包围同一原子或同一分子，密度变化可能来自轨道形状调整，跨片段净转移量可以很小。

金属表面与分子、二维层或电极接触后，表面电子云重叠产生泡利排斥。靠近金属最外层的耗尽薄层与其外侧积累薄层构成“push-back”图案，方向与功函数差驱动的长程转移可能相反。一对相邻正负薄层对应局域极化；积分平台的稳定值给出划定区域内的净电子数。

In/MoS₂与 Au/MoS₂采用相同的 1×10⁻³ e Å⁻³等值面，两组界面均在金属表面附近出现正负振荡。累计曲线在 MoS₂一侧的平台却给出相反符号：In 向 MoS₂转移约 5.0×10¹² e cm⁻²，Au 接触使 MoS₂损失约 1.7×10¹³ e cm⁻²。局域峰形和跨界面净转移属于两个量，功函数与电荷中性能级决定后者的方向。

层状异质结构中，正负区域在法向上分离便产生界面偶极，即使两侧各自含有多组抵消峰。偶极矩面密度与 ∫zΔρ(z)dz相关，原点平移在总电荷守恒时不改变结果。峰的空间间距参与偶极矩，两个体系即使拥有近似峰高，也可能因正负峰间距不同而产生不同电势阶跃。

法向正负峰怎样转成界面偶极/物性线索

三维密度沿界面平面平均可写为 Δρ̄(z) = A⁻¹∫∫Δρ(x,y,z)dxdy。横坐标必须标出原子层位置，正负峰才能归属到吸附物、表面层或次表层。平面平均保留法向电荷守恒与偶极，面内的配位差异、桥位与顶位差异会被平均掉。

覆盖度变化会同时改变吸附物数量、侧向相互作用和表面弛豫。若用峰值绝对值衡量“电荷转移强弱”，峰宽和峰位移动均被遗漏。单位面积累计电荷由曲线面积决定，表面偶极还包含位置权重。峰高、峰面积与一阶矩分别回答不同物理量，三者之间不存在通用的单调替代关系。

O/Fe/Ir(001)随氧覆盖度从 0.25 ML 增至 1 ML，O、Fe 与 Ir 原子层附近的正负峰位置和幅度均发生变化，表面偶极与 DMI 呈非单调响应。0.5 ML 和 0.75 ML 的实空间等值面还显示不同吸附位点周围的密度重排。覆盖度序列中的共变关系支持表面电荷不对称与 DMI 调制相联系，定量 DMI仍由自旋轨道耦合下的自旋螺旋能量给出。

界面偶极会改变真空能级与平均静电势阶跃，可进一步投射到功函数、肖特基势垒和能带对齐。磁性异质结还包含重元素自旋轨道耦合、交换作用和轨道杂化。差分电荷曲线用于定位界面电场的来源，物性数值来自相应哈密顿量与能量差，例如 DMI 由顺、逆手性磁构型的能量差或自旋螺旋色散提取。

CoFe/W 模型把界面增密峰 Peak1、减密峰 Peak2 与 W 厚度逐一对应。Peak1与计算 DMI 都出现起伏，但转折位置并非完全重合；Peak2的变化又具有另一种趋势。同一参数序列中的多条曲线可以排除“某个峰越高，DMI必定越大”的单指标解释，并把讨论限定在该晶格、层数和磁构型下。

轨道成键与实验量怎样检验空间图案

差分电荷密度把全部占据态对实空间密度的贡献相加，图上没有能量标签，也没有轨道标签。键轴附近的积累区可提出成键假设，PDOS用于确定参与杂化的原子轨道和能区，COHP或 COOP则区分成键态与反键态；键长、振动频率和吸附能提供结构及能量响应。

同一体系若在费米能级以下出现轨道峰重叠，仍可能同时含有成键和反键贡献。−pCOHP正值常按所用约定表示成键贡献，积分 COHP需注明能量上限。轨道重叠与键强度由分能量分析约束，差分电荷图负责给出密度重排的空间位置，两类图谱的变量不可互换。

NHC/Cu(100)模型中，Cu 吸附原子与卡宾 C 之间出现电子积累，Cu 周围出现耗尽；Bader 分区给出 Cu 吸附原子约 +0.179 e。过渡态 PDOS 中 Cu d_xy与 C p_x重叠增强，−pCOHP 和 ICOHP 显示吸附原子改变了 CO 键的反键态占据，C—O键由 1.162 Å延长至 1.175 Å。空间密度、轨道能区、键长与 0.440 eV 能垒共同限定了 NHC 修饰后的 CO 氢化机制。

界面体系还可用功函数、XPS/UPS能级、开尔文探针、势垒高度、载流子浓度和器件电流检验计算趋势。不同观测量对平均电势、局域化学环境或输运过程的响应各异。实验信号应对应差分电荷提出的具体后果：偶极对应真空能级位移，局部氧化对应芯能级变化，能带弯曲对应空间电势与载流子分布。

PVK/钙钛矿和 PPETA/钙钛矿模型中的黄色、青色区域分别对应电子积累与耗尽，PPETA界面的法向 Δρ峰更强且分布更宽。UPS测得功函数由 4.37 eV移至 4.78 eV，电离势由 5.97 eV移至 6.02 eV；C–V峰位和空穴单载流子 J–V 曲线随界面修饰同步变化。真空能级、载流子注入与界面密度重排由这些独立量分别检验，器件电流仍包含迁移率、缺陷俘获和接触电阻的贡献。