Strain Relief during Ice Growth on a Hexagonal Template

Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces

Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.

Water-Induced Bandgap Engineering in Nanoribbons of Hexagonal Boron Nitride

Different from hexagonal boron nitride (hBN) sheets, the bandgap of hBN nanoribbons (BNNRs) can be changed by spatial/electrostatic confinement. It is predicted that a transverse electric field can narrow the bandgap and even cause an insulator-metal transition in BNNRs. However, experimentally introducing an overhigh electric field across the BNNR remains challenging. Here, it is theoretically and experimentally demonstrated that water adsorption greatly reduces the bandgap of zigzag-oriented BNNRs (zBNNRs). Ab initio calculations show that water molecules can be favorably assembled within the trench between two adjacent BNNRs to form a polar ice layer, which induces a transverse equivalent electric field of over 2 V nm-1 accounting for the bandgap reduction. Field-effect transistors are successfully fabricated from zBNNRs with different widths. The conductance of water-adsorbed zBNNRs can be tuned over 3 orders in magnitude via modulation of the equivalent electrical field at room temperature. Furthermore, photocurrent response measurements are taken to determine the optical bandgaps of zBNNRs with water adsorption. The zBNNR with increased width can exhibit a bandgap down to 1.17 eV. This study offers fundamental insights into new routes toward realizing electronic/optoelectronic devices and circuits based on hexagonal boron nitride.

Orientational Water Bonding on Pt(111): Beyond the Frontier Orbital Principle

For decades, our understanding of water-metal bonding has been dominated by the frontier orbital principle in which globally stable water-metal interactions are ruled by HOMO interacting with metal surfaces. Using density functional theory calculations, herein, we have revealed that the frontier orbital principle cannot be applied to metastable water bonding on Pt(111), where the decisive role of HOMO is replaced by HOMO-1 in terms of the greatest orbital shifts and depopulations as the two different bonding indicators. Unlike the stable water configuration in which both HOMO-1 and HOMO prefer to overlap with metal states through σ-like orbital interactions, metastable configurations exhibit delicate competition or balance between σ-like and π-like orbital interactions exerted by HOMO-1 and HOMO, respectively. These findings have significantly deepened our understanding of orbital roles in water-metal bonding interactions and bridged the gap between theoretical understanding of electrified waters at electrochemical interfaces and water science on metal surfaces.

Wetting of a Stepped Platinum (211) Surface

Steps stabilize water adsorption on metal surfaces, providing favorable binding sites for water during wetting or ice nucleation, but there is limited understanding of the local water arrangements formed on such surfaces. Here we describe the structural evolution of water on the stepped Pt(211) surface using thermal desorption, low-energy electron diffraction, and scanning tunneling microscopy to probe the water structure. At low coverage water forms linear structures comprising zigzag chains along the steps that are decorated by H-bonded rings every one or two units along the terrace. Simple 2-coordinate H-bonded chains are not observed, indicating the Pt step binds too weakly to compensate entirely for a low water H-bond coordination number. As the coverage increases, water chains assemble into a disordered (2 × 1) structure, likely made up of the same narrow water chains along the steps with little or no H-bonding between adjacent structures. The chain structure disappears as water adsorption saturates the surface to form an incommensurate, disordered network of water rings of different size. Although the steps on Pt(211) clearly stabilize water adsorption and direct growth, the surface does not support the simple 1D chains previously proposed or an ordered 2D network such as seen on other surfaces. We discuss reasons for this and the factors that determine the behavior of the first water layer on stepped metal surfaces.

Interfacial Layer Breaker: A Violation of Stokes' Law in High-Speed Atomic Force Microscope Flows

Structured water near surfaces is important in nonclassical crystallization, biomineralization, and restructuring of cellular membranes. In addition to equilibrium structures, studied by atomic force microscopy (AFM), high-speed AFM (H-S AFM) can now detect piconewton forces in microseconds. With increasing speeds and decreasing tip diameters, there is a danger that continuum water models will not hold, and molecular dynamic (MD) simulations would be needed for accurate predictions. MD simulations, however, can only evolve over tens of nanoseconds due to memory and computational efficiency/speed limitations, so new methods are needed to bridge the gap. Here, we report a hybrid, multiscale simulation method, which can bridge the size and time scale gaps to existing experiments. Structured water is studied between a moving silica AFM colloidal tip and a cleaved mica surface. The computational domain includes 1,472,766 atoms. To mimic the effect of long-range hydrodynamic forces occurring in water, when moving the AFM tip at speeds from 5 × 10^-7 to 30 m/s, a hybrid multiscale method with local atomistic resolution is used, which serves as an effective open-domain boundary condition. The multiscale simulation is thus equivalent to using a macroscopically large computational domain with equilibrium boundary conditions. Quantification of the drag force shows the breaking of continuum behavior. Nonmonotonic dependence on both the tip speed and distance from the surface implies breaking of the hydration layer around the moving tip at time scales smaller than water cluster formation and strong water compressibility effects at the highest speeds.

Submolecular Insights into Interfacial Water by Hydrogen-Sensitive Scanning Probe Microscopy

ConspectusWater-solid interfaces have attracted extensive attention because of their crucial roles in a wide range of chemical and physical processes, such as ice nucleation and growth, dissolution, corrosion, heterogeneous catalysis, and electrochemistry. To understand these processes, enormous efforts have been made to obtain a molecular-level understanding of the structure and dynamics of water on various solid surfaces. By the use of scanning probe microscopy (SPM), many remarkable structures of H-bonding networks have been directly visualized, significantly advancing our understanding of the delicate competition between water-water and water-solid interactions. Moreover, the detailed dynamics of water molecules, such as diffusion, clustering, dissociation, and intermolecular and intramolecular proton transfer, have been investigated in a well-controlled manner by tip manipulation. However, resolving the submolecular structure of surface water has remained a great challenge for a long time because of the small size and light mass of protons. Discerning the position of hydrogen in water is not only crucial for the accurate determination of the structure of H-bonding networks but also indispensable in probing the proton transfer dynamics and the quantum nature of protons.In this Account, we focus on the recent advances in the H-sensitive SPM technique and its applications in probing the structures, dynamics, and nuclear quantum effects (NQEs) of surface water and ion hydrates at the submolecular level. First, we introduce the development of high-resolution scanning tunneling microscopy/spectroscopy (STM/S) and qPlus-based atomic force microscopy (qPlus-AFM), which allow access to the degrees of freedom of protons in both real and energy space. qPlus-AFM even allows imaging of interfacial water in a weakly perturbative manner by measuring the high-order electrostatic force between the CO-terminated tip and the polar water molecule, which enables the subtle difference of OH directionality to be discerned. Next we showcase the applications of H-sensitive STM/AFM in addressing several key issues related to water-solid interfaces. The surface wetting behavior and the H-bonding structure of low-dimensional ice on various hydrophilic and hydrophobic solid surfaces are characterized at the atomic scale. Then we discuss the quantitative assessment of NQEs of surface water, including proton tunneling and quantum delocalization. Moreover, the weakly perturbative and H-sensitive SPM technique can be also extended to investigations of water-ion interactions on solid surfaces, revealing the effect of hydration structure on the interfacial ion transport. Finally, we provide an outlook on the further directions and challenges for SPM studies of water-solid interfaces.

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.

Impact of surface nanostructure and wettability on interfacial ice physics

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.

Two-dimensional monolayer salt nanostructures can spontaneously aggregate rather than dissolve in dilute aqueous solutions

It is well known that NaCl salt crystals can easily dissolve in dilute aqueous solutions at room temperature. Herein, we reported the first computational evidence of a novel salt nucleation behavior at room temperature, i.e., the spontaneous formation of two-dimensional (2D) alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement. Microsecond-scale classical molecular dynamics (MD) simulations showed that NaCl or LiCl, initially fully dissolved in confined water, can spontaneously nucleate into 2D monolayer nanostructures with either ordered or disordered morphologies. Notably, the NaCl nanostructures exhibited a 2D crystalline square-unit pattern, whereas the LiCl nanostructures adopted non-crystalline 2D hexagonal ring and/or zigzag chain patterns. These structural patterns appeared to be quite generic, regardless of the water and ion models used in the MD simulations. The generic patterns formed by 2D monolayer NaCl and LiCl nanostructures were also confirmed by ab initio MD simulations. The formation of 2D salt structures in dilute aqueous solution at room temperature is counterintuitive. Free energy calculations indicated that the unexpected spontaneous salt nucleation behavior can be attributed to the nanoscale confinement and strongly compressed hydration shells of ions.

Aqueous solutions under nanoscale confinement exhibit interesting physicochemical properties. This work reports evidence on the spontaneous formation of two-dimensional alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement by computer simulations.

Advances in Atomic Force Microscopy: Imaging of Two- and Three-Dimensional Interfacial Water

Interfacial water is closely related to many core scientific and technological issues, covering a broad range of fields, such as material science, geochemistry, electrochemistry and biology. The understanding of the structure and dynamics of interfacial water is the basis of dealing with a series of issues in science and technology. In recent years, atomic force microscopy (AFM) with ultrahigh resolution has become a very powerful option for the understanding of the complex structural and dynamic properties of interfacial water on solid surfaces. In this perspective, we provide an overview of the application of AFM in the study of two dimensional (2D) or three dimensional (3D) interfacial water, and present the prospect and challenges of the AFM-related techniques in experiments and simulations, in order to gain a better understanding of the physicochemical properties of interfacial water.

Water Dissociation and Hydroxyl Formation on Ni(110)

Nickel is an active catalyst for hydrogenation and re-forming reactions, with the reactions showing a strong dependence on the surface exposed. Here, we describe the mixed hydroxyl-water phases formed during water dissociation on Ni(110) using scanning tunneling microscopy and low-current low-energy electron diffraction. Water dissociation starts between 150 and 180 K as the H-bond structure evolves from linear one-dimensional (1D) chains of intact water into a two-dimensional (2D) network containing short rows of face-sharing hexagonal rings. As further water desorbs, the hexagonal rows adopt a local (2 × 3) arrangement, forming small, disordered domains separated by strain relief features. Decomposition of this phase occurs near 220 K to form linear 1D structures consisting of flat, zigzag water chains, with each water stabilized by donating one H to hydroxyl to form a branched chain structure. The OH-H₂O chains repel each other, with the saturation layer ordering into a (2 0, 1 4) structure that decomposes to OH near 245 K as further water desorbs. The structure of the mixed OH/H₂O phases is discussed and contrasted with those found on the related Cu(110) surface, with the differences attributed to strain in the 2D H-bond network caused by the short Ni lattice spacing and strong bond to OH/H₂O.

Predicting heterogeneous ice nucleation with a data-driven approach

Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.

A Helical Monolayer Ice

Understanding water wetting layers on solid surfaces is essential for many natural and industrial processes. Here we find a helical ice monolayer with every six water molecules helically arranged along the normal of the basal plane by performing an intensive structural search based on ab initio calculations. The helical ice is more stable than all previous models of monolayer and bilayer ices in a wide range of water densities both in vacuum and on weakly interacting substrates due to a stronger network of hydrogen bonds enabled by the helical geometry. More compelling is the fact that this model adequately explains a recent experimental ice monolayer grown on graphite in terms of the lattice parameter, water density, and Moiré pattern. The helical character in the new ice model echoes previously reported helical motifs in one-dimensional ice structures and suggests an unexpected capability of hydrogen bonds in driving the surface reconstruction of ice structures.

Formation of Linear Water Chains on Ni(110)

Materials that bind strongly to water structure the contact layer, modifying its chemical and physical properties in a manner that depends on the symmetry and reactivity of the surface. Although detailed models have been developed for several inert surfaces, much less is known about reactive surfaces, particularly those with a symmetry different from that of ice. Here we investigate water adsorption on a rectangular surface, Ni(110), an active re-forming catalyst that interacts strongly with water. Instead of forming a network of H-bonded cyclic rings, water forms flat 1D water chains, leaving half the Ni atoms exposed. Second layer water also follows the surface symmetry, forming chains of alternating pentamer and heptamer rings in preference to an extended 2D structure. This behavior is different from that found on other surfaces studied previously and is driven by the short lattice spacing of the solid and the strength of the Ni-water bond.

Advances in Atomic Force Microscopy: Weakly Perturbative Imaging of the Interfacial Water

The structure and dynamics of interfacial water, determined by the water-interface interactions, are important for a wide range of applied fields and natural processes, such as water diffusion (Kim et al., 2013), electrochemistry (Markovic, 2013), heterogeneous catalysis (Over et al., 2000), and lubrication (Zilibotti et al., 2013). The precise understanding of water-interface interactions largely relies on the development of atomic-scale experimental techniques (Guo et al., 2014) and computational methods (Hapala et al., 2014b). Scanning probe microscopy has been extensively applied to probe interfacial water in many interdisciplinary fields (Ichii et al., 2012; Shiotari and Sugimoto, 2017; Peng et al., 2018a). In this perspective, we review the recent progress in the noncontact atomic force microscopy (nc-AFM) imaging and AFM simulation techniques and discuss how the newly developed techniques are applied to study the properties of interfacial water. The nc-AFM with the quadrupole-like CO-terminated tip can achieve ultrahigh-resolution imaging of the interfacial water on different surfaces, trace the reconstruction of H-bonding network and determine the intrinsic structures of the weakly bonded water clusters and even their metastable states. In the end, we present an outlook on the directions of future AFM studies of interfacial water as well as the challenges faced by this field.

Imaging and Dynamics of Water Hexamer Confined in Nanopores

Epitaxial two-dimensional (2D) nanostructures with regular patterns show great promise as templates for adsorbate confinement. Prospectively, employing 2D semiconductors with reduced density of states leads to a long excited-state lifetime that allows us to directly image the dynamics of the adsorbate. We show that epitaxial blue phosphorene (blueP) on Au(111) provides such a platform to trap water molecules in the periodic nanopores without formation of strong bonds. The trapped water aggregate is tentatively assigned to a hexamer based on our scanning tunneling microscopy studies and first-principles calculations. Real-space observation of conformational switching of the hexamer induced by inelastic electrons is achieved by using low-temperature scanning tunneling microscopy with molecular resolution. We found a localized interfacial charge rearrangement between the water hexamer and P atoms underneath that is responsible for the reversible desorption and adsorption of water molecules by changing the sample bias polarity from positive to negative, offering a promising strategy for engineering the electronic properties of blueP.