Data-driven massive reaction networks reveal mechanistic pathways underlying catalytic CO2 hydrogenation

Heterogeneous catalytic pathways for clean energy conversion involve thousands of elementary steps, but most models involve only a few dozen reactions. We combine extensive density functional theory calculations, machine learning for activation barrier prediction, and human intelligence-inspired elementary reaction enumeration and identification. This enables automated kinetic modeling of CO2 hydrogenation on copper, a key process to produce fuels and chemicals. We construct one of the largest datasets consisting of 152 elementary CO2 hydrogenation reactions and experimentally determine CO2 conversion, finding that even large networks with 100+ reactions are insufficient. In contrast, our approach reveals 9389 elementary reactions, reducing human bias in the reaction pathway. We unravel 40-fold higher CO2 conversion rates, following experimental trends of methanol and CO production. We establish the crucial role of inter-species hydrogen transfer and hydrogenation by molecular H2, a surprising data-driven discovery validated post-facto. The proposed strategy to comprehensively model complex reaction mechanisms will significantly advance catalysis research. Modeling complex heterogeneous catalytic networks with DFT alone is too expensive. Combining DFT with machine learning and elementary reaction enumeration, the proposed framework enables kinetic modeling of thousands of elementary reactions for CO2 reduction on copper, successfully validating experimentally observed trends.