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.1, pp. 200202
Section 2.1.5. Ideal probability density distributions^{a}School of Chemistry, Tel Aviv University, Tel Aviv 69 978, Israel, and ^{b}St John's College, Cambridge, England 
In applications of the centrallimit theorem, and its extensions, to intensity statistics the 's of equation (2.1.4.19) have the form (atomic scattering factor of the jth atom) times (a trigonometric expression characteristic of the space group and Wyckoff position; also known as the trigonometric structure factor). These trigonometric expressions for all the space groups, and general Wyckoff positions, are given in Tables A1.4.3.1 through A1.4.3.7 , and their first few even moments (fixedindex averaging) are given in Table 2.1.7.1. One cannot, of course, conclude that the magnitudes of the structure factor always have a normal distribution – even if the structure is homoatomic; one must look at each problem and see what components of the structure factor can be put in the form (2.1.4.19), deduce the m and to be used for each, and combine the components to obtain the asymptotic (large N, not large x) expression for the problem in question. Ordinarily the components are the real and the imaginary parts of the structure factor; the structure factor is purely real only if the structure is centrosymmetric, the spacegroup origin is chosen at a crystallographic centre and the atoms are nondispersive.
The ideal acentric distributions are obtained by applying the centrallimit theorem to the real and the imaginary parts of the structure factor, as given by equation (2.1.1.1). Consider first a crystal with no rotational symmetry (space group P1). The real part, A, of the structure factor is then given by where N is the number of atoms in the unit cell and is the phase angle of the jth atom. The centrallimit theorem then states that A tends to be normally distributed about its mean value with variance equal to its meansquare deviation from its mean. Under the assumption that the phase angles are uniformly distributed on the 0–2π range, the mean value of each cosine is zero, so that its variance is Under the same assumption, the mean value of each is onehalf, so that the variance becomes where Σ is the sum of the squares of the atomic scattering factors [cf. equation (2.1.2.4)]. The asymptotic form of the distribution of A is therefore given by A similar calculation, with sines instead of cosines, gives an analogous distribution for the imaginary part B, so that the joint probability of the real and imaginary parts of F is Ordinarily, however, we are more interested in the distribution of the magnitude, , of the structure factor than in the distribution of A and B. Using polar coordinates in equation (2.1.5.5) and integrating over the angle φ gives It is usually convenient, in structurefactor and intensity statistics, to express the results in terms of the normalized structure factor E and its magnitude . If has been put on an absolute scale (see Section 2.2.4.3 ), we have so that is the normalizedstructurefactor version of (2.1.5.6).
Distributions resulting from noncentrosymmetric crystals are known as acentric distributions; those arising from centrosymmetric crystals are known as centric. These adjectives are used to describe distributions, not crystal symmetry.
When a nondispersive crystal is centrosymmetric, and the spacegroup origin is chosen at a crystallographic centre of symmetry, the imaginary part B of its structure amplitude is zero. In the simplest case, space group , the contribution of the jth atom plus its centrosymmetric counterpart is . The calculation of goes through as before, with allowance for the fact that there are pairs instead of N independent atoms, giving or equivalently or
Additional crystallographic symmetry elements do not produce any essential alterations in the ideal centric or acentric distribution; their main effect is to replace the parameter Σ by a `distribution parameter', called S by Wilson (1950) and Rogers (1950), in certain groups of reflections. In addition, in noncentrosymmetric space groups, the distribution of certain groups of reflections becomes centric, though the general reflections remain acentric. The changes are summarized in Tables 2.1.3.1 and 2.1.3.2. The values of S are integers for lattice centring, glide planes and those screw axes that produce absences, and approximate integers for rotation axes and mirror planes; the modulations of the average intensity in reciprocal space outlined in Section 2.1.3.2 apply.
It should be noted that if intensities are normalized to the average of the group to which they belong, rather than to the general average, the distributions given in equations (2.1.5.8) and (2.1.5.11) are not affected.
The distributions just derived are asymptotic, as they are limiting values for large N. They are the only ideal distributions, in this sense, when there is only strict crystallographic symmetry and no dispersion. However, other ideal (asymptotic) distributions arise when there is noncrystallographic symmetry, or if there is dispersion. The subcentric distribution, where is a modified Bessel function of the first kind and k is the ratio of the scattering from the centrosymmetric part to the total scattering, arises when a noncentrosymmetric crystal contains centrosymmetric parts or when dispersion introduces effective noncentrosymmetry into the scattering from a centrosymmetric crystal (Srinivasan & Parthasarathy, 1976, ch. III; Wilson, 1980a,b; Shmueli & Wilson, 1983). The bicentric distribution arises, for example, when the `asymmetric unit in a centrosymmetric crystal is a centrosymmetric molecule' (Lipson & Woolfson, 1952); is a modified Bessel function of the second kind. There are higher hypercentric, hyperparallel and sesquicentric analogues (Wilson, 1952; Rogers & Wilson, 1953; Wilson, 1956). The ideal subcentric and bicentric distributions are expressed in terms of known functions, but the higher hypercentric and the sesquicentric distributions have so far been studied only through their moments and integral representations. Certain hypersymmetric distributions can be expressed in terms of Meijer's G functions (Wilson, 1987b).
When only the intrinsic probability distributions are being considered, it does not greatly matter whether the variable chosen is the intensity of reflection (I), or its positive square root, the modulus of the structure factor , since both are necessarily real and nonnegative. In an obvious notation, the relation between the intensity distribution and the structurefactor distribution is or Statistical fluctuations in counting rates, however, introduce a small but finite probability of negative observed intensities (Wilson, 1978a, 1980a) and thus of imaginary structure factors. This practical complication is treated in IT C (2004, Parts 7 and 8 ).
Both the ideal centric and acentric distributions are simple members of the family of gamma distributions, defined by where n is a parameter, not necessarily integral, and is the gamma function. Thus the ideal acentric intensity distribution is and the ideal centric intensity distribution is The properties of gamma distributions and of the related beta distributions, summarized in Table 2.1.5.1, are used in Section 2.1.6 to derive the probability density functions of sums and of ratios of intensities drawn from one of the ideal distributions.

The integral of the probability density function from the lower end of its range up to an arbitrary value x is called the cumulative probability distribution, or simply the distribution function, , of x. It can always be written if the lower end of its range is not actually one takes as identically zero between and the lower end of its range. For the distribution of A [equation (2.1.5.4) or (2.1.5.9)] the lower limit is in fact ; for the distribution of , , I and the lower end of the range is zero. In such cases, equation (2.1.5.21) becomes In crystallographic applications the cumulative distribution is usually denoted by , rather than by the capital letter corresponding to the probability density function designation. The cumulative forms of the ideal acentric and centric distributions (Howells et al., 1950) have found many applications. For the acentric distribution of [equation (2.1.5.8)] the integration is readily carried out: The integral for the centric distribution of [equation (2.1.5.11)] cannot be expressed in terms of elementary functions, but the integral required has so many important applications in statistics that it has been given a special name and symbol, the error function erf(x), defined by For the centric distribution, then The error function is extensively tabulated [see e.g. Abramowitz & Stegun (1972), pp. 310–311, and a closely related function on pp. 966–973].
The foregoing and following structurefactor and intensity statistics, presented in this chapter, are centred around the unambiguous determination of crystal symmetry. These statistics, especially the ideal ones, are routinely applied to Xray diffraction data for completing the symmetry information where needed. On the other hand, the symmetry of very fine grained, imperfectly crystalline material or individual nanostructures can be experimentally determined by the convergentbeam electron diffraction technique using an electron microscope with a nanoscale probe, as described and illustrated in Section 2.5.3 of this volume.
References
International Tables for Crystallography (2004). Vol. C, Mathematical, Physical and Chemical Tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.Abramowitz, M. & Stegun, I. A. (1972). Handbook of Mathematical Functions. New York: Dover.
Howells, E. R., Phillips, D. C. & Rogers, D. (1950). The probability distribution of Xray intensities. II. Experimental investigation and the Xray detection of centers of symmetry. Acta Cryst. 3, 210–214.
Kendall, M. & Stuart, A. (1977). The Advanced Theory of Statistics, Vol. 1, 4th ed. London: Griffin.
Lipson, H. & Woolfson, M. M. (1952). An extension of the use of intensity statistics. Acta Cryst. 5, 680–682.
Rogers, D. (1950). The probability distribution of Xray intensities. IV. New methods of determining crystal classes and space groups. Acta Cryst. 3, 455–464.
Rogers, D. & Wilson, A. J. C. (1953). The probability distribution of Xray intensities. V. A note on some hypersymmetric distributions. Acta Cryst. 6, 439–449.
Shmueli, U. & Wilson, A. J. C. (1983). Generalized intensity statistics: the subcentric distribution and effects of dispersion. Acta Cryst. A39, 225–233.
Srinivasan, R. & Parthasarathy, S. (1976). Some Statistical Applications in Xray Crystallography. Oxford: Pergamon Press.
Wilson, A. J. C. (1950). The probability distribution of Xray intensities. III. Effects of symmetry elements on zones and rows. Acta Cryst. 3, 258–261.
Wilson, A. J. C. (1952). Hypercentric and hyperparallel distributions of Xray intensities. Research (London), 5, 588–589.
Wilson, A. J. C. (1956). The probability distribution of Xray intensities. VII. Some sesquicentric distributions. Acta Cryst. 9, 143–144.
Wilson, A. J. C. (1978a). On the probability of measuring the intensity of a reflection as negative. Acta Cryst. A34, 474–475.
Wilson, A. J. C. (1980a). Relationship between `observed' and `true' intensity: effects of various counting modes. Acta Cryst. A36, 929–936.
Wilson, A. J. C. (1980b). Effect of dispersion on the probability distribution of Xray reflections. Acta Cryst. A36, 945–946.
Wilson, A. J. C. (1987b). Functional form of the ideal hypersymmetric distributions of structure factors. Acta Cryst. A43, 554–556.