International
Tables for Crystallography Volume B Reciprocal space Edited by U. Shmueli © International Union of Crystallography 2006 
International Tables for Crystallography (2006). Vol. B, ch. 3.4, pp. 391392

Consider the case of the sodium chloride crystal structure (a facecentred cubic structure) as a simple example for evaluation of the Coulombic sum. The sodium ion can be taken at the origin, and the chloride ion halfway along an edge of the unit cell. The results can easily be generalized for this structure type by using the unitcell edge length, a, as a scaling constant.
First, consider the nearest neighbours. Each sodium and each chloride ion is surrounded by six ions of opposite sign at a distance of . The Coulombic energy for the first coordination sphere is . Table 3.4.2.1 shows that the converged value of the lattice energy is . Thus the nearestneighbour energy is over three times more negative than the total lattice energy. In the second coordination sphere each ion is surrounded by 12 similar ions at a distance of . The energy contribution of the second sphere is . Thus, major cancellation occurs and the net energy for the first two coordination spheres is which actually has the wrong sign for a stable crystal. The third coordination sphere again makes a negative contribution. Each ion is surrounded by eight ions of opposite sign at a distance of . The energy contribution is , now giving a total so far of . In the fourth coordination sphere each ion is surrounded by six others of the same sign at a distance of a. The energy contribution is to yield a total of .
It is seen immediately by examining the numbers that the Coulombic sum is converging very slowly in direct space. Madelung (1918) devised a method for accurate evaluation of the sodium chloride lattice sum. However, his method is not generally applicable for more complex lattice structures. Evjen (1932) emphasized the importance of summing over a neutral domain, and replaced the sum with an integral outside of the first few shells of nearest neighbours. But the method of Ewald remained as the only completely general and accurate method of evaluating the Coulombic sum for a general lattice. Although it was derived in a somewhat different way, Ewald's method is equivalent to accelerated convergence for the special case of .
In the reciprocal lattice of sodium chloride only points with indices (hkl) all even or all odd are permitted by the facecentred symmetry. The reciprocal cell has edge length and the reciprocalaxis directions coincide with the directlattice axis directions. The closest points to the origin are the eight (111) forms at a distance of . For sodium chloride this distance is 0.3078 Å^{−1}.
Table 3.4.12.1 shows the effect of convergence acceleration on the directspace portion of the sum for the sodium chloride structure. The constant term is included in the values given. This constant term is always large if w is not zero; for instance, when this term is (Table 3.4.12.2). For the reciprocallattice sum is zero to six figures. Thus, only the direct sum (plus the constant term) is needed, evaluated out to 20 Å in direct space, to obtain sixfigure accuracy. As shown in Table 3.4.2.1 above, the same summation effort without the use of accelerated convergence gave 8% error, or only slightly better than onefigure accuracy. The acceleratedconvergence technique therefore yielded nearly five orders of magnitude improvement in accuracy, even without evaluation of the reciprocallattice sum.


The column showing shows an example of how the reciprocallattice sum can also be neglected if lower accuracy is required. Table 3.4.12.2 shows that the reciprocallattice sum is still only 0.003. But now the directlattice sum only needs to be evaluated out to 14 Å, with further savings in calculation effort. For w values larger than 0.15 the reciprocal sum is needed. For this sum must be evaluated out to 0.8 Å^{−1} to obtain sixfigure accuracy.
Table 3.4.12.3 illustrates an application for the dispersion sum. When five figures of accuracy can be obtained without consideration of the reciprocal sum. The direct sum is required out to 18 Å. If , better than fourfigure accuracy can still be obtained without evaluating the reciprocallattice sum. In this case, the direct lattice needs to be summed only to 12 Å, and there is a saving of an order of magnitude in the length of the calculation. As with the Coulombic sum, if w is greater than 0.15 the reciprocallattice summation is needed; Table 3.4.12.4 shows the values.


The time required to obtain a lattice sum of given accuracy will vary depending on the particular structure considered and of course on the computer and program which are used. An example of timing for the benzene dispersion sum is given in Table 3.4.12.5 for the PCK83 program (Williams, 1984) running on a VAX11/750 computer. In this particular case direct terms were evaluated at a rate of about 200 terms s^{−1} and reciprocal terms, being a sum themselves, were evaluated at a slower rate of about 5 terms s^{−1}.

Table 3.4.12.5 shows the time required for evaluation of the dispersion sum using various values of the convergence constant, w. The timing figures show that there is an optimum choice for w; for the PCK83 program the optimum value indicated is 0.15–0.2 Å^{−1}. In the program of Pietila & Rasmussen (1984) values in the range 0.15–0.2 Å^{−1} are also suggested. For the WMIN program (Busing, 1981) a slightly higher value of 0.25 Å^{−1} is suggested. Trial calculations can be used to determine the optimum value of w for the situation of a particular crystal structure, program and computer.
References
Busing, W. R. (1981). WMIN, a computer program to model molecules and crystals in terms of potential energy functions. Oak Ridge National Laboratory Report ORNL5747. Oak Ridge, Tennessee 37830, USA.Evjen, H. M. (1932). The stability of certain heteropolar crystals. Phys. Rev. 39, 675–694.
Madelung, E. (1918). Das elektrische Feld in Systemen von regelmässig angeordneten Punktladungen. Phys. Z. 19, 524–532.
Pietila, L.O. & Rasmussen, K. (1984). A program for calculation of crystal conformations of flexible molecules using convergence acceleration. J. Comput. Chem. 5, 252–260.
Williams, D. E. (1984). PCK83, a crystal molecular packing analysis program. Quantum Chemistry Program Exchange, Department of Chemistry, Indiana University, Bloomington, Indiana 47405, USA.