如何利用密度泛函理论（DFT）复现同步辐射 EXAFS 的配位环境？

拿到同步辐射扩展 X 射线吸收精细结构（EXAFS）谱后，常见做法是拟合出金属—配体键长，再与密度泛函理论（DFT）优化结果对照。如果两个距离接近，是否就找到了真实配位结构？问题在于，多种结构可能具有相近的第一壳层距离，却产生不同的高壳层和多重散射信号。DFT 与实验之间还缺一步：把候选结构计算成 EXAFS 谱，再比较它能否解释实验中实际可分辨的特征。

DFT 几何优化给出特定计算条件下的原子坐标、能量和电子结构；EXAFS 测量的是吸收原子附近的光电子散射干涉。常用路线是“DFT 结构—散射计算—理论 EXAFS—实验对照”。FEFF 等程序根据局域坐标计算散射振幅与相位。仅有态密度、Bader 电荷或一幅优化后的球棍图，都还没有完成 EXAFS 的正向计算。

对近似高斯分布的单次散射壳层，常用表达式为：

χ(k) = Σ_j [N_jS₀²F_j(k)/(kR_j²)]
× exp[−2k²σ_j² − 2R_j/λ(k)]
× sin[2kR_j + φ_j(k)]

这里 χ(k) 是无量纲振荡信号，k 是光电子波数，单位 Å⁻¹；R 是吸收原子到邻居的距离，单位 Å；N 是该壳层配位数。F 与 φ 是散射振幅和总相位，λ 是光电子平均自由程，S₀² 为振幅缩减因子，σ² 是原子间距离的均方相对涨落，单位 Ų。多重散射需要按实际路径处理，不能把每条路径都当成一个配位壳层。

这条式子说明了为何“距离差不多”不够：距离影响振荡相位，邻居种类影响 F 和 φ，配位数及无序共同影响振幅。未经相位校正的傅里叶变换峰位不是实际键长，不能拿图上的 R 峰顶直接与 DFT 距离相减。应使用散射模型得到的距离，并核对 k 空间振荡，而非只比较一个峰的位置。

以碳载体上的金属位点为例，可以比较不同氮配位形式、缺位结构、轴向吸附物，以及合理的金属邻近结构。候选范围由合成、元素组成、显微观察和测量状态约束，不必穷举所有想象中的模型；但只算一个预设结构，也无法检验其他解释。每个模型应写清它预计改变哪一组散射信号。

模型的电荷、自旋、配体和环境同样重要。干燥粉末、浸在电解液中的样品与工作电位下的催化剂，可能有不同的轴向配位。计算裸露位点去解释存在水或反应物的谱，失配未必来自 DFT 精度。过渡金属应检查合理的自旋初态；泛函、色散修正和必要的关联处理，要根据体系选择并检验敏感性，不能照搬另一篇论文的参数。

优化时检查截断能、k 点、超胞和力收敛对局域距离的影响；表面模型还需处理层数、固定层及周期镜像。DFT 优化出的整数邻居数，是该模型的几何属性；EXAFS 拟合配位数则受振幅标定与结构平均影响。二者数值不同，不应立即解释成“部分原子失去了配体”。

图1｜Co、Fe、Ni 单原子催化剂的结构分析。a—c 为实验近边谱，d—f 为实验与理论近边谱对照，g—i 为 EXAFS 拟合及金属—氧、金属—氮路径分量。

来源：10.1038/s41467-021-27143-5

图1的 d—f 是 X 射线吸收近边结构（XANES）的实验与模拟对照，g—i 才是 EXAFS 拟合及 M—N、M—O 路径分量。原研究比较了不同氮配位与轴向含氧模型；部分候选能改善近边特征，却与 EXAFS 距离不符。这个例子说明，XANES 的吻合不能替代 EXAFS 检验，同一结构要接受不同观测量的约束。

此外，原子数或配体数不同的模型不能直接按总电子能排序。涉及吸附、脱附或不同化学组成时，须定义相应反应与化学势；即使某模型能量较低，也不能排除合成过程中保留下来的亚稳结构。

导出结构时，把待测元素的一个原子设为吸收中心，展开周期邻居，再截取包含目标壳层的原子簇。晶胞边缘的邻居必须通过周期平移补齐，薄片的真空区则不应被误当成缺失配位。FEFF 输入中的吸收边、元素势与原子坐标，需要与实验及候选模型一致；DFT 波函数通常不会因为导出坐标就自动成为 FEFF 的散射势。

簇半径、自洽势范围、最大路径长度和多重散射阶数应分别测试。逐步增大它们，观察目标 k、R 区间的理论谱是否稳定。只想分析第一壳层时不一定需要大量长路径；但配体骨架或近直线原子排列的多重散射可能影响较远区域，不能仅凭“都是轻元素”就删掉。

FEFF 路径带有有效半路径长度 R_eff 和简并度。单次散射的一组等价路径，其简并度可对应该组邻居数；多重散射的简并度描述等价路径条数。若软件已包含四条等价 M—N 路径，振幅表达式再额外乘四，就会重复计数。允许配位数变化时，应明确定义相对原始简并度的比例，或采用单位简并度路径乘目标 N 的一致约定。

图2｜液氮温度下金属铜的实验 EXAFS（空心圆）与计算结果（实线）。原研究计算前将理论晶格参数放大 1%；纵向箭头表示振幅尺度。

来源：10.1107/S1600577524005484

图2用已知结构的金属铜检验计算路线。原研究在计算谱前将理论晶格参数放大了 1%，图注明确记录了这一调整。对自己的模型，类似修正也要保留原始结果、说明依据，并让不同候选遵守一致规则。先用已知结构检验计算，再比较未知位点，有助于判断谱的系统偏差是否被误归因于结构错误。

静态 DFT 优化通常对应某一局部能量极小值，而测量包含许多吸收原子及随时间变化的构型。对于相对刚性的结构，可用静态路径配合有物理依据的 σ² 描述阻尼；溶液配合物、柔软配体或明显重构的位点，则可能需要第一性原理分子动力学（AIMD）提供构型分布。

AIMD 应在与目标状态相适应的温度及环境下运行。舍弃平衡前轨迹，检查局域距离、配位和温度是否稳定，再从生产轨迹抽取构型。相邻帧往往高度相关，增加文件数量并不等于增加独立采样；应比较不同时间段及抽帧间隔下的平均谱是否收敛。若短轨迹未跨越配体交换或重构势垒，所得平均只代表已采样的结构区域。

对每个瞬时构型、每个相关吸收中心计算 χ(k)，按相应权重求平均：

χ_avg(k) = Σ_s w_sχ_s(k)，Σ_s w_s = 1

s 表示构型与吸收位点，w 为归一化权重。先平均 χ(k)，再傅里叶变换并取模。先把原子坐标平均成一幅结构，会丢失距离分布；先对每帧傅里叶变换的模求平均，则会丢失构型之间的相位抵消。

图3｜水溶液中钷配合物的光谱和计算结构。a 为近边谱分析；b、c 对照实验 EXAFS、路径拟合和 AIMD 计算谱；d 为水环境中的构型；e、f 为成键轨道分析。

来源：10.1038/s41586-024-07267-6

图3的 b、c 同时给出实验、路径拟合和 AIMD 计算结果，d 展示显式水环境中的构型。原文从平衡后的轨迹取样并用 FEFF 计算，因此曲线对应构型集合，不能用 d 面板那一个瞬时结构代表。还需注意：路径拟合依据配体与晶体结构先验，固定了第一壳层 N = 9；它检验了这一模型下的距离等参数，并非仅靠 EXAFS 振幅独立测出九配位。

图4｜AIMD 得到的钷周围氧（a）、碳（b）和氮（c）的径向分布。红线对应左轴的径向分布函数，蓝线对应右轴的积分配位数。

来源：10.1038/s41586-024-07267-6

图4的 a 将氧原子的径向分布与积分配位数放在一起。积分在合适的壳层边界内回答“平均有多少邻居”，峰的宽度则反映距离分布；边界应结合第一极小值和配体归属选择。这类分布可以检查静态模型漏掉了哪些涨落，但 RDF 与 EXAFS 并不是同一种信号，后者还包含元素相关的散射振幅和相位。

逐帧计算已经包含轨迹采样到的几何无序，不能再原样叠加同一份热无序阻尼。额外 σ² 只能用于解释确实缺失且有依据的贡献。低温下经典核运动的 AIMD 还可能遗漏零点振动，需要另行评估；以任意加热轨迹去压低振幅，不能替代物理分析。

实验端先核对能量标定、背景扣除、归一化、异常点和重复扫描。荧光自吸收、探测器非线性或束流损伤可能改变振幅；若把这些变化直接交给 N 和 σ²，结构参数就会吸收仪器误差。理论与实验应使用相同的 k 范围、权重、窗函数及窗宽，随后在同一 R 范围比较。

S₀² 宜由合适的已知配位标准约束，并检查其转移到目标体系的适用性。能量原点修正 ΔE₀ 与距离变化会相关；同一吸收边、同一数据集中的路径通常共享合理的 ΔE₀。不要为每条路径分别放开能量偏移，直到所有候选都能贴上实验。

两个近邻壳层能否分开，受有效 k 范围限制。作为教学估算，若 k 从 3 到 11 Å⁻¹，则 Δk = 8 Å⁻¹，径向分辨率约为 π/(2Δk) ≈ 0.20 Å。DFT 中相差 0.05 Å 的两组距离，未必能在实验中独立解析。这个分辨率也不同于在指定模型下拟合单一平均键长的统计精度。

拟合信息量常用 N_ind ≈ 2ΔkΔR/π + 2 估算，其中 ΔR 是拟合区间宽度，单位 Å；不同约定的常数项略有差别。沿用上述 Δk，若 ΔR = 1.5 Å，独立信息量只有约 10。增加插值点或同时采用多种 k 权重，都不会凭空产生同等数量的新信息。放开多个壳层的 N、R、σ² 和各自的 ΔE₀，很容易超过数据承载能力。

比较候选时，同时看 k 空间的振荡、复数傅里叶变换的实部和虚部、模及残差。只看模的峰形，可能遗漏相位错误。候选使用同样的数据区间与相近的自由参数约束，再检查参数相关性、改变窗口后的稳定性及噪声水平。拟合残差略小而 σ² 为负、振幅异常或距离修正过大时，应先排查模型与处理条件。

对有取向的薄膜、单晶或二维材料，谱强度还取决于 X 射线偏振与键方向。计算粉末平均谱去比较某一偏振方向的实验，可能把方向差异误判为配位数降低。

图5｜室温下 WSe₂ 的 Se K 边 EXAFS。空心圆为实验，实线为计算；上方对应偏振垂直表面，下方对应偏振平行表面，两组曲线作了纵向平移。

来源：10.1107/S1600577524005484

图5上下两组分别对应偏振垂直和平行样品表面的测量，计算也保留了方向条件。两组振荡的形状与振幅存在差异，说明复现局域配位还需复现测量几何；讨论“少了几个邻居”之前，应检查样品取向。

若样品含有多种金属位点，可在可比的吸收边归一化条件下写成 χ_exp(k) ≈ Σ_a p_aχ_a(k)。p 对应吸收原子在不同环境中的比例，不能直接使用载体或整体相的质量分数。混合谱应在线性 χ 层面组合。只有各候选谱具有可分辨差异时，拟合出的比例才有意义；M—N 与 M—O 等相近轻元素贡献，往往还需要化学组成、XANES 或其他表征约束。

实施时，先保存实验处理参数和可用范围；随后为每个候选记录组成、自旋、环境、优化条件与坐标。逐一生成散射路径，核对周期邻居、简并度和路径收敛，再决定采用静态无序模型还是构型平均。最后在固定的对照条件下计算与拟合，记录残差、参数不确定性及被排除候选的具体失配特征。

遇到失配时，可按特征查原因：振荡位置持续错开，检查距离与能量原点；高 k 阻尼不对，检查无序、采样与实验振幅失真；较远壳层缺峰，检查邻居范围、多重散射和模型是否漏掉结构组分。每次修改保留修改前的谱，避免几何、振幅和背景一起改变后无法判断改善来自哪里。

最终结论应对应比较的强度。如果一个候选在受约束的条件下复现了多个区间，并明显优于合理替代结构，可以说实验支持该配位模型。若几个候选在噪声内无法区分，就保留可兼容的结构范围。DFT 提供可检验的原子结构，EXAFS 决定这些结构在当前数据下能被支持到什么程度；有时可信结果是一组环境，而不是唯一的球棍图。