Joint ICTP-IAEA School on Novel Experimental Methodologies for Synchrotron Radiation Applications in Nano-science and Environmental Monitoring
Evaluation of XRF Spectra
from basics to advanced systems
Piet Van Espen [email protected]
20 Nov 2014
1 Content
1. Introduction: some concepts 2. Simple peak integration 3. Method of least squares 4. Fitting of x-ray spectra 5. Improvements to the model 6. Final remarks
2 Some basic concepts
Nature and science: a personal view
Nature = Signal + Noise
Science = Model + Statistics
3 Why spectrum evaluation?
element concentrations ⇔ net intensity of fluorescence lines
interference free continuum corrected
But: frequent peak overlap presence of a continuum
Especially in energy-dispersive spectra
4 4 Information content of a spectrum
Spectrum contains • Information: energy and intensity of x-rays the signal
• Amplitude noise: due to Poisson statistics ► fluctuations in the spectrum the noise • Energy noise: finite resolution of the detector ► nearly Gaussian peaks with a width of ~160 eV
5 Amplitude noise Counting events involves Poisson statistics
Poisson probability density function: The probability to observe N counts if the true number is µ
Property:
Poisson distribution Normal distribution ≅ µ µ and σ2 = µ
approximation is very good for µ (or N) ≥ 9
Poisson : P(N | µ=3) Normal : P(x | µ=3 σ2 = 3)
6 2 Spectrum analysis Energy Noise Resolution of ED-XRF spectrometers
2.35⇥ F E (9) Full Width at Half Maximum (FWHM) of a peak
2 2 2 FWHMPeak = FWHMElec + FWHMDet (10) 2 Spectrum analysis Electronic noise Intrinsic contribution ~100 eV 2.35⇥ F E (9) FWHM2 = FWHM2 +Det2 (10) Peakϵ energy to createElec e-h pairPealk 3.85 eV F Fano factor ~0.114 E x-ray energy in eV
Mn Kα @ 5.895 keV
FWHMDet = 120 eV
FWHMElec = 100 eV
=> FWHMPeak = 156 eV
7
2
2 Energy Noise Resolution of ED-XRF spectrometers
Cr – Mn – Fe overlap at ~20 eV Cr – Mn – Fe overlap at ~160 eV
8 Without amplitude noise (counting statistics) there would be NO PROBLEM
But it is part of the nature We can only measure longer or with a more efficient system
Without energy noise there would be LITTLE PROBLEM
The natural line width of x-rays is only a few eV!!! The observed peak width is the result of the detection process with a fundamental limitation imposed by the Fano factor
9 Information content of a spectrum
If no energy noise or no amplitude noise ► could determine the “information” unambiguously
Need methods to extract information in a optimum way These methods rely on “addition” information (knowledge) to extract the useful information
Not the method itself is important (if implemented correctly) but the correctness of the additional information.
10 10 2. Simple peak integration
11 Simple peak integration
Estimate
Uncertainty
As good as is can get We have to make assumptions if the assumptions (model) are correct! integration limits linear background no interference
12 3. Method of Least Squares
13 Method of least squares
Need to “estimate” the net peak area with highest possible • correctness (no systematic error) • precision (smallest random error)
Least-squares estimation (fitting): • unbiased • minimum variance
Limiting factors: • counting statistical fluctuations (precision) • accuracy of the fitting model
14 Method of least squares, straight line
2 SpectrumData: analysis {xi,yi}, i=1, 2, …, N 2 SpectrumModel: y(i analysis) = b + b x noise 2.35p0 1F i E (9) ⇥ ⇥ 2.35p F E (9) Fitting the model: ⇥ ⇥ estimating2 b0 and2 b1 2 FWHMPeak = FWHMElec + FWHMDet (10) 2 2 2 FWHMPeak = FWHMElec + FWHMDet (10) Criterion: SS = [y y(i)]2 = [y b b x ]2 = min (11) i i 0 1 i i SS = [y y(ii)]2 = [y b b x ]2 = min (11) X i X i 0 1 i Xi Xi ⇤SS Set of 2 equations in 2 unknowns b and b =0⇤SS yi = b0n + b1 xi (12) 0 1 ⇤b0 ! =0i yi = b0in + b1 xi Normal equations(12) ⇤b X ! X 0 i i ⇤SS X X 2 =0 xiyi = b0 xi + b1 x (13) ⇤b ⇤SS! i Direct2 analytical solution 1 =0i xiyi = b0 xii + b1 xi (13) ⇤b1 X ! X X Xi X Xi
15
2 2 4. Fitting X-ray Spectra
16 Least-squares estimate of x-ray spectrum parameters position area
1 (x x )2 Peak described by a Gaussian y(i)=b + A exp i p p2⇡ 2 2 linear parameter continuum width non-linear
n2 2 1 1 2 Criterion, agreement between model and data = [y(i) yi] (1) ⌫ wi iX=n1
y(i)=yCont(i)+ Aj RijP (i, Eij) (2) " # Minimum:Xj Xk
No direct analytical solution 2 Gain 2 1 E(i) Ejk P = G(i, Ejk)=Search χ exp for minimum (3) p2⇡ 2 jk jk " ✓ ◆ #
E(i)=Zero + Gain i (4) ⇥
17 1/2 Noise 2 = + ✏ FanoE (5) jk p jk "✓2 2 ln 2◆ #
Gain E(i) Ejk E(i) Ejk 1 T (i, Ejk)= exp erfc + 2 exp 1 p2 p2 2 2 h i
Gain E(i) Ejk S(i, Ejk)= erfc 2Ejk p2
P (i, Ejk)=G(i, Ejk)+fSS(i, Ejk)+fT T (i, Ejk) (6)
fS(Ejk)=µDet(Ejk)(a0 + a1Ejk)
fTK↵(Ejk)=b0 + b1µDet(Ejk)
fTK (Ejk)=c0 + c1µDet(Ejk)