FFLUX molecular simulations driven by atomic Gaussian process regression models

Yulian T. Manchev,Paul L. A. Popelier
DOI: https://doi.org/10.1002/jcc.27323
2024-02-14
Journal of Computational Chemistry
Abstract:Gaussian process regression (GPR) models of topological atoms, defined by Quantum chemical topology (QCT), are created via GPyTorch and interfaced with FFLUX, a next‐generation machine learning force field, for the first time. The models are rigorously benchmarked on energy and notably force prediction accuracy, as well as geometry optimization performance in FFLUX. An improved kernel function which captures the periodicity of the input descriptors is used in the models and simulations. Machine learning (ML) force fields are revolutionizing molecular dynamics (MD) simulations as they bypass the computational cost associated with ab initio methods but do not sacrifice accuracy in the process. In this work, the GPyTorch library is used to create Gaussian process regression (GPR) models that are interfaced with the next‐generation ML force field FFLUX. These models predict atomic properties of different molecular configurations that appear in a progressing MD simulation. An improved kernel function is utilized to correctly capture the periodicity of the input descriptors. The first FFLUX molecular simulations of ammonia, methanol, and malondialdehyde with the updated kernel are performed. Geometry optimizations with the GPR models result in highly accurate final structures with a maximum root‐mean‐squared deviation of 0.064 Å and sub‐kJ mol−1 total energy predictions. Additionally, the models are tested in 298 K MD simulations with FFLUX to benchmark for robustness. The resulting energy and force predictions throughout the simulation are in excellent agreement with ab initio data for ammonia and methanol but decrease in quality for malondialdehyde due to the increased system complexity. GPR model improvements are discussed, which will ensure the future scalability to larger systems.
chemistry, multidisciplinary
What problem does this paper attempt to address?