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.

Revisiting Aspirin Polymorphic Stability Using a Machine Learning Potential

Study of Entropy-Driven Polymorphic Stability for Aspirin Using Accurate Neural Network Interatomic Potential

New Insights into the Crystallographic Disorder in the Polymorphic Forms of Aspirin from Low-Frequency Vibrational Analysis

Studying Mechanical Properties and Phase Transitions of Aspirin Polymorphs with Terahertz Spectroscopy and Ab Initio Simulations

Modern Semiempirical Electronic Structure Methods and Machine Learning Potentials for Drug Discovery: Conformers, Tautomers, and Protonation States

Solid–solid Phase Transformation of Aspirin at High Pressures and Room Temperatures

Revealing the intergrowth phenomenon of aspirin polymorphs through a swift cooling crystallization process

No Headache for PIPs: A PIP Potential for Aspirin Outperforms Other Machine-Learned Potentials

Quantitative Analysis of Physical Stability Mechanisms of Amorphous Solid Dispersions by Molecular Dynamic Simulation

Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol

Mechanistic Insights Into A Classic Wonder Drug-Aspirin

Harnessing Machine Learning for Efficient Large-Scale Interatomic Potential for Sildenafil and pharmaceuticals containing H, C, N, O, and S

A critical examination of compound stability predictions from machine-learned formation energies

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

Ab initio Accuracy Neural Network Potential for Drug-like Molecules

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

A Robust Crystal Structure Prediction Method to Support Small Molecule Drug Development with Large Scale Validation and Prospective Studies

Nucleation Control and Isolation of Polymorphic Forms of Aspirin through an Efficient Template‐Assisted Swift Cooling Process

Modeling the Structure of Ketoprofen-Poly(vinylpyrrolidone) Amorphous Solid Dispersions with Empirical Potential Structure Refinements of X-ray Scattering Data

Formation Energy Prediction of Material Crystal Structures using Deep Learning

Identifying polymer states by machine learning