Stable and Accurate Atomistic Simulations of Flexible Molecules using Conformationally Generalisable Machine Learned Potentials

Christopher D Williams,Jas Kalayan,Neil A Burton,Richard A Bryce
DOI: https://doi.org/10.26434/chemrxiv-2024-r75jz
2024-02-16
Abstract:Computational simulation methods based on machine learned potentials (MLPs) promise to revolutionise shape prediction of flexible molecules in solution, but their widespread adoption has been limited by the way in which training data is generated. Here, we present an approach which allows the key conformational degrees of freedom to be properly represented in reference molecular datasets. MLPs trained on these datasets using a global descriptor scheme are generalisable in conformational space, providing quantum chemical accuracy for all conformers. These MLPs are capable of propagating long, stable molecular dynamics trajectories, an attribute that has remained a challenge for MLPs. We deploy the MLPs in obtaining converged conformational free energy surfaces for flexible molecules via well-tempered metadynamics simulations; this approach provides a hitherto inaccessible route to accurately computing the structural, dynamical and thermodynamical properties of a wide variety of flexible molecular systems.
Chemistry
What problem does this paper attempt to address?
The paper discusses the problem of using machine learning potential energy surfaces (MLPs) in molecular dynamics simulations to accurately and stably simulate flexible molecules. Currently, although MLPs have the potential to improve the accuracy of flexible molecular conformation predictions, their widespread application is limited by the way training data is generated. The researchers propose a new approach to ensure that key conformational degrees of freedom are adequately represented in the reference molecular dataset, so that MLPs trained with global descriptors have generalizability in conformational space and can provide quantum chemical accuracy for all conformations in terms of energy and forces. This method can generate long and stable molecular dynamics trajectories and accurately calculate the structural, dynamic, and thermodynamic properties of various flexible molecular systems by obtaining the converged free energy surface of flexible molecules through well-tempered metadynamics simulations. In the paper, the researchers point out the shortcomings of the current training protocol, such as the unstable trajectory problem, and propose solutions, emphasizing the importance of comprehensive conformational sampling for training MLPs.