Metadynamics simulations of calcite crystallization on self-assembled monolayers

We show that recent developments in the application of metadynamics methods to direct simulations of crystallization make it possible to predict the orientation of crystals grown on self-assembled monolayers. In contrast to previous studies, the method allows for dynamic treatment of the organic component and the inclusion of explicit surface water without the need for computationally intensive interfacial energy calculations or prior knowledge of the interfacial structure. The method is applied to calcite crystallization on carboxylate terminated alkanethiols arrayed on Au (111). We demonstrate that a dynamic treatment of the monolayer is sufficient to reproduce the experimental results without the need to impose epitaxial constraints on the system. We also observe an odd-even effect in the variation of selectivity with organic chain length, reproducing experimentally observed orientations in both cases. Analysis of the ordering process in our simulations suggests a cycle of mutual control in which both the organic and mineral components induce complementary local order across the interface, leading to the formation of a critical crystalline region. The influence of pH, together with some factors that might affect the range of applicability of our method, is discussed.