Search

The hydroxide anion plays an essential role in many chemical and biochemical reactions. But a molecular-scale description of its hydration state, and hence also its transport, in water is currently controversial. The statistical mechanical quasi-chemical theory of solutions suggests that HO[H2O]3- is the predominant species in the aqueous phase under standard conditions. This result is in close agreement with recent spectroscopic studies on hydroxide water clusters, and with the available thermodynamic hydration free energies. In contrast, a recent ab initio molecular dynamics simulation has suggested that HO[H_2O]4- is the only dominant aqueous solution species. We apply adiabatic ab initio molecular dynamics simulations, and find good agreement with both the quasi-chemical theoretical predictions and experimental results. The present results suggest a picture that is simpler, more traditional, but with additional subtlety. These coordination structures are labile but the tri-coordinate species is the prominent case. This conclusion is unaltered with changes in the electronic density functional. No evidence is found for rate-determining activated inter-conversion of a HO[H2O]4- trap structure to HO[H2O]3-, mediating hydroxide transport. The view of HO- diffusion as the hopping of a proton hole has substantial validity, the rate depending largely on the dynamic disorder of the water hydrogen-bond network., 7 pages, 5 figures, additional results included

Physical Review Letters

This paper describes a mixed direct-iterative method for boundary integral formulations of dielectric solvation models. We give an example for which a direct solution at thermal accuracy is nontrivial and for which Gauss-Seidel iteration diverges in rare but reproducible cases. This difficulty is analyzed by obtaining the eigenvalues and the spectral radius of the iteration matrix. This establishes that the nonconvergence is due to inaccuracies of the asymptotic approximations for the matrix elements for accidentally close boundary element pairs on different spheres. This difficulty is cured by checking for boundary element pairs closer than the typical spatial extent of the boundary elements and for those pairs performing an `in-line' Monte Carlo integration to evaluate the required matrix elements. This difficulty are not expected and have not been observed when only a direct solution is sought. Finally, we give an example application of these methods to deprotonation of monosilicic acid in water., to be published in proceeding of the Pacific Symposium on Biomcomputing; latex with sprocl.sty, 17 pages with 3 figures included with psfig