Efficient kinetic Monte Carlo simulations with long-range electrostatic interactions
Abstract
Charge transport in solid-state materials is often governed by charge carrier hopping processes in the presence of long-range electrostatic interactions. Kinetic Monte Carlo (kMC) simulations provide a framework for describing such rare-event dynamics over extended timescales. However, the efficient treatment of long-range interactions remains a major computational challenge, since each particle jump modifies the energy landscape globally and, in principle, requires updating all transition rates after every carrier motion. We present an efficient and generally applicable update process for these transition rates in the presence of long-range electrostatic interactions. To illustrate the method, we study charge diffusion on a simple cubic host lattice under an external electric field. The transport behavior is investigated in high states of charge (SOC), and the influence of temperature, electric field strength, and SOC is examined. Our simulations show strongly suppressed transport at 100\,\% SOC (50\,\% occupation) caused by a freezing of the charge carriers into a Coulomb superlattice. Slight deviations from this reference in terms of charge carrier concentration lead to a rapid increase in conductivity. A deeper analysis reveals that this behavior can be rationalized as transport of non-interacting defects in the Coulomb superlattice. These results demonstrate the capability of the proposed update procedure to efficiently capture non-equilibrium transport phenomena in interacting charged systems and provide a foundation for simulations of more complex charge diffusion problems.
Related
- Atomistic Insights into MOCVD Growth of GaAs using a Closed Loop Autonomous System Combining Directed Reinforcement Learning and Kinetic Monte Carlo Simulations
- SmoQyElPhQMC.jl: An open-source Julia package for efficient and scalable quantum Monte Carlo simulations of electron-phonon coupled models
- Revisiting kinetic electrostatic electron non-linear (KEEN) waves in the presence of dynamical ions
- On the description of chemical bonding using the electrostatic potential. I. Theoretical introduction and covalent interactions
- Atomic Environment Aware Lennard-Jones Parameterization for Electrostatic Embedding ML/MM Simulations
BibTeX
@misc{viand2026efficient,
title = {{Efficient kinetic Monte Carlo simulations with long-range electrostatic interactions}},
author = {Roya Ebrahimi Viand and Bat-Amgalan Bat-Erdene and Karsten Reuter and Sebastian Matera},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.12538},
eprint = {2609.12538},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.12538},
}
RIS
TY - UNPB
TI - Efficient kinetic Monte Carlo simulations with long-range electrostatic interactions
AU - Roya Ebrahimi Viand
AU - Bat-Amgalan Bat-Erdene
AU - Karsten Reuter
AU - Sebastian Matera
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2609.12538
UR - https://arxiv.org/abs/2609.12538
AB - Charge transport in solid-state materials is often governed by charge carrier hopping processes in the presence of long-range electrostatic interactions. Kinetic Monte Carlo (kMC) simulations provide a framework for describing such rare-event dynamics over extended timescales. However, the efficient treatment of long-range interactions remains a major computational challenge, since each particle jump modifies the energy landscape globally and, in principle, requires updating all transition rates after every carrier motion. We present an efficient and generally applicable update process for these transition rates in the presence of long-range electrostatic interactions. To illustrate the method, we study charge diffusion on a simple cubic host lattice under an external electric field. The transport behavior is investigated in high states of charge (SOC), and the influence of temperature, electric field strength, and SOC is examined. Our simulations show strongly suppressed transport at 100\,\% SOC (50\,\% occupation) caused by a freezing of the charge carriers into a Coulomb superlattice. Slight deviations from this reference in terms of charge carrier concentration lead to a rapid increase in conductivity. A deeper analysis reveals that this behavior can be rationalized as transport of non-interacting defects in the Coulomb superlattice. These results demonstrate the capability of the proposed update procedure to efficiently capture non-equilibrium transport phenomena in interacting charged systems and provide a foundation for simulations of more complex charge diffusion problems.
ER -