Tables for
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2010). Vol. B, ch. 5.2, pp. 647-653   | 1 | 2 |

Chapter 5.2. Dynamical theory of electron diffraction

A. F. Moodie,a J. M. Cowleyb and P. Goodmanc

aDepartment of Applied Physics, Royal Melbourne Institute of Technology, 124 La Trobe Street, Melbourne, Victoria 3000, Australia,bArizona State University, Box 871504, Department of Physics and Astronomy, Tempe, AZ 85287–1504, USA, and cSchool of Physics, University of Melbourne, Parkville, Australia

In this chapter, the main formulations of the dynamical theory of electron diffraction are outlined. These include the defining equations, forward scattering, the evolution operator, the projection approximation, semi-reciprocal space, the two-beam approximation, the eigenvalue approach, translational invariance, dispersion surfaces, the multislice formulation, the Born series and other approximations.

5.2.1. Introduction

| top | pdf |

Since electrons are charged, they interact strongly with matter, so that the single scattering approximation has a validity restricted to thin crystals composed of atoms of low atomic number. Further, at energies of above a few tens of keV, the wavelength of the electron is so short that the geometry of two-beam diffraction can be approximated in only small unit cells.

It is therefore necessary to develop a scattering theory specific to electrons and, preferably, applicable to imaging as well as to diffraction. The development, started by Born (1926[link]) and Bethe (1928[link]), and continuing into the present time, is the subject of an extensive literature, which includes reviews [for instance: Howie (1978[link]), Humphreys (1979[link])] and historical accounts (Goodman, 1981[link]), and is incorporated in Chapter 5.1[link] . Here, an attempt will be made to present only that outline of the main formulations which, it is hoped, will help the nonspecialist in the use of the tables. No attempt will be made to follow the historical development, which has been tortuous and not always logical, but rather to seek the simplest and most transparent approach that is consistent with brevity. Only key points in proofs will be sketched in an attempt to display the nature, rather than the rigorous foundations of the arguments.

5.2.2. The defining equations

| top | pdf |

No many-body effects have yet been detected in the diffraction of fast electrons, but the velocities lie well within the relativistic region. The one-body Dirac equation would therefore appear to be the appropriate starting point. Fujiwara (1962[link]), using the scattering matrix, carried through the analysis for forward scattering, and found that, to a very good approximation, the effects of spin are negligible, and that the solution is the same as that obtained from the Schrödinger equation provided that the relativistic values for wavelength and mass are used. In effect a Klein–Gordon equation (Messiah, 1965[link]) can be used in electron diffraction (Buxton, 1978[link]) in the form [\nabla^{2} \psi_{\rm b} + {8\pi^{2} m|e|\varphi\over h^{2}} \psi_{\rm b} + {8 \pi^{2} m_{0}|e|W\over h^{2}} \left(1 + {|e|W\over 2m_{0}c^{2}}\right) \psi_{\rm b} = 0.]Here, W is the accelerating voltage and [\varphi], the potential in the crystal, is defined as being positive. The relativistic values for mass and wavelength are given by [m = m_{0} (1 - v^{2} / c^{2})^{-1/2}], and taking `e' now to represent the modulus of the electronic charge, [|e|], [\lambda = h[2m_{0}eW (1 + eW / 2m_{0}c^{2})]^{-1/2},]and the wavefunction is labelled with the subscript b in order to indicate that it still includes back scattering, of central importance to LEED (low-energy electron diffraction).

In more compact notation, [[\nabla^{2} + k^{2} (1 + \varphi / W)] \psi_{\rm b} = (\nabla^{2} + k^{2} + 2 k\sigma \varphi) \psi_{\rm b} = 0. \eqno(] Here [k = |{\bf k}|] is the scalar wavenumber of magnitude [2\pi / \lambda], and the interaction constant [\sigma = 2\pi me\lambda / h^{2}]. This constant is approximately [10^{-3}] for 100 kV electrons.

For fast electrons, [\varphi / W] is a slowly varying function on a scale of wavelength, and is small compared with unity. The scattering will therefore be peaked about the direction defined by the incident beam, and further simplification is possible, leading to a forward-scattering solution appropriate to HEED (high-energy electron diffraction).

5.2.3. Forward scattering

| top | pdf |

A great deal of geometric detail can arise at this point and, further, there is no generally accepted method for approximation, the various procedures leading to numerically negligible differences and to expressions of precisely the same form. Detailed descriptions of the geometry are given in the references.

The entrance surface of the specimen, in the form of a plate, is chosen as the x, y plane, and the direction of the incident beam is taken to be close to the z axis. Components of the wavevector are labelled with suffixes in the conventional way; [{\bf K}_{0} = {\bf k}_{x} + {\bf k}_{y}] is the transverse wavevector, which will be very small compared to [{\bf k}_{z}]. In this notation, the excitation error for the reflection is given by [\zeta_{\bf h} = {K_{0}^{2} - |{\bf K}_{0} + 2\pi {\bf h}|^{2}\over 4\pi |{\bf k}_{z}|}.]

An intuitive method argues that, since [\varphi / W \ll 1], then the component of the motion along z is little changed by scattering. Hence, making the substitution [\psi_{\rm b} = \psi \exp \{ik_{z}z\}] and neglecting [\partial^{2} \psi / \partial z^{2}], equation ([link] becomes [{\partial \psi\over \partial z} = i \left[{1\over 2k_{z}} (\nabla_{x, \, y}^{2} + K_{0}^{2}) + \sigma \varphi \right] \psi, \eqno(]where [\nabla_{x, \, y}^{2} \equiv {\partial^{2}\over \partial x^{2}} + {\partial^{2}\over \partial y^{2}},]and [\psi (x, y, 0) = \exp \{i(k_{x}x + k_{y}y)\}].

Equation ([link] is of the form of a two-dimensional time-dependent Schrödinger equation, with the z coordinate replacing time. This form has been extensively discussed. For instance, Howie (1966[link]) derived what is essentially this equation using an expansion in Bloch waves, Berry (1971[link]) used a Green function in a detailed and rigorous derivation, and Goodman & Moodie (1974[link]), using methods due to Feynman, derived the equation as the limit of the multislice recurrence relation. A method due to Corones et al. (1982[link]) brings out the relationship between the HEED and LEED equations. Equation ([link] is cast in the form of a first-order system, [{\partial\over \partial z} \pmatrix{\psi_{\rm b}\cr \displaystyle{\partial \psi_{\rm b}\over \partial z}\cr} = \pmatrix{0 &1\cr -(\nabla_{x, \, y}^{2} + k^{2} + 2k\sigma \varphi) &0\cr} \pmatrix{\psi_{\rm b}\cr \displaystyle{\partial \psi_{\rm b}\over \partial z}}.]

A splitting matrix is introduced to separate the wavefunction into the forward and backward components, [\psi_{\rm b}^{\pm}], and the fast part of the phase is factored out, so that [\psi_{\rm b}^{\pm} = \psi^{\pm} \exp \{\pm ik_{z}z\}]. In the resulting matrix differential equation, the off-diagonal terms are seen to be small for fast electrons, and equation ([link] reduces to the pair of equations [{\partial \psi^{\pm}\over \partial z} = \pm i \left[{1\over 2k_{z}} (\nabla_{x, \, y}^{2} + K_{0}^{2}) + \sigma \varphi \right] \psi^{\pm}. \eqno(]The equation for [\psi^{\pm}] is the Lontovich & Fock (1946[link]) parabolic equation.

5.2.4. Evolution operator

| top | pdf |

Equation ([link] is a standard and much studied form, so that many techniques are available for the construction of solutions. One of the most direct utilizes the causal evolution operator. A recent account is given by Gratias & Portier (1983[link]).

In terms of the `Hamiltonian' of the two-dimensional system, [-{\bf H}(z) \equiv {1\over 2k_{z}} (\nabla_{x, \, y}^{2} + K_{0}^{2}) + \sigma \varphi,]the evolution operator [{\bf U}(z, z_{0})], defined by [\psi(z) = {\bf U}(z, z_{0}) \psi_{0}], satisfies [i {\partial\over \partial z} {\bf U}(z, z_{0}) = {\bf H}(z) {\bf U}(z, z_{0}), \eqno(]or [{\bf U}(z, z_{0}) = 1 - i \textstyle\int\limits_{z_{0}}^{z} {\bf U}(z, z_{1}) {\bf H}(z_{1})\,\hbox{d}z_{1}. \eqno(]

5.2.5. Projection approximation – real-space solution

| top | pdf |

Many of the features of the more general solutions are retained in the practically important projection approximation in which [\varphi (x, y, z)] is replaced by its projected mean value [\varphi_{\rm p} (x, y)], so that the corresponding Hamiltonian [{\bf H}_{\rm p}] does not depend on z. Equation ([link]) can then be solved directly by iteration to give [{\bf U}_{\rm p} (z, z_{0}) = \exp \{- i{\bf H}_{\rm p} (z - z_{0})\}, \eqno(]and the solution may be verified by substitution into equation ([link]).

Many of the results of dynamical theory can be obtained by expansion of equation ([link] as [{\bf U}_{\rm p} \equiv {\bf 1} - i{\bf H}_{\rm p} (z - z_{0}) + {i^{2}\over 2!} {\bf H}_{\rm p}^{2} (z - z_{0}) - \ldots,]followed by the direct evaluation of the differentials. Such expressions can be used, for instance, to explore symmetries in image space.

5.2.6. Semi-reciprocal space

| top | pdf |

In the derivation of electron-diffraction equations, it is more usual to work in semi-reciprocal space (Tournarie, 1962[link]). This can be achieved by transforming equation ([link] with respect to x and y but not with respect to z, to obtain Tournarie's equation [{\hbox{d}^{2} |U\rangle\over \hbox{d}z^{2}} = - {\bf M}_{\rm b}(z)|U\rangle. \eqno(]Here [|U\rangle] is the column vector of scattering amplitudes and [{\bf M}_{\rm b} (z)] is a matrix, appropriate to LEED, with k vectors as diagonal elements and Fourier coefficients of the potential as nondiagonal elements.

This equation is factorized in a manner parallel to that used on the real-space equation [equation ([link]] (Lynch & Moodie, 1972[link]) to obtain Tournarie's forward-scattering equation [{\hbox{d}|U^{\pm}\rangle\over \hbox{d}z} = \pm i{\bf M}^{\pm} (z) |U^{\pm}\rangle, \eqno(]where [\eqalign{ {\bf M}^{\pm} (z) &= \pm [{\bf K} + (1/2) {\bf K}^{-1} {V} (z)],\cr [K_{ij}] &= \delta_{ij} K_{i},}]and [[V_{ij}] = 2k_{z} {\textstyle\sum\limits_{l}} V_{i - j} \exp \{ - 2\pi ilz\},]where [V_{i}\equiv\sigma{}v_i] are the scattering coefficients and [v_{i}] are the structure amplitudes in volts. In order to simplify the electron-diffraction expression, the third crystallographic index `l' is taken to represent the periodicity along the z direction.

The double solution involving M of equation ([link]) is of interest in displaying the symmetry of reciprocity, and may be compared with the double solution obtained for the real-space equation [equation ([link]]. Normally the [{\bf M}^{+}] solution will be followed through to give the fast-electron forward-scattering equations appropriate in HEED. [{\bf M}^{-}], however, represents the equivalent set of equations corresponding to the z reversed reciprocity configuration. Reciprocity solutions will yield diffraction symmetries in the forward direction when coupled with crystal-inverting symmetries (Section 2.5.3[link] ).

Once again we set out to solve the forward-scattering equation ([link],b[link]) now in semi-reciprocal space, and define an operator [{\bf Q}(z)] [compare with equation ([link])] such that [|U_{z}\rangle = {\bf Q}_{z}| U_{0} \rangle \quad\hbox{with}\quad U_{0} = |0\rangle \semi]i.e., [{\bf Q}_{z}] is an operator that, when acting on the incident wavevector, generates the wavefunction in semi-reciprocal space.

Again, the differential equation can be transformed into an integral equation, and once again this can be iterated. In the projection approximation, with [{\bf M}] independent of z, the solution can be written as [{\bf Q}_{\rm p} = \exp \{i{\bf M}_{\rm p} (z - z_{0})\}.]A typical off-diagonal element is given by [V_{i-j}/\cos\theta_i], where [\theta_{i}] is the angle through which the beam is scattered. It is usual in the literature to find that [\cos \theta_{i}] has been approximated as unity, since even the most accurate measurements are, so far, in error by much more than this amount.

This expression for [{\bf Q}_{\rm p}] is Sturkey's (1957[link]) solution, a most useful relation, written explicitly as [|U\rangle = \exp \{i{\bf M}_{\rm p}{T}\}|0\rangle \eqno(]with T the thickness of the crystal, and [|0\rangle], the incident state, a column vector with the first entry unity and the rest zero. [{\bf S} = \exp \{i{\bf M}_{\rm p}{T}\}]is a unitary matrix, so that in this formulation scattering is described as rotation in Hilbert space.

5.2.7. Two-beam approximation

| top | pdf |

In the two-beam approximation, as an elementary example, equation ([link] takes the form [\pmatrix{u_{0}\cr u_{\bf h}\cr} = \exp \left\{i \pmatrix{0 &V^*({\bf h})\cr V({\bf h}) &K_{\bf h}\cr} T\right\} \pmatrix{0\cr 1\cr}. \eqno(]If this expression is expanded directly as a Taylor series, it proves surprisingly difficult to sum. However, the symmetries of Clifford algebra can be exploited by summing in a Pauli basis thus, [\eqalign{ &\exp \left\{i \pmatrix{0 &V^*({\bf h})\cr V({\bf h}) &K_{\bf h}\cr} T\right\}\cr &\quad = \exp \left\{i {K_{\bf h} T\over 2}\right\} {\bf E} \exp \left\{i \left({K_{\bf h}\over 2} \boldsigma_{3} + V^{\rm R} \boldsigma_{1} - V^{\rm I} \boldsigma_{2}\right) T\right\}.}]Here, the [\boldsigma_{i}] are the Pauli matrices [\displaylines{\boldsigma_{1} = \pmatrix{0 &1\cr 1 &0\cr},\,\boldsigma_{2} = \pmatrix{0 &i\cr -i &0\cr},\,\boldsigma_{3} = \pmatrix{-1 &0\cr 0 &1\cr},\cr {\bf E} = \pmatrix{1 &0\cr 0 &1\cr},}]and [V^{\rm R}], [V^{\rm I}] are the real and imaginary parts of the complex scattering coefficients appropriate to a noncentrosymmetric crystal, i.e. [V_{\bf h} = V^{\rm R} + iV^{\rm I}]. Expanding, [\eqalign{ &\exp \left \{i \left({K_{\bf h}\over 2} \boldsigma_{3} + V^{\rm R} \boldsigma_{1} - V^{\rm I} \boldsigma_{2}\right) T\right\}\cr &\quad = {\bf E} + i \left({K_{\bf h}\over 2} \boldsigma_{3} + V^{\rm R} \boldsigma_{1} - V^{\rm I} \boldsigma_{2}\right) T\cr &\qquad - {1\over 2} \left({K_{\bf h}\over 2} \boldsigma_{3} + V^{\rm R} \boldsigma_{1} - V^{\rm I} \boldsigma_{2}\right)^{2} T^{2} + \ldots,}]using the anti-commuting properties of [\boldsigma_{i}]: [\left.\matrix{\boldsigma_{i} \boldsigma_{j} + \boldsigma_{j} \boldsigma_{i} &= 0\cr \phantom{\sigma_{i} \sigma_{j}+ }\boldsigma_{i} \boldsigma_{i} &= 1\cr}\right\}]and putting [[(K_{\bf h}/2)^{2} +V({\bf h}) V^*({\bf h})] = \Omega], [{\bf M}_2=[(K_{\bf h}/2)\boldsigma_3] [+\ V^{\rm R}\boldsigma_1-V^{\rm I}\boldsigma_2]], so that [{\bf M}_{2}^{2} = \Omega {\bf E}] and [{\bf M}_{2}^{3} = \Omega {\bf M}_{2}], the powers of the matrix can easily be evaluated. They fall into odd and even series, corresponding to sine and cosine, and the classical two-beam approximation is obtained in the form [\eqalignno{ &{\bf Q}_{2} = \exp \{i (K_{\bf h}/2) T\} {\bf E} \left[(\cos \Omega^{1/2} T) {\bf E} + i \left({\sin \Omega^{1/2} T\over \Omega^{1/2}}\right) {\bf M}_{2}\right].\cr & &(}]

This result was first obtained by Blackman (1939[link]), using Bethe's dispersion formulation. Ewald and, independently, Darwin, each with different techniques, had, in establishing the theoretical foundations for X-ray diffraction, obtained analogous results (see Section 5.1.3[link] ).

The two-beam approximation, despite its simplicity, exemplifies some of the characteristics of the full dynamical theory, for instance in the coupling between beams. As Ewald pointed out, a formal analogy can be found in classical mechanics with the motion of coupled pendulums. In addition, the functional form [(\sin ax)/x], deriving from the shape function of the crystal emerges, as it does, albeit less obviously, in the N-beam theory.

This derivation of equation ([link] exhibits two-beam diffraction as a typical two-level system having analogies with, for instance, lasers and nuclear magnetic resonance and exhibiting the symmetries of the special unitary group SU(2) (Gilmore, 1974[link]).

5.2.8. Eigenvalue approach

| top | pdf |

In terms of the eigenvalues and eigenvectors, defined by [{\bf H}_{\rm p}|\,j\rangle = \gamma_{j}|\,j\rangle,]the evolution operator can be written as [{\bf U}(z, z_{0}) = \textstyle\int |\,j\rangle \exp \{\gamma_{j} (z - z_{0})\} \langle \,j|\, \hbox{d}j.]

This integration becomes a summation over discrete eigen states when an infinitely periodic potential is considered.

Despite the early developments by Bethe (1928[link]), an N-beam expression for a transmitted wavefunction in terms of the eigenvalues and eigenvectors of the problem was not obtained until Fujimoto (1959[link]) derived the expression [U_{\bf h} = {\textstyle\sum\limits_{j}} \psi_{0}^{\,j*} \psi_{\bf h}^{\,j} \exp \{-i2\pi \gamma_{j} T\}, \eqno(]where [\psi_{\bf h}^{\,j}] is the h component of the j eigenvector with eigenvalue [\gamma_{j}].

This expression can now be related to those obtained in the other formulations. For example, Sylvester's theorem (Frazer et al., 1963[link]) in the form [f({\bf M}) = {\textstyle\sum\limits_{j}} {\bf A}_{j}\,f(\gamma_{j})]when applied to Sturkey's solution yields [\boldPhi_{\bf h} = \exp (i{\bf M}_{\rm p}z) = \textstyle\sum {\bf P}_{j} \exp \{i2\pi \gamma_{j}z\}](Kainuma, 1968[link]; Hurley et al., 1978[link]). Here, the [{\bf P}_{j}] are projection operators, typically of the form [{\bf P}_{j} = \prod_{n \neq j} {({\bf M}_{\rm p} - {\bf E}\gamma_{n})\over \gamma_{j} - \gamma_{n}}.]On changing to a lattice basis, these transform to [\psi_{0}^{\,j*} \psi_{\bf h}^{\,j}].

Alternatively, the semi-reciprocal differential equation can be uncoupled by diagonalizing [{\bf M}_{\rm p}] (Goodman & Moodie, 1974[link]), a process which involves the solution of the characteristic equation [|{\bf M}_{\rm p} - \gamma_{j}{\bf E}| = 0. \eqno(]

5.2.9. Translational invariance

| top | pdf |

An important result deriving from Bethe's initial analysis, and not made explicit in the preceding formulations, is that the fundamental symmetry of a crystal, namely translational invariance, by itself imposes a specific form on wavefunctions satisfying Schrödinger's equation.

Suppose that, in a one-dimensional description, the potential in a Hamiltonian [{\bf H}_{t}(x)] is periodic, with period t. Then, [\varphi (x + t) = \varphi (x)]and [{\bf H}_{t}\psi (x) = {\bf E}\psi (x).]Now define a translation operator [\boldGamma f(x) = f(x + t),]for arbitrary [f(x)]. Then, since [\boldGamma \varphi (x) = \varphi (x)], and [\nabla^{2}] is invariant under translation, [\boldGamma {\bf H}_{t}(x) = {\bf H}_{t}(x)]and [\boldGamma {\bf H}_{t}(x)\psi (x) = {\bf H}_{t}(x + t)\psi (x + t) = {\bf H}_{t}(x) \boldGamma \psi (x).]

Thus, the translation operator and the Hamiltonian commute, and therefore have the same eigenfunctions (but not of course the same eigenvalues), i.e. [\boldGamma \psi (x) = \alpha \psi (x).]This is a simpler equation to deal with than that involving the Hamiltonian, since raising the operator to an arbitrary power simply increments the argument [\boldGamma^{m} \psi (x) = \psi (x + mt) = \alpha^{m} \psi (x).]But [\psi (x)] is bounded over the entire range of its argument, positive and negative, so that [|\alpha| = 1], and [\alpha] must be of the form [\exp\{i2\pi kt\}].

Thus, [\psi (x + t) = \boldGamma \psi (x) = \exp \{i2\pi kt\}\psi (x)], for which the solution is [\psi (x) = \exp \{i2\pi kt\} q(x)]with [q(x + t) = q(x)].

This is the result derived independently by Bethe and Bloch. Functions of this form constitute bases for the translation group, and are generally known as Bloch functions. When extended in a direct fashion into three dimensions, functions of this form ultimately embody the symmetries of the Bravais lattice; i.e. Bloch functions are the irreducible representations of the translational component of the space group.

5.2.10. Bloch-wave formulations

| top | pdf |

In developing the theory from the beginning by eigenvalue techniques, it is usual to invoke the periodicity of the crystal in order to show that the solutions to the wave equation for a given wavevector k are Bloch waves of the form [\psi = C({\bf r}) \exp \{i {\bf k} \cdot {\bf r}\},]where [C({\bf r})] has the periodicity of the lattice, and hence may be expanded in a Fourier series to give [\psi = {\textstyle\sum\limits_{\bf h}} C_{\bf h} ({\bf k}) \exp \{i ({\bf k} + 2\pi {\bf h}) \cdot {\bf r}\}. \eqno(]The [C_{\bf h} ({\bf k})] are determined by equations of consistency obtained by substitution of equation ([link] into the wave equation.

If N terms are selected in equation ([link] there will be N Bloch waves where wavevectors differ only in their components normal to the crystal surface, and the total wavefunction will consist of a linear combination of these Bloch waves. The problem is now reduced to the problem of equation ([link].

The development of solutions for particular geometries follows that for the X-ray case, Chapter 5.1[link] , with the notable differences that:

  • (1) The two-beam solution is not adequate except as a first approximation for particular orientations of crystals having small unit cells and for accelerating voltages not greater than about 100 keV. In general, many-beam solutions must be sought.

  • (2) For transmission HEED, the scattering angles are sufficiently small to allow the use of a small-angle forward-scattering approximation.

  • (3) Polarization effects are negligible except for very low energy electrons.

Humphreys (1979[link]) compares the action of the crystal, in the Bloch-wave formalism, with that of an interferometer, the incident beam being partitioned into a set of Bloch waves of different wavevectors. `As each Bloch wave propagates it becomes out of phase with its neighbours (due to its different wavevector). Hence interference occurs. For example, if the crystal thickness varies, interference fringes known as thickness fringes are formed.' For the two-beam case, these are the fringes of the pendulum solution referred to previously.

5.2.11. Dispersion surfaces

| top | pdf |

One of the important constructs of the Bloch-wave formalism is the dispersion surface, a plot of the permitted values of the z component of a Bloch wavevector against the component of the incident wavevector parallel to the crystal surface. The curve for a particular Bloch wave is called a branch. Thus, for fast electrons, the two-beam approximation has two branches, one for each eigenvalue, and the N-beam approximation has N.

A detailed treatment of the extensive and powerful theory that has grown from Bethe's initial paper is to be found, for example, in Hirsch et al. (1965[link]). Apart from its fundamental importance as a theoretical tool, this formulation provides the basis for one of the most commonly used numerical techniques, the essential step being the estimation of the eigenvalues from equation ([link] [see IT C (2004[link], Section[link] )].

5.2.12. Multislice

| top | pdf |

Multislice derives from a formulation that generates a solution in the form of a Born series (Cowley & Moodie, 1962[link]). The crystal is treated as a series of scattering planes on to which the potential from the slice between z and [z + \Delta z] is projected, separated by vacuum gaps [\Delta z], not necessarily corresponding to any planes or spacings of the material structure. The phase change in the electron beam produced by passage through a slice is given by [q = \exp \left \{- i\sigma \textstyle\int\limits_{z_{1}}^{z_{1} + \Delta z} \varphi (x, y, z)\, \hbox{d}z\right \},]and the phase distribution in the xy plane resulting from propagation between slices is given by [p = \exp \left \{{ik(x^{2} + y^{2})\over 2\Delta z}\right \},]where the wavefront has been approximated by a paraboloid. Thus, the wavefunction for the (n + 1)th slice is given by [\eqalignno{ \psi_{n + 1} &= \left[\psi_{n} * \exp \left \{{ik(x^{2} + y^{2})\over 2\Delta z}\right \}\right] \exp \{- i\sigma \varphi_{n + 1}\}\cr &= [\psi_{n} * p]q, &(}]where [\ast] is the convolution operator (Cowley, 1981[link]).

This equation can be regarded as the finite difference form of the Schrödinger equation derived by Feynman's (1948[link]) method. The calculation need be correct only to first order in [\Delta z]. Writing the convolution in equation ([link] explicitly, and expanding in a Taylor series, the integrals can be evaluated to yield equation ([link] (Goodman & Moodie, 1974[link]).

If equation ([link] is Fourier transformed with respect to x and y, the resulting recurrence relation is of the form [U_{n + 1} = [U_{n}P] * Q_{n}, \eqno(]where P and Q are obtained by Fourier transforming p and q above. This form is convenient for numerical work since, for a perfect crystal, it is: discrete, as distinct from equation ([link] which is continuous in the variables [see IT C (2004[link], Section[link] )]; numerically stable at least up to 5000 beams; fast; and only requires a computer memory proportional to the number of beams (Goodman & Moodie, 1974[link]).

5.2.13. Born series

| top | pdf |

In the impulse limit of equation ([link], the integrals can be evaluated to give the Born series (Cowley & Moodie, 1957[link])[U(h,k)=\textstyle\sum\limits_nU_n(h,k),]where[\eqalignno{U_n(h,k)&=\textstyle \sum \limits_{l}\textstyle \sum \limits_{h_1k_1l_1}\ldots\textstyle \sum \limits_{h_{n-1}k_{n-1}l_{n-1}}i^nV(h_1,k_1,l_1)&\cr &\quad\ldots V\left(h-\textstyle \sum \limits_{r=1}^{n-1}h_r, k-\textstyle \sum \limits_{r=1}^{n-1}k_r, l-\textstyle \sum \limits_{r=1}^{n-i}l_r\right)&\cr&\quad\times[\exp\{-2\pi i \zeta T\}/(2 \pi i )^n]&\cr &\quad\times\left(\exp\{i \pi \zeta T\}(\sin\pi\zeta T/\zeta)\left[(\zeta-\zeta_1)\ldots(\zeta-\zeta_{n-1})\right]^{-1} \right.&\cr&\left.\quad+\textstyle \sum \limits_{m=1}^{n-1}\exp\{i \pi \zeta_m T\}(\sin\pi\zeta_m T/\zeta_m)\left[(\zeta_m-\zeta_1)\right.\right.&\cr&\left.\left.\quad\ldots(\zeta_m-\zeta_{m-1})(\zeta_m-\zeta_{m+1})\ldots(\zeta_m-\zeta)\right]^{-1}\right)&\cr&&(\cr}]and where n is the order of interaction. Here [\zeta] is the excitation error of the reflection with index h, k, and [\zeta_i] are the excitation errors for the reflections with indices [h_i], [k_i], [l_i]. Thus each constituent process may be represented by a diagram, starting on the origin of reciprocal space, possibly looped, and ending on the point with coordinates (hk).

This solution can also be obtained by iteration of the Green-function integral equation, the integrals being evaluated by means of suitably chosen contours on the complex [k_{z}] plane (Fujiwara, 1959[link]), as well as by expansion of the scattering matrix (Fujimoto, 1959[link]).

Clearly, two or more of the [\zeta_i] will, in general, be equal in nearly all of the terms in equation ([link]). Confluence is, however, readily described, the divided differences of arbitrary order transforming into differentials of the same order (Moodie, 1972[link]).

The physical picture that emerges from equation ([link]) is that of n-fold scattering, the initial wave being turned through n − 1 intermediate states, processes that can be presented by scattering diagrams in reciprocal space (Gjønnes & Moodie, 1965[link]).

For a given scattering vector, constituent functions are evaluated for all possible paths in three dimensions, and those functions are then summed over l. There are therefore two distinct processes by which upper-layer lines can perturb wavefunctions in the zone, namely: by scattering out of the zone and then back in; and by intrusion of the effective shape function from another zone, the latter process being already operative in the first Born, or kinematical approximation.

The constituent functions to be evaluated can be transformed into many forms. One of the more readily described is that which assigns to each diagram an effective dynamical shape function. If there are no loops in the diagram of order n, this effective shape function is the (n + 1)th divided difference of the constituent phase-shifted kinematical shape transforms. For general diagrams, divided differences in loops are replaced by the corresponding differentials. The resulting function is multiplied by the convolution of the contributing structure amplitudes and diagrams of all orders summed (Moodie, 1972[link]).

While scattering diagrams have no utility in numerical work, they find application in the analysis of symmetries, for instance in the determination of the presence or absence of a centre of inversion [for a recent treatment, see Moodie & Whitfield (1995[link])] and in the detection of screw axes and glide planes (Gjønnes & Moodie, 1965[link]). Methods for the direct determination of all space groups are described by Goodman (1975[link]) and by Tanaka et al. (1983[link]) (see Section 2.5.3[link] ).

Equation ([link]) can be rewritten in a form particularly suited to the classification of approximations, and to describing the underlying symmetry of the formulation. The equation is written for compactness as[U_n(h)=E_n(h)Z_n(\zeta),]so that [E_n(h)] depends only on crystal structure and [Z_n(\zeta)] only on diffraction geometry. A transformation (Cowley & Moodie, 1962[link]) involving bialternants leads to[U_n=\textstyle \sum \limits_{r=0}^{\infty}E_n(h)\left[(2\pi i T)^{n+r}/(n+r)!\right]h_r(\zeta,\zeta_1\ldots\zeta_{n-1}),\eqno(]where [h_r] is the complete homogeneous symmetric polynomial function of n variables of order r.

Upper-layer-line effects can, of course, be calculated in any of the formulations.

5.2.14. Approximations

| top | pdf |

So far, only the familiar first Born and two-beam approximations and the projection approximation have been mentioned. Several others, however, have a considerable utility.

A high-voltage limit can be calculated in standard fashion to give [U_{\rm HVL} (h, k) = \hbox{\scr F} \exp \left\{ -i \sigma_{c} \textstyle\int\limits_{0}^{T} \varphi (x,y,z)\, \hbox{d}z\right\}, \eqno(]where [\hbox{\scr F}] is the Fourier transform operator, and [\sigma_{c} = 2\pi m_{0}e\lambda_{c} / h^{2}] with [\lambda_{c} = (h / m_{0}c)], the Compton wavelength. The phase-grating approximation, which finds application in electron microscopy, involves the assumption that equation ([link] has some range of validity when [\sigma_{c}] is replaced by [\sigma]. This is equivalent to ignoring the curvature of the Ewald sphere and can therefore apply to thin crystals [see Section 2.5.2[link] and IT C (2004[link], Section 4.3.8[link] )].

Approximations that involve curtailing the number of beams evidently have a range of validity that depends on the size of the unit cell. The most explored case is that of three-beam inter­actions. Kambe (1957[link]) has demonstrated that phase information can be obtained from the diffraction data; Gjønnes & Høier (1971[link]) analysed the confluent case, and Hurley & Moodie (1980[link]) have given an explicit inversion for the centrosymmetric case. Analyses of the symmetry of the defining differential equation, and of the geometry of the noncentrosymmetric case, have been given by Moodie et al. (1996[link], 1998[link]).

Niehrs and his co-workers (e.g. Blume, 1966[link]) have shown that, at or near zones, effective two-beam conditions can sometimes obtain, in that, for instance, the central beam and six equidistant beams of equal structure amplitude can exhibit two-beam behaviour when the excitation errors are equal. Group-theoretical treatments have been given by Fukuhara (1966[link]) and by Kogiso & Takahashi (1977[link]). Explicit reductions for all admissible noncentrosymmetric space groups have been obtained by Moodie & Whitfield (1994[link]). Extensions of such results have application in the interpretation of lattice images and convergent-beam patterns.

The approximations near the classical limit have been extensively explored [for instance, see Berry (1971[link])] but channelling has effectively become a separate subject and cannot be discussed here.


Berry, M. V. (1971). Diffraction in crystals at high voltages. J. Phys. C, 4, 697–722.
Bethe, H. A. (1928). Theorie der Beugung von Elektronen an Kristallen. Ann. Phys. (Leipzig), 87, 55–129.
Blackman, M. (1939). Intensities of electron diffraction rings. Proc. Phys. Soc. London Sect. A, 173, 68–72.
Blume, J. (1966). Die Kantenstreifung im Elektronen-Mikroskopischen Bild Wurfelformiger MgO Kristalle bei Durchstrahlung im Richtung der Raumdiagonal. Z. Phys. 191, 248–272.
Born, M. (1926). Quantenmechanik der Stossvorgange. Z. Phys. 38, 803–826.
Buxton, B. (1978). Graduate Lecture-Course Notes: Dynamical Diffraction Theory. Cambridge University, England.
Corones, J., De Facio, B. & Kreuger, R. J. (1982). Parabolic approximations to the time-independent elastic wave equation. J. Math. Phys. 23, 577–586.
Cowley, J. M. (1981). Diffraction Physics, pp. 26–30. Amsterdam: North-Holland.
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.
Cowley, J. M. & Moodie, A. F. (1962). The scattering of electrons by thin crystals. J. Phys. Soc. Jpn, 17, Suppl. B11, 86–91.
Feynman, R. (1948). Space–time approach to non-relativistic quantum mechanics. Rev. Mod. Phys. 201, 367–387.
Frazer, R. A., Duncan, W. J. & Collar, A. R. (1963). Elementary Matrices, pp. 78–79. Cambridge University Press.
Fujimoto, F. (1959). Dynamical theory of electron diffraction in the Laue case. J. Phys. Soc. Jpn, 14, 1558–1568.
Fujiwara, K. (1959). Application of higher order Born approximation to multiple elastic scattering of electrons by crystals. J. Phys. Soc. Jpn, 14, 1513–1524.
Fujiwara, K. (1962). Relativistic dynamical theory of electron diffraction. J. Phys. Soc. Jpn, 17, Suppl. B11, 118–123.
Fukuhara, A. (1966). Many-ray approximations in the dynamical theory of electron diffraction. J. Phys. Soc. Jpn, 21, 2645–2662.
Gilmore, R. (1974). Lie Groups, Lie Algebras, and Some of Their Applications. New York: Wiley–Interscience.
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.
Gjønnes, J. & Moodie, A. F. (1965). Extinction conditions in the dynamic theory of electron diffraction. Acta Cryst. 19, 65–67.
Goodman, P. (1975). A practical method of three-dimensional space-group analysis using convergent-beam electron diffraction. Acta Cryst. A31, 804–810.
Goodman, P. (1981). Editor. Fifty Years of Electron Diffraction. Dordrecht: Kluwer Academic Publishers.
Goodman, P. & Moodie, A. F. (1974). Numerical evaluation of N-beam wave functions in electron scattering by the multislice method. Acta Cryst. A30, 280–290.
Gratias, D. & Portier, R. (1983). Time-like perturbation method in high energy electron diffraction. Acta Cryst. A39, 576–584.
Hirsch, P. B., Howie, A., Nicholson, R. B., Pashley, D. W. & Whelan, M. J. (1965). Electron Microscopy of Thin Crystals. London: Butterworths.
Howie, A. (1966). Diffraction channelling of fast electrons and positrons in crystals. Philos. Mag. 14, 223–237.
Howie, A. (1978). In Electron Diffraction 1927–1977, edited by P. J. Dobson, J. B. Pendry & C. J. Humphreys, pp. 1–12. Inst. Phys. Conf. Ser. No. 41. Bristol/London: Institute of Physics.
Humphreys, C. J. (1979). The scattering of fast electrons by crystals. Rep. Prog. Phys. 42, 1825–1887.
Hurley, A. C., Johnson, A. W. S., Moodie, A. F., Rez, P. & Sellar, J. R. (1978). Algebraic approaches to N-beam theory. In Electron Diffraction 1927–1977, edited by P. J. Dobson, J. B. Pendry & C. J. Humphreys, pp. 34–40. Inst. Phys. Conf. Ser. No. 41. Bristol/London: Institute of Physics.
Hurley, A. C. & Moodie, A. F. (1980). The inversion of the three-beam intensities for scalar scattering by a general centrosymmetric crystal. Acta Cryst. A36, 737–738.
International Tables for Crystallography (2004). Vol. C. Mathematical, Physical and Chemical Tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.
Kainuma, Y. (1968). Averaged intensities in the many beam dynamical theory of electron diffraction. Part I. J. Phys. Soc. Jpn, 25, 498–510.
Kambe, K. (1957). Study of simultaneous reflection in electron diffraction by crystal. J. Phys. Soc. Jpn, 12, 13–31.
Kogiso, M. & Takahashi, H. (1977). Group-theoretical method in the many-beam theory of electron diffraction. J. Phys. Soc. Jpn, 42, 223–229.
Lontovitch, M. & Fock, R. (1946). Solution of the problem of propagation of electromagnetic waves along the Earth's surface by the method of parabolic equation. (Translated from Russian by J. Smorodinsky.) J. Phys. 10, 13–24.
Lynch, D. F. & Moodie, A. F. (1972). Numerical evaluation of low energy electron diffraction intensities. Surf. Sci. 32, 422–438.
Messiah, A. (1965). Quantum Mechanics, Vol. II, pp. 884–888. Amsterdam: North-Holland.
Moodie, A. F. (1972). Reciprocity and shape functions in multiple scattering diagrams. Z. Naturforsch. Teil A, 27, 437–440.
Moodie, A. F., Etheridge, J. & Humphreys, C. J. (1996). The symmetry of three-beam scattering equations: inversion of three-beam diffraction patterns from centrosymmetric crystals. Acta Cryst. A52, 596–605.
Moodie, A. F., Etheridge, J. & Humphreys, C. J. (1998). The Coulomb interaction and the direct measurement of structural phase. In The Electron, edited by A. Kirkland & P. Brown, pp. 235–246. IOM Communications Ltd. London: The Institute of Materials.
Moodie, A. F. & Whitfield, H. J. (1994). The reduction of N-beam scattering from noncentrosymmetric crystals to two-beam form. Acta Cryst. A50, 730–736.
Moodie, A. F. & Whitfield, H. J. (1995). Friedel's law and non-centrosymmetric space groups. Acta Cryst. A51, 198–201.
Sturkey, L. (1957). The use of electron-diffraction intensities in structure determination. Acta Cryst. 10, 858–859.
Tanaka, M., Sekii, H. & Nagasawa, T. (1983). Space-group determination by dynamic extinction in convergent-beam electron diffraction. Acta Cryst. A39, 825–837.
Tournarie, M. (1962). Recent developments of the matrical and semi-reciprocal formulation on the field of dynamical theory. J. Phys. Soc. Jpn, 17, Suppl. B11, 98–100.

to end of page
to top of page