Source author record

Torsten Keßler

Torsten Keßler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis cond-mat.other math-ph math.MP math.NT

Catalog footprint

What is connected

3works

6topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Singular Euler-Maclaurin expansion

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\times 10^{10}$ particles. A reference implementation in Mathematica is provided online.

preprint2021arXiv

Singular Euler-Maclaurin expansion on multidimensional lattices

We extend the classical Euler-Maclaurin expansion to sums over multidimensional lattices that involve functions with algebraic singularities. This offers a tool for the precise quantification of the effect of microscopic discreteness on macroscopic properties of a system. First, the Euler-Maclaurin summation formula is generalised to lattices in higher dimensions, assuming a sufficiently regular summand function. We then develop this new expansion further and construct the singular Euler-Maclaurin (SEM) expansion in higher dimensions, an extension of our previous work in one dimension, which remains applicable and useful even if the summand function includes a singular function factor. We connect our method to analytical number theory and show that all operator coefficients can be efficiently computed from derivatives of the Epstein zeta function. Finally we demonstrate the numerical performance of the expansion and efficiently compute singular lattice sums in infinite two-dimensional lattices, which are of high relevance in solid state and quantum physics. An implementation in Mathematica is provided online along with this article.

preprint2020arXiv

Vlasov-Poisson system tackled by particle simulation utilising boundary element methods

This paper presents a grid-free simulation algorithm for the fully three-dimensional Vlasov--Poisson system for collisionless electron plasmas. We employ a standard particle method for the numerical approximation of the distribution function. Whereas the advection of the particles is grid-free by its very nature, the computation of the acceleration involves the solution of the non-local Poisson equation. To circumvent a volume mesh, we utilise the Fast Boundary Element Method, which reduces the three-dimensional Poisson equation to a system of linear equations on its two-dimensional boundary. This gives rise to fully populated matrices which are approximated by the $\mathcal H^2$-technique, reducing the computational time from quadratic to linear complexity. The approximation scheme based on interpolation has shown to be robust and flexible, allowing a straightforward generalisation to vector-valued functions. In particular, the Coulomb forces acting on the particles are computed in linear complexity. In first numerical tests, we validate our approach with the help of classical non-linear plasma phenomena. Furthermore, we show that our method is able to simulate electron plasmas in complex three-dimensional domains with mixed boundary conditions in linear complexity.