EReaxFF: A Pseudoclassical Treatment of Explicit Electrons within Reactive Force Field Simulations

We present a computational tool, eReaxFF, for simulating explicit electrons within the framework of the standard ReaxFF reactive force field method. We treat electrons explicitly in a pseudoclassical manner that enables simulation several orders of magnitude faster than quantum chemistry (QC) methods, while retaining the ReaxFF transferability. We delineate here the fundamental concepts of the eReaxFF method and the integration of the Atom-condensed Kohn-Sham DFT approximated to second order (ACKS2) charge calculation scheme into the eReaxFF. We trained our force field to capture electron affinities (EA) of various species. As a proof-of-principle, we performed a set of molecular dynamics (MD) simulations with an explicit electron model for representative hydrocarbon radicals. We establish a good qualitative agreement of EAs of various species with experimental data, and MD simulations with eReaxFF agree well with the corresponding Ehrenfest dynamics simulations. The standard ReaxFF parameters available in the literature are transferrable to the eReaxFF method. The computationally economic eReaxFF method will be a useful tool for studying large-scale chemical and physical systems with explicit electrons as an alternative to computationally demanding QC methods.