International
Tables for Crystallography Volume C Mathematical, physical and chemical tables Edited by E. Prince © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. C. ch. 6.1, p. 588
|
In a cumulant expansion (Johnson & Levy, 1974), the entire series is expressed in exponential form. The cumulant expansion about S = 0 for the generalized temperature factor is where the coefficient tensor , a symmetric tensor of order p, is the pth-order cumulant. The inverse Fourier transform is the Edgeworth expansion around the Gaussian p.d.f. Cumulants can be expressed in terms of moments and vice versa. The pth moment (if it exists) of a general p.d.f., ρ(x), is a symmetric tensor defined as The relations between the lower-order moments and cumulants are and, conversely, In the Gram–Charlier and Fourier-invariant expansions, the Fourier-transform relationship between the p.d.f. and the temperature factor to given order can be made exact. Each cumulant contributes to all higher-order quasi-moment terms and vice versa. Hence, a given cumulant expansion is to an extent arbitrarily truncated (Kuhs, 1983). Care is required when interpreting the coefficients (Zucker & Schulz, 1982).
On the other hand, the cumulant expansion has the advantage of yielding tractable expressions for the one-particle potential in the quantum regime (Mair, 1980a). In that regime, equation (6.1.1.36) for the one-particle potential is invalid, and the expressions relating V(u) to ρ(u) in the Gram–Charlier and Fourier-invariant expansions are cumbersome (Mair & Wilkins, 1976).
Coefficients obtained by applying least-squares methods to structure-factor equations related to the truncated cumulant expansions do not necessarily yield non-negative p.d.f.'s nor are the linear-term coefficients necessarily faithful representations of the mean. Caution must be exercised in interpreting the results.
All the methods are satisfactory in the case of rapidly converging potential series. The methods are equivalent up to λ2 in the van Hove order parameter (Mair, 1980b). Difficulties are encountered with convergence of the series in the case of strong anharmonicity, in which case numerical or alternative analytical models may be necessary. If the anharmonicity is such that the difference between the expansions is significant, it may be preferable to evaluate the Fourier transforms directly, as recommended by Mackenzie & Mair (1985).
References
Johnson, C. K. & Levy, H. A. (1974). Thermal motion of independent atoms. International tables for X-ray crystallography. Vol. IV, pp. 317–319. Birmingham: Kynoch Press.Google ScholarKuhs, W. F. (1983). Statistical description of multimodal atomic probability densities. Acta Cryst. A39, 149–158.Google Scholar
Mackenzie, J. K. & Mair, S. L. (1985). Anharmonic temperature factors: the limitations of perturbation-theory expressions. Acta Cryst. A41, 81–85.Google Scholar
Mair, S. L. (1980a). Temperature dependence of the anharmonic Debye–Waller factor. J. Phys. C, 13, 2857–2868.Google Scholar
Mair, S. L. (1980b). The anharmonic Debye–Waller factor in the classical limit. J. Phys. C, 13, 1419–1425.Google Scholar
Mair, S. L. & Wilkins, S. W. (1976). Anharmonic Debye–Waller factor using quantum statistics. J. Phys. C, 9, 1145–1158.Google Scholar
Zucker, U. H. & Schulz, H. (1982). Statistical approaches for the treatment of anharmonic motion in crystals. I. A comparison of the most frequently used formalisms of anharmonic thermal vibrations. Acta Cryst. A38, 563–568.Google Scholar