Methane steam reforming is an important industrial process
for
hydrogen production, employing Ni as a low-cost, highly active catalyst,
which, however, suffers from coking due to methane cracking. Coking
is the accumulation of a stable poison over time, occurring at high
temperatures; thus, to a first approximation, it can be treated as
a thermodynamic problem. In this work, we developed an Ab initio kinetic
Monte Carlo (KMC) model for methane cracking on Ni(111) at steam reforming
conditions. The model captures C–H activation kinetics in detail,
while graphene sheet formation is described at the level of thermodynamics,
to obtain insights into the “terminal (poisoned) state”
of graphene/coke within reasonable computational times. We used cluster
expansions (CEs) of progressively higher fidelity to systematically
assess the influence of effective cluster interactions between adsorbed
or covalently bonded C and CH species on the “terminal state”
morphology. Moreover, we compared the predictions of KMC models incorporating
these CEs into mean-field microkinetic models in a consistent manner.
The models show that the “terminal state” changes significantly
with the level of fidelity of the CEs. Furthermore, high-fidelity
simulations predict C–CH island/rings that are largely disconnected
at low temperatures but completely encapsulate the Ni(111) surface
at high temperatures.
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’.
We extend the work of Ravipati et al.[Comput. Phys. Commun., 2022, 270, 108148] in benchmarking the performance of large-scale, distributed, on-lattice kinetic Monte Carlo (KMC) simulations. Our software package, Zacros,...
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.