

The program uses fast convergent series for lattice sums L_{s} with real exponents s>3 for LennardJones type of potentials (~r^{ s}) associated with the simple cubic, facecentred cubic, bodycentred cubic and hexagonal closedpacked lattices. Four different methods are used: 1) Direct summation (good for large exponents); 2) Dirichlet series (number theoretical treatment, good for large exponents); 3) Terras' decomposition of the Epstein zeta function for pure quadratic forms (good for small exponents); 4) Van der Hoff  Benson expansion (good for small exponents). This is an opensource code and free for distribution (except for commercial purposes as it contains a modified version for the Bessel function published in Numerical Recipes). For citations please use:  Antony Burrows and Peter Schwerdtfeger, Program Jones  A Fortran program for calculating lattices sums for the cubic and hexagonal lattices, Massey University, 2019. (pdf file with more details available upon request). For an introduction to extended LennardJones potentials and the use of corresponding lattice sums for solidstate properties see:  P. Schwerdtfeger, N. Gaston, R. P. Krawczyk, R. Tonner, and G. E. Moyano, Extension of the LennardJones potential: Theoretical investigations into raregas clusters and crystal lattices of He, Ne, Ar and Kr using manybody interaction expansions, Phys. Rev. B 73, 064112 (2006). A more detailed paper is on the way. For a general introduction into lattice sums see:  J. M. Borwein, M. Glasser, R. McPhedran, J. Wan, and I. Zucker, Lattice sums then and now, 150 (Cambridge University Press, 2013). To download the fortran program click here: For the number theoretical treatment in terms of the Dirichlet series a very small database containing the number of representations for an integer r_{3}(n) with n≤100 is included in the program. Sloane's bigger data base for the four lattices with n≤10,000 can also be used and can be downloaded here (recreated by using Program Jones): Download for the Sloane Database A very large database (533.4 MB) for n≤20,000,000 well exceeding Sloane's data can be used as well and can downloaded here: Download for Extended Sloane Database (compressed 178.1 MB) History: A detailed historical survey on lattice sums can be found in Borwein's book. For the sc, bcc, and fcc lattices, the first attempt to calculate the corresponding lattice sums to sufficient accuracy was by LennardJones and Ingham, who used expansions in terms of Bessel functions. Details can also be found in the work by M. Born R. D. Misra [On the stability of crystal lattices. iv, Math. Proc. Cambridge Phil. Soc. 36, 466478 (1940)]. In 1953 B. M. E. Van der Hoff and G. C. Benson introduced an expansion in terms of Bessel functions for lattice sums [A method for the evaluation of some lattice sums occurring in calculations of physical properties of crystals, Can. J. Phys. 31, 10871094 (1953)]. In 1973, A. A. Terras introduced an expansion of quadratic forms in terms of Bessel functions to successively reduce the dimension N in the Epstein zeta function by blockdiagonalizing the matrix S of the quadratic form [Bessel series expansions of the Epstein zeta function and the functional equation, Trans. Am. Math. Soc. 183, 477486 (1973)] shown at the top of the webpage. Further improvement came from the work by I. J. Zucker in 1975 who evaluated lattice sums and Madelung constants for cubic and tetragonal lattices to 10 significant digits by using a Mellin transformation in conjunction with Jacobian θfunctions (Madelung constants and lattice sums for invariant cubic lattice complexes and certain tetragonal structures, J. Phys. A: Math. Gen. 8, 1734 (1975); and many other papers]. The lattice sum for the hexagonal closed packed structure was given early by G. Kane and M. GoeppertMayer in 1940 [Lattice summations for hexagonal closepacked crystals, J. Chem. Phys. 8, 642 (1940)], but has never been evaluated to high precision as fast convergent expansions for this lattice were unknown so far. Note: We are currently working on lattice sums in N dimensions using Terras' decomposition of the Epstein zeta function and analytical continuations. 
Maintained by Peter Schwerdtfeger  Last updated: June 2019  Copyright 2014  Massey University
