Hydration of mineral surfaces, a critical process for many technological applications, encompasses multiple coupled chemical reactions and topological changes, challenging both experimental characterization and computational modeling. In this work, we used reactive force field simulations to understand the surface properties, hydration, and dissolution of a model mineral, tricalcium silicate. We show that the computed static quantities, i.e., surface energies and water adsorption energies, do not provide useful insight into predict mineral hydration because they do not account for major structural changes at the interface when dynamic effects are included. Upon hydration, hydrogen atoms from dissociated water molecules penetrate into the crystal, forming a disordered calcium silicate hydrate layer that is similar for most of the surfaces despite wide-ranging static properties. Furthermore, the dynamic picture of hydration reveals the hidden role of surface topology, which can lead to unexpected water tessellation that stabilizes the surface against dissolution. KEYWORDS: dissolution, hydration, molecular dynamics, surface properties, water adsorption, calcium silicate
■ INTRODUCTIONThe hydration and dissolution of minerals is of crucial importance for a broad range of natural and synthetic phenomena: from the stability of minerals and rocks in watersheds and aquifers, 1 to hydrogen production and other heterogeneous catalytic reactions, 2,3 to the production, utilization, and ultimate degradation of building materials, 4 and to the durability of biomaterials, 5 to name only a few examples. The key properties at play during dissolution are controlled by interfacial reactions between water molecules and the solid surfaces, which can be probed experimentally via techniques such as AFM, vertical scanning interferometry (VSI), and the study of isotope exchange or solute fluxes to provide general information regarding the reaction rates or qualitative information about the formation of steps, pitches, and so on. 6 However, in many applications, much more control over the hydration and dissolution is desired, such as the ability to moderate the rate 7 or preferentially favor the dissolution of one crystal facet over other.8 To achieve such control, a deeper understanding of the molecular mechanisms governing hydration and dissolution is necessary.In this work, we use atomic-scale simulations to study the hydration of a calcium orthosilicate. We take a multistep approach, which is critical to piece together the relevant physical processes underlying hydration and dissolution. First, we analyze the surface energies and water molecule adsorption, extending beyond traditional approaches in order to obtain a more accurate picture of these properties. Next, we compare these surface properties with the results obtained from simulating a realistic hydration process over a period of 2 ns. We discuss the differences between the static and dynamic pictures and show that topological changes render the properties of the f...