Protons in Supercritical Water:  A Multistate Empirical Valence Bond Study

Molecular dynamics simulations have been performed to analyze microscopic details related to aqueous solvation of excess protons along the supercritical T = 673 K isotherm, spanning a density interval from a typical liquid down to vapor environments. The simulation methodology relies on a multistate empirical valence bond Hamiltonian model that includes a proton translocation mechanism. Our results predict a gradual stabilization of the solvated Eigen cation H(3)O.(H(2)O)(3)(+) at lower densities, in detriment of the symmetric Zundel dimer H.(H(2)O)(2)(+). At all densities, the average solvation structure in the close vicinity of the hydronium is characterized by three hydrogen bond acceptor water molecules and presents minor changes in the solute water distances. Characteristic times for the proton translocation jumps have been computed using population relaxation time correlation functions. Compared to room temperature results, the rates at high densities are 4 times faster and become progressively slower in steamlike environments. Diffusion coefficients for the excess proton have also been computed. In agreement with conductometric data, our results show that contributions from the Grotthus mechanism to the overall proton transport diminish at lower densities and predict that in steamlike environments, the proton diffusion is almost 1 order of magnitude slower than that for pure water. Spectroscopic information for the solvated proton is accordant to the gradual prevalence of proton localization in Eigen-like structures at lower densities.