Researcher profile

Torsten Keßler

Torsten Keßler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
6topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

3 published item(s)

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.