Existing force field optimization methods for ReaxFF employ gradient-free blackbox optimization methods such as Genetic Algorithms (GA) and Evolutionary Algorithms (EA). These methods perform a global search in a high dimensional space and come with a high computational cost because they do not utilize any gradient information, but rather rely on error evaluations at different points in the search space.