Abstract
We model the autoionization of water by determining the free energy of hydration of the major intermediate species of water ions. We represent the smallest ions─the hydroxide ion OH(-), the hydronium ion H(3)O(+), and the Zundel ion H(5)O(2)(+)─by bonded models and the more extended ionic structures by strong nonbonded interactions (e.g., the Eigen H(9)O(4)(+) = H(3)O(+) + 3(H(2)O) and the Stoyanov H(13)O(6)(+) = H(5)O(2)(+) + 4(H(2)O)). Our models are faithful to the precise QM energies and their components to within 1% or less. Using the calculated free energies and atomization energies, we compute the pK(a) of pure water from first principles as a consistency check and arrive at a value within 1.3 log units of the experimental one. From these calculations, we conclude that the hydronium ion, and its hydrated state, the Eigen cation, are the dominant species in the water autoionization process.