Rietveld Analysis of X-ray and neutron diffraction patterns
Download
Report
Transcript Rietveld Analysis of X-ray and neutron diffraction patterns
Rietveld Analysis of X-ray and neutron
diffraction patterns
Analysis of the whole diffraction pattern
Profile fitting is included
Not only the integrated intensities
Refinement of the structure parameters from diffraction data
Quantitative phase analysis
Lattice parameters
Atomic positions and occupancies
Temperature vibrations
Grain size and micro-strain (in the recent versions)
Not intended for the structure solution
The structure model must be known before starting the Rietveld refinement
1
Non-refinable parameters in the Rietveld
method
Space group
Chemical composition
Analytical function describing the shape of the diffraction
profiles
Wavelength of the radiation (can be refined in Fullprof or in
LHRL; suitable for the synchrotron data)
Intensity ratio in Ka1, Ka2 doublet
Origin of the polynomial function describing the background
2
Rietveld analysis
History
Computer programs
H.M. Rietveld - neutron data, fixed
wavelength
D.E. Cox - X-ray data
R.B. Von Dreele - neutron data, TOF
D.B. Wiles & R.A. Young - X-ray data, 2
wavelengths, more phases
Helsinki group - spherical functions for
preferred orientation but a single
wavelength
Fullprof, LHRL - surface absorption
BGMN - automatic calculation, crystallite
size and microstrain in form of ellipsoids
P. Scardi et at - size, strain
H.M. Rietveld
DBW2.9, DBW3.2 (Wiles & Young)
University of Helsinki
Fullprof (J. Rodriguez-Carvajal)
BGMN (R. Bergmann)
LHRL (C.J. Howard & B.A. Hunter)
P. Scardi et al.
Bärlocher
GSAS
3
Integral intensity
Calculated intensity:
yic yib Gikp I k
p
k
G is the normalised profile function, I is the intensity of the k-th reflection. The
summation is performed over all phases p, and over all reflections contributing to
the respective point.
The intensity of the Bragg reflections
I k Sm k Lk Fk Pk Ak Ek
2
4
Scattering by one elementary cell
Structure factor
n
Fk N j f j exp 2ih tk r j 2 2h tk B j h k
j 1
Fk N j f j exp 2i hx j ky j z j
n
j 1
exp 11h 2 22 k 2 33 2 2 12 hk 2 13 h 2 23 k
Calculation is performed in the oblique axes (for the respective crystal
system)
5
Temperature vibrations
Atomic displacement (in Cartesian co-ordinates)
B j u j u tj
u12
u1u2
u1u3
u1u3
u 2 u3
2
u3
u1u2
u22
u 2 u3
1 a* cot * a*
1 t
B 2 F βF; F 0 1 b* sin *
2
0
0
a cos
b cosa
c
6
Crystal symmetry restrictions
Six anisotropic temperature factors per atom in a general case
(symmetrical matrix)
For an atom in a site of special symmetry the B-matrix must be invariant
to the symmetry operations (in the Cartesian axis system)
P t BP B
An example - rotation axes parallel with z
cosa
P sin a
0
sin a
cosa
0
0
0
1
7
Temperature vibrations - special cases
Isotropic atomic vibrations
sin 2
Fk N j f j exp 2i hx j ky j z j exp B j
2
j 1
n
B j 8 2 u 2
j
Overall temperature factor
2
n
2
2 sin
N j f j exp 2i hx j ky j z j
Fk exp 8 u
2
j 1
8
Scattering by one atom
Atomic scattering factor
sin 2
c Df Df
f ai exp bi
2
i 1
4
a, b, c are from the “International Tables for Crystallography”
Df’, Df” must be checked and changed for synchrotron radiation
Another possibility
Include our set of the atomic scattering factors
9
Preferred orientation of grains (texture)
Gauss-like distribution
1 G exp G sin a
1 G exp G sin a
Pk G2 1 G2 exp G1a k2
Pk G2
Pk G2
2
2
1
k
3
2
1
k
March-Dollase correction
1
Pk G12 cos2 a k sin 2 a k
G1
3
2
Spherical functions
10
Absorption correction
0.0
For flat samples - micro-absorption
and surface absorption (Hermann &
Ermrich)
1 P0 1 a 0
1
1
sin sin
Apparent decrease of the
temperature factors or even
“negative” temperature factors
ln (Intensity ratio)
Ak 1 P0 Ps ( )
-0.1
-0.2
-0.3
-0.4
# 1
# 2
-0.5
-0.6
0.0
0.1
0.2
(sin / )
0.3
2
11
Absorption correction
For thin samples (powder on glass) in
symmetrical arrangement
thick sample, high absorption
t : A 1 (2)
thin sample, low absorption
t 0: A t sin
experim ental data
absorption factor
apparent tem perature
log (Intensity ratio)
1
2t
I I0
1
exp
2
sin
0.0
-0.1
-0.2
-0.3
-0.4
0.00
0.05
0.10
0.15
0.20
0.25
(sin / )
0.30
0.35
0.40
2
12
Extinction correction
(for large crystallites)
Ek EB cos2 k EL sin 2 k
Extinction for the Bragg case ( = 90)
1
EB
1 x
Extinction for the Laue case ( = 0)
F
x D k
Ve
2
x x 2 5x3
1
for x 1
2 4 48
EL
2 1 1 3 15 for x 1
2
3
x 8 x 128x 1024x
13
Profile functions
Gauss
Lorentz (Cauchy)
C
2
exp 02 2 i 2 k ; C0 4 ln 2
k
k
C0
G
L
2 C0
1
; C0 4
k 1 C0 2 2 2
i
k
k2
C0
k
4 2 1
2
1
2
2
i
k
2
k
m
Pearson VII
PVII
Pseudo-Voigt
pV L 1 G
m
; C0
1
2
2 m 2 1
m 0.5
m
k2 U tan2 k V tan k W
14
Background
Subtraction of the background intensities
Interpolation of the background intensities
Polynomial function (six refinable parameters)
Origin of the background - improves the pivoting of the normal matrix
A special function for amorphous components
sin B2 m1Q
yib B0 B1Q B2 m
B2 m1Q
m 1
n
15
Minimisation routine
Uses the Newton-Raphson algorithm to minimise the quantity
R wi yio yic ; wi 1 yi
2
i
Normal matrix
M Dx D y
yic yic
M mn wi
i
xm xn
y
ym wi ic yio yic0
xm
i
wi yio yic
2
1
m M mm
i
N P
16
Reliability factors
The profile R-factor ………
y y
y
io
Rp
ic
i
i
The weighted Rp ………………………………………
I I
I
ko
The Bragg R-factor ………
RB
wi yio yic
i
wi yio2
i
2
io
Rwp
kc
i
ko
i
The expected Rf ………………………………………
The goodness of fit
GOF
1
2
w y
i
yic
2
io
i
N P
Rwp
R
exp
Rexp
2
N P
w y2
i io
i
1
2
17
Connecting parameters, constrains
Young - parameter coupling
Coding of variables: number of the parameter in the normal matrix +
weight for the calculated increment
Lattice parameters in the cubic system: 41.00 41.00 41.00
Fractional co-ordinates at 12k in P63/mmc, (x 2x z): 20.50 21.00 31.00
Fullprof - constrains
Inter-atomic distances may be constrained
BGMN - working with molecules
Definition of the molecule (in Cartesian co-ordinates)
Translation and rotation of the whole molecule
18
Structure of the input file
(Fullprof for anglesite)
COMM PbSO4 D1A(ILL),Rietveld Round Robin, R.J. Hill,JApC 25,589(1992)
!Job Npr Nph Nba Nex Nsc Nor Dum Iwg Ilo Ias Res Ste Nre Cry Uni Cor
1 7 1 0 2 0 0 0 0 0 0 0 0 0 0 0 0
!
!Ipr Ppl Ioc Mat Pcr Ls1 Ls2 Ls3 Syo Prf Ins Rpa Sym Hkl Fou Sho Ana
0 0 1 0 1 0 0 0 0 1 6 1 1 0 0 1 1
!
! lambda1 Lambda2 Ratio Bkpos Wdt Cthm muR AsyLim Rpolarz
1.54056 1.54430 0.5000 70.0000 6.0000 1.0000 0.0000 160.00 0.0000
!NCY Eps R_at R_an R_pr R_gl Thmin
Step
Thmax PSD Sent0
5 0.10 1.00 1.00 1.00 1.00 10.0000 0.0500 155.4500 0.000 0.000
!
! Excluded regions (LowT HighT)
0.00 10.00
154.00 180.00
!
34 !Number of refined parameters
!
! Zero Code Sycos Code Sysin Code Lambda Code MORE
-0.0805 81.00 0.0000 0.00 0.0000 0.00 0.000000 0.00 0
! Background coefficients/codes
207.37 39.798 65.624 -31.638 -90.077 47.978
21.000 31.000 41.000 51.000 61.000 71.000
! Data for PHASE number: 1 ==> Current R_Bragg: 4.16
PbSO4
!Nat Dis Mom Pr1 Pr2 Pr3 Jbt Irf Isy Str Furth ATZ Nvk Npr More
5 0 0 0.0 0.0 0.0 0 0 0 0 0
0.00 0 7 0
Pnma
<-- Space group symbol
!Atom Typ
X
Y
Z Biso Occ /Line below:Codes
Pb PB 0.18748 0.25000 0.16721 1.40433 0.50000 0 0 0
171.00 0.00 181.00 281.00 0.00
S S
0.06544 0.25000 0.68326 0.41383 0.50000 0 0 0
191.00 0.00 201.00 291.00 0.00
O1 O
0.90775 0.25000 0.59527 1.97333 0.50000 0 0 0
211.00 0.00 221.00 301.00 0.00
O2 O
0.19377 0.25000 0.54326 1.48108 0.50000 0 0 0
231.00 0.00 241.00 311.00 0.00
O3 O
0.08102 0.02713 0.80900 1.31875 1.00000 0 0 0
251.00 261.00 271.00 321.00 0.00
! Scale
Shape1 Bov Str1 Str2 Str3 Strain-Model
1.4748
0.0000 0.0000 0.0000 0.0000 0.0000
0
11.00000 0.00 0.00 0.00 0.00 0.00
! U
V
W
X
Y
GauSiz LorSiz Size-Model
0.15485 -0.46285 0.42391 0.00000 0.08979 0.00000 0.00000 0
121.00 131.00 141.00 0.00 151.00 0.00 0.00
! a
b
c
alpha beta
gamma
8.480125 5.397597 6.959482 90.000000 90.000000 90.000000
91.00000 101.00000 111.00000 0.00000 0.00000 0.00000
! Pref1 Pref2 Asy1 Asy2 Asy3 Asy4
0.00000 0.00000 0.28133 0.03679-0.09981 0.00000
0.00 0.00 161.00 331.00 341.00 0.00
19
Quantitative phase analysis
Volume fraction
SVe2
V
2
SV
e
p
p
Weight fraction
SZMVe2
m
2
SZMV
e
p
p
Use the correct occupancies : N = occupancy / max # of Wyckoff positions
20
Tips and tricks (on the course of the refinement)
Instrumental parameters
Structure parameters
Scale factor (always)
Background (1)
Line broadening and shape (3)
Zero shift (4)
Sample displacement or
transparency (5)
Preferred orientation (7)
Surface absorption (7)
Extinction (7)
Scale factor (always)
Lattice parameter (2)
Atomic co-ordinates (6)
Temperature factors (8)
Occupancies (8), N = occ/max(N)
important for quantitative phase
analysis
21
Tips and tricks (how to obtain reliable data)
Use only good adjusted diffractometer
Bad adjustment causes the line shift and broadening; the latter cannot
be corrected in the Rietveld programs
Use only fine powders
Coarse powder “randomises” the integral intensities
Coarse powder causes problems with rough surface
Use sufficient counting time
The error in intensity is proportional to sqrt(N) as for the Poisson
distribution
Apply dead-time correction
For strong diffraction lines, the use of the dead-time correction is
strongly recommended
22
Effect of the grain size
Variations in observed intensities
(bad statistics)
Figure: Effect of specimen rotation
and particle size on Si powder
intensity using conventional
diffractometer and CuKa radiation.
International Tables for
Crystallography, Vol. C,
ed. A.J.C. Wilson,
Kluwer Academic Publishers, 1992.
23
In the Rietveld refinement don’t
refine parameters which are fixed by the structure relations
(fractional co-ordinates, lattice parameters)
refine all three parameters describing the line broadening
concurrently
refine the anisotropic temperature factors from X-ray powder
diffraction data
use diffraction patterns measured in a narrow range
forget that the number of structure parameters being refined
cannot be larger than the number of lines
24
Corundum
20-140
20-120
20-100
20-80
Reference
a[Å]
4.7584(2)
4.7585(2)
4.7587(2)
4.7591(3)
4.7586(1)
c[Å]
12.9895(3)
12.9899(4)
12.9905(4)
12.9914(5)
12.9897(1)
z(Al)
0.3521(1)
0.3521(1)
0.3520(1)
0.3521(1)
0.35216
B(Al) [Å ]
0.41(1)
0.38(2)
0.51(3)
0.59(4)
x(O)
0.3060(2)
0.3062(2)
0.3061(3)
0.3063(3)
B(O) [Å ]
0.43(3)
0.34(3)
0.36(4)
0.33(5)
ZER
0.055(9)
0.075(9)
0.090(12)
0.117(25)
DSP
-0.091(8)
-0.110(8)
-0.124(11)
-0.150(23)
TSP
0.003(1)
0.002(2)
0.003(3)
0.002(5)
R(Bragg)
3.8%
3.6%
3.2%
2.7%
2
2
0.30624
0.6%
25
Auxiliary methods
and computer programs
The most critical parameters for the convergence of the Rietveld
refinement - lattice parameters
FIRESTAR
1
2 sin
2
d hkl
sin
i i 2
1
h 2 a *2 k 2b *2 2c *2 2hka* b * cos *
2
d hkl
2
2kb * c * cosa * 2hc * a * cos *
2
1
2
d hkl
min
Only the crystal system must be known (not the space group)
The diffraction pattern must be indexed
26
Problems with
positions of diffraction lines
Residual stresses in bulk materials
Anisotropic deformation of crystallites (anisotropy of mechanical
properties)
Presence of errors in the structure (stacking faults, …)
Use of the programs working with net integral intensities
(POWOW, POWLS) is recommended
How to get the net intensities?
Numerical integration (not for the overlapped lines)
Profile fitting using analytical functions (for overlapped lines) - DIFPATAN
27
Indexing of the diffraction pattern
in unknown phases
Computer program TREOR (Trials and Errors)
Requirements
A single phase in the specimen
High-quality data (particularly, the error in the positions of diffraction
lines must not exceed 0.02° in 2)
Very good alignment of the diffractometer or the use of an internal
standard (mixed to the specimen)
28