Connecting the quantum and classical mechanics simulation world: Applications of reactive step molecular dynamics simulations

J Chem Phys. 2021 May 21;154(19):194105. doi: 10.1063/5.0048618.

Abstract

This article presents the application of the reactive step molecular dynamics simulation method [M. Biedermann, D. Diddens, and A. Heuer, J. Chem. Theory Comput. 17, 1074 (2021)] toward two different atomistic, chemically reactive systems. During reactive steps, transitions from reactant to product molecules are modeled according to physically correct transition probabilities based on quantum chemical information about the reactions such as molecular reaction rates via instant exchange of the employed force field and a subsequent, short relaxation of the structure. In the first application, we study the follow-up reactions of singly reduced ethylene carbonate (EC) radicals in EC solution, first, via extensive ab initio molecular dynamics simulations and, second, with the reactive step algorithm. A direct comparison of both simulation methods shows excellent agreement. Then, we employ the reactive step algorithm to simulate the enolate formation of 2-methylcyclopropanone with the base lithium diisopropylamine. Thereby, we can demonstrate that the reactive step algorithm is also capable of capturing effects from kinetic vs thermodynamic control of chemical reactions during simulation.