Mechanistic Study of Chemical Reactions Using Path Sampling

Abstract

Molecular simulations are the ideal tool to obtain detailed information at the molecular scale which is often invisible to experiments. However, we are still far from a situation in which industries and pharmaceutical laboratories effectively design new materials and medicines based on molecular modeling. Present simulation techniques such as standard molecular dynamics can not reach time- and length-scales required to study complex chemical and biological processes, or they are based on inaccurate and oversimplified models which make them unreliable. Several methods have been proposed over the last decades in order to overcome these problems. Many of those methods are able to provide thermodynamic properties much faster than plain molecular dynamics does, but they force the system in an unnatural way and, therefore, information about the unbiased dynamics is lost. The Transition Path Sampling (TPS) pioneered another approach showing how one can harvest paths (short molecular trajectories) via a statistically sound Monte Carlo approach. Based on this principle TPS also provided an algorithm for the determination of rate constants. Several aspects of this original algorithm were improved by the Transition Interface Sampling (TIS) method. Another variation of TIS, Replica Exchange TIS (RETIS), improved the algorithmic efficiency even further. TIS and RETIS divide the phase space in subregions using interfaces and gather a collection of true dynamical unbiased trajectories connecting the reactant and product states without a priori assumption of the reaction coordinate. From these ensembles of pathways, reaction mechanisms can be extracted and rate constants can be computed.

In this work, for the first time, we use a combination of the RETIS simulations with Ab initio molecular dynamics to study two realistic and challenging reactions. These reactions are autoionization of water and silicate oligomerization reactions. Although these reactions have previously been studied experimentally and/or theoretically, the exact mechanisms of reactions and rate constants are still under discussion. The autodissociation of a water molecule in pure water or an aqueous solution is a fundamental event in acid-base chemistry and is an extremely rare event on the femtosecond time scale of molecular motions. Silicate oligomerization is also a fundamental reaction that is at the basis of sol-gel chemistry and the formation of nanoporous silicate materials. RETIS provides a quantitative and deep qualitative analysis of these fundamental processes. In addition to the RETIS simulations, we have tested a ReaxFF forcefield for the silicate dimerization reactions and have examined this process with Ab initio molecular dynamics simulations using the more conventional thermodynamic integration methods. With these molecular dynamics simulations, we can study both dynamics properties and time dependent kinetic phenomena. Finally, we introduce an approach to analyze collective variables regarding their predictive power for a reaction. The method is based on already available path sampling data produced by, for instance, transition interface sampling or forward flux sampling. ​