差分电荷密度图，等值面该怎么选？

说明：本文华算科技主要介绍差分电荷密度等值面的数学含义、阈值对空间形貌的影响，以及不同模型之间的比较条件和截面积分验证。

差分电荷密度的等值面是什么？

等值面把三维空间中具有相同电子密度差的点连成曲面，曲面穿过的位置由选定的数值决定。对于组分 A、B 接触产生的电子重排，常用 Δn(r) = n_AB(r) − n_A(r) − n_B(r) 表示。这里的 n 是电子数密度，三个密度在同一位置相减，正值记录相对参考态的积累，负值记录耗尽。

把显示阈值记为 c > 0，两种曲面分别满足Δn(r) = +c 与 Δn(r) = −c，常用单位为 e Å⁻³。在固定颜色的等值面图中，曲面上各点的密度差幅度相同，颜色只用于区分正负。曲面包围的区域可包含更大的幅度，它的外形随 c 改变，并非材料中预先存在的实体表皮。

单层、双层 1T-ZrO₂ 的原图左列绘制电子密度 n，采用约 0.60 和 0.45 e Å⁻³；右列绘制有正负值的电子密度差截面。n 与 Δn 是两种空间函数，总密度的等值面值不能照搬给密度差；ELF 的无量纲数值或单个轨道的等值面值也有各自定义。

密度差来自一次有明确参考的相减。冻结片段中的原子应保留整体结构的坐标，并采用一致的晶胞、网格和电子结构处理。若分别弛豫片段，原子位置改变本身就会贡献密度差。选 c 只能改变这些差值的显示范围，无法去掉错误参考或坐标错位产生的斑块。

电子数密度增加的符号与真实电荷密度的符号需分开：电子带负电，正的 Δn 表示电子增多。正负值与颜色的配对由图例规定，同一种黄色或蓝色在不同论文中可能表示相反含义。单位、相减顺序和正负色标共同确定这张图所画的量。

等值面阈值升高或降低会改变什么？

同一网格中，较低的 c 可绘出包围低幅度密度差的外围曲面，满足 Δn ≥ c 或 Δn ≤ −c 的空间范围随之扩大。提高 c 则只保留高幅度部分，外围壳层可能收缩，连成一片的区域可能分离。某个色团消失，说明该处未达到新阈值，原始密度差仍然存在。

AeF⁻（Ae 为 Be 至 Ba）的分子计算曾比较 0.001、0.005 和 0.01 e a₀⁻³ 等取值，a₀ 为玻尔半径。较低阈值与较高阈值显露的空间特征并不相同：较高取值能显出低阈值包络中不易辨认的键区积累。这些数值服务于该组分子的密度场，未构成其他材料的默认取值范围。

苯腈–铜离子计算中的等值面统一为 0.005 e Å⁻³，含两分子乙醇的络合物又采用两种视向展示。视向决定前后曲面如何遮挡，内侧特征可能被前景曲面盖住；两种金属模型电荷及有无乙醇则真正改变了计算对象。这幅图比较的是视向和模型，未改变显示阈值。

用于主图的 c 应让目标键区或界面两侧的正负重排可辨认，并在邻近取值下保持相同的位置和符号判断。若一个结论只在极窄取值内出现，应检查它究竟来自低幅度结构、包络相连还是视向遮挡。幅度差异很大的局域与弥散区域，可以分别补充另一阈值图和二维截面，每幅图标明自己的 c。

在分子附近电子密度梯度陡峭的位置，两个参考网格轻微错位就可产生相邻的正负残差。比较逐步加密的网格及充分自洽的密度差，可检查同一位置的幅度与符号是否稳定。低 c 会把残差和低幅度物理重排一起显示，增大 c 只能使部分区域不再显示，原始网格并未改变。

不同体系的差分电荷等值面怎样比较？

比较两个催化剂或吸附构型时，相同物理单位下的同一绝对阈值才能筛选相同幅度的密度差。0.001 e a₀⁻³ 与 0.001 e Å⁻³ 具有不同的体积单位，不能按数字相同比较。输入网格若含体积缩放，还应按文件约定恢复体积密度，避免把文件中的原始数值当作 e Å⁻³。

两种含锆前驱体片段在 ZrO₂ 上的密度差采用同一 5 × 10⁻³ e Å⁻³ 等值面值。这种对照固定了密度幅度的筛选条件，原图中键区位置与密度差形貌可以在同一标准下比较。相邻 PDOS 描述能量分辨的轨道贡献，图内以 e 标注的分区电荷数另由电荷分析得到。

有些研究为了显示幅度悬殊的响应，会给不同材料采用不同 c。Zn 吸附于 CuO 和石墨的原图分别使用0.005 与 0.001 e Å⁻³，并另列分区电子转移结果。石墨图中的大包络只表明低阈值下的空间范围，不能据此把它判为转移电子更多的体系。

逐图按各自峰值归一化可以比较重排形状，却会抹去绝对幅度差；按晶胞最大值选取相同比例，也可能受到某个局部强峰支配。若目标是比较界面整体响应，应固定绝对单位与显示值，并报告截面或积分结果。若目标是突出很弱的空间尾部，可以降低 c，但须把它与强弱排序分开。

晶胞横向重复会增加图中的色团数量，视向、裁切范围和透明度又会改变重叠程度。相同 c 还需配合相同的展示范围：比较每个吸附物附近的重排，就显示同等数量的吸附物和邻近位点；比较单位面积的响应，则明确表面面积与覆盖度。

差分电荷等值面该怎样用截面和积分核对？

三维等值面仅保留 Δn = ±c 的位置，二维截面可连续显示选定平面上的 Δn 数值。ZrO₂–MoS₂ 异质结的原图把正负等值面、取样平面和对应截面并列，三维部分取约 ±0.001 e Å⁻³。截面穿过哪排原子会改变局部图案，色标范围也要在对照中统一。

等值面包络的体积具有 ų 量纲，电子数变化来自对密度差在指定区域内积分：ΔN_Ω = ∫_Ω Δn(r) dV。将包络体积乘以 c，只把内部各点近似成同一个值，实际内部密度差并不恒定。若 Ω 随 c 改变，得到的也是随阈值改变的选区积分。

对于平面界面，可定义 Δn̄(z) = A⁻¹∫_AΔn(x,y,z) dxdy，再计算累计电子数变化 Q(z) = A∫_z₀^zΔn̄(z′) dz′。Q(z) 与显示阈值 c 无关，但依赖积分方向、起点和用于分隔片段的位置。平面平均还会抵消同一高度的正负贡献，局部键区的重排应由截面补充。

当整体与片段电子数满足 N_AB = N_A + N_B 时，全晶胞密度差积分应为零；带电模型则按实际电子数差检验。Bader 等分区方法通过指定原子域给出电子数估计，阈值扫描负责确认空间形貌。对某一键区或界面的电子数判断，应保留其积分区域、参考态和数值收敛结果。