3. 分峰拟合的数学基础与峰形函数

XPS分峰拟合的本质,是用若干解析函数去逼近实测谱线。选错峰形函数,即使拟合残差看起来很小,得到的化学态比例也可能严重偏离真实值。因此,理解高斯、洛伦兹、Voigt以及Doniach-Sunjic函数的物理来源与数学特性,是分峰工作的第一道门槛。

一、三种基础峰形的物理来源

理想谱线在现实中会被两类机制展宽。第一类是仪器展宽,来源于X射线源的自然线宽、分析器的能量分辨率以及电子透镜的像差,这类展宽在统计上近似服从高斯分布。第二类是寿命展宽,来源于光电子发射后内层空穴的有限寿命,由海森堡不确定性原理决定,其线型为洛伦兹分布。实测峰是这两者的卷积结果,即Voigt函数。

三种函数的数学形式决定了它们的尾部行为差异极大。高斯函数为 exp(-x²/2σ²),在偏离峰位约3σ处已衰减到峰高的1%以下,尾部收敛极快。洛伦兹函数为 1/(1+x²/γ²),衰减仅按 x⁻² 进行,在远离峰位处仍有可观强度。Voigt函数是两者的卷积,没有初等解析式,通常用Faddeeva函数或伪Voigt近似计算。

这一差异直接决定了拟合的实践后果:若用纯高斯拟合一个实际含有洛伦兹成分的峰,为了补偿长尾,拟合算法会额外引入虚假的宽峰或抬高背景,导致定量结果失真。反之,若用纯洛伦兹拟合以仪器展宽为主的谱,峰宽会被系统性高估。

高斯 / 洛伦兹 / Voigt 峰形对比(同一FWHM、同一峰高) 相对能量偏移 (eV) 强度 -6 -4 -2 0 2 4 6 高斯(仪器展宽主导) 洛伦兹(寿命展宽主导) Voigt(两者卷积) 尾部差异区 洛伦兹尾部强度远高于高斯 → 影响背景与相邻峰拟合

二、混合高斯-洛伦兹函数(GL)的参数选择

由于Voigt函数计算成本较高,日常分峰中更常用伪Voigt或混合高斯-洛伦兹函数(GL)来近似。GL函数的定义为:

GL(x) = η · L(x) + (1 - η) · G(x)

其中:
  G(x) = exp(-4ln2 · (x-x₀)² / FWHM²)     高斯分量
  L(x) = 1 / (1 + 4(x-x₀)² / FWHM²)       洛伦兹分量
  η    = 洛伦兹分数,取值 0~1
  x₀   = 峰位
  FWHM = 半高全宽

η=0时退化为纯高斯,η=1时退化为纯洛伦兹。实际拟合中,η是一个可自由优化的参数,但不应完全放任。对于金属单质峰(如Au 4f、Cu 2p₃/₂),仪器展宽占主导,η通常取0.1~0.3;对于氧化物或半导体中的离子峰,寿命展宽贡献增大,η可取0.4~0.7;对于含有未成对电子的过渡金属化合物,η可能更高。

一个常见的错误是让所有分峰共享同一个η值。实际上,同一元素的不同化学态,其空穴寿命可能不同,η应允许独立变化。但若谱图信噪比差、峰重叠严重,独立η会导致过拟合,此时可约束同一元素各分峰的η在合理范围内浮动,例如限制在0.2~0.6之间。

三、非对称峰形与Doniach-Sunjic函数

金属态谱峰常表现出明显的不对称性:峰位右侧(高结合能侧)拖尾比左侧更长。这一现象源于传导电子对光致空穴的屏蔽响应,即“费米边奇异效应”。此时高斯、洛伦兹或对称GL函数都无法正确描述,必须使用Doniach-Sunjic(DS)函数:

DS(x) = Γ(1-α) · cos[πα/2 + (1-α)·arctan((x-x₀)/γ)] / [ (x-x₀)² + γ² ]^((1-α)/2)

其中:
  α  = 不对称指数,α=0时退化为洛伦兹
  γ  = 半高半宽(HWHM)
  Γ  = 伽马函数

α的典型取值:清洁金属表面如Pt 4f、Pd 3d,α约为0.1~0.2;过渡金属如Ni、Fe的金属态,α可达0.2~0.35;若表面存在氧化或吸附,金属态比例下降,α应相应减小。α取值过大(如超过0.4)往往意味着拟合在“吃掉”背景或相邻峰,需要警惕。

使用DS函数的三个前提:第一,该峰确实对应金属态,而非氧化物或氮化物;第二,背景扣除方式与DS函数兼容,Shirley背景在强不对称峰下可能失效,建议改用Tougaard背景;第三,α必须与γ、峰面积联合约束,避免α与背景参数相互补偿导致虚假拟合。

四、函数选择的判断路径

面对一张实际谱图,可按以下顺序判断:先看峰形是否对称。若明显不对称且峰位处于金属态结合能范围,优先尝试DS函数。若基本对称,再看峰宽与仪器分辨率的关系:若峰宽接近仪器分辨率(通常0.3~0.8 eV),η取低值;若峰宽显著大于仪器分辨率,说明寿命展宽或化学环境分布贡献大,η取高值。

对于有机聚合物、绝缘体中的C 1s、O 1s峰,峰形通常接近高斯,η取0.1~0.2即可。对于过渡金属氧化物如Co 2p、Mn 2p,由于多重分裂和 shake-up 卫星峰的存在,主峰本身可能仍近似对称,但需要额外添加卫星峰分量,此时主峰用GL函数、卫星峰用宽高斯是常见策略。

样品类型 典型峰 推荐函数 η或α参考值 背景建议
清洁金属 Pt 4f、Au 4f DS α=0.1~0.2 Tougaard
金属氧化物 Cu 2p₃/₂、Zn 2p GL η=0.4~0.6 Shirley
有机聚合物 C 1s、O 1s GL η=0.1~0.2 Shirley
过渡金属化合物 Co 2p、Mn 2p GL + 卫星峰 η=0.3~0.5 Shirley
半导体/绝缘体 Si 2p、N 1s GL η=0.2~0.4 Shirley

五、实操中的参数约束与验证

以一份假设的Pt 4f谱图为例,仪器分辨率约0.45 eV,Pt 4f₇/₂峰位约71.1 eV,4f₅/₂与4f₇/₂间距约3.33 eV,面积比理论值为4:3。拟合时先固定双峰间距和面积比,仅放开峰位微调、FWHM和α。若使用DS函数,α初值设0.15,FWHM初值设0.9 eV,背景用Tougaard,参数B设为约1.5 eV。拟合后检查残差:若残差在峰位两侧呈系统性正负交替,说明α或背景参数不合理;若残差随机分布且幅值小于峰高的2%,可认为拟合可接受。

另一个关键验证是物理合理性检查:拟合得到的各化学态结合能是否落在文献合理区间?自旋轨道分裂值是否与理论一致?面积比是否接近理论值?若这些条件不满足,即使残差很小,也应重新审视峰形函数与约束条件。分峰拟合的目标不是“拟合得好看”,而是“拟合得正确”。

快速记忆:高斯管仪器,洛伦兹管寿命,Voigt是两者卷积,GL是工程近似,DS管金属不对称。选函数先问“这个峰是什么化学态”,再问“峰形对称吗”,最后用物理约束锁住参数。