Indexing X Ray Diffraction Patterns
Indexing X Ray Diffraction Patterns
Indexing X Ray Diffraction Patterns
LEARNING OBJECTIVES Upon completion of this module you will be able to index an X-ray diffraction pattern, identify the Bravais lattice, and calculate the lattice parameters for crystalline materials.
BACKGROUND We need to know about crystal structures because structure, to a large extent, determines properties. X-ray diffraction (XRD) is one of a number of experimental tools that are used to identify the structures of crystalline solids. The XRD patterns, the product of an XRD experiment, are somewhat like fingerprints in that they are unique to the material that is being examined. The information in an XRD pattern is a direct result of two things: (1) The size and shape of the unit cells determine the relative positions of the diffraction peaks; (2) Atomic positions within the unit cell determine the relative intensities of the diffraction peaks (remember the structure factor?). Taking these things into account, we can calculate the size and shape of a unit cell from the positions of the XRD peaks and we can determine the positions of the atoms in the unit cell from the intensities of the diffraction peaks. Full identification of crystal structures is a multi-step process that consists of: (1) Calculation of the size and shape of the unit cell from the XRD peak positions; (2) Computation of the number of atoms/unit cell from the size and shape of the cell, chemical composition, and measured density; (3) Determination of atom positions from the relative intensities of the XRD peaks We will only concern ourselves with step (1), calculation of the size and shape of the unit cell from XRD peak positions. We loosely refer to this as indexing. The laboratory module is broken down into two sections. The first addresses how to index patterns from cubic materials. The second addresses how to index patterns from non-cubic materials.
90 80 70 2,2,0 (45.30,69.2)
60
50 40 30 3,3,1 4,0,0 (72.80,13.7) (65.99,9.7) 4,2,2 (83.66,17.4) 5,1,1 3,3,3 (90.05,9.4) 3,1,1 (53.68,39.7)
20
10 0 0 5 10 15 20 25 30 35 40 45 50 55 60
6,2,0 (118.86,9.6) 2 ()
65
70
75
80
85
90
95
100
105
110
115
120
This is an example of a properly indexed diffraction pattern. All peaks are accounted for. One now needs only to assign the correct Bravais lattice and to calculate lattice parameters. How to we correctly index a pattern? The correct procedures follow.
PROCEDURE FOR INDEXING AN XRD PATTERN The procedures are standard. They work for any crystal structure regardless of whether the material is a metal, a ceramic, a semiconductor, a zeolite, etc There are two methods of analysis. You will do both. One I will refer to as the mathematical method. The second is known as the analytical method. The details are covered below. Mathematical Method Interplanar spacings in cubic crystals can be written in terms of lattice parameters using the plane spacing equation:
1 h2 + k 2 + l 2 = d2 a2 You should recall Braggs law ( = 2d sin ), which can be re-written either as:
2 = 4d 2 sin 2 OR sin 2 =
2
4d 2
Combining this relationship with the plane spacing equation gives us a new relationship: 1 h 2 + k 2 + l 2 4 sin 2 , = = 2 d2 a2 which can be rearranged to:
2 sin 2 = 2 ( h 2 + k 2 + l 2 ) 4a 2 The term in parentheses 2 is constant for any one pattern (because the X-ray 4a wavelength and the lattice parameters a do not change). Thus sin 2 is proportional to h 2 + k 2 + l 2 . This proportionality shows that planes with higher Miller indices will diffract at higher values of . 2 Since 2 is constant for any pattern, we can write the following relationship for any 4a two different planes:
2 sin 2 1 4a 2 = sin 2 2 2 2 4a 2 2 2 2 2 2 ( h1 + k1 + l1 ) sin 2 1 ( h1 + k1 + l1 ) = . or 2 sin 2 2 ( h22 + k22 + l2 ) 2 2 2 ( h2 + k2 + l2 )
The ratio of sin 2 values scales with the ratio of h 2 + k 2 + l 2 values. In cubic systems, the first XRD peak in the XRD pattern will be due to diffraction from planes with the lowest Miller indices, which interestingly enough are the close packed planes (i.e.: simple cubic, (100), h 2 + k 2 + l 2 =1; body-centered cubic, (110) h 2 + k 2 + l 2 =2; and face-centered cubic, (111) h 2 + k 2 + l 2 =3).
Since h, k, and l are always integers, we can obtain h 2 + k 2 + l 2 values by dividing the sin 2 values for the different XRD peaks with the minimum one in the pattern (i.e., the sin 2 value from the first XRD peak) and multiplying that ratio by the proper integer (either 1, 2 or 3). This should yield a list of integers that represent the various h 2 + k 2 + l 2 values. You can identify the correct Bravais lattice by recognizing the sequence of allowed reflections for cubic lattices (i.e., the sequence of allowed peaks written in terms of the quadratic form of the Miller indices). Primitive Body-centered Face-centered Diamond cubic
h2 h2 h2 h2 + k2 + k2 + k2 + k2 + l2 + l2 + l2 + l2
2 sin 2 = 2 ( h 2 + k 2 + l 2 ) 4a
which can be re-written as: a2 =
2
4sin
2
(h
+ k 2 + l 2 ) OR a =
2sin
h2 + k 2 + l 2
Worked Example Consider the following XRD pattern for Aluminum, which was collected using CuK radiation.
Intensity (%) 100 90 80 70 60 50 40 30 20 10 0 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115
(38.43,100.0)
ALUMINIUM =1.540562
(44.67,46.9)
(65.02,26.4)
(78.13,27.9)
(82.33,7.8)
Index this pattern and determine the lattice parameters. Steps: (1) Identify the peaks. (2) Determine sin 2 . (3) Calculate the ratio sin 2 / sin 2 min and multiply by the appropriate integers. (4) Select the result from (3) that yields h 2 + k 2 + l 2 as an integer. (5) Compare results with the sequences of h 2 + k 2 + l 2 values to identify the Bravais lattice. (6) Calculate lattice parameters.
Here we go!
(1) Identify the peaks and their proper 2 values. Eight peaks for this pattern. Note: most patterns will contain 1 and 2 peaks at higher angles. It is common to neglect 2 peaks. sin 2 sin 2 1 2 2 2 2 sin min sin min 3 sin sin 2 min
Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin2
h2+k2+l2
hkl
a ()
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83
sin2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859
h2+k2+l2
hkl
a ()
116.36
0.7220
(3) Calculate the ratio sin 2 / sin 2 min and multiply by the appropriate integers. 1 Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
h2+k2+l2
hkl
a ()
(4) Select the result from (3) that most closely yields h 2 + k 2 + l 2 as a series of integers. 1 Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
h2+k2+l2
hkl
a ()
(5) Compare results with the sequences of h 2 + k 2 + l 2 values to identify the miller indices for the appropriate peaks and the Bravais lattice. 1 Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
h2+k2+l2
3 4 8 11 12 16 19 20
hkl
111 200 220 311 222 400 331 420
a ()
4.0538 4.0539 4.0538 4.0538 4.0538 4.0541 4.0540 4.0541
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
h2+k2+l2
3 4 8 11 12 16 19 20
hkl
111 200 220 311 222 400 331 420
a ()
4.0538 4.0539 4.0538 4.0538 4.0538 4.0541 4.0540 4.0541
Analytical Method
This is an alternative approach that will yield the same results as the mathematical method. It will give you a nice comparison. Recall:
For any cubic system, the possible values of h 2 + k 2 + l 2 correspond to the sequence:
h 2 + k 2 + l 2 = 1,2,3,4,5,6,8,9,10,11
If we determine sin 2 for each peak and we divide the values by the integers 2,3,4,5,6,8,9,10,11, we can obtain a common quotient, which is the value of K corresponding to h 2 + k 2 + l 2 = 1. K is related to the lattice parameter as follows:
2 K = 2 OR a = 2 K 4a
If we divide the sin 2 values for each reflection by K, we get the h 2 + k 2 + l 2 values. The sequence of h 2 + k 2 + l 2 values can be used to label each XRD peak and to identify the Bravais lattice. Lets do an example for the Aluminum pattern presented above. Steps: (1) (2) (3) (4)
Identify the peaks. Determine sin 2 . Calculate the ratio sin 2 /(integers) Identify the lowest common quotient from (3) and identify the integers to which it corresponds. Let the lowest common quotient be K. (5) Divide sin 2 by K for each peak. This will give you a list of integers corresponding to h 2 + k 2 + l 2 . (6) Select the appropriate pattern of h 2 + k 2 + l 2 values and identify the Bravais lattice. (7) Calculate lattice parameters.
Here we go again!
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin 2
sin 2 2
sin 2 3
sin 2 4
sin 2 5
sin 2 6
sin 2 8
2
38.43 44.67 65.02 78.13 82.33 98.93
sin 2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776
sin 2 2
sin 2 3
sin 2 4
sin 2 5
sin 2 6
sin 2 8
7 8
111.83 116.36
0.6859 0.7220
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin 2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
sin 2 2
0.0542 0.0722 0.1444 0.1986 0.2166 0.2888 0.3430 0.3610
sin 2 3
0.0361 0.0481 0.0963 0.1324 0.1444 0.1925 0.2286 0.2407
sin 2 4
0.0271 0.0361 0.0722 0.0993 0.1083 0.1444 0.1715 0.1805
sin 2 5
0.0217 0.0289 0.0578 0.0794 0.0867 0.1155 0.1372 0.1444
sin 2 6
0.0181 0.0241 0.0481 0.0662 0.0722 0.0963 0.1143 0.1203
sin 2 8
0.0135 0.0181 0.0361 0.0496 0.0542 0.0722 0.0857 0.0903
(4) Identify the lowest common quotient from (3) and identify the integers to which it corresponds. Let the lowest common quotient be K.
Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin 2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
sin 2 2
0.0542 0.0722 0.1444 0.1986 0.2166 0.2888 0.3430 0.3610
sin 2 3
0.0361 0.0481 0.0963 0.1324 0.1444 0.1925 0.2286 0.2407
sin 2 4
0.0271 0.0361 0.0722 0.0993 0.1083 0.1444 0.1715 0.1805
sin 2 5
0.0217 0.0289 0.0578 0.0794 0.0867 0.1155 0.1372 0.1444
sin 2 6
0.0181 0.0241 0.0481 0.0662 0.0722 0.0963 0.1143 0.1203
sin 2 8
0.0135 0.0181 0.0361 0.0496 0.0542 0.0722 0.0857 0.0903
K = 0.0361 (5) Divide sin 2 by K for each peak. This will give you a list of integers corresponding to h 2 + k 2 + l 2 .
Peak No.
1 2 3 4 5 6 7
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83
sin 2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859
sin 2 K
3.000 4.000 8.001 11.001 12.002 16.000 19.001
h2 + k 2 + l 2
hkl
116.36
0.7220
20.000
(6) Select the appropriate pattern of h 2 + k 2 + l 2 values and identify the Bravais lattice.
Peak No.
1 2 3 4 5 6 7 8
2
38.43 44.67 65.02 78.13 82.33 98.93 111.83 116.36
sin 2
0.1083 0.1444 0.2888 0.3972 0.4333 0.5776 0.6859 0.7220
sin 2 K
3.000 4.000 8.001 11.001 12.002 16.000 19.001 20.000
h2 + k 2 + l 2
3 4 8 11 12 16 19 20
hkl
111 200 220 311 222 400 331 420
Sequence suggests a Face-Centered Cubic Bravais Lattice (7) Calculate lattice parameters.
a=
2 K
1.540562 A 2 0.0361
= 4.0541
These methods will work for any cubic material. This means metals, ceramics, ionic crystals, minerals, intermetallics, semiconductors, etc
The graphical method developed by Hull and Davey1 are convenient for indexing diffraction patterns, in particular for systems of lower symmetry. The reason is that this method allows one to determine structure even if lattice parameters are unknown. The mathematical methods that will be illustrated in later sections of this module require such knowledge, in particular the values of the various lattice parameter ratios (c/a, b/a, c/b etc). The steps involved in constructing and indexing patterns using Hull-Davey charts is very straightforward. First, consider the plane spacing equations for the crystal structures of interest. Some examples are shown below: Hexagonal Tetragonal Orthorhombic Etc. You should recall Braggs law ( = 2d sin ), which can be re-written either as:
1 4 h 2 + hk + k 2 l 2 = + 2 d2 3 a2 c 1 h2 + k 2 l 2 = + 2 d2 a2 c 2 2 1 h k l2 = + + d 2 a2 b2 c2
2 = 4d 2 sin 2
or sin 2 =
2
4d 2
Combining Braggs law with the plane spacing equations yields the relationship:
Hexagonal Tetragonal Orthorhombic Etc which can be rearranged in terms of sin2 to:
1 4 h 2 + hk + k 2 = + d2 3 a2 2 2 1 h +k l2 = + 2 d2 a2 c 2 2 1 h k l2 = 2+ 2 + 2 d2 a b c
l 2 4sin 2 = 2 c2
= = 4 sin 2
2
4 sin 2
A.W. Hull and W.P. Davey, Phys. Rev., vol. 17, pp. 549, 1921; W.P. Davey, Gen. Elec. Rev., vol. 25, pp. 564, 1922.
2 You should note that as unlike cubic systems where 2 is constant, your results for 4a non-cubic systems will depend upon ratios of lattice parameters (i.e., c/a, b/a, etc.) and your interaxial angles (i.e., , , ). We will illustrate this (sort of) below. This is due to the non-equivalence of indices in these systems (e.g., tetragonal 001 100; orthorhombic 001 010 100; etc).
Lets concentrate on hexagonal systems for the time being. I may ask you to derive relationships for tetragonal and orthorhombic systems in a homework assignment. As noted previously, the mathematical method requires knowledge of the c/a ratio. We dont know what it is so we need to construct a Hull-Davey chart. To accomplish this goal, we must first rewrite our revised d-spacing equations as follows:
0.1
1.0
102
103
110 200
201
2.00
1.50
1.00
10.0 0.50
c/a ratio
To determine the c/a ratio, one only needs to collect an XRD pattern, identify the peak locations in terms of the Bragg angle, calculate the d-spacing for each peak and to construct a single range d-spacing scale (2log d) that is the same size as the logarithmic [s] scale (you can use sin2 instead if you prefer). I know this is confusing so I have schematically illustrated what I mean in the next set of figures.
h1k1l1
1.0
h2k2l2 h3k3l3
1.0
1.0
log [s]
sin2-scale
d scale
h4k4l4
+ + c/a ratio
Next, you need to calculate the d-spacing or sin2 values for the observed peaks and mark them on a strip laid along side the appropriate d- or sin2 - scale.
h1k1l1
1.0
h2k2l2 h3k3l3
1.0
1.0
log [s]
sin2-scale
d scale
h4k4l4
+ + c/a ratio
The strip should be moved horizontally and vertically across the log [s] c/a plot until a position is found where each mark on your strip coincides with a line on the chart. This is illustrated schematically on the next figure.
Please keep in mind that my illustrations for the Hull-Davey method are SCHEMATIC. This method is very difficult to convey. You should consult the classical references to find out more information about this technique.
0.10
h1k1l1
10.0
1.0
h2k2l2 h3k3l3
1.0
log [s]
d scale
h4k4l4
+
10.0 1.0
c/a ratio
This is our c/a ratio for the pattern!
This method really does work as I showed you in class. Once you know your c/a ratio, you can index the XRD pattern. As we noted above, there are two ways to do this. The first is the mathematical method.
2 4 l2 sin 2 = 2 ( h 2 + hk + k 2 ) + (c / a ) 2 4a 3
Note that the lattice parameter a and the ratio of lattice parameters c/a are constant for a 2 given diffraction pattern. Thus, 2 is constant for any pattern. The pattern can now be 4a indexed in by considering the terms in brackets:
4 2 ( h + hk + k 2 ) 3
l2 (c / a ) 2 Lets start with term 1. This term only depends on the indices h and k. Thus its value can be calculated for different values of h and k. This is done below for various hk values.
The second term can be determined by substituting in the known c/a ratio. illustrated for zinc (c/a = 1.8563) in the table below. Term 2 calculated for zinc (c/a = 1.8563)
l 0 1 2 3 4 5 6 l 0 1 4 9 16 25 36
2 2 2
This is
The next step is to add the values for the two terms that are permitted by the structure factor (i.e., the values corresponding to the allowed hkl values) and to rank them in increasing order. The structure factor calculation for hexagonal systems yields the following rules: 1. When h + 2k = 3N (where N is an integer), there is no peak. 2. When l is odd, there is no peak. Both criteria must be met!
Indices (hkl) l h + 2k Peak 301 Odd 3 NO 103 Odd 1 3N YES Etc Compare with values in appendix 9 in Cullity.
Several values for the bracketed quantity are calculated below minus the peaks forbidden by the structure factor.
h 0 1 1 1 1 1 0 1 2 2 1 2
k 0 0 0 0 0 1 0 1 0 0 0 0
l 2 0 1 2 3 0 4 2 0 1 4 2
sum 1.1608 1.3333 1.6235 2.4942 3.9452 4.0000 4.6433 5.1608 5.3333 5.6235 5.9766 6.4942
The values in this table have been calculated for specific (hkl) planes. We can assign specific hkl values for each of the peaks in a hexagonal unknown by noting that the sequence of peaks will be the same as indicated in the table. Lattice parameters can be determined in two ways: We can calculate a by looking for peaks where l = 0 (i.e., hk0 peaks). If you substitute l = 0 into:
2 4 l2 sin 2 = 2 ( h 2 + hk + k 2 ) + (c / a ) 2 4a 3
you will get,
2 4 sin 2 = 2 ( h 2 + hk + k 2 ) 4a 3
OR a=
3 sin
h 2 + hk + k 2
You can now perform this calculation for every hk0 peak, which will yield values for a. Similarly, values for c can be determined by looking for 00l type peaks. In these instances, h = k = 0. Thus,
2 4 2 l2 2 sin = 2 ( h + hk + k ) + (c / a ) 2 4a 3
2
becomes
c=
2sin
Worked Example Consider the following XRD pattern for Titanium, which was collected using CuK radiation.
Titanium Pow der (-325 m esh)
Intensity (counts)
20
30
40
50
60
70
80
90
100
Two Theta
Index this pattern and determine the lattice parameters. Steps: (1) Identify the peaks. (2) Determine values of
4 2 ( h + hk + k 2 ) for reflections allowed by the structure factor. 3 l2 (3) Determine values of for the allowed reflections and the known c/a ratio (c / a ) 2 (4) Add the solutions from parts (2) and (3) together and re-arrange them in increasing order.
(5) Use this order to assign indices to the peaks in your diffraction pattern. (6) Look for hk0 type reflections and calculate a for these reflections. (7) Look for 00l type reflections. Calculate c for these reflections.
Here we go!
0 1 2 3
l2 for the allowed reflections and the known c/a ratio. (c / a ) 2 Titanium: c/a = 1.587
l
0 1 2 3 4 5 6
l2
0 1 4 9 16 25 36
(4) Add the solutions from parts (2) and (3) together and re-arrange them in increasing order.
hkl 002 100 101 102 103 110 004 112 200 201 104 202 203 105 114 210 211 204 006 212 106 213 300 205 302
Pt.1+Pt.2 1.588 1.333 1.730 2.921 4.906 4.000 6.352 5.588 5.333 5.730 7.685 6.921 8.906 11.258 10.352 9.333 9.730 11.685 14.292 10.921 15.625 12.906 12.000 15.258 13.588
hkl 100 002 101 102 110 103 200 112 201 004 202 104 203 210 211 114 212 105 204 300 213 302 006 205 106
Pt.1+Pt.2 1.333 1.588 1.730 2.921 4.000 4.906 5.333 5.588 5.730 6.352 6.921 7.685 8.906 9.333 9.730 10.352 10.921 11.258 11.685 12.000 12.906 13.588 14.292 15.258 15.625
(5) Use this order to assign indices to the peaks in your diffraction pattern.
Peak 35.275 38.545 40.320 53.115 63.095 70.765 74.250 76.365 77.500 82.360 86.940 92.900 I/Io 21 18 100 16 11 9 10 8 14 2 2 10 d (nm) sin2 0.091805 2.5423 0.108941 2.3338 0.118779 2.2351 0.199895 1.7229 0.273744 1.4723 0.335278 1.3303 0.36428 1.2763 0.382132 1.2461 0.39178 1.2307 0.433526 1.1699 0.473309 1.1197 0.525296 1.0628
hkl 100 002 101 102 110 103 200 112 201 004 202 104 AVG a c h2+hk+k2 l2
(6) Look for hk0 type reflections and calculate a for these reflections.
Peak 35.275 38.545 40.320 53.115 63.095 70.765 74.250 76.365 77.500 82.360 86.940 92.900 I/Io 21 18 100 16 11 9 10 8 14 2 2 10 d (nm) sin2 0.091805 2.5423 0.108941 2.3338 0.118779 2.2351 0.199895 1.7229 0.273744 1.4723 0.335278 1.3303 0.36428 1.2763 0.382132 1.2461 0.39178 1.2307 0.433526 1.1699 0.473309 1.1197 0.525296 1.0628
hkl 100 002 101 102 110 103 200 112 201 004 202 104 AVG a 2.936 c h2+hk+k2 1 l2
2.945 2.947
3 4
2.943
c/a:
(7) Look for 00l type reflections. Calculate c for these reflections.
Peak 35.275 38.545 40.320 53.115 63.095 70.765 74.250 76.365 77.500 82.360 86.940 92.900 I/Io 21 18 100 16 11 9 10 8 14 2 2 10 d (nm) sin2 0.091805 2.5423 0.108941 2.3338 0.118779 2.2351 0.199895 1.7229 0.273744 1.4723 0.335278 1.3303 0.36428 1.2763 0.382132 1.2461 0.39178 1.2307 0.433526 1.1699 0.473309 1.1197 0.525296 1.0628
hkl 100 002 101 102 110 103 200 112 201 004 202 104 AVG a 2.936 c h2+hk+k2 1 l2
4.668
2.945 2.947
3 4
4.680
16
2.943
4.674
c/a:
1.588
Pretty good correlation with ICDD value. Actual c/a for Titanium is 1.5871
This method, though effective for most powder XRD data, can yield the wrong results if XRD peaks are missing from your XRD pattern. In other words, missing peaks can cause you to assign the wrong hkl values to a peak. Other methods should be available.
Analytical Method for Non-Cubic Crystals
To accurately apply this technique, one must first consider our altered plane spacing equation:
2 4 l2 sin 2 = 2 ( h 2 + hk + k 2 ) + (c / a ) 2 4a 3
Since a and c/a are constants for any given pattern, we can re-arrange this equation to:
sin 2 = A ( h2 + hk + k 2 ) + Cl 2
where A = . Since h, k, and l are always integers, the term in 3a 2 4c 2 parentheses, h 2 + hk + k 2 can only have values like 0, 1, 3, 4, 7, 9, 12 and l2 can only have values like 0, 1, 4, 9,. We need to calculate sin 2 for each peak, divide each sin 2 value by the integers 3, 4, 7, 9 and look for the common quotient (i.e., the sin 2 / n value that is equal to one of the observed sin 2 values). The sin 2 values representing this common quotient refer to hk0 type peaks. Thus this common quotient can be tentatively assigned as A. We can now re-arrange terms in our modified equation to obtain C. This is done as follows:
and C =
sin 2 = A ( h 2 + hk + k 2 ) + Cl 2 Cl 2 = sin 2 A ( h 2 + hk + k 2 )
We get the value of C by subtracting from each sin 2 the values of nA (i.e., A, 3A, 4A, 7A,) where A is the common quotient that we identified above. Next, we need to look for the remainders that are in the ratio of 1, 4, 9, 16, which will be peaks of the 00l type. We can determine C from these peaks. The remaining peaks are neither hk0-type nor 00l type. Instead they are hkl-type. They can be indexed from a combination of A and C values. Lets do an example.
Worked Example Steps: (1) (2) (3) (4) (5) (6) Identify the peaks and calculate sin 2 for each peak. Divide each sin 2 value by the integers 3, 4, 7, 9. Look for the common quotient. Let the lowest common quotient represent A. Assign hk0 type indices to peaks. Calculate sin 2 - nA where n = 1, 3, 4, 7.
(7) Look for the lowest common quotient. From this we can identify 00l type peaks. Recall, that 001 is not allowed for hexagonal systems. The first 00l type peak will be 002. We can calculate C from: C l 2 = sin 2 A (h 2 + hk + k 2 ) (8) Look for values of sin 2 that increase by factors of 4, 9 (this is because l = 1, 2, 3 and l2 = 1, 4, 9). Peaks exhibiting these characteristics are 00l type peaks, which can be assigned the indices 004, 009, etc). Also note that the values of sin 2 will be some integral number times the value observed in (7) which indicates the indices of the peak (9) Peaks that are neither hk0 nor 00l can be identified using combinations of our calculated A and C values. (10) Calculate the lattice parameters from the values of A and C.
Confused yet? You could be. I was the first time I learned these things. Let me show you an example that should make all things clear.
Here we go!
Consider the diffraction pattern for Titanium as shown below. This one is a little different than the specimen that we analyzed above.
Intensity (%) 100 (40.18,100.0)
90
80
70 60
50 40
30 20
(35.10,25.5) (76.24,13.1) (102.41,4.4) (105.82,1.4) (53.02,12.8) (62.96,13.4) (70.70,13.0) (77.37,9.3) (86.79,2.1) (109.08,8.3) (114.35,5.4) (119.30,2.7) 2 () 100 105 110 115 120
Steps to success: 1. Calculate sin2 for each peak 2. Divide each sin2 value by integers 3, 4, 7 (from h2+hk+k2 allowed by the structure factor) 3. Look for lowest common quotient. 4. Let lowest common quotient = A . 5. Peaks with lowest common quotient are hk 0 type peaks. Assign allowed hk 0 indices to peaks.
Peak 35.100 38.390 40.170 53.000 62.940 70.650 74.170 76.210 77.350 82.200 86.740 92.680 102.350 105.600 109.050 114.220 119.280
I/Io 25.5 25.6 100 12.8 13.4 13 1.8 13.1 9.3 1.7 2.1 1.8 4.4 1.4 8.3 5.4 2.7
(sin2)/3 (sin2)/4 (sin2)/7 (sin2)/9 (sin2)/12 sin2 0.0909 0.0303 0.0227 0.0130 0.0101 0.0076 0.1081 0.0360 0.0270 0.0154 0.0120 0.0090 0.1179 0.0393 0.0295 0.0168 0.0131 0.0098 0.1991 0.0664 0.0498 0.0284 0.0221 0.0166 0.0908 0.2725 0.0681 0.0389 0.0303 0.0227 0.3343 0.1114 0.0836 0.0478 0.0371 0.0279 0.3636 0.1212 0.0909 0.0519 0.0404 0.0303 0.3808 0.1269 0.0952 0.0544 0.0423 0.0317 0.3905 0.1302 0.0976 0.0558 0.0434 0.0325 0.4321 0.1440 0.1080 0.0617 0.0480 0.0360 0.4716 0.1572 0.1179 0.0674 0.0524 0.0393 0.5234 0.1745 0.1308 0.0748 0.0582 0.0436 0.6069 0.2023 0.1517 0.0867 0.0674 0.0506 0.6345 0.2115 0.1586 0.0906 0.0705 0.0529 0.6632 0.2211 0.1658 0.0947 0.0737 0.0553 0.7051 0.2350 0.1763 0.1007 0.0783 0.0588 0.7445 0.2482 0.1861 0.1064 0.0827 0.0620 A= 0.0908
=7xA
210
(sin2)/LCQ 1.0 1.2 1.3 2.2 3.0 3.7 4.0 4.2 4.3 4.8 5.2 5.8 6.7 7.0 7.3 7.8 8.2
6. Subtract from each sin2 value 3A , 4A , 7A (from h 2+hk +k 2 allowed by the structure factor) 7. Look for lowest common quotient (LCQ). From this you can identify 00l -type peaks. The first allowed peak for hexagonal systems is 002. Determine C from the equation: C l 2 = sin2-A (h 2+hk +k 2) since h =0 and k =0, then: C=LCQ/l 2 = sin2/l 2 8. Look for values of sin2 that increase by factors of 4, 9, 16... (because l = 1,2,3,4..., l 2=1,4,9,16...) The peaks exhibiting these characteristics are 00l -type peaks (002...).
We identify the 4th peak as 102 because we observe the 2 LCQ for sin -1A. Recall that the 1 comes from the 2 2 quadratic form of Miller indices (i.e., h +hk +k =1). We identify the 8th peak as 112 because we observe the 2 LCQ for sin -3A. Recall that the 1 comes from the 2 2 quadratic form of Miller indices (i.e., h +hk +k =3). We identify the 11th peak as ... etc...
1.54062 Peak 35.100 38.390 40.170 53.000 62.940 70.650 74.170 76.210 77.350 82.200 86.740 92.680 102.350 105.600 109.050 114.220 119.280
I/Io 25.5 25.6 100 12.8 13.4 13 1.8 13.1 9.3 1.7 2.1 1.8 4.4 1.4 8.3 5.4 2.7
sin2-A sin2-3A sin2-4A sin2-4A sin2 0.0909 0.1081 0.0173 0.1179 0.0271 0.1083 0.1991 0.2725 0.1817 0.0001 0.3343 0.2435 0.0618 0.3636 0.2728 0.0911 0.0003 0.3808 0.2900 0.1083 0.0175 0.3905 0.2997 0.1180 0.0272 0.4321 0.3413 0.1597 0.0688 0.4716 0.3807 0.1991 0.1083 0.5234 0.4326 0.2509 0.1601 0.6069 0.5161 0.3345 0.2436 0.6345 0.5436 0.3620 0.2711 0.6632 0.5724 0.3907 0.2999 0.0274 0.7051 0.6143 0.4326 0.3418 0.0693 0.7445 0.6537 0.4721 0.3812 0.1087
h 1 0 1 1 2 1 0 2
k 0 0 0 1 0 1 0 0
l 0 2 2 0 0 2 4 2
l =LCQ/C 4.0
LCQ = 0.1083
9. Peaks that are not hk 0 or 00l can be identified using combinations of A and C values. This is accomplished by considering: sin2 = C l 2 + A (h 2+hk +k 2) Cycle through allowed values for l and hk , and compare sin2 value to labeled peaks. 10. Once A and C are known, the lattice parameters can be calculated.
2 2 2 2 sin = C l + A (h +hk +k )
A 0.0908 Peak 35.100 38.390 40.170 53.000 62.940 70.650 74.170 76.210 77.350 82.200 86.740 92.680 102.350 105.600 109.050 114.220 119.280
C 0.0270 I/Io 25.5 25.6 100 12.8 13.4 13 1.8 13.1 9.3 1.7 2.1 1.8 4.4 1.4 8.3 5.4 2.7 sin 0.0909 0.1081 0.1179 0.1991 0.2725 0.3343 0.3636 0.3808 0.3905 0.4321 0.4716 0.5234 0.6069 0.6345 0.6632 0.7051 0.7445
2
sin
2
h 1 0 1 1 1 1 2 1 2 0 2 1 2 2 2 1 2
k 0 0 0 0 1 0 0 1 0 0 0 0 0 1 1 1 1
l 0 2 1 2 0 3 0 2 1 4 2 4 3 0 1 4 2
Calculated
0.0908 0.1081 0.1179 0.1989 0.2725 0.3341 0.3633 0.3806 0.3903 0.4324 0.4714 0.5232 0.6065 0.6358 0.6628 0.7049 0.7439
10. Once A and C are known, the lattice parameters can be calculated.
a 2.951
c 4.686
c /a 1.588
HOW DO THESE VALUES COMPARE WITH THOSE FROM THE ICDD CARDS?