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.

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

Adaptive Force Field Parameter Optimization for Expanding Reaction Simulations Within Wide-Ranged Temperature

A brand new reactive potential made molecular dynamics on chemical reaction possible

Parameterization of a reactive force field using a Monte Carlo algorithm

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

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

INDEEDopt: a deep learning-based ReaxFF parameterization framework

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

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

Ab-Initio Molecular Dynamics Acceleration Scheme with an Adaptive Machine Learning Framework

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

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

Parametrization of a Reactive Force Field (ReaxFF) for Molecular Dynamics Simulations of Si Nanoparticles

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

Modeling atomic layer deposition of alumina using reactive force field molecular dynamics

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

Machine Learning for Parameter Auto-tuning in Molecular Dynamics Simulations: Efficient Dynamics of Ions near Polarizable Nanoparticles

Global optimization of parameters in the reactive force field ReaxFF for SiOH

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