Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions

Abstract

Water molecules and the associated proton transfer (PT) are prevalent in chemical and biological systems and have been a hot research topic. Spectroscopic characterization and ab initio molecular dynamics (AIMD) simulations have previously revealed insights into acidic and basic liquids. Presumably, the situation in the acidic/basic solution is not necessarily the same as in pure water; in addition, the autoionization constant for water is only 10 −14 under ambient conditions, making the study of PT in pure water challenging. To overcome this issue, we modeled periodic water box systems containing 1000 molecules for tens of nanoseconds based on a neural network potential (NNP) with quantum mechanical accuracy. The NNP was generated by training a dataset containing the energies and atomic forces of 17 075 configurations of periodic water box systems, and these data points were calculated at the MP2 level that considers electron correlation effects. We found that the size of the system and the duration of the simulation have a significant impact on the convergence of the results. With these factors considered, our simulations showed that hydronium (H 3 O + ) and hydroxide (OH − ) ions in water have distinct hydration structures, thermodynamic and kinetic properties, e.g., the longer-lasting and more stable hydrated structure of OH − ions than that of H 3 O + , as well as a significantly higher free energy barrier for the OH − -associated PT than that of H 3 O + , leading the two to exhibit completely different PT behaviors. Given these characteristics, we further found that PT via OH − ions tends not to occur multiple times or between many molecules. In contrast, PT via H 3 O + can synergistically occur among multiple molecules and prefers to adopt a cyclic pattern among three water molecules, while it occurs mostly in a chain pattern when more water molecules are involved. Therefore, our studies provide a detailed and solid microscopic explanation for the PT process in pure water.

Reshaping the QM Region On-the-Fly: Adaptive-Shape QM/MM Dynamic Simulations of a Hydrated Proton in Bulk Water

“…The curve for g D–H (Figure b), which corresponds to the pair of the hydronium oxygen and the water hydrogen, has two distinct peaks, which are located at r = 1.0 and 3.1 Å, respectively. Overall, our results are in good agreement with previous computational studies. ,,− …”