Atomic imaging of the edge structure and growth of a two-dimensional hexagonal ice

Ultrahigh resistance of hexagonal boron nitride to mineral scale formation

Formation of mineral scale on a material surface has profound impact on a wide range of natural processes as well as industrial applications. However, how specific material surface characteristics affect the mineral-surface interactions and subsequent mineral scale formation is not well understood. Here we report the superior resistance of hexagonal boron nitride (hBN) to mineral scale formation compared to not only common metal and polymer surfaces but also the highly scaling-resistant graphene, making hBN possibly the most scaling resistant material reported to date. Experimental and simulation results reveal that this ultrahigh scaling-resistance is attributed to the combination of hBN’s atomically-smooth surface, in-plane atomic energy corrugation due to the polar boron-nitrogen bond, and the close match between its interatomic spacing and the size of water molecules. The latter two properties lead to strong polar interactions with water and hence the formation of a dense hydration layer, which strongly hinders the approach of mineral ions and crystals, decreasing both surface heterogeneous nucleation and crystal attachment.

Scale formation may have detrimental effects on the properties and functions of materials’ surfaces. Here the authors report the high scaling resistance of hexagonal boron nitride and relate it to the atomic level structure and interaction with water molecules.

Effects of ice and supercooled water on the metastability of methane hydrate: DSC analysis and MD simulations

Methane hydrate (MH) has been viewed as a potential abundant clean energy resource worldwide. Its related technologies play important roles in applications of gas and energy storage, flow assurance of natural gas pipelines etc. Unlike the well-researched stability and decomposition of MH at temperatures above 273 K, the metastability of MH below the ice freezing point, i.e. the anomalous slow decomposition out of thermodynamically stable regions, remains to be unravelled. Studies regarding the influences of ice and supercooled water (SW) on the metastable properties of MH led to varied conclusions, i.e. the as-proposed self-preservation effect and metastable MH-SW-gas equilibrium. In this study, a series of DSC experiments were performed to investigate the thermal stability boundaries and the associated metastable behaviours of MH-ice-gas and MH-SW-gas samples in porous medium. The DSC analysis probed accurate thermal stabilities and characterized decomposition behaviors of the samples, contributing to the hypothesis of potential influences from SW and ice on the metastability of MH. MD simulations were also validated and performed. Active guest-host interactions by the SW layers between MH and gas phases were identified, suggesting probable microscopic configurations related to the metastability of the MH-SW-gas system. Indications of the DSC and MD simulation results call for future high-resolution in situ experimental validations.

Robustness of Bilayer Hexagonal Ice against Surface Symmetry and Corrugation

Two-dimensional (2D) bilayer hexagonal ice (BHI) is regarded as the first intrinsic 2D ice crystal. However, the robustness of such a structure or its derivatives against surface symmetry and corrugation is still unclear. Here, we report the formation of 2D BHI on gold surfaces with 1D corrugation, using noncontact atomic force microscopy. The hexagonal arrangement of the first wetting layer was visualized on the Au(110)-1×2 surface. Upon depositing more water molecules, the first layer would rearrange and shrink, resulting in the formation of buckled BHI. Such a buckled BHI is hydrophobic despite the appearance of dangling OH, due to the strong interlayer bonding. Furthermore, the BHI is also stable on the Au(100)-5×28 surface. This work reveals the unexpected generality of the BHI on corrugated surfaces with nonhexagonal symmetry, thus shedding new light on the microscopic understandings of the low-dimensional ice formation on solid surfaces or under confinement.

Visualizing Eigen/Zundel cations and their interconversion in monolayer water on metal surfaces

The nature of hydrated proton on solid surfaces is of vital importance in electrochemistry, proton channels, and hydrogen fuel cells but remains unclear because of the lack of atomic-scale characterization. We directly visualized Eigen- and Zundel-type hydrated protons within the hydrogen bonding water network on Au(111) and Pt(111) surfaces, using cryogenic qPlus-based atomic force microscopy under ultrahigh vacuum. We found that the Eigen cations self-assembled into monolayer structures with local order, and the Zundel cations formed long-range ordered structures stabilized by nuclear quantum effects. Two Eigen cations could combine into one Zundel cation accompanied with a simultaneous proton transfer to the surface. Moreover, we revealed that the Zundel configuration was preferred over the Eigen on Pt(111), and such a preference was absent on Au(111).

Seeing how ice breaks the rule

Basic defects in ice monolayers are seen using a microscope

Computational Analysis of Hydrogen Bond Vibrations of Ice III in the Far-Infrared Band

The hydrogen-disordered structure of ice III makes it difficult to analyze its vibrational spectrum theoretically. To clarify the contribution of hydrogen bonds (HBs), we constructed a 24-molecule supercell to mimic the real structure and performed first-principles density functional theory calculations. The calculated curve of phonon density of states showed good correspondence with the experimental data. Based on the theory of two kinds of HB vibrational modes, we analyzed the distributions of two-bond modes and four-bond modes. The energy splitting of these modes results in a flat vibrational band, which is a common phenomenon in high-pressure ice phases. These findings verified the general rule that there are two types of HB vibrations in ice, thereby furthering our understanding of HB interactions in water ice and their broad role in nature.

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.

A qPlus-based scanning probe microscope compatible with optical measurements

We design and develop a scanning probe microscope (SPM) system based on the qPlus sensor for atomic-scale optical experiments. The microscope operates under ultrahigh vacuum and low temperature (6.2 K). In order to obtain high efficiency of light excitation and collection, two front lenses with high numerical apertures (N.A. = 0.38) driven by compact nano-positioners are directly integrated on the scanner head without degrading its mechanical and thermal stability. The electric noise floor of the background current is 5 fA/Hz^1/2, and the maximum vibrational noise of the tip height is below 200 fm/Hz^1/2. The drift of the tip-sample spacing is smaller than 0.1 pm/min. Such a rigid scanner head yields small background noise (oscillation amplitude of ∼2 pm without excitation) and high quality factor (Q factor up to 140 000) for the qPlus sensor. Atomic-resolution imaging and inelastic electron tunneling spectroscopy are obtained under the scanning tunneling microscope mode on the Au(111) surface. The hydrogen-bonding structure of two-dimensional (2D) ice on the Au(111) surface is clearly resolved under the atomic force microscope (AFM) mode with a CO-terminated tip. Finally, the electroluminescence spectrum from a plasmonic AFM tip is demonstrated, which paves the way for future photon-assisted SPM experiments.

Anomalous Phase Behaviors of Monolayer NaCl Aqueous Solutions Induced by Effective Coulombic Interactions within Angstrom-Scale Slits

Interests in subnanofluidic devices have called for molecular dynamics (MD) simulation studies of the thermodynamic behavior of monolayer salt solution within angstrom-scale slits. However, it still remains a grand challenge to accurately describe the Coulombic interactions by incorporating the effects of charge transfer and electronic dielectric screening. Herein, by using the electronic continuum model, where the effective ion charges are fine-tuned with a scaling factor of λ, we present simulation evidence that the effective Coulombic interactions among Na⁺/Cl^- ions can strongly affect the behavior of monolayer ionic aqueous solution. Our microsecond-scale MD simulations show that only the counterions with moderate effective charges (0.3 ≤ λ ≤ 0.8) can dissolve in monolayer water, whereas the high effective charges (λ ≥ 0.85) induce ions to assemble into monolayer nanocrystals, and ions with the low effective charges (λ ≤ 0.2) exhibit gas-like nanobubble. These findings could provide deeper insights into the physical chemistry behind subnanofluidic iontronic devices.

Insights into Ionic Liquids: From Z-Bonds to Quasi-Liquids

Ionic liquids (ILs) hold great promise in the fields of green chemistry, environmental science, and sustainable technology due to their unique properties, such as a tailorable structure, the various types available, and their environmentally friendly features. On the basis of multiscale simulations and experimental characterizations, two unique features of ILs are as follows: (1) strong coupling interactions between the electrostatic forces and hydrogen bonds, namely in the Z-bond, and (2) the unique semiordered structure and properties of ultrathin films, specifically regarding the quasi-liquid. In accordance with the aforementioned theoretical findings, many cutting-edge applications have been proposed: for example, CO₂ capture and conversion, biomass conversion and utilization, and energy storage materials. Although substantial progress has been made recently in the field of ILs, considerable challenges remain in understanding the nature of and devising applications for ILs, especially in terms of e.g. in situ/real-time observation and highly precise multiscale simulations of the Z-bond and quasi-liquid. In this Perspective, we review recent developments and challenges for the IL research community and provide insights into the nature and function of ILs, which will facilitate future applications.

Local Chiral Inversion of Thymine Dimers by Manipulating Single Water Molecules

Water, as one of the most important and indispensable small molecules in vivo, plays a crucial role in driving biological self-assembly processes. Real-space detection and identification of water-induced organic structures and further capture of dynamic dehydration processes are important yet challenging, which would help to reveal the cooperation and competition mechanisms among water-involved noncovalent interactions. Herein, introduction of water molecules onto the self-assembled thymine (T) structures under ultrahigh vacuum (UHV) conditions results in the hydration of hydrogen-bonded T dimers forming a well-ordered water-involved T structure. Reversibly, a local dehydration process is achieved by in situ scanning tunneling microscopy (STM) manipulation on single water molecules, where the adjacent T dimers connected with water molecules undergo a local chiral inversion process with the hydrogen-bonding configuration preserved. Such a strategy enables real-space identification and detection of the interactions between water and organic molecules, which may also shed light on the understanding of biologically relevant self-assembly processes driven by water.

Graphitic-like Hexagonal Phase of Alkali Halides in Quasi-Two-Dimensional Confined Space under Ambient Conditions

The discovery of specific matter phases with abnormal physical properties in low-dimensional systems and/or on particular substrates, such as the hexagonal phase of ice and two-dimensional (2D) CaCl with an abnormal valence state, continuously reveals more fundamental mechanisms of the nature. Alkali halides, represented by NaCl, are one of the most common compounds and usually thought to be well-understood. In the past decades, many theoretical studies suggested the existence of one particular phase, that is, the graphitic-like hexagonal phase of alkali halides at high pressure or in low-dimension states, with the expectation of improved properties of this matter phase but lacking experimental evidence due to severe technical challenges. Here, by optimized cryo-electron microscopy, we report the direct atomic-resolution observation and in situ characterization of the prevalent and stable graphitic-like alkali halide hexagonal phases, which were spontaneously formed by unsaturated NaCl and LiCl solution, respectively, in the quasi-2D confined space between reduced graphene oxide layers under ambient conditions. Combined with a control experiment, density functional theory calculations, and previous theoretical studies, we believe that a delicate balance among the cation-π interaction of the solute and substrate, electrostatic interactions of anions and cations, solute-solvent interactions, and thermodynamics under confinement synergistically results in the formation of such hexagonal crystalline phases. These findings highlight the effects of the substrate and the confined space on the formation of specific matter phases and provide a universal scheme for the preparation of special graphitic-like hexagonal phases of alkali halides.

Patterning Configuration of Surface Hydrophilicity by Graphene Nanosheet towards the Inhibition of Ice Nucleation and Growth

Freezing of liquid water occurs in many natural phenomena and affects countless human activities. The freezing process mainly involves ice nucleation and continuous growth, which are determined by the energy and structure fluctuation in supercooled water. Herein, considering the surface hydrophilicity and crystal structure differences between metal and graphene, we proposed a kind of surface configuration design, which was realized by graphene nanosheets being alternately anchored on a metal substrate. Ice nucleation and growth were investigated by molecular dynamics simulations. The surface configuration could induce ice nucleation to occur preferentially on the metal substrate where the surface hydrophilicity was higher than the lateral graphene nanosheet. However, ice nucleation could be delayed to a certain extent under the hindering effect of the interfacial water layer formed by the high surface hydrophilicity of the metal substrate. Furthermore, the graphene nanosheets restricted lateral expansion of the ice nucleus at the clearance, leading to the formation of a curved surface of the ice nucleus as it grew. As a result, ice growth was suppressed effectively due to the Gibbs–Thomson effect, and the growth rate decreased by 71.08% compared to the pure metal surface. Meanwhile, boundary misorientation between ice crystals was an important issue, which also prejudiced the growth of the ice crystal. The present results reveal the microscopic details of ice nucleation and growth inhibition of the special surface configuration and provide guidelines for the rational design of an anti-icing surface.

Studying Ice with Environmental Scanning Electron Microscopy

Scanning electron microscopy (SEM) is a powerful imaging technique able to obtain astonishing images of the micro- and the nano-world. Unfortunately, the technique has been limited to vacuum conditions for many years. In the last decades, the ability to introduce water vapor into the SEM chamber and still collect the electrons by the detector, combined with the temperature control of the sample, has enabled the study of ice at nanoscale. Astounding images of hexagonal ice crystals suddenly became real. Since these first images were produced, several studies have been focusing their interest on using SEM to study ice nucleation, morphology, thaw, etc. In this paper, we want to review the different investigations devoted to this goal that have been conducted in recent years in the literature and the kind of information, beyond images, that was obtained. We focus our attention on studies trying to clarify the mechanisms of ice nucleation and those devoted to the study of ice dynamics. We also discuss these findings to elucidate the present and future of SEM applied to this field.

Hydrogen bonded trimesic acid networks on Cu(111) reveal how basic chemical properties are imprinted in HR-AFM images

High resolution non-contact atomic force microscopy measurements characterize assemblies of trimesic acid molecules on Cu(111) and the link group interactions, providing the first fingerprints utilizing CO-based probes for this widely studied paradigm for hydrogen bond driven molecular self assembly. The enhanced submolecular resolution offered by this technique uniquely reveals key aspects of the competing interactions. Accurate comparison between full-density-based modeled images and experiment allows to identify key structural elements in the assembly in terms of the electron-withdrawing character of the carboxylic groups, interactions of those groups with Cu atoms in the surface, and the valence electron density in the intermolecular region of the hydrogen bonds. This study of trimesic acid assemblies on Cu(111) combining high resolution atomic force microscopy measurements with theory and simulation forges clear connections between fundamental chemical properties of molecules and key features imprinted in force images with submolecular resolution.

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.

Coexistence of Superconductivity and Nontrivial Band Topology in Monolayered Cobalt Pnictides on SrTiO3

As an intrinsically layered material, FeSe has been extensively explored for potentially revealing the underlying mechanisms of high transition temperature (high-T_c) superconductivity and realizing topological superconductivity and Majorana zero modes. Here we use first-principles approaches to identify that the cobalt pnictides of CoX (X = As, Sb, Bi), none of which is a layered material in bulk form. Nevertheless, all can be stabilized as monolayered systems either in freestanding form or supported on the SrTiO₃(001) substrate. We further show that each of the cobalt pnictides may potentially harbor high-T_c superconductivity beyond the Cu- and Fe-based superconducting families, and the underlying mechanism is inherently tied to their isovalency nature with the FeSe monolayer. Most strikingly, each of the monolayered CoX's on SrTiO₃ is shown to be topologically nontrivial, and our findings provide promising new platforms for realizing topological superconductors in the two-dimensional limit.

Important Roles of Water Clusters Confined in a Nanospace as Revealed by a Synchrotron X-ray Diffraction Study

States of water molecules confined in a nanospace designed by montmorillonite (negatively charged silicate layer) and charge compensating benzylammonium were investigated. Caffeine was used as a probe because of its compatibility for the fine structure of the interlayer water. Powder synchrotron radiation X-ray diffraction (SXRD) and adsorption isotherms of the water vapor revealed a metastable structure of bimolecular water layers (2WLs) in the interlayer space. Water molecules readily penetrated to expand the interlayer space to 0.56 nm. The interlayer space did not increase further even in the presence of excess water. According to the isosteric heat of water, the expansion was limited because of moderate hydration as forming 2WLs. Caffeine molecules replaced a part of the water molecules in the 2WLs to expand the interlayer space to 0.65 nm. Time-resolved SXRD with an accumulation time of 500 ms revealed that the interlayer expansion reached a steady state within a few minutes. The caffeine intercalation proceeded, involving a change in the molecular orientation that increased the contact area of the caffeine molecules. The interlayer expansion was limited in all the solvents examined (mixtures of water with methanol, ethanol, acetone, and tetrahydrofuran), while the packing density of the incorporated caffeine was maximized in the absence of an organic solvent. The water molecules confined in the interlayer space acted as an actuator to accommodate a large quantity of amphiphilic molecules by adapting the nanostructure, which was achieved by releasing the confined water molecules.

Temperature-dependent kinetic pathways of heterogeneous ice nucleation competing between classical and non-classical nucleation

Ice nucleation on the surface plays a vital role in diverse areas, ranging from physics and cryobiology to atmospheric science. Compared to ice nucleation in the bulk, the water-surface interactions present in heterogeneous ice nucleation complicate the nucleation process, making heterogeneous ice nucleation less comprehended, especially the relationship between the kinetics and the structures of the critical ice nucleus. Here we combine Markov State Models and transition path theory to elucidate the ensemble pathways of heterogeneous ice nucleation. Our Markov State Models reveal that the classical one-step and non-classical two-step nucleation pathways can surprisingly co-exist with comparable fluxes at T = 230 K. Interestingly, we find that the disordered mixing of rhombic and hexagonal ice leads to a favorable configurational entropy that stabilizes the critical nucleus, facilitating the non-classical pathway. In contrast, the favorable energetics promotes the formation of hexagonal ice, resulting in the classical pathway. Furthermore, we discover that, at elevated temperatures, the nucleation process prefers to proceed via the classical pathway, as opposed to the non-classical pathway, since the potential energy contributions override the configurational entropy compensation. This study provides insights into the mechanisms of heterogeneous ice nucleation and sheds light on the rational designs to control crystallization processes.

Formation of porous ice frameworks at room temperature

Bulk crystalline ices with ultralow densities have been demonstrated to be thermodynamically metastable at negative pressures. However, the direct formation of these bulk porous ices from liquid water at negative pressures is extremely challenging. Inspired by approaches toward porous media based on host-guest chemistry, such as metal-organic frameworks and covalent organic frameworks, we herein demonstrate via molecular dynamics simulations that a class of ultralow-density porous ices with upright channels can be formed spontaneously from liquid water at 300 K with the assistance of carbon nanotube arrays. We refer to these porous ice polymorphs as water oxygen-vertex frameworks (WOFs). Notably, our simulations revealed that the liquid-WOF phase transition is first-order and occurs at room temperature. All the WOFs exhibited the unique structural feature that they can be regarded as assemblies of nanoribbons of hexagonal bilayer ice (2D ice I) at their armchair or zigzag edges. Based on density functional theory calculations, a comprehensive phase diagram of the WOFs was constructed considering both the thermodynamic and thermal stabilities of the porous ices at negative pressures. Like other types of porous media, these WOFs may be applicable to gas storage, purification, and separation. Moreover, these biocompatible porous ice networks may be exploited as medical-related carriers.

A high-accuracy machine-learning water model for exploring water nanocluster structures

Water, the most important molecule on the Earth, possesses many essential and unique physical properties that are far from completely understood, partly due to serious difficulties in identifying the precise microscopic structures of water. Hence, identifying the structures of water nanoclusters is a fundamental and challenging issue for studies on the relationship between the macroscopic physical properties of water and its microscopic structures. For large-scale simulations (at the level of nm and ns) of water nanoclusters, a calculation method with simultaneous accuracy at the level of quantum chemistry and efficiency at the level of an empirical potential method is in great demand. Herein, a machine-learning (ML) water model was utilized to explore the microscopic structural features at different length scales for water nanoclusters with a size up to several nm. The ML water model can be employed to efficiently predict the structures of water nanoclusters with a similar accuracy to that of density functional theory and with substantially lower computational resource demands. To validate the low-lying structure search results with experimental spectral results, an ML water model combined with velocity autocorrelation function analysis was used to simulate the vibrational spectra of water nanoclusters with up to thousands of water molecules. By comparing the simulated and experimentally recorded vibrational spectra, the atomic structures determined by a simulation based on the ML water model are all verified. To demonstrate its ability to represent water's structural evolution at large length and time scales, the ML water model was employed to model the structural evolution during the crystal-liquid transition, and the phase transition temperatures of water clusters with different sizes were precisely predicted. The ML water model provides an efficient theoretical calculation tool for exploring the structures and physical properties of water and their relationships, especially for clusters with relatively large sizes and processes with relatively long durations.

Edge premelting of two-dimensional ices

The surface of a three-dimensional ice crystal naturally has a quasi-liquid layer (QLL) at temperatures below its bulk melting point, due to a phenomenon called surface premelting. Here, we show that the edges of a two-dimensional (2D) bilayer hexagonal ice adsorbed on solid surfaces undergo premelting as well, resulting in the formation of quasi-liquid bands (QLBs) at the edges. Our extensive molecular dynamics simulations show that the QLB exhibits structure and dynamics indistinguishable from the bilayer liquid phase, acting as a lower-dimensional analog of the QLL on the bulk ice. We further find that at low temperatures, the width of the QLBs at armchair-type edges of the 2D ice is almost identical to that at zigzag-type edges but becomes far greater than the latter at temperatures near the melting point. The chirality-dependent edge premelting of 2D ices should add an important new ingredient to the heterogeneity of premelting.

Hydrogen Bonding and the Structural Properties of Glycerol-Water Mixtures with a Microwave Field: a Molecular Dynamics Study

In a microwave field, the dielectric properties, molecular structures, and hydrogen bonding dynamics of glycerol in its mixtures with water were determined by the molecular dynamics simulation method. The dipole-dipole correlation of glycerol is linked to the field intensity of microwaves. The results show that as the field intensity is increased, even glycerol in the second coordination shell can become correlated with each other. The structures of up to 35 glycerol molecules are observed. More than that, it was observed that lifetimes of glycerol-glycerol hydrogen bonds were prolonged, while the average hydrogen bond number was also increased. Besides, the structures in a strong microwave field mimic the weak C-H⋯O hydrogen bonds seen in high-glycerol concentration mixtures, yet the concentration is lower. These results indicate that with the assistance of the microwave field, glycerol molecules become concentrated and are more likely to establish stable interactions with others. As a consequence, the spherical clusters composed by glycerol molecules in our nanosheet synthesis experiment are easier to form.

First-Principles Molecular Dynamics Simulations of the Spontaneous Freezing Transition of 2D Water in a Nanoslit

As with bulk ices, two-dimensional (2D) ices exhibit diverse crystalline structures, and the majority of these 2D structures have been predicted based on classical molecular dynamics (MD) simulations. Here, the spontaneous freezing transition of 2D liquid water within hydrophobic nanoslits is demonstrated for the first time using first-principles MD simulations. Various 2D ices are observed under different lateral pressure and temperature conditions. Notably, the liquid water confined to a 6.0 Å-wide nanoslit can spontaneously freeze into a monolayer ice consisting of an array of zigzag water chains at 2.5 GPa and 250 K. Moreover, within an 8.0 Å-wide nanoslit and at 4.0 GPa and 300 K, a previously unreported bilayer ice forms spontaneously that has a structure resembling that of the double surface layers of bulk ice-VII. Both 2D crystalline ices do not obey the ice rule, suggesting first-principles simulation can access a certain phase space that is not easily approached using classical simulations.

A Computational Validation of Water Molecules Adsorption on an NaCl Surface

It was reported that a scanning tunneling microscopy (STM) study observed the adsorption geometry of a water monomer and a tetramer on NaCl(100) film. Based on first-principles density functional theory (DFT), the adsorption behavior of water on the NaCl surface was simulated with CASTEP code. The results showed that the water monomer almost lay on the NaCl(001) surface with one O–H bond tilted slightly downward. This was quite different from the STM observations. In fact, the experimental observation was influenced by the Au(111) substrate, which showed an upright form. A recent report on observations of two-dimensional ice structure on Au(111) substrate verified our simulations. However, the water tetramer formed a stable quadrate structure on the surface, which was consistent with observation. The intermolecular hydrogen bonds present more strength than surface adsorption. The simulations presented a clearer picture than experimental observations.

Gas hydrates in confined space of nanoporous materials: new frontier in gas storage technology

Gas hydrates (clathrate hydrates, clathrates, or hydrates) are crystalline inclusion compounds composed of water and gas molecules. Methane hydrates, the most well-known gas hydrates, are considered a menace in flow assurance. However, they have also been hailed as an alternative energy resource because of their high methane storage capacity. Since the formation of gas hydrates generally requires extreme conditions, developing porous material hosts to synthesize gas hydrates with less-demanding constraints is a topic of great interest to the materials and energy science communities. Though reports of modeling and experimental analysis of bulk gas hydrates are plentiful in the literature, reliable phase data for gas hydrates within confined spaces of nanoporous media have been sporadic. This review examines recent studies of both experiments and theoretical modeling of gas hydrates within four categories of nanoporous material hosts that include porous carbons, metal-organic frameworks, graphene nanoslits, and carbon nanotubes. We identify challenges associated with these porous systems and discuss the prospects of gas hydrates in confined space for potential applications.

Semiconducting MnB5monolayer as a potential photovoltaic material

Exploring new two-dimensional (2D) materials is of great significance for both basic research and practical applications. Although boron can form various 3D and 2D allotropes due to its ease of forming multi-center bonds, the coexistence of honeycomb and kagome boron structures has never been observed in any 2D material yet. In this article we apply first-principle swarm structural searches to predict the existence of a stable MnB₅structure, consisting of a sandwich of honeycomb and kagome borophenes. More interestingly, a MnB₅nanosheet is a semiconductor with a band gap of 1.07 eV and a high optical absorption in a broad band, which satisfies the requirements of a very good photovoltaic material. Upon moderate strain, MnB₅undergoes a conversion from an indirect to a direct band gap semiconductor. The power conversion efficiency of a heterostructure solar cell made of MnB₅is up to 18%. The MnB₅nanosheet shows a robust dynamical and thermal stability, stemming from the presence of intra- and interlayer multi-center σ and π bonds. These characteristics make MnB₅a promising photovoltaic material.

Mechanical hydrolysis imparts self-destruction of water molecules under steric confinement

Decoding behavioral aspects associated with the water molecules in confined spaces such as an interlayer space of two-dimensional nanosheets is key for the fundamental understanding of water-matter interactions and identifying unexpected phenomena of water molecules in chemistry and physics. Although numerous studies have been conducted on the behavior of water molecules in confined spaces, their reach stops at the properties of the planar ice-like formation, where van der Waals interactions are the predominant interactions and many questions on the confined space such as the possibility of electron exchange and excitation state remain unsettled. We used density functional theory and reactive molecular dynamics to reveal orbital overlap and induction bonding between water molecules and graphene sheets under much less pressure than graphene fractures. Our study demonstrates high amounts of charge being transferred between water and the graphene sheets, as the interlayer space becomes smaller. As a result, the inner face of the graphene nanosheets is functionalized with hydroxyl and epoxy functional groups while released hydrogen in the form of protons either stays still or traverses a short distance inside the confined space via the Grotthuss mechanism. We found signatures of a new hydrolysis mechanism in the water molecules, i.e. mechanical hydrolysis, presumably responsible for relieving water from extremely confined conditions. This phenomenon where water reacts under extreme confinement by disintegration rather than forming ice-like structures is observed for the first time, illustrating the prospect of treating ultrafine porous nanostructures as a driver for water splitting and material functionalization, potentially impacting the modern design of nanofilters, nanochannels, nano-capacitators, sensors, and so on.

Probing Specific Adsorption of Electrolytes at Kaolinite-Aqueous Interfaces by Atomic Force Microscopy

Adsorption of electrolytes (ions) at solid-liquid interfaces alters the physical and chemical properties of materials and hence plays a critical role in manufacturing and processing of nanomaterials featuring large surface or interfacial areas of desired structures and morphology. Many experiments and theoretical calculations using various electrical double layer (EDL) models have been conducted to understand how and where ions adsorb at charged surfaces in a liquid. However, conclusions from previous research remain inconclusive because of model-dependent approaches to studying ion adsorption at diverse solid-liquid interfaces. In this study, atomic force microscopy is used to image in liquids the surface lattice structure of two kaolinite basal planes in the presence and absence of monovalent and divalent cations. Distinct adsorption of ions through different mechanisms (such as electrostatic attraction and specific adsorption) is identified through atomic resolution imaging without the assumption of an EDL structure.

100 MHz large bandwidth preamplifier and record-breaking 50 kHz scanning rate quantum point contact mode probe microscopy imaging with atomic resolution

The high-bandwidth preamplifier is a vital component designed to increase the scanning speed of a high-speed scanning tunneling microscope (STM). However, the bandwidth is limited not only by the characteristic GΩ feedback resistor R_F but also by the characteristic unity-gain-stable operational amplifier (UGS-OPA) in the STM preamplifier. Here, we report that paralleling a resistor with the tunneling junction (PRTJ) can break both limitations. Then, the UGS-OPA can be replaced by a higher rate, higher antinoise ability, decompensated OPA. By doing so, a bandwidth of more than 100 MHz was achieved in the STM preamplifier with decompensated OPA657, and a higher bandwidth is possible. High-clarity atomic resolution STM images were obtained under about 10 MHz bandwidth and quantum point contact microscopy mode with a record-breaking line rate of 50 k lines/s and a record-breaking frame rate of 250 frames/s. Both the PRTJ method and the decompensated OPA will pave the way for higher scanning speeds and play a key role in the design of high-performance STMs.

A folded ice monolayer

A highly stable ice monolayer with folded structural motifs is predicted by means of a novel tiling method augmented with ab initio calculations. This ice monolayer has every two neighboring water hexamers connected by a water square yet folded into two distinct planes, and is thus coined as a folded ice model. It is in the ground state in a range of water densities from 0.08 to 0.12 Å-2, with a stronger energy preference at a lower water density. Its stability shown by ab initio molecular dynamics simulations can sustain up to a temperature of 100 K. The tiling method also enables the prediction of a family of considerably stable ice monolayers with a variety of puckered structures. These results enrich our knowledge of low-dimensional water structures and pave a way to explore more exotic ice nanostructures under confinements.

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.

Atomic structure causing an obvious difference in thermal conductance at the Pd-H2O interface: a molecular dynamics simulation

Thermal transfer across solid-liquid interfaces is influenced by multiple factors such as surface wettability, interfacial water layer density, molecular structure, and mass density depletion length. However, the dominant factors in interfacial heat transport are yet to be investigated. In this work, we explore the contributions from these factors by employing the Pd-water model for water molecules forming ordered, partially ordered, and disordered structures on Pd (100), (110) and (111) surfaces, respectively. The results revealed that the ordered water layer on the (100) surface can introduce a "phonon bridge" at the solid-liquid interface to improve thermal transfer, while the partially ordered water layer on the (110) surface can further promote thermal transfer due to the enhanced interfacial friction. On the other hand, the decreased density depletion length also makes dominant contributions to the enhancement of interfacial thermal transfer. The results are explained by the interfacial friction coefficient, surface potential energy distribution and density depletion length. We also introduce an efficient technique by tuning the vacancy defects on the solid surface to tune the atomic structure as well as the thermal transfer. Our study reveals the complex relationship between the atomic structure of the crystal face, the water layer structure and the thermal boundary conductance, which will inspire more experimental and theoretical studies toward the improvement of interfacial thermal transport by tuning the structure of the water layer.

Computational Prediction of Novel Ice Phases: A Perspective

Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.

Intrinsic reconstruction of ice-I surfaces

Understanding the precise atomic structure of ice surfaces is critical for revealing the mechanisms of physical and chemical phenomena at the surfaces, such as ice growth, melting, and chemical reactions. Nevertheless, no conclusive structure has been established. In this study, noncontact atomic force microscopy was used to address the characterization of the atomic structures of ice Ih(0001) and Ic(111) surfaces. The topmost hydrogen atoms are arranged with a short-range (2 × 2) order, independent of the ice thickness and growth substrates used. The electrostatic repulsion between non-hydrogen-bonded water molecules at the surface causes a reduction in the number of the topmost hydrogen atoms together with a distortion of the ideal honeycomb arrangement of water molecules, leading to a short-range-ordered surface reconstruction.

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.