Highly confined water: two-dimensional ice, amorphous ice, and clathrate hydrates

HTR+: a novel algorithm for identifying type and polycrystal of gas hydrates

In this work, the hierarchical topology ring (HTR+) algorithm, an extension of the HTR algorithm, was developed for identifying gas hydrate types, cage structures, and grain boundaries (GBs) within polycrystalline structures. Utilizing molecular dynamics trajectories of polycrystalline hydrates, the accuracy of the HTR+ algorithm is validated in identifying sI, sII and sH hydrate types, hydrate grains, and GBs in multi-hydrate polycrystals, as well as clathrate cages at GBs. Additionally, during the hydrate nucleation and growth processes, clathrate cages, hydrate type, hydrate grains and ice structures are accurately recognized. Significantly, this algorithm demonstrates high efficiency, particularly for large hydrate systems. HTR+ algorithm emerges a powerful tool for identifying micro/mesoscopic structures of gas hydrates, enabling an in-depth understanding of the formation mechanisms and properties of gas hydrates.

Formation of compounds with diverse polyelectrolyte morphologies and nonlinear ion conductance in a two-dimensional nanofluidic channel

Aqueous electrolytes subjected to angstrom-scale confinement have recently attracted increasing interest because of their distinctive structural and transport properties, as well as their promising applicability in bioinspired nanofluidic iontronics and ion batteries. Here, we performed microsecond-scale molecular dynamics simulations, which provided evidence of nonlinear ionic conductance under an external lateral electric field due to the self-assembly of cations and anions with diverse polyelectrolyte morphologies (e.g., extremely large ion clusters) in aqueous solutions within angstrom-scale slits. Specifically, we found that the cations and anions of Li2SO4 and CaSO4 formed chain-like polyelectrolyte structures, whereas those of Na2SO4 and MgSO4 predominantly formed a monolayer of hydrated salt. Additionally, the cations and anions of K2SO4 assembled into a hexagonal anhydrous ionic crystal. These ion-dependent diverse polyelectrolyte morphologies stemmed from the enhanced Coulomb interactions, weakened hydration and steric constraints within the angstrom-scale slits. More importantly, once the monolayer hydrated salt or ionic crystal structure was formed, the field-induced ion current exhibited an intriguing gating effect at a low field strength. This abnormal ion transport was attributed to the concerted movement of cations and anions within the solid polyelectrolytes, leading to the suppression of ion currents. When the electric field exceeded a critical strength, however, the ion current surged rapidly due to the dissolution of many cations and anions within a few nanoseconds in the aqueous solution.

Formation of a two-dimensional helical square tube ice in hydrophobic nanoslit using the TIP5P water model

In extreme and nanoconfinement conditions, the tetrahedral arrangement of water molecules is challenged, resulting in a rich and new phase behavior unseen in bulk phases. The unique phase behavior of water confined in hydrophobic nanoslits has been previously observed, such as the formation of a variety of two-dimensional (2D) ices below the freezing temperature. The primary identified 2D ice phase, termed square tube ice (STI), represents a unique arrangement of water molecules in 2D ice, which can be viewed as an array of 1D ice nanotubes stacked in the direction parallel to the confinement plane. In this study, we report the molecular dynamics (MD) simulations evidence of a novel 2D ice phase, namely, helical square tube ice (H-STI). H-STI is characterized by the stacking of helical ice nanotubes in the direction parallel to the confinement plane. Its structural specificity is evident in the presence of helical square ice nanotubes, a configuration unseen in both STI and single-walled ice nanotubes. A detailed analysis of the hydrogen bonding strength showed that H-STI is a 2D ice phase diverging from the Bernal-Fowler-Pauling ice rules by forming only two strong hydrogen bonds between adjacent molecules along its helical ice chain. This arrangement of strong hydrogen bonds along ice nanotube and weak bonds between the ice nanotube shows a similarity to quasi-one-dimensional van der Waals materials. Ab initio molecular dynamics simulations (over a 30 ps) were employed to further verify H-STI's stability at 1 GPa and temperature up to 200 K.

The performance of OPC and OPC3 water models in predictions of 2D structures under nanoconfinement

Nanoconfined water plays an important role in broad fields of science and engineering. Classical molecular dynamics (MD) simulations have been widely used to investigate water phases under nanoconfinement. The key ingredient of MD is the force field. In this study, we systematically investigated the performance of a recently introduced family of globally optimal water models, OPC and OPC3, and TIP4P/2005 in describing nanoconfined two-dimensional (2D) water ice. Our studies show that the melting points of the monolayer square ice (MSI) of all three water models are higher than the melting points of the corresponding bulk ice Ih. Under the same conditions, the melting points of MSI of OPC and TIP4P/2005 are the same and are ∼90 K lower than that of the OPC3 water model. In addition, we show that OPC and TIP4P/2005 water models are able to form a bilayer AA-stacked structure and a trilayer AAA-stacked structure, which are not the cases for the OPC3 model. Considering the available experimental data and first-principles simulations, we consider the OPC water model as a potential water model for 2D water ice MD studies.

Terahertz electric field induced melting and transport of monolayer water confined in double-walled carbon nanotubes

Understanding the structures and dynamics of confined water in nanochannels holds great promise for various applications, ranging from membrane separation to blue energy collection. A setting of particular interest is the confined monolayer water within double-walled carbon nanotubes (DWCNTs), which demonstrates rich ice morphologies; however, the dynamics of this peculiar system are still unexplored. In this work, a series of molecular dynamics (MD) simulations reveal that the two-dimensional ice in DWCNTs can be effectively melted by terahertz electric fields but not by static electric fields, exhibiting an interesting ice to vapor-like transition along with extraordinary dynamical behaviors. Specifically, under appropriate field frequency, the water flow presents a sharp increase with the increase in field strength, indicating an excellent gating behavior. These remarkable findings are attributed to the resonance effect between the terahertz electric field and inherent vibration of water hydrogen bonds, causing the water molecules to change from the frozen to super permeation states. The amount of confined water exhibits a sudden reduction, confirming the breakdown of the hydrogen bond network. The distributions of density profiles, hydrogen bond number and dipole orientation demonstrate more details of the water structural change. Furthermore, under a certain field strength, the water flow shows a peculiar maximum behavior with the increase in field frequency, implying frequency optimization for water transport. These findings not only enhance our understanding of the phase transition behavior of water molecules confined within DWCNTs under the influence of a terahertz electric field but also provide a promising avenue for designing innovative nanofluidic devices.

Hydrogen Solubility in Confined Water

Confined water has demonstrated distinct structural and dynamic properties compared to bulk water. Although many studies have explored the water structure within simple geometries using materials such as carbon and silica, studies on gas solubility in confined water and the underlying physics of water structure-solubility remain limited. Recent research has illuminated the concept of "oversolubility", wherein gases display increased solubility within liquids confined in small pores compared to their bulk form. This study focuses on zeolites, naturally abundant materials with versatile applications, to study the hydrogen solubility within confined water through careful experimentation. Our findings underscore the relationship between the pore dimension and gas solubility enhancement within confined water. Hydrogen solubility is closely associated with the rearrangement of water molecules within the porous framework of the zeolite. Our research shows that a 2 nm pore size results in the greatest increase in hydrogen solubility in the water trapped inside the zeolite framework. The double donor-double acceptor (DDAA) bonds play a critical role in hydrogen solubility. Our research provides fundamental insight into the role of the molecular bonding type on hydrogen solubility in water, paving the way for potential applications in hydrogen storage and utilization.

Assisted Self-Assembly of Nanoporous Ices via Carbon Nanomaterial Templates

Self-assembly is a widely used synthetic method in nanoscience to assemble well-organized structures. Self-assembly processes usually occur in a water solvent environment. However, the self-assembly of water molecules is rarely studied. Herein, we show a strategy to fabricate porous ice via carbon nanomaterial-assisted self-assembly. Diverse frameworks of nanoporous ice are formed by using orthorhombic and tetragonal arrays of carbon nanotubes or carbon-atom chains as templates. In contrast to many bulk ices discovered in nature, nanoporous ices are shown to be stable only under negative pressure. Hence, nanoporous ices cannot be produced through the direct nucleation of water at negative pressure. The template-assisted self-assembly method is shown to be the most effective method to fabricate nanoporous ice in quantity. Several key factors for the self-assembly of nanoporous ices are identified, including proper gap spacings in the carbon nanomaterial template and suitable interactions between water and the carbon nanomaterials.

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.

Hydrogen-bonded structures and low temperature transitions of the confined water in subnano channels

The interfacial and confined water have long been attractive objects due to their crucial roles in biological, geological processes, etc. In this paper, we investigate the hydrogen-bonded structures of water and their low temperature transitions in the subnano channels of AlPO4-11 for the first time on the basis of infrared spectroscopy. The number of the adsorbed water molecules is estimated to be 8.45 per channel in one unit cell by thermogravimetric analysis. It is found that the confined water molecules are involved in saturated and unsaturated coordination with different hydrogen bond strengths at ambient temperature. The former refers to ice-like four-coordinated water and the latter includes liquid-like structures, Al-coordinated and relatively free water molecules. Unique coordination between water molecules and framework Al sites is responsible for the ice-like structures in the channels above the ice melting point. The appearance of liquid-like structures is closely related to the strong channel confinement, which does not allow the formation of extensive tetrahedral hydrogen-bonded configuration. As temperature decreases, a structural transformation of confined water happens in the channels of AlPO4-11. Isolated small water oligomers and two new components with stronger hydrogen bonds, such as low-density amorphous ice-like structures and a kind of low-density liquid-like structures are preferred. Our results provide important insights into the structural organizations and thermal-dynamic behaviors of confined water in extreme narrow channels.

Ice structures assembled from cubic water clusters of D2d and S4 symmetry

The study of self-assembly processes is of key importance for fundamental science and modern technologies. Cubic water clusters of D2d and S4 symmetry show great potential as building blocks for self-assembly. The objective of this paper is to construct possible ice structures formed by hydrogen bonding of these very stable water clusters. A number of such structures are herein presented, including quasi-2D and 3D ices as well as spatial layered and tubular ices. The energetics and structure of many configurations differing in the arrangement of hydrogen atoms in hydrogen bonds have been studied. It was established that the proton disorder of all such ices is of island type. The residual entropy of these ices is equal to ln(3)/4 in dimensionless form. For layered structures formed by the stacking of multiple bilayers, the determining role of the van der Waals interactions is shown. Note that, for all considered ices, the lowest-energy configurations are formed only by clusters of D2d symmetry.

Puckered Zigzag Monolayer Ice: Does a Confined Flat Four-Coordinated Monolayer Ice Always Have a Corresponding Puckered Phase?

We note that a flat, four-coordinated monolayer ice under confinement always has a corresponding puckered phase. Recently, a monolayer ice consisting of an array of zigzag water chains (ZZMI) predicted by first-principles calculations of water under confinement is a flat four-coordinated monolayer ice. Herein, to investigate whether puckered ZZMI exists stably, we perform molecular dynamics simulations of two-dimensional (2D) ice formation for water constrained in graphene nanocapillaries. We find a novel monolayer ice structure that can be viewed as the ZZMI puckered along the direction perpendicular to the zigzag chain (pZZMI). Unlike ZZMI that does not satisfy the ice rule, each water molecule in pZZMI can form four hydrogen bonds (HBs) via forming two stable intersublayer HBs and two intrasublayer HBs. This work provides a fresh perspective on 2D confined ice, highlighting the intrinsic connections between 2D confined ices.

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.

Evidence of Formation of 1-10 nm Diameter Ice Nanotubes in Double-Walled Carbon Nanotube Capillaries

Water exhibits rich phase behaviors in nanoscale confinement. Since the simulation evidence of the formation of single-walled ice nanotubes (INTs) in single-walled carbon nanotubes was confirmed experimentally, INTs have been recognized as a form of low-dimensional hydrogen-bonding network. However, the single-walled INTs reported in the literature all possess subnanometer diameters (<1 nm). Herein, based on systematic and large-scale molecular dynamics simulations, we demonstrate the spontaneous freezing transition of liquid water to single-walled INTs with diameters reaching ∼10 nm when confined to capillaries of double-walled carbon nanotubes (DW-CNTs). Three distinct classes of INTs are observed, namely, INTs with flat square walls (INTs-FSW), INTs with puckered rhombic walls (INTs-PRW), and INTs with bilayer hexagonal walls (INTs-BHW). Surprisingly, when water is confined in DW-CNT (3, 3)@(13, 13), an INT-FSW freezing temperature of 380 K can be reached, which is even higher than the boiling temperature of bulk water at atmospheric pressure. The freezing temperatures of INTs-FSW decrease as their caliber increases, approaching to the freezing temperature of two-dimensional flat square ice at the large-diameter limit. In contrast, the freezing temperature of INTs-PRW is insensitive to their diameter. Ab initio molecular dynamics simulations are performed to examine the stability of the INT-FSW and INT-PRW. The highly stable INTs with diameters beyond subnanometer scale can be exploited for potential applications in nanofluidic technologies and for mass transport as bioinspired nanochannels.

Nanoporous ices: an emerging class in the water/ice family

The history of scientific research on diverse ice structures dates back to more than a century. To date, 20 three-dimensional crystalline ice phases (ice I-ice XX) have been identified in the laboratory, among which ice XVI and ice XVII belong to a class of low-density nanoporous ices. Nanoporous ices can also be viewed as a special class of porous materials or water ice, as they possess a relatively high fraction of nano-cavities and/or nano-channels built into the hydrogen-bonded water framework. As such, like the prototypical class of porous materials (e.g., MOFs and COFs), nanoporous ices can be named as water oxygen-vertex frameworks (WOFs). Because of their large surface-to-volume ratio, WOFs may be potential media for gas storage, gas purification and separation. They may be applied to the biomedical field owing to their excellent biocompatibility. The field of porous ices is still emerging, as many porous ice structures that are predicted to be stable by computer simulations require future experimental confirmation. For future theoretical/computational studies, as the machine-learning method becomes an increasingly popular research tool in the material science and chemical science fields, more reliable porous ice structures and phase diagrams will be predicted with the development of more accurate machine-learning force fields.

Differences in water and vapor transport through angstrom-scale pores in atomically thin membranes

The transport of water through nanoscale capillaries/pores plays a prominent role in biology, ionic/molecular separations, water treatment and protective applications. However, the mechanisms of water and vapor transport through nanoscale confinements remain to be fully understood. Angstrom-scale pores (~2.8-6.6 Å) introduced into the atomically thin graphene lattice represent ideal model systems to probe water transport at the molecular-length scale with short pores (aspect ratio ~1-1.9) i.e., pore diameters approach the pore length (~3.4 Å) at the theoretical limit of material thickness. Here, we report on orders of magnitude differences (~80×) between transport of water vapor (~44.2-52.4 g m^-2 day^-1 Pa^-1) and liquid water (0.6-2 g m^-2 day^-1 Pa^-1) through nanopores (~2.8-6.6 Å in diameter) in monolayer graphene and rationalize this difference via a flow resistance model in which liquid water permeation occurs near the continuum regime whereas water vapor transport occurs in the free molecular flow regime. We demonstrate centimeter-scale atomically thin graphene membranes with up to an order of magnitude higher water vapor transport rate (~5.4-6.1 × 10⁴g m^-2 day^-1) than most commercially available ultra-breathable protective materials while effectively blocking even sub-nanometer (>0.66 nm) model ions/molecules.

A Clathrate Hydrate Structure Hidden in Plain Sight

Clathrate hydrates have diverse crystal structures, and among them, the three (sI, sII, and sH) most prevalent ones cover nearly all known structures, while the norm is to consider other structures only when specific guest molecules are present. Here we report the observation of a hidden clathrate structure: the tetragonal structure (TS-I) in commonly formed gas hydrates, as evidenced from molecular dynamics simulations. We show that when two (or more) sI crystal grains with different growth directions come into contact or when the growth of a sI crystal encounters geometrical frustration, the TS-I results as a cocrystal. We give evidence that TS-I may also play an important role in the combination and/or transition between sI and sII. These results imply that this previously neglected structure may be commonly present whenever sI or sII is formed. This hidden structure must be identified, experimentally and in simulations; confining the possible structures may hinder an in-depth understanding of clathrate hydrates.

The pressure induced phase diagram of double-layer ice under confinement: a first-principles study

Here, we present double-layer ice confined within various carbon nanotubes (CNTs) using state-of-the-art pressure induced (-5 GPa to 5 GPa) dispersion corrected density functional theory (DFT) calculations. We find that the double-layer ice exhibits remarkably rich and diverse phase behaviors as a function of pressure with varying CNT diameters. The lattice cohesive energies for various pure double layer ice phases follow the order of hexagonal > pentagonal > square tube > hexagonal-close-packed (HCP) > square > buckled-rhombic (b-RH). The confinement width was found to play a crucial role in the square and square tube phases in the intermediate pressure range of about 0-1 GPa. Unlike the phase transition in pure bilayer ice structures, the relative enthalpies demonstrate that the pentagonal phase, rather than the hexagonal structure, is the most stable ice polymorph at ambient pressure as well as in a deep negative pressure region, whereas the b-RH phase dominates under high pressure. The relatively short O⋯O distance of b-RH demonstrates the presence of a strong hydrogen bonding network, which is responsible for stabilizing the system.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

Two-dimensional partitioned square ice confined in graphene/graphite nanocapillaries

As one of the most fascinating confined water/ice phenomena, two-dimensional square ice has been extensively studied and experimentally confirmed in recent years. Apart from the unidirectional homogeneous square icing patterns considered in previous studies, the multidirectional partitioned square icing patterns are discovered in this study and characterized by molecular dynamics (MD) simulations. Square icing parameters are proposed to quantitatively distinguish the partitioned patterns from the homogeneous patterns and the liquid water. The number of graphene monolayers n is varied in this study, and the results show that it is more energetically favorable to form partitioned square icing patterns when the water molecules are confined between graphite sheets (n ≥ 2) compared to graphene (n = 1). This phenomenon is insensitive to n as long as n ≥ 2 because of the short-range nature of the interaction between water molecules and the carbon substrate. Moreover, it is energetically unfavorable to form partitioned square icing patterns for a single layer of water molecules even for n ≥ 2, verifying that the interaction between layers of water molecules is another dominant factor in the formation of partitioned structures. The conversion from partitioned structure to homogeneous square patterns is investigated by changing the pressure and the temperature. Based on the comprehensive MD simulations, this study unveils the formation mechanism of the partitioned square icing patterns.

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.

2D Material Nanofiltration Membranes: From Fundamental Understandings to Rational Design

Since the discovery of 2D materials, 2D material nanofiltration (NF) membranes have attracted great attention and are being developed with a tremendously fast pace, due to their energy efficiency and cost effectiveness for water purification. The most attractive aspect for 2D material NF membranes is that, anomalous water and ion permeation phenomena have been constantly observed because of the presence of the severely confined nanocapillaries (<2 nm) in the membrane, leading to its great potential in achieving superior overall performance, e.g., high water flux, high rejection rates of ions, and high resistance to swelling. Hence, fundamental understandings of such water and ion transport behaviors are of great significance for the continuous development of 2D material NF membranes. In this work, the microscopic understandings developed up to date on 2D material NF membranes regarding the abnormal transport phenomena are reviewed, including ultrafast water and ion permeation rates with the magnitude several orders higher than that predicted by conventional diffusion behavior, ion dehydration, ionic Coulomb blockade, ion-ion correlations, etc. The state-of-the-art structural designs for 2D material NF membranes are also reviewed. Discussion and future perspectives are provided highlighting the rational design of 2D material membrane structures in the future.

Self-assembled rhomboidal ammonia monolayer confined in two vertically stacked graphene oxide/graphene nanosheets

Confined water molecules have attracted widespread research interest due to their versatile phase behaviors. Ammonia (NH₃, isoelectronic with water) molecules are also expected to realize the delicate self-assembled hydrogen-bonded network like water in confinement. Here, the structures and phase behavior of NH₃ monolayers confined in two structurally symmetrical graphene oxide (GO) or graphene (G) nanosheets are investigated using first-principles calculations and ab initio molecular dynamics simulations. A highly ordered new rhomboidal phase with all NH₃ molecules adopting a Y-shaped configuration, in which one N-H bond is parallel to the confining planes and two other N-H bonds point to the top/bottom GO/G layers, respectively, was discovered at low temperature, resulting from the symmetrical confinement and subtle interlayer/intermolecular interactions. Remarkably, this new phase is so stable that a quite large strain is needed to destroy it. At room temperature, these NH₃ monolayers behave like a liquid. These rhomboidal NH₃ monolayers confined in GO/G nanosheets not only offer diverse hydrogen-bonded networks but also possess potential piezoelectricity for future device applications.

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.

Continuous Water Filling in a Graphene Nanochannel: A Molecular Dynamics Study

Low dimensional materials especially carbon materials hold high promise in the fields of water purification, mineral separation, energy harvesting/conversion, and so on. The fluidic devices fabricated by direct synthesis, lithography, or self-assembly of low dimensional materials provide opportunities for exploring the novel properties and applications of nanoconfined transport. Here, continuous filling of water and acetone molecules into a graphene nanochannel is investigated. A stairlike nonlinear dependence of the number of filling water molecules on interlayer distance d is found when d < 1 nm due to the existence of out-plane layered and in-plane ordered monolayer structure, while near-linear dependence is found for acetone because of the freely rotating configurations along with varying d during the filling process. The entropy, potential energy, and free energy of the confined system during the continuous filling are analyzed to understand the structural evolution of water. The energy-costs are discussed depending on the structure evolution of water during the filling, which is crucial to understanding the swelling and capillary condensation widely existing in the angstrom/nanometer-scale separation membranes.

Anomalously low friction of confined monolayer water with a quadrilateral structure

In this work, we explored how the structure of monolayer water confined between two graphene sheets is coupled to its dynamic behavior. Our molecular dynamics simulations show that there is a remarkable interrelation between the friction of confined water with two walls and its structure under extreme confinement. When the water molecules formed a regular quadrilateral structure, the friction coefficient is dramatically reduced. Such a low-friction coefficient can be attributed to the formation of long-range ordered hydrogen bond network, which not only decreases the structure corrugation in the direction perpendicular to the walls but also promotes the collective motion of the confined water. The regular quadrilateral structure can be formed only if the number density of confined water falls within a certain range. Higher number density results in larger structure corrugations, which increases the friction, while smaller number density leads to an irregular hydrogen bond network in which the collective motion cannot play the role. We demonstrated that there are four distinct stages in the diagram of the friction coefficient vs the number density of confined water. This research clearly established the connection between the dynamic characteristics of confined monolayer water and its structure, which is beneficial to further understand the mechanism of the high-speed water flow through graphene nanocapillaries observed in recent experiments.

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.

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.

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.

Electric field triggered release of gas from a quasi-one-dimensional hydrate in the carbon nanotube

We systematically investigate the effects of an axial electric field on the formation and decomposition of quasi-one-dimensional nitrogen gas hydrates within a single-walled carbon nanotube (SWNT) by using molecular dynamics (MD) simulations. We find that the nitrogen hydrate in the SWNT undergoes a series of structure phase transitions with increasing electric field. Corresponding to the structure transition, the nitrogen gas releases from the carbon nanotube in the electric field range of 1 V nm^-1 to 2 V nm^-1. However, nitrogen molecules are trapped as guest molecules, forming a molecule wire, in the ice nanotube when the electric field is less than 1 V nm^-1 or larger than 2 V nm^-1. Our simulations indicate that the nanotube is an excellent tiny gas tank that can be used to trap gas molecules and control their release triggered sensitively by electric signals. The key to this phenomenon is the change in orientations of water dipoles induced by the electric field, which leads to the structural change in the hydrogen-bonding network and the change in the diffusion coefficient of the water molecules. Our findings here may help understanding the mechanism of the electrorelease of gas from hydrates confined in the nanoscale space.

Electric Field Induced Dewetting of Hydrophobic Nanocavities at Ambient Temperature

The understanding of water dewetting in nanoporous materials is of great importance in various fields of science and technology. Herein, we report molecular dynamics simulation results of dewetting of water droplet in hydrophobic nanocavities between graphene walls under the influence of electric field. At ambient temperature, the rate of dewetting induced by electric field is significantly large. Whereas, it is a very low rate of dewetting induced by high temperature (423 K) due to the strong interaction of the hydrogen-bonding networks of water droplets in nanocavities. In addition, the electric filed induced formation of a water column has been found in a vacuum chamber. When the electric field is turned off, the water column will transform into a water droplet. Importantly, the results demonstrate that the rate of electric field-induced dewetting increases with growth of the electric field. Overall, our results suggest that electric field may have a great potential application for nanomaterial dewetting.

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.

Two-dimensional hydrogen hydrates: structure and stability

The structure and stability of two-dimensional hydrogen hydrate were investigated in this work using density functional theory. The results are in line with expectations that the occupied cages are more stable after their confinement between two parallel hydrophobic sheets. The four two-dimensional hydrogen hydrate crystals - BLHH-I, BLHH-II, BLHH-III and BLHH-IV - that we predicted were much more stable in a restricted environment than in a free environment, even close to or exceeding conventional hydrogen hydrates. Besides, we found that the stability of two-dimensional hydrates is inversely related to the increase in temperature. Our work highlights that two-dimensional hydrates provide a new research idea in the field of hydrogen storage.

A general topological network criterion for exploring the structure of icy nanoribbons and monolayers

We develop intuitive metrics for quantifying complex nucleating systems under confinement.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

Direct observation of 2-dimensional ices on different surfaces near room temperature without confinement

Water-solid interfaces play important roles in a wide range of fields, including atmospheric science, geochemistry, electrochemistry, and food science. Herein, we report simulation evidence of 2-dimensional (2D) ice formation on various surfaces and the dependence of the 2D crystalline structure on the hydrophobicity and morphology of the underlying surface. Contrary to the prevailing view that nanoscale confinement is necessary for the 2D liquid-to-bilayer ice transition, we find that the liquid-to-bilayer hexagonal ice (BHI) transition can occur either on a model smooth surface or on model fcc-crystal surfaces with indices of (100), (110), and (111) near room temperature. We identify a critical parameter that characterizes the water-surface interaction, above which the BHI can form on the surface. This critical parameter increases as the temperature increases. Even at temperatures above the freezing temperature of bulk ice (Ih ), we find that BHI can also form on a superhydrophilic surface due to the strong water-surface interaction. The tendency toward the formation of BHI without confinement reflects a proper water-surface interaction that can compensate for the entropy loss during the freezing transition. Furthermore, phase diagrams of 2D ice formation are described on the plane of the adsorption energy versus the fcc lattice constant (Eads-afcc), where 4 monolayer square-like ices are also identified on the fcc model surfaces with distinct water-surface interactions.

Ice Nucleation of Confined Monolayer Water Conforms to Classical Nucleation Theory

We confirmed that monolayer water confined by parallel graphene sheets spontaneously crystallizes from a structurally and dynamically heterogeneous liquid phase under moderate supercooling via direct molecular dynamics simulation. Square-lattice-like geometric order is observed at the early stage of nucleation and is preserved during the entire nucleus growth process. The diffusion coefficient and free energy profile in the cluster space extracted from a Bayesian trajectory analysis agree well with the classical nucleation theory (CNT) prediction and yield thermodynamic quantities exhibiting linear temperature dependence. The effectiveness of maximum cluster size as the descriptor of ice nucleation dynamics in the CNT framework can be attributed to the dynamical time scale decoupling and strong structural pattern dependence of density fluctuation in the liquid phase.

Assembly of peptides in mica-graphene nanocapillaries controlled by confined water

Water in nanoscale-confined geometries has unique physicochemical properties in contrast to bulk water, and is believed to play important roles in biological processes although there is less direct information available in the literature. Here, we report the self-assembly behaviors of a neurodegenerative disease related peptide termed GAV-9 encapsulated in mica-graphene nanocapillaries interacting with water nanofilms condensed under ambient conditions, based on atomic force microscopy (AFM) imaging and molecular dynamics (MD) simulations. The results revealed that, upon increase in the humidity, the GAV-9 peptide monomers adsorbed the confined water molecules and transitioned to unexpected hydrogel-like structures. Our MD simulations also suggested that in the confined mica-graphene nanocapillaries, the GAV-9 peptide monomers would indeed form water-rich hydrogel structures instead of highly ordered nanofilaments. The interfacial water confined in the mica-graphene nanocapillary is found to be crucial for such a transition. Moreover, the distribution of confined water layers largely depended on the locations of the preformed peptide nanofilaments, and the peptide nanofilaments further assembled into nanosheets with the water layer increasing, but depolymerized to amorphous peptide assemblies with the water layer decreasing. The polymerization and depolymerization of the peptide nanofilaments could be controlled in a reversible manner. Our results have supplied a simplified model system to uncover the effects of the confined interfacial water on the biological process at the molecular level.

In Situ Study of Molecular Structure of Water and Ice Entrapped in Graphene Nanovessels

Water is ubiquitous in natural systems, ranging from the vast oceans to the nanocapillaries in the earth crust or cellular organelles. In bulk or in intimate contact with solid surfaces, water molecules arrange themselves according to their hydrogen (H) bonding, which critically affects their short- and long-range molecular structures. Formation of H-bonds among water molecules designates the energy levels of certain nonbonding molecular orbitals of water, which are quantifiable by spectroscopic techniques. While the molecular architecture of water in nanoenclosures is of particular interest to both science and industry, it requires fine spectroscopic probes with nanometer spatial resolution and sub-eV energy sensitivity. Graphene liquid cells (GLCs), which feature opposing closely spaced sheets of hydrophobic graphene, facilitate high-resolution transmission electron microscopy (TEM) and electron energy-loss spectroscopy (EELS) measurements of attoliter water volumes encapsulated tightly in the GLC nanovessels. We perform in situ TEM and EELS analysis of water encased in thin GLCs exposed to room and cryogenic temperatures to examine the nanoscale arrangement of the contained water molecules. Simultaneous quantification of GLC thickness leads to the conclusion that H-bonding strengthens under increased water confinement. The present results demonstrate the feasibility of nanoscale chemical characterization of aqueous fluids trapped in GLC nanovessels and offer insights on water molecule arrangement under high-confinement conditions.

Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopores

The freezing of water mostly proceeds via heterogeneous ice nucleation, a process in which an effective nucleation medium not only expedites ice crystallization but also may effectively direct the polymorph selection of ice. Here, we show that water confined within a hydrophobic slit nanopore exhibits a freezing behavior strongly distinguished from its bulk counterpart. Such a difference is reflected by a strong, non-monotonic pore-size dependence of freezing temperature but, more surprisingly, by an unexpected stacking ordering of crystallized two-dimensional ice containing just a few ice layers. In particular, confined trilayer ice is found to exclusively crystallize into a well-ordered, hexagonal stacking sequence despite the fact that nanopore exerts no explicit constraint on stacking order. The absence of cubic stacking sequence is found to be originated from the intrinsically lower thermodynamic stability of cubic ice over hexagonal ice at the interface, which contrasts sharply the nearly degenerated stability of bulk hexagonal and cubic ices. Detailed examination clearly reveals that the divergence is attributed to the inherent difference between the two ice polymorphs in their surface phonon modes, which is further found to generically occur at both hydrophobic and hydrophilic surfaces.

Low-Dimensional Confined Ice Has the Electronic Signature of Liquid Water

Water confined in nanomaterials demonstrates anomalous behavior. Recent experiments and simulations have established that room-temperature water inside carbon nanotubes and between graphene layers behaves as solid ice: its molecules form four hydrogen bonds in a highly organized network with long-range order and exhibit low mobility. Here, we applied a first-principle energy decomposition analysis to reveal that the strength and patterns of donor-acceptor interactions between molecules in these low-dimensional ice structures resemble those in bulk liquid water rather than those in hexagonal ice. A correlation analysis shows that this phenomenon originates from a variety of hydrogen-bond distortions, different in 1D and 2D ice, from the tetrahedral configuration due to constraints imposed by nanomaterials. We discuss the implications of the reported interplay between the electronic and geometric structure of hydrogen bonds in "room-temperature ice" for computer modeling of confined water using traditional force fields.

Thick Two-Dimensional Water Film Confined between the Atomically Thin Mica Nanosheet and Hydrophilic Substrate

The interesting properties of water molecules confined in a two-dimensional (2D) environment have aroused great attention. However, the study of 2D-confined water at the hydrophilic-hydrophilic interface is largely unexplored due to the lack of appropriate system. In this work, the behavior of water molecules confined between an atomically thin mica nanosheet and a hydrophilic SiO₂/Si substrate was investigated using an atomic force microscope in detail at ambient conditions. The confined water molecules aggregated as droplets when the relative humidity (RH) of the environment was 11%. A large-area 2D water film with a uniform thickness of ∼2 nm was observed when the mica flake was incubated at 33% RH for 1 h before being mechanically exfoliated on a SiO₂/Si substrate. Interestingly, the water film showed ordered edges with a predominant angle of 120°, which was the same with the lattice orientation of the mica nanosheet on top of it. The water film showed a fluidic behavior at the early stage and reached a stable state after 48 h under ambient conditions. The surface properties of the upper mica nanosheet and the underlying substrate played a crucial role in manipulating the behavior of confined water molecules. When the surface of the upper mica nanosheet was modified by Na⁺, Ni²⁺_, and aminopropyltriethoxysilane (APS), only some small water droplets were observed instead of a water film. The surface of the underlying SiO₂/Si substrate was functionalized by hydrophilic APS and hydrophobic octadecyltrimethoxysiliane (OTS). The small water droplets were imaged on a hydrophobic OTS-SiO₂/Si substrate, while the water film with regular edges was maintained on a hydrophilic APS-SiO₂/Si substrate. Our results might provide an alternative molecular view for investigating structures and properties of confined water molecules in 2D environments.

Ice-Liquid Oscillations in Nanoconfined Water

Nanoscale confinement has a strong effect on the phase behavior of water. Studies in the last two decades have revealed a wealth of novel crystalline and quasicrystalline structures for water confined in nanoslits. Less is known, however, about the nature of ice-liquid coexistence in extremely nanoconfined systems. Here, we use molecular simulations to investigate the ice-liquid equilibrium for water confined between two nanoscopic disks. We find that the nature of ice-liquid phase coexistence in nanoconfined water is different from coexistence in both bulk water and extended nanoslits. In highly nanoconfined systems, liquid water and ice do not coexist in space because the two-phase states are unstable. The confined ice and liquid phases coexist in time, through oscillations between all-liquid and all-crystalline states. The avoidance of spatial coexistence of ice and liquid originates on the non-negligible cost of the interface between confined ice and liquid in a small system. It is the result of the small number of water molecules between the plates and has no analogue in bulk water.

Phase behaviors of deeply supercooled bilayer water unseen in bulk water

Akin to bulk water, water confined to an isolated nanoslit can show a wealth of new 2D phases of ice and amorphous ice, as well as unusual phase behavior. Indeed, 2D water phases, such as bilayer hexagonal ice and monolayer square ice, have been detected in the laboratory, confirming earlier computational predictions. Herein, we report theoretical evidence of a hitherto unreported state, namely, bilayer very low density amorphous ice (BL-VLDA), as well as evidence of a strong first-order transition between BL-VLDA and the BL amorphous ice (BL-A), and a weak first-order transition between BL-VLDA and the BL very low density liquid (BL-VLDL) water. The diffusivity of BL-VLDA is typically in the range of 10-9 cm2/s to 10-10 cm2/s. Similar to bulk (3D) water, 2D water can exhibit two forms of liquid in the deeply supercooled state. However, unlike supercooled bulk water, for which the two forms of liquid can coexist and merge into one at a critical point, the 2D BL-VLDL and BL high-density liquid (BL-HDL) phases are separated by the highly stable solid phase of BL-A whose melting line exhibits the isochore end point (IEP) near 220 K in the temperature-pressure diagram. Above the IEP temperature, BL-VLDL and BL-HDL are indistinguishable. At negative pressures, the metastable BL-VLDL exhibits a spatially and temporally heterogeneous structure induced by dynamic changes in the nanodomains, a feature much less pronounced in the BL-HDL.