电子局域函数与差分电荷分布:计算来源、读图要点及边界判据详解(elf电子局域函数) ypxx.net

说明:本文华算科技主要介绍电子局域函数和差分电荷分布的计算来源、读图对象和判断边界。

电子密度、参考态和局域性是什么?

第一性原理计算给出的核心对象是 电子密度 ρ(r),它在空间中连续分布,告诉我们电子更常停留在哪些区域。晶体、表面、缺陷和吸附体系里的电子不会天然按原子边界切开,也不会自动写出“谁向谁转移了多少电子”;总电荷密度先回答电子在哪里,后续图谱再把同一份密度加工成不同问题的答案。

图1. KPbBr3、RbPbBr3 和 CsPbBr3 在 (0 0 1) 面上的电荷密度分布。DOI:10.1038/s41598-025-11320-3

差分电荷分布通常写成 Δρ(r),含义是把相互作用后的体系密度减去若干 参考态密度。参考态可以是干净表面、孤立吸附物或两个分离的异质结组分;参考态一变,正负区域就会改变。差分图适合讨论吸附、掺杂、界面接触后电子积累与耗尽的空间位置,不能当作总电荷密度图或原子价态表。

与 Δρ(r) 的直接相减不同,电子局域函数 ELF 由 同自旋电子局域程度、动能密度和均匀电子气标尺构造,数值常在 0 到 1 之间。ELF 接近 1 常对应局域电子对、孤对电子或成键电子更集中的空间;接近 0.5 常接近金属态或均匀电子气特征。这里的“局域”指电子对出现的相对概率和成键形态。

所以二者不该互相替代。差分电荷分布问的是“形成当前构型后,电子密度相对参考态在哪里增加、在哪里减少”;ELF 问的是“当前密度和动能密度对应怎样的 局域电子结构”。同一计算结果可以同时画出 Δρ 和 ELF,但 颜色含义、零点和参考对象 都不同。

ELF适合回答哪些成键问题?

读 ELF 前要先确认 结构对象。GeO2 熔体案例里,讨论对象不是抽象的“电子更红”,而是四配位 [GeØ4]n 网络和三配位 [GeOØ2]n 链之间的局部配位差异。有了 Ge-O 距离、配位数和链/网络 motif,ELF 才能对应成键环境。

图2. GeO2 熔化前后的结构模型,以及三配位 [GeOØ2]n 链结构。DOI:10.1038/s41467-023-42890-3

在 GeO2 的 DOS 与价电子 ELF 图中,链状 [GeOØ2]n 和网络状 [GeØ4]n 的局域电子分布不同。ELF 能把 Ge-O 键区、桥氧附近电子对和非键 p 态 的空间特征联系起来,适合判断某个区域像共价键、孤对电子、离域金属电子还是弱相互作用空隙。

图3. [GeOØ2]n 链和 [GeØ4]n 网络的 DOS 及价电子 ELF 图。DOI:10.1038/s41467-023-42890-3

ELF 的优势在于成键图像直观,但不能把颜色深浅直接翻译成键能。键强 需要结合键长、振动频率、COHP/ICOHP、反键态占据和能量变化。含 d 电子体系还受 U 值、自旋态和赝势 影响,ELF 截面对比应保持相同计算设置和切片平面。

Cu 离子与磷脂模型的对比很适合说明这一点。Cu+ 形成 3d10-3d10 相互作用 时,ELF 曲线和二维截面在金属-金属方向给出局域电子信号;Zn2+ 模型则缺少类似的 金属-金属键特征。这里 ELF 支持的是“价层电子是否在两个金属之间形成局域成键区域”。

图4. Cu 离子和 Zn 离子结合磷脂模型后的电子密度、ELF 和自然键轨道对比。DOI:10.1038/s41598-018-32322-4

差分电荷分布适合回答哪些转移问题?

差分电荷分布的价值在于比较。吸附物靠近表面、掺杂原子进入晶格、两个二维材料接触之后,电子密度会重新排布。Δρ 的正负等值面 把这些变化变成空间图像:哪里积累、哪里耗尽、界面偶极大致指向哪里;参考态写得清楚,差分图才能限定电子重排位置。

NiFe 双原子催化位点的例子里,DOS 和差分电荷密度共同指向 Fe、Ni 与配位 N 原子 周围的电子重排。橙色和水青色区域描述积累与耗尽,DOS 则说明 Fe 3d、Ni 3d 和 N pz 态 在费米能级附近怎样变化。只保留差分图,很难判断空间重排来自哪类轨道。

图5. Ni-SAC、Fe-SAC 和 NiFe-DASC 的 DOS 与差分电荷密度分布。DOI:10.1038/s41467-021-24052-5

差分图也不该直接写成“电子转移越多,反应越好”。在 CO2RR 或 OER 模型中,电荷从金属位点、配位 N 或吸附 CO 之间重新分配,只能说明 吸附构型的电子云响应。真正影响反应路径的是 *COOH 形成、*CO 脱附和自由能差;Δρ 可以解释某一步的电子重排,但不能替代吸附能和自由能路径。

图6. Fe-SAC 与 NiFe-DASC 上 CO 吸附的轨道相互作用、差分电荷密度和自由能路径。DOI:10.1038/s41467-021-24052-5

两类图谱结论不一致时怎样限定判断?

当吸附后某处 Δρ 显示电子积累,但 ELF 没有形成明显键区高值时,这可能对应静电极化、弱配位或金属态离域;反过来,ELF 在某个键区较高,也可能只是局域电子对或孤对电子明显,也不能说明该区域相对参考态获得了净电子。一个看相对变化,一个看当前局域结构,冲突往往来自问题设定不同。

GeO2 熔体案例还提醒我们,ELF 读出的成键图像需要和 结构、振动谱或能量结果 互相约束。三配位链模型能解释 Raman 肩峰,ELF 和 DOS 再补充局部成键信息。若只拿 ELF 截面判断“某种键存在”,却不看 Ge-O 距离、振动模式和相对能量,结论会超出图谱本身的支持范围。

图7. GeO2 熔体的计算 Raman 光谱与 X1、X2 振动模式。DOI:10.1038/s41467-023-42890-3

二者的分工可以写得很具体:讨论吸附或界面电荷转移时,用 Δρ 说明积累、耗尽和偶极方向,再用 Bader/Hirshfeld 给出分区积分;讨论成键类型时,用 ELF 指认键区局域电子、孤对电子或金属离域特征,再用 PDOS、COHP/ICOHP 或键长变化检验轨道来源和能量贡献。

电子局域函数和差分电荷分布都来自电子结构计算,但二者尺度不同。ELF 的关键词是局域电子结构、成键区域和电子对特征;差分电荷分布的关键词是参考态、电子积累、电子耗尽和界面重排。把这两个问题区分开,读图时就不容易把“局域”误写成“转移”。