Abstract:van der Waals dispersion-repulsion interactions, commonly represented in atomistic force fields by the Lennard-Jones (LJ) potential, play an important role in the accuracy of molecular simulations. Training the force field parameters used in the LJ potential is challenging, generally requiring adjustment based on simulations of macroscopic physical properties. The large computational expense of these simulations, especially if many parameters are trained simultaneously, limits the size of training data set and number of optimization steps that can be taken, often requiring modelers to perform optimizations within a local parameter region. To allow for more global LJ parameter optimization against large training sets, we introduce a multi-fidelity optimization technique which uses Gaussian process surrogate modeling to build inexpensive models of physical properties as a function of LJ parameters. This allows for fast evaluation of objective functions, greatly accelerating searches over parameter space and enabling the use of global optimization algorithms. We use an iterative framework which performs optimization with differential evolution at the surrogate level, followed by validation at the simulation level and surrogate refinement. Using this technique on two previously studied training sets, containing up to 195 physical property targets, we refit a subset of the LJ parameters for the OpenFF 1.0.0 (Parsley) force field. We demonstrate that this multi-fidelity technique can find improved parameter sets compared to a purely simulation-based optimization by searching more broadly and escaping local minima. In particular, this technique often finds significantly different parameter minima that have comparably accurate performance. In most cases, these parameter sets are transferable to other similar molecules in a test set. This multi-fidelity technique provides a platform for fast optimization against physical properties that can be refined and applied in multiple ways to the development of molecular models.

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

Stochastic modeling of non-linear adsorption with Gaussian kernel density estimators

An Adaptive Gaussian Process‐based Method for Efficient Bayesian Experimental Design in Groundwater Contaminant Source Identification Problems

Accelerating Black-Box Molecular Property Optimization by Adaptively Learning Sparse Subspaces

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

A Bayesian approach for parameter estimation and prediction using a computationally intensive model

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Parameter optimisation using Bayesian inference for spallation models

Fast Bayesian Force Fields from Active Learning: Study of Inter-Dimensional Transformation of Stanene

Accelerating inference for stochastic kinetic models

Fast, accurate and lightweight sequential simulation-based inference using Gaussian locally linear mappings

Additive Multi-Index Gaussian process modeling, with application to multi-physics surrogate modeling of the quark-gluon plasma

Iterative Construction of Gaussian Process Surrogate Models for Bayesian Inference

Bayesian optimization of atomic structures with prior probabilities from universal interatomic potentials

Hamiltonian Monte Carlo Acceleration Using Surrogate Functions with Random Bases

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Using physical property surrogate models to perform accelerated multi-fidelity optimization of force field parameters

Fast and accurate nonadiabatic molecular dynamics enabled through variational interpolation of correlated electron wavefunctions