Tables for
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2010). Vol. B, ch. 2.5, pp. 356-361   | 1 | 2 |

Section 2.5.4. Electron-diffraction structure analysis (EDSA)2

B. K. Vainshteinc and B. B. Zvyagind

2.5.4. Electron-diffraction structure analysis (EDSA)2

| top | pdf | Introduction

| top | pdf |

Electron-diffraction structure analysis (EDSA) (Vainshtein, 1964[link]) based on electron diffraction (Pinsker, 1953[link]) is used for the investigation of the atomic structure of matter together with X-ray and neutron diffraction analysis. The peculiarities of EDSA, as compared with X-ray 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[link] ) the electrons are scattered by the electrostatic field of an object. The values of the atomic scattering amplitudes, fe, are three orders higher than those of X-rays, fx, and neutrons, fn. Therefore, a very small quantity of a substance is sufficient to obtain a diffraction pattern. EDSA is used for the investigation of very thin single-crystal 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 hydro-oxides, various inorganic, organic, semiconducting and metallo-organic 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 [d_{hkl}] and [I_{hkl}] 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[link]; Pinsker, 1953[link]; Zvyagin, 1967[link]; Pinsker et al., 1981[link]; Dorset, 1976[link]; Zvyagin et al., 1979[link]).

There are different kinds of electron diffraction (ED) depending on the experimental conditions: high-energy (HEED) (above 30–200 kV), low-energy (LEED) (10–600 V), transmission (THEED) and reflection (RHEED). In electron-diffraction studies use is made of special apparatus – electron-diffraction 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, high-resolution electron diffraction (HRED) is obtained. ED patterns may also be observed in electron microscopes by a selected-area method (SAD). Other types of electron diffraction are: MBD (microbeam), HDD (high-dispersion), CBD (convergent-beam), SMBD (scanning-beam) and RMBD (rocking-beam) diffraction (see Sections 2.5.2[link] and 2.5.3[link]). 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[link]; Tsipursky & Drits, 1977[link]; Zhukhlistov et al., 1997[link], 1998[link]; Zvyagin et al., 1996[link]).

Electron-diffraction studies of the structure of molecules in vapours and gases is a large special field of research (Vilkov et al., 1978[link]). See also Stereochemical Applications of Gas-Phase Electron Diffraction (1988)[link]. The geometry of ED patterns

| top | pdf |

In HEED, the electron wavelength λ is about 0.05 Å or less. The Ewald sphere with radius [\lambda^{-1}] 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 [L\lambda] (Fig.[link]). The basic formula is [r = |{\bf h}| L\lambda, \hbox{ or } rd = L\lambda, \eqno(]where r is the distance from the pattern centre to the reflection, h is the reciprocal-space vector, d is the appropriate interplanar distance and L is the specimen-to-screen distance. The deviation of the Ewald sphere from a plane at distance h from the origin of the coordinates is [\delta_{h} = h^{2} \lambda / 2]. Owing to the small values of λ and to the rapid decrease of fe depending on [(\sin \theta)/\lambda], the diffracted beams are concentrated in a small angular interval ([\leq 0.1] rad).


Figure | top | pdf |

Ewald spheres in reciprocal space. Dotted line: electrons, solid line: X-rays.

Single-crystal ED patterns image one plane of the RL. They can be obtained from thin ideal crystalline plates, mosaic single-crystal 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 unit-cell determination are carried out according to the formulae relating the RL to the DL (direct lattice) (Vainshtein, 1964[link]; Pinsker, 1953[link]; Zvyagin, 1967[link]).

Under electron-diffraction 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 [{\bf c}^{*}] is normal to this plane, and the RL points are considered as being disposed along direct lines parallel to the axis [c^{*}] 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 [a^{*}b^{*}] 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 [(a \sin \gamma)^{-1}, (b \sin \gamma)^{-1} \hbox{ and angle } \gamma' = \pi - \gamma \eqno(](Fig.[link]). The points hkl along these directions hk are at distances [\eta = ha^{*} \cos \beta^{*} + kb^{*} \cos \alpha^{*} + lc^{*} \eqno(]from the ab plane.


Figure | top | pdf |

Triclinic reciprocal lattice. Points: open circles, projection net: black circles.

By changing the crystal orientation it is possible to obtain an image of the [a^{*} b^{*}] plane containing hk0 reflections, or of other RL planes, with the exception of planes making a small angle with the axis [c^{*}].

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; [(a \sin \gamma \cos \psi_{h})^{-1}, (b \sin \gamma \cos \psi_{k})^{-1}, \eqno(]with an angle γ″ between them satisfying the relation [\cos \gamma'' = \sin \psi_{h} \sin \psi_{k} - \cos \psi_{h} \cos \psi_{k} \cos \gamma, \eqno(]and by a system of parallel directions [p_{h}h + p_{k}k = l\semi \quad l = 0, \pm 1, \pm 2,\ldots. \eqno(]The angles [\psi_{h}, \psi_{k}] are formed by directions 0h, 0k in the plane of the pattern with the plane ab. The coefficients [p_{h}, p_{k}] depend on the unit-cell 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 unit-cell calculations (Vainshtein, 1964[link]; Zvyagin, 1967[link]; Zvyagin et al., 1979[link]).

In RED patterns obtained with an incident beam nearly parallel to the plane ab one can reveal all the RL planes passing through [c^{*}] 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[link] ) 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[link], Section 4.3.8[link] ). Patterns obtained with a convergent beam contain Kossel lines and are used for determining the symmetry of objects under investigation (see Section 2.5.3[link]).

Texture ED patterns are a widely used kind of ED pattern (Pinsker, 1953[link]; Vainshtein, 1964[link]; Zvyagin, 1967[link]). 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 [c^{*}] 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 X-ray 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 [[hkl]^{*}] 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 [[001]^{*}], i.e. the axis [c^{*}]. The matrices of appropriate transformations will define a transition to the general case (see IT A , 2005[link]). The RL directions [hk = \hbox{constant}], parallel to the texture axis, transform to cylindrical surfaces, the points with [\eta_{hkl} = \hbox{constant}] 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 [c^{*}] 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.[link]).


Figure | top | pdf |

Formation of ellipses on an electron-diffraction pattern from an oblique texture.

At the normal electron-beam incidence (tilting angle [\varphi = 0^{\circ}]) the ED pattern represents a cross section of cylinders perpendicular to the axis [c^{*}], i.e. a system of rings.

On tilting the specimen to an angle ϕ with respect to its normal position (usually [\varphi \simeq 60^{\circ}]) the patterns image an oblique cross section of the cylindrical RL, and are called oblique-texture (OT) ED patterns. The ellipses [(hk = \hbox{constant})] and layer lines [(l = \hbox{constant})] for orthogonal lattices are the main characteristic lines of ED patterns along which the reflections are arranged. The shortcoming of oblique-texture 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 [(90^{\circ} - \varphi)] around the axis [c^{*}] and, in particular, of reflections 00l. However, at ϕ ≲ 60–70° the set of reflections is usually sufficient for structural determination.

For unit-cell determination and reflection indexing the values d (i.e. [|{\bf h}|]) are used, and the reflection positions defined by the ellipses hk to which they belong and the values η are considered. The periods [a^{*}, b^{*}] are obtained directly from [h_{100}] and [h_{010}] values. The period [c^{*}], if it is normal to the plane [a^{*} b^{*}] ([\gamma^{*}] being arbitrary), is calculated as [c^{*} = \eta/l = (h_{hkl}^{2} - h_{hk0}^{2})^{1/2}/l. \eqno(]For oblique-angled lattices [c^{*} = [(h_{{l_1}+l}^{2} + h_{{l_1}-l}^{2} - 2h_{l}^{2})/2]^{1/2}/l. \eqno(]In the general case of oblique-angled lattices the coaxial cylinders hk have radii [\eqalignno{ b_{hk} &= (1/\sin \gamma) [(h^{2} / a^{2}) + (k^{2} / b^{2}) &\cr &\quad - (2hk \cos \gamma / ab)]^{1/2} &(\cr}]and it is always possible to use the measured or calculated values [b_{hk}] in ([link] instead of [h_{hk0}], since [\eta = (h_{hkl}^{2} - b_{hk}^{2})^{1/2}. \eqno(]In OT patterns the [b_{hk}] and η values are represented by the lengths of the small axes of the ellipses [B_{hk} = L\lambda b_{hk}] and the distances of the reflections hkl from the line of small axes (equatorial line of the pattern) [D_{hkl} = L\lambda \eta /\sin \varphi = hp + ks + lq. \eqno(]

Analysis of the [B_{hk}] values gives a, b, γ, while p, s and q are calculated from the [D_{hkl}] values. It is essential that the components of the normal projections [c_{\rm n}] of the axis c on the plane ab measured in the units of a and b are [\eqalign{ x_{\rm n} &= (c / a) (\cos \beta - \cos \alpha \cos \gamma)/\sin^{2} \gamma\cr &= - p / q,\cr y_{\rm n} &= (c / b) (\cos \alpha - \cos \beta \cos \gamma)/\sin^{2} \gamma\cr &= - s / q.\cr} \eqno(]Obtaining [x_{\rm n}, y_{\rm n}] one can calculate [c_{\rm n} = [(x_{\rm n} a)^{2} + (y_{\rm n} b)^{2} + 2x_{\rm n} y_{\rm n} ab \cos \gamma]^{1/2}. ]Since [\eqalign{ d_{001} &= L\lambda / q \sin \varphi,\cr c &= (c_{\rm n}^{2} + d_{001}^{2})^{1/2}.\cr} \eqno(]

The α, β values are then defined by the relations [\eqalign{ \cos \alpha &= (x_{\rm n} a \cos \gamma + y_{\rm n} b)/c,\cr \cos \beta &= (x_{\rm n} a + y_{\rm n} b \cos \gamma)/c.\cr} \eqno(]

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 [h_{hkl}]. The ED pattern, like an X-ray powder pattern, is a set of rings with radii [r_{hkl} = h_{hkl} L\lambda. \eqno(] Intensities of diffraction beams

| top | pdf |

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[link] ) [\eqalign{ f_{\rm e}^{\rm abs} (s) &= 4\pi K \int\varphi (r) r^{2} {\sin sr \over sr} \ \hbox{d}r\semi\cr K &= {2\pi me \over h^{2}}\semi\,f_{\rm e} = K^{-1} f_{\rm e}^{\rm abs},\cr} \eqno(]where [\varphi (r)] is the potential of an atom and [s = 4\pi (\sin \theta)/\lambda]. The absolute values of [f_{\rm e}^{\rm abs}] have the dimensionality of length L. In EDSA it is convenient to use [f_{\rm e}] without K. The dimensionality of [f_{\rm e}] is [potential [L^{3}]]. With the expression of [f_{\rm e}] in V Å3 the value [K^{-1}] in ([link] is 47.87 V Å2.

The scattering atomic amplitudes [f_{\rm e} (s)] differ from the respective [f_{\rm x} (s)] X-ray values in the following: while [f_{\rm x} (0) = Z] (electron shell charge), the atomic amplitude at [s = 0] [f_{\rm e} (0) = 4\pi {\textstyle\int} \varphi (r) r^{2} \ \hbox{d}r \eqno(]is the `full potential' of the atom. On average, [f_{\rm e} (0) \simeq Z^{1/3}], but for small atomic numbers Z, owing to the peculiarities in the filling of the electron shells, [f_{\rm e} (0)] exhibits within periods of the periodic table of elements `reverse motion', i.e. they decrease with Z increasing (Vainshtein, 1952[link], 1964[link]). At large [(\sin \theta)/\lambda], [f_{\rm e} \simeq Z]. The atomic amplitudes and, consequently, the reflection intensities, are recorded, in practice, up to values of (sin θ)/λ ≃ 0.8–1.2 Å−1, i.e. up to dmin ≃ 0.4–0.6 Å.

The structure amplitude [\Phi_{hkl}] of a crystal is determined by the Fourier integral of the unit-cell potential (see Chapter 1.2[link] ), [\Phi_{hkl} = {\textstyle\int\limits_{\Omega}} \varphi ({\bf r}) \exp \{2\pi i ({\bf r} \cdot {\bf h})\} \hbox{ d}v_{r}, \eqno(]where Ω is the unit-cell volume. The potential of the unit cell can be expressed by the potentials of the atoms of which it is composed: [\varphi ({\bf r}) = {\textstyle\sum\limits_{{\rm cell}, \, i}} \varphi_{{\rm at}\, i} ({\bf r} - {\bf r}_{i}). \eqno(]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 [w({\bf r})] describing the thermal motion: [\varphi_{\rm at} = \varphi_{\rm at} ({\bf r}) * w ({\bf r}). \eqno(]Accordingly, the atomic temperature factor of the atom in a crystal is [{f_{{\rm e}T} [(\sin \theta) / \lambda] = f_{\rm e}\,f_{T} = f_{\rm e} [(\sin \theta) / \lambda] \exp \{- B [(\sin \theta) / \lambda]^{2}\},} \eqno(]where the Debye temperature factor is written for the case of isotropic thermal vibrations. Consequently, the structure amplitude is [\Phi_{hkl} = {\textstyle\sum\limits_{{\rm cell}, \, i}} f_{{\rm e}T_{i}} \exp \{2\pi i (hx_{i} + ky_{i} + lz_{i})\}. \eqno(]This general expression is transformed (see IT I, 1952[link]) 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[link]; Vainshtein, 1964[link]) can be estimated from the formula [A = \lambda \left|{\langle \Phi_{{\bf h}} \rangle \over \Omega}\right| t {{\lower7pt\hbox{$\lt$}}\atop{\raise2pt\hbox{$\sim$}}} 1, \eqno(]where [\langle \Phi_{{\bf h}} \rangle] is the averaged absolute value of [\Phi_{{\bf h}}] (see also Section 2.5.2[link]). Since [\langle \Phi_{{\bf h}} \rangle] are proportional to [Z^{0.8}], condition ([link] is better fulfilled for crystals with light and medium atoms. Condition ([link] 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 ([link] 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 [I_{h} [(\sin \theta)/\lambda]] averaged over definite angular intervals, and sums [{\textstyle\sum}\, f_{\rm obs}^{2} [(\sin \theta)/\lambda]] calculated for the same angular intervals.

For mosaic single-crystal films the integral intensity of reflection is [I_{{\bf h}} = j_{0} S\lambda^{2} \left|{\Phi_{{\bf h}} \over \Omega}\right|^{2} {td_{{\bf h}} \over \alpha} \simeq \Phi_{{\bf h}}^{2} d_{{\bf h}}\semi \eqno(]for textures [I_{{\bf h}} = j_{0} S\lambda^{2} \left|{\Phi_{{\bf h}} \over \Omega}\right|^{2} {tL\lambda p \over 2\pi R' \sin \varphi} \simeq \Phi_{{\bf h}}^{2} \,p / R'. \eqno(]Here [j_{0}] is the incident electron-beam 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 [I_{{\bf h}} = j_{0} S\lambda^{2} \left|{\Phi_{{\bf h}} \over \Omega}\right|^{2} {td_{{\bf h}}^{2} \,p \Delta S \over 4\pi L\lambda} \simeq \Phi_{{\bf h}}^{2} d_{{\bf h}}^{2}\, p \eqno(]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 [A \geq 1] ([link]. 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 two-wave approximation (Blackman, 1939[link]) or by taking into account multibeam theories (Fujimoto, 1959[link]; Cowley, 1981[link]; Avilov et al. 1984[link]; see also Chapter 5.2[link] ).

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 [\Phi^{\rm calc}] and observed [\Phi^{\rm obs}] structure amplitudes of R = 5–15%, which is sufficient for structural determinations.

With the use of electron-diffractometry techniques, reliability factors as small as [R =] 2–3% have been reached and more detailed data on the distribution of the inner-crystalline potential field have been obtained, characterizing the state and bonds of atoms, including hydrogen (Zhukhlistov et al., 1997[link], 1998[link]; Avilov et al., 1999[link]).

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[link]). The experimental intensity measurement is made by a photo method or by direct recording (Avilov, 1979[link]). In some cases the amplitudes [\Phi_{hkl}] can be determined from dynamic scattering patterns – the bands of equal thickness from a wedge-shaped crystal (Cowley, 1981[link]), or from rocking curves. Structure analysis

| top | pdf |

The unit cell is defined on the basis of the geometric theory of electron-diffraction patterns, and the space group from extinctions. It is also possible to use the method of converging beams (Section 2.5.3[link]). The structural determination is based on experimental sets of values [|\Phi_{hkl}|^{2}] or [|\Phi_{hkl}|] (Vainshtein, 1964[link]).

The trial-and-error method may be used for the simplest structures. The main method of determination is the construction of the Patterson functions [P(xyz) = {1 \over \Omega} \left[\Phi_{000}^{2} + 2 \sum\limits_{hkl = -\infty}^{hkl = +\infty} |\Phi_{hkl}|^{2} \cos 2\pi (hx + ky + lz)\right] \eqno(]and their analysis on the basis of heavy-atom methods, superposition methods and so on (see Chapter 2.3[link] ). Direct methods are also used (Dorset et al., 1979[link]). Thus the phases of structure factors are calculated and assigned to the observed moduli [\Phi_{\bf h} = |\Phi_{{\bf h}, \, \rm obs}| \exp \{i\alpha_{\rm calc}\}. \eqno(]

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 three-dimensional Fourier series of the potential (see also Chapter 1.3[link] ) [\varphi (xyz) = {1 \over \Omega} \sum\limits_{{\bf h}} \Phi_{hkl} \exp \{-2\pi i (hx + ky + lz)\} \eqno(]or projections [\varphi' (xy) = {1 \over S} \sum\limits_{{\bf h}} \Phi_{hk0} \exp \{-2\pi i (hx + ky)\}. \eqno(]The general formulae ([link] and ([link] transform, according to known rules, to the expressions for each space group (see IT I, 1952[link]). If [\Phi_{hkl}] are expressed in V Å3 and the volume Ω or the cell area S in Å3 and Å2, respectively, then the potential ϕ is obtained directly in volts, while the projection of the potential [\varphi'] is in V Å. The amplitudes [|\Phi_{hkl}|] are reduced to an absolute scale either according to a group of strong reflections [{\textstyle\sum} |\Phi_{{\bf h}}|^{\rm calc} = {\textstyle\sum} |\Phi_{{\bf h}}|^{\rm obs} \eqno(]or using the Parseval equality [\sum\limits_{{\bf h} = -\infty}^{+\infty} |\Phi_{{\bf h}}|^{2} = \Omega^{2} \langle \varphi^{2}\rangle = \Omega \sum\limits_{i({\rm cell})} {1 \over 2\pi^{2}} \int\limits_{0}^{\infty} f_{{\rm e}T_{i}}^{2} (s)\, s^{2} \hbox{ d}s \eqno (]or Wilson's statistical method [\langle \Phi^{2} [(\sin \theta) / \lambda]\rangle = {\textstyle\sum\limits_{i}} \,f_{{\rm e}T_{i}}^{2} [(\sin \theta) / \lambda]. \eqno(]The term [\Phi_{000}] defines the mean inner potential of a crystal, and is calculated from [f_{\rm e}(0)] [([link], ([link]] [\langle \varphi_{\rm cr}\rangle = \Phi_{000}/\Omega = {1 \over \Omega} \sum f_{\rm e}(0). \eqno(]The Fourier series of the potential in EDSA possess some peculiarities (Vainshtein, 1954[link], 1964[link]) which make them different from the electron-density Fourier series in X-ray analysis. Owing to the peculiarities in the behaviour of the atomic amplitudes ([link], which decrease more rapidly with increasing [(\sin \theta) / \lambda] compared with [f_{\rm x}], the peaks of the atomic potential [\varphi_{\rm at} (r) = {1 \over 2\pi^{2}} \int f_{{\rm e}T} (s) {\sin sr \over sr} s^{2} \hbox{ d}s \eqno(]are more `blurred' and exhibit a larger half-width than the electron-density peaks [\rho_{\rm at}(r)]. On average, this half-width corresponds to the `resolution' of an electron-diffraction pattern – about 0.5 Å or better. The potential in the maximum (`peak height') does not depend as strongly on the atomic number as in X-ray analysis: [\varphi (0) = {1 \over 2\pi^{2}} \int f_{{\rm e}T}(s)\, s^{2} \hbox{ d}s \sim Z^{0.75}, \eqno(]while in X-ray diffraction [\rho (0) \sim Z^{1.2}]. In such a way, in EDSA the light atoms are more easily revealed in the presence of heavy atoms than in X-ray diffraction, permitting, in particular, hydrogen atoms to be revealed directly without resorting to difference syntheses as in X-ray diffraction. Typical values of the atomic potential [\varphi (0)] (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 electron-diffraction 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 [\varphi_{i}(0)] gives information on population percentage.

There is a relationship between the nuclear distribution, electron density and the potential as given by the Poisson equation [\nabla^{2} \varphi ({\bf r}) = -4 \pi e [\rho_{+} ({\bf r}) - \rho_{-} ({\bf r})]. \eqno(]This makes it possible to interrelate X-ray diffraction, EDSA and neutron-diffraction data. Thus for the atomic amplitudes [f_{\rm e}(s) = 4\pi Ke [Z - f_{\rm x} (s)]\,s^{-2}, \eqno(]where Z is the nuclear charge and [f_{\rm x}] the X-ray atomic scattering amplitude, and for structure amplitudes [\Phi_{hkl} = \pi Ke [Z_{hkl} - F_{hkl}]|{\bf h}|^{-2}, \eqno(]where [F_{hkl}] is the X-ray structure amplitude of the electron density of a crystal and [Z_{hkl}] is the amplitude of scattering from charges of nuclei in the cell taking into account their thermal motion. The values [Z_{hkl}] can be calculated easily from neutron-diffraction 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 high-resolution electron-microscopy 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 three-dimensional 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 [\varphi' ({\bf x}) * A({\bf x})] 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.


International Tables for Crystallography (2004). Vol. C. Mathematical, Physical and Chemical Tables, edited E. Prince, 3rd ed. Dordrecht: Kluwer Academic Publishers.
International Tables for Crystallography (2005). Vol. A. Space-Group Symmetry, edited by Th. Hahn, 5th ed. Heidelberg: Springer.
International Tables for X-ray Crystallography (1952). Vol. I. Symmetry Groups. Birmingham: Kynoch Press.
Stereochemical Applications of Gas-Phase 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 high-energy 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 multiple-beam 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: North-Holland.
Dorset, D. L. (1976). The interpretation of quasi-kinematical single-crystal 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 n-beam dynamical scattering. Acta Cryst. A35, 1001–1009.
Fujimoto, F. (1959). Dynamical theory of electron diffraction in Laue-case. I. General theory. J. Phys. Soc. Jpn, 14, 1158–1168.
Pinsker, Z. G. (1953). Electron Diffraction. London: Butterworth.
Pinsker, Z. G., Zvyagin, B. B. & Imamov, R. M. (1981). Principal results of electron-diffraction 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.
Vainshtein, B. K. (1952). Dependence of electron scattering on the atomic number. Dokl. Akad. Nauk SSSR, 85, 1239–1242.
Vainshtein, B. K. (1954). On the studies of crystal lattice potential by electron diffraction. Tr. Inst. Krist. Akad. Nauk SSSR, 9, 259–276.
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.
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). Electron-Diffraction 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). High-Voltage Electron Diffraction Investigations of Layered Minerals. Moscow: Nauka.
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.) Nauka-Physmathlit, pp. 225–234.

to end of page
to top of page