English

Singular Euler-Maclaurin expansion

Numerical Analysis 2022-01-28 v3 Other Condensed Matter Numerical Analysis

Abstract

We present the singular Euler--Maclaurin expansion, a new method for the efficient computation of large singular sums that appear in long-range interacting systems in condensed matter and quantum physics. In contrast to the traditional Euler--Maclaurin summation formula, the new method is applicable also to the product of a differentiable function and a singularity. For suitable non-singular functions, we show that the approximation error decays exponentially in the expansion order and polynomially in the characteristic length scale of the non-singular function, where precise error estimates are provided. The sum is approximated by an integral plus a differential operator acting on the non-singular function factor only. The singularity furthermore is included in a generalisation of the Bernoulli polynomials that form the coefficients of the differential operator. We demonstrate the numerical performance of the singular Euler--Maclaurin expansion by applying it to the computation of the full non-linear long-range forces inside a macroscopic one-dimensional crystal with 2×10102\times 10^{10} particles. A reference implementation in Mathematica is provided online.

Keywords

Cite

@article{arxiv.2003.12422,
  title  = {Singular Euler-Maclaurin expansion},
  author = {Andreas A. Buchheit and Torsten Keßler},
  journal= {arXiv preprint arXiv:2003.12422},
  year   = {2022}
}

Comments

Accompanying code: https://github.com/andreasbuchheit/singular_euler_maclaurin

R2 v1 2026-06-23T14:29:20.570Z