Robustness of Bilayer Hexagonal Ice against Surface Symmetry and Corrugation

Improving Molecular-Dynamics Simulations for Solid-Liquid Interfaces with Machine-Learning Interatomic Potentials

Emerging developments in artificial intelligence have opened infinite possibilities for material simulation. Depending on the powerful fitting of machine learning algorithms to first-principles data, machine learning interatomic potentials (MLIPs) can effectively balance the accuracy and efficiency problems in molecular dynamics (MD) simulations, serving as powerful tools in various complex physicochemical systems. Consequently, this brings unprecedented enthusiasm for researchers to apply such novel technology in multiple fields to revisit the major scientific problems that have remained controversial owing to the limitations of previous computational methods. Herein, we introduce the evolution of MLIPs, provide valuable application examples for solid-liquid interfaces, and present current challenges. Driven by solving multitudinous difficulties in terms of the accuracy, efficiency, and versatility of MLIPs, this booming technique, combined with molecular simulation methods, will provide an underlying and valuable understanding of interdisciplinary scientific challenges, including materials, physics, and chemistry.

Two rhombic ice phases from aqueous salt solutions under graphene confinement

Water exhibits rich ice phases depending upon its respective formation conditions, and in particular, the two-dimensional ice with nonhexagonal symmetry adsorbed on solids relates to the exceptional arrangement of water molecules. Despite extensive reporting of two-dimensional ice on various solid surfaces, the geometry and thermodynamics of ice formation from an aqueous salt solution are still unknown. In this Letter, we show the formation of single- and two-phase mixed two-dimensional rhombic ice from aqueous salt solutions with different concentrations under strong compressed confinement of graphene at ambient temperature by using classical molecular dynamics simulations and first-principles calculations. The two rhombic ice phases exhibit identical geometry and thermodynamic properties, but different projections of the oxygen atoms against solid surface symmetry, where they relate to the stable and metastable arrangements of water molecules confined between two graphene layers. A single-phase rhombic ice would grow from the confined saturated aqueous solutions since the previously stable rhombic molecular arrangement becomes an unstable high-energy state by introducing salt ions nearby. Our result reveals different rhombic ice phases growing from pure water and aqueous solutions, highlighting the deciding role of salt ions in the ice formation process due to their common presence in liquids.

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+.

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.

Visualizing the Promoting Role of Interfacial Water in the Deprotonation of Formic Acid on Cu(111)

Water plays a crucial role in various heterogeneous catalytic reactions, but the atomic-scale characterization of how water participates in these chemical processes remains a significant challenge. Here we directly visualize the promoting role of interfacial water in the deprotonation of formic acid (FA) on a metal surface, using combined scanning tunneling microscopy and qPlus-based noncontact atomic force microscopy. We find the dissociation of FA when coadsorbed with water on the Cu(111) surface, resulting in the formation of hydronium and formate ions. Interestingly, most of the hydrated proton and formate ions exhibit a phase-separated behavior on Cu(111), in which Eigen and Zundel cations assemble into a monolayer hexagonal hydrogen-bonding (H-bonding) network, and bidentate formate ions are solvated with water and aggregate into one-dimensional chains or two-dimensional H-bonding networks. This phase-separated behavior is essential for preventing the proton transfer back from hydronium to formate and the reformation of FA. Density functional theory calculations reveal that the participation of water significantly reduces the deprotonation barrier of FA on Cu(111), in which water catalyzes the decomposition of FA through the Grotthuss proton transfer mechanism. In addition, the separate solvation of hydronium and bidentate formate ions is energetically preferred due to the enhanced interaction with the copper substrate. The promoting role of water in the deprotonation of FA is further confirmed by the temperature-programmed desorption experiment, which shows that the intensity of the H2 desorption peak significantly increases and the desorption of FA declines when water and FA coadsorbed on the Cu(111) surface.

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.

Origin of the Hydrophobic Behaviour of Hydrophilic CeO2

The nature of the hydrophobicity found in rare-earth oxides is intriguing. The CeO2 (100) surface, despite its strongly hydrophilic nature, exhibits hydrophobic behaviour when immersed in water. In order to understand this puzzling and counter-intuitive effect we performed a detailed analysis of the confined water structure and dynamics. We report here an ab-initio molecular dynamics simulation (AIMD) study which demonstrates that the first adsorbed water layer, in immediate contact with the hydroxylated CeO2 surface, generates a hydrophobic interface with respect to the rest of the liquid water. The hydrophobicity is manifested in several ways: a considerable diffusion enhancement of the confined liquid water as compared with bulk water at the same thermodynamic condition, a weak adhesion energy and few H-bonds above the hydrophobic water layer, which may also sustain a water droplet. These findings introduce a new concept in water/rare-earth oxide interfaces: hydrophobicity mediated by specific water patterns on a hydrophilic surface.

Temperature-pressure phase diagram of confined monolayer water/ice at first-principles accuracy with a machine-learning force field

Understanding the phase behaviour of nanoconfined water films is of fundamental importance in broad fields of science and engineering. However, the phase behaviour of the thinnest water film - monolayer water - is still incompletely known. Here, we developed a machine-learning force field (MLFF) at first-principles accuracy to determine the phase diagram of monolayer water/ice in nanoconfinement with hydrophobic walls. We observed the spontaneous formation of two previously unreported high-density ices, namely, zigzag quasi-bilayer ice (ZZ-qBI) and branched-zigzag quasi-bilayer ice (bZZ-qBI). Unlike conventional bilayer ices, few inter-layer hydrogen bonds were observed in both quasi-bilayer ices. Notably, the bZZ-qBI entails a unique hydrogen-bonding network that consists of two distinctive types of hydrogen bonds. Moreover, we identified, for the first time, the stable region for the lowest-density [Formula: see text] monolayer ice (LD-48MI) at negative pressures (<-0.3 GPa). Overall, the MLFF enables large-scale first-principle-level molecular dynamics (MD) simulations of the spontaneous transition from the liquid water to a plethora of monolayer ices, including hexagonal, pentagonal, square, zigzag (ZZMI), and hexatic monolayer ices. These findings will enrich our understanding of the phase behaviour of the nanoconfined water/ices and provide a guide for future experimental realization of the 2D ices.