International
Tables for Crystallography Volume F Crystallography of biological macromolecules Edited by E. Arnold, D. M. Himmel and M. G. Rossmann © International Union of Crystallography 2012 
International Tables for Crystallography (2012). Vol. F, ch. 2.1, pp. 5257
Section 2.1.4. Basic diffraction physics^{a}Laboratory of Biophysical Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 
The scattering of an Xray beam by a crystal results from interaction between the electric component of the beam and the electrons in the crystal. The magnetic component has hardly any effect and can be disregarded.
If a monochromatic polarized beam hits an electron, the electron starts to oscillate in the direction of the electric vector of the incident beam (Fig. 2.1.4.1). This oscillating electron acts as the aerial of a transmitter and radiates Xrays with the same or lower frequency as the incident beam. The frequency change is due to the Compton effect: the photons of the incident beam collide with the electron and lose part of their energy. This is inelastic scattering, and the scattered radiation is incoherent with the incident beam. Compton scattering contributes to the background in a diffraction experiment. In elastic scattering, the scattered radiation has the same wavelength as the incident radiation, and this is the radiation responsible for the interference effects in diffraction. It was shown by Thomson that if the electron is completely free the following hold:
In terms of energy,The scattered energy per unit solid angle is
It was shown by Klein & Nishina (1929) [see also Heitler (1966)] that the scattering by an electron can be discussed in terms of the classical Thomson scattering if the quantum energy . This is not true for very short Xray wavelengths. For , and are exactly equal, but for , is 0.0243 times . Since wavelengths in macromolecular crystallography are usually in the range 0.8–2.5 Å, the classical approximation is allowed. It should be noted that:

This can be derived along classical lines by calculating the phase difference between the Xray beams scattered by each of the two electrons. A derivation based on quantum mechanics leads exactly to the same result by calculating the transition probability for the scattering of a primary quantum , given a secondary quantum (Heitler, 1966, p. 193). For simplification we shall give only the classical derivation here. In Fig. 2.1.4.2, a system of two electrons is drawn with the origin at electron 1 and electron 2 at position r. They scatter the incident beam in a direction given by the vector s. The direction of the incident beam is along the vector . The length of the vectors can be chosen arbitrarily, but for convenience they are given a length . The two electrons scatter completely independently of each other.

The black dots are electrons. The origin of the system is at electron 1; electron 2 is at position r. The electrons are irradiated by an Xray beam from the direction indicated by vector . The radiation scattered by the electrons is observed in the direction of vector s. Because of the path difference , scattered beam 2 will lag behind scattered beam 1 in phase. Reproduced with permission from Drenth (1999). Copyright (1999) SpringerVerlag. 
Therefore, the amplitudes of the scattered beams 1 and 2 are equal, but they have a phase difference resulting from the path difference between the beam passing through electron 2 and the beam passing through electron 1. The path difference is . Beam 2 lags behind in phase compared with beam 1, and with respect to wave 1 its phase angle iswhere .
From Fig. 2.1.4.3, it is clear that the direction of S is perpendicular to an imaginary plane reflecting the incident beam at an angle θ and that the length of S is given byThe total scattering from the twoelectron system is 1 + if the resultant amplitude of the waves from electrons 1 and 2 is set to 1. In an Argand diagram, the waves are represented by vectors in a twodimensional plane, as in Fig. 2.1.4.4(a).^{1} Thus far, the origin of the system was chosen at electron 1. Moving the origin to another position simply means an equal change of phase angle for all waves. Neither the amplitudes nor the intensities of the reflected beams change (Fig. 2.1.4.4b).

The direction of the incident wave is indicated by and that of the scattered wave by s. Both vectors are of length . A plane that makes equal angles with s and can be regarded as a mirror reflecting the incident beam. Reproduced with permission from Drenth (1999). Copyright (1999) SpringerVerlag. 
Electrons in an atom are bound by the nucleus and are – in principle – not free electrons.
However, to a good approximation, they can be regarded as such if the frequency of the incident radiation ν is greater than the natural absorption frequencies, , at the absorption edges of the scattering atom, or the wavelength of the incident radiation is shorter than the absorptionedge wavelength (Section 2.1.4.4). This is normally true for light atoms but not for heavy ones (Table 2.1.4.1).

If the electrons in an atom can be regarded as free electrons, the scattering amplitude of the atom is a real quantity, because the electron cloud has a centrosymmetric distribution, i.e. .
A small volume, , at r contains electrons, and at −r there are electrons. The combined scattering of the two volume elements, in units of the scattering of a free electron, isthis is a real quantity.
The scattering amplitude of an atom is called the atomic scattering factor f. It expresses the scattering of an atom in terms of the scattering of a single electron. f values are calculated for spherically averaged electrondensity distributions and, therefore, do not depend on the scattering direction. They are tabulated in IT C (2004) as a function of . The f values decrease appreciably as a function of (Fig. 2.1.4.5). This is due to interference effects between the scattering from the electrons in the cloud. In the direction , all electrons scatter in phase and the atomic scattering factor is equal to the number of electrons in the atom.
A plane of atoms reflects an Xray beam with a phase retardation of with respect to the scattering by a single atom. The difference is caused by the difference in path length from source (S) to atom (M) to detector (D) for the different atoms in the plane (Fig. 2.1.4.6). Suppose the plane is infinitely large. The shortest connection between S and D via the plane is S–M–D. The plane containing S, M and D is perpendicular to the reflecting plane, and the lines SM and MD form equal angles with the reflecting plane. Moving outwards from atom M in the reflecting plane, to P for instance, the path length S–P–D is longer. At the edge of the first Fresnel zone, the path is longer (Fig. 2.1.4.6). This edge is an ellipse with its centre at M and its major axis on the line of intersection between the plane SMD and the reflecting plane. Continuing outwards, many more elliptic Fresnel zones are formed. Clearly, the beams radiated by the many atoms in the plane interfere with each other. The situation is represented in the Argand diagram in Fig. 2.1.4.7. Successive Fresnel zones can be subdivided into an equal number of subzones. If the distribution of electrons is sufficiently homogeneous, it can be assumed that the subzones in one Fresnel zone give the same amplitude at D. Their phases are spaced at regular intervals and their vectors in the Argand diagram lie in a half circle. In the lower part of Fig. 2.1.4.7, this is illustrated for the first Fresnel zone. For the second Fresnel zone (upper part), the radius is slightly smaller, because the intensity radiated by more distant zones decreases (Kauzmann, 1957). Therefore, the sum of vectors pointing upwards is shorter than that of those pointing downwards, and the resulting scattered wave lags in phase behind the scattering by the atom at M.

S is the Xray source and D is the detector. The scattering is by the atoms in a plane. The shortest distance between S and D via a point in the plane is through M. Path lengths via points in the plane further out from M are longer, and when these beams reach the detector they lag behind in phase with respect to the MD beam. The plane is divided into zones, such that from one zone to the next the path difference is . 

Schematic picture of the Argand diagram for the scattering by atoms in a plane. All electrons are considered free. The vector of the incident beam points to the left. The atom at M (see Fig. 2.1.4.6) has a phase difference of π with respect to the incident beam. Subzones in the first Fresnel zone have the endpoints of their vectors on the lower half circle. For the next Fresnel zone, they are on the upper half circle, which has a smaller radius because the amplitude decreases gradually for subsequent Fresnel zones (Kauzmann, 1957). The sum of all vectors points down, indicating a phase lag of with respect to the beam scattered by the atom at M. 
In classical dispersion theory, the scattering power of an atom is derived by supposing that the atom contains dipole oscillators. In units of the scattering of a free electron, the scattering of an oscillator with eigen frequency and moderate damping factor was found to be a complex quantity:where is the frequency of the incident radiation [James, 1965; see also IT C (2004), equation (4.2.6.8) ]. When in equation (2.1.4.4), approaches unity, as is the case for scattering by a free electron; when , approaches zero, demonstrating the lack of scattering from a fixed electron. Only for does the imaginary part have an appreciable value.
Fortunately, quantum mechanics arrives at the same result by adding a rational meaning to the damping factors and interpreting as absorption frequencies of the atom (Hönl, 1933). For heavy atoms, the most important transitions are to a continuum of energy states, with or etc., where and are the frequencies of the K and L absorption edges.
In practice, the complex atomic scattering factor, , is separated into three parts: . f is the contribution to the scattering if the electrons are free electrons and it is a real number (Section 2.1.4.3). f′ is the real part of the correction to be applied and f″ is the imaginary correction; f″ is always in phase ahead of f (Fig. 2.1.4.8). is the total real part of the atomic scattering factor.

The atomic scattering factor as a vector in the Argand diagram. (a) When the electrons in the atom can be regarded as free. (b) When they are not completely free and the scattering becomes anomalous with a real anomalous contribution and an imaginary contribution . Reproduced with permission from Drenth (1999). Copyright (1999) SpringerVerlag. 
The imaginary correction is connected with absorption by oscillators having . It can be calculated from the atomic absorption coefficient of the anomalously scattering element. For each of the K, L etc. absorption edges, is virtually zero for frequencies below the edge, but it rises steeply at the edge and decreases gradually at higher frequencies.
The real correction can be derived from by means of the Kramers–Kronig transform [IT C (2004), Section 4.2.6.2.2 ]. For frequencies close to an absorption edge, becomes strongly negative.
Values for f, and are always given in units equal to the scattering by one free electron. f values are tabulated in IT C (2004) as a function of , and the anomalousscattering corrections for forward scattering as a function of the wavelength. Because the anomalous contribution to the atomic scattering factor is mainly due to the electrons close to the nucleus, the value of the corrections diminishes much more slowly than f as a function of the scattering angle.
A unit cell contains a large number of electrons, especially in the case of biological macromolecules. The waves scattered by these electrons interfere with each other, thereby reducing the effective number of electrons in the scattered wave. The exception is scattering in the forward direction, where the beams from all electrons are in phase and add to each other. The effective number of scattering electrons is called the structure factor F because it depends on the structure, i.e. the distribution of the atoms in the unit cell. It also depends on the scattering direction. If small electrondensity changes due to chemical bonding are neglected, the structure factor can be regarded as the sum of the scattering by the atoms in the unit cell, taking into consideration their positions and the corresponding phase differences between the scattered waves. For n atoms in the unit cellwhere S is a vector perpendicular to the plane reflecting the incident beam at an angle θ; the length of S is given by [equation (2.1.4.3) in Section 2.1.4.2].
The origin of the system is chosen at the origin of the selected unit cell. Atom j is at position with respect to the origin. Another unit cell has its origin at and , where t, u and v are whole numbers, and a, b and c are the basis vectors of the unit cell. With respect to the first origin, its scattering is
The wave scattered by a crystal is the sum of the waves scattered by all unit cells. Assuming that the crystal has a very large number of unit cells , the amplitude of the wave scattered by the crystal is
For an infinitely large crystal, the three summations over the exponential functions are delta functions. They have the property that they are zero unlesswhere h, k and l are whole numbers, either positive, negative or zero. These are the Laue conditions. If they are fulfilled, all unit cells scatter in phase and the amplitude of the wave scattered by the crystal is proportional to the amplitude of the structure factor F. Its intensity is proportional to .
S vectors satisfying equation (2.1.4.7) are denoted by S(hkl) or S(h), and the corresponding structure factors as or F(h).
Bragg's law for scattering by a crystal is better known than the Laue conditions:where d is the distance between reflecting lattice planes, θ is the reflecting or glancing angle and λ is the wavelength (Fig. 2.1.4.9). It can easily be shown that the Laue conditions and Bragg's law are equivalent by combining equation (2.1.4.7) with the following information:

From equation (2.1.4.9) it follows that vector S(hkl) is perpendicular to a plane determined by the points a/h, b/k and c/l, and according to conditions (3) this is a lattice plane. Therefore, scattering by a crystal can indeed be regarded as reflection by lattice planes. The projection of a/h, b/k and c/l on vector S(hkl) is (Laue condition), but it is also equal to the spacing between the lattice planes (see Fig. 2.1.1.3), and, therefore, . Combining this with equation (2.1.4.3) yields Bragg's law, [equation (2.1.4.8)].
For noncentrosymmetric structures, the structure factor,is an imaginary quantity and can also be written as^{2}
It is sometimes convenient to split the structure factor into its real part, A(S), and its imaginary part, B(S). For centrosymmetric structures, if the origin of the structure is chosen at the centre of symmetry.
The average value of the structurefactor amplitude decreases with increasing or, because , with increasing reflecting angle θ.
This is caused by two factors:

In protein crystal structures determined at high resolution, each atom is given its own individual thermal parameter B.^{3} Anisotropic thermal vibration is described by six parameters instead of one, and the evaluation of this anisotropic thermal vibration requires more data (Xray intensities) than are usually available. Only at very high resolution (better than 1.5 Å) can one consider the incorporation of anisotropic temperature factors.
The value of can be regarded as the effective number of electrons per unit cell scattering in the direction corresponding to S. This is true if the values of are on an absolute scale; this means that the unit of scattering is the scattering by one electron in a specific direction. The experimental values of are normally on an arbitrary scale. The average value of the scattered intensity, , on an absolute scale is , where is the atomic scattering factor reduced by the temperature factor. This can be understood as follows:
For a large number of reflections, S varies considerably, and assuming that the angles are evenly distributed over the range 0–2π for , the average value for the terms with will be zero and only the terms with remain, giving
Because of the thermal vibrationswhere i denotes a specific atom and is the scattering factor for the atom i at rest.
It is sometimes necessary to transform the intensities and the structure factors from an arbitrary to an absolute scale. Wilson (1942) proposed a method for estimating the required scale factor K and, as an additional bonus, the thermal parameter B averaged over the atoms:
To determine K and B, equation (2.1.4.11) is written in the form
Because depends on , average intensities, , are calculated for shells of narrow ranges. is plotted against . The result should be a straight line with slope , intersecting the vertical axis at ln K (Fig. 2.1.4.10).

The Wilson plot for phospholipase A_{2} with data to 1.7 Å resolution. Only beyond 3 Å resolution is it possible to fit the curve to a straight line. Reproduced with permission from Drenth (1999). Copyright (1999) SpringerVerlag. 
For proteins, the Wilson plot gives rather poor results because the assumption in deriving equation (2.1.4.11) that the angles, , are evenly distributed over the range 0–2π for is not quite valid, especially not in the ranges at low resolution.
As discussed above, the average value of the structure factors, F(S), decreases with the scattering angle because of two effects:
This decrease is disturbing for statistical studies of structurefactor amplitudes. It is then an advantage to eliminate these effects by working with normalized structure factors, E(S), defined by
The application of equation (2.1.4.14) to gives
The average value, , is equal to 1. The advantage of working with normalized structure factors is that the scaling is not important, because if equation (2.1.4.14) is written asa scale factor affects numerator and denominator equally.
In practice, the normalized structure factors are derived from the observed data as follows:where is a correction factor for spacegroup symmetry. For general reflections it is 1, but it is greater than 1 for reflections having h parallel to a symmetry element. This can be understood as follows. For example, if m atoms are related by this symmetry element, (with j from 1 to m) is the same in their contribution to the structure factor
They act as one atom with scattering factor rather than as m different atoms, each with scattering factor f. According to equation (2.1.4.11), this increases by a factor on average. To make the F values of all reflections statistically comparable, F(h) must be divided by . For a detailed discussion, see IT B (2008), Chapter 2.1 , by U. Shmueli and A. J. C. Wilson.
References
International Tables for Crystallography (2008). Vol. B. Reciprocal Space, edited by U. Shmueli, 3rd ed. Heidelberg: Springer.International Tables for Crystallography (2004). Vol. C. Mathematical, Physical and Chemical Tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.
Heitler, W. G. (1966). The Quantum Theory of Radiation, 3rd ed. Oxford University Press.
Hönl, H. (1933). Atomfaktor für Röntgenstrahlen als Problem der Dispersionstheorie (KSchale). Ann. Phys. 18, 625–655.
James, R. W. (1965). The Optical Principles of the Diffraction of Xrays, p. 135. London: G. Bell and Sons Ltd.
Kauzmann, W. (1957). Quantum Chemistry. New York: Academic Press.
Klein, O. & Nishina, Y. (1929). Über die Streuung von Strahlung durch freie Elektronen nach der neuen relativistischen Quantendynamik von Dirac. Z. Phys. 52, 853–868.
Wilson, A. J. C. (1942). Determination of absolute from relative Xray intensity data. Nature (London), 150, 151–152.