深入剖析EXAFS拟合：配位壳层与散射路径的界定、拆分依据及独立点数约束

说明：本文主要介绍同步辐射 EXAFS中配位壳层与散射路径的区别、同一壳层分出多条路径的结构来源，以及路径数量受距离分辨率和独立点数约束的方式。

同步辐射单色光跨过某种元素的吸收边时，吸收原子放出的光电子向外传播，再被周围原子回散射。出射波与回波在吸收中心叠加，形成随光电子波数 k 起伏的 χ(k)。傅里叶变换把不同振荡频率投到 R 空间，靠近某段径向距离的贡献聚成峰包络；配位壳层围绕吸收原子作径向分组，它描述一批距离相近的近邻。

散射路径是理论曲线求和时的最小单元。单重散射路径可写成“吸收原子→散射原子→吸收原子”，多重散射路径还经过第二个或更多原子。FEFF 为每条路径计算有效路径长度、简并度、背散射振幅和相位函数；元素种类、几何路线或距离分布发生改变，路径函数也随之改变。壳层负责按半径归组，路径负责写入原子身份与光电子行程。

多元素面心立方合金把两者差别放进具体结构。中心 Cr 原子的第一近邻位于相近半径，近邻却含 Fe、Co、Ni、Cr 与 Al。“第一壳层”在几何上仍是一层，EXAFS 模型里却出现 Cr–Fe、Cr–Co、Cr–Ni、Cr–Cr 和 Cr–Al 多条路径。Fe、Co、Ni 的背散射函数彼此接近，Al 的振幅与相位差异较大，这些分量在同一 R 区间发生相长或相消。

总 EXAFS 信号由全部路径分量相加，可写成 χ(k)＝Σχ_γ(k)，下标 γ 对应一条路径。每个分量的振幅由配位数、S₀²、背散射振幅、平均自由程和 σ²控制，相位由路径距离、能量零点与相移函数控制。一个峰包络可以容纳多条不同相位的曲线，峰顶数量与路径数量没有一一对应关系。

一层近邻若拥有相同元素、相同键长分布和相同散射拓扑，等价原子会合并成一条路径，原子个数通过简并度或配位数计数。只要其中一项分化，单路径平均便会丢失振幅或相位信息。常见来源包括散射原子种类不同、同种原子的键长分裂、单重与多重散射并存。

Ni 单原子位点的第一近邻可同时出现 N、B 与 C。三类原子即使位于相近距离，原子序数与电子结构仍会改变背散射振幅 f(k) 和相移 δ(k)；模型因而写成 Ni–N、Ni–B、Ni–C 三条路径。路径名称限定散射原子，每条路径再拥有自己的配位数、距离与无序度。

元素差异能否在数据中分开，还取决于散射对比。N、C、B 的 f(k) 较弱且接近，Ni–N 与 Ni–C 会形成强相关；晶体结构、化学计量和多吸收边信息能限制配位数组合。相邻过渡金属的散射函数也会高度相似，多元素合金常把若干近邻元素合成一类有效散射体，同时会丢失该壳层内的元素占位比例。

局域畸变会把化学身份相同的键分成不同距离。Jahn–Teller 畸变的 NiO₆八面体中，四条较短 Ni–O 与两条较长 Ni–O 同属第一配位壳，路径模型却要写成两组 Ni–O。短键与长键携带不同振荡相位，温度升高后两组距离靠拢，双路径模型的参数优势随局域畸变减弱而下降。

距离分辨率常用 ΔR≈π/(2Δk) 估算，其中 Δk取自可用 EXAFS 波数范围。若两组键长差远小于 ΔR，两个路径会以高度相关的 ΔR 和 σ²交换振幅，此时一条平均路径配合较大的 σ²更贴近数据容量。距离差达到分辨尺度且高 k区保留足够信噪比时，双路径会留下可重复的拍频与相位差。

单重散射沿 A→B→A 往返，多重散射沿 A→B→C→A 或更长路线闭合。两条路线的 R_eff接近时，傅里叶变换会把它们压进同一峰区。近线性键角还能产生聚焦效应，使多重散射振幅升高。路径按完整行程分类，路径所在的壳层仍按有效半路径长度归组。

卟啉环把多个碳原子放在相近径向范围，Fe–C_α、Fe–C_m、Fe–C_β单重散射及沿环骨架传播的多重散射会叠在高壳层区域。按距离和对称关系分组后，等价原子共享路径参数，非等价碳位点保留独立的相位函数。简并度合并对称等价路线，它不负责合并化学环境不同的路线。

路径拆分处理非等价贡献。一条路径的简并度可以大于 1，例如规则八面体中六个等价配体对应一条六重简并路径；畸变后四条短键与两条长键成为两条路径，简并度变成 4 和 2。“简并度”统计等价贡献，两者对应不同的结构含义。

路径数量增加时，每条路径无须都增加一整套自由参数。同一吸收边常共享 ΔE₀，S₀²可由参比样确定；晶体占位、化学计量和多吸收边配位守恒可约束配位数。距离改变量与 σ²是否共享，要服从具体键型和温度响应。共享参数必须对应同一物理来源，相邻距离不构成共享的充分条件。

多元素合金把参数规模问题放大到极端。若 n 种吸收原子都允许 n 种近邻，配位数、距离、σ²与 ΔE₀会随元素组合迅速增长。多边联合拟合、变温序列和占位对称关系把独立变量压缩到数据能够承载的范围，约束后的参数仍保留吸收原子与散射原子的对应关系。

可用信息量常按 N_ind≈2ΔkΔR/π＋2 估算。Δk由有效 EXAFS 波数范围给出，ΔR来自拟合窗口宽度；浮动变量数应低于独立点数。两个路径的振幅和相位函数很接近时，参数相关矩阵还会进一步压缩真实可分辨的信息。

Cr/Si-MFI 的第一峰在各反应气氛下保持为一个主包络，拟合模型仍含两条 Cr–O 路径，随后还有 Cr–Si 路径。气氛改变两类 Cr–O 的配位数和键长，主峰外形随路径权重连续变化。壳层是实验谱里的聚合坐标，多路径则把这段坐标还原成可计算的局域结构成分。

一组反应序列共享路径类型和全局能量偏移，配位数、键长与 σ²随气氛独立响应。若每个状态更换一套路径清单，序列中的结构变化会与模型变化混合；序列内采用相同的 Cr–O/Cr–Si 路径类型，参数随 CO₂/C₂H₆比例的趋势便对应局域配位变化。

单路径模型若在某段 k 区持续留下拍频残差，实部零点与实验谱错开，并且结构模型给出两组非等价近邻，多路径具有明确来源。双路径参数还应在 k¹、k²、k³加权及相邻窗口内保持相近距离；距离、配位数和 σ²应处于对应键型的物理区间。

若新增路径只让残差数值下降，却让两个配位数互相补偿、误差超过拟合值，或换一个 k 窗口后距离发生大幅漂移，这条路径缺少数据约束。信息不足时，平均路径与较大的 σ²保留的是径向分布宽度；把它改成两条距离相近的路径，会把噪声写成局域结构。

多吸收边拟合让同一原子对从两个方向接受约束，例如 Ru–Rh 在 Ru K 边和 Rh K 边拥有互易配位关系。每条边仍由多条近邻路径相加，联合约束减少了任意分配配位数的空间。R 空间限定径向窗口，回变换后的 q 空间限定相位与振幅，两种坐标对应同一组路径参数。

同一个壳层可由一条路径、两条分裂路径或若干异原子与多重散射路径组成。选择取决于局域结构中有多少非等价近邻，也取决于 Δk、噪声和联合数据能分开多少分量。最小且有结构来源的路径集合会同时复现实部、幅值和多种 k 加权下的振荡，配位数与键长则停留在对应原子对的数值范围内。