EXAFS拟合关键解读：σ²为何现负值？从物理本质到优化策略全解析

说明：本文华算科技主要介绍EXAFS拟合中 σ² 的物理含义、负值对高 k 振幅的作用、造成参数负值的模型与数据因素，以及绝对 σ²、相对 Δσ² 和约束拟合的区分方法。

单条散射路径对应持续涨落的局域几何。热振动、静态畸变、不同局域构型会让有效路径长度形成分布，EXAFS 中的σ² 是该路径长度涨落的均方量，单位为 Ų。对单散射路径，它常被称为均方相对位移；“相对”二字强调吸收原子与散射原子的相关运动。其数值也受到两端原子运动相关性的影响。

在常用的高斯近似里，振幅因子含有 exp(−2k²σ²)。σ² 越大，高 k 区衰减越快；径向分布越宽，长程振荡越容易被平均掉。这个指数项同时解释了为何升温、缺陷增多或局域结构分散时，χ(k) 的高 k 振幅常会降低。

从数学定义看，方差由平方项求平均，绝对 σ² 的物理下限是 0。零代表理想化的单一距离且没有热涨落；真实材料在有限温度下通常保留动态贡献，结构不均一还会叠加静态贡献。因而，拟合输出的负号不能被理解成“原子振动小于零”，它指向参数化、参照基线或数据—模型匹配中的问题。

图2. Fe 箔与球磨 FeMo 在 77–300 K 的第一壳层 MSRD 和原子位移量。DOI：10.1107/S160057671700022X

这里的量还需与衍射中的原子位移参数分开。EXAFS 追踪的是一条路径两端的相对位移及其相关运动，多重散射路径则对应整条有效路径长度的涨落。把晶体学位移参数原样代入某条 EXAFS 路径，数值与物理对象都未必一致。

若把 σ² 拟合成负数，阻尼项会变成 exp(+2k²|σ²|)。振荡由衰减转为随 k² 增长的非物理高 k 增益。优化器借此补偿模型振幅偏低，这与负的热振动或负的结构无序无关。这种数学补偿会随拟合上限升高而迅速扩大。

增益对 k 很敏感。以 σ²=−0.001 Ų 为例，额外因子在 k=4 Å⁻¹ 时约为 1.03，到 k=12 Å⁻¹ 时已约为 1.33；若负值扩大到 −0.003 Ų，后一个因子约为 2.37。低 k 看似平静而高 k 被显著抬高，正是检查负值来源的一条直接线索。核验时把原始 χ(k) 与路径分量并列比较。

负号的严重程度必须连同不确定度阅读。−0.0002±0.0008 Ų与零相容，常指向信息量不足或参数相关性偏强；−0.0040±0.0003 Ų则很难归入随机波动。前者宜检查置信区间与相关矩阵，后者需立即审查路径、振幅标尺和数据处理。两类数值对应不同处置强度。

优化器只负责降低目标函数，不会自动维护每个参数的物理含义。凡是让计算振幅低于数据、让某条路径分担了不属于它的信号，或让高 k 噪声获得过大权重的环节，都可能把 σ² 推向负区间。负值因此更像一枚模型失配报警器。报警应触发参数和残差的双重审查。

配位数 N、振幅缩减因子 S₀²、路径简并度和数据归一化控制整体振幅，σ² 则施加带有 k² 斜率的阻尼。若 N 或 S₀² 设得偏低，程序可以用负 σ² 把高 k 补回来。N—S₀²—σ² 的相关性必须借助标准样、共享参数或独立化学计量信息解耦。协方差矩阵能量化这种耦合程度。

归一化误差、荧光自吸收校正过度或不同谱之间的振幅标定不一致，也会制造相似补偿。未经校正的自吸收通常压低实验振幅；只有校正过头、使实验振幅偏高时，才会朝负 σ² 方向施压。核对振幅尺度后再调无序度，比直接把 σ² 锁在零以上更能暴露源头。两种改法在残差中的指纹不同。

遗漏相邻壳层、多重散射或分裂键长后，一条高斯路径会被迫吸收多组相位与振幅。错误散射原子、错误配位结构或不恰当的 ΔE₀、ΔR 关系，也能把优化带入伪极小值。此时较低 R 因子不足以保证参数具有物理意义，残差在 k 空间和 R 空间的位置比单个总分更有诊断价值。残差的周期性还能帮助定位缺失路径。

当真实距离分布偏斜、双峰或具有长尾时，只拟合二阶累积量会把复杂线形压进一个 σ²。可在数据质量允许时测试分裂路径或高阶累积量，并比较多组 k 权重下的残差。若负值仅在某一种权重或某一小段区间出现，高 k 噪声、背景扣除和窗口边缘属于排查对象。权重变化造成的参数跳变尤其值得记录。

窄 Δk、窄 ΔR 与过多变量会让协方差矩阵接近奇异。独立信息量可用N_ind≈2ΔkΔR/π+2估算；变量数逼近这一容量时，σ² 的符号容易受起始值、参数上下限和局部极小值支配。改变初值、扩大可用 k 区间、减少自由参数并重复拟合，可以检验负值是否稳定。

软件输出栏里的 “sigma2” 未必总是绝对方差。有些温变、原位或差分分析采用 σ²(样品)=σ²(参照)+Δσ²，拟合变量其实是差值。若样品的距离分布比参照更窄，Δσ² 可以小于 0，只要相加后的总 σ² 仍不小于 0。总方差才与真实路径宽度对应。

检查从参数表达式、共享变量和基准温度开始：它究竟是 sigma2，还是 sigma2_ref + delta_sigma2。随后统一 Ų 与 10⁻³ Å² 的记法，防止把缩放后的表格数值直接代回模型。绝对量为负与差值为负是两件事，参数名省略 Δ 会让二者混在一列。核对表达式可直接消除这种歧义。

空间扫描或时间序列还应采用一致的 k、R 区间、路径集合和 S₀² 标定。只有拟合定义一致，Δσ² 的正负才对应样品相对参照的展宽变化；若每个点都自由漂移 S₀²、配位数和能量偏移，参数地图可能主要记录协方差。共享基准有助于分离真实空间差异。

第一步核对单位、绝对量或差值定义以及路径所绑定的表达式；第二步固定从标准样获得的 S₀²，检查配位数、简并度与归一化；第三步在 k 空间和 R 空间定位残差，再替换散射路径或距离分布模型；第四步改变 k 权重、拟合窗口和起始值。每次只改一类条件，参数漂移才有可比性。若多项一起变化，参数响应就难以归因。

多路径模型宜用化学计量和结构对称性建立约束，并让参数数目低于独立信息量。若某条弱路径的 σ² 反复穿过零，而删去该路径后主要壳层参数稳定、残差没有成体系地增大，这条路径在当前数据分辨率下可能不可辨识；保留一个带负方差的路径只会转移误差。

给绝对 σ² 设置 min=0 能阻止非物理解，但非负约束是一项诊断操作。若参数停在零下限、误差无法可靠估计，或配位数与 S₀² 随之大幅跳动，此前的负值仍未得到处理。应比较约束前后的残差位置、相关系数和其他路径参数；“没有负号”也不足以保证拟合可靠。

经过多初值、多窗口和多 k 权重测试后，若绝对 σ² 仍稳定处于零以下，应更换路径集合、振幅标定或距离分布模型；若仅 Δσ² 为负且总方差为正，则可按参照态收窄来理解。符号、误差、相关矩阵与残差形状四项齐全，负值的来源才能被准确归类。这也是决定保留、重建或删除该路径的依据。