Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions

• Published in: Journal of chemical theory and computation 2023-07, Vol.19 (13), p.4243-4254
• Main Authors: Huo, Jun, Chen, Jianghao, Liu, Pei, Hong, Benkun, Zhang, Jian, Dong, Hao, Li, Shuhua
• Format: Article
• Language: English
• Subjects: Autoionization; Condensed Matter, Interfaces, and Materials; Data points; Free energy; Molecular dynamics; Neural networks; Protons; Quantum mechanics; Simulation; Water chemistry
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/acs.jctc.3c00244

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 (H3O+) 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 H3O+, as well as a significantly higher free energy barrier for the OH–-associated PT than that of H3O+, 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 H3O+ 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.