International
Tables for Crystallography Volume B Reciprocal space Edited by U. Shmueli © International Union of Crystallography 2010 
International Tables for Crystallography (2010). Vol. B, ch. 2.3, pp. 244247
Section 2.3.1. Introduction ^{a}Department of Biological Sciences, Purdue University, West Lafayette, Indiana 47907, USA,^{b}CABM & Rutgers University, 679 Hoes Lane, Piscataway, New Jersey 08854–5638, USA, and ^{c}Department of Biological Sciences, Columbia University, New York 10027, USA 
Historically, the Patterson has been used in a variety of ways to effect the solutions of crystal structures. While some simple structures (Ketelaar & de Vries, 1939; Hughes, 1940; Speakman, 1949; Shoemaker et al., 1950) were solved by direct analysis of Patterson syntheses, alternative methods have largely superseded this procedure. An early innovation was the heavyatom method which depends on the location of a small number of relatively strong scatterers (Harker, 1936). Imageseeking methods and Patterson superposition techniques were first contemplated in the late 1930s (Wrinch, 1939) and applied sometime later (Beevers & Robertson, 1950; Clastre & Gay, 1950; Garrido, 1950a; Buerger, 1959). This experience provided the encouragement for computerized vectorsearch methods to locate individual atoms automatically (Mighell & Jacobson, 1963; Kraut, 1961; Hamilton, 1965; Simpson et al., 1965) or to position known molecular fragments in unknown crystal structures (Nordman & Nakatsu, 1963; Huber, 1965). The Patterson function has been used extensively in conjunction with the isomorphous replacement method (Rossmann, 1960; Blow, 1958) or anomalous dispersion (Rossmann, 1961a) to determine the position of heavyatom substitution. Pattersons have been used to detect the presence and relative orientation of multiple copies of a given chemical motif in the crystallographic asymmetric unit in the same or different crystals (Rossmann & Blow, 1962). Finally, the orientation and placement of known molecular structures (`molecular replacement') into unknown crystal structures can be accomplished via Patterson techniques.
The function, introduced by Patterson in 1934 (Patterson, 1934a,b), is a convolution of electron density with itself and may be defined as where is the `Patterson' function at u, is the crystal's periodic electron density and V is the volume of the unit cell. The Patterson function, or series, can be calculated directly from the experimentally derived Xray intensities as The derivation of (2.3.1.2) from (2.3.1.1) can be found in this volume (see Section 1.3.4.2.1.6 ) along with a discussion of the physical significance and symmetry of the Patterson function, although the principal properties will be restated here.
The Patterson can be considered to be a vector map of all the pairwise interactions between the atoms in a unit cell. The vectors in a Patterson correspond to vectors in the real (direct) crystal cell but translated to the Patterson origin. Their weights are proportional to the product of densities at the tips of the vectors in the real cell. The Patterson unit cell has the same size as the real crystal cell. The symmetry of the Patterson comprises the Laue point group of the crystal cell plus any additional lattice symmetry due to Bravais centring. The reduction of the real space group to the Laue symmetry is produced by the translation of all vectors to the Patterson origin and the introduction of a centre of symmetry. The latter is a consequence of the relationship between the vectors AB and BA. The Patterson symmetries for all 230 space groups are tabulated in IT A (2005).
An analysis of Patterson peaks can be obtained by considering N atoms with form factors in the unit cell. Then Using Friedel's law, which can be decomposed to On substituting (2.3.1.3) in (2.3.1.2), we see that the Patterson consists of the sum of total interactions of which N are of weight at the origin and are of weight at .
The weight of a peak in a real cell is given by where U is the volume of the atom i. By analogy, the weight of a peak in a Patterson (form factor ) will be given by Although the maximum height of a peak will depend on the spread of the peak, it is reasonable to assume that heights of peaks in a Patterson are proportional to the products of the atomic numbers of the interacting atoms.
There are a total of interactions in a Patterson due to N atoms in the crystal cell. These can be represented as an square matrix whose elements , indicate the position and weight of the peak produced between atoms i and j (Table 2.3.1.1). The N vectors corresponding to the diagonal of this matrix are located at the Patterson origin and arise from the convolution of each atom with itself. This leaves vectors whose locations depend on the relative positions of all of the atoms in the crystal cell and whose weights depend on the atom types related by the vector. Complete specification of the unique nonorigin Patterson vectors requires description of only the elements in either the upper or the lower triangle of this matrix, since the two sets of vectors represented by the two triangles are related by a centre of symmetry . Patterson vector positions are usually represented as , where u, v and w are expressed as fractions of the Patterson cell axes.

If we assume a constant number of atoms per unit volume, the number of atoms N in a unit cell increases in direct proportion with the volume of the unit cell. Since the number of nonorigin peaks in the Patterson function is and the Patterson cell is the same size as the real cell, the problem of overlapping peaks in the Patterson function becomes severe as N increases. To make matters worse, the breadth of a Patterson peak is roughly equal to the sum of the breadth of the original atoms. The effective width of a Patterson peak will also increase with increasing thermal motion, although this effect can be artificially reduced by sharpening techniques. Naturally, a loss of attainable resolution at high scattering angles will affect the resolution of atomic peaks in the real cell as well as peaks in the Patterson cell. If U is the van der Waals volume per average atom, then the fraction of the cell occupied by atoms will be . Similarly, the fraction of the cell occupied by Patterson peaks will be or . With the reasonable assumption that for a typical organic crystal, then the cell can contain at most five atoms for there to be no overlap, other than by coincidence, of the peaks in the Patterson. As N increases there will occur a background of peaks on which are superimposed features related to systematic properties of the structure.
The contrast of selected Patterson peaks relative to the general background level can be enhanced by a variety of techniques. For instance, the presence of heavy atoms not only enhances the size of a relatively small number of peaks but ordinarily ensures a larger separation of the peaks due to the lightatom skeleton on which the heavy atoms are hung. That is, the factor f (above) is substantially reduced. Another example is the effect of systematic atomic arrangements (e.g. αhelices or aromatic rings) resulting in multiple peaks which stand out above the background. In the isomorphous replacement method, isomorphous difference Pattersons are utilized in which the contrast of the Patterson interactions between the heavy atoms is enhanced by removal of the predominant interactions which involve the rest of the structure.
A. L. Patterson, in his first indepth exposition of his newly discovered series (Patterson, 1935), introduced the major modifications to the Patterson which are still in use today. He illustrated, with onedimensional Fourier series, the techniques of removing the Patterson origin peak, sharpening the overall function and also removing peaks due to atoms in special positions. Each one of these modifications can improve the interpretability of Pattersons, especially those of simple structures. Whereas the recommended extent of such modifications is controversial (Buerger, 1966), most studies which utilize Patterson functions do incorporate some of these techniques [see, for example, Jacobson et al. (1961), Braun et al. (1969) and Nordman (1980a)]. Since Patterson's original work, other workers have suggested that the Patterson function itself might be modified; Fourier inversion of the modified Patterson then provides a new and perhaps more tractable set of structure factors (McLachlan & Harker, 1951; Simonov, 1965; Raman, 1966; Corfield & Rosenstein, 1966). Karle & Hauptman (1964) suggested that an improved set of structure factors could be obtained from an originremoved Patterson modified such that it was everywhere nonnegative and that Patterson density values less than a bonding distance from the origin were set to zero. Nixon (1978) was successful in solving a structure which had previously resisted solution by using a set of structure factors which had been obtained from a Patterson in which the largest peaks had been attenuated.
The N origin peaks [see expression (2.3.1.3)] may be removed from the Patterson by using coefficients A Patterson function using these modified coefficients will retain all interatomic vectors. However, the observed structure factors must first be placed on an absolute scale (Wilson, 1942) in order to match the scatteringfactor term. In practice, Patterson origins can also be removed by using coefficientswhere is the average reflection intensity, usually calculated in several resolution shells. This formula has the advantage that the observed structure factors do not need to be on absolute scale.
Analogous to origin removal, the vector interactions due to atoms in known positions can also be removed from the Patterson function. Patterson showed that nonorigin Patterson peaks arising from known atoms 1 and 2 may be removed by using the expression where and are the positions of atoms 1 and 2 and and are their respective thermal correction factors. Using onedimensional Fourier series, Patterson illustrated how interactions due to known atoms can obscure other information.
Patterson also introduced a means by which the peaks in a Patterson function may be artificially sharpened. He considered the effect of thermal motion on the broadening of electrondensity peaks and consequently their Patterson peaks. He suggested that the coefficients could be corrected for thermal effects by simulating the atoms as point scatterers and proposed using a modified set of coefficients where , the average scattering factor per electron, is given by A common formulation for this type of sharpening expresses the atomic scattering factors at a given angle in terms of an overall isotropic thermal parameter B as The Patterson coefficients then become The normalized structure factors, E (see Chapter 2.2 ), which are used in crystallographic direct methods, are also a common source of sharpened Patterson coefficients . Although the centre positions and total contents of Patterson peaks are unaltered by sharpening, the resolution of individual peaks may be enhanced. The degree of sharpening can be controlled by adjusting the size of the assumed B factor; Lipson & Cochran (1966, pp. 165–170) analysed the effect of such a choice on Patterson peak shape.
All methods of sharpening Patterson coefficients aim at producing a point atomic representation of the unit cell. In this quest, the highresolution terms are enhanced (Fig. 2.3.1.1). Unfortunately, this procedure must also produce a serious Fourier truncation error which will be seen as large ripples about each Patterson peak (Gibbs, 1898). Consequently, various techniques have been used (mostly unsuccessfully) in an attempt to balance the advantages of sharpening with the disadvantages of truncation errors.
Schomaker and Shoemaker [unpublished; see Lipson & Cochran (1966, p. 168)] used a function in which s is the length of the scattering vector, to produce a Patterson synthesis which is less sensitive to errors in the loworder terms. Jacobson et al. (1961) used a similar function, which they rationalize as the sum of a scaled exponentially sharpened Patterson and a gradient Patterson function (the value of k was empirically chosen as ). This approach was subsequently further developed and generalized by Wunderlich (1965).
Interpretation of any Patterson requires some assumption, such as the existence of discrete atoms. A complete interpretation might also require an assumption of the number of atoms and may require other external information (e.g. bond lengths, bond angles, van der Waals separations, hydrogen bonding, positive density etc.). To what extent is the solution of a Patterson function unique? Clearly the greater is the supply of external information, the fewer will be the number of possible solutions. Other constraints on the significance of a Patterson include the error involved in measuring the observed coefficients and the resolution limit to which they have been observed. The larger the error, the larger the number of solutions. When the error on the amplitudes is infinite, it is only the other physical constraints, such as packing, which limit the structural solutions. Alternative solutions of the same Patterson are known as homometric structures.
During their investigation of the mineral bixbyite, Pauling & Shappell (1930) made the disturbing observation that there were two solutions to the structure, with different arrangements of atoms, which yielded the same set of calculated intensities. Specifically, atoms occupying equipoint set 24d in space group can be placed at either or without changing the calculated intensities. Yet the two structures were not chemically equivalent. These authors resolved the ambiguity by placing the oxygen atoms in question at positions which gave the most acceptable bonding distances with the rest of the structure.
Patterson interpreted the above ambiguity in terms of the series: the distance vector sets or Patterson functions of the two structures were the same since each yielded the same calculated intensities (Patterson, 1939). He defined such a pair of structures a homometric pair and called the degenerate vector set which they produced a homometric set. Patterson went on to investigate the likelihood of occurrence of homometric structures and, indeed, devoted a great deal of his time to this matter. He also developed algebraic formalisms for examining the occurrence of homometric pairs and multiplets in selected onedimensional sets of points, such as cyclotomic sets, and also sets of points along a line (Patterson, 1944). Some simple homometric pairs are illustrated in Fig. 2.3.1.2.
Drawing heavily from Patterson's inquiries into the structural uniqueness allowed by the diffraction data, Hosemann, Bagchi and others have given formal definitions of the different types of homometric structures (Hosemann & Bagchi, 1954). They suggested a classification divided into pseudohomometric structures and homomorphs, and used an integral equation representing a convolution operation to express their examples of finite homometric structures. Other workers have chosen various means for describing homometric structures [Buerger (1959, pp. 41–50), Menzer (1949), Bullough (1961, 1964), Hoppe (1962)].
Since a Patterson function is centrosymmetric, the Pattersons of a crystal structure and of its mirror image are identical. Thus the enantiomeric ambiguity present in noncentrosymmetric crystal structures cannot be overcome by using the Patterson alone and represents a special case of homometric structures. Assignment of the correct enantiomorph in a crystal structure analysis is generally not possible unless a recognizable fragment of known chirality emerges (e.g. Lamino acids in proteins, Driboses in nucleic acids, the known framework of steroids and other natural products, the righthanded twist of αhelices, the lefthanded twist of successive strands in a βsheet, the fold of a known protein subunit etc.) or anomalousscattering information is available and can be used to resolve the ambiguity (Bijvoet et al., 1951).
It is frequently necessary to select arbitrarily one enantiomorph over another in the early stages of a structure solution. Structurefactor phases calculated from a single heavy atom in space group P1, P2 or (for instance) will be centrosymmetric and both enantiomorphs will be present in Fourier calculations based on these phases. In other space groups (e.g. ), the selected heavy atom is likely to be near one of the planes containing the axes and thus produce a weaker `ghost' image of the opposite enantiomorph. The mixture of the two overlapped enantiomorphic solutions can cause interpretive difficulties. As the structure solution progresses, the `ghosts' are exorcized owing to the dominance of the chosen enantiomorph in the phases.
Patterson also defined a second, less well known, function (Patterson, 1949) as This function can be computed directly only for centrosymmetric structures. It can be calculated for noncentrosymmetric structures when the phase angles are known or assumed. It will represent the degree to which the known or assumed structure has a centre of symmetry at u. That is, the product of the density at and is large when integrated over all values x within the unit cell. Since atoms themselves have a centre of symmetry, the function will contain peaks at each atomic site roughly proportional in height to the square of the number of electrons in each atom plus peaks at the midpoint between atoms proportional in height to the product of the electrons in each atom. Although this function has not been found very useful in practice, it is useful for demonstrating the presence of weak enantiomorphic images in a given tentative structure determination.
References
International Tables for Crystallography (2005). Vol. A, SpaceGroup Symmetry, edited by Th. Hahn. Heidelberg: Springer.Beevers, C. A. & Robertson, J. M. (1950). Interpretation of the Patterson synthesis. Acta Cryst. 3, 164.
Bijvoet, J. M., Peerdeman, A. F. & van Bommel, A. J. (1951). Determination of the absolute configuration of optically active compounds by means of Xrays. Nature (London), 168, 271–272.
Blow, D. M. (1958). The structure of haemoglobin. VII. Determination of phase angles in the noncentrosymmetric [100] zone. Proc. R. Soc. London Ser. A, 247, 302–336.
Braun, P. B., Hornstra, J. & Leenhouts, J. I. (1969). Automated crystalstructure determination by Patterson search using a known part of the molecule. Philips Res. Rep. 24, 85–118.
Buerger, M. J. (1959). Vector Space and its Application in CrystalStructure Investigation. New York: John Wiley.
Buerger, M. J. (1966). Background for the use of imageseeking functions. Trans. Am. Crystallogr. Assoc. 2, 1–9.
Bullough, R. K. (1961). On homometric sets. I. Some general theorems. Acta Cryst. 14, 257–269.
Bullough, R. K. (1964). On homometric sets. II. Sets obtained by singular transformations. Acta Cryst. 17, 295–308.
Clastre, J. & Gay, R. (1950). La détermination des structures cristallines à partir du diagramme de Patterson. Compt. Rend. 230, 1876–1877.
Corfield, P. W. R. & Rosenstein, R. D. (1966). Maximum information from the minimum function. Trans. Am. Crystallogr. Assoc. 2, 17–28.
Garrido, J. (1950a). Sur la détermination des structures cristallines au moyen de la transformée de Patterson. Compt. Rend. 230, 1878–1879.
Gibbs, J. W. (1898). Remarks regarding Michelson's letter. Nature (London), 59, 200.
Hamilton, W. C. (1965). The crystal structure of orthorhombic acetamide. Acta Cryst. 18, 866–870.
Harker, D. (1936). The application of the threedimensional Patterson method and the crystal structures of proustite, Ag_{3}AsS_{3}, and pyrargyrite, Ag_{3}SbS_{3}. J. Chem. Phys. 4, 381–390.
Hoppe, W. (1962). `NahezuHomometrische Lösungen' und Faltmolekülmethode. Z. Kristallogr. 117, 249–258.
Hosemann, R. & Bagchi, S. N. (1954). On homometric structures. Acta Cryst. 7, 237–241.
Huber, R. (1965). Die automatisierte Faltmolekülmethode. Acta Cryst. 19, 353–356.
Hughes, E. W. (1940). The crystal structure of dicyandiamide. J. Am. Chem. Soc. 62, 1258–1267.
Jacobson, R. A., Wunderlich, J. A. & Lipscomb, W. N. (1961). The crystal and molecular structure of cellobiose. Acta Cryst. 14, 598–607.
Karle, J. & Hauptman, H. (1964). Positivity, point atoms, and Pattersons. Acta Cryst. 17, 392–396.
Ketelaar, J. A. A. & de Vries, T. A. (1939). The crystal structure of tetra phosphonitrile chloride, P_{4}N_{4}Cl_{8}. Recl Trav. Chim. 58, 1081–1099.
Kraut, J. (1961). The crystal structure of 2aminoethanol phosphate. Acta Cryst. 14, 1146–1152.
Lipson, H. & Cochran, W. (1966). The Determination of Crystal Structures. Ithaca: Cornell University Press.
McLachlan, D. Jr & Harker, D. (1951). Finding the signs of the F's from the shifted Patterson product. Proc. Natl Acad. Sci. USA, 37, 846–849.
Menzer, G. (1949). Über die mehrdeutigkeit der Kristallstrukturbestimmung. Z. Naturforsch. Teil A, 4, 11–21.
Mighell, A. D. & Jacobson, R. A. (1963). Analysis of threedimensional Patterson maps using vector verification. Acta Cryst. 16, 443–445.
Nixon, P. E. (1978). Overlapping Patterson peaks and direct methods: the structure of prostratin. Acta Cryst. A34, 450–453.
Nordman, C. E. (1980a). Vectorspace Patterson search and other storedfunction sampling procedures. In Computing in Crystallography, edited by R. Diamond, S. Ramaseshan & K. Venkatesan, pp. 5.01–5.13. Bangalore: Indian Academy of Sciences.
Nordman, C. E. & Nakatsu, K. (1963). Interpretation of the Patterson function of crystals containing a known molecular fragment. The structure of an Alstonia alkaloid. J. Am. Chem. Soc. 85, 353–354.
Patterson, A. L. (1934a). A Fourier series representation of the average distribution of scattering power in crystals. Phys. Rev. 45, 763.
Patterson, A. L. (1934b). A Fourier series method for the determination of the components of interatomic distances in crystals. Phys. Rev. 46, 372–376.
Patterson, A. L. (1935). A direct method for the determination of the components of interatomic distances in crystals. Z. Kristallogr. 90, 517–542.
Patterson, A. L. (1939). Homometric structures. Nature (London), 143, 939–940.
Patterson, A. L. (1944). Ambiguities in the Xray analysis of crystal structures. Phys. Rev. 65, 195–201.
Patterson, A. L. (1949). An alternative interpretation for vector maps. Acta Cryst. 2, 339–340.
Pauling, L. & Shappell, M. D. (1930). The crystal structure of bixbyite and the Cmodification of the sesquioxides. Z. Kristallogr. 75, 128–142.
Raman, S. (1966). Patterson functions and vector sets. Trans. Am. Crystallogr. Assoc. 2, 10–16.
Rossmann, M. G. (1960). The accurate determination of the position and shape of heavyatom replacement groups in proteins. Acta Cryst. 13, 221–226.
Rossmann, M. G. (1961a). The position of anomalous scatterers in protein crystals. Acta Cryst. 14, 383–388.
Rossmann, M. G. & Blow, D. M. (1962). The detection of subunits within the crystallographic asymmetric unit. Acta Cryst. 15, 24–31.
Shoemaker, D. P., Donohue, J., Schomaker, V. & Corey, R. B. (1950). The crystal structure of L_{8}threonine. J. Am. Chem. Soc. 72, 2328–2349.
Simonov, V. I. (1965). Calculation of the phases of the structure amplitudes by Fourier transformation of the sum, product and minimum functions. Proc. Indian Acad. Sci. A62, 213–223.
Simpson, P. G., Dobrott, R. D. & Lipscomb, W. N. (1965). The symmetry minimum function: high order image seeking functions in Xray crystallography. Acta Cryst. 18, 169–179.
Speakman, J. C. (1949). The crystal structures of the acid salts of some monobasic acids. Part I. Potassium hydrogen bisphenyl acetate. J. Chem. Soc. pp. 3357–3365.
Wilson, A. J. C. (1942). Determination of absolute from relative Xray intensity data. Nature (London), 150, 151–152.
Wrinch, D. M. (1939). The geometry of discrete vector maps. Philos. Mag. 27, 98–122.
Wunderlich, J. A. (1965). A new expression for sharpening Patterson functions. Acta Cryst. 19, 200–202.