The contributions of electron scattering to the study of the structures of crystalline solids are many and diverse. This section will deal only with the scattering of highenergy electrons (in the energy range of 10^{4} to 10^{6} eV) in transmission through thin samples of crystalline solids and the derivation of information on crystal structures from diffraction patterns and highresolution images. The range of wavelengths considered is from about 0.122 Å (12.2 pm) for 10 kV electrons to 0.0087 Å (0.87 pm) for 1 MeV electrons. Given that the scattering amplitudes of atoms for electrons have much the same form and variation with as for Xrays, it is apparent that the angular range for strong scattering of electrons will be of the order of 10^{−2} rad. Only under special circumstances, usually involving multiple elastic and inelastic scattering from very thick specimens, are scattering angles of more than 10^{−1} rad of importance.
The strength of the interaction of electrons with matter is greater than that of Xrays by two or three orders of magnitude. The singlescattering, first Born approximation fails significantly for scattering from single heavy atoms. Diffracted beams from single crystals may attain intensities comparable with that of the incident beam for crystal thicknesses of 10^{2} Å, rather than 10^{4} Å or more. It follows that electrons may be used for the study of very thin samples, and that dynamical scattering effects, or the coherent interaction of multiply scattered electron waves, will modify the diffracted amplitudes in a significant way for all but very thin specimens containing only light atoms.
The experimental techniques for electron scattering are largely determined by the possibility of focusing electron beams by use of strong axial magnetic fields, which act as electron lenses having focal lengths as short as 1 mm or less. Electron microscopes employing such lenses have been produced with resolutions approaching 1 Å. With such instruments, images showing individual isolated atoms of moderately high atomic number may be obtained. The resolution available is sufficient to distinguish neighbouring rows of adjacent atoms in the projected structures of thin crystals viewed in favourable orientations. It is therefore possible in many cases to obtain information on the structure of crystals and of crystal defects by direct inspection of electron micrographs.
The electromagnetic electron lenses may also be used to form electron beams of very small diameter and very high intensity. In particular, by the use of cold fieldemission electron guns, it is possible to obtain a current of 10^{−10} A in an electron beam of diameter 10 Å or less with a beam divergence of less than 10^{−2} rad, i.e. a current density of 10^{4} A cm^{−2} or more. The magnitudes of the electron scattering amplitudes then imply that detectable signals may be obtained in diffraction from assemblies of fewer than 10^{2} atoms. On the other hand, electron beams may readily be collimated to better than 10^{−6} rad.
The cross sections for inelastic scattering processes are, in general, less than for the elastic scattering of electrons, but signals may be obtained by the observation of electron energy losses, or the production of secondary radiations, which allow the analysis of chemical compositions or electronic excited states for regions of the crystal 100 Å or less in diameter.
On the other hand, the transfer to the sample of large amounts of energy through inelastic scattering processes produces radiation damage which may severely limit the applicability of the imaging and diffraction techniques, especially for biological and organic materials, unless the information is gathered from large specimen volumes with low incident electron beam densities.
Structure analysis of crystals can be performed using electron diffraction in the same way as with Xray or neutron diffraction. The mathematical expressions and the procedures are much the same. However, there are peculiarities of the electrondiffraction case which should be noted.
(1) Structure analysis based on electron diffraction is possible for thin specimens for which the conditions for kinematical scattering are approached, e.g. for thin mosaic singlecrystal specimens, for thin polycrystalline films having a preferred orientation of very small crystallites or for very extensive, very thin single crystals of biological molecules such as membranes one or a few molecules thick.
(2) Dynamical diffraction effects are used explicitly in the determination of crystal symmetry (with no Friedel's law limitations) and for the measurement of structure amplitudes with high accuracy.
(3) For many radiationresistant materials, the structures of crystals and of some molecules may be determined directly by imaging atom positions in projections of the crystal with a resolution of 2 Å or better. The information on atom positions is not dependent on the periodicity of the crystal and so it is equally possible to determine the structures of individual crystal defects in favourable cases.
(4) Techniques of microanalysis may be applied to the determination of the chemical composition of regions of diameter 100 Å or less using the same instrument as for diffraction, so that the chemical information may be correlated directly with morphological and structural information.
(5) Crystalstructure information may be derived from regions containing as few as 10^{2} or 10^{3} atoms, including very small crystals and single or multiple layers of atoms on surfaces.

The material of this section is also reviewed in the text by Spence (2003).
(1) The elastic scattering of electrons results from the interaction of the charged electrons with the electrostatic potential distribution, , of the atoms or crystals. An incident electron of kinetic energy eW gains energy in the potential field. Alternatively it may be stated that an incident electron wave of wavelength is diffracted by a region of variable refractive index
(2) The most important inelastic scattering processes are:
(a) thermal diffuse scattering, with energy losses of the order of 2 × 10^{−2} eV, separable from the elastic scattering only with specially devised equipment; the angular distribution of thermal diffuse scattering shows variations with which are much the same as for the Xray case in the kinematical limit;
(b) bulk plasmon excitation, or the excitation of collective energy states of the conduction electrons, giving energy losses of 3 to 30 eV and an angular range of scattering of 10^{−4} to 10^{−3} rad;
(c) surface plasmons, or the excitation of collective energy states of the conduction electrons at discontinuities of the structure, with energy losses less than those for bulk plasmons and a similar angular range of scattering;
(d) interband or intraband excitation of valenceshell electrons giving energy losses in the range of 1 to 10^{2} eV and an angular range of scattering of 10^{−4} to 10^{−2} rad;
(e) innershell excitations, with energy losses of 10^{2} eV or more and an angular range of scattering of 10^{−3} to 10^{−2} rad, depending on the energy losses involved.
(3) In the original treatment by Bethe (1928) of the elastic scattering of electrons by crystals, the Schrödinger equation is written for electrons in the periodic potential of the crystal; i.e. where K_{0} is the wavevector in zero potential (outside the crystal) (magnitude ) and W is the accelerating voltage. The solutions of the equation are Bloch waves of the form where is the incident wavevector in the crystal and h is a reciprocallattice vector. Substitution of (2.5.2.2) and (2.5.2.3) in (2.5.2.1) gives the dispersion equations Here κ is the magnitude of the wavevector in a medium of constant potential (the `inner potential' of the crystal). The refractive index of the electron in the average crystal potential is then Since is positive and of the order of 10 V and W is to V, is positive and of the order of .
Solution of equation (2.5.2.4) gives the Fourier coefficients of the Bloch waves and application of the boundary conditions gives the amplitudes of individual Bloch waves (see Chapter 5.2
).
(4) The experimentally important case of transmission of highenergy electrons through thin specimens is treated on the assumption of a plane wave incident in a direction almost perpendicular to an infinitely extended planeparallel lamellar crystal, making use of the smallangle scattering approximation in which the forwardscattered wave is represented in the paraboloidal approximation to the sphere. The incidentbeam direction, assumed to be almost parallel to the z axis, is unique and the z component of k is factored out to give where and . [See Lynch & Moodie (1972), Portier & Gratias (1981), Tournarie (1962), and Chapter 5.2
.]
This equation is analogous to the timedependent Schrödinger equation with z replacing t. Retention of the ± signs on the righthand side is consistent with both ψ and being solutions, corresponding to propagation in opposite directions with respect to the z axis. The doublevalued solution is of importance in consideration of reciprocity relationships which provide the basis for the description of some dynamical diffraction symmetries. (See Section 2.5.3.)
(5) The integral form of the wave equation, commonly used for scattering problems, is written, for electron scattering, as
The wavefunction within the integral is approximated by using successive terms of a Born series
The first Born approximation is obtained by putting in the integral and subsequent terms are generated by putting in the integral.
For an incident plane wave, and for a point of observation at a large distance from the scattering object , the first Born approximation is generated as where or, putting and collecting the preintegral terms into a parameter μ, This is the Fouriertransform expression which is the basis for the kinematical scattering approximation. It is derived on the basis that all terms for are very much smaller than and so is a weak scattering approximation.
In this approximation, the scattered amplitude for an atom is related to the atomic structure amplitude, , by the relationship, derived from (2.5.2.8), For centrosymmetrical atom potential distributions, the are real, positive and monotonically decreasing with . A measure of the extent of the validity of the first Born approximation is given by the fact that the effect of adding the higherorder terms of the Born series may be represented by replacing in (2.5.2.10) by the complex quantities and for single heavy atoms the phase factor η may vary from 0.2 for to 4 or 5 for large , as seen from the tables of IT C (2004, Section 4.3.3
).
(6) Relativistic effects produce appreciable variations of the parameters used above for the range of electron energies considered. The relativistic values are where is the Compton wavelength, , and Values for these quantities are listed in IT C (2004, Section 4.3.2
). The variations of λ and σ with accelerating voltage are illustrated in Fig. 2.5.2.1. For high voltages, σ tends to a constant value, .


There are two alternative sets of signs for the functions describing wave optics. Both sets have been widely used in the literature. There is, however, a requirement for internal consistency within a particular analysis, independently of which set is adopted. Unfortunately, this requirement has not always been met and, in fact, it is only too easy at the outset of an analysis to make errors in this way. This problem might have come into prominence somewhat earlier were it not for the fact that, for centrosymmetric crystals (or indeed for centrosymmetric projections in the case of planar diffraction), only the signs used in the transmission and propagation functions can affect the results. It is not until the origin is set away from a centre of symmetry that there is a need to be consistent in every sign used.
Signs for electron diffraction have been chosen from two points of view: (1) defining as positive the sign of the exponent in the structurefactor expression and (2) defining the forward propagating freespace wavefunction with a positive exponent.
The second of these alternatives is the one which has been adopted in most solidstate and quantummechanical texts.
The first, or standard crystallographic convention, is the one which could most easily be adopted by crystallographers accustomed to retaining a positive exponent in the structurefactor equation. This also represents a consistent International Tables usage. It is, however, realized that both conventions will continue to be used in crystallographic computations, and that there are by now a large number of operational programs in use.
It is therefore recommended (a) that a particular sign usage be indicated as either standard crystallographic or alternative crystallographic to accord with Table 2.5.2.1, whenever there is a need for this to be explicit in publication, and (b) that either one or other of these systems be adhered to throughout an analysis in a selfconsistent way, even in those cases where, as indicated above, some of the signs appear to have no effect on one particular conclusion.
The forwardscattering approximation to the manybeam dynamical diffraction theory outlined in Chapter 5.2
provides the basis for the calculation of diffraction intensities and electronmicroscope image contrast for thin crystals. [See Cowley (1995), Chapter 5.2
and IT C (2004) Sections 4.3.6
and 4.3.8
.] On the other hand, there are various approximations which provide relatively simple analytical expressions, are useful for the determination of diffraction geometry, and allow estimates to be made of the relative intensities in diffraction patterns and electron micrographs in favourable cases.
(a) The kinematical approximation, derived in Section 2.5.2.2 from the first Born approximation, is analagous to the corresponding approximation of Xray diffraction. It assumes that the scattering amplitudes are directly proportional to the threedimensional Fourier transform of the potential distribution, . so that the potential distribution takes the place of the chargedensity distribution, , relevant for Xray scattering.
The validity of the kinematical approximation as a basis for structure analysis is severely limited. For lightatom materials, such as organic compounds, it has been shown by Jap & Glaeser (1980) that the thickness for which the approximation gives reasonable accuracy for zoneaxis patterns from single crystals is of the order of 100 Å for 100 keV electrons and increases, approximately as , for higher energies. The thickness limits quoted for polycrystalline samples, having crystallite dimensions smaller than the sample thickness, are usually greater (Vainshtein, 1956). For heavyatom materials the approximation is more limited since it may fail significantly for single heavy atoms.
(b) The phaseobject approximation (POA), or highvoltage limit, is derived from the general manybeam dynamical diffraction expression, equation (5.2.13.1), Chapter 5.2
, by assuming the Ewald sphere curvature to approach zero. Then the scattering by a thin sample can be expressed by multiplying the incoming wave amplitude by the transmission function where is the projection of the potential distribution of the sample in the z direction, the direction of the incident beam. The diffractionpattern amplitudes are then given by twodimensional Fourier transform of (2.5.2.16).
This approximation is of particular value in relation to the electron microscopy of thin crystals. The thickness for its validity for 100 keV electrons is within the range 10 to 50 Å, depending on the accuracy and spatial resolution involved, and increases with accelerating voltage approximately as . In computational work, it provides the starting point for the multislice method of dynamical diffraction calculations (IT C, 2004, Section 4.3.6.1
).
(c) The twobeam approximation for dynamical diffraction of electrons assumes that only two beams, the incident beam and one diffracted beam (or two Bloch waves, each with two component amplitudes), exist in the crystal. This approximation has been adapted, notably by Hirsch et al. (1965), for use in the electron microscopy of inorganic materials.
It forms a convenient basis for the study of defects in crystals having small unit cells (metals, semiconductors etc.) and provides good preliminary estimates for the determination of crystal thicknesses and structure amplitudes for orientations well removed from principal axes, and for electron energies up to 200–500 keV, but it has decreasing validity, even for favourable cases, for higher energies. It has been used in the past as an `extinction correction' for powderpattern intensities (Vainshtein, 1956).
(d) The Bethe second approximation, proposed by Bethe (1928) as a means for correcting the twobeam approximation for the effects of weakly excited beams, replaces the Fourier coefficients of potential by the `Bethe potentials' Use of these potentials has been shown to account well for the deviations of powderpattern intensities from the predictions of twobeam theory (Horstmann & Meyer, 1965) and to predict accurately the extinctions of Kikuchi lines at particular accelerating voltages due to relativistic effects (Watanabe et al., 1968), but they give incorrect results for the smallthickness limit.

(1) Comparison with Xray diffraction. The relations of realspace and reciprocalspace functions are analogous to those for Xray diffraction [see equations (2.5.2.2), (2.5.2.10) and (2.5.2.15)]. For diffraction by crystals where the integral of (2.5.2.18) and the summation of (2.5.2.19) are taken over one unit cell of volume (see Dawson et al., 1974).
Important differences from the Xray case arise because
(a) the wavelength is relatively small so that the Ewaldsphere curvature is small in the reciprocalspace region of appreciable scattering amplitude;
(b) the dimensions of the singlecrystal regions giving appreciable scattering amplitudes are small so that the `shape transform' regions of scattering power around the reciprocallattice points are relatively large;
(c) the spread of wavelengths is small (10^{−5} or less, with no whiteradiation background) and the degree of collimation is better (10^{−4} to 10^{−6}) than for conventional Xray sources.
As a consequence of these factors, singlecrystal diffraction patterns may show many simultaneous reflections, representing almostplanar sections of reciprocal space, and may show fine structure or intensity variations reflecting the crystal dimensions and shape.

(2) Kinematical diffractionpattern intensities are calculated in a manner analogous to that for Xrays except that
(a) no polarization factor is included because of the smallangle scattering conditions;
(b) integration over regions of scattering power around reciprocallattice points cannot be assumed unless appropriate experimental conditions are ensured.
For a thin, flat, lamellar crystal of thickness H, the observed intensity is where is the excitation error for the h reflection and Ω is the unitcell volume.
For a singlecrystal diffraction pattern obtained by rotating a crystal or from a uniformly bent crystal or for a mosaic crystal with a uniform distribution of orientations, the intensity is where is the crystal volume and is the latticeplane spacing. For a polycrystalline sample of randomly oriented small crystals, the intensity per unit length of the diffraction ring is where is the multiplicity factor for the h reflection and L is the camera length, or the distance from the specimen to the detector plane. The special cases of `oblique texture' patterns from powder patterns having preferred orientations are treated in IT C (2004, Section 4.3.5
).
(3) Twobeam dynamical diffraction formulae: complex potentials including absorption. In the twobeam dynamical diffraction approximation, the intensities of the directly transmitted and diffracted beams for transmission through a crystal of thickness H, in the absence of absorption, are where is the extinction distance, , and where is the deviation from the Bragg angle.
For the case that , with the incident beam at the Bragg angle, this reduces to the simple Pendellösung expression
The effects on the elastic Bragg scattering amplitudes of the inelastic or diffuse scattering may be introduced by adding an outofphase component to the structure amplitudes, so that for a centrosymmetric crystal, becomes complex by addition of an imaginary component. Alternatively, an absorption function , having Fourier coefficients , may be postulated so that is replaced by . The are known as phenomenological absorption coefficients and their validity in manybeam diffraction has been demonstrated by, for example, Rez (1978).
The magnitudes depend on the nature of the experiment and the extent to which the various inelastically or diffusely scattered electrons are included in the measurements being made. If measurements are made of purely elastic scattering intensities for Bragg reflections or of image intensity variations due to the interaction of the sharp Bragg reflections only, the main contributions to the absorption coefficients are as follows (Radi, 1970):
(a) from plasmon and singleelectron excitations, is of the order of 0.1 V_{0} and , for , is negligibly small;
(b) from thermal diffuse scattering; is of the order of 0.1 V_{h} and decreasing more slowly than with scattering angle.
Including absorption effects in (2.5.2.26) for the case gives The Borrmann effect is not very pronounced for electrons because , but can be important for the imaging of defects in thick crystals (Hirsch et al., 1965; Hashimoto et al., 1961).
Attempts to obtain analytical solutions for the dynamical diffraction equations for more than two beams have met with few successes. There are some situations of high symmetry, with incident beams in exact zoneaxis orientations, for which the manybeam solution can closely approach equivalent two or threebeam behaviour (Fukuhara, 1966). Explicit solutions for the threebeam case, which displays some aspects of manybeam character, have been obtained (Gjønnes & Høier, 1971; Hurley & Moodie, 1980).



Electron optics. Electrons may be focused by use of axially symmetric magnetic fields produced by electromagnetic lenses. The focal length of such a lens used as a projector lens (focal points outside the lens field) is given bywhere is the relativistically corrected accelerating voltage and is the z component of the magnetic field. An expression in terms of experimental constants was given by Liebman (1955) as where is a constant, NI is the number of ampere turns of the lens winding, S is the length of the gap between the magnet pole pieces and D is the bore of the pole pieces.
Lenses of this type have irreducible aberrations, the most important of which for the paraxial conditions of electron microscopy is the thirdorder spherical aberration, coefficient , giving a variation of focal length of for a beam at an angle α to the axis. Chromatic aberration, coefficient , gives a spread of focal lengths for variations and of the accelerating voltage and lens currents, respectively.
The objective lens of an electron microscope is the critical lens for the determination of image resolution and contrast. The action of this lens in a conventional transmission electron microscope (TEM) is described by use of the Abbe theory for coherent incident illumination transmitted through the object to produce a wavefunction (see Fig. 2.5.2.2).

Figure 2.5.2.2
 top  pdf  Diagram representing the critical components of a conventional transmission electron microscope (TEM) and a scanning transmission electron microscope (STEM). For the TEM, electrons from a source A illuminate the specimen and the objective lens forms an image of the transmitted electrons on the image plane, B. For the STEM, a source at B is imaged by the objective lens to form a small probe on the specimen and some part of the transmitted beam is collected by a detector at A.

The amplitude distribution in the back focal plane of the objective lens is written where is the Fourier transform of and T(u, v) is the transfer function of the lens, consisting of an aperture function and a phase function exp where the phase perturbation due to lens defocus and aberrations is usually approximated as and u, v are the reciprocalspace variables related to the scattering angles , by
The image amplitude distribution, referred to the object coordinates, is given by Fourier transform of (2.5.2.31) as where , given by Fourier transform of , is the spread function. The image intensity is then
In practice the coherent imaging theory provides a good approximation but limitations of the coherence of the illumination have appreciable effects under highresolution imaging conditions.
The variation of focal lengths according to (2.5.2.30) is described by a function . Illumination from a finite incoherent source gives a distribution of incidentbeam angles . Then the image intensity is found by integrating incoherently over and :where denotes the Fouriertransform operation.
In the scanning transmission electron microscope (STEM), the objective lens focuses a small bright source of electrons on the object and directly transmitted or scattered electrons are detected to form an image as the incident beam is scanned over the object (see Fig. 2.5.2.2). Ideally the image amplitude can be related to that of the conventional transmission electron microscope by use of the `reciprocity relationship' which refers to point sources and detectors for scalar radiation in scalar fields with elastic scattering processes only. It may be stated: `The amplitude at a point B due to a point source at A is identical to that which would be produced at A for the identical source placed at B'.
For an axial point source, the amplitude distribution produced by the objective lens on the specimen is If this is translated by the scan to X, Y, the transmitted wave is
The amplitude on the plane of observation following the specimen is then and the image signal produced by a detector having a sensitivity function H(u, v) is If H(u, v) represents a small detector, approximated by a delta function, this becomes which is identical to the result (2.5.2.35) for a plane incident wave in the conventional transmission electron microscope.
(a) The weakphaseobject approximation. For sufficiently thin objects, the effect of the object on the incidentbeam amplitude may be represented by the transmission function (2.5.2.16) given by the phaseobject approximation. If the fluctuations, , about the mean value of the projected potential are sufficiently small so that , it is possible to use the weakphaseobject approximation (WPOA) where is referred to the average value, . The assumption that only firstorder terms in need be considered is the equivalent of a singlescattering, or kinematical, approximation applied to the twodimensional function, the projected potential of (2.5.2.16). From (2.5.2.42), the image intensity (2.5.2.35) becomes where the spread function s(xy) is the Fourier transform of the imaginary part of T(uv), namely .
The optimum imaging condition is then found, following Scherzer (1949), by specifying that the defocus should be such that is close to unity for as large a range of as possible. This is so for a negative defocus such that decreases to a minimum of about before increasing to zero and higher as a result of the fourthorder term of (2.5.2.33) (see Fig. 2.5.2.3). This optimum, `Scherzer defocus' value is given by or
The resolution limit is then taken as corresponding to the value of when becomes zero, before it begins to oscillate rapidly with U. The resolution limit is then For example, for mm and Å (200 keV), Å.
Within the limits of the WPOA, the image intensity can be written simply for a number of other imaging modes in terms of the Fourier transforms and of the real and imaginary parts of the objectivelens transfer function , where r and u are twodimensional vectors in real and reciprocal space, respectively.
For darkfield TEM images, obtained by introducing a central stop to block out the central beam in the diffraction pattern in the backfocal plane of the objective lens, Here, as in (2.5.2.42), should be taken to imply the difference from the mean potential value, .
For brightfield STEM imaging with a very small detector placed axially in the central beam of the diffraction pattern (2.5.2.39) on the detector plane, the intensity, from (2.5.2.41), is given by (2.5.2.43).
For a finite axially symmetric detector, described by , the image intensity is where is the Fourier transform of (Cowley & Au, 1978).
For STEM with an annular darkfield detector which collects all electrons scattered outside the central spot of the diffraction pattern in the detector plane, it can be shown that, to a good approximation (valid except near the resolution limit) Since is the intensity distribution of the electron probe incident on the specimen, (2.5.2.48) is equivalent to the incoherent imaging of the function .
Within the range of validity of the WPOA or, in general, whenever the zero beam of the diffraction pattern is very much stronger than any diffracted beam, the general expression (2.5.2.36) for the modifications of image intensities due to limited coherence may be conveniently approximated. The effect of integrating over the variables , may be represented by multiplying the transfer function T(u, v) by socalled `envelope functions' which involve the Fourier transforms of the functions and .
For example, if is approximated by a Gaussian of width (at of the maximum) centred at and is a circular aperture function the transfer function for coherent radiation is multiplied by where
(b) The projected chargedensity approximation. For very thin specimens composed of moderately heavy atoms, the WPOA is inadequate. Within the region of validity of the phaseobject approximation (POA), more complicated relations analagous to (2.5.2.43) to (2.5.2.47) may be written. A simpler expression may be obtained by use of the twodimensional form of Poisson's equation, relating the projected potential distribution to the projected chargedensity distribution . This is the PCDA (projected chargedensity approximation) (Cowley & Moodie, 1960),
This is valid for sufficiently small values of the defocus , provided that the effects of the spherical aberration may be neglected, i.e. for image resolutions not too close to the Scherzer resolution limit (Lynch et al., 1975). The function includes contributions from both the positive atomic nuclei and the negative electron clouds. For underfocus ( negative), single atoms give dark spots in the image. The contrast reverses with defocus.

(a) Introduction. It follows from (2.5.2.43) and (2.5.2.42) that, within the severe limitations of validity of the WPOA or the PCDA, images of very thin crystals, viewed with the incident beam parallel to a principal axis, will show dark spots at the positions of rows of atoms parallel to the incident beam. Provided that the resolution limit is less than the projected distances between atom rows (1–3 Å), the projection of the crystal structure may be seen directly.
In practice, theoretical and experimental results suggest that images may give a direct, although nonlinear, representation of the projected potential or chargedensity distribution for thicknesses much greater than the thicknesses for validity of these approximations, e.g. for thicknesses which may be 50 to 100 Å for 100 keV electrons for 3 Å resolutions and which increase for comparable resolutions at higher voltage but decrease with improved resolutions.
The use of highresolution imaging as a means for determining the structures of crystals and for investigating the form of the defects in crystals in terms of the arrangement of the atoms has become a widely used and important branch of crystallography with applications in many areas of solidstate science. It must be emphasized, however, that image intensities are strongly dependent on the crystal thickness and orientation and also on the instrumental parameters (defocus, aberrations, alignment etc.). It is only when all of these parameters are correctly adjusted to lie within strictly defined limits that interpretation of images in terms of atom positions can be attempted [see IT C (2004, Section 4.3.8
)].
Reliable interpretations of highresolution images of crystals (`crystal structure images') may be made, under even the most favourable circumstances, only by the comparison of observed image intensities with intensities calculated by use of an adequate approximation to manybeam dynamical diffraction theory [see IT C (2004, Section 4.3.6
)]. Most calculations for moderate or large unit cells are currently made by the multislice method based on formulation of the dynamical diffraction theory due to Cowley & Moodie (1957). For smaller unit cells, the matrix method based on the Bethe (1928) formulation is also frequently used (Hirsch et al., 1965).
(b) Fourier images. For periodic objects in general, and crystals in particular, the amplitudes of the diffracted waves in the back focal plane are given from (2.5.2.31) by For rectangular unit cells of the projected unit cell, the vector h has components and . Then the set of amplitudes (2.5.2.34), and hence the image intensities, will be identical for two different sets of defocus and spherical aberration values and if, for an integer N, i.e.
This relationship is satisfied for all h, k if is an integer and and where m, n are integers (Kuwabara, 1978). The relationship for is an expression of the Fourier image phenomenon, namely that for a planewave incidence, the intensity distribution for the image of a periodic object repeats periodically with defocus (Cowley & Moodie, 1960). Hence it is often necessary to define the defocus value by observation of a nonperiodic component of the specimen such as a crystal edge (Spence et al., 1977).
For the special case , the image intensity is also reproduced exactly for except that in this case the image is translated by a distance parallel to each of the axes.

(1) Definition and measurement. The `resolution' of an electron microscope or, more correctly, the `least resolvable distance', is usually defined by reference to the transfer function for the coherent imaging of a weak phase object under the Scherzer optimum defocus condition (2.5.2.44). The resolution figure is taken as the inverse of the U value for which first crosses the axis and is given, as in (2.5.2.45), by
It is assumed that an objective aperture is used to eliminate the contribution to the image for U values greater than the first zero crossing, since for these contributions the relative phases are distorted by the rapid oscillations of and the corresponding detail of the image is not directly interpretable in terms of the projection of the potential distribution of the object.
The resolution of the microscope in practice may be limited by the incoherent factors which have the effect of multiplying the WPOA transfer function by envelope functions as in (2.5.2.49).
The resolution, as defined above, and the effects of the envelope functions may be determined by Fourier transform of the image of a suitable thin, weakly scattering amorphous specimen. The Fouriertransform operation may be carried out by use of an optical diffractometer. A more satisfactory practice is to digitize the image directly by use of a twodimensional detector system in the microscope or from a photographic recording, and perform the Fourier transform numerically.
For the optical diffractometer method, the intensity distribution obtained is given from (2.5.2.43) as a radially symmetric function of U, where is the product of the envelope functions.
In deriving (2.5.2.54) it has been assumed that:
(a) the WPOA applies;
(b) the optical transmission function of the photographic record is linearly related to the image intensity, ;
(c) the diffraction intensity is a radially symmetric, smoothly varying function such as is normally produced by a sufficiently large area of the image of an amorphous material;
(d) there is no astigmatism present and no drift of the specimen; either of these factors would remove the radial symmetry.
From the form of (2.5.2.54) and a preknowledge of , the zero crossings of and the form of may be deduced. Analysis of a throughfocus series of images provides more complete and reliable information.

(2) Detail on a scale much smaller than the resolution of the electron microscope, as defined above, is commonly seen in electron micrographs, especially for crystalline samples. For example, lattice fringes, having the periodicity of the crystal lattice planes, with spacings as small as 0.6 Å in one direction, have been observed using a microscope having a resolution of about 2.5 Å (Matsuda et al., 1978), and twodimensionally periodic images showing detail on the scale of 0.5 to 1 Å have been observed with a similar microscope (Hashimoto et al., 1977).
Such observations are possible because
(a) for periodic objects the diffraction amplitude in (2.5.2.31) is a set of delta functions which may be multiplied by the corresponding values of the transfer function that will allow strong interference effects between the diffracted beams and the zero beam, or between different diffracted beams;
(b) the envelope functions for the WPOA, arising from incoherent imaging effects, do not apply for strongly scattering crystals; the more general expression (2.5.2.36) provides that the incoherent imaging factors will have much less effect on the interference of some sets of diffracted beams.
The observation of finely spaced lattice fringes provides a measure of some important factors affecting the microscope performance, such as the presence of mechanical vibrations, electrical interference or thermal drift of the specimen. A measure of the fineness of the detail observable in this type of image may therefore be taken as a measure of `instrumental resolution'.


Currently most electrondiffraction patterns are obtained in conjunction with images, in electron microscopes of one form or another, as follows.
(a) Selectedarea electrondiffraction (SAED) patterns are obtained by using intermediate and projector lenses to form an image of the diffraction pattern in the backfocal plane of the objective lens (Fig. 2.5.2.2). The area of the specimen from which the diffraction pattern is obtained is defined by inserting an aperture in the image plane of the objective lens. For parallel illumination of the specimen, sharp diffraction spots are produced by perfect crystals.
A limitation to the area of the specimen from which the diffraction pattern can be obtained is imposed by the spherical aberration of the objective lens. For a diffracted beam scattered through an angle α, the spread of positions in the object for which the diffracted beam passes through a small axial aperture in the image plane is , e.g. for mm, rad (10,0,0 reflection from gold for 100 keV electrons), Å, so that a selectedarea diameter of less than about 2000 Å is not feasible. For higher voltages, the minimum selectedarea diameter decreases with if the usual assumption is made that increases for highervoltage microscopes so that is a constant.
(b) Convergentbeam electrondiffraction (CBED) patterns are obtained when an incident convergent beam is focused on the specimen, as in an STEM instrument or an STEM attachment for a conventional TEM instrument.
For a large, effectively incoherent source, such as a conventional hotfilament electron gun, the intensities are added for each incidentbeam direction. The resulting CBED pattern has an intensity distributionwhere is the Fourier transform of the exit wave at the specimen for an incidentbeam direction .
(c) Coherent illumination from a small bright source such as a field emission gun may be focused on the specimen to give an electron probe having an intensity distribution and a diameter equal to the STEM darkfield image resolution [equation (2.5.2.47)] of a few Å. The intensity distribution of the resulting microdiffraction pattern is then where is the Fourier transform of the exit wave at the specimen. Interference occurs between waves scattered from the various incidentbeam directions. The diffraction pattern is thus an inline hologram as envisaged by Gabor (1949).
(d) Diffraction patterns may be obtained by using an optical diffractometer (or computer) to produce the Fourier transform squared of a small selected region of a recorded image. The optical diffractionpattern intensity obtained under the ideal conditions specified under equation (2.5.2.54) is given, in the case of weak phase objects, by or, more generally, by where is the Fourier transform of the wavefunction at the exit face of the specimen and c is a constant depending on the characteristics of the photographic recording medium.
