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.

The automated optimisation of a coarse-grained force field using free energy data

Utilizing Machine Learning for Efficient Parameterization of Coarse Grained Molecular Force Fields

Machine Learning of coarse-grained Molecular Dynamics Force Fields

Machine Learning Directed Optimization of Classical Molecular Modeling Force Fields

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

Empirical optimization of molecular simulation force fields by Bayesian inference

Force field optimization by end-to-end differentiable atomistic simulation

Optimizing Force Fields with Experimental Data Using Ensemble Reweighting and Potential Contrasting

Machine learned coarse-grained protein force-fields: Are we there yet?

Differentiable molecular simulation can learn all the parameters in a coarse-grained force field for proteins

A Neural-Network-Based Mapping and Optimization Framework for High-Precision Coarse-Grained Simulation

Open Force Field BespokeFit: Automating Bespoke Torsion Parametrization at Scale

A Scalable Graph Neural Network Method For Developing An Accurate Force Field Of Large Flexible Organic Molecules

Top-down machine learning of coarse-grained protein force-fields

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field

Statistically Optimal Force Aggregation for Coarse-Graining Molecular Dynamics

Fine-Tuning Property Domain Weighting Factors and the Objective Function in Force Field Parameter Optimization

Two for One: Diffusion Models and Force Fields for Coarse-Grained Molecular Dynamics

Synthetic Force-Field Database for Training Machine Learning Models to Predict Mobility-Preserving Coarse-Grained Molecular-Simulation Potentials

Optimization of Reactive Force Field Simulation: Refactor, Parallelization, and Vectorization for Interactions

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