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.

Advancing X-ray scattering metrology using inverse genetic algorithms

Fitting X-ray diffraction peaks by a Genetic Algorithm and Monte Carlo Markov Chain procedure

Using the Genetic Algorithm to Determine the Parameters of X-Ray Diffractograms Through Different Functions

Automated Materials Spectroscopy Analysis using Genetic Algorithms

Hybrid approach to reconstruct nanoscale grating dimensions using scattering and fluorescence with soft X-rays

Optimization of synchrotron radiation parameters using swarm intelligence and evolutionary algorithms

Obtaining 3D Atomic Reconstructions from Electron Microscopy Images Using a Bayesian Genetic Algorithm: Possibilities, Insights, and Limitations

Reproducible x-ray reflectometry optimization: statistical analysis of differential evolution fitting of multilayer structural models

Powder X-Ray Diffraction Assisted Evolutionary Algorithm for Crystal Structure Prediction

Advances in Kriging-Based Autonomous X-Ray Scattering Experiments

Shape- and Element-Sensitive Reconstruction of Periodic Nanostructures with Grazing Incidence X-ray Fluorescence Analysis and Machine Learning

A genetic algorithm for the atomistic design and global optimisation of substitutionally disordered materials

Small angle x-ray scattering metrology for sidewall angle and cross section of nanometer scale line gratings

Data-driven approach for synchrotron X-ray Laue microdiffraction scan analysis

Robustness of "cut and splice" genetic algorithms in the structural optimization of atomic clusters

A Comparison of Genetic Algorithms and Simulated Annealing in Maximizing the Thermal Conductivity of Discrete Massive Chains

Systematic Comparison of Genetic Algorithm and Basin Hopping Approaches to the Global Optimization of Si(111) Surface Reconstructions

Machine learning for rapid inference of critical dimensions in optical metrology of nanopatterned surfaces

Towards scanning nanostructure X‐ray microscopy

An Improved Real--Space Genetic Algorithm for Crystal Structure and Polymorph Prediction

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