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.

System-specific parameter optimization for non-polarizable and polarizable force fields

System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

Construction of atomistic force field model for perfluorosulfonic acid polymer

Broadening access to small-molecule parameterization with the force field toolkit

FFParam-v2.0: A Comprehensive Tool for CHARMM Additive and Drude Polarizable Force-Field Parameter Optimization and Validation

Polarizable Force Fields Based On Physical Models And Quantum Chemical Calculations

Development of a Polarizable Force Field Using Multiple Fluctuating Charges per Atom

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

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

A Scalable Molecular Force Field Parameterization Method Based on Density Functional Theory and Quantum-Level Machine Learning

Efficient and Precise Force Field Optimization for Biomolecules Using DPA-2

PCMRESP: A Method for Polarizable Force Field Parameter Development and Transferability of the Polarizable Gaussian Multipole Models Across Multiple Solvents

Balancing Group 1 Monoatomic Ion-Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems

Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation

Polarizable Force Field of Intrinsically Disordered Proteins with CMAP and Reweighting Optimization