Molecular modeling of reactive systems with REACTER

From batteries to biology, many important technologies and physical phenomena operate as out-of-equilibrium reactive systems. Accurately modeling the nanoscale dynamics of non-equilibrium reactive systems and how they respond to external stimuli is challenging, especially if both atomistic resolution and large scales (>105 atoms) are required. REACTER is a protocol for modeling chemical reactions during classical molecular dynamics (MD) simulations. Coupling traditional fixed-valence force fields with heuristic reactive MD is advantageous for large-scale simulations of dynamic systems that can include the complex reaction mechanisms common in organic chemistry. This paper details the current features of the LAMMPS implementation of REACTER, known as fix bond/react, and surveys recent applications of the protocol in a variety of fields, including photopolymers, high-performance composites, and membranes. Conceived as a tool for modeling polymerization processes, the scope of REACTER is expanding as it is applied to new materials and supporting features are implemented. Three new case studies are presented that highlight the capabilities of REACTER, including modeling hierarchical materials, the mechanics of molecular machines, and large-scale dynamics of heterogeneous catalysis.