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

International Tables for Crystallography (2010). Vol. B, ch. 1.3, pp. 97-98   | 1 | 2 |

Section 1.3.4.4.7.5. Cruickshank's modified Fourier method

G. Bricognea

aGlobal Phasing Ltd, Sheraton House, Suites 14–16, Castle Park, Cambridge CB3 0AX, England, and LURE, Bâtiment 209D, Université Paris-Sud, 91405 Orsay, France

1.3.4.4.7.5. Cruickshank's modified Fourier method

| top | pdf |

Cruickshank consolidated and extended Cochran's derivations in a series of classic papers (Cruickshank, 1949b[link], 1950[link], 1952[link], 1956[link]). He was able to show that all the coefficients involved in the right-hand side and normal matrix of the least-squares method could be calculated by means of suitable differential Fourier syntheses even when the atoms overlap. This remarkable achievement lay essentially dormant until its independent rediscovery by Agarwal in 1978 (Section 1.3.4.4.7.6[link]).

To ensure rigorous equivalence between the summations over [{\bf h} \in {\scr H}] (in the expressions of least-squares right-hand side and normal matrix elements) and genuine Fourier summations, multiplicity-corrected weights were introduced by:[\eqalign{ w_{{\bf h}} &= {1 \over |G_{{\bf h}}|} W_{{\bf h}}\quad \hbox{ if } {\bf h} \in G{\bf h}\quad \hbox{ with } {\bf h} \in {\scr H},\cr w_{\bf h} &= 0\phantom{{241 \over |G_{{\bf h}}{12 \over 1}|}}\quad \hbox{otherwise},}]where Gh denotes the orbit of h and [G_{{\bf h}}] its isotropy subgroup (Section 1.3.4.2.2.5[link]). Similarly, derivatives with respect to parameters of symmetry-unique atoms were expressed, via the chain rule, as sums over the orbits of these atoms.

Let [p = 1, \ldots, n] be the label of a parameter [u_{p}] belonging to atoms with label j. Then Cruickshank showed that the pth element of the right-hand side of the normal equations can be obtained as [D_{p, \, j} ({\bf x}_{j})], where [D_{p, \, j}] is a differential synthesis of the form[\eqalign{ D_{p, \, j} ({\bf x}) &= {\textstyle\sum\limits_{{\bf h}}} P_{p} ({\bf h}) g_{j} ({\bf h}) w_{{\bf h}} (|F_{{\bf h}}^{\,\rm calc}| - |F_{{\bf h}}|^{\rm obs})\cr &\quad \times \exp (i\varphi_{{\bf h}}^{\rm calc}) \exp (-2 \pi i{\bf h} \cdot {\bf x})}]with [P_{p}({\bf h})] a polynomial in (h, k, l) depending on the type of parameter p. The correspondence between parameter type and the associated polynomial extends Booth's original range of differential syntheses, and is recapitulated in the following table.[\eqalign{ &\quad\hbox{Parameter type}\qquad P (h, k, l)\cr &\overline{\hskip13pc}\cr &\quad x \hbox{ coordinate} \quad \qquad - 2 \pi ih\cr &\quad y \hbox{ coordinate} \quad \qquad - 2 \pi ik\cr &\quad z \hbox{ coordinate} \quad \qquad - 2 \pi il\cr &\quad B \hbox{ isotropic} \quad \qquad \,- {\textstyle{1 \over 4}} (d_{{\bf h}}^{*})^{2}\cr &\quad B^{11} \hbox{ anisotropic} \qquad - h^{2}\cr &\quad B^{12} \hbox{ anisotropic} \qquad - hk\cr &\quad B^{13} \hbox{ anisotropic} \qquad - hl\cr &\quad B^{22} \hbox{ anisotropic} \qquad - k^{2}\cr &\quad B^{23} \hbox{ anisotropic} \qquad - kl\cr &\quad B^{33} \hbox{ anisotropic} \qquad - l^{2}.}]Unlike Cochran's original heuristic argument, this result does not depend on the atoms being resolved.

Cruickshank (1952)[link] also considered the elements of the normal matrix, of the form[\sum\limits_{{\bf h}} w_{{\bf h}} {\partial |F_{{\bf h}}^{\,\rm calc}| \over \partial u_{p}} {\partial |F_{{\bf h}}^{\,\rm calc}| \over \partial u_{q}}]associated with positional parameters. The [3 \times 3] block for parameters [{\bf x}_{j}] and [{\bf x}_{k}] may be written[\eqalign{&{\textstyle\sum\limits_{{\bf h}}} w_{{\bf h}} ({\bf hh}^{T}) {\scr Re} [(-2 \pi i) g_{j} ({\bf h}) \exp (-2 \pi i{\bf h} \cdot {\bf x}_{j}) \exp (i \varphi_{{\bf h}}^{\rm calc})]\cr &\quad \times {\scr Re} [(-2 \pi i) g_{k} ({\bf h}) \exp (-2 \pi i{\bf h} \cdot {\bf x}_{k}) \exp (i \varphi_{{\bf h}}^{\rm calc})]}]which, using the identity[{\scr Re}(z_{1}) {\scr Re}(z_{2}) = {\textstyle{1 \over 2}}[{\scr Re} (z_{1} z_{2}) + {\scr Re} (z_{1} \overline{z_{2}})],]becomes[\eqalign{&2\pi^{2} {\textstyle\sum\limits_{{\bf h}}} w_{{\bf h}} ({\bf hh}^{T}) g_{j}({\bf h})g_{k}({\bf h})\cr &\quad \times \{\exp [-2\pi i{\bf h} \cdot ({\bf x}_{j} - {\bf x}_{k})]\cr &\quad - \exp (2i\varphi_{{\bf h}}^{\rm calc}) \exp [ - 2\pi i{\bf h} \cdot ({\bf x}_{j} + {\bf x}_{k})]\}}](Friedel's symmetry makes [{\scr Re}] redundant on the last line). Cruickshank argued that the first term would give a good approximation to the diagonal blocks of the normal matrix and to those off-diagonal blocks for which [{\bf x}_{j}] and [{\bf x}_{k}] are close. On this basis he was able to justify the `n-shift rule' of Shoemaker et al. (1950)[link]. Cruickshank gave this derivation in a general space group, but using a very terse notation which somewhat obscures it. Using the symmetrized trigonometric structure-factor kernel [\Xi^{-}] of Section 1.3.4.2.2.9[link] and its multiplication formula, the above expression is seen to involve the values of a Fourier synthesis at points of the form [{\bf x}_{j} \pm S_{g} ({\bf x}_{k})].

Cruickshank (1956)[link] showed that this analysis could also be applied to the refinement of temperature factors.

These two results made it possible to obtain all coefficients involved in the normal equations by looking up the values of certain differential Fourier syntheses at [{\bf x}_{j}] or at [{\bf x}_{j} \pm S_{g} ({\bf x}_{k})]. At the time this did not confer any superiority over the standard form of the least-squares procedure, because the accurate computation of Fourier syntheses was an expensive operation. The modified Fourier method was used by Truter (1954)[link] and by Ahmed & Cruickshank (1953a)[link], and was incorporated into the program system described by Cruickshank et al. (1961)[link]. A more recent comparison with the least-squares method was made by Dietrich (1972)[link].

There persisted, however, some confusion about the nature of the relationship between Fourier and least-squares methods, caused by the extra factors [g_{j} ({\bf h})] which make it necessary to compute a differential synthesis for each type of atom. This led Cruickshank to conclude that `in spite of their remarkable similarities the least-squares and modified-Fourier methods are fundamentally distinct'.

References

Ahmed, F. R. & Cruickshank, D. W. J. (1953a). A refinement of the crystal structure analysis of oxalic acid dihydrate. Acta Cryst. 6, 385–392.
Cruickshank, D. W. J. (1949b). The accuracy of atomic co-ordinates derived by least-squares or Fourier methods. Acta Cryst. 2, 154–157.
Cruickshank, D. W. J. (1950). The convergence of the least-squares and Fourier refinement methods. Acta Cryst. 3, 10–13.
Cruickshank, D. W. J. (1952). On the relations between Fourier and least-squares methods of structure determination. Acta Cryst. 5, 511–518.
Cruickshank, D. W. J. (1956). The determination of the anisotropic thermal motion of atoms in crystals. Acta Cryst. 9, 747–753.
Cruickshank, D. W. J., Pilling, D. E., Bujosa, A., Lovell, F. M. & Truter, M. R. (1961). Crystallographic calculations on the Ferranti Pegasus and mark I computers. In Computing Methods and the Phase Problem in X-Ray Crystal Analysis, edited by R. Pepinsky, J. M. Robertson & J. C. Speakman, pp. 32–78. Oxford: Pergamon Press.
Dietrich, H. (1972). A reconsideration of Fourier methods for the refinement of crystal structures. Acta Cryst. B28, 2807–2814.
Shoemaker, D. P., Donohue, J., Schomaker, V. & Corey, R. B. (1950). The crystal structure of LS-threonine. J. Am. Chem. Soc. 72, 2328–2349.
Truter, M. R. (1954). Refinement of a non-centrosymmetrical structure: sodium nitrite. Acta Cryst. 7, 73–77.








































to end of page
to top of page