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

International Tables for Crystallography (2006). Vol. B, ch. 3.4, pp. 385-386   | 1 | 2 |

Section 3.4.3. Preliminary description of the method

D. E. Williamsa

aDepartment of Chemistry, University of Louisville, Louisville, Kentucky 40292, USA

3.4.3. Preliminary description of the method

| top | pdf |

Ewald (1921[link]) developed a method which modified the mathematical representation of the Coulombic lattice sum to improve the rate of convergence. This method was based on partially transforming the lattice sum into reciprocal space. Bertaut (1952[link]) presented another method for derivation of the Ewald result which used the concept of the crystallographic structure factor. His formula extended the Ewald treatment to a composite lattice with more than one atom per lattice point. Nijboer & DeWette (1957[link]) developed a general Fourier transform method for the evaluation of [R^{-n}] sums in simple lattices. Williams (1971[link]) extended this treatment to a composite lattice and gave general formulae for the [R^{-n}] sums for any crystal. A review article, on which this chapter is based, appeared later (Williams, 1989[link]).

Consider a function, W(R), which is unity at [R = 0] and smoothly declines to zero as R approaches infinity. If each term of the lattice sum is multiplied by W(R), the rate of convergence is increased. However, the rate of convergence of the remainder of the original sum, which contains the difference terms, is not increased.[\eqalign{ V_{n} &= (1/2) {\textstyle\sum\limits_{j}^{\rm one\; cell}} \ {\textstyle\sum\limits_{k}^{\rm all\; cells}\!\!\!'}\ Q_{jk} R_{jk}^{-n} W(R)\cr &\quad + (1/2) {\textstyle\sum\limits_{j}^{\rm one\; cell}} \ {\textstyle\sum\limits_{k}^{\rm all\; cells}\!\!\!'}\ Q_{jk} R_{jk}^{-n} [1 - W(R)].}]

In the accelerated-convergence method the difference terms are expressed as an integral of the product of two functions. According to Parseval's theorem (described below) this integral is equal to an integral of the product of the two Fourier transforms of the functions. Finally, the integral over the Fourier transforms of the functions is converted to a sum in reciprocal (or Fourier-transform) space. The choice of the convergence function W(R) is not unique; an obvious requirement is that the relevant Fourier transforms must exist and have correct limiting behaviour. Nijboer and DeWette suggested using the incomplete gamma function for W(R). More recently, Fortuin (1977[link]) showed that this choice of convergence function leads to optimal convergence of the sums in both direct and reciprocal space: [W(R) = \Gamma (n/2, \pi w^{2} R^{2})/\Gamma (n/2),] where [\Gamma (n/2)] and [\Gamma (n/2, \pi w^{2} R^{2})] are the gamma function and the incomplete gamma function, respectively: [\Gamma (n/2, \pi w^{2} R^{2}) = \textstyle\int\limits_{\pi w^{2} R^{2}}^{\infty} t^{(n/2)-1} \exp (-t)\;\hbox{d}t] and [\Gamma (n/2) = \Gamma (n/2, 0).]

The complement of the incomplete gamma function is [\gamma (n/2, \pi w^{2} R^{2}) = \Gamma (n/2) - \Gamma (n/2, \pi w^{2} R^{2}).]

References

Bertaut, E. F. (1952). L'énergie électrostatique de réseaux ioniques. J. Phys. (Paris), 13, 499–505.
Ewald, P. P. (1921). Die Berechnung optischer und elektrostatischer Gitterpotentiale. Ann. Phys. (Leipzig), 64, 253–287.
Fortuin, C. M. (1977). Note on the calculation of electrostatic lattice potentials. Physica (Utrecht), 86A, 574–586.
Nijboer, B. R. A. & DeWette, F. W. (1957). On the calculation of lattice sums. Physica (Utrecht), 23, 309–321.
Williams, D. E. (1971). Accelerated convergence of crystal lattice potential sums. Acta Cryst. A27, 452–455.
Williams, D. E. (1989). Accelerated convergence treatment of [{\bi R}^{-n}] lattice sums. Crystallogr. Rev. 2, 3–23. Corrections: Crystallogr. Rev. 2, 163–166.








































to end of page
to top of page