New Hexagonal-rhombic Trilayer Ice Structure Confined between Hydrophobic Plates

Two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement

Icing plays an important role in various physical-chemical process. Although the formation of two-dimensional ice requires nanoscale confinement, two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement remains poorly understood. Here, a critical value of a surface energy parameter is identified to characterize the liquid-solid interface interaction, above which two-dimensional and three-dimensional coexisting ice can surprisingly form on the surface. The two-dimensional ice growth mechanisms could be revealed by capturing the growth and merged of the metastable edge structures. The phase diagram about temperature and pressure vs energy parameters is predicted to distinguish liquid water, two-dimensional ice and three-dimensional ice. Furthermore, the deicing characteristics of coexisting ice demonstrate that the ice adhesion strength is linearly related to the ratio of ice-surface interaction energy to ice temperature. In addition, for gas-solid phase transition, the phase diagram about temperature and energy parameters is predicted to distinguish gas, liquid water, two-dimensional ice and three-dimensional ice. This work gives a perspective for studying the singular structure and dynamics of ice in nanoscale and provides a guide for future experimental realization of the coexisting ice.

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.

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.

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.

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.

CO Separation from H2 via Hydrate Formation in Single-Walled Carbon Nanotubes

Hydrogen is an alternative fuel without generating greenhouse gas or other harmful emissions. Industrial hydrogen production, however, always contains a small fraction of carbon monoxide (CO) (∼0.5-2%) that must be removed for use in fuel cells. Here, we present molecular dynamics simulation evidence on facile separation of CO from H₂ at ambient pressure via the formation of quasi-one-dimensional (Q1D) clathrate hydrates within single-walled carbon nanotubes (SW-CNTs). At ambient pressure, Q1D CO (or H₂) clathrates in SW-CNTs are formed spontaneously when the SW-CNTs are immersed in CO (or H₂) aqueous solution. More interestingly, for the CO/H₂ aqueous solution, highly preferential adsorption of CO over H₂ occurs within the octagonal or nonagonal ice nanotubes inside of SW-CNTs. These results suggest that the formation of Q1D hydrates within SW-CNTs can be a viable and safe method for the separation of CO from H₂, which can be exploited for hydrogen purification in fuel cells.

Buckling failure of square ice-nanotube arrays constrained in graphene nanocapillaries

Graphene confinement provides a new physical and mechanical environment with ultrahigh van der Waals pressure, resulting in new quasi-two-dimensional phases of few-layer ice. Polymorphic transition can occur in bilayer constrained water/ice system. Here, we perform a comprehensive study of the phase transition of AA-stacked bilayer water constrained within a graphene nanocapillary. The compression-limit and superheating-limit (phase) diagrams are obtained, based on the extensive molecular-dynamics simulations at numerous thermodynamic states. Liquid-to-solid, solid-to-solid, and solid-to-liquid-to-solid phase transitions are observed in the compression and superheating of bilayer water. Interestingly, there is a temperature threshold (∼275 K) in the compression-limit diagram, which indicates that the first-order and continuous-like phase transitions of bilayer water depend on the temperature. Two obviously different physical processes, compression and superheating, display similar structural evolution; that is, square ice-nanotube arrays (BL-VHDI) will bend first and then transform into bilayer triangular AA stacking ice (BL-AAI). The superheating limit of BL-VHDI exhibits local maxima, while that of BL-AAI increases monotonically. More importantly, from a mechanics point of view, we propose a novel mechanism of the transformation from BL-VHDI to BL-AAI, both for the compression and superheating limits. This structural transformation can be regarded as the "buckling failure" of the square-ice-nanotube columns, which is dominated by the lateral pressure.

Two-dimensional interlocked pentagonal bilayer ice: how do water molecules form a hydrogen bonding network?

The tradeoff between the conditions of an ideal hydrogen bonding network can serve as a generic guidance to understand the rich phase behaviors of nanoconfined water.

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

Understanding phase behavior of highly confined water, ice, amorphous ice, and clathrate hydrates (or gas hydrates), not only enriches our view of phase transitions and structures of quasi-two-dimensional (Q2D) solids not seen in the bulk phases but also has important implications for diverse phenomena at the intersection between physical chemistry, cell biology, chemical engineering, and nanoscience. Relevant examples include, among others, boundary lubrication in nanofluidic and lab-on-a-chip devices, synthesis of antifreeze proteins for ice-growth inhibition, rapid cooling of biological suspensions or quenching emulsified water under high pressure, and storage of H2 and CO2 in gas hydrates. Classical molecular simulation (MD) is an indispensable tool to explore states and properties of highly confined water and ice. It also has the advantage of precisely monitoring the time and spatial domains in the sub-picosecond and sub-nanometer scales, which are difficult to control in laboratory experiments, and yet allows relatively long simulation at the 10(2) ns time scale that is impractical with ab initio molecular dynamics simulations. In this Account, we present an overview of our MD simulation studies of the structures and phase behaviors of highly confined water, ice, amorphous ice, and clathrate, in slit graphene nanopores. We survey six crystalline phases of monolayer (ML) ice revealed from MD simulations, including one low-density, one mid-density, and four high-density ML ices. We show additional supporting evidence on the structural stabilities of the four high-density ML ices in the vacuum (without the graphene confinement), for the first time, through quantum density-functional theory optimization of their free-standing structures at zero temperature. In addition, we summarize various low-density, high-density, and very-high-density Q2D bilayer (BL) ice and amorphous ice structures revealed from MD simulations. These simulations reinforce the notion that the nanoscale confinement not only can disrupt the hydrogen bonding network in bulk water but also can allow satisfaction of the ice rule for low-density and high-density Q2D crystalline structures. Highly confined water can serve as a generic model system for understanding a variety of Q2D materials science phenomena, for example, liquid-solid, solid-solid, solid-amorphous, and amorphous-amorphous transitions in real time, as well as the Ostwald staging during these transitions. Our simulations also bring new molecular insights into the formation of gas hydrate from a gas and water mixture at low temperature.