CatWalk: Reaction Pathway Generation and Transition State Sampling via Score-based Diffusion Models
Abstract
Computing transition states for heterogeneous catalyst surface reactions is computationally demanding due to the relatively high cost in precisely locating a saddle point using ab initio electronic structure calculations. In this work, we present a new generative score-based diffusion model—CatWalk—that predicts the initial structures for reaction pathways of heterogeneous catalysts with surface adsorbates. Given only the desired reaction type and reactant structure, CatWalk produces nine pathway structures by iteratively perturbing the coordinates of the previous output. CatWalk was trained on 731 nudged elastic band (NEB) calculations from the CatTSunami dataset, which contains diverse substrate elements with adsorbates undergoing three types of chemical reactions: desorption, dissociation, and atom transfer. The model can be directed to target one of these three reaction types during generation, and due to its stochastic nature can generate various products via multiple pathways from the same initial state. Optimizing CatWalkgenerated pathways using a pretrained foundation machine-learned interatomic potential (MLIP) with the nudged elastic band (NEB) method leads to several lower-energy pathways compared to what is obtained by the traditional linear interpolation initialization method for NEBs. Additionally, repeated intermediate structure generations using CatWalk enabled sampling of different pathways and final states which can lead to multiple saddle points with variation in energy by up to 0.2 eV. Overall, CatWalk provides a general and scalable framework for accelerating reaction-pathway discovery, enabling broader exploration of catalytic mechanisms at significantly reduced computational cost.
Related
- A reaction volume bias Monte Carlo trial for sampling chemisorption in confinement
- A Scoping Review of Diffusion-Based De Novo Molecular Generation: Architectures, Evaluation, and Open Challenges for Computational Drug Discovery
- Learning continuous reaction paths for transition-state prediction
- Transition states and dynamical bottlenecks in the Li-NC <-> Li-CN isomerizing reaction with the Lagrangian betweenness
- Decoupling CO2 Availability from Bulk Solubility via Biohybrid Reaction-Diffusion Nanoconfinement.
BibTeX
@misc{pakornchote2026catwalk,
title = {{CatWalk: Reaction Pathway Generation and Transition State Sampling via Score-based Diffusion Models}},
author = {Teerachote Pakornchote and Nima Karimitari and Jacob Clary and Ravishankar Sundararaman and Graeme Henkelman and Christopher Sutton},
howpublished = {ChemRxiv},
year = {2026},
month = aug,
doi = {10.26434/chemrxiv.15008101/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15008101/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - CatWalk: Reaction Pathway Generation and Transition State Sampling via Score-based Diffusion Models
AU - Teerachote Pakornchote
AU - Nima Karimitari
AU - Jacob Clary
AU - Ravishankar Sundararaman
AU - Graeme Henkelman
AU - Christopher Sutton
PY - 2026
DA - 2026/08/31
PB - ChemRxiv
DO - 10.26434/chemrxiv.15008101/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15008101/v1
AB - Computing transition states for heterogeneous catalyst surface reactions is computationally demanding due to the relatively high cost in precisely locating a saddle point using ab initio electronic structure calculations. In this work, we present a new generative score-based diffusion model—CatWalk—that predicts the initial structures for reaction pathways of heterogeneous catalysts with surface adsorbates. Given only the desired reaction type and reactant structure, CatWalk produces nine pathway structures by iteratively perturbing the coordinates of the previous output. CatWalk was trained on 731 nudged elastic band (NEB) calculations from the CatTSunami dataset, which contains diverse substrate elements with adsorbates undergoing three types of chemical reactions: desorption, dissociation, and atom transfer. The model can be directed to target one of these three reaction types during generation, and due to its stochastic nature can generate various products via multiple pathways from the same initial state. Optimizing CatWalkgenerated pathways using a pretrained foundation machine-learned interatomic potential (MLIP) with the nudged elastic band (NEB) method leads to several lower-energy pathways compared to what is obtained by the traditional linear interpolation initialization method for NEBs. Additionally, repeated intermediate structure generations using CatWalk enabled sampling of different pathways and final states which can lead to multiple saddle points with variation in energy by up to 0.2 eV. Overall, CatWalk provides a general and scalable framework for accelerating reaction-pathway discovery, enabling broader exploration of catalytic mechanisms at significantly reduced computational cost.
ER -