Water autoionization reaction 2H2O 3 H3O ؊ ؉ OH ؊ is a textbook process of basic importance, resulting in pH ؍ 7 for pure water. However, pH of pure water surface is shown to be significantly lower, the reduction being caused by proton stabilization at the surface. The evidence presented here includes ab initio and classical molecular dynamics simulations of water slabs with solvated H3O ؉ and OH ؊ ions, density functional studies of (H2O)48H ؉ clusters, and spectroscopic isotopic-exchange data for D2O substitutional impurities at the surface and in the interior of ice nanocrystals. Because H3O ؉ does, but OH ؊ does not, display preference for surface sites, the H2O surface is predicted to be acidic with pH < 4.8. For similar reasons, the strength of some weak acids, such as carbonic acid, is expected to increase at the surface. Enhanced surface acidity can have a significant impact on aqueous surface chemistry, e.g., in the atmosphere.density functional theory ͉ IR spectroscopy ͉ molecular dynamics ͉ water autoionization ͉ ice nanocrystals I n room-temperature liquid, one in 6 ϫ 10 8 water molecules is autoionized, yielding the standard value of pH ϭ 7. Autoionization in crystal ice should be less favorable, because, in contrast to water, ice is a very poor solvent of ionic and polar substances (1). As recently realized (2-5), the chemistry and composition of aqueous surfaces are quite distinct from that of the bulk; therefore, autoionization behavior should be reexamined at the surface.A number of recent computations (6-8) indicated the preference of hydronium cations for surface positions. Surface propensity of H 3 O ϩ was also deduced from vibrational spectroscopy of large protonated water clusters (6), as well as vibrational sum frequency generation (8, 9) and second harmonic generation (10) spectroscopic experiments probing extended aqueous interfaces. Interestingly, zeta potential measurements and titration experiments on oil droplets dispersed in water indicated the presence of negative charges at the interface, interpreted as adsorbed OH Ϫ ions (11). Similar conclusions have also been drawn from zeta potentials of air bubbles in water (12). More work is clearly needed to reconcile this apparent discord between predictions of surface-selective spectroscopies and molecular simulations on one side and electrochemical measurements on the other side.H 3 O ϩ forms three strong proton-donor bonds to H 2 O, but acts as a poor proton acceptor. A surface position with only H atoms hydrogen-bonded is preferred to interior positions, because the latter are associated with disruption of the approximately tetrahedral hydrogen-bond network in water (10). The present work focuses on the effect of surface stabilization of hydronium on water autoionization and surface pH.
CalculationsOverview. Modeling of proton-tranfer systems is a nontrivial problem, because standard (empirical) potential energy surfaces do not include a possibility of proton hopping between different water molecules or transitions between the two li...