(142b) Machine-Learning Enabled Optimization of Force Fields
AIChE Spring Meeting and Global Congress on Process Safety
Thursday, April 22, 2021 - 1:50pm to 2:10pm
In this work we propose a machine learning-enabled automated force field optimization framework. Specifically, we show integrating Gaussian Process (GP) regression (e.g., surrogate) models and support vector machine (SVM) classifiers facilitates rapid tuning of force fields and provides a quick and efficient route to highly accurate, physics-based molecular models. We train a GP surrogate model on the results of molecular simulations so that for a given set of parameters, our GP model predicts the resulting simulation experimental property prediction. We explore different methods of constructing our GP model, including various kernel and mean functions, to determine how to best harness a GP modelâs ability to identify optimal regions of force field model parameter space and obtain the most reliable simulation result predictions. We also examine the use of a SVM classifier to capture discontinuities within experimental properties which a GP model cannot replicate.
As a demonstration case, we optimize force fields for two hydrofluorocarbons (HFCs), HFC-32 and HFC-125, for properties including liquid and vapor densities, vapor pressure, and enthalpy of vaporization. Results show we can find at least 26 HFC-32 and 45 HFC-125 force field parameter sets in a timeframe of weeks which give mean absolute percent error in all of the properties of interest of at most 5%. Additionally, we find that these parameter sets are able to predict transport and critical properties accurately for HFC-32 and HFC-125 without the need for further tuning. We also applied our framework to tuning an ammonium perchlorate force field to predict solid properties, including lattice parameters, unit cell structure, and Hydrogen bond distances, angles, and symmetry. Multiple parameter sets have been found that outperform existing force fields in reproducing experimental observations of the listed quantities. Future work involves expanding this tuning method to encompass additional thermodynamic and transport properties of interest in an automated parameterization workflow. As further capabilities are developed, we envision this tool will be integrated within computer aided molecular and process design schemes to facilitate rapid multiscale design and optimization.
Wang, J., & Kollman, P. A. (2001). Automatic parameterization of force field by systematic search and genetic algorithms. Journal of Computational Chemistry, 22(12), 1219-1228.