XAS是一种常见的用于表征物质电子结构的方法,其常见的一种测量方法是,利用X射线照射样品,产生光电子,样品荷电,产生电流转导至地。产生光电子越多,样品荷电越多,产生电流越大。而在实际测量的XAS中,常见低能端峰位更高的现象,这便是本底信号所致。本底产生于深度激发的电子在样品中的能量损耗。为了扣除本底,出现了种种数学模型,如此处的Shirley。
Shirley本底扣除方法中,其基本假设是:高能端的电子将部分地被非弹散射,形成低能端电子的本底噪声。定量地说便是,对于初态能量为E的光电子,在到达探测器之前:
受到非弹散射而跃迁至任意能量为ϵ < E \epsilon<E ϵ < E 的末态的概率(密度)是不依赖于始末态的常数: ∀ ϵ < E , p ( E → ϵ ) = p 0 \forall\epsilon<E,p(E\rightarrow\epsilon)=p_0 ∀ ϵ < E , p ( E → ϵ ) = p 0 损失能量的电子数d n ′ \text{d}{n'} d n ′ 正比于总的电子数(d n + d n ′ \text{d}{n}+\text{d}{n'} d n + d n ′ ),其比例系数是一个不依赖于始末态的常数: d n ′ ( E ) = α ( d n ( E ) + d n ′ ( E ) ) \text{d}{n'}(E)=\alpha(\text{d}{n}(E)+\text{d}{n'}(E)) d n ′ ( E ) = α ( d n ( E ) + d n ′ ( E )) 记测量的光谱曲线为J ( E ) J(E) J ( E ) ,本底为S ( E ) S(E) S ( E ) ,原谱减去本底得到的校正谱为J ( E ) − S ( E ) J(E)-S(E) J ( E ) − S ( E ) 。校正谱为弹性散射的电子信号,处于[ E , E + d E ] [E,E+\text{d}{E}] [ E , E + d E ] 区间内的电子数d n \text{d}{n} d n 正比于该区间的信号强度(假设比例系数在很窄的能量区间内不依赖于能量,实际上,光电子强度与入射光强度、跃迁概率等均有关):
d n ( E ) = β [ J ( E ) − S ( E ) ] d E \text{d}{n}(E)=\beta[J(E)-S(E)]\text{d}{E} d n ( E ) = β [ J ( E ) − S ( E )] d E 代入基本假设2,可以解出d n ′ \text{d}{n'} d n ′ :
d n ′ = α β 1 − α [ J ( E ) − S ( E ) ] d E \text{d}{n'}=\frac{\alpha\beta}{1-\alpha}[J(E)-S(E)]\text{d}{E} d n ′ = 1 − α α β [ J ( E ) − S ( E )] d E 这些电子将按照p 0 d ϵ p_0\text{d}{\epsilon} p 0 d ϵ 为概率,为[ ϵ , ϵ + d ϵ ] [\epsilon,\epsilon+\text{d}{\epsilon}] [ ϵ , ϵ + d ϵ ] 区间内提供本底信号强度:
S ( ϵ , E ) d ϵ = α β p 0 d ϵ 1 − α [ J ( E ) − S ( E ) ] d E S(\epsilon,E)\text{d}{\epsilon}=\frac{\alpha\beta p_0\text{d}{\epsilon}}{1-\alpha}[J(E)-S(E)]\text{d}{E} S ( ϵ , E ) d ϵ = 1 − α α β p 0 d ϵ [ J ( E ) − S ( E )] d E 其中,S ( ϵ , E ) d ϵ S(\epsilon,E)\text{d}{\epsilon} S ( ϵ , E ) d ϵ 表示初态为能量E E E 的电子在能量ϵ \epsilon ϵ 处形成的那一部分本底信号强度。因而,总的本底信号是对所有E > ϵ E>\epsilon E > ϵ 的求和:
S ( ϵ ) = ∑ E = ϵ + ∞ S ( ϵ , E ) = S ( ε ) = k ∫ ε + ∞ [ J ( E ) − S ( E ) ] d E S(\epsilon)=\sum_{E=\epsilon}^{+\infty}S(\epsilon,E) =S\left( \varepsilon \right) =k\int_{\varepsilon}^{+\infty}{\left[ J\left( E \right) -S\left( E \right) \right] \text{d}E} S ( ϵ ) = E = ϵ ∑ + ∞ S ( ϵ , E ) = S ( ε ) = k ∫ ε + ∞ [ J ( E ) − S ( E ) ] d E 其中,k = α β p 0 1 − α k=\frac{\alpha\beta p_0}{1-\alpha} k = 1 − α α β p 0 。注意到若能量为ϵ \epsilon ϵ 处并没有峰,则该处的信号均为本底,J ( E ) − S ( E ) = 0 J(E)-S(E)=0 J ( E ) − S ( E ) = 0 ,因而:
若ϵ \epsilon ϵ 处有峰。假设E b E_b E b 为峰的高能边(边界),则: S ( ε ) = k { ∫ ε E b [ J ( E ) − S ( E ) ] d E + ∑ E j > ε ∫ E j E j + Δ E j [ J ( E ) − S ( E ) ] d E } S\left( \varepsilon \right) =k\left\{ \int_{\varepsilon}^{E_b}{\left[ J\left( E \right) -S\left( E \right) \right] \text{d}E}+\sum_{E_j>\varepsilon}{\int_{E_j}^{E_j+\Delta E_j}{\left[ J\left( E \right) -S\left( E \right) \right] \text{d}E}} \right\} S ( ε ) = k ⎩ ⎨ ⎧ ∫ ε E b [ J ( E ) − S ( E ) ] d E + E j > ε ∑ ∫ E j E j + Δ E j [ J ( E ) − S ( E ) ] d E ⎭ ⎬ ⎫ 其中,E j E_j E j 为其他峰的低能边。上式不难化简为:
S ( ϵ ) = k Q ( ϵ ) + b ′ S(\epsilon) =kQ(\epsilon)+b' S ( ϵ ) = k Q ( ϵ ) + b ′ 式中,Q ( ϵ ) = ∫ ϵ E b [ J − S ] d E Q(\epsilon)=\int_\epsilon^{E_b}[J-S]\text{d}{E} Q ( ϵ ) = ∫ ϵ E b [ J − S ] d E 为校正谱处于ϵ \epsilon ϵ 与E b E_b E b 之间的面积,而b ′ b' b ′ 则取决于该峰与其他峰的相对位置。在该峰区域内,b ′ b' b ′ 是不依赖于ϵ \epsilon ϵ 的常数。
若ϵ \epsilon ϵ 处无峰。则: S ( ε ) = k ∑ E j > ε ∫ E j E j + Δ E j [ J ( E ) − S ( E ) ] d E S\left( \varepsilon \right) =k\sum_{E_j>\varepsilon}{\int_{E_j}^{E_j+\Delta E_j}{\left[ J\left( E \right) -S\left( E \right) \right] \text{d}E}} S ( ε ) = k E j > ε ∑ ∫ E j E j + Δ E j [ J ( E ) − S ( E ) ] d E 这是一个不依赖于ϵ \epsilon ϵ 的常数(只要ϵ \epsilon ϵ 在变化时不跨越任何有峰的区域)。这表明,按照Shirley的假设,两个峰位之间应该均是水平的区域,这些水平的区域会随着能量的降低,像台阶一样逐步升高(每跨越一个峰,均是上一个台阶)。
至此,我们描绘了Shirley假设下,光电子谱应该有的峰形。在使用Shirley时,我们或许需要尽量保证原谱符合其描述的情况:即峰位左右侧为水平、且高能边低于低能边。否则,贸然采用Shirley或许并不合理。
下面是求解Shirley的算法部分。首先假设某峰的高能边(E b E_b E b )和低能边(E a E_a E a )的高度J ( E a ) = S ( E a ) = a , J ( E b ) = S ( E b ) = b J(E_a)=S(E_a)=a,J(E_b)=S(E_b) = b J ( E a ) = S ( E a ) = a , J ( E b ) = S ( E b ) = b 。结合S ( ϵ ) = k Q ( ϵ ) + b ′ S(\epsilon)=kQ(\epsilon)+b' S ( ϵ ) = k Q ( ϵ ) + b ′ ,可得:
b ′ = b , k = a − b I , I = ∫ E a E b [ J ( ϵ ) − S ( ϵ ) ] d ϵ b'=b, k=\frac{a-b}{I},I=\int_{E_a}^{E_b}[J(\epsilon)-S(\epsilon)]\text{d}{\epsilon} b ′ = b , k = I a − b , I = ∫ E a E b [ J ( ϵ ) − S ( ϵ )] d ϵ 其中I I I 便是校正谱中,该峰的总面积。不妨定义泛函:
ψ [ S ( ϵ ) , E ] = ∫ E E b [ J ( ϵ ) − S ( ϵ ) ] d ϵ \psi[S(\epsilon),E]=\int_{E}^{E_b}[J(\epsilon)-S(\epsilon)]\text{d}{\epsilon} ψ [ S ( ϵ ) , E ] = ∫ E E b [ J ( ϵ ) − S ( ϵ )] d ϵ 则有:
S ( ϵ ) = b + ( a − b ) ψ [ S ( ϵ ) , E ] ψ [ S ( ϵ ) , E a ] S(\epsilon)=b+(a-b)\frac{\psi[S(\epsilon),E]}{\psi[S(\epsilon),E_a]} S ( ϵ ) = b + ( a − b ) ψ [ S ( ϵ ) , E a ] ψ [ S ( ϵ ) , E ] 上式的右侧可以视为S ( ϵ ) S(\epsilon) S ( ϵ ) 的泛函,它将S ( ϵ ) S(\epsilon) S ( ϵ ) 映射为另一函数。而这一映射的结果得到了S ( ϵ ) S(\epsilon) S ( ϵ ) 自身,这表明,上式的解便是这一泛函的不动点。在某些情况下,我们可以通过不断迭代函数来得到函数的不动点,如:
a 1 = f ( x ) , a 2 = f ( f ( x ) ) , ⋯ , f ( a ∞ ) = a ∞ a_1=f(x),a_2=f(f(x)),\cdots,f(a_{\infty})=a_\infty a 1 = f ( x ) , a 2 = f ( f ( x )) , ⋯ , f ( a ∞ ) = a ∞ 算法求解Shirley时,便是采用这一迭代法。