Track papers across labs, topics, and trends.
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 associated with manual pre-processing steps. The workflow integrates machine learning potentials (PET-MAD) with the eOn saddle point search software, automating endpoint minimization, structural alignment, and initial path generation. Validation on the HCN to HNC isomerization demonstrates the pipeline's ability to recover known energy profiles without manual intervention, ensuring identical execution across platforms.
Say goodbye to ad-hoc scripts: this automated workflow slashes manual intervention in NEB calculations, ensuring reproducible reaction path optimization 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.