Water-Ice Analogues of Polycyclic Aromatic Hydrocarbons: Water Nanoclusters on Cu(111)

Quasi-one-dimensional hydrogen bonding in nanoconfined ice

The Bernal-Fowler ice rules stipulate that each water molecule in an ice crystal should form four hydrogen bonds. However, in extreme or constrained conditions, the arrangement of water molecules deviates from conventional ice rules, resulting in properties significantly different from bulk water. In this study, we employ machine learning-driven first-principles simulations to identify a new stabilization mechanism in nanoconfined ice phases. Instead of forming four hydrogen bonds, nanoconfined crystalline ice can form a quasi-one-dimensional hydrogen-bonded structure that exhibits only two hydrogen bonds per water molecule. These structures consist of strongly hydrogen-bonded linear chains of water molecules that zig-zag along one dimension, stabilized by van der Waals interactions that stack these chains along the other dimension. The unusual interplay of hydrogen bonding and van der Waals interactions in nanoconfined ice results in atypical proton behavior such as potential ferroelectric behavior, low dielectric response, and long-range proton dynamics.

The structure-giving role of Rb+ ions for water-ice nanoislands supported on Cu(111)

We characterize the effect of rubidium ions on water-ice nanoislands in terms of area, fractal dimension, and apparent height by low-temperature scanning tunneling microscopy. Water nanoislands on the pristine Cu(111) surface are compared to those at similar coverage on a Rb+ pre-covered Cu(111) surface to reveal the structure-giving effect of Rb+. The presence of Rb+ induces changes in the island shape, and hence, the water network, without affecting the nanoisland volume. The broad area distribution shifts to larger values while the height decreases from three bilayers to one or two bilayers. The nanoislands on the Rb+ pre-covered surface are also more compact, reflected in a shift in the fractal dimension distribution. We relate the changes to a weakening of the hydrogen-bond network by Rb+.

Structure Discovery in Atomic Force Microscopy Imaging of Ice

The interaction of water with surfaces is crucially important in a wide range of natural and technological settings. In particular, at low temperatures, unveiling the atomistic structure of adsorbed water clusters would provide valuable data for understanding the ice nucleation process. Using high-resolution atomic force microscopy (AFM) and scanning tunneling microscopy, several studies have demonstrated the presence of water pentamers, hexamers, and heptamers (and of their combinations) on a variety of metallic surfaces, as well as the initial stages of 2D ice growth on an insulating surface. However, in all of these cases, the observed structures were completely flat, providing a relatively straightforward path to interpretation. Here, we present high-resolution AFM measurements of several water clusters on Au(111) and Cu(111), whose understanding presents significant challenges due to both their highly 3D configuration and their large size. For each of them, we use a combination of machine learning, atomistic modeling with neural network potentials, and statistical sampling to propose an underlying atomic structure, finally comparing its AFM simulated images to the experimental ones. These results provide insights into the early phases of ice formation, which is a ubiquitous phenomenon ranging from biology to astrophysics.

Identification of a common ice nucleus on hydrophilic and hydrophobic close-packed metal surfaces

Establishing a general model of heterogeneous ice nucleation has long been challenging because of the surface water structures found on different substrates. Identifying common water clusters, regardless of the underlying substrate, is one of the key steps toward solving this problem. Here, we demonstrate the presence of a common water cluster found on both hydrophilic Pt(111) and hydrophobic Cu(111) surfaces using scanning tunneling microscopy and non-contact atomic force microscopy. Water molecules self-assemble into a structure with a central flat-lying hexagon and three fused pentagonal rings, forming a cluster consisting of 15 individual water molecules. This cluster serves as a critical nucleus during ice nucleation on both surfaces: ice growth beyond this cluster bifurcates to form two-dimensional (three-dimensional) layers on hydrophilic (hydrophobic) surfaces. Our results reveal the inherent similarity and distinction at the initial stage of ice growth on hydrophilic and hydrophobic close-packed metal surfaces; thus, these observations provide initial evidence toward a general model for water-substrate interaction.

Hydrogen Exchange through Hydrogen Bonding between Methanol and Water in the Adsorbed State on Cu(111)

The interaction between submonolayers of methanol and water on Cu(111) is studied at 95-160 K temperature range with surface-sensitive infrared spectroscopy using isotopically labeled molecules. The initial interaction of methanol with the preadsorbed amorphous solid water at 95 K is through hydrogen-bonding with the dangling hydroxyl groups of water. Upon increasing the temperature up to 140 K, methanol and deuterated water form H-bonded structures which allow hydrogen-deuterium exchange between the hydroxyl group of methanol and the deuterated water. The evolution of the O-D and O-H stretching bands indicate that the hydrogen transfer is dominant at around 120-130 K, slightly below the desorption temperature of methanol. Above 140 K, methanol desorbs and a mixture of hydrogen-related water isotopologues remains on the surface. The isotopic composition of this mixture versus the initial D₂O:CH₃OH ratio supports a potential exchange mechanism via hydrogen hopping between alternating methanol and water molecules in a hydrogen-bonded network.

Electrostatic Discovery Atomic Force Microscopy

While offering high resolution atomic and electronic structure, scanning probe microscopy techniques have found greater challenges in providing reliable electrostatic characterization on the same scale. In this work, we offer electrostatic discovery atomic force microscopy, a machine learning based method which provides immediate maps of the electrostatic potential directly from atomic force microscopy images with functionalized tips. We apply this to characterize the electrostatic properties of a variety of molecular systems and compare directly to reference simulations, demonstrating good agreement. This approach offers reliable atomic scale electrostatic maps on any system with minimal computational overhead.

Insights into lithium manganese oxide-water interfaces using machine learning potentials

Unraveling the atomistic and the electronic structure of solid-liquid interfaces is the key to the design of new materials for many important applications, from heterogeneous catalysis to battery technology. Density functional theory (DFT) calculations can, in principle, provide a reliable description of such interfaces, but the high computational costs severely restrict the accessible time and length scales. Here, we report machine learning-driven simulations of various interfaces between water and lithium manganese oxide (Li_xMn₂O₄), an important electrode material in lithium ion batteries and a catalyst for the oxygen evolution reaction. We employ a high-dimensional neural network potential to compute the energies and forces several orders of magnitude faster than DFT without loss in accuracy. In addition, a high-dimensional neural network for spin prediction is utilized to analyze the electronic structure of the manganese ions. Combining these methods, a series of interfaces is investigated by large-scale molecular dynamics. The simulations allow us to gain insights into a variety of properties, such as the dissociation of water molecules, proton transfer processes, and hydrogen bonds, as well as the geometric and electronic structure of the solid surfaces, including the manganese oxidation state distribution, Jahn-Teller distortions, and electron hopping.

Origins of fast diffusion of water dimers on surfaces

The diffusion of water molecules and clusters across the surfaces of materials is important to a wide range of processes. Interestingly, experiments have shown that on certain substrates, water dimers can diffuse more rapidly than water monomers. Whilst explanations for anomalously fast diffusion have been presented for specific systems, the general underlying physical principles are not yet established. We investigate this through a systematic ab initio study of water monomer and dimer diffusion on a range of surfaces. Calculations reveal different mechanisms for fast water dimer diffusion, which is found to be more widespread than previously anticipated. The key factors affecting diffusion are the balance of water-water versus water-surface bonding and the ease with which hydrogen-bond exchange can occur (either through a classical over-the-barrier process or through quantum-mechanical tunnelling). We anticipate that the insights gained will be useful for understanding future experiments on the diffusion and clustering of hydrogen-bonded adsorbates.

The experimental observation that water dimers diffuse more rapidly than monomers across materials’ surfaces is yet to be clarified. Here the authors show by ab initio calculations classical and quantum mechanical mechanisms for faster water dimer diffusion on a broad range of metal and non-metal surfaces.

Rich topologies of monolayer ices via unconventional electrowetting

Accurate manipulation of a substance on the nanoscale and ultimately down to the level of a single atom or molecule is an ongoing subject of frontier research. Herein, we show that topologies of water monolayers on substrates, in the complete wetting condition, can be manipulated into rich forms of ordered structures via electrowetting. Notably, two new topologies of monolayer ices were identified from our molecular dynamics simulations: one stable below room temperature and the other one having the ability to be stable at room temperature. Moreover, the wettability of the substrate can be tuned from superhydrophobic to superhydrophilic by uniformly changing the charge of each atomic site of the dipole or quadrupole distributed in an orderly manner on the model substrate. At a certain threshold value of the atomic charge, water droplets on the substrate can spread out spontaneously, achieving a complete electrowetting. Importantly, unlike the conventional electrowetting, which involves application of a uniform external electric field, we proposed non-conventional electrowetting, for the first time, by invoking the electric field of dipoles and quadrupoles embedded in the substrate. Moreover, different topologies of water monolayers can be achieved by using the non-conventional electrowetting. A major advantage of the non-conventional electrowetting is that the contact-angle saturation, a long-standing and known limitation in the field of electrowetting, can be overcome by tuning uniformly the lattice atomic charge at the surface, thereby offering a new way to mitigate the contact-angle saturation for various electrowetting applications.

Strain Relief during Ice Growth on a Hexagonal Template

Heterogeneous ice nucleation at solid surfaces impacts many areas of science, from environmental processes, such as precipitation, to microbiological systems and food processing, but the microscopic mechanisms underpinning nucleation remain unclear. Discussion of ice growth has often focused around the role of the surface in templating the structure of water, forcing the first layer to adopt the registry of the underlying substrate rather than that of ice. To grow a thick ice film, water in the first few ice layers must accommodate this strain, but understanding how this occurs requires detailed molecular-scale information that is lacking. Here we combine scanning tunneling microscopy, low-energy electron diffraction, and work-function measurements with electronic structure calculations to investigate the initial stages of ice growth on a Pt alloy surface, having a lattice spacing 6% larger than ice. Although the first layer of water forms a strictly commensurate hexagonal network, this behavior does not extend to the second layer. Instead, water forms a 2D structure containing extended defect rows made from face-sharing pentamer and octamer rings. The defect rows allow the majority of second-layer water to remain commensurate with the solid surface while compensating lateral strain by increasing the water density close to that of an ice surface. The observation of octamer-pentamer rows in ice films formed on several surfaces suggests that the octamer-pentamer defect motif acts as a flexible strain relief mechanism in thin ice films, providing a mechanism that is not available during the growth of strained films in other materials, such as semiconductors.

Identifying Few-Molecule Water Clusters with High Precision on Au(111) Surface

Revealing the nature of a hydrogen-bond network in water structures is one of the imperative objectives of science. With the use of a low-temperature scanning tunneling microscope, water clusters on a Au(111) surface were directly imaged with molecular resolution by a functionalized tip. The internal structures of the water clusters as well as the geometry variations with the increase of size were identified. In contrast to a buckled water hexamer predicted by previous theoretical calculations, our results present deterministic evidence for a flat configuration of water hexamers on Au(111), corroborated by density functional theory calculations with properly implemented van der Waals corrections. The consistency between the experimental observations and improved theoretical calculations not only renders the internal structures of absorbed water clusters unambiguously, but also directly manifests the crucial role of van der Waals interactions in constructing water-solid interfaces.

Chirality at two-dimensional surfaces: A perspective from small molecule alcohol assembly on Au(111)

The delicate balance between hydrogen bonding and van der Waals interactions determines the stability, structure, and chirality of many molecular and supramolecular aggregates weakly adsorbed on solid surfaces. Yet the inherent complexity of these systems makes their experimental study at the molecular level very challenging. In this quest, small alcohols adsorbed on metal surfaces have become a useful model system to gain fundamental insight into the interplay of such molecule-surface and molecule-molecule interactions. Here, through a combination of scanning tunneling microscopy and density functional theory, we compare and contrast the adsorption and self-assembly of a range of small alcohols from methanol to butanol on Au(111). We find that longer chained alcohols prefer to form zigzag chains held together by extended hydrogen bonded networks between adjacent molecules. When alcohols bind to a metal surface datively via one of the two lone electron pairs of the oxygen atom, they become chiral. Therefore, the chain structures are formed by a hydrogen-bonded network between adjacent molecules with alternating adsorbed chirality. These chain structures accommodate longer alkyl tails through larger unit cells, while the position of the hydroxyl group within the alcohol molecule can produce denser unit cells that maximize intermolecular interactions. Interestingly, when intrinsic chirality is introduced into the molecule as in the case of 2-butanol, the assembly changes completely and square packing structures with chiral pockets are observed. This is rationalized by the fact that the intrinsic chirality of the molecule directs the chirality of the adsorbed hydroxyl group meaning that heterochiral chain structures cannot form. Overall this study provides a general framework for understanding the effect of simple alcohol molecular adstructures on hydrogen bonded aggregates and paves the way for rationalizing 2D chiral supramolecular assembly.

Water agglomerates on Fe3O4(001)

Determining the structure of water adsorbed on solid surfaces is a notoriously difficult task and pushes the limits of experimental and theoretical techniques. Here, we follow the evolution of water agglomerates on Fe3O4(001); a complex mineral surface relevant in both modern technology and the natural environment. Strong OH-H2O bonds drive the formation of partially dissociated water dimers at low coverage, but a surface reconstruction restricts the density of such species to one per unit cell. The dimers act as an anchor for further water molecules as the coverage increases, leading first to partially dissociated water trimers, and then to a ring-like, hydrogen-bonded network that covers the entire surface. Unraveling this complexity requires the concerted application of several state-of-the-art methods. Quantitative temperature-programmed desorption (TPD) reveals the coverage of stable structures, monochromatic X-ray photoelectron spectroscopy (XPS) shows the extent of partial dissociation, and noncontact atomic force microscopy (AFM) using a CO-functionalized tip provides a direct view of the agglomerate structure. Together, these data provide a stringent test of the minimum-energy configurations determined via a van der Waals density functional theory (DFT)-based genetic search.

The effects of different heterocycles and solvents on the ESIPT mechanisms of three novel photoactive mono-formylated benzoxazole derivatives

The detailed effects of different heterocycles and solvents on the dynamical ESIPT mechanisms of three novel mono-formylated benzoxazole derivatives A–C in two different surroundings have been expounded by the TDDFT method at the B3LYP/TZVP level.