A new interatomic potential for dissociative water was developed for use in molecular dynamics simulations. The simulations use a multibody potential, with both pair and three-body terms, and the Wolf summation method for the long-range Coulomb interactions. A major feature in the potential is the change in the shortrange O-H repulsive interaction as a function of temperature and/or pressure in order to reproduce the densitytemperature curve between 273 K and 373 at 1 atm, as well as high-pressure data at various temperatures. Using only the change in this one parameter, the simulations also reproduce room-temperature properties of water, such as the structure, cohesive energy, diffusion constant, and vibrational spectrum, as well as the liquid-vapor coexistence curve. Although the water molecules could dissociate, no dissociation is observed at room temperature. However, behavior of the hydronium ion was studied by introduction of an extra H + into a cluster of water molecules. Both Eigen and Zundel configurations, as well as more complex configurations, are observed in the migration of the hydronium.

Molecular dynamics (MD) computer simulation of the adsorption of water molecules onto the vitreous silica surface was performed using a new dissociative water potential. 58 The simulations showed dissociative chemisorption of water molecules onto the silica surface, forming silanol (SiOH) groups at a concentration consistent with experimental data. Water penetration and silanol formation ∼7-8 Å below the outermost oxygen are observed. Because of the dissociative nature of the water potential, formation of hydronium ions is allowed, and, whereas seldom observed in the simulations of bulk water, hydronium ions are formed during the reactions causing the formation of the silanols. The formation of hydronium ions has also been observed in ab initio calculations of water adsorption onto silica surfaces. The time evolution of the reactions involving hydronium ions in our MD simulations is similar to that observed in first-principles MD calculations. Hydronium ions offer a mechanism by which initially singly coordinated terminal oxygen (Si-O -) receives a H + ion from a relatively distant chemisorbed H 2 O molecule via multiple H + ion transfer, creating two SiOH sites.

Understanding the dissolution of silicate glasses and minerals from atomic to macroscopic levels is a challenge with major implications in geoscience and industry. One of the main uncertainties limiting the development of predictive models lies in the formation of an amorphous surface layer––called gel––that can in some circumstances control the reactivity of the buried interface. Here, we report experimental and simulation results deciphering the mechanisms by which the gel becomes passivating. The study conducted on a six-oxide borosilicate glass shows that gel reorganization involving high exchange rate of oxygen and low exchange rate of silicon is the key mechanism accounting for extremely low apparent water diffusivity (∼10−21 m2 s−1), which could be rate-limiting for the overall reaction. These findings could be used to improve kinetic models, and inspire the development of new molecular sieve materials with tailored properties as well as highly durable glass for application in extreme environments.

Table 1c.

Dilatometric measurement of the thermal expansion of water in porous silica shows that the expansion coefficient increases systematically as the pore size decreases below about 15 nm. This behavior is quantitatively reproduced by molecular dynamics (MD) simulations based on a new dissociative potential. According to MD, the structure of the water is modified within approximately 6 A of the pore wall, so that it resembles bulk water at a higher pressure. On the basis of this observation, it is possible to account for the measured expansion, as the thermal expansion coefficient of bulk water increases with temperature over the range considered in this study.