  Classical molecular dynamics (MD) simulations enable modeling of materials
and examination of microscopic details that are not accessible experimentally.
The predictive capability of MD relies on the force field (FF) used to describe
interatomic interactions. FF parameters are typically determined to reproduce
selected material properties computed from density functional theory (DFT)
and/or measured experimentally. A common practice in parameterizing FFs is to
use least-squares local minimization algorithms. Genetic algorithms (GAs) have
also been demonstrated as a viable global optimization approach, even for
complex FFs. However, an understanding of the relative effectiveness and
efficiency of different optimization techniques for the determination of FF
parameters is still lacking. In this work, we evaluate various FF parameter
optimization schemes, using as example a training data set calculated from DFT
for different polymorphs of Ir$O_2$. The Morse functional form is chosen for
the pairwise interactions and the optimization of the parameters against the
training data is carried out using (1) multi-start local optimization
algorithms: Simplex, Levenberg-Marquardt, and POUNDERS, (2) single-objective
GA, and (3) multi-objective GA. Using random search as a baseline, we compare
the algorithms in terms of reaching the lowest error, and number of function
evaluations. We also compare the effectiveness of different approaches for FF
parameterization using a test data set with known ground truth (i.e generated
from a specific Morse FF). We find that the performance of optimization
approaches differs when using the Test data vs. the DFT data. Overall, this
study provides insight for selecting a suitable optimization method for FF
parameterization, which in turn can enable more accurate prediction of material
properties and chemical phenomena.
