Theory and practice of X-ray diffraction experiment

Theory and practice of X-ray
diffraction experiment
Power Point presentation from lecture # 2 is available at:
http://dl.dropbox.com/u/23622306/UIUC/Lecture%202.ppt
Power Point presentation from lecture # 4 is available at:
http://dl.dropbox.com/u/23622306/UIUC/lecture%204.ppt
Grain size and types of diffraction experiment
•
•
•
•
Single crystal
Bulk powder
Nanopowder
Amorphous/glass
Reciprocal space
•
•
•
•
z
•
hkl
X-ray
x
•
•
•
y
Direct space relates to atoms in the unit cell, while reciprocal
space relates to peaks in diffraction experiment.
In diffraction experiment planes of atoms in the crystal act as
a “selective” mirror and reflect (in the same way as in optics)
incoming beam of x-rays if specific geometric conditions are
satisfied.
Vectors in reciprocal space correspond to families of planes
in the crystal (the direction of the rec. vector is normal to the
family of planes).
In a conventional crystal (as opposed to incommensurately
modulated crystal or quasi-crystal) vectors in reciprocal
space can occupy only points on a 3-dimensional grid.
Grid coordinates of reciprocal vectors are known as Miller
indices hkl
Geometry of the reciprocal space grid is determined by the
unit cell of the crystal. And can be directly measured.
D-spacing is equal to inverse of the reciprocal vector length.
Coordinates of vectors in reciprocal space are described in
laboratory (instrument-related) reference coordinate system.
•
UB matrix relates the Miller indices of reciprocal vector (hkl) with
its Cartesian coordinates in lab reference system (xyz) at zero
goniometer position.
Orientation matrix
xyz=UB hkl
•
•
•
•
Columns in the orientation matrix are coordinates of the principal
vectors in reciprocal space 1 0 0 , 0 1 0 and 0 0 1.
UB matrix is composed of two sub-matrices, U and B
U describes the orientation of the crystal axes with respect to
the laboratory reference system, while B stores information
about the unit cell parameters.
By inverting the above equation one can calculate what are the
Miller indices of a measured reciprocal vector xyz
hkl=UB-1 xyz
z
hkl
X-ray
x
y
Diffraction condition and Bragg equation for a single crystal
hkl
1q
diffracted X-ray
nl=2dsinq
2q
incident X-ray
1q
2q
Diffraction geometry is analogous to reflection in a selective
mirror:
hkl
Diffraction is effective only at selected incidence angles with
respect to the “reflecting” plane, determined by the Bragg
equation.
When diffraction is effective “reflection” angle equals to the
incidence angle, with the reciprocal vector bisecting the two.
Intensity of the “reflected” beam is not equal to the intensity of
the incident beam, it is related to the mean electron density
within the direct lattice plane.
Explanation of Bragg derivation
nl=2dsinq
•
•
•
In diffraction experiment we measure the direction and intensity of the diffracted
beams as well as the orientation of the crystal corresponding to each of the
diffraction peaks.
Diffraction does not occur at any arbitrary crystal orientation. In general to bring a
given reciprocal vector to a diffraction position one may need to rotate the sample
using goniometer or change the incident wavelength (an exception is
polychromatic experiment in which a range of incident wavelengths is available).
When the diffraction event occurs, the following relation between the incident and
diffracted beam vectors, incident wavelength and the reciprocal vector is satisfied:
Ewald construction
1/l(S0+Sd)=R xyz
•
•
R is the goniometer rotation matrix bringing the crystal from the zero-goniometer
position to the position at which maximum peak intensity is observed.
In the Ewald construction the center of the crystal and the center of reciprocal
space do not coincide.
S0
Sd
Rxyz
Radius of the Evald sphere is 1/l
Experimental approach to
single crystal diffraction
Effective diffraction vs. observed diffraction
•
Cause effective diffraction and find (observe) diffraction signal
•
Measure Sd and R for a number of diffraction peaks
•
Calculate xyz using the Ewald equation
•
Find orientation matrix/index the peaks (assign Miller indices)
•
Measure/calculate peak intensities I(hkl)
Wide oscillation ± 13°
Monochromatic SXD experiment
Goniometry in matrix representation
Atomic scattering factor
•
The incoming X-rays are scattered by the electrons of the atoms. As the wavelength of the X-rays (1.5
to 0.5 A) is of the order of the atom diameter, most of the scattering is in the forward direction. For
neutrons of the same wavelength the scattering factor is not angle dependent due to the fact that the
atomic nucleus is magnitudes smaller than the electron cloud. It is also obvious that the X-ray
scattering power will depend on the number of electrons in the particular atom. The X-ray scattering
power of an atom decreases with increasing scattering angle and is higher for heavier atoms. A plot of
scattering factor f in units of electrons vs. sin(theta)/lambda shows this behavior. Note that for zero
scattering angle the value of f equals the number of electrons.
Electron density around an isolated atom (r) is spherically
symmetric, with max at nucleus position, and falls off smoothly
with distance.
f s  

  r e
2i r s 
dr  f s 

Cromer and Mann 9-parameter equation
Anomalous dispersion
The scattering factor contains additional (complex) contributions from anomalous dispersion effects
(essentially resonance absorption) which become substantial in the vicinity of the X-ray absorption
edge of the scattering atom. These anomalous contributions can be calculated as well and their
presence can be exploited in the MAD phasing technique.
Excitation Scans
We can observe the Δf” by measuring
the absorption of the x-rays by the
atom. Often we us the fluorescence of
the absorbing atom as a measure of
absorptivity. That is, we measure an
“excitation” spectrum.
How to get Δf ’?
The real, “dispersive” component is
calculated from Δf” by the KramersKronig relationship. Very roughly, Δf’ is
the negative first derivative of Δf”.
From Ramakrishnan’s study of GH5
Thermal vibrations and Debye-Waller factor
•
There is an additional weakening of the scattering power of the atoms by the so
called Temperature-, B- , or Debye-Waller factor. This exponential factor is also
angle dependent and effects the high angle reflections substantially (one of the
reasons for cryo-cooling crystals is to reduce the attenuation of the high angle
reflections due to this B-factor).
Structure factor
1
 r  
Vc
F e
 2i h r 
h
h
Fh    r e
V
N
Fh   f i sh e
i 1
I h  Fh Fh
*
2i hx i 
2i hr 
dr
Lorenz correction
Accounts for the different speed with which the reciprocal vectors move
through the Ewald sphere
L
cos q
sin 2 2q
Polarization correction
105
100
95
I
90
85
80
75
70
-200
-150
-100
-50
0
azimuth (deg)
50
100
150
200
Other intensity corrections
•
•
•
•
•
Absorption correction
Extinction correction
Preferred orientation
Illuminated volume
Incident intensity
Symmetry of peak intensities. Friedel pairs and Laue classes
Friedel's law, named after Georges Friedel, is a property of Fourier
transforms of real functions.1
Given a real function , its Fourier transform
Laue
classes
Non-centrosymmetric
groups
having the same
Laue class
1
has the following properties.
2/m
2, m
mmm
222, 2mm
-3
-3m
where is the complex conjugate of .
Centrosymmetric points are called Friedel's pairs.
The Friedel pair symmetry is broken by anomalous dispersion. The
effect is strongly pronounced for incident energies close to the
absorption edge.
3
32, 3m
4/m
4,
4/mmm
422,
6/m
6,
6/mmm
622,
m-3
m3m
, 42m
, 62m
23
432,
32
Systematic absences and space group determination
A centered hkl
B centered
C centered
F centered
I centered
R (obverse)
R (reverse)
k + l = 2n
h + l = 2n
h + k = 2n
k + l = 2n, h + l = 2n, h + k = 2n
h + k + l = 2n
-h + k + l = 3n
h - k + l = 3n
Screw || [100]
21, 42
41, 43
Screw || [010]
21, 42
41, 43
Screw || [001]
21, 42, 63
31, 32, 62, 64
41, 43
61, 65
h00
h = 2n
h = 4n
0k0
k = 2n
k = 4n
00l
l = 2n
l = 3n
l = 4n
l = 6n
Glide reflecting in a
b glide
c glide
n glide
d glide
Glide reflecting in b
a glide
c glide
n glide
d glide
Glide reflecting in c
b glide
a glide
n glide
d glide
0kl
k = 2n
l = 2n
k + l = 2n
k + l = 4n
h0l
h = 2n
l = 2n
h + l = 2n
h + l = 4n
hk0
k = 2n
h = 2n
k + h = 2n
k + h = 4n
Point detector experiment
•
•
•
•
•
•
•
•
•
•
Scintillator detector converts x-ray to electric signal
Center the sample on rotation axis and with the beam
Search for peaks (usually 10-20) at different 2theta and sample orientations
Center the peaks that were found
Determine orientation matrix
Find more peaks to constrain and refine the orientation matrix better
Calculate position for a list of peaks that need to be measured based on orientation
matrix
Position each peak individually, record a rocking curve/peak profile
Integrate, scale and correct peak intensities
Solve/refine the structure
Area detector experiment
•
•
•
•
•
•
•
CCD detector detects visible light. Phosphor screen in front of the detector converts
x-ray to visible light.
Center the sample on rotation axis and with the beam
Collect diffraction images while rotating the sample
Determine detector coordinates and sample orientations for each diffraction peak
Reconstruct the reciprocal space in 3-d and determine the orientation matrix (index)
Predict peak positions in recorded diffraction images and retrieve peak intensities
(structure factor amplitudes)
Solve/refine the structure
Area detector and step
scan approach
With wide rotation image we determine the direction of the
diffracted beam Sd but the rotation angles (necessary to calculate
R) at which each peak occurs are unknown.
Step scan allows to determine R for each peak by finding a step
image at which the peak has maximum intensity.
With R and Sd we can calculate xyz.
Step scan images
General algorithm for
analysis of single-crystal
XRD data
zoom
fitting
Residual
Image
Initial peak list (Sd)
Step scan analysis (R)
Orientation matrix (UB)
Predicted peak list
Integrated intensities I(hkl)
Scaled and corrected intensities
Figures of merit
Rint
F  F


F
2
o
2
o
2
o
Rsigma
 F 


F
2
o
2
o
Demonstration of single-crystal data processing with
GSE_ADA
Powder diffraction
•
•
•
•
•
Sample is composed of a very large (preferably >106)
number of single crystal specimens with random
distribution of crystal orientations.
Diffraction for all reciprocal vectors occurs (is effective) at
any arbitrary orientation of the sample. Sample rotation is
not necessary.
Single crystal diffraction peaks (directional beams) become
cones of radiation.
Peaks corresponding to different reciprocal vectors
(different hkls) that have similar length (d-spacing) overlap.
Diffraction signal carries only information about the length
of the reciprocal vectors, but not their orientation.
Different experimental approaches to powder diffraction
•
Polychromatic EDX approach
–
–
–
–
–
–
–
•
Monochromatic approach with point detector
–
–
–
–
–
•
White beam
Point energy-dispersive detector at fixed angle
Non-scanning signal accumulation
Can provide access to large d-spacing range w/o much of angular access
Suffers from preferred orientation problems
Peak usually quite broad
Peak intensity interpretation difficult (correction)
Mono beam
Point detector at variable angle
Scanning signal accumulation
Suffers from preferred orientation problems
Can provide patterns with very narrow peak profiles
Monochromatic approach with area detector
–
–
–
–
Mono beam
Area detector at fixed or variable angle/position
Non-scanning signal accumulation
Does not suffer from preferred orientation problems
High resolution powder diffraction with analyzer
crystal
Bragg-Brentano powder diffractometer
Integration of 2-dimensional diffraction pattern
Calibration of detector geometry
–
–
–
–
•
Goniometer geometry calibration
–
–
•
•
•
Detector coordinates of the point of intersection of the
incident beam and the detector surface
Sample-to-detector distance
Detector non-orthogonality with respect to the beam
Pixel size
Goniometer zeros
Goniometer axis alignment with the detector orientation
Instrumental function
Calibrating with a diffraction standard
Cake transform
0.07
0.065
0.06
peak width [deg]
•
0.055
0.05
0.045
0.04
190
210
230
250
270
detector distance [mm]
290
310
Energy dispersive powder diffraction
Ge solid state semiconductor detector
Multiple wavelengths
Powder pattern fitting techniques
•
Individual peak fitting
– Peak positions are not constrained – they are individually refined
– Peak profiles are not related with each other
– Closely overlapping peaks are very hard to deal with
•
LeBail (unit-cell constrained) refinement
– Individual peak positions are not refined – they are calculated from the unit cell
parameters which are refined
– Individual peak profiles are not refined – there is a global function (e.g. Cagliotti
function) that connects all peak profiles.
– Peak intensities are refined free (not constrained by the structure model)
– Gives much more reliable way of refining positions and intensities of closely
overlapping peaks.
•
Rietveld (model-biased) refinement
– Individual peak positions are not refined – they are calculated from the unit cell
parameters which are refined
– Individual peak profiles are not refined – there is a global function (e.g. Cagliotti
function) that connects all peak profiles.
– Individual peak intensities are not refined free, they are calculated from the
structure model which is refined.
Peak width analysis
•
Instrumental function
–
–
–
•
Sample-related factors
–
–
–
•
Focused beam divergence away from the sample
Detector pixel size
Phosphor point spread function
Scherrer formula
Hall-Williamson plot
Cagliotti function
Deconvolution
size = K l / FW(S) cos(q)
FW(S) cos(q) = K l / size + 4 strain sin(q)
FW(S)D = FWHMD - FW(I)D
The dimensionless shape factor K has a typical value of about 0.9, but varies with the actual
shape of the crystallite.
Scherrer formula is not applicable to grains larger than about 0.1 μm, which precludes those
observed in most metallographic and ceramographic microstructures.
Intensity(Counts)
Peak width analysis
13nm nano powder
FW(S) cos(q) = K l / size + 4 strain sin(q)
[CNMDAC3_023.chi]
5000
4500
4000
Bulk fine powder
3500
3000
2500
2000
1500
1000
500
0
4
5
6
7
8
9
10
11
12
13
T wo-T heta (deg)
14
15
16
17
18
19
20
21
S&M and Quantitative analysis
demonstration
•
Image integration
–
–
•
Powder pattern analysis
–
–
•
American Mineralogist Crystal Structure Database
Project RRUFF
Crystallography Open Database
ICSD
PDF
Peak and pattern fitting, refinement of structure
–
–
–
–
–
•
MDI Jade
Treor (index permutation method)
Ito (zone indexing method)
Dicvol (dichotomy algorithm)
Database and S&M
–
–
–
–
–
•
Jade
GSE_Shell
Indexing
–
–
–
–
•
Fit2d
Powder3d
GSAS
Fullprof
PowderCell
Endeavour
Jade
Structure determination ab initio
–
–
–
Endeavour
Fox
Expo
Tools of the trade
Powder diffraction
Tools of the trade
Single-crystal diffraction
•
Commercial
– Bruker Saint/SMART, APEX
– HKL Research Denzo/HKL2000
– Oxford (Agilent) Crysalis
•
Academic
– XDS (Wofgang Kabsch) Linux
– GSE_ADA (Przemek Dera) Windows
Structure determination
Structure identification
“Logical” choices
Search and match
Indexing-based
Structure refinement (requires the structure model to be approximately known)
Rietveld method
Single-crystal refinement
Structure solution (ab initio)
direct methods
charge flipping
simulated annealing
High pressure apparatus
Incident x-ray beam
Incident x-ray beam
Beryllium seats
Diamond anvil
Diamond anvil
Metal
Metal
gasket
gasket
Steel frame
Steel frame
Diffracted x-ray beam
Diffracted x-ray beam
• Sample is immersed in hydrostatic
liquid, which freezes at some point
during compression (usually below 10
GPa).
• Diffraction pressure calibrant is
placed in the sample chamber along
with the sample.
• Both incident and diffracted beam
travel though diamonds, Be disks,
pressure medium, and sample.
High-pressure crystallography challenges
Experimental challenges:
–
–
–
–
–
–
–
–
–
–
Very small sample (<0.01mm)
Angular access restricted (low completeness)
Absorption limits the incident energy to >15 keV
Absorption and extinction affect intensity measurement
High background (scattering, Compton, etc.)
Multiple SXD diffraction signal
Contamination by PXD signal
Poor sample quality (strain, multi-grain assemblages)
More challenging sample centering
Beam size vs. sample size
Homework (detailed instructions will be sent by
Friday 9/16), due on Tuesday 9/27
1.
2.
3.
4.
5.
Download and install Rosetta. Feel free to use any other program
of your preference to complete the tasks described below (in
which case you can skip 1).
Download the example unknown powder patterns xx1.chi, xx2.chi,
xx3.chi
Identify minerals present in each sample
Download the example QA powder patterns of quartz-albite
mixture QA25.chi, QA50.chi. The number in file name stands for
the albite content (as prepared).
Verify how accurate is the quantitative analysis done with
Rosetta.