Ab initio molecular dynamics simulation of the solvation and transport of hydronium and hydroxyl ions in water

Charge defects in water created by excess or missing protons appear in the form of solvated hydronium H3O+ and hydroxyl OH− ions. Using the method of ab initio molecular dynamics, we have investigated the structure and proton transfer dynamics of the solvation complexes, which embed the ions in the network of hydrogen bonds in the liquid. In our ab initio molecular dynamics approach, the interatomic forces are calculated each time step from the instantaneous electronic structure using density functional methods. All hydrogen atoms, including the excess proton, are treated as classical particles with the mass of a deuterium atom. For the H3O+ ion we find a dynamic solvation complex, which continuously fluctuates between a (H5O2)+ and a (H9O4)+ structure as a result of proton transfer. The OH− has a predominantly planar fourfold coordination forming a (H9O5)− complex. Occasionally this complex is transformed in a more open tetrahedral (H7O4)− structure. Proton transfer is observed only for the more waterlike (H7O4)− complex. Transport of the charge defects is a concerted dynamical process coupling proton transfer along hydrogen bonds and reorganization of the local environment. The simulation results strongly support the structural diffusion mechanism for charge transport. In this model, the entire structure—and not the constituent particles—of the charged complex migrates through the hydrogen bond network. For H3O+, we propose that transport of the excess proton is driven by coordination fluctuations in the first solvation shell (i.e., second solvation shell dynamics). The rate-limiting step for OH− diffusion is the formation of the (H7O4)− structure, which is the solvation state showing proton transfer activity.