Hydrogen Bond Dynamics in the Solvation Shell on Proton Transfer in Aqueous Solution

Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions

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₃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₃O⁺, as well as a significantly higher free energy barrier for the OH^--associated PT than that of H₃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₃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.

Subpicosecond Molecular Rearrangements Affect Local Electric Fields and Auto-Dissociation in Water

Molecular simulations of auto-dissociation of water molecules in an 81,000 atom bulk water system show that the electric field variations caused by local bond length and angle variations enhance proton transfer within ∼600 fs prior to auto-dissociation. In this paper, auto-dissociation relates to the initial separation of a proton from a water molecule to another, forming the H33O+ and OH- ions. Only transfers for which a proton's initial nearest covalently bonded oxygen remained the same for at least 1 ps prior to the transfer and for which that proton's new nearest acceptor oxygen remained the same for at least 1 ps after the transfer were evaluated. Electric fields from solvent atoms within 6 Å of a transferring proton (H*) are dominant, with little contribution from farther molecules. However, exclusion of the accepting oxygen in such electric field calculations shows that the field on H* from the other solvent atoms weakens as the time to transfer becomes less than 600 fs, indicating the primary importance of the accepting oxygen on enabling auto-dissociation. All resultant OH- and H3O+ ion pairs recombined at times greater than 1 ps after auto-dissociation. A concentration of 8.01 × 1017 cm-3 for these ion pairs was observed. The simulations indicate that transient auto-dissociation in water is more common than that inferred from dc-conductivity experiments (10-5 vs 10-7) and is consistent with the results of calculations that include nuclear quantum effects. The conductivity experiments require the rearrangement of farther water molecules to form hydrogen-bonded "water wires" that afford long-range and measurable proton transport away from the reaction site. Nonetheless, the relatively large number of picosecond-lived auto-dissociation products might be engineered within 2D layers and oriented external fields to offer new energy-related systems.