Multivalent Surface Cations Enhance Heterogeneous Freezing of Water on Muscovite Mica

Insights into the structure and dynamics of K+ ions at the muscovite-water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. In this work, we make a first step in this direction by investigating the behavior of the native K+ ions at the interface between neat water and the muscovite mica (001) surface with ab-initio-based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K+ ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K+ diffusion between mica's ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K+ ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Interaction of surface cations of cleaved mica with water in vapor and liquid forms

Natural minerals contain ions that become hydrated when they come into contact with water in vapor and liquid forms. Muscovite mica - a common phyllosilicate with perfect cleavage planes - is an ideal system to investigate the details of ion hydration. The cleaved mica surface is decorated by an array of K+ ions that can be easily exchanged with other ions or protons when immersed in an aqueous solution. Despite the vast interest in the atomic-scale hydration processes of these K+ ions, experimental data under controlled conditions have remained elusive. Here, atomically resolved non-contact atomic force microscopy (nc-AFM) is combined with X-ray photoelectron spectroscopy (XPS) to investigate the cation hydration upon dosing water vapor at 100 K in ultra-high vacuum (UHV). The cleaved surface is further exposed to ultra-clean liquid water at room temperature, which promotes ion mobility and partial ion-to-proton substitution. The results offer the first direct experimental views of the interaction of water with muscovite mica under UHV. The findings are in line with previous theoretical predictions.

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation

Heterogeneous nucleation is the dominant form of liquid-to-solid transition in nature. Although molecular simulations are most uniquely suited to studying nucleation, the waiting time to observe even a single nucleation event can easily exceed the current computational capabilities. Therefore, there exists an imminent need for methods that enable computationally fast and feasible studies of heterogeneous nucleation. Seeding is a technique that has proven to be successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of homogeneous crystal nucleation. In this article, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudorigid bodies and simulated on a surface with metastable liquid above its melting temperature. This allows the seeds to adapt to the surface and identify favorable seed-surface configurations, which is necessary for reliable predictions of crystal polymorphs that form and the corresponding heterogeneous nucleation rates. We demonstrate and validate RSeeds for heterogeneous ice nucleation on a flexible self-assembled monolayer surface, a mineral surface based on kaolinite, and two model surfaces. RSeeds predicts the correct ice polymorph, exposed crystal plane, and rotation on the surface. RSeeds is semiquantitative and can be used to estimate the critical nucleus size and nucleation rate when combined with classical nucleation theory. We demonstrate that RSeeds can be used to evaluate nucleation rates spanning many orders of magnitude.

Resolving the intrinsic short-range ordering of K+ ions on cleaved muscovite mica

Muscovite mica, KAl₂(Si₃Al)O₁₀(OH)₂, is a common layered phyllosilicate with perfect cleavage planes. The atomically flat surfaces obtained through cleaving lend themselves to scanning probe techniques with atomic resolution and are ideal to model minerals and clays. Despite the importance of the cleaved mica surfaces, several questions remain unresolved. It is established that K⁺ ions decorate the cleaved surface, but their intrinsic ordering - unaffected by the interaction with the environment - is not known. This work presents clear images of the K⁺ distribution of cleaved mica obtained with low-temperature non-contact atomic force microscopy (AFM) under ultra-high vacuum (UHV) conditions. The data unveil the presence of short-range ordering, contrasting previous assumptions of random or fully ordered distributions. Density functional theory (DFT) calculations and Monte Carlo simulations show that the substitutional subsurface Al³⁺ ions have an important role for the surface K⁺ ion arrangement.

Ion-dependent protein-surface interactions from intrinsic solvent response

The phyllosilicate mineral muscovite mica is widely used as a surface template for the patterning of macromolecules, yet a molecular understanding of its surface chemistry under varying solution conditions, required to predict and control the self-assembly of adsorbed species, is lacking. We utilize all-atom molecular dynamics simulations in conjunction with an electrostatic analysis based in local molecular field theory that affords a clean separation of long-range and short-range electrostatics. Using water polarization response as a measure of the electric fields that arise from patterned, surface-bound ions that direct the adsorption of charged macromolecules, we apply a Landau theory of forces induced by asymmetrically polarized surfaces to compute protein-surface interactions for two muscovite-binding proteins (DHR10-mica6 and ^C98RhuA). Comparison of the pressure between surface and protein in high-concentration KCl and NaCl aqueous solutions reveals ion-specific differences in far-field protein-surface interactions, neatly capturing the ability of ions to modulate the surface charge of muscovite that in turn selectively attracts one binding face of each protein over all others.

Examining the Role of Different Molecular Interactions on Activation Energies and Activation Volumes in Liquid Water

There are a large number of force fields available to model water in molecular dynamics simulations, which each have their own strengths and weaknesses in describing the behavior of the liquid. One particular weakness in many of these models is their description of dynamics away from ambient conditions, where their ability to reproduce measurements is mixed. To investigate this issue, we use the recently developed fluctuation theory for dynamics to directly evaluate measures of the local temperature and pressure dependence: the activation energy and the activation volume. We examine these activation parameters for hydrogen-bond jump exchange times, OH reorientation times, and diffusion coefficients calculated from the SPC/E, SPC/Fw, TIP3P-PME, TIP3P-PME/Fw, OPC3, TIP4P/2005, TIP4P/Ew, E3B2, and E3B3 water models. Activation energy decompositions available through the fluctuation theory approach provide mechanistic insight into the origins of different temperature dependences between the various models, as well as the influence of three-body effects and flexibility.

Brief Overview of Ice Nucleation

The nucleation of ice is vital in cloud physics and impacts on a broad range of matters from the cryopreservation of food, tissues, organs, and stem cells to the prevention of icing on aircraft wings, bridge cables, wind turbines, and other structures. Ice nucleation thus has broad implications in medicine, food engineering, mineralogy, biology, and other fields. Nowadays, the growing threat of global warming has led to intense research activities on the feasibility of artificially modifying clouds to shift the Earth’s radiation balance. For these reasons, nucleation of ice has been extensively studied over many decades and rightfully so. It is thus not quite possible to cover the whole subject of ice nucleation in a single review. Rather, this feature article provides a brief overview of ice nucleation that focuses on several major outstanding fundamental issues. The author’s wish is to aid early researchers in ice nucleation and those who wish to get into the field of ice nucleation from other disciplines by concisely summarizing the outstanding issues in this important field. Two unresolved challenges stood out from the review, namely the lack of a molecular-level picture of ice nucleation at an interface and the limitations of classical nucleation theory.