International
Tables for Crystallography Volume C Mathematical, physical and chemical tables Edited by E. Prince © International Union of Crystallography 2006 
International Tables for Crystallography (2006). Vol. C, ch. 5.4, pp. 537540
https://doi.org/10.1107/97809553602060000598 Chapter 5.4. Electrondiffraction methodsLatticeparameter determination by electrondiffraction methods is reviewed. Obtaining cell axes and angles from a single pattern with suitable Laue zones and two patterns with different zone axes is discussed. Methods based on Kikuchi or highorder Laue zone (HOLZ) lines are also described: these have the advantage that the camera constant does not need to be known. 
This article treats the recovery of cell axes and angles from (a) a single pattern with suitable Laue zones and (b) two patterns with different zone axes. It is assumed that instrument distortions, if significant, are corrected and that the patterns are free of artefacts such as twinning, double diffraction etc. (Edington, 1975). The treatment is valid for convergentbeam, micro and selectedarea electrondiffraction patterns and accelerating voltages above approximately 30 kV. Relevant papers are by LePage (1992) and Zuo (1993), and background reading is contained in Edington (1975), Gard (1976), and Hirsch, Howie, Nicholson, Pashley & Whelan (1965).
The basic requirement in the determination of the unit cell of a crystal is to find, from one or more diffraction patterns, the basis vector set, a*, b*, c*, of a primitive reciprocal unit cell. The Cartesian components of these vectors form an orientation matrixwhich, when inverted, gives the vector components of the corresponding realspace cell. The elements of UB can be measured directly from the diffraction pattern in millimetres. Define axes x and y to be in the recording plane and z in the beam direction. A point in the diffraction pattern x, y, z is then related to the indices h, k, l by Note that points with nonzero z are observed on the plane z = 0, see Fig. 5.4.1.1.
The metric M of UB^{−1} is used to find the unitcell edges and angles as where T means the transpose. Then, gives and where L is the effective distance between the diffracting crystal and the recording plane and λ is the wavelength. These quantities are defined in Fig. 5.4.1.1 together with the nomenclature and geometrical relationships required in this article.
If necessary, the cell is reduced to the Bravais cell according to the procedures given in IT A (2005, Chapter 9.2 ), before calculating the metric.
In practice, there may be a difficulty in choosing a vector set that describes a primitive reciprocal cell. Although a record of any reasonably dense plane of reciprocal space immediately exposes two basis vectors of a cell, the third vector lies out of the plane of the diffraction pattern containing the first two vectors and may not be directly measurable. Hence, some care must be taken to ensure that the third vector chosen makes the cell primitive. This presents no difficulty in the nonzerozone analysis given in Subsection 5.4.1.3 but needs consideration when two patterns are used as described in Subsection 5.4.1.2.
Unitcell parameters may be partly or fully determined depending on the extent of knowledge of the effective camera length L and wavelength λ, and upon the type and number of patterns. The different situations are distinguished in Table 5.4.1.1 for photographic recording. Note that L is simply a magnification factor.

The accuracy of the nonzerozone analysis described in Subsection 5.4.1.3 will depend on the influence of spiral, radial, and elliptical distortions. The camera length and wavelength need to be known to better than 1%.
For measurements made using microscope deflector systems, a knowledge of L may not be required depending on the location of the deflectors in the microscope column. Instead, the calibration factor of the deflectors, in suitable units, will replace L.
Two patterns are required that represent different sections through the reciprocal lattice. The angle of rotation p between these sections must be known, as well as the trace of the rotation axis in the plane of the pattern. Define the plane of the first pattern as the xy plane with the x axis, for convenience, coincident with the trace of the rotation axis. The coordinates of reflections in the first pattern are then measured. The coordinates of spots in the patterns rotated by p relative to the first pattern are then measured, care being taken to align the trace of the rotation axis with the x axis of the measuring equipment for each pattern. The coordinates of these reflections are then reduced to the coordinate system of the first pattern by the relations , , . The coordinates of all reflections measured are placed in a table that is scanned to extract the three shortest noncoplanar vectors. If the patterns come from dense, neighbouring zones, it is likely that these vectors define a primitive cell.
One pattern with well defined Laue zones taken at an arbitrary zone allows the recovery of a third vector not in the plane of the pattern. First find the components for the two shortest zerozone vectors a* and b* (see Fig. 5.4.1.2 ). For the third vector, measure the x and y coordinates of a nonzerozone vector. These are then transformed to components in the coordinate system defined by the two shortest vectors of the zero zone. Then the fractional part of each transformed component is taken. If the vector comes from a Laue zone of order n, these fractional parts must be divided by n. A transformation back to Cartesian coordinates gives the x and y components of the third vector c*, point O′ in Fig. 5.4.1.2. In practice, it is best to measure a set of vectors at roughly equal angles around the zone so that an average can be taken to improve accuracy.

A diffraction pattern with the crystal oriented at a zone axis. The lower diagram shows the zerozone vectors a* and b*, and OO′, the projection on the xy plane of c*. 
For certain space groups and special orientations, it is possible for half of the zerozone reflections to be absent. If sufficient nonzerozone vectors have been measured, two different vectors OO′ should be found. New axes (a* + b*)/2 and (a* − b*)/2 must then be chosen and the x and y coordinates of c* recalculated.
The z component of c* is obtained (Fig. 5.4.1.1) from where r is the radius of the Laue circle of order n.
We now have the orientation matrix UB as and the measurement of the pattern is complete.
After transformation of the axes defined by UB to the Bravais axes, the inverse of the resulting UB_{Bravais} will index the pattern using
Latticeparameter determination based on selectedarea electrondiffraction patterns requires accurate calibration of the camera constant K. This constant depends on the electron wavelength λ and the camera length L and is given by K = λL. Because the camera constant cannot be determined with sufficient accuracy in a transmission electron microscope, a number of methods for determination of lattice parameters (or electron wavelength) have been developed based on Kikuchi or HOLZ (highorder Laue zone) lines (Uyeda, Nonoyama & Kogiso, 1965; Høier, 1969; Olsen, 1976a; Jones, Rackham & Steeds, 1977). In these methods, the lattice parameters can be determined, provided the electron wavelength is known (or vice versa) without knowing the value of K.
In the method of Uyeda, Nonoyama & Kogiso (1965), the electron wavelength can be determined from a single Kikuchi pattern provided that the lattice parameters of the crystal are known. Fig. 5.4.2.1 illustrates the geometry involved in the method. A and B are two zone axes of the specimen. The pairs of Kikuchi lines g, −g and h, −h belong to the zones A and B, respectively. The points P and Q are the intersections of the Kikuchi lines g and h with the line AB.
As a first approximation, the wavelength of the electrons is given by where D is the measured distance between the two Kikuchi lines in either of the pairs and d is the corresponding interplanar spacing. is the angle between the zone axes A and B. When the Kikuchi pattern is indexed, and d can be calculated from the known lattice parameters of the specimen. and are effective interplanar spacings and are given by and are the interplanar spacings corresponding to g and h, respectively. The angle (or ) is the angle between the lattice plane g (or h) and the plane defined by A and B. These angles are given by The values of and to be used in (5.4.2.1) can be calculated from (5.4.2.2) and (5.4.2.3).
If Δ and D are measured on the photographic plate or a print of any magnification, the electron wavelength can be calculated by using (5.4.2.1). Only the ratio (Δ/D) is required for the determination of λ. The distance AB is not required. The points A and B are used only to fix the points P and Q. A and B do not need to fall inside the photographic plate when a Kikuchi pattern is symmetrical across the line AB, because in such cases P and Q can be determined from intersections of Kikuchi lines that are symmetrical with each other. The expression for λ in (5.4.2.1) is only a first approximation. More accurate expressions can be found in the paper by Uyeda, Nonoyama & Kogiso (1965).
A simpler and more accurate method has been developed by Høier (1969). He showed that in the cubic case it is possible to determine the ratio between the lattice parameter a and the electron wavelength λ with a relative accuracy of 0.1%. Only two quantities have to be measured on the photographic plate or a print at any magnification: the height of a triangle formed by three Kikuchi lines and one separation between a defect–excess line pair (Fig. 5.4.2.2 ). If three indexed Kikuchi lines intersect at the same point on a photographic plate, the following equations can be derived from Bragg's law: In addition, the length of K is equal to the electron wavelength where the wavelength is written for generality (for electrons, λ_{i} = λ). If the Kikuchi lines do not belong to the same zone, (5.4.2.4) and (5.4.2.5) can be solved and a/λ determined.
Exact triple intersections are rare. A practical method as proposed by Høier (1969) is therefore based on three lines that nearly intersect at the same point (Fig. 5.4.2.2). The dimensions of the triangle ABC in Fig. 5.4.2.2 change with λ. Let us assume that only the wavelength in the beam giving the reflection varies. The increment necessary to shift the line to A is then determined by By measuring the height in the triangle ABC and the line separation , the wavelength to be used in (5.4.2.5) can be calculated and the ratio a/λ determined. Care must be taken to avoid areas where the lines are displaced from kinematical positions owing to dynamical interactions.
The method of Høier (1969) for cubic crystals was later extended to lowersymmetry cases (triclinic) by Olsen (1976a), who also developed computer programs for leastsquares refinement of the lattice parameters. The derivation of the equations for the procedure is carried out in the following in a way slightly different from that described by Olsen (1976a).
If three Kikuchi lines not belonging to the same zone intersect at the same point on the photographic plate, the direction K from the origin of the Ewald sphere to the intersection point of the Kikuchi lines is given by Bragg's law: The length of the vector K is given by where λ is the electron wavelength.
Because triple intersections are rare, a practical method is therefore based on three lines that nearly intersect at the same point, as proposed by Høier (1969). In order to obtain an exact intersection, one of the lines can be changed from to , where is a vector approximately parallel to . For this Kikuchi line (in the following assumed to be line no. 3), the Bragg condition (5.4.2.7) gives For a pair of Kikuchi lines, the simple relation holds, where is the line separation measured on the photographic plate of the centres of the two Kikuchi lines. is the shift (measured on the plate) in the position of one of the lines in order to obtain an exact triple intersection.
Substitution of (5.4.2.10) into (5.4.2.9) gives
From n measurements of intersections, the following equations are obtained: where are the indices (given in reciprocal space) of the Kikuchi lines, is the Kronecker delta, are the indices (given in real space) of the direction to the intersection number j, and is the interplanar spacing corresponding to the reflection and can be expressed in terms of the cell dimensions in real space. The length of the vectors can also be expressed in terms of the lattice parameters in real space a, b, c, α, β, and γ as Equations (5.4.2.12) can be solved with respect to , and can then be expressed in terms of the lattice parameters by substituting the expressions for from (5.4.2.12) into (5.4.2.13).
If the number of intersections is greater than the number of unknown lattice parameters, a leastsquaresrefinement procedure can be used. This involves minimizing the function
Because the expression for Q is nonlinear in the lattice parameters, a refinement procedure can be used to solve the equations derived from (5.4.2.14) using the leastsquares procedure. The accuracy in the lattice parameters or wavelength can be found by statistical methods. A computer program for refinement of the lattice parameters (or electron wavelength) is available (Olsen, 1976b). In this program, some of the lattice parameters can be held constant or kept equal during the refinement.
Methods based on measurements of distances between zone axes (poles) in the Kikuchi patterns (Thomas, 1970) are not very accurate because optical distortions make longdistance measurements inaccurate on the photographic plate or a magnified print. The methods by Høier (1969) and Olsen (1976a) are based on `near intersections'. Better accuracy may be obtained if the high voltage can be varied in order to obtain exact intersections. A simple method for determination of lattice parameters can be applied for crystals with symmetry as low as orthorhombic (Gjønnes & Olsen, 1984). The method is a simplified version of the Uyeda, Nonoyama & Kogiso (1965) method. If two parallel Kikuchi lines belonging to different zones overlap, then:
These equations can be solved to give the ratio between lattice parameter and electron wavelength.
A simple, rapid procedure for acceleratingvoltage (or electronwavelength) determination of a transmission electron microscope has been described by FitzGerald & Johnson (1984). In their method, it is necessary to measure the ratio of two easily found distances between points defined by the intersections of Kikuchi lines near the zone of a silicon crystal. It is not necessary to index the Kikuchi lines, because the method is based on (a) a particular crystal and (b) an easily recognizable crystal orientation, and because the points between which distances need to be measured are specified in their paper and can be easily found. Polynominals for converting the distance measurements both to electron wavelength and to accelerating voltage are given for the range from 100 to 200 kV. A 300 K temperature change must occur (temperature coefficient of K) before the error in the lattice parameter of silicon due to thermal expansion becomes significant.
A method for measuring small local changes in the lattice parameter of bulk specimens based on selectedarea channelling patterns has been described by Walker & Booker (1982). The method utilizes the small changes in the position of highindex channelling lines due to small changes in the lattice parameters, and is based on scanning electron microscopy (SEM). The method is rapid and convenient, is suitable for bulk specimens, and can be applied to areas only a few micrometres across. The method has been used for Si and GaP and is in this case based on the intersection of pairs of , , type lines (the points A, B, and C in Fig. 5.4.2.3 ). Selectedarea channelling patterns from Si and GaP obtained under precisely the same conditions have the same characteristic array of lines, but with slightly different distances AB and BC. The detection limit of this method is at present 3 parts in 10^{4}.

Schematic diagram showing the indexing of the most prominent lines in the selectedarea channelling pattern near the [111] zone of Si. Accelerating voltage: 25 kV. 
When a small electron probe is used to illuminate a crystal, as in convergentbeam electron diffraction (CBED), very fine lines are often found in the diffraction discs (Steeds, 1979). These lines are called HOLZ lines. They are due to upperlayer diffraction effects and can be used in latticeparameter determination (Rackham, Jones & Steeds, 1974; Jones, Rackham & Steeds, 1977). For highest accuracy, relatively thick crystals are required. The limitation to the accuracy is set either by the weakness of the lines or by energy losses in the specimen. Relative changes in the lattice parameters can be measured to an accuracy of 1 part in 10^{4}, whereas the accuracy in the absolute determination of lattice parameters is typically of an order of magnitude worse. The lattice parameters can be measured from crystal regions as small as 20 nm in diameter with an accuracy better than 1 part in 10^{3}. If more accurate results are desired, it is necessary to make measurements on lines that are not affected by interactions between HOLZ reflections.
References
Edington, J. W. (1975). Electron diffraction in the electron microscope. Monographs in practical electron microscopy in materials science, No. 2. Eindhoven: N. V. Philips Gloeilampenfabrieken.FitzGerald, J. D. & Johnson, A. W. S. (1984). A simplified method of electron microscope voltage measurement. Ultramicroscopy, 12, 231–236.
Gard, J. A. (1976). Electron microscopy in mineralogy, p. 52. Berlin: Springer.
Gjønnes, J. & Olsen, A. (1984). Analytical electron microscopy. JEOL News, 22E, 13–18.
Hirsch, P. B., Howie, A., Nicholson, R. B., Pashley, D. W. & Whelan, M. J. (1965). Electron microscopy of thin crystals. London: Butterworth.
Høier, R. (1969). A method to determine the ratio between lattice parameter and electron wavelength from Kikuchi line intersections. Acta Cryst. A25, 516–518.
International Tables for Crystallography (2005). Vol. A. Heidelberg: Springer.
Jones, P. M., Rackham, G. M. & Steeds, J. W. (1977). Higher order Laue zone effects in electron diffraction and their use in lattice parameter determination. Proc. R. Soc. London Ser. A, 354, 197–222.
LePage, Y. (1992). Ab initio primitive cell parameters from single convergent beam patterns: a converse route to the identification of microcrystals with electrons. Microsc. Res. Tech. 21, 158–165.
Olsen, A. (1976a). Lattice parameter determination using Kikuchiline intersections: application to olivine and feldspar. J. Appl. Cryst. 9, 9–13.
Olsen, A. (1976b). Determination of lattice constants using Kikuchi line intersections. Solid State Group Report Series. Institute of Physics, University of Oslo, Norway.
Rackham, G. M., Jones, P. M. & Steeds, J. W. (1974). Upper layer diffraction effects in zone axis patterns. Proceedings of the Eighth International Congress on Electron Microscopy, Canberra, Australia, pp. 336–337.
Steeds, J. W. (1979). Convergent beam electron diffraction. Introduction to analytical electron microscopy, edited by J. J. Hren, J. I. Goldstein & D. C. Joy, pp. 387–422. New York: Plenum.
Thomas, G. (1970). Kikuchi electron diffraction and applications. Modern diffraction and imaging techniques in material science, edited by S. Amelinckx, S. Gevers, G. Remaut & J. Van Landuyt, pp. 131–185. Amsterdam: NorthHolland.
Uyeda, R., Nonoyama, M. & Kogiso, M. (1965). Determination of the wavelength of electrons from a Kikuchi pattern. J. Electron Microsc. 14, 296–300.
Walker, A. R. & Booker, G. R. (1982). A selectedarea channelling pattern (SACP) method for measuring small local changes in lattice parameter with bulk specimens. Electron microscopy 1982, Vol. 1, pp. 651–652. Hamburg: Elsevier.
Zuo, J. M. (1993). New method of Bravais lattice determination. Ultramicroscopy, 52, 459–464.