On the Efficient Computation of Large Scale Singular Sums with Applications to Long-Range Forces in Crystal Lattices

Abstract

We develop a new expansion for representing singular sums in terms of integrals and vice versa. This method provides a powerful tool for the efficient computation of large singular sums that appear in long-range interacting systems in condensed matter and quantum physics. It also offers a generalised trapezoidal rule for the precise computation of singular integrals. In both cases, the difference between sum and integral is approximated by derivatives of the non-singular factor of the summand function, where the coefficients in turn depend on the singularity. We show that for a physically meaningful set of functions, the error decays exponentially with the expansion order. For a fixed expansion order, the error decays algebraically both with the grid size, if the method is used for quadrature, or the characteristic length scale of the summand function in case the sum over a fixed grid is approximated by an integral. In absence of a singularity, the method reduces to the Euler–Maclaurin summation formula. We demonstrate the numerical performance of our new expansion by applying it to the computation of the full nonlinear long-range forces inside a domain wall in a macroscopic one-dimensional crystal with 2times 10^{10} particles. The code of our implementation in Mathematica is provided online. For particles that interact via the Coulomb repulsion, we demonstrate that finite size effects remain relevant even in the thermodynamic limit of macroscopic particle numbers. Our results show that widely-used continuum limits in condensed matter physics are not applicable for quantitative predictions in this case.