Abstract:The growing interest in machine learning (ML) tools within chemistry and material science stems from their novelty and ability to predict properties almost as accurately as underlying electronic structure calculations or experiments. Transition path sampling (TPS) offers a practical way to explore transition routes between metastable minima such as conformers and isomers on the multidimensional potential energy surface. However, TPS has historically suffered from the computational cost vs. accuracy trade-off between affordable force-field simulations and expensive high-fidelity quantum mechanical calculations. ML interatomic potentials combined with TPS offer a new approach for the exploration of transition pathways at near-quantum mechanical accuracy, while keeping the computational cost comparable to classical force fields. In this study, we employ the HIP-NN-TS and ANI-1x neural network-based ML potentials, both trained on the ANI-1x dataset of 5 million HCNO structures. We first verify the correctness of our approach by applying it to alanine dipeptide and compare the resulting energy surface and transition paths to the literature. Our findings suggest that proposed approach holds promise for conformational searches, as evidenced by the chemical accuracy (errors ≲ 1 kcal/mol) for thermal molecular dynamics trajectories of alanine dipeptide. While we were able to successfully reconstruct alanine dipeptide’s potential landscape using both HIP-NN-TS and ANI-1x frameworks, we observed that ML models with lower accuracy may still recover additional important conformations. We also find that active learning, augmenting the training data by structures taken from TPS trajectories, improved the accuracy by ~30% with small amounts of additional data. Finally, we evaluated a more intricate case, azobenzene, and reiterated that seemingly simple torsions may bear a challenge for ML potentials and limit their applications in TPS. Inability of HIP-NN-TS to correctly describe the energetics of major rotational pathway in azobenzene isomerization highlights deficiencies of the reference method in describing the electronic degrees of freedom. Our study underscores the importance of domain expertise in selecting physically meaningful pathways for benchmarking ML potentials, especially considering the intricacies of electronic structure in chemical dynamics and non-equilibrium processes.

Investigating the hyperparameter space of deep neural network models for reaction coordinates: Revisiting the solvent coordinate in alanine dipeptide isomerization

Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide

Unveiling interatomic distances influencing the reaction coordinates in alanine dipeptide isomerization: An explainable deep learning approach

Capturing Dichotomic Solvent Behavior in Solute--Solvent Reactions with Neural Network Potentials

Development of a New Parameter Optimization Scheme for a Reactive Force Field (ReaxFF) Based on a Machine Learning Approach

Reactive Dynamics and Spectroscopy of Hydrogen Transfer from Neural Network-Based Reactive Potential Energy Surfaces

Machine Learning-Guided Adaptive Parametrization for Coupling Terms in a Mixed United-Atom/Coarse-Grained Model for Diphenylalanine Self-Assembly in Aqueous Ionic Liquids

Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach

Deep Learning Reaction Network: a machine learning framework for modeling time resolved data

Artificial Neural Network-Based Hamiltonian Monte Carlo for High-Dimensional Bayesian Inference of Reaction Kinetics Models

Multitask Machine Learning of Collective Variables for Enhanced Sampling of Rare Events

ReactAIvate: A Deep Learning Approach to Predicting Reaction Mechanisms and Unmasking Reactivity Hotspots

Modeling of water gas shift reaction using neural network trained on detailed kinetic mechanisms

Machine learning based implicit solvent model for aqueous-solution alanine dipeptide molecular dynamics simulations

Low-dimensional projection of reactivity classes in chemical reaction dynamics using supervised dimensionality reduction

Neural Network Based ∆-Machine learning approach efficiently brings the DFT potential energy surface to the CCSD(T) quality: a case for the OH + CH3OH reaction

Reaction Coordinates by Nonlinear Dimensionality Reduction

The Neural Network Based Δ-machine Learning Approach Efficiently Brings the DFT Potential Energy Surface to the CCSD(T) Quality: a Case for the OH + CH3OH Reaction.

Challenges and Opportunities for Machine Learning Potentials in Transition Path Sampling: Alanine Dipeptide and Azobenzene Studies

Machine Learning Identifies Chemical Characteristics That Promote Enzyme Catalysis

Generic Interpretable Reaction Condition Predictions with Open Reaction Condition Datasets and Unsupervised Learning of Reaction Center