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. 276-345
https://doi.org/10.1107/97809553602060000558 Chapter 2.5. Electron diffraction and electron microscopy in structure determinationaArizona State University, Box 871504, Department of Physics and Astronomy, Tempe, AZ 85287-1504, USA, bSchool of Physics, University of Melbourne, Parkville, Australia 3052, cInstitute of Crystallography, Academy of Sciences of Russia, Leninsky prospekt 59, Moscow B-117333, Russia, dInstitute of Ore Mineralogy (IGEM), Academy of Sciences of Russia, Staromonetny 35, 109017 Moscow, Russia, and eExxonMobil Research and Engineering Co., 1545 Route 22 East, Clinton Township, Annandale, New Jersey 08801, USA This chapter deals with electron diffraction and imaging in the transmission geometry at kilovolt and higher beam energies for the purposes of crystal structure determination. Techniques such as low-energy electron diffraction, gas electron diffraction and reflection electron diffraction are not considered. Following a brief review of the history of the subject and its relationship to X-ray diffraction, the fundamental theory is presented by J. M. Cowley in Section 2.5.2. Unlike X-rays, which diffract from the electron density in a crystal, electrons are scattered elastically by the Coulomb potential, which is related to the density through Poisson's equation and includes the nuclear contribution. Electrons are much more strongly scattered than X-rays and have much smaller wavelengths. (This leads to a very ‚flat’ Ewald sphere and the simultaneous excitation of many Bragg beams.) Electron sources have comparable brightness to third-generation synchrotrons fitted with an undulator, so that Bragg intensities are very high. Because of their limited penetration, samples for transmission electron microscopes (TEMs) usually have thicknesses of less than a micron, while for interpretable atomic resolution images, thicknesses of tens of nanometres are preferred. Multiple scattering complicates TEM image interpretation and diffraction for larger thicknesses. The current (2006) spatial resolution of the best aberration-corrected TEM instruments is about 0.8 nm. Microdiffraction patterns may be obtained using a beam width of sub-nanometre dimensions, while the analysis of characteristic X-rays excited by such a probe provides chemical information. For organic monolayers (such as proteins), the amount of structural information obtained per unit of radiation damage (which is a strong function of resolution) exceeds that obtainable by X-ray diffraction. Cowley's treatment includes sections on electron scattering factors, Bethe's 1928 multiple scattering theory, Born's series, sign conventions, two-beam dynamical theory and single-scattering theory. This is followed by the theory of electron microscope imaging at high resolution, including the weak-phase object and other approximations, and the Scherzer focusing condition. The section ends with treatments of atomic resolution imaging in crystals and the factors which limit it, with Fourier or Talbot self-imaging, and with a brief discussion of coherent nanodiffraction. Section 2.5.3 by P. Goodman describes the convergent-beam electron diffraction (CBED) method used for space-group determination. This information is obtained in the microdiffraction mode, allowing nanocrystals to be analysed. A cone of illumination is used, which broadens the Bragg spots into discs, whose internal symmetries are analysed together with the symmetry of the whole pattern of discs. This section reviews the history of the subject, and provides tables which allow a systematic deduction of point and space groups from these data in most cases. Multiple scattering effects, in particular, can be used to distinguish between centrosymmetric and noncentrosymmetric crystals, while cancellation along symmetry-related multiple scattering paths allows identification of screw and glide elements. Examples from the literature are cited and sources of useful multiple-scattering software for simulating CBED patterns are noted. Section 2.5.4 by B. K. Vainshtein and B. B. Zvyagin is devoted to the use of transmission electron diffraction patterns to solve the structures of thin crystal structures using the kinematic theory. Texture and polycrystalline patterns are considered, as is the relationship to X-ray work. Section 2.5.5 by B. K. Vainshtein again outlines the theory of high-resolution electron imaging, extending this to include image processing, image cross-correlation and alignment, and image filtering and enhancement. Section 2.5.6 by B. K. Vainshtein discusses the problem of three-dimensional image reconstruction from projections, widely used in cryo-electron microscopy of proteins, and increasingly also now in materials science. The real-space methods of Radon, the method of back-projection, iterative methods and reciprocal-space techniques are described in full. Section 2.5.7 by D. L. Dorset summarizes solutions to the phase problem which may be applied to electron diffraction data. These include many of those currently in use for X-ray diffraction, including Patterson maps, direct methods and trial-and-error search techniques. Much of the section concerns electron diffraction data from thin organic films, analysed using the three-phase invariants of the direct methods approach to phasing. The tangent formula is reviewed and useful software is suggested. Density modification and maximum entropy methods are outlined. Because of the strong dependence of multiply scattered electron diffraction intensities on the thickness of the sample, and on local orientation changes (effects which are not accounted for in the structure analysis of the data), the thinnest possible samples must be used, and phasing measures are sought which are robust with respect to multiple scattering perturbations. Conversely, the ability to obtain data from nanometre-sized regions can greatly assist the effort to obtain high-quality perfect-crystal data free of defects, bending or thickness variation. Keywords: electron crystallography; electron diffraction; electron microscopy; EDSA; electron-diffraction structure analysis; kinematical diffraction; weak phase objects; crystal structure imaging; image reconstruction; convergent-beam electron diffraction; CBED; space-group determination; zone-axis patterns; Friedel's law; in-disc symmetries; symmetry elements; crystal defects; twins; Ewald sphere; optimal defocus; non-classical crystallography; image enhancement; three-dimensional reconstruction; back-projection method; Fourier transformation; discretization; direct phase determination; phase invariant sums; tangent formula; density modification; convolution techniques; maximum entropy; maximum likelihood; absorption in electron diffraction; image resolution. |
2.5.1. Foreword1
Given that electrons have wave properties and the wavelengths lie in a suitable range, the diffraction of electrons by matter is completely analogous to the diffraction of X-rays. While for X-rays the scattering function is the electron-density distribution, for electrons it is the potential distribution which is similarly peaked at the atomic sites. Hence, in principle, electron diffraction may be used as the basis for crystal structure determination. In practice it is used much less widely than X-ray diffraction for the determination of crystal structures but is receiving increasing attention as a means for obtaining structural information not readily accessible with X-ray- or neutron-diffraction techniques.
Electrons having wavelengths comparable with those of the X-rays commonly used in diffraction experiments have energies of the order of 100 eV. For such electrons, the interactions with matter are so strong that they can penetrate only a few layers of atoms on the surfaces of solids. They are used extensively for the study of surface structures by low-energy electron diffraction (LEED) and associated techniques. These techniques are not covered in this series of volumes, which include the principles and practice of only those diffraction and imaging techniques making use of high-energy electrons, having energies in the range of 20 keV to 1 MeV or more, in transmission through thin specimens.
For the most commonly used energy ranges of high-energy electrons, 100 to 400 keV, the wavelengths are about 50 times smaller than for X-rays. Hence the scattering angles are much smaller, of the order of rad, the recording geometry is relatively simple and the diffraction pattern represents, to a useful first approximation, a planar section of reciprocal space.
The elastic scattering of electrons by atoms is several orders of magnitude greater than for X-rays. This fact has profound consequences, which in some cases are highly favourable and in other cases are serious hindrances to structure analysis work. On the one hand it implies that electron-diffraction patterns can be obtained from very small single-crystal regions having thicknesses equal to only a few layers of atoms and, with recently developed techniques, having diameters equivalent to only a few interatomic distances. Hence single-crystal patterns can be obtained from microcrystalline phases.
However, the strong scattering of electrons implies that the simple kinematical single-scattering approximation, on which most X-ray diffraction structure analysis is based, fails for electrons except for very thin crystals composed of light-atom materials. Strong dynamical diffraction effects occur for crystals which may be 100 Å thick, or less for heavy-atom materials. As a consequence, the theory of dynamical diffraction for electrons has been well developed, particularly for the particular special diffracting conditions relevant to the transmission of fast electrons (see Chapter 5.2 ), and observations of dynamical diffraction effects are commonly made and quantitatively interpreted. The possibility has thus arisen of using the observation of dynamical diffraction effects as the basis for obtaining crystal structure information. The fact that dynamical diffraction is dependent on the relative phases of the diffracted waves then implies that relative phase information can be deduced from the diffraction intensities and the limitations of kinematical diffraction, such as Friedel's law, do not apply. The most immediately practicable method for making use of this possibility is convergent-beam electron diffraction (CBED) as described in Section 2.5.3.
A further important factor, determining the methods for observing electron diffraction, is that, being charged particles, electrons can be focused by electromagnetic lenses. The irreducible aberrations of cylindrical magnetic lenses have, to date, limited the resolution of electron microscopes to the extent that the least resolvable distances (or `resolutions') are about 100 times the electron wavelength. However, with microscopes having a resolution of better than 2 Å it is possible to distinguish the individual rows of atoms, parallel to the incident electron beam, in the principal orientations of many crystalline phases. Thus `structure images' can be obtained, sometimes showing direct representation of projections of crystal structures [see IT C (2004), Section 4.3.8 ]. However, the complications of dynamical scattering and of the coherent imaging processes are such that the image intensities vary strongly with crystal thickness and tilt, and with the defocus or other parameters of the imaging system, making the interpretation of images difficult except in special circumstances. Fortunately, computer programs are readily available whereby image intensities can be calculated for model structures [see IT C (2004), Section 4.3.6 ] Hence the means exist for deriving the projection of the structure if only by a process of trial and error and not, as would be desirable, from a direct interpretation of the observations.
The accuracy with which the projection of a structure can be deduced from an image, or series of images, improves as the resolution of the microscope improves but is not at all comparable with the accuracy attainable with X-ray diffraction methods. A particular virtue of high-resolution electron microscopy as a structural tool is that it may give information on individual small regions of the sample. Structures can be determined of `phases' existing over distances of only a few unit cells and the defects and local disorders can be examined, one by one.
The observation of electron-diffraction patterns forms an essential part of the technique of structure imaging in high-resolution electron microscopy, because the diffraction patterns are used to align the crystals to appropriate axial orientations. More generally, for all electron microscopy of crystalline materials the image interpretation depends on knowledge of the diffraction conditions. Fortunately, the diffraction pattern and image of any specimen region can be obtained in rapid succession by a simple switching of lens currents. The ready comparison of the image and diffraction data has become an essential component of the electron microscopy of crystalline materials but has also been of fundamental importance for the development of electron-diffraction theory and techniques.
The individual specimen regions giving single-crystal electron-diffraction patterns are, with few exceptions, so small that they can be seen only by use of an electron microscope. Hence, historically, it was only after electron microscopes were commonly available that the direct correlations of diffraction intensities with crystal size and shape could be made, and a proper basis was available for the development of the adequate dynamical diffraction theory.
For the complete description of a diffraction pattern or image intensities obtained with electrons, it is necessary to include the effects of inelastic scattering as well as elastic scattering. In contrast to the X-ray diffraction case, the inelastic scattering does not produce just a broad and generally negligible background. The average energy loss for an inelastic scattering event is about 20 eV, which is small compared with the energy of about 100 keV for the incident electrons. The inelastically scattered electrons have a narrow angular distribution and are diffracted in much the same way as the incident or elastically scattered electrons in a crystal. They therefore produce a highly modulated contribution to the diffraction pattern, strongly peaked about the Bragg spot positions (see Chapter 4.3 ). Also, as a result of the inelastic scattering processes, including thermal diffuse scattering, an effective absorption function must be added in the calculation of intensities for elastically scattered electrons.
The inelastic scattering processes in themselves give information about the specimen in that they provide a measure of the excitations of both the valence-shell and the inner-shell electrons of the solid. The inner-shell electron excitations are characteristic of the type of atom, so that microanalysis of small volumes of specimen material (a few hundreds or thousands of atoms) may be achieved by detecting either the energy losses of the transmitted electrons or the emission of the characteristic X-ray [see IT C (2004), Section 4.3.4 ].
An adverse effect of the inelastic scattering processes, however, is that the transfer of energy to the specimen material results in radiation damage; this is a serious limitation of the application of electron-scattering methods to radiation-sensitive materials such as organic, biological and many inorganic compounds. The amount of radiation damage increases rapidly as the amount of information per unit volume, derived from the elastic scattering, is increased, i.e. as the microscope resolution is improved or as the specimen volume irradiated during a diffraction experiment is decreased. At the current limits of microscopic resolution, radiation damage is a significant factor even for the radiation-resistant materials such as semiconductors and alloys.
In the historical development of electron-diffraction techniques the progress has depended to an important extent on the level of understanding of the dynamical diffraction processes and this understanding has followed, to a considerable degree, from the availability of electron microscopes. For the first 20 years of the development, with few exceptions, the lack of a precise knowledge of the specimen morphology meant that diffraction intensities were influenced to an unpredictable degree by dynamical scattering and the impression grew that electron-diffraction intensities could not meaningfully be interpreted.
It was the group in the Soviet Union, led initially by Dr Z. G. Pinsker and later by Dr B. K. Vainshtein and others, which showed that patterns from thin layers of a powder of microcrystals could be interpreted reliably by use of the kinematical approximation. The averaging over crystal orientation reduced the dynamical diffraction effects to the extent that practical structure analysis was feasible. The development of the techniques of using films of crystallites having strongly preferred orientations, to give patterns somewhat analogous to the X-ray rotation patterns, provided the basis for the collection of three-dimensional diffraction data on which many structure analyses have been based [see Section 2.5.4 and IT C (2004), Section 4.3.5 ].
In recent years improvements in the techniques of specimen preparation and in the knowledge of the conditions under which dynamical diffraction effects become significant have allowed progress to be made with the use of high-energy electron diffraction patterns from thin single crystals for crystal structure analysis. Particularly for crystals of light-atom materials, including biological and organic compounds, the methods of structure analysis developed for X-ray diffraction, including the direct methods (see Section 2.5.7), have been successfully applied in an increasing number of cases. Often it is possible to deduce some structural information from high-resolution electron-microscope images and this information may be combined with that from the diffraction intensities to assist the structure analysis process [see IT C (2004), Section 4.3.8.8 ].
The determination of crystal symmetry by use of CBED (Section 2.5.3) and the accurate determination of structure amplitudes by use of methods depending on the observation of dynamical diffraction effects [IT C (2004), Section 4.3.7 ] came later, after the information on morphologies of crystals, and the precision electron optics associated with electron microscopes, became available.
In spite of the problem of radiation damage, a great deal of progress has been made in the study of organic and biological materials by electron-scattering methods. In some respects these materials are very favourable because, with only light atoms present, the scattering from thin films can be treated using the kinematical approximation without serious error. Because of the problem of radiation damage, however, special techniques have been evolved to maximize the information on the required structural aspects with minimum irradiation of the specimen. Image-processing techniques have been evolved to take advantage of the redundancy of information from a periodic structure and the means have been devised for combining information from multiple images and diffraction data to reconstruct specimen structure in three dimensions. These techniques are outlined in Sections 2.5.5 and 2.5.6. They are based essentially on the application of the kinematical approximation and have been used very effectively within that limitation.
For most inorganic materials the complications of many-beam dynamical diffraction processes prevent the direct application of these techniques of image analysis, which depend on having a linear relationship between the image intensity and the value of the projected potential distribution of the sample. The smaller sensitivities to radiation damage can, to some extent, remove the need for the application of such methods by allowing direct visualization of structure with ultra-high-resolution images and the use of microdiffraction techniques.
2.5.2. Electron diffraction and electron microscopy1
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 high-energy electrons (in the energy range of to eV) in transmission through thin samples of crystalline solids and the derivation of information on crystal structures from diffraction patterns and high-resolution 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 X-rays, it is apparent that the angular range for strong scattering of electrons will be of the order of rad. Only under special circumstances, usually involving multiple elastic and inelastic scattering from very thick specimens, are scattering angles of more than rad of importance.
The strength of the interaction of electrons with matter is greater than that of X-rays by two or three orders of magnitude. The single-scattering, 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 Å, rather than Å 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 field-emission electron guns, it is possible to obtain a current of A in an electron beam of diameter 10 Å or less with a beam divergence of less than rad, i.e. a current density of 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 atoms. On the other hand, electron beams may readily be collimated to better than 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 X-ray or neutron diffraction. The mathematical expressions and the procedures are much the same. However, there are peculiarities of the electron-diffraction case which should be noted.
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 structure-factor expression and (2) defining the forward propagating free-space wavefunction with a positive exponent.
The second of these alternatives is the one which has been adopted in most solid-state and quantum-mechanical 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 structure-factor 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 self-consistent way, even in those cases where, as indicated above, some of the signs appear to have no effect on one particular conclusion.
electron interaction constant ; (relativistic) electron mass; electron wavelength; (magnitude of) electron charge; Planck's constant; ; volume of the unit cell; continuous reciprocal-space vector, components u, v; discrete reciprocal-space coordinate; crystal potential averaged along beam direction (positive); slice thickness; absorption potential [positive; typically ]; defocus (defined as negative for underfocus); spherical aberration coefficient; excitation error relative to the incident-beam direction and defined as negative when the point h lies outside the Ewald sphere; atomic scattering factor for electrons, , related to the atomic scattering factor for X-rays, , by the Mott formula . Fourier transform of periodic slice transmission function.
|
The forward-scattering approximation to the many-beam dynamical diffraction theory outlined in Chapter 5.2 provides the basis for the calculation of diffraction intensities and electron-microscope 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.
|
|
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 by where 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 third-order 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).
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 reciprocal-space 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 high-resolution 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 incident-beam angles . Then the image intensity is found by integrating incoherently over and : where denotes the Fourier-transform 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.
|
|
|
Currently most electron-diffraction patterns are obtained in conjunction with images, in electron microscopes of one form or another, as follows.
|
2.5.3. Space-group determination by convergent-beam electron diffraction2
Convergent-beam electron diffraction, originating in the experiments of Kossel and Möllenstedt (Kossel & Möllenstedt, 1938) has been established over the past two decades as a powerful technique for the determination of space group in inorganic materials, with particular application when only microscopic samples are available. Relatively recently, with the introduction of the analytical electron microscope, this technique – abbreviated as CBED – has become available as a routine, so that there is now a considerable accumulation of data from a wide range of materials. A significant extension of the technique in recent times has been the introduction of LACBED (large-angle CBED) by Tanaka & Terauchi (1985). This technique allows an extensive angular range of single diffraction orders to be recorded and, although this method cannot be used for microdiffraction (since it requires an extensive single-crystal area), new LACBED applications appear regularly, particularly in the field of semiconductor research (see Section 2.5.3.6).
The CBED method relies essentially on two basic properties of transmission electron diffraction, namely the radical departure from Friedel's law and the formation of characteristic extinction bands within space-group-forbidden reflections. Departure from Friedel's law in electron diffraction was first noted experimentally by Miyake & Uyeda (1950). The prediction of space-group-forbidden bands (within space-group-forbidden reflections) by Cowley & Moodie (1959), on the other hand, was one of the first successes of N-beam theory. A detailed explanation was later given by Gjønnes & Moodie (1965). These are known variously as `GM' bands (Tanaka et al., 1983), or more simply and definitively as `GS' (glide–screw) bands (this section). These extinctions have a close parallel with space-group extinctions in X-ray diffraction, with the reservation that only screw axes of order two are accurately extinctive under N-beam conditions. This arises from the property that only those operations which lead to identical projections of the asymmetric unit can have N-beam dynamical symmetries (Cowley et al., 1961).
Additionally, CBED from perfect crystals produces high-order defect lines in the zero-order pattern, analogous to the defect Kikuchi lines of inelastic scattering, which provide a sensitive measurement of unit-cell parameters (Jones et al., 1977; Fraser et al., 1985; Tanaka & Terauchi, 1985).
The significant differences between X-ray and electron diffraction, which may be exploited in analysis, arise as a consequence of a much stronger interaction in the case of electrons (Section 2.5.2). Hence, thin, approximately parallel-sided crystal regions must be used in high-energy (100 kV–1 MV) electron transmission work, so that diffraction is produced from crystals effectively infinitely periodic in only two dimensions, leading to the relaxation of three-dimensional diffraction conditions known as `excitation error' (Chapter 5.2 ). Also, there is the ability in CBED to obtain data from microscopic crystal regions of around 50 Å in diameter, with corresponding exposure times of several seconds, allowing a survey of a material to be carried out in a relatively short time.
In contrast, single-crystal X-ray diffraction provides much more limited symmetry information in a direct fashion [although statistical analysis of intensities (Wilson, 1949) will considerably supplement this information], but correspondingly gives much more direct three-dimensional geometric data, including the determination of unit-cell parameters and three-dimensional extinctions.
The relative strengths and weaknesses of the two techniques make it useful where possible to collect both convergent-beam and X-ray single-crystal data in a combined study. However, all parameters can be obtained from convergent-beam and electron-diffraction data, even if in a somewhat less direct form, making possible space-group determination from microscopic crystals and microscopic regions of polygranular material. Several reviews of the subject are available (Tanaka, 1994; Steeds & Vincent, 1983; Steeds, 1979). In addition, an atlas of characteristic CBED patterns for direct phase identification of metal alloys has been published (Mansfield, 1984), and it is likely that this type of procedure, allowing N-beam analysis by comparison with standard simulations, will be expanded in the near future.
Symmetry analysis is necessarily tied to examination of patterns near relevant zone axes, since the most intense N-beam interaction occurs amongst the zero-layer zone-axis reflections, with in addition a limited degree of upper-layer (higher-order Laue zone) interaction. There will generally be several useful zone axes accessible for a given parallel-sided single crystal, with the regions between axes being of little use for symmetry analysis. Only one such zone axis can be parallel to a crystal surface normal, and a microcrystal is usually chosen at least initially to have this as the principal symmetry axis. Other zone axes from that crystal may suffer mild symmetry degradation because the N-beam lattice component (`excitation error' extension) will not have the symmetry of the structure (Goodman, 1974; Eades et al., 1983).
Upper-layer interactions, responsible for imparting three-dimensional information to the zero layer, are of two types: the first arising from `overlap' of dynamic shape transforms and causing smoothly varying modulations of the zero-layer reflections, and the second, caused by direct interactions with the upper-layer, or higher-order Laue zone lines, leading to a sharply defined fine-line structure. These latter interactions are especially useful in increasing the accuracy of space-group determination (Tanaka et al., 1983), and may be enhanced by the use of low-temperature specimen stages. The presence of these defect lines in convergent-beam discs, occurring especially in low-symmetry zone-axis patterns, allows symmetry elements to be related to the three-dimensional structure (Section 2.5.3.5; Fig. 2.5.3.4c).
To the extent that such three-dimensional effects can be ignored or are absent in the zero-layer pattern the projection approximation (Chapter 5.2 ) can be applied. This situation most commonly occurs in zone-axis patterns taken from relatively thin crystals and provides a useful starting point for many analyses, by identifying the projected symmetry.
Convergent-beam diffraction symmetries are those of Schrödinger's equation, i.e. of crystal potential, plus the diffracting electron. The appropriate equation is given in Section 2.5.2 [equation (2.5.2.6)] and Chapter 5.2 [equation (5.2.2.1)] in terms of the real-space wavefunction ψ. The symmetry elements of the crystal responsible for generating pattern symmetries may be conveniently classified as of two types (I and II) as follows.
|
A minimal summary of basic theoretical points, otherwise found in Chapter 5.2 and numerous referenced articles, is given here.
For a specific zero-layer diffraction order the incident and diffracted vectors are and . Then the three-dimensional vector has the pattern-space projection, . The point gives the symmetrical Bragg condition for the associated diffraction disc, and is identifiable with the angular deviation of from the vertical z axis in three-dimensional space (see Fig. 2.5.3.1). also defines the symmetry centre within the two-dimensional disc diagram (Fig. 2.5.3.2); namely, the intersection of the lines S and G, given by the trace of excitation error, , and the perpendicular line directed towards the reciprocal-space origin, respectively. To be definitive it is necessary to index diffracted amplitudes relating to a fixed crystal thickness and wavelength, with both crystallographic and momentum coordinates, as , to handle the continuous variation of (for a particular diffraction order), with angles of incidence as determined by , and registered in the diffraction plane as the projection of .
Vector diagram in semi-reciprocal space, using Ewald-sphere constructions to show the `incident', `reciprocity' and `reciprocity × centrosymmetry' sets of vectors. Dashed lines connect the full vectors to their projections in the plane of observation. |
Reciprocity was introduced into the subject of electron diffraction in stages, the essential theoretical basis, through Schrödinger's equation, being given by Bilhorn et al. (1964), and the N-beam diffraction applications being derived successively by von Laue (1935), Cowley (1969), Pogany & Turner (1968), Moodie (1972), Buxton et al. (1976), and Gunning & Goodman (1992).
Reciprocity represents a reverse-incidence configuration reached with the reversed wavevectors and , so that the scattering vector is unchanged, but is changed in sign and hence reversed (Moodie, 1972). The reciprocity equation, is valid independently of crystal symmetry, but cannot contribute symmetry to the pattern unless a crystal-inverting symmetry element is present (since belongs to a reversed wavevector). The simplest case is centrosymmetry, which permits the right-hand side of (2.5.3.1) to be complex-conjugated giving the useful CBED pattern equation Since K is common to both sides there is a point-by-point identity between the related distributions, separated by 2g (the distance between g and reflections). This invites an obvious analogy with Friedel's law, , with the reservation that (2.5.3.2) holds only for centrosymmetric crystals. This condition (2.5.3.2) constitutes what has become known as the ±H symmetry and, incidentally, is the only reciprocity-induced symmetry so general as to not depend upon a disc symmetry-point or line, nor on a particular zone axis (i.e. it is not a point symmetry but a translational symmetry of the pattern intensity).
|
Horizontal glides, a′, n′ (diperiodic, primed notation), generate zero-layer absent rows, or centring, rather than GS bands (see Fig. 2.5.3.3). This is an example of the projection approximation in its most universally held form, i.e. in application to absences. Other examples of this are: (a) appearance of both G and S extinction bands near their intersection irrespective of whether glide or screw axes are involved; and (b) suppression of the influence of vertical, non-primitive translations with respect to observations in the zero layer. It is generally assumed as a working rule that the zero-layer or ZOLZ pattern will have the rotational symmetry of the point-group component of the vertical screw axis (so that ). Elements included in Table 2.5.3.1 on this pretext are given in parentheses. However, the presence of rather than 2 ( rather than 3 etc.) should be detectable as a departure from accurate twofold symmetry in the first-order-Laue-zone (FOLZ) reflection circle (depicted in Fig. 2.5.3.3). This has been observed in the cubic structure of Ba2Fe2O5Cl2, permitting the space groups I23 and to be distinguished (Schwartzman et al., 1996). A summary of all the symmetry components described in this section is given diagrammatically in Table 2.5.3.2.
Diagrammatic representation of the influence of non-symmorphic elements: (i) Alternate rows of the zero-layer pattern are absent owing to the horizontal glide plane. The pattern is indexed as for an `a' glide; the alternative indices (in parentheses) apply for a `b' glide. (ii) GS bands are shown along the central row of the zero layer, for odd-order reflections. |
The following guidelines, the result of accumulated experience from several laboratories, are given in an experimentally based sequence, and approximately in order of value and reliability.
|
These results are illustrated in Table 2.5.3.2 and by actual examples in Section 2.5.3.5.
Space groups may very well be identified using CBED patterns from an understanding of the diffraction properties of real-space symmetry elements, displayed for example in Table 2.5.3.2. It is, however, of great assistance to have the symmetries tabulated in reciprocal space, to allow direct comparison with the pattern symmetries.
There are three generally useful ways in which this can be done, and these are set out in Tables 2.5.3.3 to 2.5.3.5. The simplest of these is by means of point group, following the procedures of Buxton et al. (1976). Next, the CBED pattern symmetries can be listed as diperiodic groups which are space groups in two dimensions, allowing identification with a restricted set of three-dimensional space groups (Goodman, 1984b). Finally, the dynamic extinctions (GS bands and zero-layer absences) can be listed for each non-symmorphic space group, together with the diffraction conditions for their observation (Tanaka et al., 1983; Tanaka & Terauchi, 1985). Descriptions for these tables are given below.
|
|
|
|
The list of examples given here must necessarily be regarded as unsatisfactory considering the vastness of the subject, although some attempt has been made to choose a diverse range of problems which will illustrate the principles involved. Some particular aspects, however, need further mention.
One of these concerns the problem of examining large-unit-cell materials with a high diffraction-pattern density. This limits the possible convergence angle, if overlap is to be avoided, and leaves numerous but featureless discs [for example Goodman (1984b)]. Technical advances which have been made to overcome this problem include the beam-rocking technique (Eades, 1980) and LACBED (Tanaka et al., 1980), both of which are reviewed by Tanaka & Terauchi (1985) and Eades et al. (1983). The disadvantage of these latter methods is that they both require a significantly larger area of specimen than does the conventional technique, and it may be that more sophisticated methods of handling the crowded conventional patterns are still needed.
Next, the matter of accuracy must be considered. There are two aspects of the subject where this is of concern. Firstly, there is a very definite limit to the sensitivity with which symmetry can be detected. In a simple structure of medium-light atoms, displacements of say 0.1 Å or less from a pseudomirror plane could easily be overlooked. An important aspect of CBED analysis, not mentioned above, is the N-beam computation of patterns which is required when something approaching a refinement (in the context of electron diffraction) is being attempted. Although this quantitative aspect has a long history [for example see Johnson (1972)], it has only recently been incorporated into symmetry studies as a routine (Creek & Spargo, 1985; Tanaka, 1994). Multi-slice programs which have been developed to produce computer-simulated pattern output are available (Section 2.5.3.8).
Next there is concern as to the allocation of a space group to structures which microscopically have a much lower symmetry (Goodman et al., 1984). This arises because the volume sampled by the electron probe necessarily contains a large number of unit cells. Reliable microscopic interpretation of certain nonstoichiometric materials requires that investigations be accompanied by high-resolution microscopy. Frequently (especially in mineralogical samples), nonstoichiometry implies that a space group exists only on average, and that the concept of absolute symmetry elements is inapplicable.
From earlier and concluding remarks it will be clear that combined X-ray/CBED and CBED/electron-microscopy studies of inorganic materials represents the standard ideal approach to space-group analysis at present; given this approach, all the space-group problems of classical crystallography appear soluble. As has been noted earlier, it is important that HREM be considered jointly with CBED in determining space group by electron crystallography, and that only by this joint study can the so-called `phase problem' be completely overcome. The example of the space-group pairs and has already been cited. Using CBED, it might be expected that FOLZ lines would show a break from twofold symmetry with the incident beam aligned with a axis. However, a direct distinction should be made apparent from high-resolution electron micrographs. Other less clear-cut cases occur where the HREM images allow a space-group distinction to be made between possible space groups of the same arithmetic class, especially when only one morphology is readily obtained (e.g. , , ).
The slightly more subtle problem of distinguishing enantiomorphic space-group pairs can be solved by one of two approaches: either the crystal must be rotated around an axis by a known amount to obtain two projections, or the required three-dimensional phase information can be deduced from specific three-beam-interaction data. This problem is part of the more general problem of solving handedness in an asymmetric structure, and is discussed in detail by Johnson & Preston (1994).
|
2.5.4. Electron-diffraction structure analysis (EDSA)4
Electron-diffraction structure analysis (EDSA) (Vainshtein, 1964) based on electron diffraction (Pinsker, 1953) 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 ) 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 X-rays, , 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 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 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: 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 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).
Electron-diffraction 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 Gas-Phase 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 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 . 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).
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; Pinsker, 1953; Zvyagin, 1967).
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 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 non-zero 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 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; 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 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 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 electron-beam 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 oblique-texture (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 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 around the axis 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. ) 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 oblique-angled lattices In the general case of oblique-angled 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 X-ray 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 X-ray 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 unit-cell potential (see Chapter 1.2 ), 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: 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 single-crystal films the integral intensity of reflection is for textures Here 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 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 two-wave 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 inner-crystalline 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 wedge-shaped crystal (Cowley, 1981), or from rocking curves.
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 5.2.2 ). The structural determination is based on experimental sets of values or (Vainshtein, 1964).
The trial-and-error 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 heavy-atom 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 three-dimensional 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 electron-density Fourier series in X-ray 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 half-width than the electron-density peaks . 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: while in X-ray diffraction . 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 (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 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 X-ray diffraction, EDSA and neutron-diffraction data. Thus for the atomic amplitudes where Z is the nuclear charge and the X-ray atomic scattering amplitude, and for structure amplitudes where is the X-ray 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 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 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.
2.5.5. Image reconstruction 5
In many fields of physical measurements, instrumental and informative techniques, including electron microscopy and computational or analogue methods for processing and transforming signals from objects investigated, find a wide application in obtaining the most accurate structural data. The signal may be radiation from an object, or radiation transmitted through the object, or reflected by it, which is transformed and recorded by a detector.
The image is the two-dimensional signal on the observation plane recorded from the whole three-dimensional volume of the object, or from its surface, which provides information on its structure. In an object this information may change owing to transformation of the scattered wave inside an instrument. The real image is composed of and noise from signal disturbances:
Image-reconstruction methods are aimed at obtaining the most accurate information on the structure of the object; they are subdivided into two types (Picture Processing and Digital Filtering, 1975; Rozenfeld, 1969):
|
These two methods may be used separately or in combination.
The image should be represented in the form convenient for perception and analysis, e.g. in digital form, in lines of equal density, in points of different density, in half-tones or colour form and using, if necessary, a change or reversal of contrast.
Reconstructed images may be used for the three-dimensional reconstruction of the spatial structure of an object, e.g. of the density distribution in it (see Section 2.5.6).
This section is connected with an application of the methods of image processing in transmission electron microscopy (TEM). In TEM (see Section 2.5.2), the source-emitted electrons are transmitted through an object and, with the aid of a system of lenses, form a two-dimensional image subject to processing.
Another possibility for obtaining information on the structure of an object is structural analysis with the aid of electron diffraction – EDSA. This method makes use of information in reciprocal space – observation and measurement of electron-diffraction patterns and calculation from them of a two-dimensional projection or three-dimensional structure of an object using the Fourier synthesis. To do this, one has to find the relative phases of the scattered beams.
The wavefunction of an electron-microscopic image is written as Here is the incident plane wave. When the wave is transmitted through an object, it interacts with the electrostatic potential [ is the three-dimensional vector in the space of the object]; this process is described by the Schrödinger equation (Section 2.5.2.1). As a result, on the exit surface of an object the wave takes the form where q is the transmission function and x is the two-dimensional vector . The diffraction of the wave is described by the two-dimensional Fourier operator:
Here, we assume the initial wave amplitude to be equal to unity and the initial phase to be zero, so that , which defines, in this case, the wavefunction in the back focal plane of an objective lens with the reciprocal-space coordinates . The function Q is modified in reciprocal space by the lens transfer function . The scattered wave transformation into an image is described by the inverse Fourier operator .
The process of the diffraction , as seen from (2.5.5.1), is the same in both TEM and EDSA. Thus, in TEM under the lens actions the image formation from a diffraction pattern takes place with an account of the phases, but these phases are modified by the objective-lens transfer function. In EDSA, on the other hand, there is no distorting action of the transfer function and the `image' is obtained by computing the operation .
The computation of projections, images and Fourier transformation is made by discretization of two-dimensional functions on a two-dimensional network of points – pixels in real space and in reciprocal space .
The intensity distribution of an electron wave in the image plane depends not only on the coherent and inelastic scattering, but also on the instrumental functions. The electron wave transmitted through an object interacts with the electrostatic potential which is produced by the nuclei charges and the electronic shells of the atoms. The scattering and absorption of electrons depend on the structure and thickness of a specimen, and the atomic numbers of the atoms of which it is composed. If an object with the three-dimensional distribution of potential is sufficiently thin, then the interaction of a plane electron wave with it can be described as the interaction with a two-dimensional distribution of potential projection , where b is the specimen thickness. It should be noted that, unlike the three-dimensional function of potential with dimension , the two-dimensional function of potential projection has the potential-length dimension which, formally, coincides with the charge dimension. The transmission function, in the general case, has the form (2.5.2.42), and for weak phase objects the approximation is valid.
In the back focal plane of the objective lens the wave has the form where is the Scherzer phase function (Scherzer, 1949) of an objective lens (Fig. 2.5.5.1), is the aperture function, the spherical aberration coefficient, and Δf the defocus value [(2.5.2.32) –(2.5.2.35)].
The bright-field image intensity (in object coordinates) is where . The phase function (2.5.5.7) depends on defocus, and for a weak phase object (Cowley, 1981) where , which includes only an imaginary part of function (2.5.5.6). While selecting defocus in such a way that under the Scherzer defocus conditions [(2.5.2.44), (2.5.2.45)] , one could obtain In this very simple case the image reflects directly the structure of the object – the two-dimensional distribution of the projection of the potential convoluted with the spread function . In this case, no image restoration is necessary. Contrast reversal may be achieved by a change of defocus.
At high resolution, this method enables one to obtain an image of projections of the atomic structure of crystals and defects in the atomic arrangement – vacancies, replacements by foreign atoms, amorphous structures and so on; at resolution worse than atomic one obtains images of dislocations as continuous lines, inserted phases, inclusions etc. (Cowley, 1981). It is also possible to obtain images of thin biological crystals, individual molecules, biological macromolecules and their associations.
Image restoration. In the case just considered (2.5.5.10), the projection of potential , convoluted with the spread function, can be directly observed. In the general case (2.5.5.9), when the aperture becomes larger, the contribution to image formation is made by large values of spatial frequencies U, in which the function sin χ oscillates, changing its sign. Naturally, this distorts the image just in the region of appropriate high resolution. However, if one knows the form of the function sin χ (2.5.5.7), the true function can be restored.
This could be carried out experimentally if one were to place in the back focal plane of an objective lens a zone plate transmitting only one-sign regions of sin χ (Hoppe, 1971). In this case, the information on is partly lost, but not distorted. To perform such a filtration in an electron microscope is a rather complicated task.
Another method is used (Erickson & Klug, 1971). It consists of a Fourier transformation of the measured intensity distribution TQ (2.5.5.6) and division of this transform, according to (2.5.5.7a,b), by the phase function sin χ. This gives Then, the new Fourier transformation yields (in the weak-phase-object approximation) the true distribution The function sin χ depending on defocus Δf should be known to perform this procedure. The transfer function can also be found from an electron micrograph (Thon, 1966). It manifests itself in a circular image intensity modulation of an amorphous substrate or, if the specimen is crystalline, in the `noise' component of the image. The analogue method (optical Fourier transformation for obtaining the image ) can be used (optical diffraction, see below); digitization and Fourier transformation can also be applied (Hoppe et al., 1973).
The thin crystalline specimen implies that in the back focal objective lens plane the discrete kinematic amplitudes are arranged and, by the above method, they are corrected and released from phase distortions introduced by the function sin χ (see below) (Unwin & Henderson, 1975).
For the three-dimensional reconstruction (see Section 2.5.6) it is necessary to have the projections of potential of the specimen tilted at different angles α to the beam direction (normal beam incidence corresponds to ). In this case, the defocus Δf changes linearly with increase of the distance l of specimen points from the rotation axis . Following the above procedure for passing on to reciprocal space and correction of sin χ, one can find (Henderson & Unwin, 1975).
Elastic interaction of an incident wave with a weak phase object is defined on its exit surface by the distribution of potential projection ; however, in the general case, the electron scattering amplitude is a complex one (Glauber & Schomaker, 1953). In such a way, the image itself has the phase and amplitude contrast. This may be taken into account if one considers not only the potential projection , but also the `imaginary potential' which describes phenomenologically the absorption in thin specimens. Then, instead of (2.5.5.5), the wave on the exit surface of a specimen can be written as and in the back focal plane if and Usually, μ is small, but it can, nevertheless, make a certain contribution to an image. In a sufficiently good linear approximation, it may be assumed that the real part cos χ of the phase function (2.5.5.7a) affects , while , as we know, is under the action of the imaginary part sin χ.
Thus, instead of (2.5.5.6), one can write and as the result, instead of (2.5.5.10),
The functions and can be separated by object imaging using the through-focus series method. In this case, using the Fourier transformation, one passes from the intensity distribution (2.5.5.15) in real space to reciprocal space. Now, at two different defocus values and [(2.5.5.6), (2.5.5.7a,b)] the values and can be found from the two linear equations (2.5.5.14). Using the inverse Fourier transformation, one can pass on again to real space which gives and (Schiske, 1968). In practice, it is possible to use several through-focus series and to solve a set of equations by the least-squares method.
Another method for processing takes into account the simultaneous presence of noise and transfer function zeros (Kirkland et al., 1980). In this method the space frequencies corresponding to small values of the transfer function modulus are suppressed, while the regions where such a modulus is large are found to be reinforced.
When the specimen thickness exceeds a certain critical value (50–100 Å), the kinematic approximation does not hold true and the scattering is dynamic. This means that on the exit surface of a specimen the wave is not defined as yet by the projection of potential (2.5.5.3), but one has to take into account the interaction of the incident wave and of all the secondary waves arising in the whole volume of a specimen.
The dynamic scattering calculation can be made by various methods. One is the multi-slice (or phase-grating) method based on a recurrent application of formulae (2.5.5.3) for n thin layers thick, and successive construction of the transmission functions (2.5.5.4), phase functions , and propagation function (Cowley & Moodie, 1957).
Another method – the scattering matrix method – is based on the solution of equations of the dynamic theory (Chapter 5.2 ). The emerging wave on the exit surface of a crystal is then found to diffract and experience the transfer function action [(2.5.5.6), (2.5.5.7a,b)].
The dynamic scattering in crystals may be interpreted using Bloch waves: It turns out that only a few (bound and valence Bloch waves) have strong excitation amplitudes. Depending on the thickness of a crystal, only one of these waves or their linear combinations (Kambe, 1982) emerges on the exit surface. An electron-microscopic image can be interpreted, at certain thicknesses, as an image of one of these waves [with a correction for the transfer function action (2.5.5.6), (2.5.5.7a,b)]; in this case, the identical images repeat with increasing thickness, while, at a certain thickness, the contrast reversal can be observed. Only the first Bloch wave which arises at small thickness, and also repeats with increasing thickness, corresponds to the projection of potential , i.e. the atom projection distribution in a thin crystal layer.
An image of other Bloch waves is defined by the function , but their maxima or minima do not coincide, in the general case, with the atomic positions and cannot be interpreted as the projection of potential. It is difficult to reconstruct from these images, especially when the crystal is not ideal and contains imperfections. In these cases one resorts to computer modelling of images at different thicknesses and defocus values, and to comparison with an experimentally observed pattern.
The imaging can be performed directly in an electron microscope not by a photo plate, but using fast-response detectors with digitized intensity output on line. The computer contains the necessary algorithms for Fourier transformation, image calculation, transfer function computing, averaging, and correction for the observed and calculated data. This makes possible the interpretation of the pattern observed directly in experiment (Herrmann et al., 1980).
The real electron-microscope image is subdivided into two components: The main of these, , is a two-dimensional image of the `ideal' object obtained in an electron microscope with instrumental functions inherent to it. However, in the process of object imaging and transfer of this information to the detector there are various sources of noise. In an electron microscope, these arise owing to emission-current and accelerating-voltage fluctuations, lens-supplying current (temporal fluctuations), or mechanical instabilities in a device, specimen or detector (spatial shifts). The two-dimensional detector (e.g. a photographic plate) has structural inhomogeneities affecting a response to the signal. In addition, the specimen is also unstable; during preparation or imaging it may change owing to chemical or some other transformations in its structure, thermal effects and so on. Biological specimens scatter electrons very weakly and their natural state is moist, while in the electron-microscope column they are under vacuum conditions. The methods of staining (negative or positive), e.g. of introducing into specimens substances containing heavy atoms, as well as the freeze-etching method, somewhat distort the structure of a specimen. Another source of structure perturbation is radiation damage, which can be eliminated at small radiation doses or by using the cryogenic technique. The structure of stained specimens is affected by stain graininess. We assume that all the deviations of a specimen image from the `ideal' image are included in the noise term . The substrate may also be inhomogeneous. All kinds of perturbations cannot be separated and they appear on an electron microscope image as the full noise content .
The image enhancement involves maximum noise suppression and hence the most accurate separation of a useful signal from the real image (2.5.5.1). At the signal/noise ratio such a separation appears to be rather complicated. But in some cases the real image reflects the structure sufficiently well, e.g. during the atomic structure imaging of some crystals . In other cases, especially of biological specimen imaging, the noise N distorts substantially the image, . Here one should use the methods of enhancement. This problem is usually solved by the methods of statistical processing of sets of images . If one assumes that the informative signal is always the same, then the noise error may be reduced.
The image enhancement methods are subdivided into two classes:
These methods can be used separately or in combination. The enhancement can be applied to both the original and the restored images; there are also methods of simultaneous restoration and enhancement.
The image can be enhanced by analogue (mainly optical and photographic) methods or by computational methods for processing digitized functions in real and reciprocal space.
The cases where the image has translational symmetry, rotational symmetry, and where the image is asymmetric will be considered.
Periodic images. An image of the crystal structure with atomic or molecular resolution may be brought to self-alignment by a shift by a and b periods in a structure projection. This can be performed photographically by printing the shifted image on the same photographic paper or, vice versa, by shifting the paper (McLachlan, 1958).
The Fourier filtration method for a periodic image with noise N is based on the fact that in Fourier space the components and are separated. Let us carry out the Fourier transformation of the periodic signal with the periods a, b and noise N: The left part of (2.5.5.18) represents the Fourier coefficients distributed discretely with periods and in the plane . This is the two-dimensional reciprocal lattice. The right-hand side of (2.5.5.18) is the Fourier transform distributed continuously in the plane. Thus these parts are separated. Let us `cut out' from distribution (2.5.5.18) only values using the `window' function . The window should match each of the real peaks which, owing to the finite dimensions of the initial periodic image, are not points, as this is written in an idealized form in (2.5.5.18) with the aid of δ functions. In reality, the `windows' may be squares of about , in size, or a circle. Performing the Fourier transformation of product (2.5.5.18) without , and set of windows , we obtain: the periodic component without the background, . The zero coefficient in (2.5.5.19) should be decreased, since it is due, in part, to the noise. When the window w is sufficiently small, in (2.5.5.19) represents the periodic distribution (average over all the unit cells of the projection) included in (2.5.5.18). Nevertheless, some error from noise in an image does exist, since with we also introduced into the inverse Fourier transformation the background transform values which are within the `windows'.
This approach is realized by an analogue method [optical diffraction and filtering of electron micrographs in a laser beam (Klug & Berger, 1964)] and can also be carried out by computing.
As an example, Fig. 2.5.5.2(b) shows an electron micrograph of the periodic structure of a two-dimensional protein crystal, while Fig. 2.5.5.2(c) represents optical diffraction from this layer. In order to dissect the aperiodic component in a diffraction plane, according to the scheme in Fig. 2.5.5.2(a), one places a mask with windows covering reciprocal-lattice points. After such a filtration, only the component makes a contribution during the image formation by means of a lens, while the component diffracted by the background is delayed. As a result, an optical pattern of the periodic structure is obtained (Fig. 2.5.5.2d).
(a) Diagram of an optical diffractometer. D is the object (an electron micrograph), is the diffraction plane and a mask that transmits only , is the plane of the (filtered) image; (b) an electron micrograph of a crystalline layer of the protein phosphorylase b; (c) its optical diffraction pattern (the circles correspond to the windows in the mask that transmits only the diffracted beams from the periodic component of the image); (d) the filtered image. Parts (b)–(d) are based on the article by Kiselev et al. (1971). |
Optical diffractometry also assists in determining the parameters of a two-dimensional lattice and its symmetry.
Using the same method, one can separate the superimposed images of two-dimensional structures with different periodicity and in different orientation, the images of the `near' and `far' sides of tubular periodic structures with monomolecular walls (Klug & DeRosier, 1966; Kiselev et al., 1971), and so on.
Computer filtering involves measuring the image optical density , digitization, and Fourier transformation (Crowther & Amos, 1971). The sampling distance usually corresponds to one-third of the image resolution. When periodic weak phase objects are investigated, the transformation (2.5.5.18) yields the Fourier coefficients. If necessary, we can immediately make corrections in them using the microscope transfer function according to (2.5.5.6), (2.5.5.7a,b) and (2.5.5.11a), and thereby obtain the true kinematic amplitudes . The inverse transformation (2.5.5.16) gives a projection of the structure (Unwin & Henderson, 1975; Henderson & Unwin, 1975).
Sometimes, an observed image is `noised' by the to a great extent. Then, one may combine data on real and reciprocal space to construct a sufficiently accurate image. In this case, the electron-diffraction pattern is measured and structure-factor moduli from diffraction reflection intensities are obtained: At the same time, the structure factors are calculated from the processed structure projection image by means of the Fourier transformation. However, owing to poor image quality we take from these data only the values of phases since they are less sensitive to scattering density distortions than the moduli, and construct the Fourier synthesis
Here the possibilities of combining various methods open up, e.g. for obtaining the structure-factor moduli from X-ray diffraction, and phases from electron microscopy, and so on (Gurskaya et al., 1971).
Images with point symmetry. If a projection of an object (and consequently, the object itself) has a rotational N-fold axis of symmetry, the structure coincides with itself on rotation through the angle . If the image is rotated through arbitrary angles and is aligned photographically with the initial image, then the best density coincidence will take place at a rotation through which defines N. The pattern averaging over all the rotations will give the enhanced structure image with an times reduced background (Markham et al., 1963).
Rotational filtering can be performed on the basis of the Fourier expansion of an image in polar coordinates over the angles (Crowther & Amos, 1971). The integral over the radius from azimuthal components gives their power where a is the maximum radius of the particle. A set forms a spectrum, the least common multiple N of strong peaks defining the N-fold symmetry. The two-dimensional reconstructed image of a particle with rotational symmetry is defined by the synthesis (2.5.5.24) with .
Asymmetric images. In this case, a set of images is processed by computational or analogue methods. The initial selection of images involves the fulfillment of the maximum similarity condition.
The averaging of n images in real space gives The signal/noise ratio on an average image is times enhanced.
The degree of similarity and accuracy of superposition of two images with an account both of translational and angular shifts is estimated by a cross-correlation function 3 of two selected images and (Frank, 1975, 1980). The value is the measure of image similarity, the x′ coordinate of the maximum indicates the shift of the images relative to each other. The first term of the resultant expression (2.5.5.26) is the cross-correlation function of noise-corrected images being compared, the second and third terms are approximately equal to zero, since the noise does not correlate with the signal; the last term is the autocorrelation function of the noise (Cramér, 1954; Frank, 1975, 1980).
The calculation of a correlation function is performed by means of Fourier transformation on the basis of the convolution theorem, since the Fourier transformation of the product of the Fourier transform of function and the conjugated Fourier transform function gives the cross-correlation function of the initial functions:
The probability density of samples for images has the form
Here J is the tentative image (as such, a certain `best' image can first be selected, while at the repeated cycle an average image is obtained), is the image investigated, σ is the standard deviation of the normal distribution of noises and the relative shift of the image. This function is called a likelihood function; it has maxima relative to the parameters , . The average image and dispersion are This method is called the maximum-likelihood method (Cramér, 1954; Kosykh et al., 1983).
It is convenient to carry out the image alignment, in turn, with respect to translational and angular coordinates. If we start with an angular alignment we first use autocorrelation functions or power spectra, which have the maximum and the symmetry centre at the origin of the coordinates. The angular correlation maximum gives the mutual angle of rotation of two images.
Then we carry out the translational alignment of rotationally aligned images using the translational correlation function (2.5.5.26) (Langer et al., 1970).
In the iteration alignment method, the images are first translationally aligned and then an angular shift is determined in image space in polar coordinates with the centre at the point of the best translational alignment. After the angular alignment the whole procedure may be repeated (Steinkilberg & Schramm, 1980).
The average image obtained may have false high-frequency components. They can be excluded by multiplying its Fourier components by some function and suppressing high-space frequencies, for instance by an `artificial temperature factor' .
For a set of similar images the Fourier filtration method can also be used (Ottensmeyer et al., 1977). To do this, one should prepare from these images an artificial `two-dimensional crystal', i.e. place them in the same orientation at the points of the two-dimensional lattice with periods a, b. The processing is then performed according to (2.5.5.18), (2.5.5.19); as a result one obtains with reduced background. Some translational and angular errors in the arrangement of the images at the artificial lattice points act as an artificial temperature factor. The method can be realized by computing or by optical diffraction.
2.5.6. Three-dimensional reconstruction 6
In electron microscopy we obtain a two-dimensional image – a projection of a three-dimensional object (Fig. 2.5.6.1): The projection direction is defined by a unit vector and the projection is formed on the plane x perpendicular to The set of various projections may be assigned by a discrete or continuous set of points on a unit sphere (Fig. 2.5.6.2). The function reflects the structure of an object, but gives information only on coordinates of points of its projected density. However, a set of projections makes it possible to reconstruct from them the three-dimensional (3D) distribution (Radon, 1917; DeRosier & Klug, 1968; Vainshtein et al., 1968; Crowther, DeRosier & Klug, 1970; Gordon et al., 1970; Vainshtein, 1971a; Ramachandran & Lakshminarayanan, 1971; Vainshtein & Orlov, 1972, 1974; Gilbert, 1972a; Herman, 1980). This is the task of the three-dimensional reconstruction of the structure of an object:
The projection sphere and projection of along τ onto the plane . The case represents orthoaxial projection. Points indicate a random distribution of τ. |
Besides electron microscopy, the methods of reconstruction of a structure from its projections are also widely used in various fields, e.g. in X-ray and NMR tomography, in radioastronomy, and in various other investigations of objects with the aid of penetrating, back-scattered or their own radiations (Bracewell, 1956; Deans, 1983; Mersereau & Oppenheim, 1974).
In the general case, the function (2.5.6.1) (the subscript indicates dimension) means the distribution of a certain scattering density in the object. The function is the two-dimensional projection density; one can also consider one-dimensional projections of two- (or three-) dimensional distributions. In electron microscopy, under certain experimental conditions, by functions and we mean the potential and the projection of the potential, respectively [the electron absorption function μ (see Section 2.5.4) may also be considered as `density']. Owing to a very large depth of focus and practical parallelism of the electron beam passing through an object, in electron microscopy the vector τ is the same over the whole area of the irradiated specimen – this is the case of parallel projection.
The 3D reconstruction (2.5.6.2) can be made in the real space of an object – the corresponding methods are called the methods of direct three-dimensional reconstruction (Radon, 1917; Vainshtein et al., 1968; Gordon et al., 1970; Vainshtein, 1971a; Ramachandran & Lakshminarayanan, 1971; Vainshtein & Orlov, 1972, 1974; Gilbert, 1972a).
On the other hand, three-dimensional reconstruction can be carried out using the Fourier transformation, i.e. by transition to reciprocal space. The Fourier reconstruction is based on the well known theorem: the Fourier transformation of projection of a three-dimensional object is the central (i.e. passing through the origin of reciprocal space) two-dimensional plane cross section of a three-dimensional transform perpendicular to the projection vector (DeRosier & Klug, 1968; Crowther, DeRosier & Klug, 1970; Bracewell, 1956). In Cartesian coordinates a three-dimensional transform is The transform of projection along z is In the general case (2.5.6.1) of projecting the plane along the vector τ Reconstruction with Fourier transformation involves transition from projections at various to cross sections , then to construction of the three-dimensional transform by means of interpolation between in reciprocal space, and transition by the inverse Fourier transformation to the three-dimensional distribution : Transition (2.5.6.2) or (2.5.6.6) from two-dimensional electron-microscope images (projections) to a three-dimensional structure allows one to consider the complex of methods of 3D reconstruction as three-dimensional electron microscopy. In this sense, electron microscopy is an analogue of methods of structure analysis of crystals and molecules providing their three-dimensional spatial structure. But in structure analysis with the use of X-rays, electrons, or neutrons the initial data are the data in reciprocal space in (2.5.6.6), while in electron microscopy this role is played by two-dimensional images [(2.5.6.2), (2.5.6.6)] in real space.
In electron microscopy the 3D reconstruction methods are, mainly, used for studying biological structures (symmetric or asymmetric associations of biomacromolecules), the quaternary structure of proteins, the structures of muscles, spherical and rod-like viruses, bacteriophages, and ribosomes.
An exact reconstruction is possible if there is a continuous set of projections corresponding to the motion of the vector over any continuous line connecting the opposite points on the unit sphere (Fig. 2.5.6.2). This is evidenced by the fact that, in this case, the cross sections which are perpendicular to τ in Fourier space (2.5.6.4) continuously fill the whole of its volume, i.e. give (2.5.6.3) and thereby determine .
In reality, we always have a discrete (but not continuous) set of projections . The set of is, practically, obtained by the rotation of the specimen under the beam through various angles (Hoppe & Typke, 1979) or by imaging of the objects which are randomly oriented on the substrate at different angles (Kam, 1980; Van Heel, 1984). If the object has symmetry, one of its projections is equivalent to a certain number of different projections.
The object is finite in space. For function and any of its projections there holds the normalization condition where Ω is the total `weight' of the object described by the density distribution . If one assumes that the density of an object is constant and that inside the object = constant = 1, and outside it , then Ω is the volume of an object. The volume of an object, say, of molecules, viruses and so on, is usually known from data on the density or molecular mass.
In practice, an important case is where all the projection directions are orthogonal to a certain straight line: (Fig. 2.5.6.3). Here the axis of rotation or the axis of symmetry of an object is perpendicular to an electron beam. Then the three-dimensional problem is reduced to the two-dimensional one, since each cross section is represented by its one-dimensional projections. The direction of vector τ is defined by the rotational angle ψ of a specimen: In this case, the reconstruction is carried out separately for each level : and the three-dimensional structure is obtained by superposition of layers (Vainshtein et al., 1968; Vainshtein, 1978).
In direct methods of reconstruction as well as in Fourier methods the space is represented as a discrete set of points on a two-dimensional net or on a three-dimensional lattice. It is sometimes expedient to use cylindrical or spherical coordinates. In two-dimensional reconstruction the one-dimensional projections are represented as a set of discrete values , at a certain spacing in . The reconstruction (2.5.6.9) is carried out over the discrete net with nodes . The net side A should exceed the diameter of an object ; the spacing . Then (2.5.6.8) transforms into the sum For oblique projections the above sum is taken over all the points within the strips of width a along the axis (Fig. 2.5.6.4).
The resolution δ of the reconstructed function depends on the number h of the available projections. At approximately uniform angular distribution of projections, and diameter equal to D, the resolution at reconstruction is estimated as The reconstruction resolution δ should be equal to or somewhat better than the instrumental resolution d of electron micrographs , the real resolution of the reconstructed structure being d. If the number of projections h is not sufficient, i.e. , then the resolution of the reconstructed structure is δ (Crowther, DeRosier & Klug, 1970; Vainshtein, 1978).
In electron microscopy the typical instrumental resolution d of biological macromolecules for stained specimens is about 20 Å; at the object with diameter Å the sufficient number h of projections is about 20. If the projections are not uniformly distributed in projection angles, the resolution decreases towards for such τ in which the number of projections is small.
Properties of projections of symmetric objects . If the object has an N-fold axis of rotation, its projection has the same symmetry. At orthoaxial projection perpendicular to the N-fold axis the projections which differ in angle at are identical: This means that one of its projections is equivalent to N projections. If we have h independent projections of such a structure, the real number of projections is hN (Vainshtein, 1978). For a structure with cylindrical symmetry one of its projections fully determines the three-dimensional structure.
Many biological objects possess helical symmetry – they transform into themselves by the screw displacement operation , where p is the number of packing units in the helical structure per q turns of the continuous helix. In addition, the helical structures may also have the axis of symmetry N defining the pitch of the helix. In this case, a single projection is equivalent to projections (Cochran et al., 1952).
Individual protein molecules are described by point groups of symmetry of type N or . Spherical viruses have icosahedral symmetry 532 with two-, three- and fivefold axes of symmetry. The relationship between vectors τ of projections is determined by the transformation matrix of the corresponding point group (Crowther, Amos et al., 1970).
Modelling. If several projections are available, and, especially, if the object is symmetric, one can, on the basis of spatial imagination, recreate approximately the three-dimensional model of the object under investigation. Then one can compare the projections of such a model with the observed projections, trying to draw them as near as possible. In early works on electron microscopy of biomolecules the tentative models of spatial structure were constructed in just this way; these models provide, in the case of the quaternary structure of protein molecules or the structure of viruses, schemes for the arrangement of protein subunits. Useful subsidiary information in this case can be obtained by the method of optical diffraction and filtration.
This method is also called the synthesis of projection functions. Let us consider a two-dimensional case and stretch along each one-dimensional projection (Fig. 2.5.6.5) by a certain length b; thus, we obtain the projection function Let us now superimpose h functions The continuous sum over the angles of projection synthesis is this is the convolution of the initial function with a rapidly falling function (Vainshtein, 1971b). In (2.5.6.15), the approximation for a discrete set of h projections is also written. Since the function approaches infinity at , the convolution with it will reproduce the initial function , but with some background B decreasing around each point according to the law . At orthoaxial projection the superposition of cross sections arranged in a pile gives the three-dimensional structure .
Radon operator. Radon (1917; see also Deans, 1983) gave the exact solution of the problem of reconstruction. However, his mathematical work was for a long time unknown to investigators engaged in reconstruction of a structure from images; only in the early 1970s did some authors obtain results analogous to Radon's (Ramachandran & Lakshminarayanan, 1971; Vainshtein & Orlov, 1972, 1974; Gilbert, 1972a).
The convolution in (2.5.6.15) may be eliminated using the Radon integral operator, which modifies projections by introducing around each point the negative values which annihilate on superposition the positive background values. The one-dimensional projection modified with the aid of the Radon operator has the form Now is calculated analogously to (2.5.6.14), not from the initial projections L but from the modified projection :
The reconstruction of high-symmetry structures, in particular helical ones, by the direct method is carried out from one projection making use of its equivalence to many projections. The Radon formula in discrete form can be obtained using the double Fourier transformation and convolution (Ramachandran & Lakshminarayanan, 1971).
These methods have been derived for the two-dimensional case; consequently, they can also be applied to three-dimensional reconstruction in the case of orthoaxial projection.
Let us discretize by a net of points ; then we can construct the system of equations (2.5.6.10).
When h projections are available the condition of unambiguous solution of system (2.5.6.10) is: . At we can, in practice, obtain sufficiently good results (Vainshtein, 1978).
Methods of reconstruction by iteration have also been derived that cause some initial distribution to approach one satisfying the condition that its projection will resemble the set . Let us assign on a discrete net as a zero-order approximation the uniform distribution of mean values (2.5.6.7) The projection of the qth approximation at the angle (used to account for discreteness) is .
The next approximation for each point jk is given in the method of `summation' by the formula where is the number of points in a strip of the projection . One cycle of iterations involves running around all of the angles (Gordon et al., 1970).
When carrying out iterations, we may take into account the contribution not only of the given projection, but also of all others. In this method the process of convergence improves. Some other iteration methods have been elaborated (Gordon & Herman, 1971; Gilbert, 1972b; Crowther & Klug, 1974; Gordon, 1974).
This method is based on the Fourier projection theorem [(2.5.6.3) –(2.5.6.5)]. The reconstruction is carried out according to scheme (2.5.6.6) (DeRosier & Klug, 1968; Crowther, DeRosier & Klug, 1970; Crowther, Amos et al. 1970; DeRosier & Moore, 1970; Orlov, 1975). The three-dimensional Fourier transform is found from a set of two-dimensional cross sections on the basis of the Whittaker–Shannon interpolation. If the object has helical symmetry (which often occurs in electron microscopy of biological objects, e.g. on investigating bacteriophage tails, muscle proteins) cylindrical coordinates are used. Diffraction from such structures with c periodicity and scattering density is defined by the Fourier–Bessel transform: The inverse transformation has the form so that and are the mutual Bessel transforms
Owing to helical symmetry, (2.5.6.22), (2.5.6.23) contain only those of the Bessel functions which satisfy the selection rule (Cochran et al., 1952) where N, q and p are the helix symmetry parameters, . Each layer l is practically determined by the single function with the lowest n; the contribution of other functions is neglected. Thus, the Fourier transformation of one projection of a helical structure, with an account of symmetry and phases, gives the three-dimensional transform (2.5.6.23). We can introduce into this transform the function of temperature-factor type filtering the `noise' from large spatial frequencies.
In the general case of 3D reconstruction from projections the projection vector τ occupies arbitrary positions on the projection sphere (Fig. 2.5.6.2). Then, as in (2.5.6.15), we can construct the three-dimensional spatial synthesis. To do this, let us transform the two-dimensional projections by extending them along τ as in (2.5.6.13) into three-dimensional projection functions .
Analogously to (2.5.6.15), such a three-dimensional synthesis is the integral over the hemisphere (Fig. 2.5.6.2) this is the convolution of the initial function with (Vainshtein, 1971b).
To obtain the exact reconstruction of we find, from each , the modified projection (Vainshtein & Orlov, 1974; Orlov, 1975)
By extending along τ we transform them into . Now the synthesis over the angles gives the three-dimensional function
The approximation for a discrete set of angles is written on the right. In this case we are not bound by the coaxial projection condition which endows the experiment with greater possibilities; the use of object symmetry also profits from this. To carry out the 3D reconstruction (2.5.6.25) or (2.5.6.27) one should know all three Euler's angles ψ, θ, α (Fig. 2.5.6.2).
The projection vectors should be distributed more or less uniformly over the sphere (Fig. 2.5.6.2). This can be achieved by using special goniometric devices.
Another possibility is the investigation of particles which, during the specimen preparation, are randomly oriented on the substrate. This, in particular, refers to asymmetric ribosomal particles. In this case the problem of determining these orientations arises.
The method of spatial correlation functions may be applied if a large number of projections with uniformly distributed projection directions is available (Kam, 1980). The space correlation function is the averaged characteristic of projections over all possible directions which is calculated from the initial projections or the corresponding sections of the Fourier transform. It can be used to find the coefficients of the object density function expansion over spherical harmonics, as well as to carry out the 3D reconstruction in spherical coordinates.
Another method (Van Heel, 1984) involves the statistical analysis of image types, subdivision of images into several classes and image averaging inside the classes. Then, if the object is rotated around some axis, the 3D reconstruction is carried out by the iteration method.
If such a specimen is inclined at a certain angle with respect to the beam, then the images of particles in the preferred orientation make a series of projections inclined at an angle β and having a random azimuth. The azimuthal rotation is determined from the image having zero inclination.
If particles on the substrate have a characteristic shape, they may acquire a preferable orientation with respect to the substrate, their azimuthal orientation α being random (Radermacher et al., 1987).
In the general case, the problem of determining the spatial orientations of randomly distributed identical three-dimensional particles with an unknown structure may be solved by measuring their two-dimensional projections (Fig. 2.5.6.1) if the number i of such projections is not less than three, (Vainshtein & Goncharov, 1986a,b; Goncharov et al. 1987; Goncharov, 1987). The direction of the vector along which the projection is obtained is set by the angle (Fig. 2.5.6.2).
The method is based on the analysis of one-dimensional projections of two-dimensional projections where α is the angle of the rotation about vector τ in the p plane.
Lemma 1. Any two projections and (Fig. 2.5.6.6) have common (identical) one-dimensional projections : Vectors and (Fig. 2.5.6.3) determine plane h in which they are both lying. Vector is normal to plane h and parallel to axis of the one-dimensional projection ; both and axes along which the projections and are constructed are perpendicular to .
The corresponding lemma in the Fourier space states:
Lemma 2. Any two plane transforms, and intersect along the straight line (Fig. 2.5.6.7); the one-dimensional transform is the transform of .
Thus in order to determine the orientations of a three-dimensional particle it is necessary either to use projections in real space or else to pass to the Fourier space (2.5.6.5).
Now consider real space. The projections are known and can be measured but angles of their rotation about vector (Fig. 2.5.6.8) are unknown and should be determined. Let us choose any two projections and and construct a set of one-dimensional projections and by varying angles and . In accordance with Lemma 1, there exists a one-dimensional projection, common for both and , which determines angles and along which and should be projected for obtaining the identical projection (Fig. 2.5.6.5). Comparing and and using the minimizing function it is possible to find such a common projection . (A similar consideration in Fourier space yields .)
Plane projections of a three-dimensional body. The systems of coordinates in planes (a) and (b) are chosen independently of one another. |
The mutual spatial orientations of any three non-coplanar projection vectors , , can be found from three different two-dimensional projections , and by comparing the following pairs of projections: and , and , and and , and by determining the corresponding , and . The determination of angles , and reduces to the construction of a trihedral angle formed by planes , and . Then the projections with the known can be complemented with other projections and the corresponding values of ω can be determined. Having a sufficient number of projections and knowing the orientations , it is possible to carry out the 3D reconstruction of the object [see (2.5.6.27); Orlov, 1975; Vainshtein & Goncharov, 1986a; Goncharov et al., 1987].
The concept of using experimental electron-diffraction intensities for quantitative crystal structure analyses has already been presented in Section 2.5.4. Another aspect of quantitative structure analysis, employing high-resolution images, has been presented in Sections 2.5.5 and 2.5.6. That is to say, electron micrographs can be regarded as an independent source of crystallographic phases.
Before direct methods (Chapter 2.2 ) were developed as the standard technique for structure determination in small-molecule X-ray crystallography, there were two principal approaches to solving the crystallographic phase problem. First, `trial and error' was used, finding some means to construct a reasonable model for the crystal structure a priori, e.g. by matching symmetry properties shared by the point group of the molecule or atomic cluster and the unit-cell space group. Secondly, the autocorrelation function of the crystal, known as the Patterson function (Chapter 2.3 ), was calculated (by the direct Fourier transform of the available intensity data) to locate salient interatomic vectors within the unit cell.
The same techniques had been used for electron-diffraction structure analysis (nowadays known as electron crystallography ). In fact, advocacy of the first method persists. Because of the perturbations of diffracted intensities by multiple-beam dynamical scattering (Chapter 5.2 ), it has often been suggested that trial and error be used to construct the scattering model for the unit crystal in order to test its convergence to observed data after simulation of the scattering events through the crystal. This indirect approach assumes that no information about the crystal structure can be obtained directly from observed intensity data. Under more favourable scattering conditions nearer to the kinematical approximation, i.e. for experimental data from thin crystals made up of light atoms, trial and error modelling, simultaneously minimizing an atom–atom nonbonded potential function with the crystallographic residual, has enjoyed widespread use in electron crystallography, especially for the determination of linear polymer structures (Brisse, 1989; Pérez & Chanzy, 1989).
Interpretation of Patterson maps has also been important for structure analysis in electron crystallography. Applications have been discussed by Vainshtein (1964), Zvyagin (1967) and Dorset (1994a). In face of the dynamical scattering effects for electron scattering from heavy-atom crystals realized later (e.g. Cowley & Moodie, 1959), attempts had also been made to modify this autocorrelation function by using a power series in to sharpen the peaks (Cowley, 1956). (Here , replacing the notation for the kinematical electron-diffraction structure factor employed in Section 2.5.4.) More recently, Vincent and co-workers have selected first-order-Laue-zone data from inorganics to minimize the effect of dynamical scattering on the interpretability of their Patterson maps (Vincent & Exelby, 1991, 1993; Vincent & Midgley, 1994). Vainshtein & Klechkovskaya (1993) have also reported use of the Patterson function to solve the crystal structure of a lead soap from texture electron-diffraction intensity data.
It is apparent that trial-and-error techniques are most appropriate for ab initio structure analysis when the underlying crystal structures are reasonably easy to model. The requisite positioning of molecular (or atomic) groups within the unit cell may be facilitated by finding atoms that fit a special symmetry position [see IT A (2005)]. Alternatively, it is helpful to know the molecular orientation within the unit cell (e.g. provided by the Patterson function) to allow the model to be positioned for a conformational or translational search. [Examples would include the polymer-structure analyses cited above, as well as the layer-packing analysis of some phospholipids (Dorset, 1987).] While attempts at ab initio modelling of three-dimensional crystal structures, by searching an n-dimensional parameter space and seeking a global internal energy minimum, has remained an active research area, most success so far seems to have been realized with the prediction of two-dimensional layers (Scaringe, 1992). In general, for complicated unit cells, determination of a structure by trial and error is very difficult unless adequate constraints can be placed on the search.
Although Patterson techniques have been very useful in electron crystallography, there are also inherent difficulties in their use, particularly for locating heavy atoms. As will be appreciated from comparison of scattering-factor tables for X-rays [IT C (2004) Chapter 6.1 ] with those for electrons, [IT C (2004) Chapter 4.3 ] the relative values of the electron form factors are more compressed with respect to atomic number than are those for X-ray scattering. As discussed in Chapter 2.3 , it is desirable that the ratio of summed scattering factor terms, , where Z is the scattering factor value at sin , be near unity. A practical comparison would be the value of r for copper (DL-alaninate) solved from electron-diffraction data by Vainshtein et al. (1971). For electron diffraction, compared to the value 2.36 for X-ray diffraction. Orientation of salient structural features, such as chains and rings, would be equally useful for light-atom moieties in electron or X-ray crystallography with Patterson techniques. As structures become more complicated, interpretation of Patterson maps becomes more and more difficult unless an automated search can be carried out against a known structural fragment (Chapter 2.3 ).
The `direct method' most familiar to the electron microscopist is the high-resolution electron micrograph of a crystalline lattice. Retrieval of an average structure from such a micrograph assumes that the experimental image conforms adequately to the `weak phase object' approximation, as discussed in Section 2.5.5. If this is so, the use of image-averaging techniques, e.g. Fourier filtration or correlational alignment, will allow the unit-cell contents to be visualized after the electron-microscope phase contrast transfer function is deconvoluted from the average image, also discussed in Section 2.5.5. Image analyses can also be extended to three dimensions, as discussed in Section 2.5.6, basically by employing tomographic reconstruction techniques to combine information from the several tilt projections taken from the crystalline object. The potential distribution of the unit cell to the resolution of the imaging experiment can then be used, via the Fourier transform, to obtain crystallographic phases for the electron-diffraction amplitudes recorded at the same resolution. This method for phase determination has been the mainstay of protein electron crystallography.
Once a set of phases is obtained from the Fourier transform of the deconvoluted image, they must, however, be referred to an allowed crystallographic origin. For many crystallographic space groups, this choice of origin may coincide with the location of a major symmetry element in the unit cell [see IT A (2005)]. Hence, since the Fourier transform of translation is a phase term, if an image shift is required to translate the origin of the repeating mass unit from the arbitrary position in the image to a specific site allowed by the space group, where the operation `' denotes convolution. The Fourier transform of this shifted density function will be In addition to the crystallographic phases , it will, therefore, be necessary to find the additional phase-shift term that will access an allowed unit-cell origin. Such origin searches are carried out automatically by some commercial image-averaging computer-software packages.
In addition to applications to thin protein crystals (e.g. Henderson et al., 1990; Jap et al., 1991; Kühlbrandt et al., 1994), there are numerous examples of molecular crystals that have been imaged to a resolution of 3–4 Å, many of which have been discussed by Fryer (1993). For π-delocalized compounds, which are the most stable in the electron beam against radiation damage, the best results (2 Å resolution) have been obtained at 500 kV from copper perchlorophthalocyanine epitaxically crystallized onto KCl. As shown by Uyeda et al. (1978–1979), the averaged images clearly depict the positions of the heavy Cu and Cl atoms, while the positions of the light atoms in the organic residue are not resolved. (The utility of image-derived phases as a basis set for phase extension will be discussed below.) A number of aromatic polymer crystals have also been imaged to about 3 Å resolution, as reviewed recently (Tsuji, 1989; Dorset, 1994b).
Aliphatic molecular crystals are much more difficult to study because of their increased radiation sensitivity. Nevertheless, monolamellar crystals of the paraffin n-C44H90 have been imaged to 2.5 Å resolution with a liquid-helium cryomicroscope (Zemlin et al., 1985). Similar images have been obtained at room temperature from polyethylene (Revol & Manley, 1986) and also a number of other aliphatic polymer crystals (Revol, 1991).
As noted by J. M. Cowley in Section 2.5.1, dynamical scattering can pose a significant barrier to the direct interpretation of high-resolution images from many inorganic materials. Nevertheless, with adequate control of experimental conditions (limiting crystal thickness, use of high-voltage electrons) some progress has been made. Pan & Crozier (1993) have described 2.0 Å images from zeolites in terms of the phase-grating approximation. A three-dimensional structural study has been carried out on an aluminosilicate by Wenk et al. (1992) with thin samples that conform to the weak-phase-object approximation at the 800 kV used for the imaging experiment. Heavy and light (e.g. oxygen) atoms were located in the micrographs in good agreement with an X-ray crystal structure. Heavy-atom positions from electron microscopic and X-ray structure analyses have also been favourably compared for two heavy-metal oxides (Hovmöller et al., 1984; Li & Hovmöller, 1988).
Conventional direct phasing techniques, as commonly employed in X-ray crystallography (e.g. see Chapter 2.2 ), have also been used for ab initio electron-crystallographic analyses. As in X-ray crystallography, probabilistic estimates of a linear combination of phases (Hauptman & Karle, 1953; Hauptman, 1972) are made after normalized structure factors are calculated via electron form factors, i.e. (Here, an overall temperature factor can be found from a Wilson plot. Because of multiple scattering, the value of B may be found occasionally to lie close to 0.0 Å2.) The phase invariant sums can be particularly effective for structure analysis. Of particular importance historically have been the -triple invariants where and . The probability of predicting is directly related to the value of where and Z is the value of the scattering factor at . Thus, the values of the phases are related to the measured structure factors, just as they are found to be in X-ray crystallography. The normalization described above imposes the point-atom structure (compensating for the fall-off of an approximately Gaussian form factor) often assumed in deriving the joint probability distributions. Especially for van der Waals structures, the constraint of positivity also holds in electron crystallography. (It is also quite useful for charged atoms so long as the reflections are not measured at very low angles.) Other useful phase invariant sums are the triples, where , and the quartets, where and . The prediction of a correct phase for an invariant is related in each case to the normalized structure-factor magnitudes.
The procedure for phase determination, therefore, is identical to the one used in X-ray crystallography (see Chapter 2.2 ). Using vectorial combinations of Miller indices, one generates triple and quartet invariants from available measured data and ranks them according to parameters such as A, defined above, which, as shown in Chapter 2.2 , are arguments of the Cochran formula. The invariants are thus listed in order of their reliability. This, in fact, generates a set of simultaneous equations in crystallographic phase. In order to begin solving these equations, it is permissible to define arbitrarily the phase values of a limited number of reflections (three for a three-dimensional primitive unit cell) for reflections with Miller-index parity and , where g is an even number. This defines the origin of a unit cell. For noncentrosymmetric unit cells, the condition for defining the origin, which depends on the space group, is somewhat more complicated and an enantiomorph-defining reflection must be added.
In the evaluation of phase-invariant sums above a certain probability threshold, phase values are determined algebraically after origin (and enantiomorph) definition until a large enough set is obtained to permit calculation of an interpretable potential map (i.e. where atomic positions can be seen). There may be a few invariant phase sums above this threshold probability value which are incorrectly predicted, leading either to false phase assignments or at least to phase assignments inconsistent with those found from other invariants. A small number of such errors can generally be tolerated. Another problem arises when an insufficient quantity of new phase values is assigned directly from the phase invariants after the origin-defining phases are defined. This difficulty may occur for small data sets, for example. If this is the case, it is possible that a new reflection of proper index parity can be used to define the origin. Alternatively, algebraic unknowns can be used to establish the phase linkage among certain reflections. If the structure is centrosymmetric, and when enough reflections are given at least symbolic phase assignments, maps are calculated and the correct structure is identified by inspection of the potential maps. When all goes well in this so-called `symbolic addition' procedure, the symbols are uniquely determined and there is no need to calculate more than a single map. If algebraic values are retained for certain phases because of limited vectorial connections in the data set, then a few maps may need to be generated so that the correct structure can be identified using the chemical knowledge of the investigator. The atomic positions identified can then be used to calculate phases for all observed data (via the structure-factor calculation) and the structure can be refined by Fourier (or, sometimes, least-squares) techniques to minimize the crystallographic R factor.
The first actual application of direct phasing techniques to experimental electron-diffraction data, based on symbolic addition procedures, was to two methylene subcell structures (an n-paraffin and a phospholipid; Dorset & Hauptman, 1976). Since then, evaluation of phase invariants has led to numerous other structures. For example, early texture electron-diffraction data sets obtained in Moscow (Vainshtein, 1964) were shown to be suitable for direct analysis. The structure of diketopiperazine (Dorset, 1991a) was determined from these electron-diffraction data (Vainshtein, 1955) when directly determined phases allowed computation of potential maps such as the one shown in Fig. 2.5.7.1. Bond distances and angles are in good agreement with the X-ray structure, particularly after least-squares refinement (Dorset & McCourt, 1994a). In addition, the structures of urea (Dorset, 1991b), using data published by Lobachev & Vainshtein (1961), paraelectric thiourea (Dorset, 1991b), using data published by Dvoryankin & Vainshtein (1960), and three mineral structures (Dorset, 1992a), from data published by Zvyagin (1967), have been determined, all using the original texture (or mosaic single-crystal) diffraction data. The most recent determination based on such texture diffraction data is that of basic copper chloride (Voronova & Vainshtein, 1958; Dorset, 1994c).
Potential map for diketopiperazine ([001] projection) after a direct phase determination with texture electron-diffraction intensity data obtained originally by Vainshtein (1955). |
Symbolic addition has also been used to assign phases to selected-area diffraction data. The crystal structure of boric acid (Cowley, 1953) has been redetermined, adding an independent low-temperature analysis (Dorset, 1992b). Additionally, a direct structure analysis has been reported for graphite, based on high-voltage intensity data (Ogawa et al., 1994). Two-dimensional data from several polymer structures have also been analysed successfully (Dorset, 1992c) as have three-dimensional intensity data (Dorset, 1991c,d; Dorset & McCourt, 1993).
Phase information from electron micrographs has also been used to aid phase determination by symbolic addition. Examples include the epitaxically oriented paraffins n-hexatriacontane (Dorset & Zemlin, 1990), n-tritriacontane (Dorset & Zhang, 1991) and a 1:1 solid solution of n-C32H66/n-C36H74 (Dorset, 1990a). Similarly, lamellar electron-diffraction data to ca 3 Å resolution from epitaxically oriented phospholipids have been phased by analysis of and -triplet invariants (Dorset, 1990b, 1991e,f), in one case combined with values from a 6 Å resolution electron microscope image (Dorset et al., 1990, 1993). Most recently, such data have been used to determine the layer packing of a phospholipid binary solid solution (Dorset, 1994d).
An ab initio direct phase analysis was carried out with zonal electron-diffraction data from copper perchlorophthalocyanine. Using intensities from a ca 100 Å thick sample collected at 1.2 MeV, the best map from a phase set with symbolic unknowns retrieves the positions of all the heavy atoms, equivalent to the results of the best images (Uyeda et al., 1978–1979). Using these positions to calculate an initial phase set, the positions of the remaining light C, N atoms were found by Fourier refinement so that the final bond distances and angles were in good agreement with those from X-ray structures of similar compounds (Dorset et al., 1991). A similar analysis has been carried out for the perbromo analogue (Dorset et al., 1992). Although dynamical scattering and secondary scattering significantly perturb the observed intensity data, the total molecular structure can be visualized after a Fourier refinement. Most recently, a three-dimensional structure determination was reported for C60 buckminsterfullerene based on symbolic addition with results most in accord with a rotationally disordered molecular packing (Dorset & McCourt, 1994b).
Given a triple phase relationship where h, k and form a vector sum, it is often possible to find a more reliable estimate of when all the possible vectorial contributions to it within the observed data set are considered as an average, viz: For actual phase determination, this can be formalized as follows. After calculating normalized structure-factor magnitudes from the observed to generate all possible phase triples within a reasonably high threshold, new phase values can be estimated after origin definition by use of the tangent formula (Karle & Hauptman, 1956): The reliability of the phase estimate depends on the variance , which is directly related to the magnitude of , i.e. is identical to the A value defined in the previous section. In the initial stages of phase determination is replaced by an expectation value until enough phases are available to permit its calculation.
The phase solutions indicated by the tangent formula can thus be ranked according to the phase variance and the determination of phases can be made symbolically from the most probable triple-product relationships. This procedure is equivalent to the one described above for the evaluation of phase-invariant sums by symbolic addition. This procedure may allow determination of a large enough basis phase set to produce an interpretable map.
An alternative procedure is to use an automatic version of the tangent formula in a multisolution process. This procedure is described in Chapter 2.2 . After origin definition, enough algebraic unknowns are defined (two values if centrosymmetric and four values, cycling through phase quadrants, if noncentrosymmetric) to access as many of the unknown phases as possible. These are used to generate a number of trial phase sets and the likelihood of identifying the correct solution is based on the use of some figure of merit.
Multisolution approaches employing the tangent formula include MULTAN (Germain et al., 1971), QTAN (Langs & DeTitta, 1975) and RANTAN (Yao, 1981). RANTAN is a version of MULTAN that allows for a larger initial random phase set (with suitable control of weights in the tangent formula). QTAN utilizes the definition, where for evaluating the phase variance. (Here , are modified Bessel functions.) After multiple solutions are generated, it is desirable to locate the structurally most relevant phase sets by some figure of merit. There are many that have been suggested (Chapter 2.2 ). The most useful figure of merit in QTAN has been the NQEST (De Titta et al., 1975) estimate of negative quartet invariants (see Chapter 2.2 ). More recently, this has been superseded by the minimal function (Hauptman, 1993): where and .
In the first application (Dorset et al., 1979) of multisolution phasing to electron-diffraction data (using the program QTAN), n-beam dynamical structure factors generated for cytosine and disodium 4-oxypyrimidine-2-sulfinate were used to assess the effect of increasing crystal thickness and electron accelerating voltage on the success of the structure determination. At 100 kV samples at least 80 Å thickness were usable for data collection and at 1000 kV this sample thickness limit could be pushed to 300 Å – or, perhaps, 610 Å if a partial structure were accepted for later Fourier refinement. NQEST identified the best phase solutions. Later QTAN was used to evaluate the effect of elastic crystal bend on the structure analysis of cytosine (Moss & Dorset, 1982).
In actual experimental applications, two forms of thiourea were investigated with QTAN (Dorset, 1992d), using published three-dimensional electron-diffraction intensities (Dvoryankin & Vainshtein, 1960, 1962). Analysis of the centrosymmetric paraelectric structure yielded results equivalent to those found earlier by symbolic addition (Dorset, 1991b). Analysis of the noncentrosymmetric ferroelectric form was also quite successful (Dorset, 1992d). In both cases, the correct structure was found at the lowest value of NQEST. Re-analysis of the diketopiperazine structure with QTAN also found the correct structure (Dorset & McCourt, 1994a) within the four lowest values of NQEST, but not the one at the lowest value. The effectiveness of this figure of merit became more questionable when QTAN was used to solve the noncentrosymmetric crystal structure of a polymer (Dorset, McCourt, Kopp et al., 1994). The solution could not be found readily when NQEST was used but was easily identified when the minimal function was employed instead.
MULTAN has been used to phase simulated data from copper perchlorophthalocyanine (Fan et al., 1985). Within the 2 Å resolution of the electron-microscope image, if one seeks phases for diffraction data in reciprocal-space regions where the objective lens phase contrast transfer function , the method proves to be successful. The method is also quite effective for phase extension from 2 Å to 1 Å diffraction resolution, where the low-angle data serve as a large initial phase set for the tangent formula. However, no useful results were found from an ab initio phase determination carried out solely with the electron-diffraction structure-factor magnitudes. Similar results were obtained when RANTAN was used to phase experimental data from this compound (Fan et al., 1991), i.e. the multisolution approach worked well for phase extension but not for ab initio phase determination. Additional tests were subsequently carried out with QTAN on an experimental hk0 electron-diffraction data set collected at 1200 kV (Dorset, McCourt, Fryer et al., 1994). Again, ab initio phase determination is not possible by this technique. However, if a basis set was constructed from the Fourier transform of a 2.4 Å image, a correct solution could be found, but not at the lowest value of NQEST. This figure of merit was useful, however, when the basis set was taken from the symbolic addition determination mentioned in the previous section.
Another method of phase determination, which is best suited to refining or extending a partial phase set, is the Hoppe–Gassmann density modification procedure (Hoppe & Gassmann, 1968; Gassmann & Zechmeister, 1972; Gassmann, 1976). The procedure is very simple but also very computer-intensive. Starting with a small set of (phased) , an initial potential map is calculated by Fourier transformation. This map is then modified by some real-space function, which restricts peak sizes to a maximum value and removes all negative density regions. The modified map is then Fourier-transformed to produce a set of phased structure factors. Phase values are accepted via another modification function in reciprocal space, e.g. , where p is a threshold quantity. The new set is then transformed to obtain a new and the phase refinement continues iteratively until the phase solution converges (judged by lower crystallographic R values).
The application of density modification procedures to electron-crystallographic problems was assessed by Ishizuka et al. (1982), who used simulated data from copper perchlorophthalocyanine within the resolution of the electron-microscope image. The method was useful for finding phase values in reciprocal-space regions where the transfer function . As a technique for phase extension, density modification was acceptable for test cases where the resolution was extended from 1.67 to 1.0 Å, or 2.01 to 1.21 Å, but it was not very satisfactory for a resolution enhancement from 2.5 to 1.67 Å. There appear to have been no tests of this method yet with experimental data. However, the philosophy of this technique will be met again below in the description of the the maximum entropy and likelihood procedure.
One of the first relationships ever derived for phase determination is the Sayre (1952) equation: which is a simple convolution of phased structure factors multiplied by a function of the atomic scattering factors. For structures with non-overlapping atoms, consisting of one atomic species, it is an exact expression. Although the convolution term resembles part of the tangent formula above, no statistical averaging is implied (Sayre, 1980). In X-ray crystallography this relationship has not been used very often, despite its accuracy. Part of the reason for this is that it requires relatively high resolution data for it to be useful. It can also fail for structures comprised of different atomic species.
Since, relative to X-ray scattering factors, electron scattering factors span a narrower range of magnitudes at , it might be thought that the Sayre equation would be particularly useful in electron crystallography. In fact, Liu et al. (1988) were able to extend phases for simulated data from copper perchlorophthalocyanine starting at the image resolution of 2 Å and reaching the 1 Å limit of an electron-diffraction data set. This analysis has been improved with a 2.4 Å basis set obtained from the Fourier transform of an electron micrograph of this material at 500 kV and extended to the 1.0 Å limit of a 1200 kV electron-diffraction pattern (Dorset et al., 1995). Using the partial phase sets for zonal diffraction data from several polymers by symbolic addition (see above), the Sayre equation has been useful for extending into the whole hk0 set, often with great accuracy. The size of the basis set is critical but the connectivity to access all reflections is more so. Fan and co-workers have had considerable success with the analysis of incommensurately modulated structures. The average structure (basis set) is found by high-resolution electron microscopy and the `superlattice' reflections, corresponding to the incommensurate modulation, are assigned phases in hyperspace by the Sayre convolution. Examples include a high superconductor (Mo et al., 1992) and the mineral ankangite (Xiang et al., 1990). Phases of regular inorganic crystals have also been extended from the electron micrograph to the electron-diffraction resolution by this technique (Hu et al., 1992).
In an investigation of how direct methods might be used for phase extension in protein electron crystallography, low-resolution phases from two proteins, bacteriorhodopsin (Henderson et al., 1986) and halorhodopsin (Havelka et al., 1993) were extended to higher resolution with the Sayre equation (Dorset et al., 1995). For the noncentrosymmetric bacteriorhodopsin hk0 projection a 10 Å basis set was used, whereas a 15 Å set was accepted for the centrosymmetric halorhodopsin projection. In both cases, extensions to 6 Å resolution were reasonably successful. For bacteriorhodopsin, for which data were available to 3.5 Å, problems with the extension were encountered near 5 Å, corresponding to a minimum in a plot of average intensity versus resolution. Suggestions were made on how a multisolution procedure might be successful beyond this point.
Maximum entropy has been applied to electron crystallography in several ways. In the sense that images are optimized, the entropy term where and is a pixel density, has been evaluated for various test electron-microscope images. For crystals, the true projected potential distribution function is thought to have the maximum value of S. If the phase contrast transfer function used to obtain a micrograph is unknown, test images (i.e. trial potential maps) can be calculated for different values of . The value that corresponds to the maximum entropy would be near the true defocus. In this way, the actual objective lens transfer function can be found for a single image (Li, 1991) in addition to the other techniques suggested by this group.
Another use of the maximum-entropy concept is to guide the progress of a direct phase determination (Bricogne & Gilmore, 1990; Gilmore et al., 1990). Suppose that there is a small set H of known phases (corresponding either to origin definition, or the Fourier transform of an electron micrograph, or both) with associated unitary structure-factor amplitudes . [The unitary structure factor is defined as .] As usual, the task is to expand into the unknown phase set K to solve the crystal structure. From Bayes' theorem, the procedure is based on an operation where . This means that the probability of successfully deriving a potential map, given diffraction data, is estimated. This so-called posterior probability is approximately proportional to the product of the probability of generating the map (known as the prior) and the probability of generating the data, given the map (known as the likelihood). The latter probability consults the observed data and can be used as a figure of merit.
Beginning with the basis set H, a trial map is generated from the limited number of phased structure factors. As discussed above, the map can be immediately improved by removing all negative density. The map can be improved further if its entropy is maximized using the equation given above for S. This produces the so-called maximum-entropy prior .
So far, it has been assumed that all . If large reflections from the K set are now added and their phase values are permuted, then a number of new maps can be generated and their entropies can be maximized as before. This creates a phasing `tree' with many possible solutions; individual branch points can have further reflections added via permutations to produce further sub-branches, and so on. Obviously, some figure of merit is needed to `prune' the tree, i.e. to find likely paths to a solution.
The desired figure of merit is the likelihood . First a quantity where (the calculated unitary structure factors) and (the observed unitary structure factors), is defined. From this one can calculate The null hypothesis can also be calculated from the above when , so that the likelihood gain ranks the nodes of the phasing tree in order of the best solutions.
Applications have been made to experimental electron-crystallographic data. A small-molecule structure starting with phases from an electron micrograph and extending to electron-diffraction resolution has been reported (Dong et al., 1992). Other experimental electron-diffraction data sets used in other direct phasing approaches (see above) also have been assigned phases by this technique (Gilmore, Shankland & Bricogne, 1993). These include intensities from diketopiperazine and basic copper chloride. An application of this procedure in protein structure analysis has been published by Gilmore et al. (1992) and Gilmore, Shankland & Fryer (1993). Starting with 15 Å phases, it was possible to extend phases for bacteriorhodopsin to the limits of the electron-diffraction pattern, apparently with greater accuracy than possible with the Sayre equation (see above).
The aim of electron-crystallographic data collection is to minimize the effect of dynamical scattering, so that the unit-cell potential distribution or its Fourier transform is represented significantly in the recorded signal. It would be a mistake, however, to presume that these data ever conform strictly to the kinematical approximation, for there is always some deviation from this ideal scattering condition that can affect the structure analysis. Despite this fact, some direct phasing procedures have been particularly `robust', even when multiple scattering perturbations to the data are quite obvious (e.g. as evidenced by large crystallographic residuals).
The most effective direct phasing procedures seem to be those based on the triple invariants. These phase relationships will not only include the symbolic addition procedure, as it is normally carried out, but also the tangent formula and the Sayre equation (since it is well known that this convolution can be used to derive the functional form of the three-phase invariant). The strict ordering of magnitudes is, therefore, not critically important so long as there are no major changes from large to small values (or vice versa). This was demonstrated in direct phase determinations of simulated n-beam dynamical diffraction data from a sulfur-containing polymer (Dorset & McCourt, 1992). Nevertheless, there is a point where measured data cannot be used. For example, intensities from ca 100 Å-thick epitaxically oriented copper perchlorophthalocyanine crystals become less and less representative of the unit-cell transform at lower electron-beam energies (Tivol et al., 1993) and, accordingly, the success of the phase determination is compromised (Dorset, McCourt, Fryer et al., 1994). The similarity between the Sayre convolution and the interactions of structure-factor terms in, e.g., the multislice formulation of n-beam dynamical scattering was noted by Moodie (1965). It is interesting to note that dynamical scattering interactions observed by direct excitation of and triples in convergent-beam diffraction experiments can actually be exploited to determine crystallographic phases to very high precision (Spence & Zuo, 1992, pp. 56–63).
While the evaluation of positive quartet invariant sums (see Chapter 2.2 ) seems to be almost as favourable in the electron diffraction case as is the evaluation of triples, negative quartet invariants seem to be particularly sensitive to dynamical diffraction. If dynamical scattering can be modelled crudely by a convolutional smearing of the diffraction intensities, then the lowest structure-factor amplitudes, and hence the estimates of lowest values, will be the ones most compromised. Since the negative-quartet relationships require an accurate prediction of small `cross-term' values, multiple scattering can, therefore, limit the efficacy of this invariant for phase determination. In initial work, negative quartets have been mostly employed in the NQEST figure of merit, and analyses (Dorset, McCourt, Fryer et al., 1994; Dorset & McCourt, 1994a) have shown how the degradation of weak kinematical terms effectively reduced its effectiveness for locating correct structure solutions via the tangent formula, even though the tangent formula itself (based on triple phase estimates) was quite effective for phase determination. Substitution of the minimal function for NQEST seems to have overcome this difficulty. [It should be pointed out, though, that only the -triple contribution to is considered.]
Structure refinement is another area where the effects of dynamical scattering are also problematic. For example, in the analysis of the paraelectric thiourea structure (Dorset, 1991b) from published texture diffraction data (Dvoryankin & Vainshtein, 1960), it was virtually impossible to find a chemically reasonable structure geometry by Fourier refinement, even though the direct phase determination itself was quite successful. The best structure was found only when higher-angle intensities (i.e. those least affected by dynamical scattering) were used to generate the potential map. Later analyses on heavy-atom-containing organics (Dorset et al., 1992) found that the lowest kinematical R-factor value did not correspond to the chemically correct structure geometry. This observation was also made in the least-squares refinement of diketopiperazine (Dorset & McCourt, 1994a). It is obvious that, if a global minimum is sought for the crystallographic residual, then dynamical structure factors, rather than kinematical values, should be compared to the observed values (Dorset et al., 1992). Ways of integrating such calculations into the refinement process have been suggested (Sha et al., 1993). Otherwise one must constrain the refinement to chemically reasonable bonding geometry in a search for a local R-factor minimum.
Corrections for such deviations from the kinematical approximation are complicated by the presence of other possible data perturbations, especially if microareas are being sampled, e.g. in typical selected-area diffraction experiments. Significant complications can arise from the diffraction incoherence observed from elastically deformed crystals (Cowley, 1961) as well as secondary scattering (Cowley et al., 1951). These complications were also considered for the larger (e.g. millimeter diameter) areas sampled in an electron-diffraction camera when recording texture diffraction patterns (Turner & Cowley, 1969), but, because of the crystallite distributions, it is sometimes found that the two-beam dynamical approximation is useful (accounting for a number of successful structure analyses carried out in Moscow).
References
Avilov, A. S. (1979). Electrical measurement of reflection intensities on electron diffraction from mosaic single crystals. Sov. Phys. Crystallogr. 24, 103–104.Google ScholarAvilov, 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.Google Scholar
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.Google Scholar
Bando, Y. (1981). Weak asymmetry in β-Si3N4 as revealed by convergent beam electron diffraction. Acta Cryst. B39, 185–189.Google Scholar
Bethe, H. A. (1928). Theorie der Beugung von Elektronen an Kristallen. Ann. Phys. (Leipzig), 87, 55–129.Google Scholar
Bilhorn, D. E., Foldy, L. L., Thaler, R. M. & Tobacman, W. (1964). Remarks concerning reciprocity in quantum mechanics. J. Math. Phys. 5, 435–441.Google Scholar
Blackman, M. (1939). On the intensities of electron diffraction rings. Proc. R. Soc. London Ser. A, 173, 68–82.Google Scholar
Bracewell, B. N. (1956). Strip integration in radio astronomy. Austr. J. Phys. 9, 198–217.Google Scholar
Bricogne, G. & Gilmore, C. J. (1990). A multisolution method of phase determination by combined maximization of entropy and likelihood. I. Theory, algorithms and strategy. Acta Cryst. A46, 284–297.Google Scholar
Brisse, F. (1989). Electron diffraction of synthetic polymers: the model compound approach to polymer structure. J. Electron Microsc. Tech. 11, 272–279.Google Scholar
Buxton, B., Eades, J. A., Steeds, J. W. & Rackham, G. M. (1976). The symmetry of electron diffraction zone axis patterns. Philos. Trans. R. Soc. London Ser. A, 181, 171–193.Google Scholar
Carpenter, R. W. & Spence, J. C. H. (1982). Three-dimensional strain-field information in convergent-beam electron diffraction patterns. Acta Cryst. A38, 55–61.Google Scholar
Chou, C. T., Anderson, S. C., Cockayne, D. J. H., Sikorski, A. Z. & Vaughan, M. R. (1994). Ultramicroscopy, 55, 334–347.Google Scholar
Cochran, W., Crick, F. H. C. & Vand, V. (1952). The structure of synthetic polypeptides. 1. The transform of atoms on a helix. Acta Cryst. 5, 581–586.Google Scholar
Cowley, J. M. (1953). Structure analysis of single crystals by electron diffraction. II. Disordered boric acid structure. Acta Cryst. 6, 522–529.Google Scholar
Cowley, J. M. (1956). A modified Patterson function. Acta Cryst. 9, 397–398.Google Scholar
Cowley, J. M. (1961). Diffraction intensities from bent crystals. Acta Cryst. 14, 920–927.Google Scholar
Cowley, J. M. (1969). Image contrast in transmission scanning electron microscopy. Appl. Phys. Lett. 15, 58–59.Google Scholar
Cowley, J. M. (1981). Diffraction physics, 2nd ed. Amsterdam: North-Holland.Google Scholar
Cowley, J. M. (1995). Diffraction physics, 3rd ed. Amsterdam: North-Holland.Google Scholar
Cowley, J. M. & Au, A. Y. (1978). Image signals and detector configurations for STEM. In Scanning electron microscopy, Vol. 1, pp. 53–60. AMF O'Hare, Illinois: SEM Inc.Google Scholar
Cowley, J. M. & Moodie, A. F. (1957). The scattering of electrons by atoms and crystals. I. A new theoretical approach. Acta Cryst. 10, 609–619.Google Scholar
Cowley, J. M. & Moodie, A. F. (1959). The scattering of electrons by atoms and crystals. III. Single-crystal diffraction patterns. Acta Cryst. 12, 360–367.Google Scholar
Cowley, J. M. & Moodie, A. F. (1960). Fourier images. IV. The phase grating. Proc. Phys. Soc. London, 76, 378–384.Google Scholar
Cowley, J. M., Moodie, A. F., Miyake, S., Takagi, S. & Fujimoto, F. (1961). The extinction rules for reflections in symmetrical electron diffraction spot patterns. Acta Cryst. 14, 87–88.Google Scholar
Cowley, J. M., Rees, A. L. G. & Spink, J. A. (1951). Secondary elastic scattering in electron diffraction. Proc. Phys. Soc. London Sect. A, 64, 609–619.Google Scholar
Cramér, H. (1954). Mathematical methods of statistics. University of Princeton.Google Scholar
Creek, R. C. & Spargo, A. E. C. (1985). Electron optical study of rutile. J. Appl. Cryst. 18, 197–204.Google Scholar
Crowther, R. A. & Amos, L. A. (1971). Harmonic analysis of electron microscope images with rotational symmetry. J. Mol. Biol. 60, 123–130.Google Scholar
Crowther, R. A., Amos, L. A., Finch, J. T., DeRosier, D. J. & Klug, A. (1970). Three dimensional reconstruction of spherical viruses by Fourier synthesis from electron micrographs. Nature (London), 226, 421–425.Google Scholar
Crowther, R. A., DeRosier, D. J. & Klug, A. (1970). The reconstruction of a three-dimensional structure from projections and its application to electron microscopy. Proc. R. Soc. London Ser. A, 317, 319–340.Google Scholar
Crowther, R. A. & Klug, A. (1974). Three dimensional image reconstruction on an extended field – a fast, stable algorithm. Nature (London), 251, 490–492.Google Scholar
Dawson, B., Goodman, P., Johnson, A. W. S., Lynch, D. F. & Moodie, A. F. (1974). Some definitions and units in electron diffraction. Acta Cryst. A30, 297–298.Google Scholar
Deans, S. R. (1983). The Radon transform and some of its applications. New York: John Wiley.Google Scholar
DeRosier, D. J. & Klug, A. (1968). Reconstruction of three dimensional structures from electron micrographs. Nature (London), 217, 130–134.Google Scholar
DeRosier, D. J. & Moore, P. B. (1970). Reconstruction of three-dimensional images from electron micrographs of structure with helical symmetry. J. Mol. Biol. 52, 355–369.Google Scholar
De Titta, G. T., Edmonds, J. W., Langs, D. A. & Hauptman, H. (1975). Use of negative quartet cosine invariants as a phasing figure of merit: NQEST. Acta Cryst. A31, 472–479.Google Scholar
Dong, W., Baird, T., Fryer, J. R., Gilmore, C. J., MacNicol, D. D., Bricogne, G., Smith, D. J., O'Keefe, M. A. & Hovmöller, S. (1992). Electron microscopy at 1 Å resolution by entropy maximization and likelihood ranking. Nature (London), 355, 605–609.Google Scholar
Dorset, D. L. (1976). The interpretation of quasi-kinematical single-crystal electron diffraction intensity data from paraffins. Acta Cryst. A32, 207–215.Google Scholar
Dorset, D. L. (1987). Electron diffraction structure analysis of phospholipids. J. Electron Microsc. Tech. 7, 35–46.Google Scholar
Dorset, D. L. (1990a). Direct structure analysis of a paraffin solid solution. Proc. Natl Acad. Sci. USA, 87, 8541–8544.Google Scholar
Dorset, D. L. (1990b). Direct determination of crystallographic phases for diffraction data from phospholipid multilamellar arrays. Biophys. J. 58, 1077–1087.Google Scholar
Dorset, D. L. (1991a). Electron diffraction structure analysis of diketopiperazine – a direct phase determination. Acta Cryst. A47, 510–515.Google Scholar
Dorset, D. L. (1991b). Is electron crystallography possible? The direct determination of organic crystal structures. Ultramicroscopy, 38, 23–40.Google Scholar
Dorset, D. L. (1991c). Electron diffraction structure analysis of polyethylene. A direct phase determination. Macromolecules, 24, 1175–1178.Google Scholar
Dorset, D. L. (1991d). Electron crystallography of linear polymers: direct structure analysis of poly(ɛ-caprolactone). Proc. Natl Acad. Sci. USA, 88, 5499–5502.Google Scholar
Dorset, D. L. (1991e). Direct determination of crystallographic phases for diffraction data from lipid bilayers. I. Reliability and phase refinement. Biophys. J. 60, 1356–1365.Google Scholar
Dorset, D. L. (1991f). Direct determination of crystallographic phases for diffraction data from lipid bilayers. II. Refinement of phospholipid structures. Biophys. J. 60, 1366–1373.Google Scholar
Dorset, D. L. (1992a). Direct phasing in electron crystallography: determination of layer silicate structures. Ultramicroscopy, 45, 5–14.Google Scholar
Dorset, D. L. (1992b). Direct methods in electron crystallography – structure analysis of boric acid. Acta Cryst. A48, 568–574.Google Scholar
Dorset, D. L. (1992c). Electron crystallography of linear polymers: direct phase determination for zonal data sets. Macromolecules, 25, 4425–4430.Google Scholar
Dorset, D. L. (1992d). Automated phase determination in electron crystallography: thermotropic phases of thiourea. Ultramicroscopy, 45, 357–364.Google Scholar
Dorset, D. L. (1994a). Electron crystallography of organic molecules. Adv. Electron. Electron Phys. 88, 111–197.Google Scholar
Dorset, D. L. (1994b). Electron crystallography of linear polymers. In Characterization of solid polymers. New techniques and developments, edited by S. J. Spells, pp.1–16. London: Chapman and Hall.Google Scholar
Dorset, D. L. (1994c). Electron crystallography of inorganic compounds. Direct determination of the basic copper chloride structure CuCl2·3Cu(OH)2. J. Chem. Crystallogr. 24, 219–224.Google Scholar
Dorset, D. L. (1994d). Direct determination of layer packing for a phospholipid solid solution at 0.32 nm resolution. Proc. Natl Acad. Sci. USA, 91, 4920–4924.Google Scholar
Dorset, D. L., Beckmann, E. & Zemlin, F. (1990). Direct determination of a phospholipid lamellar structure at 0.34 nm resolution. Proc. Natl Acad. Sci. USA, 87, 7570–7573.Google Scholar
Dorset, D. L. & Hauptman, H. A. (1976). Direct phase determination for quasi-kinematical electron diffraction intensity data from organic microcrystals. Ultramicroscopy, 1, 195–201.Google Scholar
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.Google Scholar
Dorset, D. L., Kopp, S., Fryer, J. R. & Tivol, W. F. (1995). The Sayre equation in electron crystallography. Ultramicroscopy, 57, 59–89.Google Scholar
Dorset, D. L. & McCourt, M. P. (1992). Effect of dynamical scattering on successful direct phase determination in electron crystallography – a model study. Trans. Am. Crystallogr. Assoc. 28, 105–113.Google Scholar
Dorset, D. L. & McCourt, M. P. (1993). Electron crystallographic analysis of a polysaccharide structure – direct phase determination and model refinement for mannan I. J. Struct. Biol. 111, 118–124.Google Scholar
Dorset, D. L. & McCourt, M. P. (1994a). Automated structure analysis in electron crystallography: phase determination with the tangent formula and least-squares refinement. Acta Cryst. A50, 287–292.Google Scholar
Dorset, D. L. & McCourt, M. P. (1994b). Disorder and molecular packing of C60 buckminsterfullerene: a direct electron-crystallographic analysis. Acta Cryst. A50, 344–351.Google Scholar
Dorset, D. L., McCourt, M. P., Fryer, J. R., Tivol, W. F. & Turner, J. N. (1994). The tangent formula in electron crystallography: phase determination of copper perchlorophthalocyanine. Microsc. Soc. Am. Bull. 24, 398–404.Google Scholar
Dorset, D. L., McCourt, M. P., Kopp, S., Wittmann, J.-C. & Lotz, B. (1994). Direct determination of polymer crystal structures by electron crystallography – isotactic poly(1-butene), form III. Acta Cryst. B50, 201–208.Google Scholar
Dorset, D. L., McCourt, M. P., Tivol, W. F. & Turner, J. N. (1993). Electron diffraction from phospholipids – an approximate correction for dynamical scattering and tests for a correct phase determination. J. Appl. Cryst. 26, 778–786.Google Scholar
Dorset, D. L., Tivol, W. F. & Turner, J. N. (1991). Electron crystallography at atomic resolution: ab initio structure analysis of copper perchlorophthalocyanine. Ultramicroscopy, 38, 41–45.Google Scholar
Dorset, D. L., Tivol, W. F. & Turner, J. N. (1992). Dynamical scattering and electron crystallography – ab initio structure analysis of copper perbromophthalocyanine. Acta Cryst. A48, 562–568.Google Scholar
Dorset, D. L. & Zemlin, F. (1990). Direct phase determination in electron crystallography: the crystal structure of an n-paraffin. Ultramicroscopy, 33, 227–236.Google Scholar
Dorset, D. L. & Zhang, W. P. (1991). Electron crystallography at atomic resolution: the structure of the odd-chain paraffin n-tritriacontane. J. Electron Microsc. Tech. 18, 142–147.Google Scholar
Dvoryankin, V. F. & Vainshtein, B. K. (1960). An electron diffraction study of thiourea. Sov. Phys. Crystallogr. 5, 564–574.Google Scholar
Dvoryankin, V. F. & Vainshtein, B. K. (1962). An electron diffraction study of the low-temperature ferroelectric form of thiourea. Sov. Phys. Crystallogr. 6, 765–772.Google Scholar
Eades, J. A. (1980). Another way to form zone axis patterns. Inst. Phys. Conf. Ser. 52, 9–12.Google Scholar
Eades, J. A., Shannon, M. D. & Buxton, B. F. (1983). Crystal symmetry from electron diffraction. In Scanning electron microscopy, 1983/III, pp. 1051–1060. Chicago: SEM Inc.Google Scholar
Erickson, H. P. & Klug, A. (1971). Measurements and compensation of defocusing and aberrations by Fourier processing of electron micrographs. Philos. Trans. R. Soc. London Ser. B, 261, 105–118.Google Scholar
Fan, H. F., Xiang, S. B., Li, F. H., Pan, Q., Uyeda, N. & Fujiyoshi, Y. (1991). Image resolution enhancement by combining information from electron diffraction pattern and micrograph. Ultramicroscopy, 36, 361–365.Google Scholar
Fan, H.-F., Zhong, Z.-Y., Zheng, C.-D. & Li, F.-H. (1985). Image processing in high-resolution electron microscopy using the direct method. I. Phase extension. Acta Cryst. A41, 163–165.Google Scholar
Frank, J. (1975). Averaging of low exposure electron micrographs of non-periodic objects. Ultramicroscopy, 1, 159–162.Google Scholar
Frank, J. (1980). The role of correlation techniques in computer image processing. In Computer processing of electron microscopy images, edited by P. W. Hawkes, pp. 187–222. Berlin: Springer-Verlag.Google Scholar
Fraser, H. L., Maher, D. M., Humphreys, C. J., Hetherington, C. J. D., Knoell, R. V. & Bean, J. C. (1985). The detection of local strains in strained superlattices. In Microscopy of semiconducting materials, pp. 1–5. London: Institute of Physics.Google Scholar
Fryer, J. R. (1993). Electron crystallography of small organic molecules. Microsc. Soc. Am. Bull. 23, 44–56.Google Scholar
Fujimoto, F. (1959). Dynamical theory of electron diffraction in Laue-case. I. General theory. J. Phys. Soc. Jpn, 14(11), 1158–1168.Google Scholar
Fujiwara, K. (1961). Relativistic dynamical theory of electron diffraction. J. Phys. Soc. Jpn, 16, 2226–2238.Google Scholar
Fukuhara, A. (1966). Many-ray approximations in the dynamical theory of electron diffraction. J. Phys. Soc. Jpn, 21, 2645–2662.Google Scholar
Gabor, D. (1949). Microscopy by reconstructed wavefronts. Proc. R. Soc. London Ser. A, 197, 454–487.Google Scholar
Gassmann, J. (1976). Improvement and extension of approximate phase sets in structure determination. In Crystallographic computing techniques, edited by F. R. Ahmed, pp. 144–154. Copenhagen: Munksgaard.Google Scholar
Gassmann, J. & Zechmeister, K. (1972). Limits of phase expansion in direct methods. Acta Cryst. A28, 270–280.Google Scholar
Germain, G., Main, P. & Woolfson, M. M. (1971). The application of phase relationships to complex structures. III. The optimum use of phase relationships. Acta Cryst. A27, 368–376.Google Scholar
Gilbert, P. F. C. (1972a). The reconstruction of a three-dimensional structure from projections and its application to electron microscopy. II. Direct methods. Proc. R. Soc. London Ser. B, 182, 89–102.Google Scholar
Gilbert, P. F. C. (1972b). Iterative methods for the three-dimensional reconstruction of an object from projections. J. Theor. Biol. 36, 105–117.Google Scholar
Gilmore, C. J., Bricogne, G. & Bannister, C. (1990). A multisolution method of phase determination by combined maximization of entropy and likelihood. II. Application to small molecules. Acta Cryst. A46, 297–308.Google Scholar
Gilmore, C. J., Shankland, K. & Bricogne, G. (1993). Applications of the maximum entropy method to powder diffraction and electron crystallography. Proc. R. Soc. London Ser. A, 442, 97–111.Google Scholar
Gilmore, C. J., Shankland, K. & Fryer, J. R. (1992). The application of the maximum entropy method to electron microscopy data for purple membrane. Trans. Am. Crystallogr. Assoc. 28, 129–139.Google Scholar
Gilmore, C. J., Shankland, K. & Fryer, J. R. (1993). Phase extension in electron crystallography using the maximum entropy method and its application to two-dimensional purple membrane data from Halobacterium halobium. Ultramicroscopy, 49, 132–146.Google Scholar
Gjønnes, J. & Høier, R. (1971). The application of non-systematic many-beam dynamic effects to structure-factor determination. Acta Cryst. A27, 313–316.Google Scholar
Gjønnes, J. & Moodie, A. F. (1965). Extinction conditions in dynamic theory of electron diffraction patterns. Acta Cryst. 19, 65–67.Google Scholar
Glauber, R. & Schomaker, V. (1953). The theory of electron diffraction. Phys. Rev. 89, 667–670.Google Scholar
Goncharov, A. B. (1987). Integral geometry and 3D-reconstruction of arbitrarily oriented identical particles from their electron micrographs. Sov. Phys. Crystallogr. 32, 663–666.Google Scholar
Goncharov, A. B., Vainshtein, B. K., Ryskin, A. I. & Vagin, A. A. (1987). Three-dimensional reconstruction of arbitrarily oriented identical particles from their electron photomicrographs. Sov. Phys. Crystallogr. 32, 504–509.Google Scholar
Goodman, P. (1974). The role of upper layer interactions in electron diffraction. Nature (London), 251, 698–701.Google Scholar
Goodman, P. (1984a). A matrix basis for CBED pattern analysis. Acta Cryst. A40, 522–526.Google Scholar
Goodman, P. (1984b). A retabulation of the 80 layer groups for electron diffraction usage. Acta Cryst. A40, 633–642.Google Scholar
Goodman, P., McLean, J. D., Wilson, I. J. & Olsen, A. (1984). Optical microdiffraction and image analysis of subsymmetries in Nb2O5 tunnel structures. In Analytical electron microscopy–1984, pp. 130–134. San Francisco Press.Google Scholar
Goodman, P. & Miller, P. (1993). Reassessment of the symmetry of the 221 PbBiSrCaCuO structure using LACBED and high-resolution SAD: the relevance of Cowley's theory of disorder scattering to a real-space structural analysis. Ultramicroscopy, 52, 549–556.Google Scholar
Goodman, P., Miller, P., White, T. J. & Withers, R. L. (1992). Symmetry determination and Pb-site ordering analysis for the n = 1,2 PbxBi2 − xSr2Can − 1CunO4 + 2n + δ compounds by convergent-beam and selected-area electron diffraction. Acta Cryst. B48, 376–387.Google Scholar
Goodman, P. & Whitfield, H. J. (1980). The space group determination of GaS and Cu3As2S3I by convergent beam electron diffraction. Acta Cryst. A36, 219–228.Google Scholar
Gordon, R. (1974). A tutorial on ART (algebraic reconstruction techniques). IEEE Trans. Nucl. Sci. NS-21, 78–93.Google Scholar
Gordon, R., Bender, R. & Herman, G. T. (1970). Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography. J. Theor. Biol. 29, 471–481.Google Scholar
Gordon, R. & Herman, G. T. (1971). Reconstruction of pictures from their projections. Commun. ACM, 14, 759–768.Google Scholar
Grzinic, G. (1985). Calculation of incommensurate diffraction intensities from disordered crystals. Philos. Mag. A, 52, 161–187.Google Scholar
Gunning, J. & Goodman, P. (1992). Reciprocity in electron diffraction. Acta Cryst. A48, 591–595.Google Scholar
Gurskaya, G. V., Lobanova, G. M. & Vainshtein, B. K. (1971). X-ray diffraction and electron-microscope study of hexagonal catalase crystal. Sov. Phys. Crystallogr. 16, 662–669.Google Scholar
Hashimoto, H., Endoh, H., Tanji, T., Ono, A. & Watanabe, E. (1977). Direct observation of fine structure within images of atoms in crystals by transmission electron microscopy. J. Phys. Soc. Jpn, 42, 1073–1074.Google Scholar
Hashimoto, H., Mannami, M. & Naiki, T. (1961). Dynamical theory of electron diffraction for the electron microscope image of crystal lattices. I. Image of single crystals. II. Image of superposed crystals (moiré pattern). Philos. Trans. R. Soc. London, 253, 459–516.Google Scholar
Hauptman, H. (1972). Crystal structure determination. The role of the cosine seminvariants. NY: Plenum Press.Google Scholar
Hauptman, H. (1993). A minimal principle in X-ray crystallography: starting in a small way. Proc. R. Soc. London Ser. A, 442, 3–12.Google Scholar
Hauptman, H. & Karle, J. (1953). Solution of the phase problem. I. The centrosymmetric crystal. American Crystallographic Association Monograph No. 3. Ann Arbor, MI: Edwards Brothers.Google Scholar
Havelka, W., Henderson, R., Heymann, J. A. W. & Oesterhelt, D. (1993). Projection structure of halorhodopsin from Halobacterium halobium at 6 Å resolution obtained by electron cryomicroscopy. J. Mol. Biol. 234, 837–846.Google Scholar
Henderson, R., Baldwin, J. M., Ceska, T. A., Zemlin, F., Beckmann, E. & Downing, K. H. (1990). Model for the structure of bacteriorhodopsin based on high-resolution electron cryomicroscopy. J. Mol. Biol. 213, 899–929.Google Scholar
Henderson, R., Baldwin, J. M., Downing, K. H., Lepault, J. & Zemlin, F. (1986). Structure of purple membrane from Halobacterium halobium: recording, measurement and evaluation of electron micrographs at 3.5 Å resolution. Ultramicroscopy, 19, 147–178.Google Scholar
Henderson, R. & Unwin, P. N. T. (1975). Three-dimensional model of purple membrane obtained by electron microscopy. Nature (London), 257, 28–32.Google Scholar
Herman, G. T. (1980). Image reconstruction from projection: the fundamentals of computerized tomography. New York: Academic Press.Google Scholar
Herrmann, K. H., Krahl, D. & Rust, H.-P. (1980). Low-dose image recording by TV techniques. In Electron microscopy at molecular dimensions, edited by W. Baumeister & W. Vogell, pp. 186–193. Berlin: Springer-Verlag.Google Scholar
Hirsch, P. B., Howie, A., Nicholson, R. B., Pashley, D. W. & Whelan, M. J. (1965). Electron microscopy of thin crystals. London: Butterworths.Google Scholar
Hoppe, W. (1971). Use of zone correction plate and other techniques for structure determination of aperiodic objects at atomic resolution using a conventional electron microscope. Philos. Trans. R. Soc. London Ser. B, 261, 71–94.Google Scholar
Hoppe, W., Bussler, P., Feltynowski, A., Hunsmann, N. & Hirt, A. (1973). Some experience with computerized image reconstruction methods. In Image processing and computer-aided design in electron optics, edited by R. W. Hawkes, pp. 92–126. London: Academic Press.Google Scholar
Hoppe, W. & Gassmann, J. (1968). Phase correction, a new method to solve partially known structures. Acta Cryst. B24, 97–107.Google Scholar
Hoppe, W. & Typke, D. (1979). Three-dimensional reconstruction of aperiodic objects in electron microscopy. In Advances in structure research by diffraction method. Oxford: Pergamon Press.Google Scholar
Horstmann, M. & Meyer, G. (1965). Messung der Elektronenbeugungsintensitaten polykristalliner Aluminium schichten bei tiefer Temperatur und Vergleich mit der dynamischen Theorie. Z. Phys. 182, 380–397.Google Scholar
Hovmöller, S., Sjögren, A., Farrants, G., Sundberg, M. & Marinder, B. O. (1984). Accurate atomic positions from electron microscopy. Nature (London), 311, 238–241.Google Scholar
Hu, H. H., Li, F. H. & Fan, H. F. (1992). Crystal structure determination of K2O·7Nb2O5 by combining high resolution electron microscopy and electron diffraction. Ultramicroscopy, 41, 387–397.Google Scholar
Hurley, A. C. & Moodie, A. F. (1980). The inversion of three-beam intensities for scalar scattering by a general centrosymmetric crystal. Acta Cryst. A36, 737–738.Google Scholar
International Tables for Crystallography (2005). Vol. A. Space-group symmetry, edited by Th. Hahn, 5th ed. Heidelberg: Springer.Google Scholar
International Tables for Crystallography (2004). Vol. C. Mathematical, physical and chemical tables, edited E. Prince, 3rd ed. Dordrecht: Kluwer Academic Publishers.Google Scholar
International Tables for X-ray Crystallography (1952). Vol. I. Symmetry groups. Birmingham: Kynoch Press. (Present distributor Kluwer Academic Publishers, Dordrecht.)Google Scholar
Ishizuka, K., Miyazaki, M. & Uyeda, N. (1982). Improvement of electron microscope images by the direct phasing method. Acta Cryst. A38, 408–413.Google Scholar
Ishizuka, K. & Taftø, J. (1982). Kinematically allowed reflections caused by scattering via HOLZ. Proc. Electron Microsc. Soc. Am. pp. 688–689.Google Scholar
Jap, B. K. & Glaeser, R. M. (1980). The scattering of high-energy electrons. II. Quantitative validity domains of the single-scattering approximations for organic crystals. Acta Cryst. A36, 57–67.Google Scholar
Jap, B. K., Walian, P. J. & Gehring, K. (1991). Structural architecture of an outer membrane channel as determined by electron crystallography. Nature (London), 350, 167–170.Google Scholar
Johnson, A. W. S. (1972). Stacking faults in graphite. Acta Cryst. A28, 89–93.Google Scholar
Johnson, A. W. S. & Preston, A. R. (1994). Some notes on the selection of structural chirality by CBED. Ultramicroscopy, 55, 348–355.Google Scholar
Jones, P. M., Rackham, G. M. & Steeds, J. W. (1977). Higher order Laue zone effects in electron diffraction and their use in lattice parameter determination. Proc. R. Soc. London Ser. A, 354, 197–222.Google Scholar
Kam, Z. (1980). Three-dimensional reconstruction of aperiodic objects. J. Theor. Biol. 82, 15–32.Google Scholar
Kambe, K. (1982). Visualization of Bloch waves of high energy electrons in high resolution electron microscopy. Ultramicroscopy, 10, 223–228.Google Scholar
Karle, J. & Hauptman, H. (1956). A theory of phase determination for the four types of non-centrosymmetric space groups , , , . Acta Cryst. 9, 635–651.Google Scholar
Kirkland, E. J., Siegel, B. M., Uyeda, N. & Fujiyoshi, Y. (1980). Digital reconstruction of bright field phase contrast images from high resolution electron micrographs. Ultra-microscopy, 5, 479–503.Google Scholar
Kiselev, N. A., Lerner, F. Ya. & Livanova, N. B. (1971). Electron microscopy of muscle phosphorylase B. J. Mol. Biol. 62, 537–549.Google Scholar
Klug, A. & Berger, J. E. (1964). An optical method for the analysis of periodicities in electron micrographs and some observations on the mechanism of negative staining. J. Mol. Biol. 10, 565–569.Google Scholar
Klug, A. & DeRosier, D. J. (1966). Optical filtering of electron micrographs: reconstruction of one-sided images. Nature (London), 212, 29–32.Google Scholar
Kossel, W. & Möllenstedt, G. (1938). Electron interference in a convergent beam. Nature (London), 26, 660.Google Scholar
Kosykh, V. P., Pustovskikh, A. I., Kirichuk, V. S., Kühne, T., Orlova, E. V., Tsuprun, V. L. & Kiselev, N. A. (1983). Use of digital storage methods to recover images of monocrystalline layers of virus particles. Sov. Phys. Crystallogr. 28, 637–643.Google Scholar
Kühlbrandt, W., Wang, D. N. & Fujiyoshi, Y. (1994). Atomic model of plant light-harvesting complex by electron crystallography. Nature (London), 367, 614–621.Google Scholar
Kuwabara, S. (1978). Nearly aberration-free crystal images in high voltage electron microscopy. J. Electron Microsc. 27, 161–169.Google Scholar
Langer, R., Frank, J., Feltynowski, A. & Hoppe, W. (1970). Anwendung des Bilddifferenzverfahrens auf die Untersuchung von Strukturänderungen dünner Kohlefolien bei Elektronenbestrahlung. Ber. Bunsenges Phys. Chem. 74(11), 1120–1126.Google Scholar
Langs, D. A. & DeTitta, G. T. (1975). A flexible and rapid phase determination and refinement procedure. Acta Cryst. A31, S16.Google Scholar
Laue, M. von (1935). Die Fluoreszenzrontgenstrahlung von Einkristallen. Ann. Phys. (Leipzig), 23, 703–726.Google Scholar
Li, D. X. & Hovmöller, S. (1988). The crystal structure of Na3Nb12O31F determined by HREM and image processing. J. Solid State Chem. 73, 5–10.Google Scholar
Li, F. H. (1991). Crystal structures from high-resolution electron microscopy. In Electron crystallography of organic molecules, edited by J. R. Fryer & D. L. Dorset, pp. 153–167. Dordrecht: Kluwer Academic Publishers.Google Scholar
Liebman, G. (1955). A unified representation of magnetic electron lens properties. Proc. Phys. Soc. London Sect. B, 68, 737–745.Google Scholar
Liu, Y.-W., Fan, H.-F. & Zheng, C.-D. (1988). Image processing in high-resolution electron microscopy using the direct method. III. Structure-factor extrapolation. Acta Cryst. A44, 61–63.Google Scholar
Lobachev, A. N. & Vainshtein, B. K. (1961). An electron diffraction study of urea. Sov. Phys. Crystallogr. 6, 313–317.Google Scholar
Lynch, D. F. & Moodie, A. F. (1972). Numerical evaluation of low energy electron diffraction intensity. I. The perfect crystal with no upper layer lines and no absorption. Surf. Sci. 32, 422–438.Google Scholar
Lynch, D. F., Moodie, A. F. & O'Keefe, M. A. (1975). n-Beam lattice images. V. The use of the charge-density approximation in the interpretation of lattice images. Acta Cryst. A31, 300–307.Google Scholar
McLachlan, D. (1958). Crystal structure and information theory. Proc. Natl Acad. Sci. USA, 44, 948–956.Google Scholar
Mansfield, J. (1984). Convergent beam electron diffraction of alloy phases by the Bristol Group under the direction of John Steeds. Bristol: Adam Hilger.Google Scholar
Markham, R., Frey, S. & Hills, G. J. (1963). Methods for the enhancement of image detail and accentuation of structure in electron microscopy. Virology, 20, 88–102.Google Scholar
Matsuda, T., Tonomura, A. & Komada, T. (1978). Observation of lattice images with a field emission electron microscope. Jpn. J. Appl. Phys. 17, 2073–2074.Google Scholar
Mermin, N. D. (1992). The space groups of icosahedral quasicrystals and cubic, orthorhombic, monoclinic and triclinic crystals. Rev. Mod. Phys. 64, 3–49.Google Scholar
Mersereau, R. M. & Oppenheim, A. V. (1974). Digital reconstruction of multi-dimensional signals from their projections. Proc. IEEE, 62(10), 1319–1338.Google Scholar
Miyake, S. & Uyeda, R. (1950). An exception to Friedel's law in electron diffraction. Acta Cryst. 3, 314.Google Scholar
Mo, Y. D., Cheng, T. Z., Fan, H. F., Li, J. Q., Sha, B. D., Zheng, C. D., Li, F. H. & Zhao, Z. X. (1992). Structural features of the incommensurate modulation in the Pb-doped Bi-2223 high-Tc phase by defect method electron diffraction analysis. Supercond. Sci. Technol. 5, 69–72.Google Scholar
Moodie, A. F. (1965). Some structural implications of n-beam interactions. International Conference on Electron Diffraction and Crystal Defects, Melbourne, Australia, paper ID-1.Google Scholar
Moodie, A. F. (1972). Reciprocity and shape function in multiple scattering diagrams. Z. Naturforsch. Teil A, 27, 437–440.Google Scholar
Moodie, A. F. & Whitfield, H. J. (1984). CBED and HREM in the electron microscope. Ultramicroscopy, 13, 265–278.Google Scholar
Moss, B. & Dorset, D. L. (1982). Effect of crystal bending on direct phasing of electron diffraction data from cytosine. Acta Cryst. A38, 207–211.Google Scholar
Ogawa, T., Moriguchi, S., Isoda, S. & Kobayashi, T. (1994). Application of an imaging plate to electron crystallography at atomic resolution. Polymer, 35, 1132–1136.Google Scholar
Orlov, S. S. (1975). Theory of three-dimensional reconstruction. II. The recovery operator. Sov. Phys. Crystallogr. 20, 429–433.Google Scholar
Ottensmeyer, F. P., Andrews, J. W., Basett-Jones, D. P., Chan, A. S. & Hewitt, J. (1977). Signal to noise enhancement in dark field electron micrographs of vasopressin: filtering of arrays of images in reciprocal space. J. Microsc. 109, 256–268.Google Scholar
Pan, M. & Crozier, P. A. (1993). Quantitative imaging and diffraction of zeolites using a slow-scan CCD camera. Ultramicroscopy, 52, 487–498.Google Scholar
Pérez, S. & Chanzy, H. (1989). Electron crystallography of linear polysaccharides. J. Electron Microsc. Tech. 11, 280–285.Google Scholar
Picture Processing and Digital Filtering (1975). Edited by T. S. Huang. Berlin: Springer-Verlag.Google Scholar
Pinsker, Z. G. (1953). Electron diffraction. London: Butterworth.Google Scholar
Pinsker, Z. G., Zvyagin, B. B. & Imamov, R. M. (1981). Principal results of electron-diffraction structural investigations. Sov. Phys. Crystallogr. 26, 669–674.Google Scholar
Pogany, A. P. & Turner, P. S. (1968). Reciprocity in electron diffraction and microscopy. Acta Cryst. A24, 103–109.Google Scholar
Pond, R. C. & Vlachavas, D. S. (1983). Bicrystallography. Proc. R. Soc. London Ser. A, 386, 95–143.Google Scholar
Portier, R. & Gratias, D. (1981). Diffraction symmetries for elastic scattering. In Electron microscopy and analysis. Inst. Phys. Conf. Ser. No. 61, pp. 275–278. Bristol, London: Institute of Physics.Google Scholar
Radermacher, M., McEwen, B. & Frank, J. (1987). Three-dimensional reconstruction of asymmetrical object in standard and high voltage electron microscopy. Proc. Microscop. Soc. Canada, XII Annual Meet., pp. 4–5.Google Scholar
Radi, G. (1970). Complex lattice potentials in electron diffraction calculated for a number of crystals. Acta Cryst. A26, 41–56.Google Scholar
Radon, J. (1917). Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten. (On the determination of functions from their integrals along certain manifolds). Ber. Verh. Saechs. Akad. Wiss. Leipzig Math. Phys. Kl. 69, 262–277.Google Scholar
Ramachandran, G. N. & Lakshminarayanan, A. V. (1971). Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. Proc. Natl Acad. Sci. USA, 68(9), 2236–2240.Google Scholar
Revol, J. F. (1991). Electron crystallography of radiation-sensitive polymer crystals. In Electron crystallography of organic molecules, edited by J. R. Fryer & D. L. Dorset, pp. 169–187. Dordrecht: Kluwer Academic Publishers. Google Scholar
Revol, J. F. & Manley, R. St. J. (1986). Lattice imaging in polyethylene single crystals. J. Mater. Sci. Lett. 5, 249–251.Google Scholar
Rez, P. (1978). In Electron diffraction 1927–1977, edited by P. J. Dobson, J. B. Pendry & C. J. Humphreys, pp. 61–67. Inst. Phys. Conf. Ser. No. 41. Bristol, London: Institute of Physics.Google Scholar
Rozenfeld, A. (1969). Picture processing by computer. New York: Academic Press.Google Scholar
Sayre, D. (1952). The squaring method: a new method for phase determination. Acta Cryst. 5, 60–65.Google Scholar
Sayre, D. (1980). Phase extension and refinement using convolutional and related equation systems. In Theory and practice of direct methods in crystallography, edited by M. F. C. Ladd & R. A. Palmer, pp. 271–286. NY: Plenum Press.Google Scholar
Scaringe, R. P. (1992). Crystallography in two dimensions: comparison of theory and experiment for molecular layers. Trans. Am. Crystallogr. Assoc. 28, 11–23.Google Scholar
Schapink, F. W., Forgany, S. K. E. & Buxton, B. F. (1983). The symmetry of convergent-beam electron diffraction patterns from bicrystals. Acta Cryst. A39, 805–813.Google Scholar
Scherzer, O. (1949). The theoretical resolution limit of the electron microscope. J. Appl. Phys. 20, 20–29.Google Scholar
Schiske, P. (1968). Zur Frage der Bildrekonstruktion durch Fokusreihen. 1 Y Eur. Reg. Conf. Electron Microsc. Rome, 1, 145–146.Google Scholar
Schwartzman, A., Goodman, P. & Johnson, A. W. S. (1996). IUCr XVII Congress and General Assembly, Seattle, Washington, USA, August 8–16, Collected Abstracts, p. C-54, Abstract PS02.03.18.Google Scholar
Sha, B.-D., Fan, H.-F. & Li, F.-H. (1993). Correction for the dynamical electron diffraction effect in crystal structure analysis. Acta Cryst. A49, 877–880.Google Scholar
Shechtman, D., Blech, I., Gratias, D. & Cahn, J. W. (1984). Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 53, 1951–1953.Google Scholar
Shoemaker, V. & Glauber, R. (1952). The Born approximation in electron diffraction. Nature (London), 170, 290–291.Google Scholar
Spence, J. C. H., O'Keefe, M. A. & Kolar, H. (1977). Image interpretation in crystalline germanium. Optik (Stuttgart), 49, 307–323.Google Scholar
Spence, J. C. H. & Zuo, J. M. (1992). Electron microdiffraction. New York: Plenum Press.Google Scholar
Steeds, J. W. (1979). Convergent beam electron diffraction. In Introduction to analytical electron microscopy, edited by J. J. Hren, J. I. Goldstein & D. C. Joy, pp. 387–422. New York: Plenum.Google Scholar
Steeds, J. W. (1983). Developments in convergent beam electron diffraction. Report to the Commission on Electron Diffraction of the International Union of Crystallography.Google Scholar
Steeds, J. W. & Evans, N. S. (1980). Practical examples of point and space group determination in convergent beam diffraction. Proc. Electron Microsc. Soc. Am. pp. 188–191.Google Scholar
Steeds, J. W., Rackham, G. M. & Shannon, M. D. (1978). On the observation of dynamically forbidden lines in two and three dimensional electron diffraction. In Electron diffraction 1927–1977. Inst. Phys. Conf. Ser. No. 41, pp. 135–139.Google Scholar
Steeds, J. W. & Vincent, R. (1983). Use of high symmetry zone axes in electron diffraction in determining crystal point and space groups. J. Appl. Cryst. 16, 317–324.Google Scholar
Steinkilberg, M. & Schramm, H. J. (1980). Eine verbesserte Drehkorrelations Methode für die Strukturbestimmung biologischer Macromoleküle durch Mittelung elektronenmikroskopischer Bilder. Hoppe–Seyler's Z. Physiol. Chem. 361, 1363–1369.Google Scholar
Stereochemical Applications of Gas-Phase Electron Diffraction (1988). Part A, edited by I. Hargittai & M. Hargittai. New York: VCH.Google Scholar
Tanaka, M. (1994). Convergent-beam electron diffraction. Acta Cryst. A50, 261–286.Google Scholar
Tanaka, M., Saito, P., Ueno, K. & Harada, Y. (1980). Large angle convergent-beam electron diffraction. J. Electron. Microsc. 29, 408–412.Google Scholar
Tanaka, M., Sekii, H. & Nagasawa, T. (1983). Space group determination by dynamic extinction in convergent beam electron diffraction. Acta Cryst. A39, 825–837.Google Scholar
Tanaka, M. & Terauchi, M. (1985). Convergent-beam electron diffraction. Tokyo: JEOL Ltd.Google Scholar
Tanaka, M., Terauchi, M. & Tsuda, K. (1994). Convergent-beam electron diffraction III. Tokyo: JEOL–Maruzen.Google Scholar
Thon, F. (1966). On the defocusing dependence of phase contrast in electron microscopical images. Z. Naturforsch. Teil A, 21, 476–478.Google Scholar
Tivol, W. F., Dorset, D. L., McCourt, M. P. & Turner, J. N. (1993). Voltage-dependent effect on dynamical scattering and the electron diffraction structure analysis of organic crystals: copper perchlorophthalocyanine. Microsc. Soc. Am. Bull. 23, 91–98.Google Scholar
Tournaire, M. (1962). Recent developments of the matrical and semi-reciprocal formulation in the field of dynamical theory. J. Phys. Soc. Jpn, 17, Suppl. B11, 98–100.Google Scholar
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.)Google Scholar
Tsuji, M. (1989). Electron microscopy. In Comprehensive polymer science, Vol. 1. Polymer characterization, edited by C. Booth & C. Price, pp. 785–840. Oxford: Pergamon PressGoogle Scholar
Turner, P. S. & Cowley, J. M. (1969). The effects of n-beam dynamical diffraction on electron diffraction intensities from polycrystalline materials. Acta Cryst. A25, 475–481.Google Scholar
Unwin, P. N. T. & Henderson, R. (1975). Molecular structure determination by electron microscopy of unstained crystalline specimens. J. Mol. Biol. 94, 425–440.Google Scholar
Uyeda, N., Kobayashi, T., Ishizuka, K. & Fujiyoshi, Y. (1978–1979). High voltage electron microscopy for image discrimination of constituent atoms in crystals and molecules. Chem. Scr. 14, 47–61.Google Scholar
Vainshtein, B. K. (1952). Dependence of electron scattering on the atomic number. Dokl. Akad. Nauk SSSR, 85, 1239–1242. (In Russian.)Google Scholar
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.)Google Scholar
Vainshtein, B. K. (1955). Elektronograficheskoe issledovanie diketopiperazina. Zh. Fiz. Khim. 29, 327–344.Google Scholar
Vainshtein, B. K. (1956). Structure analysis by electron diffraction. Moscow: Akad. Sci. USSR. [English edition (1964): Oxford: Pergamon Press.]Google Scholar
Vainshtein, B. K. (1964). Structure analysis by electron diffraction. Oxford: Pergamon Press.Google Scholar
Vainshtein, B. K. (1971a). The synthesis of projecting functions. Sov. Phys. Dokl. 16, 66–69.Google Scholar
Vainshtein, B. K. (1971b). Finding the structure of objects from projections. Sov. Phys. Crystallogr. 15, 781–787.Google Scholar
Vainshtein, B. K. (1978). Electron microscopical analysis of the three-dimensional structure of biological macromolecules. In Advances in optical and electron microscopy, Vol. 7, edited by V. E. Cosslett & R. Barer, pp. 281–377. London: Academic Press.Google Scholar
Vainshtein, B. K., Barynin, V. V. & Gurskaya, G. V. (1968). The hexagonal crystalline structure of catalase and its molecular structure. Sov. Phys. Dokl. 13, 838–841.Google Scholar
Vainshtein, B. K., D'yakon, I. A. & Ablov, A. V. (1971). Electron diffraction determination of the structure of copper DL-alaninate. Sov. Phys. Dokl. 15, 645–647.Google Scholar
Vainshtein, B. K. & Goncharov, A. B. (1986a). Determination of the spatial orientation of arbitrarily arranged identical particles of unknown structure from their projections. Sov. Phys. Dokl. 287, 278–283.Google Scholar
Vainshtein, B. K. & Goncharov, A. B. (1986b). Proceedings of the 11th International Congress on Electron Microscopy, Kyoto, Vol. 1, pp. 459–460.Google Scholar
Vainshtein, B. K. & Klechkovskaya, V. V. (1993). Electron diffraction by Langmuir–Blodgett films. Proc. R. Soc. London Ser. A, 442, 73–84.Google Scholar
Vainshtein, B. K. & Orlov, S. S. (1972). Theory of the recovery of functions from their projections. Sov. Phys. Crystallogr. 17, 213–216.Google Scholar
Vainshtein, B. K. & Orlov, S. S. (1974). General theory of direct 3D reconstruction. Proceedings of International Workshop, Brookhaven National Laboratory, pp. 158–164.Google Scholar
Van Heel, M. (1984). Multivariate statistical classification of noisy images (randomly oriented biological macromolecules). Ultramicroscopy, 13, 165–184.Google Scholar
Vilkov, L. V., Mastryukov, V. S. & Sadova, N. I. (1978). Determination of geometrical structure of free molecules. Leningrad: Khimiya. (In Russian.)Google Scholar
Vincent, R. & Exelby, D. R. (1991). Structure of metastable Al–Ge phases determined from HOLZ Patterson transforms. Philos. Mag. Lett. 63, 31–38.Google Scholar
Vincent, R. & Exelby, D. R. (1993). Structure of a metastable Al–Ge phase determined from large angle CBED patterns. Philos. Mag. B, 68, 513–528.Google Scholar
Vincent, R. & Midgley, P. A. (1994). Double conical beam rocking system for measurement of integrated electron diffraction intensities. Ultramicroscopy, 53, 271–282.Google Scholar
Voronova, A. A. & Vainshtein, B. K. (1958). An electron diffraction study of CuCl2·3Cu(OH)2. Sov. Phys. Crystallogr. 3, 445–451.Google Scholar
Watanabe, D., Uyeda, R. & Kogiso, M. (1968). An apparent variation of structure factors for electrons with accelerating voltage. An observation through Kikuchi patterns. Acta Cryst. A24, 249–250.Google Scholar
Wenk, H.-R., Downing, K. H., Ho, M.-S. & O'Keefe, M. A. (1992). 3D structure determination from electron-microscope images: electron crystallography of staurolite. Acta Cryst. A48, 700–716.Google Scholar
Wilson, A. J. C. (1949). The probability distribution of X-ray intensities. Acta Cryst. 2, 318–321.Google Scholar
Withers, R. L., Schmid, S. & Thompson, J. G. (1993). A composite modulated structure approach to the lanthanide oxide fluoride, uranium nitride fluoride and zirconium nitride fluoride solid-solution fields. Acta Cryst. B49, 941–951.Google Scholar
Wolff, P. M. de, Janssen, T. & Janner, A. (1981). The superspace groups for incommensurate crystal structures with a one-dimensional modulation. Acta Cryst. A37, 625–636.Google Scholar
Xiang, S.-B., Fan, H.-F., Wu, X.-J., Li, F.-H. & Pan, Q. (1990). Direct methods in superspace. II. The first application to an unknown incommensurate modulated structure. Acta Cryst. A46, 929–934.Google Scholar
Yao, J.-X. (1981). On the application of phase relationships to complex structures. XVIII. RANTAN – random MULTAN. Acta Cryst. A37, 642–644.Google Scholar
Zemlin, F., Reuber, E., Beckmann, E., Zeitler, E. & Dorset, D. L. (1985). Molecular resolution electron micrographs of monolamellar paraffin crystal. Science, 229, 461–462.Google Scholar
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.Google Scholar
Zhukhlistov, A. P. & Zvyagin, B. B. (1998). Crystal structure of lizardite 1T from electron diffractometry data. Crystallogr. Rep. 43, 950–955.Google Scholar
Zuo, J. M., Gjønnes, K. & Spence, J. C. H. (1989). A FORTRAN source listing for simulating three-dimensional CBED patterns with absorption by the Bloch wave method. J. Electron Microsc. Tech. 12, 29–55.Google Scholar
Zvyagin, B. B. (1967). Electron-diffraction analysis of clay mineral structures. New York: Plenum.Google Scholar
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. (In Russian).Google Scholar
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. (In Russian).Google Scholar