PNNL

Abstract

New insights into the structure of the calcite-water interface are obtained through direct model-independent comparison of multiple classical molecular dynamics (MD) simulations with high-resolution specular X-ray reflectivity (XR) data. This set of comparisons, with four different state-of-the-art force fields (including two non-polarizable, one polarizable, and one reactive force field), reveal new insights into the absolute accuracy of the simulated structures and the uniqueness of the XR-derived structural results. These four simulations, while qualitatively similar, have visibly distinct interfacial structure, and are distinguished through a quantitative comparison of the XR signals calculated from these simulations with experimental XR data. The results demonstrate that the simulated calcite-water interface structures, as a whole, are not consistent with the XR data (i.e., within their precision and accuracy). This disagreement is largely due to the simulation of the calcite lattice. The simulated interfacial water profiles show substantially different levels of agreement with the XR data. Of these, the rigid-ion model (RIM) simulations show the best consistency with the experimental XR data. Further model-dependent comparisons of the structural parameters that describe the interfacial structure (derived from both the MD simulations and the XR data) provide further insight into the sources of differences between these two approaches. Using the new insights from the RIM simulations, new structures of the calcite-water interface consistent with both the experimental data and the simulation are identified and compared to recent results.