Accelerating Transition State Search and Ligand Screening for Organometallic Catalysis with Reactive Machine Learning Potential

Organometallic catalysis lies at the heart of numerous industrial processes that produce bulk and fine chemicals. The search for transition states and screening for organic ligands are vital in designing highly active organometallic catalysts with efficient reaction kinetics. However, identifying accurate transition states necessitates computationally intensive quantum chemistry calculations. In this work, a reactive machine learning potential (RMLP) model is developed to accelerate transition state optimizations and ligand screening for organometallic catalysis based on an automated transition state database construction method and a higher-order equivariant message passing neural network. In case studies involving the ethylene hydrogenation reaction catalyzed by organometallic catalysts, RMLP rapidly predicts potential energy surfaces along intrinsic reaction coordinate paths, achieving speeds nearly 3 orders of magnitude faster than those of rigorous quantum chemistry calculations. Meanwhile, it maintains comparable accuracy with a root-mean-square deviation of 0.307 Å for transition state geometries and a mean absolute error of 0.871 kJ·mol^-1 for reaction barriers on the external test set, significantly outperforming semiempirical quantum chemistry methods. Our RMLP model offers an effective alternative to both rigorous and semiempirical quantum chemistry approaches for rapid and precise transition state optimizations, facilitating high-throughput screening of advanced organometallic catalyst ligands.