Exact distributed kinetic Monte Carlo simulations for on-lattice chemical kinetics: lessons learnt from medium- and large-scale benchmarks

• Published in: Philosophical transactions of the Royal Society of London. Series A: Mathematical, physical, and engineering sciences Vol. 381; no. 2250; p. 20220235
• Main Authors: Savva, Giannis D, Benson, Raz L, Christidi, Ilektra-Athanasia, Stamatakis, Michail
• Format: Journal Article
• Language: English
• Published: England The Royal Society 10.07.2023
• Subjects: catalysis; distributed simulation; materials science; time-warp algorithm; lattice; kinetic Monte Carlo
• ISSN: 1364-503X; 1471-2962; 1471-2962
• DOI: 10.1098/rsta.2022.0235

Kinetic Monte Carlo (KMC) simulations have been instrumental in multiscale catalysis studies, enabling the elucidation of the complex dynamics of heterogeneous catalysts and the prediction of macroscopic performance metrics, such as activity and selectivity. However, the accessible length- and time-scales have been a limiting factor in such simulations. For instance, handling lattices containing millions of sites with 'traditional' sequential KMC implementations is prohibitive owing to large memory requirements and long simulation times. We have recently established an approach for exact, distributed, lattice-based simulations of catalytic kinetics which couples the Time-Warp algorithm with the Graph-Theoretical KMC framework, enabling the handling of complex adsorbate lateral interactions and reaction events within large lattices. In this work, we develop a lattice-based variant of the Brusselator system, a prototype chemical oscillator pioneered by Prigogine and Lefever in the late 60s, to benchmark and demonstrate our approach. This system can form spiral wave patterns, which would be computationally intractable with sequential KMC, while our distributed KMC approach can simulate such patterns 15 and 36 times faster with 625 and 1600 processors, respectively. The medium- and large-scale benchmarks thus conducted, demonstrate the robustness of the approach, and reveal computational bottlenecks that could be targeted in further development efforts. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.