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 timescales are still limited, and handling lattices containing millions of sites with “traditional” sequential KMC implementations becomes prohibitive due to large memory requirements and long simulation times. [more]