什么是电子局域函数ELF？

说明：本文华算科技主要介绍电子局域函数的概率基础、拓扑盆、结构演化信号及定量证据范围。

为什么由同自旋条件概率构造？

把两个自旋相同的电子放到很近的位置，泡利不相容原理会压低第二个电子出现的概率。电子局域函数（electron localization function，ELF）追踪的正是这段短程概率下降得有多快。某处的同自旋电子彼此回避越显著，电子在该处形成局域电子对或局域电子团的倾向越强。这里的“局域”描述同自旋条件概率的空间曲率，和某一点含有多少电子属于两个问题。

单行列式波函数中，这个曲率可写成泡利动能密度。对自旋分量 σ，常用表达为 D_σ=τ_σ−|∇ρ_σ|²/(4ρ_σ)；再除以同密度均匀电子气的参考值 D_σ⁰，得到无量纲比值 χ_σ。不同程序对动能密度 τ 的定义可能相差 1/2，公式前因子也随之改变，读取软件输出时应核对程序采用的约定，避免把符号差异当成物理差异。

ELF_σ(r) = 1 / [1 + χ_σ(r)²]，　χ_σ = D_σ / D_σ⁰

这个变换把数值压到 0—1。ELF 接近 1，表示相对于同密度均匀电子气，同自旋电子在近距离出现的概率受到更强抑制；ELF=0.5 对应均匀电子气参考。接近 0 的区域表示泡利动能相对较高，它不能单独充当“缺电子区”的标签。电子密度很低的真空区还会放大数值噪声，绘图时通常配合密度截断。

图1上排用 ELF=0.75 的等值面显示局域域，下排标出吸引子和盆电子数。C—C 单键、双键与羰基附近出现不同的局域域，氧原子旁还能看到孤对电子对应的非成键域。路易斯酸配位后，O—Al 方向增加新的局域结构，原有羰基盆的电子数也发生改变。图中的等值面负责显示形状，红色数值则来自对拓扑盆内电子密度的积分，两者承担的任务不同。

ELF吸引子和盆怎样对应电子对？

ELF 标量场中的局部极大值称为吸引子。从空间任一点沿 ∇ELF 上升，轨迹会汇入某个吸引子；具有同一终点的轨迹集合构成一个盆，相邻盆由零通量面分开。吸引子给出局域中心，盆给出这片局域结构占据的空间。固定阈值绘出的彩色块只是盆的一部分，阈值改变会让表面膨胀、收缩甚至连在一起。

围绕原子核的盆记作 C(A)，主要对应内层电子；只与一个原子相连的价层盆 V(A) 常见于孤对电子、自由基或原子附近的非键电子；位于两个原子之间的 V(A,B) 常与共享电子对相联系。多中心键可出现 V(A,B,C) 等多中心盆。盆的命名依据拓扑连接和核的位置，化学归属还要核对盆的位置、电子数以及邻近原子的轨道组成。

盆划定后，可计算 N(Ω)=∫_Ωρ(r)dr，得到盆内平均电子数，并由盆人口涨落和盆间协方差分析电子共享。盆电子数会随极化、相关效应和零通量面的位置变化而偏离整数。它记录某个实空间区域收纳了多少电子，数值本身不直接给出键级或解离能。

图2a在表面羟基与水合氢离子之间画出 ELF 截面，氢原子两侧的局域区反映 O—H 键和界面氢键环境。图2b固定了这段空间对应的原子构型，图2c给出 Bader 电荷，图2d用双电层电容最低点估计零电荷电位。ELF 在这里定位局域电子对的空间关系；电荷转移方向和界面电化学响应分别由后两类量承担，单张 ELF 色彩图不包含这些定量结果。

BiPbTi₂O₆ 给出 V(A) 型盆的晶体实例。图3a—c中，Bi 和 Pb 一侧出现偏心的 6s² 孤对电子叶瓣，极性位移由此获得明确的实空间方向；中心对称参考结构中的叶瓣恢复近似球形。下排能量随归一化位移 λ 单调降低，才把孤对电子形状与极性结构的能量稳定联系起来。ELF 负责辨认偏心孤对电子，能量曲线负责检验位移是否有利。

盆还可能位于原子核之外。图4中的 YCl、Y₂Cl₃、Sc₅Cl₈ 和 Sc₇Cl₁₀ 由金属 R₆ 八面体组成层、单链或双链，空隙几何预先限定了非核局域中心可能出现的位置。几何空腔本身没有电子信息；只有 ELF 吸引子确实落入空腔，并且后续电子结构量支持该归属时，才能讨论间隙电子。

结构变化时ELF拓扑会怎样演化？

比较一组结构时，等值面的体积变化容易受绘图阈值影响。拓扑分析关注另一组事件：吸引子的出现或消失、盆的分裂或合并，以及连接两个盆的鞍点临界值。连续改变键长、压力或反应坐标时，这些事件标记局域电子组织方式发生改变。采样间隔过大可能越过短寿命盆，网格过粗也可能制造或抹去小盆。

图5把 BP 的闪锌矿、C2/m 和 P4₂/mnm 三种相放在一起。左列电子密度截面给出 B、P 周围的空间分布，右列 ELF 等值面显示压力相变后 P 原子之间形成锯齿链及后续三维连接。三幅 ELF 使用的阈值分别为 0.8、0.65 和 0.7，黄色区域的体积不能横向当作局域电子多少的排序；可信的比较对象是连接位置、吸引子拓扑和盆积分。

化学反应路径中，成键盆往往在两个反应中心靠近后出现，随后盆人口增长；断键过程则可能出现相反序列。若两条键不同步形成，两个新盆的出现位置也会错开。拓扑事件标记电子对重新分区的几何位置，反应势垒仍来自能量沿路径的变化。

图6沿内禀反应坐标列出过渡态前后多个结构。紫色吸引子及其盆人口记录 C₁—C₆ 与 C₄—C₅ 两条新键的电子重排；两组数值没有在同一结构处同步进入成键状态，对应异步环加成。这里应同时保留结构序列、盆标号和盆电子数，单独截取某个过渡态的等值面会丢失“何时出现新盆”这一信息。

哪些结论要交给盆积分/轨道/能量？

非核 ELF 极大值常见于电子化合物，也可能来自邻近原子 d 轨道的相长杂化。两类来源在彩色截面上都能形成高值斑块。区分它们时，应检查吸引子所在空隙、盆人口、相关能带的部分电荷密度以及 PDOS；若要确认电子密度自身存在非核吸引子，还可增加 QTAIM 拓扑分析。

图7中，位点 A 位于金属 R₄ 四面体中心，位点 B 位于 R₆ 八面体中心，两处都出现非核吸引子。截面中的有限 ELF 连接还显示 A、B 相距很近。论文随后发现 A 的空间太小，其局域峰主要来自 Y/Sc–d 轨道杂化；B 与较大的八面体空隙和间隙能带对应。“远离原子核”只给出候选位置，轨道来源决定它属于金属键局域还是阴离子型间隙电子。

图8分别显示费米能级附近三条能带。带 1 和带 2 的部分电荷密度呈 Y–d_xz、Y–d_yz 轨道形成的 Y—Y π 键，对应图7的 A 位点；带 3 的电荷密度沿 R₆ 层或链的空隙铺展，对应 B 位点的间隙电子。能带色散还区分二维层状和一维链状传输方向，ELF 截面没有携带这部分动量空间信息。

成键强弱可由总能差、键解离能、COHP/COOP 或力常数度量；电荷得失可由差分电荷密度与 Bader 电荷度量；局域电子对、孤对电子方向和成键盆演化则由 ELF 拓扑描述。每个量只回答它的定义所覆盖的问题，把盆电子数写成键级、把等值面大小写成键能，都会越过证据范围。

不同结构之间作 ELF 比较，还应统一赝势、价电子配置、交换关联泛函、自旋设置、实空间网格和绘图阈值。赝势计算往往不含显式芯电子，核附近的 ELF 只代表所选价电子空间；强磁性体系宜查看自旋分辨结果；低密度区应设置电子密度掩膜。完成这些控制后，ELF 才能把同自旋概率、实空间盆和化学结构放入同一套可复核的坐标中。