说明:本文华算科技主要介绍差分电荷密度的数学定义、参考密度选择、常见图形来源,以及它在成键、界面极化和缺陷分析中的使用条件。


差分电荷密度基础定义/计算来源
密度泛函计算给出基态电子数密度 n(r),它记录单位体积内的电子数随空间坐标的分布。差分电荷密度是在同一位置对两组或多组密度做减法:Δn(r) = ntarget(r) − nref(r)。减法保留了三维空间信息,单位通常写作 e Å−3 或 e bohr−3;正负值分别对应目标态相对参考态的电子数增加和减少。

许多论文把 n(r) 直接称作“电荷密度”,此时正的 Δn 表示电子积累。若程序输出或作者采用带电荷符号的 ρe(r) = −en(r),正负号的物理指向会翻转。颜色只能在图例、等值面数值和密度定义确定后解释,同一种暖色在不同论文里可能对应相反符号。
差分场还保留电子数账本。目标态与参考态包含相同电子数时,整个计算区域的积分满足 ∫Δn(r)dr = 0,局部积累必须由别处的耗尽补偿。带电缺陷、加减电子或组分发生变化时,体积分等于两套密度的电子数差;周期性带电模型使用的补偿背景也要写进解释条件。
图1把三种容易混淆的量放在同一体系中:总电子密度覆盖全部占据电子,净自旋密度取自旋向上与向下密度之差,形变电子密度则扣除了选定参考。差分电荷密度这个名称只规定了“做减法”,具体减去谁仍由研究对象和公式决定。


参考密度怎样改变减法的物理含义?
形变电荷密度常用自洽总密度减去中性自由原子密度的叠加,写成 Δndef = nSCF − ΣInIatom。它记录原子从彼此孤立到形成晶体后发生的轨道重排,成键区电子聚集、离子极化和孤对电子取向均计入同一差分场。赝势价电子密度与全电子密度在原子核附近差别很大,跨程序比较时必须核对所用密度类型。

研究吸附或异质界面时,参考项常换成冻结片段:Δnint = nAB − nAfrozen − nBfrozen。片段 A、B 的离子坐标沿用复合体系 AB,三项密度必须采用同一晶胞坐标与同一实空间网格。这样得到的是片段接触后新增的电子重排;若片段分别充分弛豫,几何形变引起的密度变化将混入结果。

外场、应变、缺陷浓度或磁态变化还可采用状态差:Δnstate = n2 − n1。两态共用离子构型时,差值偏向电子响应;各自弛豫后的密度相减则包含原子位移带来的变化。参考态写法决定差分场回答哪一个物理过程,论文中只有“charge density difference”而没有减法公式时,正负区域无法获得唯一解释。


三维函数怎转换成等值面/切片/曲线?
三维等值面选取阈值 ±c,画出满足 Δn(r) = ±c 的空间表面。较小的 c 会展开低幅度区域,较大的 c 只保留变化强的位置。等值面体积受阈值控制,不同材料、不同图之间若阈值或单位不一致,色块大小没有直接可比性。
二维切片直接保留指定平面上的 Δn(x,y,z0) 数值,适合检查键轴、层间空隙或缺陷邻域。切面未穿过相关轨道时,真实的三维重排会在图上显得很弱。切片位置属于图谱定义的一部分,晶面指数或截面坐标应与色标一同记录。

图4的 a–i 面板把同一类 Δρ 投影到 (110) 面、[X–Sb] 层的 (001) 面及 Sr 层的 (001) 面。X–Sb 键附近与 Sr 层周围的分布由各自切面呈现,后面的 XPS 和各向异性位移参数对应化学环境与晶格振动。截面差分只回答该平面上的空间分布,未经过的键或轨道仍要查询其他晶面或三维数据。
界面体系常沿平行方向积分,得到平面平均曲线 Δn̄(z) = A−1∫∫Δn(x,y,z)dxdy。继续对 z 积分可计算某段区域的累计电子数,对 zΔn 加权积分可联系界面偶极;若起始量是电子数密度,换算电荷偶极时应乘以 −e。

图5中的 Δρ(z) 把横向起伏压缩为沿界面法向的一条曲线,可定位积累与耗尽出现在哪一侧;H* 吸附等值面保留了位点周围的三维形状。两种画法来自不同空间操作,曲线峰高和等值面体积属于不同量,应回到各自的积分范围与阈值作比较。


差分场怎进入成键/极化/缺陷分析?
原子间区域出现电子积累时,可以检验成键电荷是否位于键轴或多中心配位区;积累偏向某一原子时,可以追踪极化方向。空间形状描述电子重排发生在哪里,电子转移数则来自指定区域积分、Bader 原子盆积分或其他明确分区。Bader 原子盆由总密度的零通量面划分,其数值不能由某个等值面的体积直接换算。

图6把邻界面硫族原子的平均 Bader 电荷与金属功函数列在同一材料组中。功函数差可驱动接触后的界面偶极,Bader 分区给出原子盆内的电子数,差分等值面定位积累和耗尽区域。三个量依次记录驱动力、原子盆电子数和重排位置;三类数据取值趋势相容时,界面极化解释可以返回原始数值复核。
成键强弱还涉及轨道能量和占据。差分场在两核之间积累,配套投影态密度出现能量重叠时,可讨论轨道杂化;COHP 把相互作用进一步分成成键与反键能区。局域电荷重排、轨道杂化和能量贡献回答的是三个相邻问题,它们的空间、能量和积分定义各不相同。

图7中,形变电荷密度定位 K、Na、O、Mn 邻域的电荷重排,部分电荷密度与投影态密度给出相关轨道的空间和能量分布,COHP 对代表性配位位点的键相互作用作能量分解。处理新的差分结果时,计算记录应包含参考密度公式、离子构型、电子数、密度单位、等值面阈值、切面或积分方向;缺少其中一项,对颜色范围或转移数的解释都会失去对应的计算条件。
