Phase diagram of water between hydrophobic surfaces

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.

How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

Germination of bean seeds (Vigna unguiculata L. Walp.) in strong electric fields

This paper presents a tutorial for the germination of bean seeds (Vigna unguiculata L. Walp.) in strong electrostatic fields up to 1240 V/cm. The seeds were allowed to germinate under different electric field strengths for 48 h. Although most of such germination experiments did not show any visible effect, the field strength of 945 V/cm strongly increased the seedling's vigor during the early growth stage. In the end, 30 % more yield was obtained from stimulated seeds when compared to the control group. This article postulates for the first time a hypothesis of the mechanism of action of the electric field during germination. In biological cells of any species, water confined between narrow surfaces can undergo a phase transition that shifts its melting point to higher temperatures when an external electric field is applied. This effect has already been known as electrofreezing, and has been confirmed by several experimental and molecular modeling studies. As a consequence, the transport kinetics of molecules across cell organelle membranes might be altered, which in turn leads to different plant properties. With emphasis on the presented method, this work reports: •An inexpensive electric circuit for the generation of strong electric fields•Instructions regarding the setup and operation of an adequate germination chamber.

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.

Generalization of Young-Laplace, Kelvin, and Gibbs-Thomson equations for arbitrarily curved surfaces

The Young-Laplace, Kelvin, and Gibbs-Thomson equations form a cornerstone of colloidal and surface sciences and have found successful applications in many subfields of physics, chemistry, and biology. The Gibbs-Thomson effect, for example, predicts that small crystals are in equilibrium with their liquid melt at a lower temperature than large crystals and the positive interfacial energy increases the energy required to form small particles with a high curvature interface. In cases of liquids contained within porous media (confined geometry), the effect indicates decreasing the freezing/melting temperatures and the increment of the temperature is inversely proportional to the pore size. These phenomena can be reformulated for Gaussian maps of macromolecules and can be asked the following question: can one use the equations for predicting the melting temperature and shape of polymer chains in confined geometries? The answer is no, mainly because macromolecules form highly curved surfaces (Gaussian maps), and the equations hold only for simple geometries (sphere, plane, or cylinder). Here, we present general Young-Laplace, Kelvin, and Gibbs-Thomson equations for arbitrarily curved surfaces and apply them to predict temperature distribution on a few protein surfaces. Also, after increased interest toward liquid/liquid phase separation in biology, we derive generic Ostwald ripening and show that for shape-changing condensates, instead of a monotonic growing mechanism, a variety of processes are possible. Due to the generality of equations, we clarify that at appropriate internal/external pressure conditions systems, bounded by surfaces, may adopt any shape and thermal stability is strongly influenced by the geometries of confined spaces.

Adjustable high-speed and directional diffusion of water nanodroplets confined by graphene sheets

Diffusion of confined water is important in nanofluidic and other water transport systems. In this study, the diffusion of water nanodroplets confined by graphene sheets is investigated based on molecular dynamics simulations. We find that the confined water nanodroplets can achieve a high-speed and directional motion. The impact of the size of water nanodroplets and distance of graphene sheets on diffusion is studied. The results show that the diffusion of confined water nanodroplets is adjustable and the speed is about 3 orders of magnitude faster than that of the self-diffusing water molecules in liquid water. Subsequently, the most suitable morphology of confined nanodroplets for rapid movement is found. We also find that the direction of diffusion of confined water nanodroplets is affected by the thermal vibrations of carbon atoms. Finally, the interaction energy and friction coefficient between confined nanodroplets and graphene sheets are analyzed to give an insight into the fast and directional diffusion behaviors of water nanodroplets. Our results reveal that a variation in the structure of interfacial water molecules with the distance of graphene sheets is the key to the rapid movement of confined water nanodroplets. The phenomena reported here can enrich the knowledge of molecular mechanisms for nanoconfined water systems, and may stimulate more ideas for the rapid removal of confined water.

Hydrogel Nanosheets Confined 2D Rhombic Ice: A New Platform Enhancing Chondrogenesis

Nanoconfinement within flexible interfaces is a key step towards exploiting confinement effects in several biological and technological systems wherein flexible 2D materials are frequently utilized but are arduous to prepare. Hitherto unreported, the synthesis of 2D Hydrogel nanosheets (HNS) using a template- and catalyst-free process is developed representing a fertile ground for fundamental structure-property investigations. In due course of time, nucleating folds propagating along the edges trigger co-operative deformations of HNS generating regions of nanoconfinement within trapped water islands. These severely constricting surfaces force water molecules to pack within the nanoscale regime of HNS almost parallel to the surface bringing about phase transition into puckered rhombic ice with AA and AB Bernal stacking pattern, which was mostly restricted to Molecular dynamics (MD) studies so far. Interestingly, under high lateral pressure and spatial inhomogeneity within nanoscale confinement, bilayer rhombic ice structures were formed with an in-plane lattice spacing of 0.31 nm. In this work, a systematic exploration of rhombic ice formation within HNS has been delineated using High-resolution transmission electron microscopy (HRTEM), and its ultrathin morphology was examined using Atomic Force Microscopy (AFM). Scanning Electron Microscopy (SEM) images revealed high porosity while mechanical testing presented young's modulus of 155 kPa with ~84% deformation, whereas contact angle suggested high hydrophilicity. The combinations of nanosheets, porosity, nanoconfinement, hydrophilicity, and mechanical strength, motivated us to explore their application as a scaffold for cartilage regeneration, by inducing chondrogenesis of human Wharton Jelly derived mesenchymal stem cells (hWJ MSCs). HNS promoted the formation of cell aggregates giving higher number of spheroid formation and a marked expression of chondrogenic markers (ColI, ColII, ColX, ACAN and S-100), thereby providing some cues for guiding chondrogenic differentiation.

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.

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.

Critical In-Plane Density of Polyelectrolyte Brush for the Ordered Hydrogen-Bonded Structure of Incorporated Water

A polymer electrolyte brush is a reasonable platform to confine water molecules within a nanoscopic area to study their role in the function of interacting media because of their adjustable nanospace and charge by changing the in-plane density and side chains of the brush. Here, we demonstrate how the in-plane spacing of cationic polymer brush chains, poly[2-(methacryloyloxy)ethyltrimethylammonium chloride] (PMTAC), affects the hydrogen bond configuration of incorporated water using soft X-ray emission spectroscopy. At the critical in-plane density σ = 0.30 chains/nm² of PMTAC, tetrahedrally coordinated water molecules started to melt into distorted or broken hydrogen-bonded configurations. Considering the charge on the quaternary ammonium cations, the electric field required to form a tetrahedrally coordinated hydrogen-bonded configuration was estimated as ∼500 kV cm^-1 and is effective up to ∼1 nm from the surface of the polymer chain. These findings are useful for designing specific interface properties and the resultant surface function of polyelectrolyte-based materials.

Supercooled nano-droplets of water confined in hydrophobic rubber

Hydrophobic elastomers are capable of absorbing a small amount of water that forms droplets around hydrophilic sites. These systems allow the study of confinement effects by a hydrophobic environment on the dynamics and thermodynamic behaviour of water molecules. The freezing-melting properties and the dynamics of water inside nano-droplets in butyl rubber are affected, as revealed by differential scanning calorimetry (DSC) and deuterium nuclear magnetic resonance (²H-NMR). Upon cooling down, all water crystalizes with a bimodal droplet population (d_a = 3.4 nm and d_b = 4.4 nm) in a temperature range associated with the droplet size distribution. However, the melting temperature is not shifted according to the Gibbs-Thomson equation. The relative decrease of the ²H-NMR longitudinal magnetization is not a single exponential and, by inverse Laplace transformation, it was deduced to be bimodal in agreement with the DSC measurements (T_1,a ∼ 10 ms and T_1,b ∼ 200 ms). The deduced correlation time of molecular reorientation is longer than that of bulk water and the behaviour with temperature follows the Vogel-Fulcher-Tammann (VFT) equations with a changing fragility as the droplet size is reduced when reducing hydration.

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.

On the Gibbs-Thomson equation for the crystallization of confined fluids

The Gibbs-Thomson (GT) equation describes the shift of the crystallization temperature for a confined fluid with respect to the bulk as a function of pore size. While this century old relation is successfully used to analyze experiments, its derivations found in the literature often rely on nucleation theory arguments (i.e., kinetics instead of thermodynamics) or fail to state their assumptions, therefore leading to similar but different expressions. Here, we revisit the derivation of the GT equation to clarify the system definition, corresponding thermodynamic ensemble, and assumptions made along the way. We also discuss the role of the thermodynamic conditions in the external reservoir on the final result. We then turn to numerical simulations of a model system to compute independently the various terms entering in the GT equation and compare the predictions of the latter with the melting temperatures determined under confinement by means of hyper-parallel tempering grand canonical Monte Carlo simulations. We highlight some difficulties related to the sampling of crystallization under confinement in simulations. Overall, despite its limitations, the GT equation may provide an interesting alternative route to predict the melting temperature in large pores using molecular simulations to evaluate the relevant quantities entering in this equation. This approach could, for example, be used to investigate the nanoscale capillary freezing of ionic liquids recently observed experimentally between the tip of an atomic force microscope and a substrate.

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.

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.

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.

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.

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.

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 transition-like behavior of the water monolayer close to the polarized surface of a nanotube

By molecular dynamics simulations, we have investigated effects of temperature on the dynamical behavior of water layers at the charged surface of a nanotube. The behavior of the first water monolayer at the charged surface is very different from that of bulk water. There are three different temperature regions for the axial diffusion coefficient and they increase in different ways (linearly or exponentially) with temperature. The dipole distribution of water molecules was chosen as the order parameter to analyze the phase transition-like behavior. The simulation results indicate that the transition from ordered water to disordered water is continuous, which has not been found in the bulk counterpart. The mechanism behind the unexpected phenomenon was also investigated.

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.

Thermodynamic, structural, and dynamical properties of nano-confined water using SPC/E and TIP4P models by molecular dynamics simulations

Different morphologies of water molecules are confined between two parallel graphene surfaces.

Dynamics of Water Monolayers Confined by Chemically Heterogeneous Surfaces: Observation of Surface-Induced Anisotropic Diffusion

Water present in confining geometries plays key roles in many systems of scientific and technological relevance. Prominent examples are living cells and nanofluidic devices. Despite its importance, a complete understanding of the dynamics of water in nanoscale confinement remains elusive. In this work, we use molecular dynamics (MD) simulation to investigate the diffusive dynamics of water monolayers confined in chemically heterogeneous silica slit pores. The effect of chemical heterogeneity is systematically investigated through the fraction f_SiOH of randomly distributed surface sites that possess hydroxyl functional groups. Partial hydroxylation results in heterogeneous surfaces comprising nanoscale hydrophobic and hydrophilic regions. We find that the in-plane diffusivity of water increases monotonically with f_SiOH; at low surface hydroxylation (f_SiOH ≤ 50%), slow water dynamics arise due to the formation of icelike structures in the hydrophobic regions, while at f_SiOH ≥ 75%, surface-water H-bonds in the hydrophilic regions result in faster dynamics. We show that surface patterning with ordered hydrophobic and hydrophilic "stripes" can be used to induce one-dimensional diffusion, with water diffusing through the slit pore preferentially along the direction of the hydrophilic surface patterns.

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.

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.

Enhanced Configurational Entropy in High-Density Nanoconfined Bilayer Ice

A novel kind of crystal order in high-density nanoconfined bilayer ice is proposed from molecular dynamics and density-functional theory simulations. A first-order transition is observed between a low-temperature proton-ordered solid and a high-temperature proton-disordered solid. The latter is shown to possess crystalline order for the oxygen positions, arranged on a close-packed triangular lattice with AA stacking. Uniquely among the ice phases, the triangular bilayer is characterized by two levels of disorder (for the bonding network and for the protons) which results in a configurational entropy twice that of bulk ice.

Two Dimensional Ice from First Principles: Structures and Phase Transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2  GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Structural and configurational properties of nanoconfined monolayer ice from first principles

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

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.

Compression Limit of Two-Dimensional Water Constrained in Graphene Nanocapillaries

Evaluation of the tensile/compression limit of a solid under conditions of tension or compression is often performed to provide mechanical properties that are critical for structure design and assessment. Algara-Siller et al. recently demonstrated that when water is constrained between two sheets of graphene, it becomes a two-dimensional (2D) liquid and then is turned into an intriguing monolayer solid with a square pattern under high lateral pressure [ Nature , 2015 , 519 , 443 - 445 ]. From a mechanics point of view, this liquid-to-solid transformation characterizes the compression limit (or metastability limit) of the 2D monolayer water. Here, we perform a simulation study of the compression limit of 2D monolayer, bilayer, and trilayer water constrained in graphene nanocapillaries. At 300 K, a myriad of 2D ice polymorphs (both crystalline-like and amorphous) are formed from the liquid water at different widths of the nanocapillaries, ranging from 6.0 to11.6 Å. For monolayer water, the compression limit is typically a few hundred MPa, while for the bilayer and trilayer water, the compression limit is 1.5 GPa or higher, reflecting the ultrahigh van der Waals pressure within the graphene nanocapillaries. The compression-limit (phase) diagram is obtained at the nanocapillary width versus pressure (h-P) plane, based on the comprehensive molecular dynamics simulations at numerous thermodynamic states as well as on the Clapeyron equation. Interestingly, the compression-limit curves exhibit multiple local minima.

Water in Inhomogeneous Nanoconfinement: Coexistence of Multilayered Liquid and Transition to Ice Nanoribbons

Phase behavior and the associated phase transition of water within inhomogeneous nanoconfinement are investigated using molecular dynamics simulations. The nanoconfinement is constructed by a flat bottom plate and a convex top plate. At 300 K, the confined water can be viewed as a coexistence of monolayer, bilayer, and trilayer liquid domains to accommodate the inhomogeneous confinement. With increasing liquid density, the confined water with uneven layers transforms separately into two-dimensional ice crystals with unchanged layer number and rhombic in-plane symmetry for oxygen atoms. The monolayer water undergoes the transition first into a puckered ice nanoribbon, and the bilayer water transforms into a rhombic ice nanoribbon next, followed by the transition of trilayer water into a trilayer ice nanoribbon. The sequential localized liquid-to-solid transition within the inhomogeneous confinement can also be achieved by gradually decreasing the temperature at low liquid densities. These findings of phase behaviors of water under the inhomogeneous nanoconfinement not only extend the phase diagram of confined water but also have implications for realistic nanofluidic systems and microporous materials.

Liquid-solid and solid-solid phase transition of monolayer water: high-density rhombic monolayer ice

Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.

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.

Effect of confinement on the solid-liquid coexistence of Lennard-Jones fluid

The solid-liquid coexistence of a Lennard-Jones fluid confined in slit pores of variable pore size, H, is studied using molecular dynamics simulations. Three-stage pseudo-supercritical transformation path of Grochola [J. Chem. Phys. 120(5), 2122 (2004)] and multiple histogram reweighting are employed for the confined system, for various pore sizes ranging from 20 to 5 molecular diameters, to compute the solid-liquid coexistence. The Gibbs free energy difference is evaluated using thermodynamic integration method by connecting solid-liquid phases under confinement via one or more intermediate states without any first order phase transition among them. Thermodynamic melting temperature is found to oscillate with wall separation, which is in agreement with the behavior seen for kinetic melting temperature evaluated in an earlier study. However, thermodynamic melting temperature for almost all wall separations is higher than the bulk case, which is contrary to the behavior seen for the kinetic melting temperature. The oscillation founds to decay at around H = 12, and beyond that pore size dependency of the shift in melting point is well represented by the Gibbs-Thompson equation.