A computational investigation of the phase behavior and capillary sublimation of water confined between nanoscale hydrophobic plates

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.

Connecting Non-Gaussian Water Density Fluctuations to the Lengthscale Dependent Crossover in Hydrophobic Hydration

We connect density fluctuations in liquid water to lengthscale dependent crossover in hydrophobic hydration. Specifically, we employ indirect umbrella sampling (INDUS) simulations to characterize density fluctuations in observation volumes of various sizes and shapes in water and as a function of temperature and salt concentration. Consistent with previous observations, density fluctuations are Gaussian in small molecular scale volumes, but they display non-Gaussian "low-density fat tails" in larger volumes. These non-Gaussian tails are indicative of the proximity of water to its liquid to vapor phase transition and have implications on biomolecular interactions and function. We show that the onset of non-Gaussian fluctuations in large volumes is accompanied by the formation of a cavity in the observation volume. We develop a model that uses the physics of cavity-water interface formation as a key ingredient and show that it captures the nature of non-Gaussian density fluctuations over a broad region in water and in salt solutions. We discuss the limitations of this model in the very low density region of the distribution. Our calculations provide new insights into the origins of non-Gaussian density fluctuations in water and their connections to lengthscale dependent crossover in hydrophobic hydration.

The first-principles phase diagram of monolayer nanoconfined water

Water in nanoscale cavities is ubiquitous and of central importance to everyday phenomena in geology and biology. However, the properties of nanoscale water can be substantially different from those of bulk water, as shown, for example, by the anomalously low dielectric constant of water in nanochannels¹, near frictionless water flow² or the possible existence of a square ice phase³. Such properties suggest that nanoconfined water could be engineered for technological applications in nanofluidics⁴, electrolyte materials⁵ and water desalination⁶. Unfortunately, challenges in experimentally characterizing water at the nanoscale and the high cost of first-principles simulations have prevented the molecular-level understanding required to control the behaviour of water. Here we combine a range of computational approaches to enable a first-principles-level investigation of a single layer of water within a graphene-like channel. We find that monolayer water exhibits surprisingly rich and diverse phase behaviour that is highly sensitive to temperature and the van der Waals pressure acting within the nanochannel. In addition to multiple molecular phases with melting temperatures varying non-monotonically by more than 400 kelvins with pressure, we predict a hexatic phase, which is an intermediate between a solid and a liquid, and a superionic phase with a high electrical conductivity exceeding that of battery materials. Notably, this suggests that nanoconfinement could be a promising route towards superionic behaviour under easily accessible conditions.

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.

Nanoconfined Fluids: Uniqueness of Water Compared to Other Liquids

Nanoconfinement can drastically change the behavior of liquids, puzzling us with counterintuitive properties. It is relevant in applications, including decontamination and crystallization control. However, it still lacks a systematic analysis for fluids with different bulk properties. Here we address this gap. We compare, by molecular dynamics simulations, three different liquids in a graphene slit pore: (1) A simple fluid, such as argon, described by a Lennard-Jones potential; (2) an anomalous fluid, such as a liquid metal, modeled with an isotropic core-softened potential; and (3) water, the prototypical anomalous liquid, with directional HBs. We study how the slit-pore width affects the structure, thermodynamics, and dynamics of the fluids. All the fluids show similar oscillating properties by changing the pore size. However, their free-energy minima are quite different in nature: (i) are energy-driven for the simple liquid; (ii) are entropy-driven for the isotropic core-softened potential; and (iii) have a changing nature for water. Indeed, for water, the monolayer minimum is entropy driven, at variance with the simple liquid, while the bilayer minimum is energy driven, at variance with the other anomalous liquid. Also, water has a large increase in diffusion for subnm slit pores, becoming faster than bulk. Instead, the other two fluids have diffusion oscillations much smaller than water, slowing down for decreasing slit-pore width. Our results, clarifying that water confined at the subnm scale behaves differently from other (simple or anomalous) fluids under similar confinement, are possibly relevant in nanopores applications, for example, in water purification from contaminants.

Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore

A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.

Rich topologies of monolayer ices via unconventional electrowetting

Accurate manipulation of a substance on the nanoscale and ultimately down to the level of a single atom or molecule is an ongoing subject of frontier research. Herein, we show that topologies of water monolayers on substrates, in the complete wetting condition, can be manipulated into rich forms of ordered structures via electrowetting. Notably, two new topologies of monolayer ices were identified from our molecular dynamics simulations: one stable below room temperature and the other one having the ability to be stable at room temperature. Moreover, the wettability of the substrate can be tuned from superhydrophobic to superhydrophilic by uniformly changing the charge of each atomic site of the dipole or quadrupole distributed in an orderly manner on the model substrate. At a certain threshold value of the atomic charge, water droplets on the substrate can spread out spontaneously, achieving a complete electrowetting. Importantly, unlike the conventional electrowetting, which involves application of a uniform external electric field, we proposed non-conventional electrowetting, for the first time, by invoking the electric field of dipoles and quadrupoles embedded in the substrate. Moreover, different topologies of water monolayers can be achieved by using the non-conventional electrowetting. A major advantage of the non-conventional electrowetting is that the contact-angle saturation, a long-standing and known limitation in the field of electrowetting, can be overcome by tuning uniformly the lattice atomic charge at the surface, thereby offering a new way to mitigate the contact-angle saturation for various electrowetting applications.

Nucleation of Capillary Bridges and Bubbles in Nanoconfined CO2

Using molecular simulation, we examine the capillary condensation and the capillary evaporation of CO₂ in cylindrical nanopores. More specifically, we employ the recently developed μV T-S method to determine the microscopic mechanism associated with these processes and the corresponding free energy profiles. We calculate the free energy barrier for capillary condensation and identify that the key step consists in the nucleation of a liquid bridge of a critical size. Similarly, the free energy maximum found for the capillary evaporation process is found to correspond to the nucleation of a vapor bubble of a critical size. In addition, we assess the impact of the strength of the wall-fluid on the height of the free energy barrier and on the critical size of liquid bridges (condensation process) and vapor bubbles (evaporation process). We observe that the height of the free energy barrier increases with the strength of the wall-fluid interactions. Finally, we build a theoretical model, based on capillary theory, to rationalize our findings. In particular, the simulation results reveal a linear scaling of the free energy barrier with the critical size, in excellent agreement with the theoretical predictions.

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.

Study of distance dependence of hydrophobic force between two graphene-like walls and a signature of pressure induced structure formation in the confined water

We examine the separation distance dependence of the hydrophobic force by systematically varying the distance (d) between two walls. The hydrophobic force exhibits a distance mediated crossover from a liquid-like to a gas-like behavior at around d ∼ 12 Å for 1 atm pressure. The distance dependence can be fitted to a bi-exponential form, with the longer distance part displaying a correlation length of 20 Å. In addition, the crossover is found to be accompanied by a divergent-like growth of the local relative number fluctuation of the water molecules confined between the two surfaces. Furthermore, at a fixed separation (d = 20 Å), we observe a pressure induced structural modification of confined water at high pressure. The confined water is found to form an ordered structure at high pressure (10 000 atm) and room temperature, in agreement with the experimental study [G. Algara-Siller et al. Nature 519(7544), 443 (2015)].

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.

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.

Effects of confinement between attractive and repulsive walls on the thermodynamics of an anomalous fluid

We study via molecular-dynamics simulations the thermodynamics of an anomalous fluid confined in a slit pore with one wall structured and attractive and another unstructured and repulsive. We find that the phase diagram of the homogeneous part of the confined fluid is shifted to higher temperatures, densities, and pressures with respect to the bulk, but it can be rescaled on the bulk case. We calculate a moderate increase of mobility of the homogeneous confined fluid that we interpret as a consequence of the layering due to confinement and the collective modes due to long-range correlations. We show that, as in bulk, the confined fluid has structural, diffusion, and density anomalies that order in the waterlike hierarchy, and a liquid-liquid critical point (LLCP). The overall anomalous region moves to higher temperatures, densities, and pressure, and the LLCP displaces to higher temperature compared to bulk. Motivated by experiments, we calculate also the phase diagram not just for the homogeneous part of the confined fluid but for the entire fluid in the pore, and we show that it is shifted toward higher pressures but preserves the thermodynamics, including the LLCP. Because our model has waterlike properties, we argue that in experiments with supercooled water confined in slit pores with a width of >3 nm if hydrophilic and of >1.5 nm if hydrophobic, the existence of the LLCP could be easier to test than in bulk, where it is not directly accessible.

The impacts of surface polarity on the solubility of nanoparticle

In order to study the dependence of water solubility and hydration behavior of nanoparticles on their surface polarity, we designed polar nanoparticles with varying surface polarity by assigning atomic partial charge to the surface of C60. The water solubility of the nanoparticle is enhanced by several orders of magnitude after the introduction of surface polarity. Nevertheless, when the atomic partial charge grows beyond a certain value (qM), the solubility continuously decreases to the level of nonpolar nanoparticle. It should be noted that such qM is comparable with atomic partial charge of a variety of functional groups. The hydration behaviors of nanoparticles were then studied to investigate the non-monotonic dependence of solubility on the surface polarity. The interaction between the polar nanoparticle and the hydration water is stronger than the nonpolar counterpart, which should facilitate the dissolution of the nanoparticles. On the other hand, the surface polarity also reduces the interaction of hydration water with the other water molecules and enhances the interaction between the nanoparticles which may hinder their dispersion. Besides, the introduction of surface polarity disturbs and even rearranges the hydration structure of nonpolar nanoparticle. Interestingly, the polar nanoparticle with less ordered hydration structure tends to have higher water solubility.

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.

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.

The role of material flexibility on the drying transition of water between hydrophobic objects: a thermodynamic analysis

Liquid water confined between hydrophobic objects of sufficient size becomes metastable with respect to its vapor at separations smaller than a critical drying distance. Macroscopic thermodynamic arguments predicting this distance have been restricted to the limit of perfectly rigid confining materials. However, no material is perfectly rigid and it is of interest to account for this fact in the thermodynamic analysis. We present a theory that combines the current macroscopic theory with the thermodynamics of elasticity to derive an expression for the critical drying distance for liquids confined between flexible materials. The resulting expression is the sum of the well-known drying distance for perfectly rigid confining materials and a new term that accounts for flexibility. Thermodynamic arguments show that this new term is necessarily positive, meaning that flexibility increases the critical drying distance. To study the expected magnitude and scaling behavior of the flexible term, we consider the specific case of water and present an example of drying between thin square elastic plates that are simply supported along two opposite edges and free at the remaining two. We find that the flexible term can be the same order of magnitude or greater than the rigid solution for materials of biological interest at ambient conditions. In addition, we find that when the rigid solution scales with the characteristic size of the immersed objects, the flexible term is independent of size and vice versa. Thus, the scaling behavior of the overall drying distance will depend on the relative weights of the rigid and flexible contributions.

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.

Surface effect on the electromelting behavior of nanoconfined water

Lateral diffusivity of water under confinement as a function of electric field.

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.

Sensitivity of local hydration behaviour and conformational preferences of peptides to choice of water model

Secondary structural preferences of the beta-hairpin of the 2GB1 protein in the folded and unfolded ensembles are shown to be sensitive to the choice of water model.

Ferroelectric hexagonal and rhombic monolayer ice phases

Two new monolayer ice phases are predicted from molecular dynamics simulations, both proven to be ferroelectric.

Water and water-like liquids: relationships between structure, entropy and mobility

Liquids with very diverse underlying interactions share the thermodynamic and transport anomalies of water, including metalloids, ionic melts and mesoscopic fluids. The generic feature that characterises such water-like liquids is a density-driven shift in the nature of local order in the condensed phases. The key semiquantitative relationships between structural order, thermodynamics and transport that are necessary in order to map out the consequences of this common qualitative feature for liquid-state properties and phase transformations of such systems are reviewed here. The application of these ideas to understand and model tetrahedral liquids, especially water, is discussed and possible extensions to other complex fluids are considered.

Polymorphism and polyamorphism in bilayer water confined to slit nanopore under high pressure

A distinctive physical property of bulk water is its rich solid-state phase behavior, which includes 15 crystalline (ice I-ice XIV) and at least 3 glassy forms of water, namely, low-density amorphous, high-density amorphous, and very-high-density amorphous (VHDA). Nanoscale confinement adds a new physical variable that can result in a wealth of new quasi-2D phases of ice and amorphous ice. Previous computer simulations have revealed that when water is confined between two flat hydrophobic plates about 7-9 Å apart, numerous bilayer (BL) ices (or polymorphs) can arise [e.g., BL-hexagonal ice (BL-ice I)]. Indeed, growth of the BL-ice I through vapor deposition on graphene/Pt(111) substrate has been achieved experimentally. Herein, we report computer simulation evidence of pressure-induced amorphization from BL-ice I to BL-amorphous and then to BL-VHDA(2) at 250 K and 3 GPa. In particular, BL-VHDA(2) can transform into BL-VHDA(1) via decompression from 3 to 1.5 GPa at 250 K. This phenomenon of 2D polyamorphic transition is akin to the pressure-induced amorphization in 3D ice (e.g., from hexagonal ice to HDA and then to VHDA via isobaric annealing). Moreover, when the BL-ice I is compressed instantly to 6 GPa, a new very-high-density BL ice is formed. This new phase of BL ice can be viewed as an array of square ice nanotubes. Insights obtained from pressure-induced amorphization and crystallization of confined water offer a guide with which to seek a thermodynamic path to grow a new form of methane clathrate whose BL ice framework exhibits the Archimedean 4⋅8(2) (square-octagon) pattern.