Fast electrostatic solvers for kinetic Monte Carlo simulations

05/10/2019
by   William Robert Saunders, et al.
0

Kinetic Monte Carlo (KMC) is an important computational tool in theoretical physics and chemistry. In contrast to standard Monte Carlo, KMC permits the description of time dependent dynamical processes and is not restricted to systems in equilibrium. Compared to Molecular Dynamics, it allows simulations over significantly longer timescales. Recently KMC has been applied successfully in modelling of novel energy materials such as Lithium-ion batteries and organic/perovskite solar cells. Motivated by this, we consider general solid state systems which contain free, interacting particles which can hop between localised sites in the material. The KMC transition rates for those hops depend on the change in total potential energy of the system. For charged particles this requires the frequent calculation of electrostatic interactions, which is usually the bottleneck of the simulation. To avoid this issue and obtain results in reasonable times, many studies replace the long-range potential by a phenomenological short range approximation, which leads to systematic errors and unphysical results. On the other hand standard electrostatic solvers such as Ewald summation or fast Poisson solvers are highly inefficient in the KMC setup or introduce uncontrollable systematic errors at high resolution. In this paper we describe a new variant of the Fast Multipole Method by Greengard and Rokhlin which overcomes this issue by dramatically reducing computational costs. We construct an algorithm which scales linearly in the number of charges for each KMC step, something which had not been deemed to be possible before. We demonstrate the performance and parallel scalability of the method by implementing it in a performance portable software library. We describe the high-level Python interface of the code, which makes it easy to adapt to specific use cases.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
06/30/2020

Monte Carlo simulations of interacting particles with fast and accurate electrostatics

To minimise systematic errors in Monte Carlo simulations of charged part...
research
08/29/2020

Random Batch Algorithms for Quantum Monte Carlo simulations

Random batch algorithms are constructed for quantum Monte Carlo simulati...
research
12/21/2019

Simulating sticky particles: A Monte Carlo method to sample a Stratification

Many problems in materials science and biology involve particles interac...
research
06/17/2022

Accelerated Kinetic Monte Carlo methods for general nonlocal traffic flow models

This paper presents a class of one-dimensional cellular automata (CA) mo...
research
02/08/2021

Constrained Ensemble Langevin Monte Carlo

The classical Langevin Monte Carlo method looks for i.i.d. samples from ...
research
08/03/2017

Long range forces in a performance portable Molecular Dynamics framework

Molecular Dynamics (MD) codes predict the fundamental properties of matt...
research
02/09/2022

Optimising hadronic collider simulations using amplitude neural networks

Precision phenomenological studies of high-multiplicity scattering proce...

Please sign up or login with your details

Forgot password? Click here to reset