Search papers, labs, and topics across Lattice.
This paper introduces an automated Snakemake workflow for reaction path optimization using the nudged elastic band (NEB) method, addressing the reproducibility issues caused by manual intervention in traditional NEB calculations. The workflow integrates machine learning potentials (PET-MAD) with the eOn saddle point search software, automating endpoint minimization, structural alignment, and path generation. Validation on HCN to HNC isomerization demonstrates the pipeline's ability to reproduce known energy profiles without manual intervention, ensuring identical execution across platforms.
Achieve fully reproducible reaction path optimization with an automated workflow that eliminates manual intervention and ensures consistent results across platforms.
The nudged elastic band (NEB) method is the standard approach for finding minimum energy paths and transition states on potential energy surfaces. Practical NEB calculations require several pre-processing steps: endpoint minimization, structural alignment, and initial path generation. These steps are typically handled by ad-hoc scripts or manual intervention, introducing errors and hindering reproducibility. We present a fully automated, open-source Snakemake workflow for small gas phase molecules that couples modern machine learning potentials (PET-MAD) to the eOn saddle point search software. Each step of the calculation lifecycle is encoded as an explicit dependency graph, from model retrieval and endpoint preparation through path initialization and band optimization. The workflow resolves all software dependencies from conda-forge, ensuring identical execution across platforms. Validation on the HCN to HNC isomerization demonstrates that the automated pipeline recovers the known single-barrier energy profile and product energy without manual intervention.
