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. 386388

3.4.4. Preliminary derivation to obtain a formula which accelerates the convergence of an sum over lattice points X(d)
The threedimensional directspace crystal lattice is specified by the origin vectors , and . A general vector in direct space is defined as where are the fractional cell coordinates of X. A lattice vector in direct space is defined as where are integers (specifying particular values of ) designating a lattice point. is the directcell volume which is equal to . A general point in the direct lattice is X(x); the contents of the lattice are by definition identical as the components of x are increased or decreased by integer amounts.
The reciprocallattice vectors are defined by the relations A general vector in reciprocal space H(r) is defined as A reciprocallattice vector H(h) is defined by the integer triplet (specifying particular values of ) so that In other sections of this volume a shortened notation h is used for the reciprocallattice vector. In this section the symbol H(h) is used to indicate that it is a particular value of H(r).
The threedimensional Fourier transform of a function is defined by The Fourier transform of the set of points defining the direct lattice is the set of points defining the reciprocal lattice, scaled by the directcell volume. It is useful for our purpose to express the lattice transform in terms of the Dirac delta function which is defined so that for any function We then write First consider the lattice sum over the directlattice points X(d), relative to a particular point , with omission of the origin lattice point. The special case with will also be needed: Now define a sum of Dirac delta functions Then S′ can be represented as an integral in which a term is contributed to S′ whenever the directspace vector X coincides with the lattice vector X(d), except for . Now apply the convergence function to S′:
The first integral is shown here only for the purpose of giving a consistent representation of S′; in fact, the first integral will be reconverted back into a sum and evaluated in direct space. The second integral will be transformed to reciprocal space using Parseval's theorem [see, for example, Arfken (1970)], which states that Considering only the second integral in the formula for S′ and explicitly introducing the term we have where the unprimed f includes the term which was earlier omitted from f′: Using Parseval's theorem, and evaluating the origin term, we have The Fourier transform of the complement of the incomplete gamma function divided by is (Nijboer & DeWette, 1957) If there is a change of origin and the point is used instead of X the transform is Evaluation of the two Fourier transforms in the first term gives Because of the presence of the Dirac delta function in each integral, we can convert the integrals with h unequal to zero into a sum The term needs to be evaluated in the limit. Clearly, the complex exponential goes to unity. If n is greater than 3 the limit of the indeterminate form infinity/infinity is needed: The limit can be found by L'Hospital's rule [see, for example, Widder (1961)] which states that if and both approach infinity as x approaches a constant, c, and the limit of the ratio of the first derivatives and exists, that limit is also true for the limit of the ratio of the functions: To differentiate the definite integral function, Leibnitz's formula may be used [see, for example, Arfken (1970)]. This formula states that In our case, x becomes ; f becomes which is independent of ; g becomes ; and h is infinite. Thus only the last term of Leibnitz's formula is nonzero and we obtain for the ratio of the first derivatives so that the limiting value for the term for n greater than 3 is The final result for S′ is
The significance of the terms is as follows. The first term represents the convergenceaccelerated direct sum, which does not include the origin term; the next term, also in direct space, corrects for the remainder resulting from the subtraction of the origin term; the third term comes from Parseval's theorem and is a sum over the nonzero h reciprocallattice points; and the last term is the reciprocallattice term.
If the second term becomes an indeterminate form 0/0. The limit can be found with use of L'Hospital's rule again, this time for the 0/0 form. We need the limit of , where and . To differentiate the incomplete gamma function, we can again use Leibnitz's formula. In this case only the second term of the formula is nonzero and we obtain for the ratio of the first derivatives so that the limiting value for this term as approaches zero is Therefore, the value of the sum when is
References
Arfken, G. (1970). Mathematical methods for physicists, 2nd ed. New York: Academic Press.Nijboer, B. R. A. & DeWette, F. W. (1957). On the calculation of lattice sums. Physica (Utrecht), 23, 309–321.
Widder, D. V. (1961). Advanced calculus, 2nd ed. New York: PrenticeHall.