A Multipole-Based Reactive Force Field for Hydrocarbons

J Chem Theory Comput. 2024 Nov 26;20(22):10045-10058. doi: 10.1021/acs.jctc.4c01285. Epub 2024 Nov 4.

Abstract

The computational complexity of quantum chemistry methods has prompted the development of reactive force fields, facilitating practical applications of molecular dynamics simulations for large-scale reactive systems. Current reactive force fields typically employ intricate corrections based on prior chemical knowledge, which severely impedes their further advancement. This study presents a new atomic multipole-based reactive model with bond free (OPERATOR). The force field is constructed on a simple, physically motivated model within the AMOEBA framework that closely resembles the physical representation of the chemical reaction processes. In the force field, the atomic multipoles are generated dynamically according to the atomic environments, aiming to effectively capture significant changes in the electrostatic environments during chemical reactions. Subsequently, atomic multipole-based charge penetration, polarization, and charge transfer effects are incorporated into the force field to describe the complex electrostatic interactions in the system. The force field also includes van der Waals interactions and three-body potentials. In addition, to extend these nonreactive interactions to chemical reactions, the atom distribution multipole moments are used to characterize different chemical environments. The force field has been optimized using the dataset of potential energy surfaces (PESs) of hydrocarbons derived from DFT results of millions of conformations with six degrees of freedom (DOFs). The results demonstrate that the new force field effectively replicates both the monopoles and the energies. In comparison to ReaxFF, the new force field exhibits comparable or superior performance. Furthermore, molecular dynamics simulations of n-heptane decomposition effectively reproduce the primary products and reactions observed in the experiments. Given the simplicity and physically motivated nature of the model, it is expected that the new force field will be utilized in future studies to investigate chemical reaction mechanisms involving more elements.