Abstract:ReaxFF is a computationally efficient force field to simulate complex reactive dynamics in extended molecular models with diverse chemistries, if reliable force-field parameters are available for the chemistry of interest. If not, they must be calibrated by minimizing the error ReaxFF makes on a relevant training set. Because this optimization is far from trivial, many methods, in particular genetic algorithms (GAs), have been developed to search for the global optimum in parameter space. Recently, two alternative parameter calibration techniques were proposed, i.e.\ Monte-Carlo Force Field optimizer (MCFF) and Covariance Matrix Adaptation Evolutionary Strategy (CMA-ES), which have the potential to find good parameters at a relatively low computational cost. In this work, these two methods are tested, as implemented in ADF2018, using three ReaxFF training sets, which have previously been used to benchmark GAs. Even though MCFF and CMA-ES should not be considered as exhaustive global optimizers, they can find parameters that are comparable in quality to those obtained with GAs. We observe that CMA-ES leads to slightly better results and is less sensitive to the initial guess of the parameters. Concrete recipes are provided for obtaining similar results with new training sets.Besides optimization recipes, a successful ReaxFF parameterization requires the design of a good training set. At every trial set of parameters, ReaxFF is used to optimize molecular geometries in the training set. When the optimization of some geometries fails easily, it becomes increasingly difficult to find the optimal parameters. We have addressed this issue by fixing several bugs in the ReaxFF forces and by improving the robustness of the geometry optimization. These improvements cannot eliminate all geometry convergence issues and we recommend to avoid very flexible geometries in the training set.Both MCFF and CMA-ES are still liable to converge to sub- or near-optimal parameters, which we detected by repeating the calibration with different random seeds. The existence of distinct near-optimal parameter vectors is a general pattern throughout our study and provides opportunities to improve the training set or to detect overfitting artifacts.

An Initial Design-enhanced Deep Learning-based Optimization Framework to Parameterize Multicomponent ReaxFF Force Fields

INDEEDopt: a deep learning-based ReaxFF parameterization framework

Parameterization of a reactive force field using a Monte Carlo algorithm

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

Development of a New Parameter Optimization Scheme for a Reactive Force Field (ReaxFF) Based on a Machine Learning Approach

Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach

ReaxFF Parameter Optimization with Monte-Carlo and Evolutionary Algorithms: Guidelines and Insights

JAX-ReaxFF: A Gradient-Based Framework for Fast Optimization of Reactive Force Fields.

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

Machine Learning Directed Optimization of Classical Molecular Modeling Force Fields

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

Systematic Comparison of Monte Carlo Annealing and Covariance Matrix Adaptation for the Optimization of ReaxFF Parameters

General Multiobjective Force Field Optimization Framework, with Application to Reactive Force Fields for Silicon Carbide

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

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

Enhancing ReaxFF for Molecular Dynamics Simulations of Lithium-Ion Batteries: An interactive reparameterization protocol

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

Managing Expectations and Imbalanced Training Data in Reactive Force Field Development: an Application to Water Adsorption on Alumina

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

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