International
Tables for Crystallography Volume B Reciprocal space Edited by U. Shmueli © International Union of Crystallography 2006 
International Tables for Crystallography (2006). Vol. B, ch. 2.5, pp. 306310

2.5.4. Electrondiffraction structure analysis (EDSA)^{4}
Electrondiffraction structure analysis (EDSA) (Vainshtein, 1964) based on electron diffraction (Pinsker, 1953) is used for the investigation of the atomic structure of matter together with Xray and neutron diffraction analysis. The peculiarities of EDSA, as compared with Xray structure analysis, are defined by a strong interaction of electrons with the substance and by a short wavelength λ. According to the Schrödinger equation (see Section 5.2.2 ) the electrons are scattered by the electrostatic field of an object. The values of the atomic scattering amplitudes, , are three orders higher than those of Xrays, , and neutrons, . Therefore, a very small quantity of a substance is sufficient to obtain a diffraction pattern. EDSA is used for the investigation of very thin singlecrystal films, of 5–50 nm polycrystalline and textured films, and of deposits of finely grained materials and surface layers of bulk specimens. The structures of many ionic crystals, crystal hydrates and hydrooxides, various inorganic, organic, semiconducting and metalloorganic compounds, of various minerals, especially layer silicates, and of biological structures have been investigated by means of EDSA; it has also been used in the study of polymers, amorphous solids and liquids.
Special areas of EDSA application are: determination of unit cells; establishing orientational and other geometrical relationships between related crystalline phases; phase analysis on the basis of and sets; analysis of the distribution of crystallite dimensions in a specimen and inner strains in crystallites as determined from line profiles; investigation of the surface structure of single crystals; structure analysis of crystals, including atomic position determination; precise determination of lattice potential distribution and chemical bonds between atoms; and investigation of crystals of biological origin in combination with electron microscopy (Vainshtein, 1964; Pinsker, 1953; Zvyagin, 1967; Pinsker et al., 1981; Dorset, 1976; Zvyagin et al., 1979).
There are different kinds of electron diffraction (ED) depending on the experimental conditions: highenergy (HEED) (above 30–200 kV), lowenergy (LEED) (10–600 V), transmission (THEED), and reflection (RHEED). In electrondiffraction studies use is made of special apparatus – electrondiffraction cameras in which the lens system located between the electron source and the specimen forms the primary electron beam, and the diffracted beams reach the detector without aberration distortions. In this case, highresolution electron diffraction (HRED) is obtained. ED patterns may also be observed in electron microscopes by a selectedarea method (SAD). Other types of electron diffraction are: MBD (microbeam), HDD (highdispersion), CBD (convergentbeam), SMBD (scanningbeam) and RMBD (rockingbeam) diffraction (see Sections 2.5.2 and 2.5.3). The recent development of electron diffractometry, based on direct intensity registration and measurement by scanning the diffraction pattern against a fixed detector (scintillator followed by photomultiplier), presents a new improved level of EDSA which provides higher precision and reliability of structural data (Avilov et al., 1999; Tsipursky & Drits, 1977; Zhukhlistov et al., 1997, 1998; Zvyagin et al., 1996).
Electrondiffraction studies of the structure of molecules in vapours and gases is a large special field of research (Vilkov et al., 1978). See also Stereochemical Applications of GasPhase Electron Diffraction (1988).
In HEED, the electron wavelength λ is about 0.05 Å or less. The Ewald sphere with radius has a very small curvature and is approximated by a plane. The ED patterns are, therefore, considered as plane cross sections of the reciprocal lattice (RL) passing normal to the incident beam through the point 000, to scale (Fig. 2.5.4.1). The basic formula is where r is the distance from the pattern centre to the reflection, h is the reciprocalspace vector, d is the appropriate interplanar distance and L is the specimentoscreen distance. The deviation of the Ewald sphere from a plane at distance h from the origin of the coordinates is . Owing to the small values of λ and to the rapid decrease of depending on , the diffracted beams are concentrated in a small angular interval ( rad).
Singlecrystal ED patterns image one plane of the RL. They can be obtained from thin ideal crystalline plates, mosaic singlecrystal films, or, in the RHEED case, from the faces of bulk single crystals. Point ED patterns can be obtained more easily owing to the following factors: the small size of the crystals (increase in the dimension of RL nodes) and mosaicity – the small spread of crystallite orientations in a specimen (tangential tension of the RL nodes). The crystal system, the parameters of the unit cell and the Laue symmetry are determined from point ED patterns; the probable space group is found from extinctions. Point ED patterns may be used for intensity measurements if the kinematic approximation holds true or if the contributions of the dynamic and secondary scattering are not too large.
The indexing of reflections and the unitcell determination are carried out according to the formulae relating the RL to the DL (direct lattice) (Vainshtein, 1964; Pinsker, 1953; Zvyagin, 1967).
Under electrondiffraction conditions crystals usually show a tendency to lie down on the substrate plane on the most developed face. Let us take this as (001). The vectors a and b are then parallel, while vector is normal to this plane, and the RL points are considered as being disposed along direct lines parallel to the axis with constant hk and variable l.
The interpretation of the point patterns as respective RL planes is quite simple in the case of orthogonal lattices. If the lattice is triclinic or monoclinic the pattern of the crystal in the position with the face (001) normal to the incident beam does not have to contain hk0 reflections with nonzero h and k because, in general, the planes ab and do not coincide. However, the intersection traces of direct lines hk with the plane normal to them (plane ab) always form a net with periods (Fig. 2.5.4.2). The points hkl along these directions hk are at distances from the ab plane.
By changing the crystal orientation it is possible to obtain an image of the plane containing hk0 reflections, or of other RL planes, with the exception of planes making a small angle with the axis .
In the general case of an arbitrary crystal orientation, the pattern is considered as a plane section of the system of directions hk which makes an angle ϕ with the plane ab, intersecting it along a direction [uv]. It is described by two periods along directions 0h, 0k; with an angle γ″ between them satisfying the relation and by a system of parallel directions The angles are formed by directions 0h, 0k in the plane of the pattern with the plane ab. The coefficients depend on the unitcell parameters, angle ϕ and direction [uv]. These relations are used for the indexing of reflections revealed near the integer positions hkl in the pattern and for unitcell calculations (Vainshtein, 1964; Zvyagin, 1967; Zvyagin et al., 1979).
In RED patterns obtained with an incident beam nearly parallel to the plane ab one can reveal all the RL planes passing through which become normal to the beam at different azimuthal orientations of the crystal.
With the increase of the thickness of crystals (see below, Chapter 5.1 ) the scattering becomes dynamical and Kikuchi lines and bands appear. Kikuchi ED patterns are used for the estimation of the degree of perfection of the structure of the surface layers of single crystals for specimen orientation in HREM (IT C, 2004, Section 4.3.8 ). Patterns obtained with a convergent beam contain Kossel lines and are used for determining the symmetry of objects under investigation (see Section 5.1.2 ).
Texture ED patterns are a widely used kind of ED pattern (Pinsker, 1953; Vainshtein, 1964; Zvyagin, 1967). Textured specimens are prepared by substance precipitation on the substrate, from solutions and suspensions, or from gas phase in vacuum. The microcrystals are found to be oriented with a common (developed) face parallel to the substrate, but they have random azimuthal orientations. Correspondingly, the RL also takes random azimuthal orientations, having as the common axis, i.e. it is a rotational body of the point RL of a single crystal. Thus, the ED patterns from textures bear a resemblance, from the viewpoint of their geometry, to Xray rotation patterns, but they are less complicated, since they represent a plane cross section of reciprocal space.
If the crystallites are oriented by the plane (hkl), then the axis is the texture axis. For the sake of simplicity, let us assume that the basic plane is the plane (001) containing the axes a and b, so that the texture axis is , i.e. the axis . The matrices of appropriate transformations will define a transition to the general case (see IT A , 2005). The RL directions , parallel to the texture axis, transform to cylindrical surfaces, the points with are in planes perpendicular to the texture axis, while any `tilted' lines transform to cones or hyperboloids of rotation. Each point hkl transforms to a ring lying on these surfaces. In practice, owing to a certain spread of axes of single crystals, the rings are blurred into small band sections of a spherical surface with the centre at the point 000; the oblique cross section of such bands produces reflections in the form of arcs. The main interference curves for texture patterns are ellipses imaging oblique plane cross sections of the cylinders hk (Fig. 2.5.4.3).
At the normal electronbeam incidence (tilting angle ) the ED pattern represents a cross section of cylinders perpendicular to the axis , i.e. a system of rings.
On tilting the specimen to an angle ϕ with respect to its normal position (usually ) the patterns image an oblique cross section of the cylindrical RL, and are called obliquetexture (OT) ED patterns. The ellipses and layer lines for orthogonal lattices are the main characteristic lines of ED patterns along which the reflections are arranged. The shortcoming of obliquetexture ED patterns is the absence of reflections lying inside the cone formed by rotation of the straight line coming from the point 000 at an angle around the axis and, in particular, of reflections 00l. However, at ϕ ≲ 60–70° the set of reflections is usually sufficient for structural determination.
For unitcell determination and reflection indexing the values d (i.e. ) are used, and the reflection positions defined by the ellipses hk to which they belong and the values η are considered. The periods are obtained directly from and values. The period , if it is normal to the plane ( being arbitrary), is calculated as For obliqueangled lattices In the general case of obliqueangled lattices the coaxial cylinders hk have radii and it is always possible to use the measured or calculated values in (2.5.4.5a) instead of , since In OT patterns the and η values are represented by the lengths of the small axes of the ellipses and the distances of the reflections hkl from the line of small axes (equatorial line of the pattern)
Analysis of the values gives a, b, γ, while p, s and q are calculated from the values. It is essential that the components of the normal projections of the axis c on the plane ab measured in the units of a and b are Obtaining one can calculate Since
The α, β values are then defined by the relations
Because of the small particle dimensions in textured specimens, the kinematic approximation is more reliable for OT patterns, enabling a more precise calculation of the structure amplitudes from the intensities of reflections.
Polycrystal ED patterns. In this case, the RL is a set of concentric spheres with radii . The ED pattern, like an Xray powder pattern, is a set of rings with radii
The intensities of scattering by a crystal are determined by the scattering amplitudes of atoms in the crystal, given by (see also Section 5.2.1 ) where is the potential of an atom and . The absolute values of have the dimensionality of length L. In EDSA it is convenient to use without K. The dimensionality of is [potential ]. With the expression of in V Å^{3} the value in (2.5.4.13) is 47.87 V Å^{2}.
The scattering atomic amplitudes differ from the respective Xray values in the following: while (electron shell charge), the atomic amplitude at is the `full potential' of the atom. On average, , but for small atomic numbers Z, owing to the peculiarities in the filling of the electron shells, exhibits within periods of the periodic table of elements `reverse motion', i.e. they decrease with Z increasing (Vainshtein, 1952, 1964). At large , . The atomic amplitudes and, consequently, the reflection intensities, are recorded, in practice, up to values of , i.e. up to .
The structure amplitude of a crystal is determined by the Fourier integral of the unitcell potential (see Chapter 1.2 ), where Ω is the unitcell volume. The potential of the unit cell can be expressed by the potentials of the atoms of which it is composed: The thermal motion of atoms in a crystal is taken into account by the convolution of the potential of an atom at rest with the probability function describing the thermal motion: Accordingly, the atomic temperature factor of the atom in a crystal is where the Debye temperature factor is written for the case of isotropic thermal vibrations. Consequently, the structure amplitude is This general expression is transformed (see IT I, 1952) according to the space group of a given crystal.
To determine the structure amplitudes in EDSA experimentally, one has to use specimens satisfying the kinematic scattering condition, i.e. those consisting of extremely thin crystallites. The limit of the applicability of the kinematic approximation (Blackman, 1939; Vainshtein, 1964) can be estimated from the formula where is the averaged absolute value of (see also Section 5.2.1 ). Since are proportional to , condition (2.5.4.20) is better fulfilled for crystals with light and medium atoms. Condition (2.5.4.20) is usually satisfied for textured and polycrystalline specimens. But for mosaic single crystals as well, the kinematic approximation limit is, in view of their real structure, substantially wider than estimated by (2.5.4.20) for ideal crystals. The fulfillment of the kinematic law for scattering can be, to a greater or lesser extent, estimated by comparing the decrease of experimental intensity averaged over definite angular intervals, and sums calculated for the same angular intervals.
For mosaic singlecrystal films the integral intensity of reflection is for textures Here is the incident electronbeam density, S is the irradiated specimen area, t is the thickness of the specimen, α is the average angular spread of mosaic blocks, R′ is the horizontal coordinate of the reflection in the diffraction pattern and p is the multiplicity factor. In the case of polycrystalline specimens the local intensity in the maximum of the ring reflection is measured, where ΔS is the measured area of the ring.
The transition from kinematic to dynamic scattering occurs at critical thicknesses of crystals when (2.5.4.20). Mosaic or polycrystalline specimens then result in an uneven contribution of various crystallites to the intensity of the reflections. It is possible to introduce corrections to the experimental structure amplitudes of the first strong reflections most influenced by dynamic scattering by applying in simple cases the twowave approximation (Blackman, 1939) or by taking into account multibeam theories (Fujimoto, 1959; Cowley, 1981; Avilov et al. 1984; see also Chapter 5.2 ).
The application of kinematic scattering formulae to specimens of thin crystals (5–20 nm) or dynamic corrections to thicker specimens (20–50 nm) permits one to obtain reliability factors between the calculated and observed structure amplitudes of , which is sufficient for structural determinations.
With the use of electron diffractometry techniques, reliability factors as small as 2–3% have been reached and more detailed data on the distribution of the innercrystalline potential field have been obtained, characterizing the state and bonds of atoms, including hydrogen (Zhukhlistov et al., 1997, 1998; Avilov et al., 1999).
The applicability of kinematics formulae becomes poorer in the case of structures with many heavy atoms for which the atomic amplitudes also contain an imaginary component (Shoemaker & Glauber, 1952). The experimental intensity measurement is made by a photo method or by direct recording (Avilov, 1979). In some cases the amplitudes can be determined from dynamic scattering patterns – the bands of equal thickness from a wedgeshaped crystal (Cowley, 1981), or from rocking curves.
The unit cell is defined on the basis of the geometric theory of electrondiffraction patterns, and the space group from extinctions. It is also possible to use the method of converging beams (Section 5.2.2 ). The structural determination is based on experimental sets of values or (Vainshtein, 1964).
The trialanderror method may be used for the simplest structures. The main method of determination is the construction of the Patterson functions and their analysis on the basis of heavyatom methods, superposition methods and so on (see Chapter 2.3 ). Direct methods are also used (Dorset et al., 1979). Thus the phases of structure factors are calculated and assigned to the observed moduli
The distribution of the potential in the unit cell, and, thereby, the arrangement in it of atoms (peaks of the potential) are revealed by the construction of threedimensional Fourier series of the potential (see also Chapter 1.3 ) or projections The general formulae (2.5.4.26a) and (2.5.4.26b) transform, according to known rules, to the expressions for each space group (see IT I, 1952). If are expressed in V and the volume Ω or the cell area S in and , respectively, then the potential ϕ is obtained directly in volts, while the projection of the potential is in V Å. The amplitudes are reduced to an absolute scale either according to a group of strong reflections or using the Parseval equality or Wilson's statistical method The term defines the mean inner potential of a crystal, and is calculated from [(2.5.4.13), (2.5.4.19)] The Fourier series of the potential in EDSA possess some peculiarities (Vainshtein, 1954, 1964) which make them different from the electrondensity Fourier series in Xray analysis. Owing to the peculiarities in the behaviour of the atomic amplitudes (2.5.4.13), which decrease more rapidly with increasing compared with , the peaks of the atomic potential are more `blurred' and exhibit a larger halfwidth than the electrondensity peaks . On average, this halfwidth corresponds to the `resolution' of an electrondiffraction pattern – about 0.5 Å or better. The potential in the maximum (`peak height') does not depend as strongly on the atomic number as in Xray analysis: while in Xray diffraction . In such a way, in EDSA the light atoms are more easily revealed in the presence of heavy atoms than in Xray diffraction, permitting, in particular, hydrogen atoms to be revealed directly without resorting to difference syntheses as in Xray diffraction. Typical values of the atomic potential (which depend on thermal motion) in organic crystals are: H ∼ 35, C ∼ 165, O 215 V; in Al crystals 330 V, in Cu crystals 750 V.
The EDSA method may be used for crystal structure determination, depending on the types of electrondiffraction patterns, for crystals containing up to several tens of atoms in the unit cell. The accuracy in determination of atomic coordinates in EDSA is about 0.01–0.005 Å on average. The precision of EDSA makes it possible to determine accurately the potential distribution, to investigate atomic ionization, to obtain values for the potential between the atoms and, thereby, to obtain data on the nature of the chemical bond.
If the positions in the cell are occupied only partly, then the measurement of gives information on population percentage.
There is a relationship between the nuclear distribution, electron density and the potential as given by the Poisson equation This makes it possible to interrelate Xray diffraction, EDSA and neutrondiffraction data. Thus for the atomic amplitudes where Z is the nuclear charge and the Xray atomic scattering amplitude, and for structure amplitudes where is the Xray structure amplitude of the electron density of a crystal and is the amplitude of scattering from charges of nuclei in the cell taking into account their thermal motion. The values can be calculated easily from neutrondiffraction data, since the charges of the nuclei are known and the experiment gives the parameters of their thermal motion.
In connection with the development of highresolution electronmicroscopy methods (HREM) it has been found possible to combine the data from direct observations with EDSA methods. However, EDSA permits one to determine the atomic positions to a greater accuracy, since practically the whole of reciprocal space with 1.0–0.4 Å resolution is used and the threedimensional arrangement of atoms is calculated. At the same time, in electron microscopy, owing to the peculiarities of electron optics and the necessity for an objective aperture, the image of the atoms in a crystal is a convolution, with the aperture function blurring the image up to 1.5–2 Å resolution. In practice, in TEM one obtains only the images of the heaviest atoms of an object. However, the possibility of obtaining a direct image of a structure with all the defects in the atomic arrangement is the undoubted merit of TEM.
References
International Tables for Crystallography (2005). Vol. A. Spacegroup symmetry, edited by Th. Hahn, 5th ed. Heidelberg: Springer.International Tables for Crystallography (2004). Vol. C. Mathematical, physical and chemical tables, edited E. Prince, 3rd ed. Dordrecht: Kluwer Academic Publishers.
International Tables for Xray Crystallography (1952). Vol. I. Symmetry groups. Birmingham: Kynoch Press. (Present distributor Kluwer Academic Publishers, Dordrecht.)
Stereochemical Applications of GasPhase Electron Diffraction (1988). Part A, edited by I. Hargittai & M. Hargittai. New York: VCH.
Avilov, A. S. (1979). Electrical measurement of reflection intensities on electron diffraction from mosaic single crystals. Sov. Phys. Crystallogr. 24, 103–104.
Avilov, A. S., Kuligin, A. K., Pietsch, U., Spence, J. C. H., Tsirelson, V. G.& Zuo, J. M. (1999). Scanning system for highenergy electron diffractometry. J. Appl. Cryst. 32, 1033–1038.
Avilov, A. S., Parmon, V. S., Semiletov, S. A. & Sirota, M. I. (1984). Calculation of reflected intensities in multiplebeam diffraction of fast electrons by polycrystalline specimens. Sov. Phys. Crystallogr. 29, 5–7.
Blackman, M. (1939). On the intensities of electron diffraction rings. Proc. R. Soc. London Ser. A, 173, 68–82.
Cowley, J. M. (1981). Diffraction physics, 2nd ed. Amsterdam: NorthHolland.
Dorset, D. L. (1976). The interpretation of quasikinematical singlecrystal electron diffraction intensity data from paraffins. Acta Cryst. A32, 207–215.
Dorset, D. L., Jap, B. K., Ho, M.H. & Glaeser, R. M. (1979). Direct phasing of electron diffraction data from organic crystals: the effect of nbeam dynamical scattering. Acta Cryst. A35, 1001–1009.
Fujimoto, F. (1959). Dynamical theory of electron diffraction in Lauecase. I. General theory. J. Phys. Soc. Jpn, 14(11), 1158–1168.
Pinsker, Z. G. (1953). Electron diffraction. London: Butterworth.
Pinsker, Z. G., Zvyagin, B. B. & Imamov, R. M. (1981). Principal results of electrondiffraction structural investigations. Sov. Phys. Crystallogr. 26, 669–674.
Shoemaker, V. & Glauber, R. (1952). The Born approximation in electron diffraction. Nature (London), 170, 290–291.
Tsipursky, S. I. & Drits, V. A. (1977). Efficiency of electronometric intensity registration at electron diffraction structural studies. Izv. Akad. Nauk SSSR Ser. Fiz. 41, 2263–2271. (In Russian.)
Vainshtein, B. K. (1952). Dependence of electron scattering on the atomic number. Dokl. Akad. Nauk SSSR, 85, 1239–1242. (In Russian.)
Vainshtein, B. K. (1954). On the studies of crystal lattice potential by electron diffraction. Tr. Inst. Krist. Akad. Nauk SSSR, 9, 259–276. (In Russian.)
Vainshtein, B. K. (1964). Structure analysis by electron diffraction. Oxford: Pergamon Press.
Vilkov, L. V., Mastryukov, V. S. & Sadova, N. I. (1978). Determination of geometrical structure of free molecules. Leningrad: Khimiya. (In Russian.)
Zhukhlistov, A. P., Avilov, A. S., Ferraris, G., Zvyagin, B. B. & Plotnikov, V. P. (1997). Statistical distribution of hydrogen over three positions in the brucite Mg(OH)_{2} structure from electron diffractometry data. Crystallogr. Rep. 42, 774–777.
Zhukhlistov, A. P. & Zvyagin, B. B. (1998). Crystal structure of lizardite 1T from electron diffractometry data. Crystallogr. Rep. 43, 950–955.
Zvyagin, B. B. (1967). Electrondiffraction analysis of clay mineral structures. New York: Plenum.
Zvyagin, B. B., Vrublevskaya, Z. V., Zhukhlistov, A. P., Sidorenko, S. V., Soboleva, A. F. & Fedotov, A. F. (1979). Highvoltage electron diffraction investigations of layered minerals. Moscow: Nauka. (In Russian).
Zvyagin, B. B., Zhukhlistov, A. P. & Plotnikov, A. P. (1996). Development of the electron diffractometry of minerals. Structural studies of crystals. (Coll. Works 75th Anniversary Acad. B. K. Vainshtein.) NaukaPhysmathlit, pp. 225–234. (In Russian).