Compression Limit of Two-Dimensional Water Constrained in Graphene Nanocapillaries

Molecular dynamics simulation-based study to analyse the properties of entrapped water between gold and graphene 2D interfaces

Heterostructures based on graphene and other 2D materials have received significant attention in recent years. However, it is challenging to fabricate them with an ultra-clean interface due to unwanted foreign molecules, which usually get introduced during their transfer to a desired substrate. Clean nanofabrication is critical for the utilization of these materials in 2D nanoelectronics devices and circuits, and therefore, it is important to understand the influence of the "non-ideal" interface. Inspired by the wet-transfer process of the CVD-grown graphene, herein, we present an atomistic simulation of the graphene-Au interface, where water molecules often get trapped during the transfer process. By using molecular dynamics (MD) simulations, we investigated the structural variations of the trapped water and the traction-separation curve derived from the graphene-Au interface at 300 K. We observed the formation of an ice-like structure with square-ice patterns when the thickness of the water film was <5 Å. This could cause undesirable strain in the graphene layer and hence affect the performance of devices developed from it. We also observed that at higher thicknesses the water film is predominantly present in the liquid state. The traction separation curve showed that the adhesion of graphene is better in the presence of an ice-like structure. This study explains the behaviour of water confined at the nanoscale region and advances our understanding of the graphene-Au interface in 2D nanoelectronics devices and circuits.

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.

Ultrahigh Lithium Selective Transport in Two-Dimensional Confined Ice

Inspired by selective ion transport in biological membrane proteins, researchers developed artificial ion channels that sieve monovalent cations, catering to the increasing lithium demand. In this work, we engineered an ion transport channel based on the confined ice within two-dimensional (2D) capillaries and found that the permselectivity of monovalent cations depends on the anisotropy of the confined ice. Particularly, the 2D confined ice showed an anomalous lithium selective transport along the (002) direction in the vermiculite capillary, with the Li+/Na+ and Li+/K+ permselectivity reaching up to 556 ± 86 and 901 ± 172, respectively, superior to most ion-selective channels. However, the 2D confined ice along the (100) direction showed less Li+ permselectivity. Additionally, the anisotropy of 2D confined ice can be tuned by adjusting the interlayer spacing. By providing insights into the ion transport in the 2D confined ice, our work may inspire more design of monovalent ion-selective channels for efficient lithium separation.

Water-induced strong isotropic MXene-bridged graphene sheets for electrochemical energy storage

Graphene and two-dimensional transition metal carbides and/or nitrides (MXenes) are important materials for making flexible energy storage devices because of their electrical and mechanical properties. It remains a challenge to assemble nanoplatelets of these materials at room temperature into in-plane isotropic, free-standing sheets. Using nanoconfined water-induced basal-plane alignment and covalent and π-π interplatelet bridging, we fabricated Ti3C2Tx MXene-bridged graphene sheets at room temperature with isotropic in-plane tensile strength of 1.87 gigapascals and moduli of 98.7 gigapascals. The in-plane room temperature electrical conductivity reached 1423 siemens per centimeter, and volumetric specific capacity reached 828 coulombs per cubic centimeter. This nanoconfined water-induced alignment likely provides an important approach for making other aligned macroscopic assemblies of two-dimensional nanoplatelets.

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.

Theoretical Prediction of the Anti-Icing Activity of Two-Dimensional Ice I

Two-dimensional (2D) ice I is atomic-level ice that is composed of two interlocked atomic layers saturated with hydrogen bonds. It has recently been experimentally observed, but its properties have yet to be clarified. Accordingly, we theoretically studied the hydrophobic properties of 2D ice I. On the contrary, a simulation of a hydrogen fluoride molecule on a 2D ice surface manifested that it destroyed the 2D ice structure and connected new hydrogen bonds with water molecules. Investigations of the interfacial effect between 2D and three-dimensional (3D) ice films indicated that the network structure of 2D ice was not destroyed by a 3D ice surface, as the former was saturated with hydrogen bonds. However, the surface of 3D ice reorganized to form as many hydrogen bonds as possible. Thus, the 2D ice film was hydrophobic and inhibited the growth of 3D ice. This shows that if 2D ice can be produced on an industrial scale, it can be used as an anti-3D-icing agent under low temperatures.

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.

Unexpected Behavior in Thermal Conductivity of Confined Monolayer Water

Monolayer water can be formed under extreme confinement and will present distinctive thermodynamic properties compared with bulk water. In this work, we perform molecular dynamics simulations to study the thermal conductivity of monolayer water confined in graphene channels, finding an unexpected way of thermal conductivity of monolayer water dependent on its number density, which has a close correlation with the structure of water. The monolayer water is in an amorphous state, and its thermal conductivity increases linearly with the area density when the water density is low at first. Then, the thermal conductivity increases as the number density of water rises, which is attributed to the formation of a crystal structure and the reduction of crystal defects as the number of water molecules increases. After reaching the zenith, the thermal conductivity decreases rapidly owing to the formation of a wrinkle structure of monolayer water with excessive water molecules, which weakens the phonon dispersion. Moreover, we further investigate the remarkable effects of the channel height on both the structure and thermal conductivity of monolayer water. In summary, this study demonstrates the close connection between the thermal conductivity of monolayer water and its structure, contributing to not only expanding the understanding of the thermodynamic property of nanoconfined water but also benefiting the engineering applications for nanofluidics.

Insight into the Structure and Dynamics of Ethanol-Water Binary Mixture Confined in Nanochannel by Mica and Graphene

Investigation of the structural properties and dynamics of fluid mixture confined in nanochannels has become an essential topic in many fields due to potential applications in nanofluidic devices and biological systems. Here, we study the ethanol-water blend confined between the mica and single or multilayer graphene for different slit pore widths, ethanol content, and temperatures. Our molecular dynamics simulation indicates that water molecules are adsorbed at the mica surface, while ethanol molecules prefer to be adsorbed near the graphene surface. We find that distinct layers of ethanol molecules form as the channel width and ethanol content in the mixture are increased. The diffusion of confined ethanol and water molecules depends on the nanopore widths, concentrations, and temperatures. Interestingly, at nanopore widths of 1.0 and 1.3 nm, the mobility of confined ethanol molecules is greater than that of water molecules for all ethanol concentrations. In contrast, at pore width of 0.7 nm, the opposite behavior is observed at lower concentrations of ethanol (x_EtOH = 0.1 and 0.3) in the mixture. Furthermore, the diffusivity of ethanol and water in the mixtures increases with increasing the temperatures. The hydrogen bond and cluster analysis imply the segregation of water molecules near the mica surface, while ethanol molecules are near the opposite pore wall (graphene).

Molecular-Level Insights into Unique Behavior of Water Molecules Confined in the Heterojunction between One- and Two-Dimensional Nanochannels

With the increasing importance of nanoconfined water in various heterostructures, it is quite essential to clarify the influence of nanoconfinement on the unique properties of water molecules in the pivotal heterojunction. In this work, we reported a series of classical molecular dynamics (MD) simulations to explore nanoconfined water in the subnanometer-sized and nanometer-sized heterostructures by adjusting one-dimensional (1-D) carbon nanotubes with different diameters and two-dimensional (2-D) graphene sheets with different interlayer distances. Our simulation results demonstrated that water molecules in the 1-D/2-D heterojunction show an obvious structural rearrangement associated with the remarkable breaking and formation of hydrogen bonds (HBs), and such rearrangements in the subnanometer-sized systems are much more pronounced than those in the nanometer-sized ones. When water molecules in the 1-D/2-D heterojunctions migrate from 2-D to 1-D confinements, the ordered multi-layer structure in the 2-D confinement are completely destroyed and then transform into different circular HB networks near the nanotube orifice for better connecting to the single-file or helical HB network in the 1-D nanotubes. Furthermore, water molecules in the 1-D/2-D heterojunctions can form stronger HBs with those water molecules further away from the 1-D confinement, leading to an asymmetrical orientational distribution near the orifice. More importantly, our comparison results revealed that the 1-D confinement plays a more important role than the 2-D confinement in determining both the structures and dynamics of water molecules in the 1-D/2-D heterojunction.

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.

A molecular simulation study into the stability of hydrated graphene nanochannels used in nanofluidics devices

Graphene-based nanochannels are a popular choice in emerging nanofluidics applications because of their tunable and nanometer-scale channels. In this work, molecular dynamics (MD) simulations were employed both to (i) assess the stability of dry and hydrated graphene nanochannels and (ii) elucidate the properties of water confined in these channels, using replica-scale models with 0.66-2.38 nm channel heights. The use of flexible nanochannel walls allows the nanochannel height to relax in response to the solvation forces arising from the confined fluid and the forces between the confining surfaces, without the need for application of arbitrarily high external pressures. Dry nanochannels were found to completely collapse if the initial nanochannel height was less than 2 nm, due to attractive van der Waals interactions between the confining graphene surfaces. However, the presence of water was found to prevent total nanochannel collapse, due to repulsive hydration forces opposing the attractive van der Waals force. For nanochannel heights less than ∼1.7 nm, the confining surfaces must be relaxed to obtain accurate hydration pressures and water diffusion coefficients, by ensuring commensurability between the number of confined water layers and the channel height. For very small (∼0.7 nm), hydrated channels a pressure of 231 MPa due to the van der Waals forces was obtained. In the same system, the confined water forms a mobile, liquid monolayer with a diffusion coefficient of 4.0 × 10^-5 cm² s^-1, much higher than bulk liquid water. Although this finding conflicts with most classical MD simulations, which predict in-plane order and arrested dynamics, it is supported by experiments and recently published first-principles MD simulations. Classical simulations can therefore be used to predict the properties of water confined in sub-nanometre graphene channels, providing sufficiently realistic molecular models and accurate intermolecular potentials are employed.

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.

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.

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.

Phase-dependent friction of nanoconfined water meniscus

A water meniscus naturally forms under ambient conditions at the point of contact between a nanoscale tip and an atomically flat substrate. Here, we study the effect of the phase state of this nanoscale meniscus-consisting of coexisting monolayer, bilayer and trilayer phase domains-on the frictional behavior during tip sliding by means of molecular dynamics simulations. While the meniscus experiences a domain-by-domain liquid-to-solid phase transition induced by lateral compression, we observe an evident transition in measured friction curves from continuous sliding to stick-slip and meanwhile a gradual increase in friction forces. Moreover, the stick-slip friction can be modulated by varying lattice orientation of the monolayer ice domain in the meniscus, choosing the sliding direction or applying in-plane strains to the substrate. Our results shed light on the rational design of high-performance micro- and nano-electromechanical systems relying on hydration lubrication.

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.

Dynamics of water confined in a graphene nanochannel: dependence of friction on graphene chirality

Numerous recent studies reveal that water transport through graphene-based nanochannels exhibits unconventional behavior. Theoretical and experimental works have reported that the dynamic behavior of the constrained water could be affected by the atomic structure, surface curvature, confinement height, pressure, and mechanical strain of the confining channel. However, few studies concern the role of graphene chirality in the dynamics of confined water. In the present study, using equilibrium molecular dynamics (EMD) simulations, we simulated the water flow through different graphene-based channels to investigate the influence of graphene chirality on the dynamic behavior of the confined water. The friction coefficient at the water/graphene interface was found to be dependent on the crystallographic orientation of the graphene wall. The results also indicated that the chirality-dependent friction behavior was affected by the confinement height of the channel. Detailed analyses on the physical origin of such a chirality effect suggested that the chirality-dependent friction was ascribed to the water-graphene interaction energy barrier variation when the water flow was driven along the armchair edge and zigzag edge. These findings are beneficial to the design of graphene-based nanofluidic devices.

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.

Graphynes for Water Desalination and Gas Separation

Selective transport of mass through membranes, so-called separation, is fundamental to many industrial applications, e.g., water desalination and gas separation. Graphynes, graphene analogs yet containing intrinsic uniformly distributed pores, are excellent candidates for highly permeable and selective membranes owing to their extreme thinness and high porosity. Graphynes exhibit computationally determined separation performance far beyond experimentally measured values of commercial state-of-the-art polyamide membranes; they also offer advantages over other atomically thin membranes like porous graphene in terms of controllability in pore geometry. Here, recent progress in proof-of-concept computational research into various graphynes for water desalination and gas separation is discussed, and their theoretically predicted outstanding permeability and selectivity are highlighted. Challenges associated with the future development of graphyne-based membranes are further analyzed, concentrating on controlled synthesis of graphyne, maintenance of high structural stability to withstand loading pressures, as well asthe demand for accurate computational characterization of separation performance. Finally, possible directions are discussed to align future efforts in order to push graphynes and other 2D material membranes toward practical separation applications.

Molecular Dynamics Simulation of the Influence of Nanoscale Structure on Water Wetting and Condensation

Recent advances in the microfabrication technology have made it possible to control surface properties at micro- and nanoscale levels. Functional surfaces drastically change wettability and condensation processes that are essential for controlling of heat transfer. However, the direct observation of condensation on micro- and nanostructure surfaces is difficult, and further understanding of the effects of the microstructure on the phase change is required. In this research, the contact angle of droplets with a wall surface and the initial condensation process were analyzed using a molecular dynamics simulation to investigate the impact of nanoscale structures and their adhesion force on condensation. The results demonstrated the dependence of the contact angle of the droplets and condensation dynamics on the wall structure and attractive force of the wall surface. Condensed water droplets were adsorbed into the nanostructures and formed a water film in case of a hydrophilic surface.

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.

Molecular Origin of Superlubricity between Graphene and a Highly Hydrophobic Surface in Water

Graphene is efficient to provide ultralow friction after the formation of an incommensurate interface but is limited to dry contact conditions and specific lattice structures. In this Letter, a new strategy is proposed to achieve the superlubricity of graphene through the creation of a sliding interface between graphene and a highly hydrophobic surface of self-assembled fluoroalkyl monolayers (SAFMs) in water. A superlow friction coefficient of μ = 0.0003 was obtained, demonstrating the extremely low shear stress between graphene and hydrophobic SAFMs in water. Molecular dynamics (MD) simulation shows that a nanometer-thick water layer is intercalated between graphene and hydrophobic SAFMs, and the weak interactions between water molecules and graphene provide a small energy barrier for water molecules sliding on graphene, which contributes to superlubricity. This finding reveals how to form a superlubricity interface by water intercalation, which has implications for minimizing the friction of layered materials and hydrophobic surfaces in water.

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.

Spontaneous directional motion of water molecules in single-walled carbon nanotubes with a stiffness gradient

Controlling water molecular motion at the nanoscale is critical for many important applications, such as water splitting to produce hydrogen and oxygen, biological and chemical cell reactions, nanofluidics, drug delivery, water treatment, etc. In this paper, we propose a new nanoscale device based on carbon nanotubes (CNTs) with a stiffness gradient to create a spontaneous directional motion of water molecules, and perform molecular dynamics simulations to analyze its transport characteristics. We find that the (6, 6) CNT possesses an optimal water transport rate. In the thinner CNTs, the water molecules are strongly confined by the CNT wall, resulting in a higher friction force; while in the thicker CNTs, the driving force is lower, and the water molecules tend to form ring-like configurations, resulting in a slower motion. For the (6, 6) CNT, water molecules tend to favor a chain-like configuration, through which the molecules are able to move synergistically along the stiffness gradient, and the transportation efficiency increases with the stiffness gradient but decreases with temperature. Both energetic and kinetic analyses are performed to elucidate this fascinating directional motion. Our work demonstrates a new strategy for controlling water molecular motion at the nanoscale without resorting to any active driving source, such as electric field, temperature or pressure difference.

Complete steric exclusion of ions and proton transport through confined monolayer water

It has long been an aspirational goal to create artificial structures that allow fast permeation of water but reject even the smallest hydrated ions, replicating the feat achieved by nature in protein channels (e.g., aquaporins). Despite recent progress in creating nanoscale pores and capillaries, these structures still remain distinctly larger than protein channels. We report capillaries made by effectively extracting one atomic plane from bulk crystals, which leaves a two-dimensional slit of a few angstroms in height. Water moves through these capillaries with little resistance, whereas no permeation could be detected even for such small ions as Na+and Cl−. Only protons (H+) can diffuse through monolayer water inside the capillaries. These observations improve our understanding of molecular transport at the atomic scale.

Why different water models predict different structures under 2D confinement

Experiments of nanoconfined water between graphene sheets at high pressure suggest that it forms a square ice structure (Algara-Siller et al., Nature, 2015, 519, 443). Molecular dynamics (MD) simulations have been used to attempt to recreate this structure, but there have been discrepancies in the structure formed by the confined water depending on the simulation set-up that was employed and particularly on the choice of water model. Here, using classical molecular dynamics simulations, we have systematically investigated the effect that three different water models (SPC/E, TIP4P/2005 and TIP5P) have on the structure of water confined between two rigid graphene sheets with a 0.9 nm separation. We show that the TIP4P/2005 and the TIP5P water models form a hexagonal AA-stacked structure, whereas the SPC/E model forms a rhombic AB-stacked structure. Our work demonstrates that the formation of these structures is driven by differences in the strength of hydrogen bonds predicted by the three water models, and that the nature of the graphene/water interaction only mildly affects the phase diagram. Considering the available experimental data and first-principle simulations we conclude that, among the models tested, the TIP4P/2005 and TIP5P force fields are for now the most reliable when simulating water under confinement. © 2018 Wiley Periodicals, Inc.

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.

Phase Diagram of Water Confined by Graphene

The behavior of water confined at the nanoscale plays a fundamental role in biological processes and technological applications, including protein folding, translocation of water across membranes, and filtration and desalination. Remarkably, nanoscale confinement drastically alters the properties of water. Using molecular dynamics simulations, we determine the phase diagram of water confined by graphene sheets in slab geometry, at T = 300 K and for a wide range of pressures. We find that, depending on the confining dimension D and density σ, water can exist in liquid and vapor phases, or crystallize into monolayer and bilayer square ices, as observed in experiments. Interestingly, depending on D and σ, the crystal-liquid transformation can be a first-order phase transition, or smooth, reminiscent of a supercritical liquid-gas transformation. We also focus on the limit of stability of the liquid relative to the vapor and obtain the cavitation pressure perpendicular to the graphene sheets. Perpendicular cavitation pressure varies non-monotonically with increasing D and exhibits a maximum at D ≈ 0.90 nm (equivalent to three water layers). The effect of nanoconfinement on the cavitation pressure can have an impact on water transport in technological and biological systems. Our study emphasizes the rich and apparently unpredictable behavior of nanoconfined water, which is complex even for graphene.

Coexistence of Multilayered Phases of Confined Water: The Importance of Flexible Confining Surfaces

Flexible nanoscale confinement is critical to understanding the role that bending fluctuations play on biological processes where soft interfaces are ubiquitous or to exploit confinement effects in engineered systems where inherently flexible 2D materials are pervasively employed. Here, using molecular dynamics simulations, we compare the phase behavior of water confined between flexible and rigid graphene sheets as a function of the in-plane density, ρ_2D. We find that both cases show commensurate mono-, bi-, and trilayered states; however, the water phase in those states and the transitions between them are qualitatively different for the rigid and flexible cases. The rigid systems exhibit discontinuous transitions between an (n)-layer and an (n+1)-layer state at particular values of ρ_2D, whereas under flexible confinement, the graphene sheets bend to accommodate an (n)-layer and an (n+1)-layer state coexisting in equilibrium at the same density. We show that the flexible walls introduce a very different sequence of ice phases and their phase coexistence with vapor and liquid phases than that observed with rigid walls. We discuss the applicability of these results to real experimental systems to shed light on the role of flexible confinement and its interplay with commensurability effects.

Ultra-fast vapor generation by a graphene nano-ratchet: a theoretical and simulation study

Vapor generation is of prime importance for a broad range of applications: domestic water heating, desalination and wastewater treatment, etc. However, slow and inefficient evaporation limits its development. In this study, a nano-ratchet, a multilayer graphene with cone-shaped nanopores (MGCN), to accelerate vapor generation has been proposed. By performing molecular dynamics simulation, we found that air molecules were spontaneously transported across MGCN and resulted in a remarkable pressure difference, 21 kPa, between the two sides of MGCN. We studied the dependence of the pressure difference on the ambient temperature and geometry of MGCN in detail. Through further analysis of the diffusive transport, we found that pressure difference depended on the competition between ratchet transport and Knudsen diffusion and it was further found that ratchet transport is dominant. The significant pressure difference could lead to a 15-fold or greater enhancement of vapor generation, which shows the wide applications of this nano-ratchet.

Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.

Nanoscale Insight into the Mechanism of a Highly Oriented Pyrolytic Graphite Edge Surface Wetting by "Interferencing" Water

The new molecular dynamics simulation results showing the influence of the edge carbon surface atoms on the wettability of a highly oriented pyrolytic graphite (HOPG) surface with water nanodroplets are reported. The conditions for the occurrence of the Wenzel effect are discussed, and the Cassie-to-Wenzel transition (CTWT) mechanism in the nanoscale is explored. This transition is detected by the application of a new procedure showing that the CTWT point shifts toward larger values of carbon-oxygen potential well depth with the decrease in the HOPG side angle. It is concluded that the Wenzel effect significantly contributes to the contact angles (CAs) measured for the HOPG surfaces. The Wenzel effect is also very important for the "HOPG" structures possessing the disturbed C-C interlayer distance, and its influence on the water nanodroplet CAs is strongly pronounced. The structure of water confined inside slits and on a HOPG surface is studied using the analysis of the density profiles, the number of hydrogen bonds, and, modified for the purpose of this study, structure factor. The detailed analysis of all parameters describing confined water leads to the conclusion about the presence of characteristic interference patterns revealed as a result of long-term simulation. A simple model describing this effect is proposed as the starting point for further considerations.

Structures and thermodynamics of water encapsulated by graphene

Understanding phase behaviors of nanoconfined water has driven notable research interests recently. In this work, we examine water encapsulated under a graphene cover that offers an ideal testbed to explore its molecular structures and thermodynamics. We find layered water structures for up to ~1000 trapped water molecules, which is stabilized by the spatial confinement and pressure induced by interfacial adhesion. For monolayer encapsulations, we identify representative two-dimensional crystalline lattices as well as defects therein. Free energy analysis shows that the structural orders with low entropy are compensated by high formation energies due to the pressurized confinement. There exists an order-to-disorder transition for this condensed phase at ~480-490 K, with a sharp reduction in the number of hydrogen bonds and increase in the entropy. Fast diffusion of the encapsulated water demonstrates anomalous temperature dependence, indicating the solid-to-fluid nature of this structural transition. These findings offer fundamental understandings of the encapsulated water that can be used as a pressurized cell with trapped molecular species, and provide guidance for practical applications with its presence, for example, in the design of nanodevices and nanoconfined reactive cells.

Structure and dynamics of water confined in a graphene nanochannel under gigapascal high pressure: dependence of friction on pressure and confinement

The friction coefficient at the water/graphene interface is dependent on the lateral pressure and nanochannel height under gigapascal high-pressure.

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.

Transformation between divacancy defects induced by an energy pulse in graphene

The mutual transformations among the four typical divacancy defects induced by a high-energy pulse were studied via molecular dynamics simulation. Our study revealed all six possible mutual transformations and found that defects transformed by absorbing energy to overcome the energy barrier with bonding, debonding, and bond rotations. The reversibility of defect transformations was also investigated by potential energy analysis. The energy difference was found to greatly influence the transformation reversibility. The direct transformation path was irreversible if the energy difference was too large. We also studied the correlation between the transformation probability and the input energy. It was found that the transformation probability had a local maxima at an optimal input energy. The introduction of defects and their structural evolutions are important for tailoring the exceptional properties and thereby performances of graphene-based devices, such as nanoporous membranes for the filtration and desalination of water.

AB-stacked square-like bilayer ice in graphene nanocapillaries

Water, when constrained between two graphene sheets and under ultrahigh pressure, can manifest dramatic differences from its bulk counterparts such as the van der Waals pressure induced water-to-ice transformation, known as the metastability limit of two-dimensional (2D) liquid.

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.