Can Ice-Like Structures Form on Non-Ice-Like Substrates? The Example of the K-feldspar Microcline

Feldspar minerals are the most common rock formers in Earth’s crust. As such they play an important role in subjects ranging from geology to climate science. An atomistic understanding of the feldspar structure and its interaction with water is therefore desirable, not least because feldspar has been shown to dominate ice nucleation by mineral dusts in Earth’s atmosphere. The complexity of the ice/feldspar interface arising from the numerous chemical motifs expressed on the surface makes it a challenging system. Here we report a comprehensive study of this challenging system with ab initio density functional theory calculations. We show that the distribution of Al atoms, which is crucial for the dissolution kinetics of tectosilicate minerals, differs significantly between the bulk environment and on the surface. Furthermore, we demonstrate that water does not form ice-like overlayers in the contact layer on the most easily cleaved (001) surface of K-feldspar. We do, however, identify contact layer structures of water that induce ice-like ordering in the second overlayer. This suggests that even substrates without an apparent match with the ice structure may still act as excellent ice nucleating agents

Ice formation on kaolinite: Insights from molecular dynamics simulations

The formation of ice affects many aspects of our everyday life as well as important technologies such as cryotherapy and cryopreservation. Foreign substances almost always aid water freezing through heterogeneous ice nucleation, but the molecular details of this process remain largely unknown. In fact, insight into the microscopic mechanism of ice formation on different substrates is difficult to obtain even if state-of-the-art experimental techniques are used. At the same time, atomistic simulations of heterogeneous ice nucleation frequently face extraordinary challenges due to the complexity of the water-substrate interaction and the long time scales that characterize nucleation events. Here, we have investigated several aspects of molecular dynamics simulations of heterogeneous ice nucleation considering as a prototypical ice nucleating material the clay mineral kaolinite, which is of relevance in atmospheric science. We show via seeded molecular dynamics simulations that ice nucleation on the hydroxylated (001) face of kaolinite proceeds exclusively via the formation of the hexagonal ice polytype. The critical nucleus size is two times smaller than that obtained for homogeneous nucleation at the same supercooling. Previous findings suggested that the flexibility of the kaolinite surface can alter the time scale for ice nucleation within molecular dynamics simulations. However, we here demonstrate that equally flexible (or non flexible) kaolinite surfaces can lead to very different outcomes in terms of ice formation, according to whether or not the surface relaxation of the clay is taken into account. We show that very small structural changes upon relaxation dramatically alter the ability of kaolinite to provide a template for the formation of a hexagonal overlayer of water molecules at the water-kaolinite interface, and that this relaxation therefore determines the nucleation ability of this mineral

Active sites in heterogeneous ice nucleation-the example of K-rich feldspars

Ice formation on aerosol particles is a process of crucial importance to Earth's climate and the environmental sciences but it is not understood at the molecular level. This is so partly because the nature of active sites, local surface features where ice growth commences, is still unclear. Here we report direct electron-microscopic observations of deposition growth of aligned ice crystals on feldspar, an atmospherically important component of mineral dust. Our molecular-scale computer simulations indicate that this alignment arises from the preferential nucleation of prismatic crystal planes of ice on high-energy (100) surface planes of feldspar. The microscopic patches of (100) surface, exposed at surface defects such as steps, cracks, and cavities, are thought to be responsible for the high ice nucleation efficacy of K-feldspar particles

Atomic-Scale Picture of the Composition, Decay, and Oxidation of Two-Dimensional Radioactive Films

Two-dimensional radioactive (125)I monolayers are a recent development that combines the fields of radiochemistry and nanoscience. These Au-supported monolayers show great promise for understanding the local interaction of radiation with 2D molecular layers, offer different directions for surface patterning, and enhance the emission of chemically and biologically relevant low-energy electrons. However, the elemental composition of these monolayers is in constant flux due to the nuclear transmutation of (125)I to (125)Te, and their precise composition and stability under ambient conditions has yet to be elucidated. Unlike I, which is stable and unreactive when bound to Au, the newly formed Te atoms would be expected to be more reactive. We have used electron emission and X-ray photoelectron spectroscopy (XPS) to quantify the emitted electron energies and to track the film composition in vacuum and the effect of exposure to ambient conditions. Our results reveal that the Auger electrons emitted during the ultrafast radioactive decay process have a kinetic energy corresponding to neutral Te. By combining XPS and scanning tunneling microscopy experiments with density functional theory, we are able to identify the reaction of newly formed Te to TeO2 and its subsequent dimerization. The fact that the Te2O4 units stay intact during major lateral rearrangement of the monolayer illustrates their stability. These results provide an atomic-scale picture of the composition and mobility of surface species in a radioactive monolayer as well as an understanding of the stability of the films under ambient conditions, which is a critical aspect in their future applications

Enhancement of low-energy electron emission in 2D radioactive films

High-energy radiation has been used for decades; however, the role of low-energy electrons created during irradiation has only recently begun to be appreciated. Low-energy electrons are the most important component of radiation damage in biological environments because they have subcellular ranges, interact destructively with chemical bonds, and are the most abundant product of ionizing particles in tissue. However, methods for generating them locally without external stimulation do not exist. Here, we synthesize one-atom-thick films of the radioactive isotope (125)I on gold that are stable under ambient conditions. Scanning tunnelling microscopy, supported by electronic structure simulations, allows us to directly observe nuclear transmutation of individual (125)I atoms into (125)Te, and explain the surprising stability of the 2D film as it underwent radioactive decay. The metal interface geometry induces a 600% amplification of low-energy electron emission (<10 eV; ref. ) compared with atomic (125)I. This enhancement of biologically active low-energy electrons might offer a new direction for highly targeted nanoparticle therapies

From gold to dust: a computational journey

The behaviour of atoms and molecules on different substrates is of paramount importance to an immense variety of reactions and processes in our everyday life. Catalytic reactions, the formation of ice in clouds and processes happening in fuel cells, solar cells and nanomaterials are just few such examples. The ubiquity of these systems across numerous scientific disciplines is caused, to some extent, by the sheer number of interfaces that are possible and by the complex interactions that can arise at the interface between substrate and adsorbate. Not only can the interactions at these interfaces often be tuned to optimise a particular reaction by changing substrate, but they can also be studied at the atomic level experimentally and computationally and hence can provide fascinating new insights into nature. In this thesis, a variety of interfaces, ranging from dust particles to gold surfaces, has been studied with computational techniques. Some chapters were directly motivated by experimental measurements from collaborators and others were aimed at furthering our fundamental knowledge and obtaining general insight. Exploring and understanding the structures that a substrate can induce in its adsorbates were often crucial steps in achieving this. Unravelling the interface between feldspar, a common dust particle in Earth's atmosphere, and ice, for example allowed the identification of an active site for ice nucleation at the molecular level for the first time. Similarly, discovering and characterizing the behaviour of radioactive two-dimensional films on gold surfaces opened the door to a potentially new research field, radiochemistry at the nano scale

Anomalously Low Barrier for Water Dimer Diffusion on Cu(111)

A molecular-scale description of water and ice is important in fields as diverse as atmospheric chemistry, astrophysics, and biology. Despite a detailed understanding of water and ice structures on a multitude of surfaces, relatively little is known about the kinetics of water motion on surfaces. Here, we report a detailed study on the diffusion of water monomers and the formation and diffusion of water dimers through a combination of time-lapse low-temperature scanning tunnelling microscopy experiments and first-principles electronic structure calculations on the atomically flat Cu(111) surface. On the basis of an unprecedented long-time study of individual water monomers and dimers over days, we establish rates and mechanisms of water monomer and dimer diffusion. Interestingly, we find that the monomer and the dimer diffusion barriers are similar, despite the significantly larger adsorption energy of the dimer. This is thus a violation of the rule of thumb that relates diffusion barriers to adsorption energies, an effect that arises because of the directional and flexible hydrogen bond within the dimer. This flexibility during diffusion should also be relevant for larger water clusters and other hydrogen-bonded adsorbates. Our study stresses that a molecular-scale understanding of the initial stages of ice nanocluster formation is not possible on the basis of static structure investigations alone

Ice formation on kaolinite: Insights from molecular dynamics simulations

The formation of ice affects many aspects of our everyday life as well as important technologies such as cryotherapy and cryopreservation. Foreign substances almost always aid water freezing through heterogeneous ice nucleation, but the molecular details of this process remain largely unknown. In fact, insight into the microscopic mechanism of ice formation on different substrates is difficult to obtain even if state-of-the-art experimental techniques are used. At the same time, atomistic simulations of heterogeneous ice nucleation frequently face extraordinary challenges due to the complexity of the water-substrate interaction and the long time scales that characterize nucleation events. Here, we have investigated several aspects of molecular dynamics simulations of heterogeneous ice nucleation considering as a prototypical ice nucleating material the clay mineral kaolinite, which is of relevance in atmospheric science. We show via seeded molecular dynamics simulations that ice nucleation on the hydroxylated (001) face of kaolinite proceeds exclusively via the formation of the hexagonal ice polytype. The critical nucleus size is two times smaller than that obtained for homogeneous nucleation at the same supercooling. Previous findings suggested that the flexibility of the kaolinite surface can alter the time scale for ice nucleation within molecular dynamics simulations. However, we here demonstrate that equally flexible (or non flexible) kaolinite surfaces can lead to very different outcomes in terms of ice formation, according to whether or not the surface relaxation of the clay is taken into account. We show that very small structural changes upon relaxation dramatically alter the ability of kaolinite to provide a template for the formation of a hexagonal overlayer of water molecules at the water-kaolinite interface, and that this relaxation therefore determines the nucleation ability of this mineral

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments