Price Trevor, Sivakumar Saurabh, Johnson Matthew S, Zádor Judit, Kulkarni Ambarish
Department of Chemical Engineering, University of California, Davis, California 95616, United States.
Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States.
J Phys Chem C Nanomater Interfaces. 2025 Apr 13;129(16):7751-7761. doi: 10.1021/acs.jpcc.5c00305. eCollection 2025 Apr 24.
Microkinetic models (MKMs) are widely used within the computational heterogeneous catalysis community to investigate complex reaction mechanisms, to rationalize experimental trends, and to accelerate the rational design of novel catalysts. However, constructing these models requires computationally expensive and manually tedious density functional theory (DFT) calculations for identifying transition states for each elementary reaction within the MKM. To address these challenges, we demonstrate a novel protocol that uses the open-source kinetics workflow tool Pynta to automate the iterative training of a reactive machine learning potential (rMLP). Specifically, using the silver-catalyzed partial oxidation of methanol as a prototypical example, we first demonstrate our workflow by training an rMLP to accelerate the parallel calculation of DFT-quality transition states for all 53 reactions, achieving a 7× speedup compared to a DFT-only strategy. Detailed analysis of our training curriculum reveals the shortcomings of using an adaptive sampling scheme with a single rMLP model to describe all reactions within the MKM simultaneously. We show that these limitations can be overcome using a balanced "reaction class" approach that uses multiple rMLP models, each describing a single class of similar transition states. Finally, we demonstrate that our Pynta-based workflow is also compatible with large pretrained foundational models. For example, by fine-tuning a top-performing graph neural network potential trained on the OC20 dataset, we observe an impressive 20× speedup with an 89% success rate in identifying transition states. This work highlights the synergistic potential of integrating automated tools with machine learning to advance catalysis research.
微观动力学模型(MKMs)在计算多相催化领域中被广泛应用,用于研究复杂的反应机理、解释实验趋势以及加速新型催化剂的合理设计。然而,构建这些模型需要进行计算成本高昂且人工繁琐的密度泛函理论(DFT)计算,以确定MKMs中每个基元反应的过渡态。为应对这些挑战,我们展示了一种新颖的协议,该协议使用开源动力学工作流程工具Pynta来自动迭代训练反应性机器学习势(rMLP)。具体而言,以银催化甲醇部分氧化作为典型示例,我们首先通过训练一个rMLP来加速对所有53个反应的DFT质量过渡态的并行计算,展示了我们的工作流程,与仅使用DFT的策略相比,实现了7倍的加速。对我们训练过程的详细分析揭示了使用单个rMLP模型的自适应采样方案来同时描述MKMs中的所有反应的缺点。我们表明,使用平衡的“反应类”方法,即使用多个rMLP模型,每个模型描述一类相似的过渡态,可以克服这些限制。最后,我们证明基于Pynta的工作流程也与大型预训练基础模型兼容。例如,通过微调在OC20数据集上训练的性能最佳的图神经网络势,我们在识别过渡态时观察到了令人印象深刻的20倍加速,成功率达到89%。这项工作突出了将自动化工具与机器学习相结合以推进催化研究的协同潜力。