International
Tables for
Crystallography
Volume C
Mathematical, physical and chemical tables
Edited by E. Prince

International Tables for Crystallography (2006). Vol. C. ch. 4.3, p. 259

Section 4.3.1.1. Elastic scattering from a perfect crystal

J. M. Cowleyb

4.3.1.1. Elastic scattering from a perfect crystal

| top | pdf |

The most important interaction of electrons with crystalline matter is the interaction with the electrostatic potential field. The scattering into sharp, Bragg reflections is considered in terms of the interaction of an incident plane wave with a time-independent, averaged, periodic potential field which may be written [\varphi({\bf r}) ={1\over \Omega}\sum_{\bf h}V({\bf h})\exp\{-2\pi i {\bf h}\cdot{\bf r}\}, \eqno (4.3.1.1)]where [\Omega] is the unit-cell volume and the Fourier coefficients, [V({\bf h})], may be referred to as the structure amplitudes corresponding to the reciprocal-lattice vectors h. In conformity with the crystallographic sign convention used throughout this volume [see also Volume B (IT B, 2001[link])], we choose a free-electron approximation for the incident electron beam of the form [\exp(-i{\bf k}\cdot {\bf r})] and the interaction is represented by inserting the potential (4.3.1.1)[link] in the Schrödinger wave equation [\nabla^2\psi(r)+2k\sigma\{E+\varphi(r)\}\psi(r)=0, \eqno (4.3.1.2)]where eE is the kinetic energy of the incident beam, [k\,(=2\pi/\lambda)] is the magnitude of the wavevector for the incident electrons, and σ is an `interaction constant' defined by [\sigma=2\pi me\lambda/h^2, \eqno (4.3.1.3)]where h is Planck's constant. Relativistic values of m and λ are assumed (see Subsection 4.3.1.4[link]).

The solution of equation (4.3.1.2)[link], subject to the boundary conditions imposed by the need to fit the waves in the crystal with the incoming and outgoing waves in vacuum at the crystal surfaces, then allows the directions and amplitudes of the diffracted beams to be obtained in terms of the crystal periodicities and the Fourier coefficients, [V({\bf h})], of [\varphi({\bf r})] by the eigenvalue or Bloch-wave method (Bethe, 1928[link]).

Alternatively, the scattered amplitudes may be obtained from the integral form of (4.3.1.2)[link], [\eqalignno{ \psi({\bf r})&=\exp\{-i{\bf k}_0\cdot{\bf r}\} \cr&\quad +K\int\displaystyle {\exp\{-ik|{\bf r}-{\bf r}'|\}\over |{\bf r}-{\bf r}'|}\varphi({\bf r}')\psi({\bf r}')\,{\rm d}\tau_{\bf r'}, &(4.3.1.4)}]where [\exp\{-i{\bf k}_0\cdot{\bf r}\}] represents the incident beam, K = σ/λ, and the integral is taken over the space of the variable, [{\bf r}']. An iterative solution of (4.3.1.4)[link] leads to the Born series, [\psi=\psi_0+\psi_1+\psi_2+\ldots,]where [\psi_0= \exp\{-i{\bf k}_0\cdot{\bf r}\}]and [\psi_n({\bf r})=K \int{\exp\{-ik|{\bf r}-{\bf r}'|\}\over|{\bf r}-{\bf r}'|}\varphi({\bf r}')\psi_{n-1}({\bf r}')\,{\rm d}\tau_{\bf r'}, \eqno (4.3.1.5)]for [n\ge1]. Terms of the series for [n=1,2,\ldots] may be considered to represent the contributions from single, double and multiple scattering of the incident electron beam. This method has been applied to the diffraction from crystals by Fujiwara (1959[link]).

A further formulation of the scattering problem in integral form is that due to Cowley & Moodie (1957[link]) who considered the progressive modification of an incident plane wave as it passed through successive thin slices of a crystal. The effect of the nth slice on the incident electron wave is that of a phase-grating so that the wavefunction is modified by multiplication by a transmission function, [q_n(xy)=\exp\{-i\sigma\varphi_n(xy)\}, \eqno (4.3.1.6)]where [\varphi_n(xy)] is the projection of the potential distribution within the slice in the direction of the incident beam, taken to be the z axis; [\varphi_n(x,y)=\textstyle\int\limits^{z_n+\Delta z}_{z_n}\,\varphi(x,y,z)\,{\rm d} z. \eqno (4.3.1.7)]Propagation of the wave between the centres of slices is represented by convolution with a propagation function, p(xy), so that the wave entering the (n + 1)th slice may be written [\psi_{n+1}(xy)=[\psi_n(xy)\cdot q_n(xy)]*p_n(xy). \eqno (4.3.1.8)]In the small-angle approximation, the function [p_n(xy)] is given by the usual Fresnel diffraction theory as [p(xy)=(i/\lambda\Delta z)\exp\{-ik(x^2+y^2)/2\Delta z\}. \eqno (4.3.1.9)]

In the limit that the slice thickness, [\Delta z], tends to zero, the iteration of (4.3.1.8)[link] gives an exact account of the diffraction by the crystal.

On the basis of the above-mentioned and other related formulations of the diffraction problem, several computing methods have been devised for calculation of the amplitudes and intensities of the many diffracted beams of appreciable intensity that may emerge simultaneously from a crystal (see Section 4.3.6[link]). In this way, a degree of accuracy may be achieved in the calculation of the intensities of spots in diffraction patterns or of the contrast in electron-microscope images of crystals (Section 4.3.8[link]).

References

First citation Bethe, H. A. (1928). Theorie der Beugung von Elektronen an Kristallen. Ann. Phys. Leipzig, 87, 55–129.Google Scholar
First citation 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
First citation Fujiwara, K. (1959). Application of higher-order Born approximation to multiple elastic scattering of electrons in crystals. J. Phys. Soc. Jpn, 14, 1513–1524.Google Scholar
First citation International Tables for Crystallography (2001). Vol. B, 2nd ed. Dordrecht: Kluwer Academic Publishers.Google Scholar








































to end of page
to top of page