International
Tables for
Crystallography
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2010). Vol. B, ch. 2.1, pp. 195-214   | 1 | 2 |
https://doi.org/10.1107/97809553602060000763

Chapter 2.1. Statistical properties of the weighted reciprocal lattice

U. Shmuelia* and A. J. C. Wilsonb

aSchool of Chemistry, Tel Aviv University, Tel Aviv 69 978, Israel, and bSt John's College, Cambridge, England
Correspondence e-mail:  ushmueli@post.tau.ac.il

This chapter has two purposes: (i) it gives an introduction to the principles of probability, which plays an important role in most methods of structure determination, and (ii) it describes the application of such techniques to the first phase of structure determination – the resolution of space-group ambiguities, where present. The introduction in Section 2.1.1[link] justifies briefly the application of statistical methods to the structure-factor function. In Sections 2.1.2[link] and 2.1.3[link], a discussion is presented of the average intensity of general reflections and of zones and rows, with particular attention to the relation of these averages to symmetry elements present in the crystal. This is followed in Section 2.1.4[link] by mathematical preliminaries to the basics of the calculus of probabilities and the central-limit theorem is introduced. In Section 2.1.5[link], probability density functions (p.d.f.'s) are derived which allow one to determine whether or not the crystal is centrosymmetric. These p.d.f.'s are widely applicable to structures with a large number of not too dissimilar atoms in the asymmetric unit. Such functions for structures containing dispersive scatterers and noncrystallographic centres of symmetry are also presented, and distributions of sums, averages and ratios of intensities are discussed in Section 2.1.6[link]. All these p.d.f.'s are based on the central-limit theorem and are termed ideal p.d.f.'s. However, these ideal p.d.f.'s are no longer applicable when an outstandingly heavy atom is present in the asymmetric unit. Two approaches are discussed in this chapter which may resolve this difficulty: (i) the correction-factor approach (Section 2.1.7[link]), well known from classical probability as Gram–Charlier or Edgeworth p.d.f.'s, and (ii) the Fourier method (Section 2.1.8[link]), only recently introduced to crystallography. The Gram–Charlier p.d.f.'s depend on even moments of the magnitude of the structure factor. General expressions for such moments of all orders are given for space groups of low symmetry, and the first four even moments are tabulated for all 230 space groups. The Fourier p.d.f.'s depend on their characteristic functions, or their Fourier transforms. Expressions for the atomic contributions to the characteristic function are tabulated for space groups up to and including the cubic space group [Fd\bar{3}]. Higher cubic space groups can be satisfactorily treated with the correction-factor p.d.f.'s. A comparison of these two non-ideal methods indicates that in the instances in which departures from the central-limit theorem predictions are large, the Fourier method is definitely superior to the correction-factor method.

2.1.1. Introduction

| top | pdf |

The structure factor of the hkl reflection is given by[F(hkl) = \textstyle\sum\limits_{j = 1}^{N}f_{j}\exp[2\pi i(hx_{j}+ky_{j}+lz_{j})],\eqno(2.1.1.1)]where [f_{j}] is the atomic scattering factor [complex if there is appreciable dispersion; see Chapter 1.2[link] and IT C (2004, Section 4.2.6[link] [link])], [x_{j}y_{j}z_{j}] are the fractional coordinates of the jth atom and N is the number of atoms in the unit cell. The present chapter is concerned with the statistical properties of the structure factor F and the intensity [I = FF^{*}], such as their average values, variances, higher moments and their probability density distributions.

Equation (2.1.1.1)[link] expresses F as a function of two conceptually different sets of variables: hkl taking on integral values in reciprocal space and xyz in general having non-integral values in direct space, although the special positions tabulated for each space group in IT A (2005)[link] may include the integers 0 and 1. In special positions, the non-integers often include rational fractions, but in general positions they are in principle irrational. Although hkl and xyz appear to be symmetrical variables in (2.1.1.1)[link], these limitations on their values mean that one can consider two different sets of statistical properties. In the first we seek, for example, the average intensity of the hkl reflection (indices fixed) as the positional parameters of the N atoms are distributed with equal probability over the continuous range 0–1. In the second, we seek, for example, the average intensity of the observable reflections (or of a subgroup of them having about the same value of [\sin \theta]) with the values of xyz held constant at the values they have, or are postulated to have, in a crystal structure. Other examples are obtained by substituting the words `probability density' for `average intensity'. For brevity, we may call the statistics resulting from the first process fixed-index (continuously variable parameters being understood), and those resulting from the second process fixed-parameter (integral indices being understood). Theory based on the first process is (comparatively) easy; theory based on the second hardly exists, although there is a good deal of theory concerning the conditions under which the two processes will lead to the same result (Hauptman & Karle, 1953[link]; Giacovazzo, 1977[link], 1980[link]). Mathematically, of course, the condition is that the phase angle[\vartheta = 2\pi (hx+ky+lz) \eqno(2.1.1.2)]should be distributed with uniform probability over the range 0–2π, whichever set of variables is regarded as fixed, but it is not clear when this distribution can be expected in practice for fixed-parameter averaging. The usual conclusion is that the uniform distribution will be realized if there are enough atoms, if the atomic coordinates do not approximate to rational fractions, if there are enough reflections and if stereochemical effects are negligible (Shmueli et al., 1984[link]).

Obviously, the second process (fixed parameters, varying integral indices) corresponds to the observable reality, and various approximations to it have been attempted, in preference to assuming its equivalence with the first. For example, a third (approximate) method of averaging has been used (Wilson, 1949[link], 1981[link]): xyz are held fixed and hkl are treated as continuous variables.

2.1.2. The average intensity of general reflections

| top | pdf |

2.1.2.1. Mathematical background

| top | pdf |

The process may be illustrated by evaluating, or attempting to evaluate, the average intensity of reflection by the three processes. The intensity of reflection is given by multiplying equation (2.1.1.1)[link] by its complex conjugate: [\eqalignno{I &= FF^{*} &(2.1.2.1) \cr &= \textstyle\sum\limits_{j, \, k}f_{j}\,f_{k}^{*}\exp\{2\pi i[h(x_{j}-x_{k})+\ldots]\} &(2.1.2.2) \cr &= \Sigma + \textstyle\sum\limits_{j \neq k}f_{j}\,f_{k}^{*}\exp\{2\pi i[h(x_{j}-x_{k})+\ldots]\}, \qquad &(2.1.2.3)}%2.1.2.3]where [\Sigma = \textstyle\sum\limits_{j}f_{j}\,f_{j}^{*} \eqno(2.1.2.4)]is the sum of the squares of the moduli of the atomic scattering factors. Wilson (1942[link]) argued, without detailed calculation, that the average value of the exponential term would be zero and hence that [\langle I \rangle = \Sigma. \eqno(2.1.2.5)]Averaging equation (2.1.2.3)[link] for hkl fixed, xyz ranging uniformly over the unit cell – the first process described above – gives this result identically, without complication or approximation. Ordinarily the second process cannot be carried out. We can, however, postulate a special case in which it is possible. We take a homoatomic structure and before averaging we correct the f's for temperature effects and the fall-off with [\sin \theta], so that [ff^{*}] is the same for all the atoms and is independent of hkl. If the range of hkl over which the expression for I has to be averaged is taken as a parallelepiped in reciprocal space with h ranging from −H to +H, k from −K to +K, l from −L to +L, equation (2.1.2.2)[link] can be factorized into the product of the sums of three geometrical progressions. Algebraic manipulation then easily leads to [\eqalignno{\langle I \rangle & = ff^{*} \sum _{j}\sum _{k}{{\sin \pi N_{H}(x_{j}-x_{k})}\over{N_{H}\sin \pi (x_{j}-x_{k})}} \cr &\quad\times {{\sin \pi N_{K}(y_{j}-y_{k})}\over{N_{K}\sin \pi (y_{j}-y_{k})}} {{\sin \pi N_{L}(z_{j}-z_{k})}\over{N_{L}\sin \pi (z_{j}-z_{k})}}, \cr & &(2.1.2.6)}]where [N_{H} = 2H+1], [N_{K} = 2K+1] and [N_{L} = 2L+1]. The terms with j = k give Σ, but the remaining terms are not zero. Because of the periodic nature of the trigonometric terms, the effective coordinate differences are never greater than 0.5 and in a structure of any complexity there will be many much less than 0.5. For HKL = 000, in fact, [\langle I \rangle] becomes the square of the modulus of the sum of the atomic scattering factors, [\langle I \rangle = \Phi\Phi^{*}, \eqno(2.1.2.7)]where [\Phi = \textstyle\sum\limits_{j = 1}^{N}f_{j}, \eqno(2.1.2.8)]and not the sum of the squares of their moduli; for larger HKL, [\langle I \rangle] rapidly decreases to Σ and then oscillates about that value. Wilson (1949[link], especially Section 2.1.1) suggested that the regions of averaging should be chosen so that at least one index of every reflection is [\geq 2] if [\langle I \rangle] is to be identified with Σ, and this has proven to be a useful rule-of-thumb.

The third process of averaging replaces the sum over integral values of the indices by an integration over continuous values, the appropriate values of the limits in this example being [-(H+1/2)] to [+(H+1/2)]. The effect is to replace the sines in the denominators, but not in the numerators, of equation (2.1.2.6)[link] by their arguments, and this is equivalent to the approximation [\sin x \simeq x] in the denominators only. This is a good approximation for atoms close together in the structure and thus giving the largest terms in the sums in equation (2.1.2.6)[link], and gives the correct sign and order of magnitude even for x having its maximum value of π/2.

2.1.2.2. Physical background

| top | pdf |

The preceding section has used mathematical arguments. From a physical point of view, the radiation diffracted by atoms that are resolved will interfere destructively, so that the resulting intensity will be the sum of the intensities diffracted by individual atoms, whereas that from completely unresolved atoms will interfere constructively, so that amplitudes rather than intensities add. In intermediate cases there will be partial constructive interference. Resolution in accordance with the Rayleigh (1879[link]) criterion requires that [s = (2\sin \theta)/\lambda] should be greater than half the reciprocal of the minimum interatomic distance in the crystal (Wilson, 1979[link]); full resolution requires a substantial multiple of this. This criterion is essentially equivalent to that proposed from the study of a special case of the second process in the preceding section.

2.1.2.3. An approximation for organic compounds

| top | pdf |

In organic compounds there are very many interatomic distances of about 1.5 or 1.4 Å. Adoption of the preceding criterion would mean that the inner portion of the region of reciprocal space accessible by the use of copper Kα radiation is not within the sphere of intensity statistics based on fixed-index (first process) averaging. No substantial results are available for fixed-parameter (second process) averaging, and very few from the approximation to it (third process).

To the extent to which the third process is acceptable, an approximation to the variation of [\langle I \rangle] with [\sin\theta] is obtainable. The exponent in equation (2.1.2.2)[link] can be written as [2\pi i sr_{jk} \cos\psi, \eqno(2.1.2.9)]where s is the radial distance in reciprocal space, [r_{jk}] is the distance from the jth to the kth atom and [\psi] is the angle between the vectors s and r. Averaging over a sphere of radius s, with [\psi] treated as the colatitude, gives [\langle I \rangle = \sum_{j}\sum_{k}f_{j}\,f_{k}{{\sin 2\pi sr_{jk}}\over {2\pi sr_{jk}}}. \eqno(2.1.2.10)]This is the familiar Debye expression. It has the correct limits for s zero and s large, and is in accord with the argument from resolution.

2.1.2.4. Effect of centring

| top | pdf |

In the preceding discussion there has been a tacit assumption that the lattice is primitive. A centred crystal can always be referred to a primitive lattice and if this is done no change is required. If the centred lattice is retained, many reflections are identically zero and the intensity of the nonzero reflections is enhanced by a factor of two (I and C lattices) or four (F lattice), so that the average intensity of all the reflections, zero and nonzero taken together, is unchanged.

Other symmetry elements affect only zones and rows of reflections, and so do not affect the general average when the total number of reflections is large. Their effect on zones and rows is discussed in Section 2.1.3[link].

2.1.3. The average intensity of zones and rows

| top | pdf |

2.1.3.1. Symmetry elements producing systematic absences

| top | pdf |

Symmetry elements can be divided into two types: those that cause systematic absences and those that do not. Those producing systematic absences (glide planes and screw axes) produce at the same time groups of reflections (confined to zones and rows in reciprocal space, respectively) with an average intensity an integral1 multiple of the general average. The effects for single symmetry elements of this type are given in Table 2.1.3.1[link] for the general reflections hkl and separately for any zones and rows that are affected. The `average multipliers' are given in the column headed [\langle I \rangle/\Sigma]; `distribution' and `distribution parameters' are treated in Section 2.1.5[link]. As for the centring, the fraction of reflections missing and the integer multiplying the average are related in such a way that the overall intensity is unchanged. The mechanism for compensation for the reflections with enhanced intensity is obvious.

Table 2.1.3.1| top | pdf |
Intensity-distribution effects of symmetry elements and centred lattices causing systematic absences

Abbreviations and orientation of axes: A = acentric distribution, C = centric distribution, Z = systematically zero, S = distribution parameter, <I> = average intensity. Axes are parallel to c, planes are perpendicular to c.

ElementReflectionsDistribution[S/\Sigma][\langle I\rangle/\Sigma]
[2_{1}] hkl A 1 1
hk0 C 1 1
00l [(Z + A)/2] 1 2
[3_{1}, 3_{2}] hkl A 1 1
hk0 A 1 1
00l [(2Z + A)/3] 1 3
[4_{1}, 4_{3}] hkl A 1 1
hk0 C 1 1
00l [(3Z + A)/4] 1 4
[4_{2}] hkl A 1 1
hk0 C 1 1
00l [(Z + A)/2] 2 4
[6_{1}, 6_{5}] hkl A 1 1
hk0 C 1 1
00l [(5Z + A)/6] 1 6
[6_{2}, 6_{4}] hkl A 1 1
hk0 C 1 1
00l [(2Z + A)/3] 2 6
[6_{3}] hkl A 1 1
hk0 C 1 1
00l [(Z + A)/2] 3 6
a hkl A 1 1
hk0 [(Z + A)/2] 1 2
00l C 1 1
0k0 A 2 2
C, I All [(Z + A)/2] 1 2
F All [(3Z + A)/2] 1 4

2.1.3.2. Symmetry elements not producing systematic absences

| top | pdf |

Certain symmetry elements not producing absences (mirror planes and rotation axes) cause equivalent atoms to coincide in a plane or a line projection and hence produce a zone or row in reciprocal space for which the average intensity is an integral multiple of the general average (Wilson, 1950[link]); the effects of single such symmetry elements are given in Table 2.1.3.2[link]. There is, however, no obvious mechanism for compensation for this enhancement. When reflections are few this may be an important matter in assigning an approximate absolute scale by comparing observed and calculated intensities. Wilson (1964[link]), Nigam (1972[link]) and Nigam & Wilson (1980[link]), noting that in such cases the finite size of atoms results in forbidden ranges of positional parameters, have shown that there is a diminution of the intensity of layers (rows) in the immediate neighbourhood of the enhanced zones (rows), just sufficient to compensate for the enhancement. In forming general averages, therefore, reflections from enhanced zones or rows should be included at their full intensity, not divided by the multiplier; the matter is discussed in more detail by Wilson (1987[link]a). It should be noted, however, that organic structures containing molecules related by rotation axes are rare, and such structures related by mirror planes are even rarer (Wilson, 1993[link]).

Table 2.1.3.2| top | pdf |
Intensity-distribution effects of symmetry elements not causing systematic absences

Abbreviations and orientation of axes: A = acentric distribution, C = centric distribution, S = distribution parameter, <I> = average intensity. Axes are parallel to c, planes are perpendicular to c.

ElementReflectionsDistribution[S/\Sigma = \langle I\rangle / \Sigma]
1 All A 1
[\bar{1}] All C 1
2 hkl A 1
hk0 C 1
00l A 2
[\bar{2} = m] hkl A 1
hk0 A 2
00l C 1
3 hkl A 1
hk0 A 1
00l A 3
[\bar{3}] hkl C 1
hk0 C 1
00l C 3
4 hkl A 1
hk0 C 1
00l A 4
[\bar{4}] hkl A 1
hk0 C 1
00l C 2
6 hkl A 1
hk0 C 1
00l A 6
[\bar{6} = 3/m] hkl A 1
hk0 A 2
00l C 3

2.1.3.3. More than one symmetry element

| top | pdf |

Further alterations of the intensities occur if two or more such symmetry elements are present in the space group. The effects were treated in detail by Rogers (1950[link]), who used them to construct a table for the determination of space groups by supplementing the usual knowledge of Laue group with statistical information. Only two pairs of space groups, the orthorhombic I222 and [I2_{1}2_{1}2_{1}], and their cubic supergroups I23 and [I2_{1}3_{1}], remained unresolved. Examination of this table shows that what statistical information does is to resolve the Laue group into point groups; the further resolution into space groups is equivalent to the use of Table 3.1.4.1[link] in IT A (2005)[link]. The statistical consequences of each point group, as given by Rogers, are reproduced in Table 2.1.3.3[link].

Table 2.1.3.3| top | pdf |
Average multiples for the 32 point groups (modified from Rogers, 1950[link])

The multiple gives S/Σ for the row and zone corresponding to the principal axis of the point-group symbol; those for the secondary and tertiary axes are given when the symbol contains such axes.

Point groupPrincipalSecondaryTertiary
RowZoneRowZoneRowZone
1 1 1        
[\bar{1}] 1 1        
2 2 1        
m 1 2        
[2/m] 2 2        
222 2 1 2 1 2 1
mm2 2 2 2 2 4 1
mmm 4 2 4 2 4 2
4 4 1        
[\bar{4}] 2 1        
[4/m] 4 2        
422 4 1 2 1 2 1
4mm 8 1 2 2 2 2
[\bar{4}2m] 4 1 2 1 2 2
[4/mmm] 8 2 4 2 4 2
3 3 1        
[\bar{3}] 3 1        
321 3 1 2 1 1 1
3m1 6 1 1 2 2  
31m 6 1 2 2 2 1
6 6 1        
[\bar{6}] 3 2        
[6/m] 6 2        
622 6 1 2 1 2 1
6mm 12 1 2 2 2 2
[\bar{6}m2] 6 2 2 2 4 1
[6/mmm] 12 2 4 2 4 2
231 2 1 3 1 1 1
[m\bar{3}1] 4 2 3 1 1 1
432 4 1 3 1 2 1
[\bar{4}3m] 4 1 6 1 2 2
[m\bar{3}m] 8 2 6 2 4 2

Note. The pairs of point groups, 1 and [\bar{1}] and 3 and [\bar{3}], not distinguished by average multiples, may be distinguished by their centric and acentric probability density functions.
The entry for the principal zone for the point group 422 was given incorrectly as 2 in the first edition of this volume.

2.1.4. Probability density distributions – mathematical preliminaries

| top | pdf |

For the purpose of this chapter, `ideal' probability distributions or probability density functions are the asymptotic forms obtained by the use of the central-limit theorem when the number of atoms in the unit cell, N, is sufficiently large. In order to derive them it is necessary to outline the properties of characteristic functions and to state alternative conditions for the validity of the central-limit theorem; the distributions themselves are derived in Section 2.1.5[link].

2.1.4.1. Characteristic functions

| top | pdf |

The average value of [\exp(itx)] is very important in probability theory; it is called the characteristic function of the distribution [f(x)] and is denoted by [C_{x}(t)] or, when no confusion can arise, by [C(t)]. It exists for all legitimate distributions, whether discrete or continuous. In the continuous case it is given by [C(t) = \textstyle\int\limits_{-\infty}^{\infty}\exp(itx)f(x)\,{\rm d}x, \eqno(2.1.4.1)]and is thus the Fourier transform of [f(x)]. In many cases it can be obtained from known integrals. For example, for the Cauchy distribution, [\eqalignno{C(t)& = {{a}\over{\pi}}\int\limits_{-\infty}^{\infty}{{\exp(itx)}\over{a^{2}+x^{2}}} \,{\rm d}x &(2.1.4.2) \cr & = \exp(-a|t|), &(2.1.4.3)}%2.1.4.3]and for the normal distribution, [\eqalignno{C(t)& = (2\pi\sigma^{2})^{-1/2}\int\limits_{-\infty}^{\infty}\exp\left(-{(x-m)^{2}}\over{2\sigma^{2}}\right)\exp(itx)\,{\rm d}x \cr &&(2.1.4.4) \cr & = \exp\left(imt-{{\sigma^{2}t^{2}}\over{2}}\right). &(2.1.4.5)}%2.1.4.5]Since the characteristic function is the Fourier transform of the distribution function, the converse is true, and if the characteristic function is known the probability distribution function can be obtained by the use of Fourier inversion theorem, [f(x) = (1/2\pi)\textstyle\int\limits_{-\infty}^{\infty}\exp(-itx)C(t)\,{\rm d}t. \eqno(2.1.4.6)]An alternative approach to the derivation of the distribution from a known characteristic function will be discussed below.

The most important property of characteristic functions in crystallography is the following: if x and y are independent random variables with characteristic functions [C_{x}(t)] and [C_{y}(t)], the characteristic function of their sum [z = x+y \eqno(2.1.4.7)]is the product [C_{z}(t) = C_{x}(t)C_{y}(t). \eqno(2.1.4.8)]Obviously this can be extended to any number of independent random variables.

When the moments exist, the characteristic function can be expanded in a power series in which the kth term is [m_{k}(it){}^{k}/k!]. If the power series [\exp(itx) = 1+itx+{{(it)^{2}x^{2}}\over{2!}}+{{(it)^{3}x^{3}}\over{3!}}+\ldots \eqno(2.1.4.9)]is substituted in equation (2.1.4.1)[link], one obtains [C(t) = 1+itm_{1}'+{{(it)^{2}m_{2}'}\over{2!}}+{{(it)^{3}m_{3}'}\over{3!}}+ \ldots. \eqno(2.1.4.10)]The moments are written with primes in order to indicate that equation (2.1.4.10)[link] is valid for moments about an arbitrary origin as well as for moments about the mean. If the random variable is transformed by a change of origin and scale, say [y = {{x-a}\over{b}}, \eqno(2.1.4.11)]the characteristic function for y becomes [C_{y}(t) = b\exp(-iat/b)C_{x}(t). \eqno(2.1.4.12)]

2.1.4.2. The cumulant-generating function

| top | pdf |

A function that is often more useful than the characteristic function is its logarithm, the cumulant-generating function: [K(t) = \log C(t) = k_{1}+{{k_{2}(it)^{2}}\over{2!}}+{{k_{3}(it)^{3}}\over{3!}}+ \ldots, \eqno(2.1.4.13)]where the k's are called the cumulants and may be regarded as being defined by the equation. They can be evaluated in terms of the moments by combining the series (2.1.4.10)[link] for [C(t)] with the ordinary series for the logarithm and equating the coefficients of [t^{r}]. In most cases the process as described is tedious, but it can be shortened by use of a general method [Stuart & Ord (1994[link]), Section 3.14, pp. 87–88; Exercise 3.19, p. 119]. Obviously, the cumulants exist only if the moments exist. The first few relations are [\eqalignno{k_{0} & = 0 \cr k_{1} & = m_{1}' \cr k_{2} & = m_{2} = m_{2}' - (m_{1}')^{2} &(2.1.4.14)\cr k_{3} & = m_{3} = m_{3}' - 3m_{2}'m_{1}' + 2(m_{1}')^{2} \cr k_{4} & = m_{4} - 3(m_{2})^{2} \cr & = m_{4}'-m_{3}'m_{1}'-3(m_{2}')^{2}+12m_{2}'(m_{1})^{2}-6(m_{1}')^{4}. }]Such expressions and their converses up to [k_{10}] are given by Stuart & Ord (1994[link], pp. 88–91). Since all the cumulants except [k_{1}] can be expressed in terms of the central moments only (i.e., those unprimed), only [k_{1}] is changed by a change of the origin. Because of this property, they are sometimes called the semi-invariants (or seminvariants) of the distribution. Since addition of random variables is equivalent to the multiplication of their characteristic functions [equation (2.1.4.8)[link]] and multiplication of functions is equivalent to the addition of their logarithms, each cumulant of the distribution of the sum of a number of random variables is equal to the sum of the cumulants of the distribution functions of the individual variables – hence the name cumulants. Although the cumulants (except [k_{1}]) are independent of a change of origin, they are not independent of a change of scale. As for the moments, a change of scale simply multiplies them by a power of the scale factor; if [y = x/b] [(k_{y})_{r} = (k_{x})_{r}/b^{r}. \eqno(2.1.4.15)]The cumulants of the normal distribution are particularly simple. From equation (2.1.4.5)[link], the cumulant-generating function of a normal distribution is [\eqalignno{K(t) & = imt-{{\sigma^{2}t^{2}}/{2}} &(2.1.4.16) \cr k_{1} & = m &(2.1.4.17) \cr k_{2} & = \sigma^{2}, &(2.1.4.18)}%2.1.4.18]all cumulants with [r \,\gt\, 2] are identically zero.

2.1.4.3. The central-limit theorem

| top | pdf |

A simple form of this important theorem can be stated as follows:

If [x_{1}, x_{2}, \ldots ,x_{n}] are independent and identically distributed random variables, each of them having the same mean m and variance [\sigma^{2}], then the sum [S_{n}=\textstyle\sum\limits_{j=1}^{n}x_{j} \eqno{(2.1.4.19)}]tends to be normally distributed – independently of the distribution(s) of the individual random variables – with mean [nm] and variance [n\sigma^{2}], provided n is sufficiently large.

In order to prove this theorem, let us define a standardized random variable corresponding to the sum [S_n], i.e., such that its mean is zero and its variance is unity:[\hat{S}_n={{S_n-nm}\over{\sigma \sqrt{n}}}={{\textstyle\sum_{j=1}^n(x_j-m )}\over{\sigma \sqrt{n}}}\equiv \sum\limits_{j=1}^n{{W_j}\over{\sqrt{n}}}, \eqno(2.1.4.20)]where [W_j=(x_j-m)/\sigma] is a standardized single random variable. The characteristic function of [\hat{S}_n] is therefore given by [\eqalignno{C_n(\hat{S}_n,t) &=\langle \exp (it\hat{S}_n)\rangle &\cr &=\left\langle \exp \left[ it\sum_{j=1}^n{{W_j}\over{\sqrt{n}}}\right] \right\rangle &(2.1.4.21)\cr &=\prod_{j=1}^n\left\langle \exp \left[ it{{W_j}\over{\sqrt{n}}}\right] \right\rangle &(2.1.4.22)\cr &=\left\{ \left\langle \exp \left[ it{{W_{_1}}\over{\sqrt{n}}}\right] \right\rangle \right\} ^n, &(2.1.4.23)\cr}%fd2.1.2.43]where the brackets [\langle\ \rangle] denote the operation of averaging with respect to the appropriate probability density function (p.d.f.) [cf. equation (2.1.4.1)[link]]. Equation (2.1.4.22)[link] follows from equation (2.1.4.21)[link] by the assumption of independence, while the assumption of identically distributed variables leads to the identity of the characteristic functions of the individual variables – as seen in equation (2.1.4.23)[link].

On the assumption that moments of all the orders exist – a most plausible assumption in situations usually encountered in structure-factor statistics – we can now expand the characteristic function of a single variable in a power series [cf. equation (2.1.4.10)[link]]: [\eqalignno{\left\langle \exp \left[ it{{W_{_1}}\over{\sqrt{n}}}\right] \right\rangle &=\left\langle \sum_{r=0}^\infty {{(it)^r}\over{r!}}{{W_1^r}\over{n^{r/2}}} \right\rangle &\cr &=\sum_{r=0}^\infty {{(it)^r}\over{r!}}{{\langle W_1^r\rangle }\over{n^{r/2}}} &\cr &\equiv 1-{{t^2}\over{2n}}+{{\zeta (t,n)}\over{n}}, &(2.1.4.24)\cr}]since [\langle W_1\rangle =0], [\langle W_1^2\rangle =1,] and the quantity denoted by [\zeta (t,n)] in (2.1.4.24)[link] is given by[\zeta (t,n)=\sum_{r=3}^\infty {{(it)^r}\over{r!}}{{\langle W_1^r\rangle }\over{n^{(r/2)-1}}}. \eqno(2.1.4.25)]The characteristic function of [\hat{S}_n] is therefore [\langle \exp (it\hat{S}_n)\rangle =\left[ 1-{{t^2}\over{2n}}+{{\zeta (t,n)}\over n}\right] ^n. \eqno(2.1.4.26)]Now, as is seen from (2.1.4.25)[link], for every fixed t the quantity [\zeta (t,n)] tends to zero as n tends to infinity. The cumulant-generating function of the standardized sum then becomes[\log C_n(\hat{S}_n,t)=n\log \left[ 1-{1 \over n}\left( {{t^2} \over 2}-\zeta(t,n)\right) \right] \eqno(2.1.4.27)]and the logarithm on the right-hand side of equation (2.1.4.27)[link] has the form [\log (1-z)] with [|z|\rightarrow 0] as [n\rightarrow \infty]. We may therefore use the expansion[\log(1-z)=-\left(z+{z^2\over 2}+{z^3\over 3} + \ldots\right),]which is valid for [|z| \,\lt \,1]. We then obtain[\eqalign{\log C_n(\hat{S}_n,t) &=-n\left[{1\over n}\left({t^2\over 2}-\zeta (t,n)\right) +{1\over 2n^2}\left( {t^2\over 2}-\zeta (t,n)\right) ^2\right.\cr&\left.\quad+{1\over 3n^3}\left( {t^2\over 2}-\zeta (t,n)\right) ^3+\cdots \right] \cr&=-{t^2\over 2}+\zeta (t,n)-{1\over 2n}\left( {t^2\over 2}-\zeta (t,n)\right) ^2\cr&\quad-{1\over 3n^2}\left( {t^2\over 2}-\zeta (t,n)\right) ^3-\cdots}]and finally, for every fixed t, [\lim_{n\rightarrow \infty }\log C_n(\hat{S_n},t)=-{{t^2}\over{2}}. \eqno(2.1.4.28)]Since the logarithm is a continuous function of t, it follows directly that [\lim_{n\rightarrow \infty }C_n(\hat{S_n},t)=\exp \left( -{{t^2}\over2}\right) .\eqno(2.1.4.29)]The right-hand side of (2.1.4.29)[link] is just the characteristic function of a standardized normal p.d.f., i.e., a normal p.d.f. with zero mean and unit variance [cf. equation (2.1.4.5)[link]]. The asymptotic expression for the p.d.f. of the standardized sum is therefore obtained as[p(\hat{S})={1 \over{ \sqrt{2\pi}}}\exp\left(-{{\hat{S}^2}\over 2}\right),]which proves the above version of the central-limit theorem.

Surprisingly, this theorem has a very wide applicability and values of n as low as 30 are often large enough for the theorem to be useful. Situations in which the normal p.d.f. must be modified or replaced by an altogether different one are dealt with in Sections 2.1.7[link] and 2.1.8[link] of this chapter.

2.1.4.4. Conditions of validity

| top | pdf |

The above outline of a proof of the central-limit theorem depended on the existence of moments of all orders. The components of structure factors always possess finite moments of all orders, but the existence of moments beyond the second is not necessary for the validity of the theorem and it can be proved under much less stringent conditions. In fact, if all the random variables in equation (2.1.4.19)[link] have the same distribution – as in a homoatomic structure – the only requirement is that the second moments of the distributions should exist [the Lindeberg–Lévy theorem (e.g. Cramér, 1951[link])]. If the distributions are not the same – as in a heteroatomic structure – some further condition is necessary to ensure that no individual random variable dominates the sum. The Liapounoff proof requires the existence of third absolute moments, but this is regarded as aesthetically displeasing; a theorem that ultimately involves only means and variances should require only means and variances in the proof. The Lindeberg–Cramér conditions meet this aesthetic criterion. Roughly, the conditions are that [S^{2}], the variance of the sum, should tend to infinity and [\sigma_{j}^{2}/S^{2}], where [\sigma_{j}^{2}] is the variance of the jth random variable, should tend to zero for all j as n tends to infinity. The precise formulation is quoted by Kendall & Stuart (1977[link], p. 207).

2.1.4.5. Non-independent variables

| top | pdf |

The central-limit theorem, under certain conditions, remains valid even when the variables summed in equation (2.1.4.19)[link] are not independent. The conditions have been investigated by Bernstein (1922[link], 1927[link]); roughly they amount to requiring that the variables should not be too closely correlated. The theorem applies, in particular, when each [x_{r}] is related to a finite number, [f(n)], of its neighbours, when the x's are said to be [f(n)] dependent. The [f(n)] dependence seems plausible for crystallographic applications, since the positions of atoms close together in a structure are closely correlated by interatomic forces, whereas those far apart will show little correlation if there is any flexibility in the asymmetric unit when unconstrained. Harker's (1953[link]) idea of `globs' seems equivalent to [f(n)] dependence. Long-range stereochemical effects, as in pseudo-graphitic aromatic hydrocarbons, would presumably produce long-range correlations and make [f(n)] dependence less plausible. If Bernstein's conditions are satisfied, the central-limit theorem would apply, but the actual value of [\langle x^{2} \rangle - \langle x \rangle^{2}] would have to be used for the variance, instead of the sum of the variances of the random variables in (2.1.4.19)[link]. Because of the correlations the two values are no longer equal.

French & Wilson (1978[link]) seem to have been the first to appeal explicitly to the central-limit theorem extended to non-independent variables, but many previous workers [for typical references, see Wilson (1981[link])] tacitly made the replacement – in the X-ray case substituting the local mean intensity for the sum of the squares of the atomic scattering factors.

2.1.5. Ideal probability density distributions

| top | pdf |

In applications of the central-limit theorem, and its extensions, to intensity statistics the [x_{j}]'s of equation (2.1.4.19)[link] 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[link] through [link] [link] [link] [link] [link] A1.4.3.7[link] , and their first few even moments (fixed-index averaging) are given in Table 2.1.7.1[link]. 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)[link], deduce the m and [\sigma^{2}] 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 space-group origin is chosen at a crystallographic centre and the atoms are nondispersive.

2.1.5.1. Ideal acentric distributions

| top | pdf |

The ideal acentric distributions are obtained by applying the central-limit theorem to the real and the imaginary parts of the structure factor, as given by equation (2.1.1.1)[link]. Consider first a crystal with no rotational symmetry (space group P1). The real part, A, of the structure factor is then given by [A = \textstyle\sum\limits_{j = 1}^{N}f_{j}\cos \vartheta_{j}, \eqno(2.1.5.1)]where N is the number of atoms in the unit cell and [\vartheta_{j}] is the phase angle of the jth atom. The central-limit theorem then states that A tends to be normally distributed about its mean value with variance equal to its mean-square deviation from its mean. Under the assumption that the phase angles [\vartheta_{j}] are uniformly distributed on the 0–2π range, the mean value of each cosine is zero, so that its variance is [\sigma^{2} = \textstyle\sum\limits_{j = 1}^{N}f_{j}^{2}\langle \cos^{2} \vartheta_{j} \rangle. \eqno(2.1.5.2)]Under the same assumption, the mean value of each [\cos^{2} \vartheta] is one-half, so that the variance becomes [\sigma^{2} = (1/2)\textstyle\sum\limits_{j = 1}^{N}f_{j}^{2} = (1/2)\Sigma, \eqno(2.1.5.3)]where Σ is the sum of the squares of the atomic scattering factors [cf. equation (2.1.2.4)[link]]. The asymptotic form of the distribution of A is therefore given by [p(A)\,{\rm d}A = (\pi\Sigma)^{-1/2}\exp(-A^{2}/\Sigma)\,{\rm d}A. \eqno(2.1.5.4)]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 [p(A,B)\,{\rm d}A\,{\rm d}B = (\pi\Sigma)^{-1}\exp[-(A^{2}+B^{2})/\Sigma]\,{\rm d}A\,{\rm d}B. \eqno(2.1.5.5)]Ordinarily, however, we are more interested in the distribution of the magnitude, [|F|], of the structure factor than in the distribution of A and B. Using polar coordinates in equation (2.1.5.5)[link] [[A = |F|\cos\phi,B = |F|\sin\phi]] and integrating over the angle φ gives [p(|F|)\,{\rm d}|F| = (2|F|/\Sigma)\exp(-|F|^{2}/\Sigma)\,{\rm d}|F|. \eqno(2.1.5.6)]It is usually convenient, in structure-factor and intensity statistics, to express the results in terms of the normalized structure factor E and its magnitude [|E|]. If [|F|] has been put on an absolute scale (see Section 2.2.4.3[link] ), we have [E = {{F}\over{\sqrt{\Sigma}}}\quad{\rm and}\quad |E| = {{|F|}\over{\sqrt{\Sigma}}}, \eqno(2.1.5.7)]so that [p(|E|)\,{\rm d}|E| = 2|E|\exp(-|E|^{2})\,{\rm d}|E| \eqno(2.1.5.8)]is the normalized-structure-factor version of (2.1.5.6)[link].

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.

2.1.5.2. Ideal centric distributions

| top | pdf |

When a nondispersive crystal is centrosymmetric, and the space-group 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 [P\bar{1}], the contribution of the jth atom plus its centrosymmetric counterpart is [2f_{j}\cos\vartheta_{j}]. The calculation of [p(A)] goes through as before, with allowance for the fact that there are [N/2] pairs instead of N independent atoms, giving [p(A)\,{\rm d}A = (2\pi\Sigma)^{-1/2}\exp[-A^{2}/(2\Sigma)]\,{\rm d}A \eqno(2.1.5.9)]or equivalently [p(|F|)\,{\rm d}|F| = [2/(\pi\Sigma)]^{1/2}\exp[-|F|^{2}/(2\Sigma)]\,{\rm d}|F| \eqno(2.1.5.10)]or [p(|E|)\,{\rm d}|E| = (2/\pi)^{1/2}\exp(-|E|^{2}/2)\,{\rm d}|E|. \eqno(2.1.5.11)]

2.1.5.3. Effect of other symmetry elements on the ideal acentric and centric distributions

| top | pdf |

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[link]) and Rogers (1950[link]), 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[link] and 2.1.3.2[link]. 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[link] 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)[link] and (2.1.5.11)[link] are not affected.

2.1.5.4. Other ideal distributions

| top | pdf |

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, [\eqalignno{p(|E|)\,{\rm d}|E| &= {{2|E|}\over{(1-k^{2})^{1/2}}}\exp[-|E|^{2}/(1-k^{2})] & \cr &\quad\times{} I_{0}\left({{k|E|^{2}}\over{1-k^{2}}}\right)\,{\rm d}|E|, &(2.1.5.12) \cr}]where [I_{0}(x)] 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[link], ch. III; Wilson, 1980a[link],b[link]; Shmueli & Wilson, 1983[link]). The bicentric distribution [p(|E|)\,{\rm d}|E| = \pi^{-3/2}\exp(-|E|^{2}/8)K_{0}(|E|^{2}/8)\,{\rm d}|E| \eqno(2.1.5.13)]arises, for example, when the `asymmetric unit in a centrosymmetric crystal is a centrosymmetric molecule' (Lipson & Woolfson, 1952[link]); [K_{0}(x)] is a modified Bessel function of the second kind. There are higher hypercentric, hyperparallel and sesquicentric analogues (Wilson, 1952[link]; Rogers & Wilson, 1953[link]; Wilson, 1956[link]). 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[link]).

2.1.5.5. Relation to distributions of I

| top | pdf |

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 [(|F|)], since both are necessarily real and non-negative. In an obvious notation, the relation between the intensity distribution and the structure-factor distribution is [p_{I}(I) = (1/2)I^{-1/2}p_{|F|}(I^{1/2}) \eqno(2.1.5.14)]or [p_{|F|}(|F|) = 2|F|p_{I}(|F|^{2}). \eqno(2.1.5.15)]Statistical fluctuations in counting rates, however, introduce a small but finite probability of negative observed intensities (Wilson, 1978a[link], 1980a[link]) and thus of imaginary structure factors. This practical complication is treated in IT C (2004[link], Parts 7[link] and 8[link] ).

Both the ideal centric and acentric distributions are simple members of the family of gamma distributions, defined by [\gamma_{n}(x)\,{\rm d}x = [\Gamma(n)]^{-1}x^{n-1}\exp(-x)\,{\rm d}x, \eqno(2.1.5.16)]where n is a parameter, not necessarily integral, and [\Gamma(n)] is the gamma function. Thus the ideal acentric intensity distribution is [\eqalignno{p(I)\,{\rm d}I & = \exp(-I/\Sigma)\,{\rm d}(I/\Sigma) &(2.1.5.17) \cr & = \gamma_{1}(I/\Sigma)\,{\rm d}(I/\Sigma) &(2.1.5.18)}%2.1.5.18]and the ideal centric intensity distribution is [\eqalignno{p(I)\,{\rm d}I& = (2\Sigma/\pi)^{1/2}\exp[-I/(2\Sigma)]\,{\rm d}[I/(2\Sigma)] &(2.1.5.19) \cr & = \gamma_{1/2}[I/(2\Sigma)]\,{\rm d}[I/(2\Sigma)]. &(2.1.5.20)}%2.1.5.20]The properties of gamma distributions and of the related beta distributions, summarized in Table 2.1.5.1[link], are used in Section 2.1.6[link] to derive the probability density functions of sums and of ratios of intensities drawn from one of the ideal distributions.

Table 2.1.5.1| top | pdf |
Some properties of gamma and beta distributions

If [x_{1}, x_{2}, \ldots, x_{n}] are independent gamma-distributed variables with parameters [p_{1}, p_{2}, \ldots, p_{n}], their sum is a gamma-distributed variable with [p\ =] [ p_{1} + p_{2} + \ldots + p_{n}].

If x and y are independent gamma-distributed variables with parameters p and q, then the ratio [u = x/y] has the distribution [\beta_{2} (u\semi\ p, q)].

With the same notation, the ratio [v = x/(x + y)] has the distribution [\beta_{1} (v\semi \ p, q)].

Differences and products of gamma-distributed variables do not lead to simple results. For proofs, details and references see Kendall & Stuart (1977)[link].

Name of the distribution, its functional form, mean and variance
Gamma distribution with parameter p: [\gamma_{p} (x) = [\Gamma (x)]^{-1} x^{p-1} \exp (-x)\semi \quad p \leq x \leq \infty,\quad p > 0][\hbox{mean: }\langle x\rangle = p\semi \quad \hbox{variance: } \langle (x - \langle x\rangle)^{2}\rangle = p.]
Beta distribution of first kind with parameters p and q: [\beta_{1} (x\semi\ p, q) = {\Gamma (p + q) \over \Gamma (p) \Gamma (q)} x^{p - 1} (1 - x)^{q - 1}\semi \quad 0 \leq x \leq \infty,\quad p, q > 0][\hbox{mean: }\langle x\rangle = p/(p + q)\semi ][\hbox{variance: }\langle (x - \langle x \rangle)^{2}\rangle = pq/[(p + q)^{2} (p + q + 1)].]
Beta distribution of second kind with parameters p and q: [\beta_{2} (x\semi\ p, q) = {\Gamma (p + q) \over \Gamma (p) \Gamma (q)} x^{p - 1} (1 + x)^{-p -q}\semi \quad 0 \leq x \leq \infty,\quad p, q > 0][\hbox{mean: }\langle x \rangle = p/(q - 1)\semi][\hbox{variance: }\langle (x - \langle x \rangle)^{2}\rangle = p(p + q - 1)/[(q - 1) (q - 2)].]

2.1.5.6. Cumulative distribution functions

| top | pdf |

The integral of the probability density function [f(x)] from the lower end of its range up to an arbitrary value x is called the cumulative probability distribution, or simply the distribution function, [F(x)], of x. It can always be written [F(x) = \textstyle\int\limits_{-\infty}^{x}f(u)\,{\rm d}u\semi \eqno(2.1.5.21)]if the lower end of its range is not actually [-\infty] one takes [f(x)] as identically zero between [-\infty] and the lower end of its range. For the distribution of A [equation (2.1.5.4)[link] or (2.1.5.9)[link]] the lower limit is in fact [-\infty]; for the distribution of [|F|], [|E|], I and [I/\Sigma] the lower end of the range is zero. In such cases, equation (2.1.5.21)[link] becomes [F(x) = \textstyle\int\limits_{0}^{x}f(x)\,{\rm d}x. \eqno(2.1.5.22)]In crystallographic applications the cumulative distribution is usually denoted by [N(x)], 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[link]) have found many applications. For the acentric distribution of [|E|] [equation (2.1.5.8)[link]] the integration is readily carried out: [N(|E|) = 2\textstyle\int\limits_{0}^{|E|} y\exp(-y^{2})\,{\rm d}y = 1 - \exp(-|E|^{2}). \eqno(2.1.5.23)]The integral for the centric distribution of [|E|] [equation (2.1.5.11)[link]] 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 [{\rm erf}(x) = (2/\pi^{1/2})\textstyle\int\limits_{0}^{x}\exp(-t^{2})\,{\rm d}t. \eqno(2.1.5.24)]For the centric distribution, then [\eqalignno{N(|E|) & = (2/\pi)^{1/2}\textstyle\int\limits_{0}^{|E|}y\exp(-y^{2}/2)\,{\rm d}y &(2.1.5.25) \cr & = {\rm erf}(|E|/2^{1/2}). &(2.1.5.26)}%2.1.5.26]The error function is extensively tabulated [see e.g. Abramowitz & Stegun (1972[link]), pp. 310–311, and a closely related function on pp. 966–973].

The foregoing and following structure-factor 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 X-ray 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 convergent-beam electron diffraction technique using an electron microscope with a nanoscale probe, as described and illustrated in Section 2.5.3[link] of this volume.

2.1.6. Distributions of sums, averages and ratios

| top | pdf |

2.1.6.1. Distributions of sums and averages

| top | pdf |

In Section 2.1.2.1[link], it was shown that the average intensity of a sufficient number of reflections is [\Sigma] [equation (2.1.2.4)[link]]. When the number of reflections is not `sufficient', their mean value will show statistical fluctuations about [\Sigma]; such statistical fluctuations are in addition to any systematic variation resulting from non-independence of atomic positions, as discussed in Sections 2.1.2.1–2.1.2.3[link][link][link]. We thus need to consider the probability density functions of sums like [J_{n} = \textstyle\sum\limits_{i = 1}^{n}G_{i}, \eqno(2.1.6.1)]and averages like [Y = J_{n}/n, \eqno(2.1.6.2)]where [G_{i}] is the intensity of the ith reflection. The probability density distributions are easily obtained from a property of gamma distributions: If [x_{1}, x_{2}, \ldots, x_{n}] are independent gamma-distributed variables with parameters [p_{1}, p_{2}, \ldots, p_{n}], their sum is a gamma-distributed variable with parameter p equal to the sum of the parameters. The sum of n intensities drawn from an acentric distribution thus has the distribution [p(J_{n})\,{\rm d}J_{n} = \gamma_{n}(J_{n}/\Sigma)\,{\rm d}(J_{n}/\Sigma)\semi \eqno(2.1.6.3)]the parameters of the variables added are all equal to unity, so that their sum is p. Similarly, the sum of n intensities drawn from a centric distribution has the distribution [p(J_{n})\,{\rm d}J_{n} = \gamma_{n/2}[J_{n}/(2\Sigma)]\,{\rm d}[J_{n}/(2\Sigma)]\semi \eqno(2.1.6.4)]each parameter has the value of one-half. The corresponding distributions of the averages of n intensities are then [p(Y)\,{\rm d}Y = \gamma_{n}(nY/\Sigma)\,{\rm d}(nY/\Sigma) \eqno(2.1.6.5)]for the acentric case, and [p(Y)\,{\rm d}Y = \gamma_{n/2}[nY/(2\Sigma)]\,{\rm d}[nY/(2\Sigma)] \eqno(2.1.6.6)]for the centric. In both cases the expected value of Y is [\Sigma] and the variances are [\Sigma^{2}/n] and [2\Sigma^{2}/n], respectively, just as would be expected.

2.1.6.2. Distribution of ratios

| top | pdf |

Ratios like [S_{n, \, m} = J_{n}/K_{m}, \eqno(2.1.6.7)]where [J_{n}] is given by equation (2.1.6.1)[link], [K_{m} = \textstyle\sum\limits_{j = 1}^{m}H_{j}, \eqno(2.1.6.8)]and the [H_{j}]'s are the intensities of a set of reflections (which may or may not overlap with those included in [J_{n}]), are used in correlating intensities measured under different conditions. They arise in correlating reflections on different layer lines from the same or different specimens, in correlating the same reflections from different crystals, in normalizing intensities to the local average or to [\Sigma], and in certain systematic trial-and-error methods of structure determination (see Rabinovich & Shakked, 1984[link], and references therein). There are three main cases:

  • (i) [G_{i}] and [H_{i}] refer to the same reflection; for example, they might be the observed and calculated quantities for the [hkl] reflection measured under different conditions or for different crystals of the same substance; or

  • (ii) [G_{i}] and [H_{i}] are unrelated; for example, the observed and calculated values for the [hkl] reflection for a completely wrong trial structure, of values for entirely different reflections, as in reducing photographic measurements on different layer lines to the same scale; or

  • (iii) the [G_{i}]'s are a subset of the [H_{i}]'s, so that [G_{i} = H_{i}] for [i \,\lt\, n] and [m \,\gt\, n].

Aside from the scale factor, in case (i)[link] [G_{i}] and [H_{i}] will differ chiefly through relatively small statistical fluctuations and uncorrected systematic errors, whereas in case (ii)[link] the differences will be relatively large because of the inherent differences in the intensities. Here we are concerned only with cases (ii)[link] and (iii)[link]; the practical problems of case (i)[link] are postponed to IT C (2004[link]), Chapter 7.5[link] .

There is little in the crystallographic literature concerning the probability distribution of sums like (2.1.6.1)[link] or ratios like (2.1.6.7)[link]; certain results are reviewed by Srinivasan & Parthasarathy (1976[link], ch. 5), but with a bias toward partially related structures that makes it difficult to apply them to the immediate problem.

In case (ii)[link] ([G_{i}] and [H_{i}] independent), acentric distribution, Table 2.1.5.1[link] gives the distribution of the ratio [u = nY/(mZ) \eqno(2.1.6.9)][p(u)\,{\rm d}u = \beta_{2}[nY/(mZ)\hbox{; } n, m]\,{\rm d}[nY/(mZ)], \eqno(2.1.6.10)]where [\beta_{2}] is a beta distribution of the second kind, Y is given by equation (2.1.6.2)[link] and Z by [Z = K_{m}/m, \eqno(2.1.6.11)]where n is the number of intensities included in the numerator and m is the number in the denominator. The expected value of [Y/Z] is then [\langle Y/Z \rangle = {{m}\over{m-1}} = 1+{{1}\over{m}}+\ldots \eqno(2.1.6.12)]with variance [\sigma^{2} = {{(n+m-1)m^{2}}\over{(m-1)^{2}(m-2)n}}. \eqno(2.1.6.13)]One sees that [Y/Z] is a biased estimate of the scaling factor between two sets of intensities and the bias, of the order of [m^{-1}], depends only on the number of intensities averaged in the denominator. This may seem odd at first sight, but it becomes plausible when one remembers that the mean of a quantity is an unbiased estimator of itself, but the reciprocal of a mean is not an unbiased estimator of the mean of a reciprocal. The mean exists only if [m \,\gt\, 1] and the variance only for [m \,\gt\, 2].

In the centric case, the expression for the distribution of the ratio of the two means Y and Z becomes [p(u)\,{\rm d}u = \beta_{2}[nY/(mZ)\semi \ n/2, m/2]\,{\rm d}[nY/(mZ)] \eqno(2.1.6.14)]with the expected value of [Y/Z] equal to [\langle Y/Z \rangle = {{m}\over{m-2}} = 1+{{2}\over{m}}+\ldots \eqno(2.1.6.15)]and with its variance equal to [\sigma^{2} = {{2(n+m-2)m^{2}}\over{(m-2)^{2}(m-4)n}}. \eqno(2.1.6.16)]For the same number of reflections, the bias in [\langle Y/Z \rangle] and the variance for the centric distribution are considerably larger than for the acentric. For both distributions the variance of the scaling factor approaches zero when n and m become large. The variances are large for m small, in fact `infinite' if the number of terms averaged in the denominator is sufficiently small. These biases are readily removed by multiplying [Y/Z] by [(m-1)/m] or [(m-2)/m]. Many methods of estimating scaling factors – perhaps most – also introduce bias (Wilson, 1975[link]; Lomer & Wilson, 1975[link]; Wilson, 1976[link], 1978c[link]) that is not so easily removed. Wilson (1986a[link]) has given reasons for supposing that the bias of the ratio (2.1.6.7)[link] approximates to [1+{{\sigma^{2}(I)}\over{m\langle I \rangle^{2}}}, \eqno(2.1.6.17)]whatever the intensity distribution. Equations (2.1.6.12)[link] and (2.1.6.15)[link] are consistent with this.

2.1.6.3. Intensities scaled to the local average

| top | pdf |

When the [G_{i}]'s are a subset of the [H_{i}]'s, the beta distributions of the second kind are replaced by beta distributions of the first kind, with means and variances readily found from Table 2.1.5.1[link]. The distribution of such a ratio is chiefly of interest when Y relates to a single reflection and Z relates to a group of m intensities including Y. This corresponds to normalizing intensities to the local average. Its distribution is [p(I/\langle I \rangle)\,{\rm d}(I/\langle I \rangle) = \beta_{1}(I/n\langle I \rangle\semi\ 1,n-1)\,{\rm d}(I/n\langle I \rangle) \eqno(2.1.6.18)]in the acentric case, with an expected value of [I/\langle I \rangle] of unity; there is no bias, as is obvious a priori. The variance of [I/\langle I \rangle] is [\sigma^{2} = {{n-1}\over{n+1}}, \eqno(2.1.6.19)]which is less than the variance of the intensities normalized to an `infinite' population by a fraction of the order of [2/n]. Unlike the variance of the scaling factor, the variance of the normalized intensity approaches unity as n becomes large. For intensities having a centric distribution, the distribution normalized to the local average is given by [{p(I/\langle I \rangle)\,{\rm d}(I/\langle I \rangle) = \beta_{1}[I/n\langle I \rangle\semi\ 1/2,(n-1)/2]\,{\rm d}(I/n\langle I \rangle),} \eqno(2.1.6.20)]with an expected value of [I/\langle I \rangle] of unity and with variance [\sigma^{2} = {{2(n-1)}\over{n+2}}, \eqno(2.1.6.21)]less than that for an `infinite' population by a fraction of about [3/n].

Similar considerations apply to intensities normalized to Σ in the usual way, since they are equal to those normalized to [\langle I \rangle] multiplied by [\langle I \rangle/\Sigma].

2.1.6.4. The use of normal approximations

| top | pdf |

Since [J_{n}] and [K_{m}] [equations (2.1.6.1)[link] and (2.1.6.8)[link]] are sums of identically distributed variables conforming to the conditions of the central-limit theorem, it is tempting to approximate their distributions by normal distributions with the correct mean and variance. This would be reasonably satisfactory for the distributions of [J_{n}] and [K_{m}] themselves for quite small values of n and m, but unsatisfactory for the distribution of their ratio for any values of n and m, even large. The ratio of two variables with normal distributions is notorious for its rather indeterminate mean and `infinite' variance, resulting from the `tail' of the denominator distributions extending through zero to negative values. The leading terms of the ratio distribution are given by Kendall & Stuart (1977[link], p. 288).

2.1.7. Non-ideal distributions: the correction-factor approach

| top | pdf |

2.1.7.1. Introduction

| top | pdf |

The probability density functions (p.d.f.'s) of the magnitude of the structure factor, presented in Section 2.1.5[link], are based on the central-limit theorem discussed above. In particular, the centric and acentric p.d.f.'s given by equations (2.1.5.11)[link] and (2.1.5.8)[link], respectively, are expected to account for the statistical properties of diffraction patterns obtained from crystals consisting of nearly equal atoms, which obey the fundamental assumptions of uniformity and independence of the atomic contributions and are not affected by noncrystallographic symmetry and dispersion. It is also assumed there that the number of atoms in the asymmetric unit is large. Distributions of structure-factor magnitudes which are based on the central-limit theorem, and thus obey the above assumptions, have been termed `ideal', and the subjects of the following sections are those distributions for which some of the above assumptions/restrictions are not fulfilled; the latter distributions will be called `non-ideal'.

We recall that the assumption of uniformity consists of the requirement that the fractional part of the scalar product [hx+ky+lz] be uniformly distributed over the [0, 1] interval, which holds well if [x, y, z] are rationally independent (Hauptman & Karle, 1953[link]), and permits one to regard the atomic contribution to the structure factor as a random variable. This is of course a necessary requirement for any statistical treatment. If, however, the atomic composition of the asymmetric unit is widely heterogeneous, the structure factor is then a sum of unequally distributed random variables and the Lindeberg–Lévy version of the central-limit theorem (cf. Section 2.1.4.4[link]) cannot be expected to apply. Other versions of this theorem might still predict a normal p.d.f. of the sum, but at the expense of a correspondingly large number of terms/atoms. It is well known that atomic heterogeneity gives rise to severe deviations from ideal behaviour (e.g. Howells et al., 1950[link]) and one of the aims of crystallographic statistics has been the introduction of a correct dependence on the atomic composition into the non-ideal p.d.f.'s [for a review of the early work on non-ideal distributions see Srinivasan & Parthasarathy (1976[link])]. A somewhat less well known fact is that the dependence of the p.d.f.'s of [|E|] on space-group symmetry becomes more conspicuous as the composition becomes more heterogeneous (e.g. Shmueli, 1979[link]; Shmueli & Wilson, 1981[link]). Hence both the composition and the symmetry dependence of the intensity statistics are of interest. Other problems, which likewise give rise to non-ideal p.d.f.'s, are the presence of heavy atoms in (variable) special positions, heterogeneous structures with complete or partial noncrystallographic symmetry, and the presence of outstandingly heavy dispersive scatterers.

The need for theoretical representations of non-ideal p.d.f.'s is exemplified in Fig. 2.1.7.1[link], which shows the ideal centric and acentric p.d.f.'s together with a frequency histogram of [|E|] values, recalculated for a centrosymmetric structure containing a platinum atom in the asymmetric unit of [P\bar{1}] (Faggiani et al., 1980[link]). Clearly, the deviation from the Gaussian p.d.f., predicted by the central-limit theorem, is here very large and a comparison with the possible ideal distributions can (in this case) lead to wrong conclusions.

[Figure 2.1.7.1]

Figure 2.1.7.1 | top | pdf |

Atomic heterogeneity and intensity statistics. The histogram appearing in this figure was constructed from [|E|] values which were recalculated from atomic parameters published for the centrosymmetric structure of C6H18Cl2N4O4Pt (Faggiani et al., 1980[link]). The space group of the crystal is [P\bar{1}], [Z=2], i.e. all the atoms are located in general positions. The figure shows a comparison of the recalculated distribution of [|E|] with the ideal centric [equation (2.1.5.11)[link]] and acentric [equation (2.1.5.8)[link]] p.d.f.'s, denoted by [\bar{1}] and 1, respectively.

Two general approaches have so far been employed in derivations of non-ideal p.d.f.'s which account for the above-mentioned problems: the correction-factor approach, to be dealt with in the following sections, and the more recently introduced Fourier method, to which Section 2.1.8[link] is dedicated. In what follows, we introduce briefly the mathematical background of the correction-factor approach, apply this formalism to centric and acentric non-ideal p.d.f.'s, and present the numerical values of the moments of the trigonometric structure factor which permit an approximate evaluation of such p.d.f.'s for all the three-dimensional space groups.

2.1.7.2. Mathematical background

| top | pdf |

Suppose that [p(x)] is a p.d.f. which accurately describes the experimental distribution of the random variable x, where x is related to a sum of random variables and can be assumed to obey (to some approximation) an ideal p.d.f., say [p^{(0)}(x)], based on the central-limit theorem. In the correction-factor approach we seek to represent [p(x)] as [p(x) = p^{(0)}(x)\textstyle\sum\limits_{k}d_{k}\, f_{k}(x), \eqno(2.1.7.1)]where [d_{k}] are coefficients which depend on the cause of the deviation of [p(x)] from the central-limit theorem approximation and [f_{k}(x)] are suitably chosen functions of x. A choice of the set [\{f_{k}\}] is deemed suitable, if only from a practical point of view, if it allows the convenient introduction of the cause of the above deviation of [p(x)] into the expansion coefficients [d_{k}]. This requirement is satisfied – also from a theoretical point of view – by taking [f_{k}(x)] as a set of polynomials which are orthogonal with respect to the ideal p.d.f., taken as their weight function (e.g. Cramér, 1951[link]). That is, the functions [f_{k}(x)] so chosen have to obey the relationship[\textstyle\int\limits_{a}^{b}f_{k}(x)f_{m}(x)p^{(0)}(x)\,{\rm d}x = \delta_{km} = \cases{ 1, & if\quad $k = m$\cr 0, & if\quad $k \neq m $}\,, \eqno(2.1.7.2)]where [[a,b]] is the range of existence of all the functions involved. It can be readily shown that the coefficients [d_{k}] are given by [d_{k} = \textstyle\int\limits_{a}^{b}f_{k}(x)p(x)\,{\rm d}x = \langle f_{k}(x) \rangle = \textstyle\sum\limits_{n = 0}^{k}c_{n}^{(k)}\langle x^{n} \rangle, \eqno(2.1.7.3)]where the brackets [\langle \, \rangle] in equation (2.1.7.3)[link] denote averaging with respect to the unknown p.d.f. [p(x)] and [c_{n}^{(k)}] is the coefficient of the nth power of x in the polynomial [f_{k}(x)]. The coefficients [d_{k}] are thus directly related to the moments of the non-ideal distribution and the coefficients of the powers of x in the orthogonal polynomials. The latter coefficients can be obtained by the Gram–Schmidt procedure (e.g. Spiegel, 1974[link]), or by direct use of the Szegö determinants (e.g. Cramér, 1951[link]), for any weight function that has finite moments. However, the feasibility of the present approach depends on our ability to obtain the moments [\langle x^{n} \rangle] without the knowledge of the non-ideal p.d.f., [p(x)].

2.1.7.3. Application to centric and acentric distributions

| top | pdf |

We shall summarize here the non-ideal centric and acentric distributions of the magnitude of the normalized structure factor E (e.g. Shmueli & Wilson, 1981[link]; Shmueli, 1982[link]). We assume that (i) all the atoms are located in general positions and have rationally independent coordinates, (ii) all the scatterers are dispersionless, and (iii) there is no noncrystallographic symmetry. Arbitrary atomic composition and space-group symmetry are admitted. The appropriate weight functions and the corresponding orthogonal polynomials are [\matrix {& & \hbox{Non-ideal} \cr p^{(0)}(|E|) & f_k(x) & \hbox{distribution} \cr\noalign{\vskip -7pt}\cr\noalign{\hrule}\cr & & \cr (2/\pi)^{1/2}\exp(-|E|^2/2) & He_{2k}(|E|)/[(2k)!]^{1/2} & \hbox{Centric} \cr & & \cr 2|E|\exp(-|E|^2) & L_k(|E|^2) & \hbox{Acentric}\cr \noalign{\vskip -7pt}\cr\noalign{\hrule} \cr} \eqno(2.1.7.4)]where [He_{k}] and [L_{k}] are Hermite and Laguerre polynomials, respectively, as defined, for example, by Abramowitz & Stegun (1972[link]). Equations (2.1.7.2)[link], (2.1.7.3)[link] and (2.1.7.4)[link] suffice for the general formulation of the above non-ideal p.d.f.'s of [|E|]. Their full derivation entails (i) the expression of a sufficient number of moments of [|E|] in terms of absolute moments of the trigonometric structure factor (e.g. Shmueli & Wilson, 1981[link]; Shmueli, 1982[link]) and (ii) calculation of the latter moments for the various symmetries (Wilson, 1978b[link]; Shmueli & Kaldor, 1981[link], 1983[link]). The notation below is similar to that employed by Shmueli (1982[link]).

These non-ideal p.d.f.'s of [|E|], for which the first five expansion terms are available, are given by [p_{\rm c}(|E|) = p_{\rm c}^{(0)}(|E|)\left[1+\sum_{k = 2}^{\infty}{{A_{2k}}\over{(2k)!}} He_{2k}(|E|)\right] \eqno(2.1.7.5)]and [p_{\rm a}(|E|) = p_{\rm a}^{(0)}(|E|)\left[1+\sum_{k = 2}^{\infty}{{(-1)^{k}B_{2k}}\over{k!}} L_{k}(|E|^{2})\right] \eqno(2.1.7.6)]for centrosymmetric and noncentrosymmetric space groups, respectively, where [p_{\rm c}^{(0)}(|E|)] and [p_{\rm a}^{(0)}(|E|)] are the ideal centric and acentric p.d.f.'s [see (2.1.7.4)[link]] and the unified form of the coefficients [A_{2k}] and [B_{2k}], for k = 2, 3, 4 and 5, is [\eqalign{&A_{4}\cr &A_{6}\cr &A_{8}\cr &A_{10}\cr &\cr}\eqalign{&\hbox{ or } \cr &\hbox{ or } \cr &\hbox{ or } \cr &\hbox{ or }\cr &\cr}\eqalign{B_{4} & = a_{4}Q_{4} \cr B_{6} & = a_{6}Q_{6} \cr B_{8} & = a_{8}Q_{8} + U(a_{4}^{2}Q_{4}^{2}-\gamma_{4}^{2})\cr B_{10} & = a_{10}Q_{10} + V(a_{4}a_{6}Q_{4}Q_{6}-\gamma_{4}\gamma_{6}Q_{10})\cr &\quad{}+ W\gamma_{4}^{2}Q_{10} }\eqno(2.1.7.7)](Shmueli, 1982[link]), where U = 35 or 18, V = 210 or 100 and W = 3150 or 900 according as [A_{2k}] or [B_{2k}] is required, respectively, and the other quantities in equation (2.1.7.7)[link] are given below. The composition-dependent terms in equations (2.1.7.7)[link] are [Q_{2k} = {{\sum_{j = 1}^{m}f_{j}^{2k}}\over{\left(\sum_{n = 1}^{m}f_{n}^{2}\right)^{k}}}, \eqno(2.1.7.8)]where m is the number of atoms in the asymmetric unit, [f_{j},\, j = 1,\ldots,m] are their scattering factors, and the symmetry dependence is expressed by the coefficients [a_{2k}] in equation (2.1.7.7)[link], as follows: [{a_{2k} = (-1)^{k-1}(k-1)!\alpha_{k0} +\textstyle\sum\limits_{p = 2}^{k}(-1)^{k-p}(k-p)!\alpha_{kp} \gamma_{2p},} \eqno(2.1.7.9)]where [\alpha_{kp} = \left(\matrix{ k \cr p } \right) {{(2k-1)!!}\over{(2p-1)!!}} \quad{\rm or}\quad\left(\matrix{ k \cr p } \right) {{k!}\over{p!}} \eqno(2.1.7.10)]according as the space group is centrosymmetric or noncentrosymmetric, respectively, and [\gamma_{2p}] in equation (2.1.7.9)[link] is given by [\gamma_{2p} = {{\langle |T|^{2p} \rangle}\over{\langle |T|^{2} \rangle^{p}}}, \eqno(2.1.7.11)]where [\langle |T|^{k} \rangle] is the kth absolute moment of the trigonometric structure factor [T({\bf h}) = \textstyle\sum\limits_{s = 1}^{g}\exp[2\pi i{\bf h}^{T}({\bf P}_{s}{\bf r}+{\bf t}_{s})] \equiv \xi({\bf h})+i\eta({\bf h}). \eqno(2.1.7.12)]In equation (2.1.7.12)[link], g is the number of general equivalent positions listed in IT A (2005)[link] for the space group in question, times the multiplicity of the Bravais lattice, [({\bf P}_{s},{\bf t}_{s})] is the sth space-group operator and [{\bf r}] is an atomic position vector.

The cumulative distribution functions, obtained by integrating equations (2.1.7.5)[link] and (2.1.7.6)[link], are given by [\eqalignno{N_{\rm c}(|E|)& = {\rm erf}\left({{|E|}\over{\sqrt{2}}}\right)-{{2}\over{\sqrt{\pi}}} \exp\left(-{{|E|^{2}}\over{2}}\right) \cr &\quad\times \left[\sum_{k = 2}^{\infty}{{A_{2k}}\over{(2k)!}} He_{2k-1}(|E|)\right] &(2.1.7.13)}]and [\eqalignno{N_{\rm a}(|E|) & = 1 - \exp(-|E|^{2}) + \exp(-|E|^{2}) \cr &\quad\times \left\{\sum_{k = 2}^{\infty}{{(-1)^{k}B_{2k}}\over{k!}}[L_{k-1}(|E|^{2}) - L_{k}(|E|^2)]\right\}\cr &&(2.1.7.14)}]for centrosymmetric and noncentrosymmetric space groups, respectively, where the coefficients are defined in equations (2.1.7.7)[link]–(2.1.7.12)[link][link][link][link][link]. Note that the first term on the right-hand side of equation (2.1.7.13)[link] and the first two terms on the right-hand side of equation (2.1.7.14)[link] are just the cumulative distributions derived from the ideal centric and acentric p.d.f.'s in Section 2.1.5.6[link].

The moments [\langle |T|^{2k} \rangle] were compiled for all the space groups by Wilson (1978b[link]) for [k =] 1 and 2, and by Shmueli & Kaldor (1981[link], 1983[link]) for [k =] 1, 2, 3 and 4. These results are presented in Table 2.1.7.1[link]. Closed expressions for the normalized moments [\gamma_{2p}] were obtained by Shmueli (1982[link]) for the triclinic, monoclinic and orthorhombic space groups except [Fdd2] and [Fddd] (see Table 2.1.7.2[link]). The composition-dependent terms, [Q_{2k}], are most conveniently computed as weighted averages over the ranges of [(\sin\theta)/\lambda] which were used in the construction of the Wilson plot for the computation of the [|E|] values.

Table 2.1.7.1| top | pdf |
Some even absolute moments of the trigonometric structure factor

The symbols p, q, r and s denote the second, fourth, sixth and eighth absolute moments of the trigonometric structure factor T [equation (2.1.7.12)[link]], respectively, and the columns of the table contain (for some conciseness) [p, q, r/p] and [s/p^{2}]. The numbers in parentheses, appearing beside some space-group entries, refer to hkl subsets which are defined in the note at the end of the table. These subset references are identical with those given by Shmueli & Kaldor (1981[link], 1983[link]). The symbols q, r and s are also equivalent to [\gamma_{4}P^{2}], [\gamma_{6}P^{3}] and [\gamma_{8}P^{4}], respectively, where [\gamma_{2n}] are the normalized absolute moments given by equation (2.1.7.11)[link].

Space groups(s)pq[r/p][s/p^{2}]
Point group: 1        
P1 1 1 1 1
Point group: [\bar{\bf 1}]        
[P\bar{1}] 2 6 10 17½
Point groups: 2, m        
All P 2 6 10 17½
All C 4 48 160 560
Point group: [{\bf 2/}{\bi m}]        
All P 4 36 100 306¼
All C 8 288 1600 9800
Point group: 222        
All P 4 28 64 169¾
All C and I 8 224 1024 5432
F222 16 1792 16384 173824
Point group: mm2        
All P 4 36 100 306¼
All A, C and I 8 288 1600 9800
Fmm2 16 2304 25600 313600
Fdd2 (1) 16 2304 25600 313600
Fdd2 (2) 16 1280 7168 43264
Point group: mmm        
All P 8 216 1000 5359[3\over8]
All C and I 16 1728 16000 171500
Fmmm 32 13824 256000 5488000
Fddd (1) 32 13824 256000 5488000
Fddd (2) 32 7680 71680 757120
Point group: 4        
[P4, P4_{2}] 4 36 100 306¼
[P4_{1}] (3) 4 36 100 306¼
[P4_{1}] (4) 4 20 28 42¼
[I4] 8 288 1600 9800
[I4_{1}] (5) 8 288 1600 9800
[I4_{1}] (6) 8 160 448 1352
Point group: [\bar{\bf 4}]        
[P\bar{4}] 4 28 64 169¾
[I\bar{4}] 8 224 1024 5432
Point group: [{\bf 4}/{\bi m}]        
All P 8 216 1000 5359[3\over8]
[I4/m] 16 1728 16000 171500
[I4_{1}/a] (7) 16 1728 16000 171500
[I4_{1}/a] (8) 16 960 4480 23660
Point group: 422        
[P422], [P42_{1}2], [P4_{2}22], [P4_{2}2_{1}2] 8 136 424 1682[1\over8]
[P4_{1}22], [P4_{1}2_{1}2] (3) 8 136 424 1682[1\over8]
[P4_{1}22], [P4_{1}2_{1}2] (4) 8 104 208 470[1\over8]
I422 16 1088 6784 53828
[I4_{1}22] (7) 16 1088 6784 53828
[I4_{1}22] (8) 16 832 3328 15044
Point group: 4mm        
All P 8 168 640 2970[5\over8]
I4mm, I4cm 16 1344 10240 95060
[I4_{1}md, I4_{1}cd] (7) 16 1344 10240 95060
[I4_{1}md, I4_{1}cd] (8) 16 832 3328 15188
Point groups: [\bar{\bf 4}{\bf 2}{\bi m},\bar{\bf 4}{\bi m}{\bf 2}]        
All P 8 136 424 1682[1\over8]
[I\bar{4}m2, I\bar{4}2m, I\bar{4}c2] 16 1088 6784 53828
[I\bar{4}2d] (5) 16 1088 6784 53828
[I\bar{4}2d] (6) 16 832 3328 15044
Point group: 4/mmm        
All P 16 1008 6400 51985[15\over16]
[I4/mmm], [I4/mcm] 32 8064 102400 1663550
[I4_{1}/amd, I4_{1}/acd] (5) 32 8064 102400 1663550
[I4_{1}/amd, I4_{1}/acd] (6) 32 4992 33280 265790
Point group: 3        
All P and R 3 15 31 71
Point group: [\bar{\bf 3}]        
All P and R 6 90 310 1242½
Point group: 32        
All P and R 6 66 166 508½
Point group: 3m        
P3m1, P31m, R3m 6 66 178 604½
P3c1, P31c (3), R3c (1) 6 66 178 604½
P3c1, P31c (4), R3c (2) 6 66 154 412½
Point group: [\bar{\bf 3}{\bi m}]        
[P\bar{3}1m, P\bar{3}m1, R\bar{3}m] 12 396 1780 10578¾
[P\bar{3}1c, P\bar{3}c1] (3), [R\bar{3}c] (1) 12 396 1780 10578¾
[P\bar{3}1c, P\bar{3}c1] (4), [R\bar{3}c] (2) 12 396 1540 7218¾
Point group: 6        
P6 6 90 340 1522½
[P6_{1}] (9) 6 90 340 1522½
[P6_{1}] (10) 6 54 91 161½
[P6_{1}] (11) 6 54 97 193½
[P6_{1}] (12) 6 90 280 962½
[P6_{2}] (13) 6 90 340 1522½
[P6_{2}] (14) 6 54 97 193½
[P6_{3}] (3) 6 90 340 1522½
[P6_{3}] (4) 6 90 280 962½
Point group: [\bar{\bf 6}]        
[P\bar{6}] 6 90 310 1242½
Point group: [{\bf 6/{\bi m}}]        
[P6/m] 12 540 3400 26643¾
[P6_{3}/m] (3) 12 540 3400 26643¾
[P6_{3}/m] (4) 12 540 2800 16843¾
Point group: 622        
P622 12 324 1150 5506¼
[P6_{1}22] (9) 12 324 1150 5506¼
[P6_{1}22] (10) 12 252 577 1537¾
[P6_{1}22] (11) 12 252 583 1601¾
[P6_{1}22] (12) 12 324 1090 4746¼
[P6_{2}22] (13) 12 324 1150 5506¼
[P6_{2}22] (14) 12 252 583 1601¾
[P6_{3}22] (3) 12 324 1150 5506¼
[P6_{3}22] (4) 12 324 1090 4746¼
Point group: 6mm        
P6mm 12 396 1930 12818¾
P6cc (3) 12 396 1930 12818¾
P6cc (4) 12 396 1450 6098¾
[P6_{3}cm, P6_{3}mc] (3) 12 396 1930 12818¾
[P6_{3}cm, P6_{3}mc] (4) 12 396 1630 8338¾
Point groups: [\bar{\bf 6}{\bi m}{\bf 2}, \bar{\bf 6}{\bf 2}{\bi m}]        
[P\bar{6}m2, P\bar{6}2m] 12 396 1780 10578¾
[P\bar{6}c2, P\bar{6}2c] (3) 12 396 1780 10578¾
[P\bar{6}c2, P\bar{6}2c] (4) 12 396 1540 7218¾
Point group: 6/mmm        
P6/mmm 24 2376 19300 224328[1\over8]
P6/mcc (3) 24 2376 19300 224328[1\over8]
P6/mcc (4) 24 2376 14500 106728[1\over8]
P6/mcm, P6/mmc (3) 24 2376 19300 224328[1\over8]
[P6/mcm], [P6/mmc] (4) 24 2376 16300 145928[1\over8]
Point group: 23        
P23, [P2_{1}3] 12 276 760 2695¼
I23, [I2_{1}3] 24 2208 12160 86248
F23 48 17664 194560 2759936
Point group: [{\bi m}\bar{\bf 3}]        
[Pm\bar{3}, Pn\bar{3}, Pa3] 24 1800 9400 67703[1\over8]
[Im\bar{3}, Ia\bar{3}] 48 14400 150400 2166500
[Fm\bar{3}] 96 115200 2406400 69328000
[Fd\bar{3}] (1) 96 115200 2406400 69328000
[Fd\bar{3}] (2) 96 96768 1484800 28183680
Point group: 432        
[P432, P4_{2}32] 24 1272 4648 25216[7\over8]
[P4_{1}32] (15) 24 1272 4648 25216[7\over8]
[P4_{1}32] (16) 24 1176 3568 13916[7\over8]
[P4_{1}32] (17) 24 1080 2776 8664[7\over8]
[P4_{1}32] (18) 24 984 2272 6580[7\over8]
I432 48 10176 74368 806940
[I4_{1}32] (15) 48 10176 74368 806940
[I4_{1}32] (17) 48 8640 44416 277276
F432 96 81408 1189888 25822080
[F4_{1}32] (15) 96 81408 1189888 25822080
[F4_{1}32] (18) 96 62976 581632 6738816
Point group: [\bar{\bf 4}{\bf 3}{\bi m}]        
[P\bar{4}3m] 24 1272 5128 32896[7\over8]
[P\bar{4}3n] (1) 24 1272 5128 32896[7\over8]
[P\bar{4}3n] (2) 24 1272 4168 17536[7\over8]
[I\bar{4}3m] 48 10176 82048 1052700
[I\bar{4}3d] (15); (20) 48 10176 82048 1052700
[I\bar{4}3d] (15); (21) 48 10176 66688 561180
[I\bar{4}3d] (17) 48 8640 44416 277276
[F\bar{4}3m] 96 81408 1312768 33686400
[F\bar{4}3c] (15) 96 81408 1312768 33686400
[F\bar{4}3c] (18) 96 81408 1067008 17957760
Point group: [{\bi m}\bar{\bf 3}{\bi m}]        
[Pm\bar{3}m, Pn\bar{3}m] 48 8784 72160 972717[13\over16]
[Pn\bar{3}n, Pm\bar{3}n] (1) 48 8784 72160 972717[13\over16]
[Pn\bar{3}n, Pm\bar{3}n] (2) 48 8784 56800 488877[13\over16]
[Im\bar{3}m] 96 70272 1154560 31126970
[Ia\bar{3}d] (15); (20) 96 70272 1154560 31126970
[Ia\bar{3}d] (15); (21) 96 51840 432640 4497850
[Ia\bar{3}d] (17) 96 70272 908800 15644090
[Fm\bar{3}m] 192 562176 18472960 996063040
[Fm\bar{3}c] (1) 192 562176 18472960 996063040
[Fm\bar{3}c] (2) 192 562176 14540800 500610880
[Fd\bar{3}m] (1) 192 562176 18472960 996063040
[Fd\bar{3}m] (2) 192 414720 7782400 205432640
[Fd\bar{3}c] (1) 192 562176 18472960 996063040
[Fd\bar{3}c] (2) 192 414720 6799360 136619840

Note. hkl subsets: (1) [h + k + l = 2n]; (2) [h + k + l = 2n + 1]; (3) [l = 2n]; (4) [l = 2n + 1]; (5) [2h + l = 2n]; (6) [2h + l = 2n + 1]; (7) [2k + l = 2n]; (8) [2k + l = 2n + 1]; (9) [l = 6n]; (10) [l = 6n + 1, 6n + 5]; (11) [l = 6n + 2, 6n + 4]; (12) [l = 6n + 3]; (13) [l = 3n]; (14) [l = 3n + 1], [3n + 2]; (15) hkl all even; (16) only one index odd; (17) only one index even; (18) hkl all odd; (19) two indices odd; (20) [h + k + l = 4n]; (21) [h + k + l = 4n + 2].
And the enantiomorphous space group.

Table 2.1.7.2| top | pdf |
Closed expressions for [\gamma_{2k}] [equation (2.1.7.11)[link]] for space groups of low symmetry

The normalized moments [\gamma_{2k}] are expressed in terms of [M_{k}], where [M_{k} = {(2k)! \over 2^{k}(k!)^{2}} = {(2k - 1)!! \over k!},]and [l'], which takes on the values 1, 2 or 4 according as the Bravais lattice is of type P, one of the types A, B, C or I, or type F, respectively. The expressions for [\gamma_{2k}] are identical for all the space groups based on a given point group, except Fdd2 and Fddd. The expressions are valid for general reflections and under the restrictions given in the text.

Point group(s)Expression for [\gamma_{2k}]
1 1
[\bar{1}, 2, m] [l'^{k - 1} M_{k}]
[2/m, mm2] [l'^{k - 1} M_{k}^{2}]
mmm [l'^{k - 1} M_{k}^{3}]
222 [{l'^{k - 1} \over 2^{k} (k!)^{2}} \sum\limits_{p=0}^{k} (M_{p}M_{k - p})^{3}[p! (k - p)!]^{2}]

2.1.7.4. Fourier versus Hermite approximations

| top | pdf |

As noted in Section 2.1.8.7[link] below, the Fourier representation of the probability distribution of [|F|] is usually much better than the particular orthogonal-function representation discussed in Section 2.1.7.3[link]. Many, perhaps most, non-ideal centric distributions look like slight distortions of the ideal (Gaussian) distribution and have no resemblance to a cosine function. The empirical observation thus seems paradoxical. The probable explanation has been pointed out by Wilson (1986b[link]). A truncated Fourier series is a best approximation, in the least-squares sense, to the function represented. The particular orthogonal-function approach used in equation (2.1.7.5)[link], on the other hand, is not a least-squares approximation to [p_{\rm c}(|E|)], but is a least-squares approximation to [p_{\rm c}(|E|)\exp(|E|^{2}/4). \eqno(2.1.7.15)]The usual expansions (often known as Gram–Charlier or Edgeworth) thus give great weight to fitting the distribution of the (compararively few) strong reflections, at the expense of a poor fit for the (much more numerous) weak-to-medium ones. Presumably, a similar situation exists for the representation of acentric distributions, but this has not been investigated in detail. Since the centric distributions [p_{\rm c}(|E|)] often look nearly Gaussian, one is led to ask if there is an expansion in orthogonal functions that (i) has the leading term [p_{\rm c}(|E|)] and (ii) is a least-squares (as well as an orthogonal-function)2 fit to [p_{\rm c}(|E|)]. One does exist, based on the orthogonal functions [f_{k} = n(x)He_{k}(2^{1/2}x), \eqno(2.1.7.16)]where [n(x)] is the Gaussian distribution (Myller-Lebedeff, 1907[link]). Unfortunately, no reasonably simple relationship between the coefficients [d_{k}] and readily evaluated properties of [p_{\rm c}(|E|)] has been found, and the Myller-Lebedeff expansion has not, as yet, been applied in crystallography. Although Stuart & Ord (1994[link], p. 112) dismiss it in a three-line footnote, it does have important applications in astronomy (van der Marel & Franx, 1993[link]; Gerhard, 1993[link]).

2.1.8. Non-ideal distributions: the Fourier method

| top | pdf |

The starting point of the method described in the previous section is the central-limit theorem approximation, and the method consists of finding correction factors which result in better approximations to the actual p.d.f. Conceptually, this is equivalent to improving the approximation of the characteristic function [cf. equation (2.1.4.10)[link]] over that which led to the central-limit theorem result.

The method to be described in this section does not depend on any initial approximation and will be shown to utilize the dependence of the exact value of the characteristic function on the space-group symmetry, atomic composition and other factors. This approach has its origin in a simple but ingenious observation by Barakat (1974[link]), who noted that if a random variable has lower and upper bounds then the corresponding p.d.f. can be nonzero only within these bounds and can therefore be expanded in an ordinary Fourier series and set to zero (identically) outside the bounded interval. Barakat's (1974[link]) work dealt with intensity statistics of laser speckle, where sinusoidal waves are involved, as in the present problem. This method was applied by Weiss & Kiefer (1983[link]) to testing the accuracy of a steepest-descents approximation to the exact solution of the problem of random walk, and its first application to crystallographic intensity statistics soon followed (Shmueli et al., 1984[link]). Crystallographic (e.g. Shmueli & Weiss, 1987[link]; Rabinovich et al., 1991a[link],b[link]) and noncrystallographic (Shmueli et al., 1985[link]; Shmueli & Weiss, 1985a[link]; Shmueli, Weiss & Wilson, 1989[link]; Shmueli et al., 1990[link]) symmetry was found to be tractable by this approach, as well as joint conditional p.d.f.'s of several structure factors (Shmueli & Weiss, 1985b[link], 1986[link]; Shmueli, Rabinovich & Weiss, 1989[link]). The Fourier method is illustrated below by deriving the exact counterparts of equations (2.1.7.5)[link] and (2.1.7.6)[link] and specifying them for some simple symmetries. We shall then indicate a method of treating higher symmetries and present results which will suffice for evaluation of Fourier p.d.f.'s of [|E|] for a wide range of space groups.

2.1.8.1. General representations of p.d.f.'s of [|E|] by Fourier series

| top | pdf |

We assume, as before, that (i) the atomic phase factors [\vartheta_{j} = 2\pi {\bf h}^{T}{\bf r}_{j}] [cf. equation (2.1.1.2)[link]] are uniformly distributed on (0–2π) and (ii) the atomic contributions to the structure factor are independent. For a centrosymmetric space group, with the origin chosen at a centre of symmetry, the random variable is the (real) normalized structure factor E and its bounds are [-E_{\rm M}] and [E_{\rm M}], where [E_{\rm M} = {\textstyle\sum\limits_{j = 1}^{N}}n_{j}, \hbox{ with } n_{j} = {{f_{j}}\over{\left({\textstyle\sum_{k = 1}^{N}}\,f_{k}^{2}\right)^{1/2}}}. \eqno(2.1.8.1)]Here, [E_{\rm M}] is the maximum possible value of E and [f_{j}] is the conventional scattering factor of the jth atom, including its temperature factor. The p.d.f., [p(E)], can be nonzero in the range ([-E_{\rm M},E_{\rm M}]) only and can thus be expanded in the Fourier series [p(E) = (\alpha/2)\textstyle\sum\limits_{k = -\infty}^{\infty}C_{k}\exp(-\pi ik\alpha E), \eqno(2.1.8.2)]where [\alpha = 1/E_{\rm M}]. Only the real part of [p(E)] is relevant. The Fourier coefficients can be obtained in the conventional manner by integrating over the range [(-E_{\rm M},E_{\rm M})],[C_{k} = \textstyle\int\limits_{-E_{\rm M}}^{E_{\rm M}}p(E)\exp(\pi ik\alpha E)\,{\rm d}E. \eqno(2.1.8.3)]Since, however, [p(E) = 0] for [E \,\lt\, {-}E_{\rm M}] and [E \,\gt\, E_{\rm M}], it is possible and convenient to replace the limits of integration in equation (2.1.8.3)[link] by infinity. Thus [C_{k} = \textstyle\int\limits_{-\infty}^{\infty}p(E)\exp(\pi ik\alpha E)\,{\rm d}E = \langle \exp(\pi ik\alpha E) \rangle. \eqno(2.1.8.4)]Equation (2.1.8.4)[link] shows that [C_{k}] is a Fourier transform of the p.d.f. [p(E)] and, as such, it is the value of the corresponding characteristic function at the point [t_{k} = \pi\alpha k] [i.e., [C_{k} = C(\pi\alpha k)], where the characteristic function [C(t)] is defined by equation (2.1.4.1)[link]]. It is also seen that [C_{k}] is the expected value of the exponential [\exp(\pi ik\alpha E)]. It follows that the feasibility of the present approach depends on one's ability to evaluate the characteristic function in closed form without the knowledge of the p.d.f.; this is analogous to the problem of evaluating absolute moments of the structure factor for the correction-factor approach, discussed in Section 2.1.7[link]. Fortunately, in crystallographic applications these calculations are feasible, provided individual isotropic motion is assumed. The formal expression for the p.d.f. of [|E|], for any centrosymmetric space group, is therefore [p(|E|) = \alpha\left[1+2\textstyle\sum\limits_{k = 1}^{\infty}C_{k}\cos(\pi k\alpha |E|) \right], \eqno(2.1.8.5)]where use is made of the assumption that [p(E) = p(-E)], and the Fourier coefficients are evaluated from equation (2.1.8.4)[link].

The p.d.f. of [|E|] for a noncentrosymmetric space group is obtained by first deriving the joint p.d.f. of the real and imaginary parts of E and then integrating out its phase. The general expression for E is [E = A + iB = |E|\cos\varphi + i|E|\sin\varphi, \eqno(2.1.8.6)]where [\varphi] is the phase of E. The required joint p.d.f. is [p(A,B) = (\alpha^{2}/4)\textstyle\sum\limits_{m}\textstyle\sum\limits_{n}C_{mn}\exp[-\pi i\alpha(mA+nB)], \eqno(2.1.8.7)]and introducing polar coordinates [m = r\sin\Delta] and [n = r\cos\Delta], where [r = \sqrt{m^{2}+n^{2}}] and [\Delta = \tan^{-1}(m/n)], we have [\eqalignno{p(|E|,\varphi) &= (\alpha^{2}/4)|E|\textstyle\sum\limits_{m}\textstyle\sum\limits_{n}C_{mn}\exp[-\pi i \alpha |E| & \cr &\quad\times{}\sqrt{m^{2}+n^{2}}\sin(\varphi + \Delta)]. &(2.1.8.8) \cr}]Integrating out the phase [\varphi], we obtain [p(|E|) = (\pi\alpha^{2} |E|/2)\textstyle\sum\limits_{m}\textstyle\sum\limits_{n}C_{mn}J_{0}(\pi\alpha |E|\sqrt{m^{2}+n^{2}}), \eqno(2.1.8.9)]where [J_{0}(x)] is the Bessel function of the first kind (e.g. Abramowitz & Stegun, 1972[link]). This is a general form of the p.d.f. of [|E|] for a noncentrosymmetric space group. The Fourier coefficients are obtained, similarly to the above, as [C_{mn} = \langle \exp[\pi i\alpha(mA+nB)] \rangle \eqno(2.1.8.10)]and the average in equation (2.1.8.10)[link], just as that in equation (2.1.8.4)[link], is evaluated in terms of integrals over the appropriate trigonometric structure factors. In terms of the characteristic function for a joint p.d.f. of A and B, the Fourier coefficient in equation (2.1.8.10)[link] is given by [C_{mn} = C(\pi\alpha m,\pi\alpha n)].

We shall denote the characteristic function by [C(t_{1})] if it corresponds to a Fourier coefficient of a Fourier series for a centrosymmetric space group and by [C(t_{1},t_{2})] or by [C(t,\Delta)], where [t = (t_{1}^{2}+t_{2}^{2})^ {1/2}] and [\Delta = \tan^{-1}(t_{1}/t_{2})], if it corresponds to a Fourier series for a noncentrosymmetric space group.

2.1.8.2. Fourier–Bessel series

| top | pdf |

Equations (2.1.8.5)[link] and (2.1.8.9)[link] are the exact counterparts of equations (2.1.7.5)[link] and (2.1.7.6)[link], respectively. The computational effort required to evaluate equation (2.1.8.9)[link] is somewhat greater than that for (2.1.8.5)[link], because a double Fourier series has to be summed. The p.d.f. for any noncentrosymmetric space group can be expressed by a double Fourier series, but this can be simplified if the characteristic function depends on [t = (t_{1}^{2}+t_{2}^{2})^{1/2}] alone, rather than on [t_{1}] and [t_{2}] separately. In such cases the p.d.f. of [|E|] for a noncentrosymmetric space group can be expanded in a single Fourier–Bessel series (Barakat, 1974[link]; Weiss & Kiefer, 1983[link]; Shmueli et al., 1984[link]). The general form of this expansion is [p(|E|) = 2\alpha^{2}|E|\textstyle\sum\limits_{u = 1}^{\infty}D_{u}J_{0}(\alpha\lambda_{u}|E|), \eqno(2.1.8.11)]where [D_{u} = {{C(\alpha\lambda_{u})}\over{J_{1}^{2}(\lambda_{u})}} \eqno(2.1.8.12)]and [C(\alpha\lambda_{u}) = \textstyle\prod\limits_{j = 1}^{N/g}C_{ju}, \eqno(2.1.8.13)]where [J_{1}(x)] is the Bessel function of the first kind, and [\lambda_{u}] is the uth root of the equation [J_{0}(x) = 0]; the atomic contribution [C_{ju}] to equation (2.1.8.13)[link] is computed as [C_{ju} = C(\alpha n_{j}\lambda_{u}). \eqno(2.1.8.14)]The roots [\lambda_{u}] are tabulated in the literature (e.g. Abramowitz & Stegun, 1972[link]), but can be most conveniently computed as follows. The first five roots are given by [\eqalign{ \lambda_1& = 2.4048255577\cr \lambda_2& = 5.5200781103 \cr \lambda_3& = 8.6537279129\cr \lambda_4& = 11.7915344390\cr \lambda_5& = 14.9309177085 }]and the higher ones can be obtained from McMahon's approximation (cf. Abramowitz & Stegun, 1972[link]) [{\lambda_{u} = \beta+{{1}\over{8\beta}}-{{124}\over{3(8\beta)^{3}}}+{{120928}\over{15(8\lambda)^{5}}}-{{401743168}\over{105(8\lambda)^{7}}}+\ldots,} \eqno(2.1.8.15)]where [\beta = (u-{\textstyle {1\over4}})\pi]. For [u \,\gt\, 5] the values given by equation (2.1.8.15)[link] have a relative error less than 10−11 so that no refinement of roots of higher orders is needed (Shmueli et al., 1984[link]). Numerical computations of single Fourier–Bessel series are of course faster than those of the double Fourier series, but both representations converge fairly rapidly.

2.1.8.3. Simple examples

| top | pdf |

Consider the Fourier coefficient of the p.d.f. of [|E|] for the centrosymmetric space group [P\bar{1}]. The normalized structure factor is given by [E = 2\textstyle\sum\limits_{j = 1}^{N/2}n_{j}\cos\vartheta_{j}, \quad{\rm with}\quad \vartheta_{j} = 2\pi {\bf h}^{T}\cdot{\bf r}_{j}, \eqno(2.1.8.16)]and the Fourier coefficient is [\eqalignno{C_{k} & = \langle \exp(\pi ik\alpha E) \rangle &(2.1.8.17) \cr& = \left\langle \exp\left[2\pi ik\alpha\textstyle\sum\limits_{j = 1}^{N/2}n_{j}\cos\vartheta_{j}\right] \right\rangle &(2.1.8.18) \cr & = \left\langle \textstyle\prod\limits_{j = 1}^{N/2} \exp(2\pi ik\alpha n_{j}\cos\vartheta_{j}) \right\rangle &(2.1.8.19) \cr & = \textstyle \prod\limits_{j = 1}^{N/2} \langle \exp(2\pi ik\alpha n_{j}\cos\vartheta_{j}) \rangle &(2.1.8.20) \cr & = \textstyle\prod\limits_{j = 1}^{N/2}\left\{{({1}/{2\pi})}\int\limits_{-\pi}^{\pi}\exp(2\pi ik\alpha n_{j}\cos\vartheta)\,{\rm d}\vartheta \right\}\cr &&(2.1.8.21) \cr & = \textstyle\prod\limits_{j = 1}^{N/2}J_{0}(2\pi k\alpha n_{j}). &(2.1.8.22)}%2.1.8.22]Equation (2.1.8.20)[link] is obtained from equation (2.1.8.19)[link] if we make use of the assumption of independence, the assumption of uniformity allows us to rewrite equation (2.1.8.20)[link] as (2.1.8.21)[link], and the expression in the braces in the latter equation is just a definition of the Bessel function [J_{0}(2\pi k\alpha n_{j})] (e.g. Abramowitz & Stegun, 1972[link]).

Let us now consider the Fourier coefficient of the p.d.f. of [|E|] for the noncentrosymmetric space group [P1]. We have [A = \textstyle\sum\limits_{j = 1}^{N}n_{j}\cos\vartheta_{j} \quad {\rm and}\quad B = \textstyle\sum\limits_{j = 1}^{N}n_{j}\sin\vartheta_{j}. \eqno(2.1.8.23)]These expressions for A and B are substituted in equation (2.1.8.10)[link], resulting in [\eqalignno{C_{mn} & = \left\langle \textstyle\prod\limits_{j = 1}^{N} \exp[\pi i\alpha n_{j}(m\cos\vartheta_{j} + n\sin\vartheta_{j})] \right\rangle & \cr & & (2.1.8.24) \cr & = \left\langle \textstyle\prod\limits_{j = 1}^{N} \exp[\pi i\alpha n_{j}\sqrt{m^{2}+n^{2}} \sin(\vartheta_{j}+\Delta)] \right\rangle & \cr & & (2.1.8.25)\cr &= \textstyle \prod\limits_{j = 1}^{N} J_{0}(\pi\alpha n_{j}\sqrt{m^{2}+n^{2}}). &(2.1.8.26)}]Equation (2.1.8.24)[link] leads to (2.1.8.25)[link] by introducing polar coordinates analogous to those leading to equation (2.1.8.8)[link], and equation (2.1.8.26)[link] is then obtained by making use of the assumptions of independence and uniformity in an analogous manner to that detailed in equations (2.1.8.12)[link]–(2.1.8.22)[link] above.

The right-hand side of equation (2.1.8.26)[link] is to be used as a Fourier coefficient of the double Fourier series given by (2.1.8.9)[link]. Since, however, this coefficient depends on [(m^{2}+n^{2})^{1/2}] alone rather than on m and n separately, the p.d.f. of [|E|] for P1 can also be represented by a Fourier–Bessel series [cf. equation (2.1.8.11)[link]] with coefficient [D_{u} = {{1}\over{J_{1}^{2}(\lambda_{u})}}\prod_{j = 1}^{N}J_{0}(\alpha n_{j}\lambda_ {u}), \eqno(2.1.8.27)]where [\lambda_{u}] is the uth root of the equation [J_{0}(x) = 0].

2.1.8.4. A more complicated example

| top | pdf |

We now illustrate the methodology of deriving characteristic functions for space groups of higher symmetries, following the method of Rabinovich et al. (1991a[link],b[link]). The derivation is performed for the space group [P\bar{6}] (No. 174). According to Table A1.4.3.6[link] , the real and imaginary parts of the normalized structure factor are given by [\eqalignno{A & = 2\textstyle\sum\limits_{j = 1}^{N/6}n_{j}[C(hki)c(lz)]_{j} \cr & = 2\textstyle\sum\limits_{j = 1}^{N/6}n_{j}\cos\tau_{j}\textstyle\sum\limits_{k = 1}^{3}\cos\alpha_{jk} &(2.1.8.28)}]and [\eqalignno{B & = 2\textstyle\sum\limits_{j = 1}^{N/6}n_{j}[S(hki)c(lz)]_{j} \cr & = 2\textstyle\sum\limits_{j = 1}^{N/6}n_{j}\cos\tau_{j}\textstyle\sum\limits_{k = 1}^{3}\sin\alpha_{jk}, &(2.1.8.29)}]where[\eqalign{ \alpha_{j1}& = 2\pi (hx_j+ky_j),\cr \alpha_{j2}& = 2\pi (kx_j+iy_j),\cr \alpha_{j3}& = 2\pi (ix_j+hy_j),\cr \tau_j& = 2\pi lz_j. }]Note that [\alpha_{j1}+\alpha_{j2}+\alpha_{j3} = 0], i.e., one of these contributions depends on the other two; this is a recurring problem in calculations pertaining to trigonal and hexagonal systems. For brevity, we write directly the general form of the characteristic function from which the functional form of the Fourier coefficient can be readily obtained. The characteristic function is given by [\eqalignno{C(t_{1},t_{2})& = \langle \exp[i(t_{1}A+t_{2}B)] \rangle &(2.1.8.30) \cr & = \textstyle\prod\limits_{j = 1}^{N/6}\left\langle \exp\left[2in_{j}\cos\tau_{j}\textstyle\sum\limits_{k = 1}^{3}(t_{1} \cos\alpha_{jk}+t_{2}\sin\alpha_{jk})\right]\right\rangle & \cr & & (2.1.8.31) \cr & = \textstyle \prod\limits_{j = 1}^{N/6}\biggl\langle \exp\biggl[2in_{j}t\cos\tau_{j}\textstyle\sum\limits_{k = 1}^{3}(\sin \Delta\cos\alpha_{jk} & \cr & \quad+ \cos\Delta\sin\alpha_{jk})\biggr]\biggr\rangle & (2.1.8.32) \cr & = \textstyle\prod\limits_{j = 1}^{N/6}\left\langle \exp\left[2in_{j}t\cos\tau_{j}\textstyle\sum\limits_{k = 1}^{3} \sin(\alpha_{jk}+\Delta)\right]\right\rangle , & \cr & &(2.1.8.33)}%2.1.8.33]where [\Delta = \tan^{-1}(t_{1}/t_{2})], [t = ({t_{1}^{2}+t_{2}^{2}}){}^{1/2}] and the assumption of independence was used. If we further employ the assumption of uniformity, while remembering that the angular variables [\alpha_{jk}] are not independent, the characteristic function can be written as [\eqalignno{C(t_{1},t_{2}) & = \textstyle\prod\limits_{j = 1}^{N/6}\left\{(1/2\pi)\textstyle\int\limits_{-\pi}^{\pi}\,{\rm d}\tau \biggl([1/(2\pi)^{2}] \right. \cr &\quad \times\textstyle\int\limits_{-\pi}^{\pi}\textstyle\int\limits_{-\pi}^{\pi}\textstyle\int\limits_{-\pi}^{\pi} {\rm d}\alpha_{1}\,{\rm d}\alpha_{2}\,{\rm d}\alpha_{3}\delta_{2\pi}(\alpha_{1}+\alpha_{2}+ \alpha_{3}) \cr &\quad\times \left.\left. \exp\left[2in_{j}t\cos\tau\textstyle\sum\limits_{k = 1}^{3}\sin(\alpha_{k} +\Delta)\right] \right) \right\}, & \cr & &(2.1.8.34)}]where [\delta_{2\pi}(\alpha) = {{1}\over{2\pi}}\sum_{k = -\infty}^{\infty}\exp(-ik\alpha) \eqno(2.1.8.35)]is the Fourier representation of the periodic delta function. Equation (2.1.8.34)[link] then becomes [\eqalignno{C(t_{1},t_{2}) &= \textstyle\prod\limits_{j = 1}^{N/6}\biggl\{(1/2\pi)\textstyle\int\limits_{-\pi}^{\pi}\,{\rm d}\tau \textstyle\sum\limits_{k = -\infty}^{\infty}\biggl[(1/2\pi) & \cr &\quad\times \textstyle\int\limits_{-\pi}^{\pi}\exp\left(-ik \alpha+2in_{j}t\cos\tau\sin(\alpha+\Delta)\right)\,{\rm d}\alpha\biggr]^{3}\biggr\}. & \cr & &(2.1.8.36) \cr}]If we change the variable [\alpha] to [\alpha'-\Delta], [\sin(\alpha+\Delta)] becomes [\sin\alpha'] and [-ik\alpha = -ik\alpha'+ik\Delta]. Hence [C(t_{1},t_{2}) = \textstyle\prod\limits_{j = 1}^{N/6}\biggl\{(1/2\pi)\textstyle\int\limits_{-\pi}^{\pi}\,{\rm d}\tau\textstyle\sum\limits_{k = -\infty}^{\infty}\exp(3ik\Delta) J_{k}^{3}(2n_{j}t\cos\tau)\biggr\}. \eqno(2.1.8.37)]The imaginary part of the summation, involving Bessel functions of odd orders, vanishes upon integration and the latter is restricted to the positive quadrant in [\tau]. Thus, upon replacing cosines by sines (this is permissible at this stage) the atomic contribution to the characteristic function becomes [\eqalignno{C_{j}(t,\Delta) &= (2/\pi)\textstyle\int\limits_{0}^{\pi/2}\biggl[J_{0}^{3}(2n_{j}t \sin\tau) & \cr &\quad + 2\textstyle\sum\limits_{k = 1}^{\infty}\cos(6k\Delta)J_{k}^{3}(2n_{j}t\sin\tau) \biggr]\,{\rm d}\tau & \cr & & (2.1.8.38) \cr}]and a double Fourier series must be used for the p.d.f.

2.1.8.5. Atomic characteristic functions

| top | pdf |

Expressions for the atomic contributions to the characteristic functions were obtained by Rabinovich et al. (1991a[link]) for a wide range of space groups, by methods similar to those described above. These expressions are collected in Table 2.1.8.1[link] in terms of symbols which are defined below. The following abbreviations are used in the subsequent definitions of the symbols: [\eqalign{s_{\pm} &= 2an_{j}\sin(\tau\pm\rho), \cr c_{\pm} &= 2an_{j}\cos(\tau\pm \rho)\quad {\rm and}\cr \sigma_{\pm} &= 2an_{j}\sin(\tau\pm 2\pi/3 + \rho),} ]and the symbols appearing in Table 2.1.8.1[link] are given below:[\eqalign{{}^{(a)}L_j(a,\rho) & = \langle J_0(s_+)J_0(s_-) \rangle_\tau \cr & = \textstyle\sum\limits_{k = -\infty}^\infty\cos(4k\rho)J_k^4 \cr & = J_0^4(an_j) + 2\textstyle\sum\limits_{k = 1}^\infty\cos(4k\rho)J_k^4(an_j),\cr {}^{(b)}Q_j^{(1)} (a,\rho) & = \langle J_0^2(s_+)J_0^2(s_-) \rangle_\tau,\cr {}^{(c)}Q_j^{(2)} (a,\rho) & = \langle J_0(s_+)J_0(s_-)J_0(c_+)J_0 (c_-) \rangle_\tau,\cr {}^{(d)}T_j(a,\rho) & = \textstyle\sum\limits_{k = -\infty}^\infty\exp(6ik\rho)J_k^6(an_j) \cr & = J_0^6(an_j)+2\textstyle\sum\limits_{k = 1}^\infty\cos(6k\rho)J_k^6(an_j),\cr {}^{(e)}H_j^{(1)} (a,\mu) & = \left\langle {\cal R} \left [S_j^{(1)} (\tau\semi\ a,\mu,0) \right] \right\rangle _\tau, \cr {}^{(f)}H_j^{(2)} (a,\mu) & = \left\langle {\cal R} \left [S_j^{(2)} (\tau\semi \ a,\mu,0) \right] \right\rangle _\tau,\cr {}^{(g)}\tilde{H}_j^{(1)} (a,\mu_1,\mu_2,\rho) & = \biggl\langle {{\cal R}} \biggl[S_j^{(1)} (\tau\semi \ a,\mu_1,\rho) \cr &\quad\times{}S_j^{(1)} (\tau\semi \ a,\mu_2,-\rho) \biggr] \biggr\rangle _\tau, \cr {}^{(h)}\tilde{H}_j^{(2)} (a,\mu_1,\mu_2,\rho) & = \biggl\langle {\cal R} \biggl[S_j^{(2)} (\tau\semi \ a,\mu_1,\rho) \cr &\quad\times{}S_j^{(2)} (\tau\semi \ a,\mu_2,-\rho) \biggr] \biggr\rangle _\tau, }]where [S_{j}^{(1)}(\tau\semi\ a,\mu,\rho) = \textstyle\sum\limits_{k = -\infty}^{\infty}e^{3ik\mu}J_{k}^{3} (s_{+})]and [S_{j}^{(2)}(\tau\semi\ a,\mu,\rho) = \textstyle\sum\limits_{k = -\infty}^{\infty}e^{3ik\mu}J_{k} (s_{+})J_{k}(\sigma_{+})J_{k}(\sigma_{-}).]

Table 2.1.8.1| top | pdf |
Atomic contributions to characteristic functions for [p(|E|)]

The table lists symbolic expressions for the atomic contributions to exact characteristic functions (abbreviated as c.f.) for [p(|E|)], to be computed as single Fourier series (centric), double Fourier series (acentric) and single Fourier–Bessel series (acentric), as defined in Sections 2.1.8.1[link] and 2.1.8.2[link]. The symbolic expressions are defined in Section 2.1.8.5[link]. The table is arranged by point groups, space groups and parities of the reflection indices analogously to the table of moments, Table 2.1.7.1[link], and covers all the space groups and statistically different parities of hkl up to and including space group [Fd\bar{3}]. The expressions are valid for atoms in general positions, for general reflections and presume the absence of noncrystallographic symmetry and of dispersive scatterers. The p.d.f.'s of [|E|] for any centrosymmetric space group are computed from a single Fourier series [cf. (2.1.8.5)[link]]. For a noncentrosymmetric space group the p.d.f. of [|E|] must be computed from a double Fourier series [cf. (2.1.8.9)[link]] – if the characteristic function in the table depends on Δ – and may be computed from a single Fourier–Bessel series [cf. (2.1.8.11)[link]] if it does not depend on Δ. However, the p.d.f. of [|E|] for any noncentrosymmetric space group may also be computed from a double Fourier series (cf. Section 2.1.8.1[link]).

Space group(s)gAtomic c.f.Remarks
Point group: 1      
P1 1 [ J_{0}(t n_{j})]  
Point group: [\bar{\bf 1}]      
[P\bar{1}] 2 [J_{0}(2t_{1}n_{j})]  
Point groups: 2, m      
All P 2 [J_{0}^{2}(t n_{j})]  
All C 4 [J_{0}^{2}(2t n_{j})]  
Point group: [{\bf 2/}{\bi m}]      
All P 4 [J_{0}^{2}(2t_{1} n_{j})]  
All C 8 [J_{0}^{2}(4t_{1} n_{j})]  
Point group: 222      
All P 4 [ L_{j}(t,\Delta)^{(a)}]  
All C and I 8 [ L_{j}(2t,\Delta)]  
F222 16 [ L_{j}(4t,\Delta)]  
Point group: mm2      
All P 4 [ L_{j}(t,0)]  
All C and I 8 [ L_{j}(2t,0)]  
Fmm2 16 [ L_{j}(4t,0)]  
Fdd2 16 [L_{j}(4t,0)] [ h+k+l=2n]
  16 [ L_{j}(4t,\pi/4)] [ h+k+l=2n+1]
Point group: mmm      
All P 8 [ L_{j}(2t_{1},0)]  
All C and I 16 [ L_{j}(4t_{1},0)]  
Fmmm 32 [ L_{j}(8t_{1},0)]  
Fddd 32 [ L_{j}(8t_{1},0)] [ h+k+l=2n]
  32 [ L_{j}(8t_{1},\pi/4)] [ h+k+l=2n+1]
Point group: 4      
P4, P42 4 [ L_{j}(t,0)]  
[P4_1] 4 [ L_{j}(t,0)] [ l=2n]
  4 [ L_{j}(t,\pi/4)] [ l=2n+1]
I4 8 [ L_{j}(2t,0)]  
[ I4_{1}] 8 [ L_{j}(2t,0)] [2h+l=2n]
  8 [ L_{j}(2t,\pi/4)] [ 2h+l=2n+1]
Point group: [\bar{\bf 4}]      
[P\bar{4}] 4 [L_{j}(t,\Delta)]  
[I\bar{4}] 8 [ L_{j}(2t,\Delta)]  
Point group: [{\bf 4/}{\bi m}]      
All P 8 [ L_{j}(2t_{1},0)]  
[ I4/m] 16 [ L_{j}(4t_{1},0)]  
[ I4_1/a] 16 [ L_{j}(4t_{1},0)] [ l=2n]
  16 [L_{j}(4t_{1},\pi/4)] [ l=2n+1]
Point group: 422      
P422, P4212, P4222, P42212 8 [ Q_{j}^{(1)}(t,\Delta)^{(b)}]  
P4122, P41212 8 [ Q_{j}^{(1)}(t,\Delta)] [ l=2n]
  8 [ Q_{j}^{(2)}(t,\Delta)^{(c)}] [ l=2n+1]
I422 16 [ Q_{j}^{(1)}(2t,\Delta)]  
I4122 16 [ Q_{j}^{(1)}(2t,\Delta)] [ 2k+l=2n]
  16 [ Q_{j}^{(2)}(2t,\Delta)] [ 2k+l=2n+1]
Point group: 4mm      
All P 8 [ Q_{j}^{(1)}(t,0)]  
I4mm, I4cm 16 [ Q_{j}^{(1)}(2t,0)]  
I41md, I41cd 16 [Q_{j}^{(1)}(2t,0)] [2k+l=2n]
  16 [ Q_{j}^{(1)}(2t,\pi/4)] [2k+l=2n+1]
Point groups: [\bar{\bf 4}]2m, [\bar{\bf 4}]m2      
All P 8 [ Q_{j}^{(1)}(t,\Delta)]  
I[\bar{4}]2m, I[\bar{4}]m2, I[\bar{4}]c2 16 [Q_{j}^{(1)}(2t,\Delta)]  
I[\bar{4}]2d 16 [Q_{j}^{(1)}(2t,\Delta)] [2h+l=2n]
  16 [Q_{j}^{(2)}(2t,\Delta)] [2h+l=2n+1]
Point group: [{\bf 4/}{\bi mmm}]      
All P 16 [ Q_{j}^{(1)}(2t_{1},0)]  
[I4/mmm], [I4/mcm] 32 [Q_{j}^{(1)}(4t_{1},0)]  
[ I4_{1}/amd], [I4_{1}/acd] 32 [Q_{j}^{(1)}(4t_{1},0)] [ l=2n]
  32 [ Q_{j}^{(1)}(4t_{1},\pi/4)] [ l=2n+1]
Point group: 3      
All P and R 3 [ J_{0}^{3}(t n_{j})]  
Point group: [\bar{\bf 3}]      
All P and R 6 [ J_{0}^{3}(2t_{1} n_{j})]  
Point group: 32      
All P and R 6 [T_{j}(t,\Delta)^{(d)}]  
Point group: 3m      
P3m1, P31m, R3m 6 [ T_{j}(t,\pi/2)]  
P3c1, P31c, R3c 6 [ T_{j}(t,\pi/2)] [ l=2n\,(P)], [h+k+l= 2n\,(R)]
  6 [ T_{j}(t,0)] [ l=2n+1\,(P)], [h+k+l =2n+1\,(R)]
Point group: [\bar{\bf 3}]m      
P[\bar{3}]m1, P[\bar{3}]1m, R[\bar{3}]m 12 [ T_{j}(2t_{1},\pi/2)]  
P[\bar{3}]c1, P[\bar{3}]1c, R[\bar{3}]c 12 [ T_{j}(2t_{1},\pi/2)] [ l=2n\,(P)], [h+k+l=2n\,(R)]
  12 [ T_{j}(2t_{1},0)] [ l=2n+1\,(P)], [h+k+l=2n+1\,(R)]
Point group: 6      
P6 6 [H_{j}^{(1)}(t, \pi/2)^{(e)}]  
[ P6_{1}] 6 [H_{j}^{(1)}(t, \pi/2)] [ l=6n]
  6 [ H_{j}^{(2)}(t, 0)^{(f)}] [ l=6n+1], [6n+5]
  6 [H_{j}^{(2)}(t, \pi/2)] [l=6n+2], [6n+4]
  6 [ H_{j}^{(1)}(t, 0)] [ l=6n+3]
[ P6_{2}] 6 [H_{j}^{(1)}(t, \pi/2)] [ l=3n]
  6 [H_{j}^{(2)}(t, \pi/2)] [l=3n \pm 1]
[ P6_{3}] 6 [ H_{j}^{(1)}(t, \pi/2)] [ l=2n]
  6 [ H_{j}^{(1)}(t, 0)] [ l=2n+1]
Point group: [\bar{\bf 6}]      
[P\bar{6}] 6 [H_{j}^{(1)}(t, \Delta)]  
Point group: [{\bf 6/}{\bi m}]      
[ P6/m] 12 [H_{j}^{(1)}(2t_{1}, \pi/2)]  
[ P6_{3}/m] 12 [ H_{j}^{(1)}(2t_{1}, \pi/2)] [ l=2n]
  12 [ H_{j}^{(1)}(2t_{1}, 0)] [l=2n+1]
Point group: 622      
P622 12 [\tilde{H}_{j}^{(1)} (t,\pi/2], [-\pi/2,\Delta)^{(g)}]  
[ P6_{1}22] 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2], [-\pi/2,\Delta)] [l=6n]
  12 [ \tilde{H}_{j}^{(2)}(t,0,0,\Delta)^{(h)}] [ l=6n+1], [6n+5]
  12 [ \tilde{H}_{j}^{(2)}(t,\pi/2,\pi/2,\Delta)] [ l=6n+2], [6n+4]
  12 [ \tilde{H}_{j}^{(1)}(t,0,0,\Delta)] [ l=6n+3]
[ P6_{2}22] 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2], [-\pi/2,\Delta)] [ l=3n]
  12 [ \tilde{H}_{j}^{(2)}(t,\pi/2,\pi/2,\Delta)] [ l=3n \pm 1]
[ P6_{3}22] 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2], [-\pi/2,\Delta)] [ l=2n]
  12 [ \tilde{H}_{j}^{(1)}(t,0,0,\Delta)] [ l=2n+1]
Point group: 6mm      
P6mm 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2,\pi/2,0)]  
P6cc 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2,\pi/2,0)] [ l=2n]
  12 [ \tilde{H}_{j}^{(1)}(t,\pi/2], [-\pi/2,0)] [ l=2n+1]
[P6_{3}cm], [P6_{3}mc] 12 [ \tilde{H}_{j}^{(1)}(t,\pi/2,\pi/2,0)] [l=2n]
  12 [ \tilde{H}_{j}^{(1)}(t,0,0,0)] [l=2n+1]
Point groups: [\bar{\bf 6}]2m, [\bar{\bf 6}]m2      
P[\bar{6}]2m, P[\bar{6}]m2 12 [ \tilde{H}_{j}^{(1)}(t,\Delta,\Delta,0)]  
P[\bar{6}]2c, P[\bar{6}]c2 12 [ \tilde{H}_{j}^{(1)}(t,\Delta,\Delta,0)] [ l=2n]
  12 [ \tilde{H}_{j}^{(1)}(t,\Delta+\pi/2], [-\Delta-\pi/2,0)] [ l=2n+1]
Point group: [{\bf 6/}{\bi mmm}]      
[P6/mmm] 24 [\tilde{H}_{j}^{(1)}(2t_{1},\pi/2,\pi/2,0)]  
[ P6/mcc] 24 [\tilde{H}_{j}^{(1)}(2t_{1},\pi/2,\pi/2,0)] [ l=2n]
  24 [ \tilde{H}_{j}^{(1)}(2t_{1}], [\pi/2], [-\pi/2,0)] [ l=2n+1]
[P6_{3}/mcm], [P6_{3}/mmc] 24 [ \tilde{H}_{j}^{(1)}(2t_{1},\pi/2,\pi/2,0)] [l=2n]
  24 [ \tilde{H}_{j}^{(1)}(2t_{1},0,0,0)] [ l=2n+1]
Point group: 23      
P23, P213 12 [ L_{j}^{3}(t,\Delta)]  
I23, [I2_{1}3] 24 [ L_{j}^{3}(2t,\Delta)]  
F23 48 [ L_{j}^{3}(4t,\Delta)]  
Point group: m[\bar{\bf 3}]      
Pm[\bar{3}], Pn[\bar{3}], Pa[\bar{3}] 24 [ L_{j}^{3}(2t_{1},0)]  
Im[\bar{3}], Ia[\bar{3}] 48 [ L_{j}^{3}(4t_{1},0)]  
Fm[\bar{3}] 96 [ L_{j}^{3}(8t_{1},0)]  
Fd[\bar{3}] 96 [ L_{j}^{3}(8t_{1},0)] [h+k+l=2n]
  96 [ L_{j}^{3}(8t_{1},\pi/4)] [h+k+l=2n+1]
And the enantiomorphous space group.

The averages appearing in the above summary are, in general, computed as[\left\langle f(\tau)\right\rangle = ({{2}/{\pi}})\textstyle\int\limits_{0}^{\pi/2}f(\tau)\,{\rm d}\tau, \eqno(2.1.8.39)]except [H_{j}^{(2)}] and [\tilde{H}_{j}^{(2)}] which are computed as [\left\langle f(\tau)\right\rangle = ({{3}/{\pi}})\textstyle\int\limits_{0}^{\pi/3}f(\tau)\,{\rm d}\tau, \eqno(2.1.8.40)]where [f(\tau)] is any of the atomic characteristic functions indicated above. The superscripts preceding the symbols in the above summary are appended to the corresponding symbols in Table 2.1.8.1[link] on their first occurrence.

2.1.8.6. Other non-ideal Fourier p.d.f.'s

| top | pdf |

As pointed out above, the representation of the p.d.f.'s by Fourier series is also applicable to effects of noncrystallographic symmetry. Thus, Shmueli et al. (1985[link]) obtained the following Fourier coefficient for the bicentric distribution in the space group [P\bar{1}]:[C_{k} = (2/\pi)\textstyle\int\limits_{0}^{\pi/2}\left[\textstyle\prod\limits_{j = 1}^{N/4}J_{0}(4\pi k\alpha n_{j}\cos\vartheta)\right]\,{\rm d}\vartheta, \eqno(2.1.8.41)]to be used with equation (2.1.8.5)[link]. Furthermore, if we use the important property of the characteristic function as outlined in Section 2.1.4.1[link], it is easy to write down the Fourier coefficient for a [P\bar{1}] asymmetric unit containing a centrosymmetric fragment centred at a noncrystallographic centre and a number of atoms not related by symmetry. This Fourier for the above partially bicentric arrangement is a product of expressions (2.1.8.17)[link] and (2.1.8.41)[link], with the appropriate number of atoms in each factor (Shmueli & Weiss, 1985a[link]). While the purely bicentric p.d.f. obtained by using (2.1.8.41)[link] with (2.1.8.5)[link] is significantly different from the ideal bicentric p.d.f. given by equation (2.1.5.13)[link] only when the atomic composition is sufficiently heterogeneous, the above partially bicentric p.d.f. appears to be a useful development even for an equal-atom structure.

The problem of the coexistence of several noncrystallographic centres of symmetry within the asymmetric unit of [P\bar{1}], and its effect on the p.d.f. of [|E|], was examined by Shmueli, Weiss & Wilson (1989[link]) by the Fourier method. The latter study indicates that the strongest effect is produced by the presence of a single noncrystallographic centre.

Another kind of noncrystallographic symmetry is that arising from the presence of centrosymmetric fragments in a noncentrosymmetric structure – the subcentric arrangement already discussed in Section 2.1.5.4[link]. A Fourier-series representation of a non-ideal p.d.f. corresponding to this case was developed by Shmueli, Rabinovich & Weiss (1989[link]), and was also applied to the mathematically equivalent effects of dispersion and presence of heavy scatterers in centrosymmetric special positions in a noncentrosymmetric space group.

A variety of other non-ideal p.d.f.'s occur when heavy atoms are present in special positions (Shmueli & Weiss, 1988[link]). Without going into the details of this development, it can be noted that if the atoms are distributed among k types of Wyckoff positions, the characteristic function corresponding to the p.d.f. of [|E|] is a product of the k characteristic functions, each of which is related to one of these special positions; the same property of the characteristic function as that in Section 2.1.4.1[link] is here utilized.

2.1.8.7. Comparison of the correction-factor and Fourier approaches

| top | pdf |

The need for theoretical non-ideal distributions was exemplified by Fig. 2.1.7.1[link], referred to above, and the performance of the two approaches described above, for this particular example, is shown in Fig. 2.1.8.1[link]. Briefly, the Fourier p.d.f. shows an excellent agreement with the histogram of recalculated [|E|] values, while the agreement attained by the Hermite correction factor is much less satisfactory, even for the (longest available to us) five-term expansion. It must be pointed out that (i) the inadequacy of `short' correction factors, in the example shown, is due to the large deviation from the ideal behaviour and (ii) the number of terms used there in the Fourier summation is twenty, whereafter the summation is terminated. Obviously, the computation of twenty (or more) Fourier coefficients is easier than that of five terms in the correction factor. The convergence of the Fourier series is very satisfactory. It appears that the (analytically) exact Fourier approach is the preferred one in cases of large or intermediate deviations, while the correction-factor approach may cope well with small ones. As far as the availability of symmetry-dependent centric and acentric p.d.f.'s is concerned, correction factors are available for all the space groups (see Table 2.1.7.1[link]), while Fourier coefficients of p.d.f.'s are available for the first 206 space groups (see Table 2.1.8.1[link]). It should be pointed out that p.d.f.'s based on the correction-factor method cope very well with cubic symmetries higher than [Fd\bar{3}], even if the asymmetric unit of the space group is strongly heterogeneous (Rabinovich et al., 1991b[link]).

[Figure 2.1.8.1]

Figure 2.1.8.1 | top | pdf |

The same recalculated histogram as in Fig. 2.1.7.1[link] along with the centric correction-factor p.d.f. [equation (2.1.7.5)[link]], truncated after two, three, four and five terms (dashed lines), and with that accurately computed for the correct space-group Fourier p.d.f. [equations (2.1.8.5)[link] and (2.1.8.22)[link]] (solid line).

Both approaches described in this section are related to the characteristic function of the required p.d.f. The correction-factor p.d.f.'s (2.1.7.5)[link] and (2.1.7.6)[link] can be obtained by expanding the logarithm of the appropriate characteristic function in a series of cumulants [e.g. equation (2.1.4.13)[link]; see also Shmueli & Wilson (1982[link])], truncating the series and performing its term-by-term Fourier inversion. The Fourier p.d.f., on the other hand, is computed by forming a Fourier series whose coefficients are exact analytical forms of the characteristic function at points related to the summation indices [e.g. equations (2.1.8.5)[link], (2.1.8.9)[link] and (2.1.8.11)[link], and Table 2.1.8.1[link]] and truncating the series when the terms become small enough.

References

Abramowitz, M. & Stegun, I. A. (1972). Handbook of Mathematical Functions. New York: Dover.
Barakat, R. (1974). First-order statistics of combined random sinusoidal waves with application to laser speckle patterns. Opt. Acta, 21, 903–921.
Bernstein, S. (1922). Sur la théorème limite du calcul des probabilités. Math. Ann. 85, 237–241.
Bernstein, S. (1927). Sur l'extension du théorème limite du calcul des probabilités aux sommes de quantités dépendantes. Math. Ann. 97, 1–59.
Cramér, H. (1951). Mathematical Methods of Statistics. Princeton University Press.
Faggiani, R., Lippert, B. & Lock, C. J. L. (1980). Heavy transition metal complexes of biologically important molecules. 4. Crystal and molecular structure of pentahydroxonium chloro(uracilato-N(1))(ethylene­diamine)platinum(II)chloride (H5O2)[PtCl­(NH2­CH2­CH2­NH2)­(C4­H5­N2O2)]Cl, and chloro(thyminato-N(1))­(ethyl­ene­di­amine)­platinum(II), PtCl(NH2CH2CH2NH2)­(C5H5­N2O2). Inorg. Chem. 19, 295–300.
French, S. & Wilson, K. (1978). On the treatment of negative intensity observations. Acta Cryst. A34, 517–525.
Gerhard, O. E. (1993). Line-of-sight velocity profiles in spherical galaxies: breaking the degeneracy between anisotropy and mass. Mon. Not. R. Astron. Soc. 265, 213–230.
Giacovazzo, C. (1977). On different probabilistic approaches to quartet theory. Acta Cryst. A33, 50–54.
Giacovazzo, C. (1980). Direct Methods in Crystallography. London: Academic Press.
Harker, D. (1953). The meaning of the average of [|F|^{2}] for large values of interplanar spacing. Acta Cryst. 6, 731–736.
Hauptman, H. & Karle, J. (1953). Solution of the Phase Problem. I. The Centrosymmetric Crystal. Am. Crystallogr. Assoc. Monograph No. 3. Dayton, Ohio: Polycrystal Book Service.
Howells, E. R., Phillips, D. C. & Rogers, D. (1950). The probability distribution of X-ray intensities. II. Experimental investigation and the X-ray detection of centers of symmetry. Acta Cryst. 3, 210–214.
International Tables for Crystallography (2005). Vol. A, Space-Group Symmetry, edited by Th. Hahn. Heidelberg: Springer.
International Tables for Crystallography (2004). Vol. C, Mathematical, Physical and Chemical Tables, edited by E. Prince. Dordrecht: Kluwer Academic Publishers.
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.
Lomer, T. R. & Wilson, A. J. C. (1975). Scaling of intensities. Acta Cryst. B31, 646–647.
Marel, R. P. van der & Franx, M. (1993). A new method for the identification of non-Gaussian line profiles in elliptical galaxies. Astrophys. J. 407, 525–539.
Myller-Lebedeff, W. (1907). Die Theorie der Integralgleichungen in Anwendung auf einige Reihenentwicklungen. Math. Ann. 64, 388–416.
Nigam, G. D. (1972). On the compensation of X-ray intensity. Indian J. Pure Appl. Phys. 10, 655–656.
Nigam, G. D. & Wilson, A. J. C. (1980). Compensation of excess intensity in space group P2. Acta Cryst. A36, 832–833.
Rabinovich, D. & Shakked, Z. (1984). A new approach to structure determination of large molecules by multi-dimensional search methods. Acta Cryst. A40, 195–200.
Rabinovich, S., Shmueli, U., Stein, Z., Shashua, R. & Weiss, G. H. (1991a). Exact random-walk models in crystallographic statistics. VI. P.d.f.'s of [|E|] for all plane groups and most space groups. Acta Cryst. A47, 328–335.
Rabinovich, S., Shmueli, U., Stein, Z., Shashua, R. & Weiss, G. H. (1991b). Exact random-walk models in crystallographic statistics. VII. An all-space-group study of the effects of atomic heterogeneity on the p.d.f.'s of [|E|]. Acta Cryst. A47, 336–340.
Rayleigh, Lord (1879). Investigations in optics with special reference to the spectroscope. Philos. Mag. 8, 261–274.
Rogers, D. (1950). The probability distribution of X-ray 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 X-ray intensities. V. A note on some hypersymmetric distributions. Acta Cryst. 6, 439–449.
Shmueli, U. (1979). Symmetry- and composition-dependent cumulative distributions of the normalized structure amplitude for use in intensity statistics. Acta Cryst. A35, 282–286.
Shmueli, U. (1982). A study of generalized intensity statistics: extension of the theory and practical examples. Acta Cryst. A38, 362–371.
Shmueli, U. & Kaldor, U. (1981). Calculation of even moments of the trigonometric structure factor. Methods and results. Acta Cryst. A37, 76–80.
Shmueli, U. & Kaldor, U. (1983). Moments of the trigonometric structure factor. Acta Cryst. A39, 615–621.
Shmueli, U., Rabinovich, S. & Weiss, G. H. (1989). Exact conditional distribution of a three-phase invariant in the space group P1. I. Derivation and simplification of the Fourier series. Acta Cryst. A45, 361–367.
Shmueli, U., Rabinovich, S. & Weiss, G. H. (1990). Exact random-walk models in crystallographic statistics. V. Non-symmetrically bounded distributions of structure-factor magnitudes. Acta Cryst. A46, 241–246.
Shmueli, U. & Weiss, G. H. (1985a). Centric, bicentric and partially bicentric intensity statistics. In Structure and Statistics in Crystallography, edited by A. J. C. Wilson, pp. 53–66. Guilderland: Adenine Press.
Shmueli, U. & Weiss, G. H. (1985b). Exact joint probability distributions for centrosymmetric structure factors. Derivation and application to the Σ1 relationship in the space group [P\bar{1}]. Acta Cryst. A41, 401–408.
Shmueli, U. & Weiss, G. H. (1986). Exact joint distribution of Eh , Ek and Eh+k, and the probability for the positive sign of the triple product in the space group [P\bar{1}]. Acta Cryst. A42, 240–246.
Shmueli, U. & Weiss, G. H. (1987). Exact random-walk models in crystallographic statistics. III. Distributions of [|E|] for space groups of low symmetry. Acta Cryst. A43, 93–98.
Shmueli, U. & Weiss, G. H. (1988). Exact random-walk models in crystallographic statistics. IV. P.d.f.'s of [|E|] allowing for atoms in special positions. Acta Cryst. A44, 413–417.
Shmueli, U., Weiss, G. H. & Kiefer, J. E. (1985). Exact random-walk models in crystallographic statistics. II. The bicentric distribution in the space group [P\bar{1}]. Acta Cryst. A41, 55–59.
Shmueli, U., Weiss, G. H., Kiefer, J. E. & Wilson, A. J. C. (1984). Exact random-walk models in crystallographic statistics. I. Space groups [P\bar{1}] and [P1]. Acta Cryst. A40, 651–660.
Shmueli, U., Weiss, G. H. & Wilson, A. J. C. (1989). Explicit Fourier representations of non-ideal hypercentric p.d.f.'s of [|E|]. Acta Cryst. A45, 213–217.
Shmueli, U. & Wilson, A. J. C. (1981). Effects of space-group symmetry and atomic heterogeneity on intensity statistics. Acta Cryst. A37, 342–353.
Shmueli, U. & Wilson, A. J. C. (1982). Intensity statistics: non-ideal distributions in theory and practice. In Crystallographic Statistics: Progress and Problems, edited by S. Ramaseshan, M. F. Richardson & A. J. C. Wilson, pp. 83–97. Bangalore: Indian Academy of Sciences.
Shmueli, U. & Wilson, A. J. C. (1983). Generalized intensity statistics: the subcentric distribution and effects of dispersion. Acta Cryst. A39, 225–233.
Spiegel, M. R. (1974). Theory and Problems of Fourier Analysis. Schaum's Outline Series. New York: McGraw-Hill.
Srinivasan, R. & Parthasarathy, S. (1976). Some Statistical Applications in X-ray Crystallography. Oxford: Pergamon Press.
Stuart, A. & Ord, K. (1994). Kendall's Advanced Theory of Statistics, Vol. 1, Distribution Theory, 6th ed. London: Edward Arnold.
Weiss, G. H. & Kiefer, J. E. (1983). The Pearson random walk with unequal step sizes. J. Phys. A, 16, 489–495.
Wilson, A. J. C. (1942). Determination of absolute from relative intensity data. Nature (London), 150, 151–152.
Wilson, A. J. C. (1949). The probability distribution of X-ray intensities. Acta Cryst. 2, 318–320.
Wilson, A. J. C. (1950). The probability distribution of X-ray 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 X-ray intensities. Research (London), 5, 588–589.
Wilson, A. J. C. (1956). The probability distribution of X-ray intensities. VII. Some sesquicentric distributions. Acta Cryst. 9, 143–144.
Wilson, A. J. C. (1964). The probability distribution of X-ray intensities. VIII. A note on compensation for excess average intensity. Acta Cryst. 17, 1591–1592.
Wilson, A. J. C. (1975). Effect of neglect of dispersion on apparent scale and temperature factors. In Anomalous Scattering, edited by S. Ramaseshan & S. C. Abrahams, pp. 325–332. Copenhagen: Munksgaard.
Wilson, A. J. C. (1976). Statistical bias in least-squares refinement. Acta Cryst. A32, 994–996.
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. (1978b). Variance of X-ray intensities: effect of dispersion and higher symmetries. Acta Cryst. A34, 986–994.
Wilson, A. J. C. (1978c). Statistical bias in scaling factors: Erratum. Acta Cryst. B34, 1749.
Wilson, A. J. C. (1979). Problems of resolution and bias in the experimental determination of the electron density and other densities in crystals. Acta Cryst. A35, 122–130.
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 X-ray reflections. Acta Cryst. A36, 945–946.
Wilson, A. J. C. (1981). Can intensity statistics accommodate stereochemistry? Acta Cryst. A37, 808–810.
Wilson, A. J. C. (1986a). Distributions of sums and ratios of sums of intensities. Acta Cryst. A42, 334–339.
Wilson, A. J. C. (1986b). Fourier versus Hermite representations of probability distributions. Acta Cryst. A42, 81–83.
Wilson, A. J. C. (1987a). Treatment of enhanced zones and rows in normalizing intensities. Acta Cryst. A43, 250–252.
Wilson, A. J. C. (1987b). Functional form of the ideal hypersymmetric distributions of structure factors. Acta Cryst. A43, 554–556.
Wilson, A. J. C. (1993). Space groups rare for organic structures. III. Symmorphism and inherent symmetry. Acta Cryst. A49, 795–806.








































to end of page
to top of page