Multiobjective genetic training and uncertainty quantification of reactive force fields

The ReaxFF reactive force-field approach has significantly extended the applicability of reactive molecular dynamics simulations to a wide range of material properties and processes. ReaxFF parameters are commonly trained to fit a predefined set of quantum-mechanical data, but it remains uncertain how accurately the quantities of interest are described when applied to complex chemical reactions. Here, we present a dynamic approach based on multiobjective genetic algorithm for the training of ReaxFF parameters and uncertainty quantification of simulated quantities of interest. ReaxFF parameters are trained by directly fitting reactive molecular dynamics trajectories against quantum molecular dynamics trajectories on the fly, where the Pareto optimal front for the multiple quantities of interest provides an ensemble of ReaxFF models for uncertainty quantification. Our in situ multiobjective genetic algorithm workflow achieves scalability by eliminating the file I/O bottleneck using interprocess communications. The in situ multiobjective genetic algorithm workflow has been applied to high-temperature sulfidation of MoO3 by H2S precursor, which is an essential reaction step for chemical vapor deposition synthesis of MoS2 layers. Our work suggests a new reactive molecular dynamics simulation approach for far-from-equilibrium chemical processes, which quantitatively reproduces quantum molecular dynamics simulations while providing error bars.