Aqueous systems are omnipresent in nature and technology. They show complex behaviors, which often originate in the existence of hydrogen-bond networks. Prominent examples are the anomalies of water and the non-ideal behaviors of aqueous solutions. The phenomenology becomes even richer when aqueous liquids are subject to confinement. To this day, many properties of water and its mixtures, in particular, under confinement, are not understood. In recent years, molecular dynamics simulations developed into a powerful tool to improve our knowledge in this field. Here, our simulation results for water and aqueous mixtures in the bulk and in various confinements are reviewed and some new simulation data are added to improve our knowledge about the role of interfaces. Moreover, findings for water are compared with results for silica, exploiting that both systems form tetrahedral networks.

We investigate anomalies in liquid silica with molecular dynamics simulations and present evidence for a fragile-to-strong transition at around 3100 K-3300 K. To this purpose, we studied the structure and dynamical properties of silica over a wide temperature range, finding four indicators of a fragile-to-strong transition. First, there is a density minimum at around 3000 K and a density maximum at 4700 K. The turning point is at 3400 K. Second, the local structure characterized by the tetrahedral order parameter changes dramatically around 3000 K from a higher-ordered, lower-density phase to a less ordered, higher-density phase. Third, the correlation time τ changes from an Arrhenius behavior below 3300 K to a Vogel-Fulcher-Tammann behavior at higher temperatures. Fourth, the Stokes-Einstein relation holds for temperatures below 3000 K, but is replaced by a fractional relation above this temperature. Furthermore, our data indicate that dynamics become again simple above 5000 K, with Arrhenius behavior and a classical Stokes-Einstein relation.

No abstract

We analyze the effects of spatial confinement on viscous silica using molecular dynamics simulations. For this purpose, we prepare a silica melt in a cylindrical pore, which is produced by pinning appropriate fractions of silicon and oxygen atoms in a bulk system after an equilibration period. In this way, the structure of the confined silica melt remains unaffected, while the confinement has a strong impact on the dynamics. We find that the structural relaxation of viscous silica is slowed down according to a double exponential law when approaching the pore wall. Moreover, we observe that static density correlations exist in the vicinity of the pore wall. Based on these effects, we determine dynamical and structural length scales of the silica melt. Both length scales show a similar increase upon cooling, with values on the order of the next-neighbor distances in the studied temperature range. Interestingly, we find no evidence that the growth of the length scales is affected by a fragile-to-strong transition of the silica melt. This observation casts serious doubts on the relevance of these length scales for the structural relaxation, at least for the studied glass former.

We study the wetting properties of water on silica surfaces using molecular dynamics (MD) simulations. To describe the intermolecular interaction between water and silica atoms, two types of interaction potential models are used: the standard Bródka and Zerda (BZ) model, and the Gulmen and Thompson (GT) model. We perform an in-depth analysis of the influence of the choice of the potential on the arrangement of the water molecules in partially filled pores and on top of silica slabs. We find that at moderate pore filling ratios, the GT silica surface is completely wetted by water molecules, which agrees well with experimental findings, while the commonly used BZ surface is less hydrophilic and is only partially wetted. We interpret our simulation results using an analytical calculation of the phase diagram of water in partially filled pores. Moreover, an evaluation of the contact angle of the water droplet on top of the silica slab reveals that the interaction becomes more hydrophilic with increasing slab thickness and saturates around 2.5 − 3 nm, in agreement with the experimentally found value. Our analysis also shows that the hydroaffinity of the surface is mainly determined by the electrostatic interaction, but that the van der Waals interaction nevertheless is strong enough that it can turn a hydrophobic surface into a hydrophilic surface.