The anomalously high melting temperature of bilayer ice

Effect of confinement on water properties in super-hydrophilic pores using MD simulations with the mW model

CONTEXT

We explore the influence of strongly hydrophilic confinement on various properties of water, such as density, enthalpy, potential energy, radial distribution function, entropy, specific heat capacity, structural dynamics, and transition temperatures (freezing and melting temperatures), using monatomic water (mW) model. The properties of water are found to be dependent on confinement and the wall-fluid surface interaction. Hysteresis loops are observed for density, enthalpy, potential energy, and entropy around the transition temperatures, while the size of hysteresis loops varies with confinement and surface interaction. In smaller pore sizes (H ≤ 20), the solid phase displays a higher density compared to the liquid phase, which is unconventional behavior compared to bulk water systems due to the pronounced hydrophilic properties of the confinement surface. Specific heat capacity exhibits more oscillations in the confined system compared to bulk water, stemming from uneven enthalpy differences across equal temperature intervals. During phase transformation in both heating and quenching processes, there is an abrupt change observed in specific heat capacity. Confinement exerts a notable impact on entropy in the solid phase, but its influence is negligible in the liquid phase. At lower pore sizes (H < 25 Å), there is more fluctuation in freezing temperature for all wall-fluid interactions, which diminishes beyond pore sizes of H > 25 Å. Similarly, more oscillatory behavior is observed in melting temperatures at lower pore sizes (H < 40 Å), which diminishes at higher pore sizes (H > 40 Å). During the quenching process, a sudden jump in the in-plane orientational and tetrahedral order parameters indicates the formation of an ordered phase, specifically a diamond crystalline structure. The percentages of different crystalline structures (cubic diamond, hexagonal diamond, and 2D hexagonal) vary with both the confinement size and the wall-fluid interaction strength.

METHODS

Cooling and heating simulations are conducted with the mW water model using LAMMPS for different nanoscale confinement separation sizes ranging from 8.5 to 70 Å within the temperature range of 100-350 K. The water is modeled using two-body and three-body interaction potential (Stillinger-Weber potential) and the confinement is introduced using LJ 9-3 water-wall interaction potential. Entropy is calculated using RDF data obtained from the simulation experiments for each temperature point with increments or decrements of 2.5 K. The transition temperatures are estimated using the specific heat capacity analysis.

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

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

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.

Effect of sodium chloride adsorption on the surface premelting of ice

We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na⁺/Cl^-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na⁺/Cl^- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na⁺/Cl^- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.

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.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

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.

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

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

Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopores

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

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.

Dynamics and structural behavior of water in large confinement with planar amorphous walls

We study the structure and dynamics of liquid water confined between planar amorphous walls using molecular dynamics (MD) simulations. We report MD results for systems of more than 23 000 SPC/E water molecules confined between two hydrophilic or hydrophobic walls, separated by distances of about 15 nm. We find that the walls induce ordering of the liquid and slow down the dynamics, affecting the properties of the confined water up to distances of about 8 nm at 275 K. We quantify this influence by computing dynamic and static penetration lengths and studying their temperature dependence. Our results indicate that in the temperature range considered, hydrophobic walls perturb static properties over larger lengths compared to hydrophilic walls. We also find opposite temperature trends in the dynamic penetration lengths, with hydrophobic walls increasing their range of influence on increasing the temperature.

Ice Nucleation Efficiency of Hydroxylated Organic Surfaces Is Controlled by Their Structural Fluctuations and Mismatch to Ice

Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.

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.

What Determines the Ice Polymorph in Clouds?

Ice crystals in the atmosphere nucleate from supercooled liquid water and grow by vapor uptake. The structure of the ice polymorph grown has strong impact on the morphology and light scattering of the ice crystals, modulates the amount of water vapor in ice clouds, and can impact the molecular uptake and reactivity of atmospheric aerosols. Experiments and molecular simulations indicate that ice nucleated and grown from deeply supercooled liquid water is metastable stacking disordered ice. The ice polymorph grown from vapor has not yet been determined. Here we use large-scale molecular simulations to determine the structure of ice that grows as a result of uptake of water vapor in the temperature range relevant to cirrus and mixed-phase clouds, elucidate the molecular mechanism of the formation of ice at the vapor interface, and compute the free energy difference between cubic and hexagonal ice interfaces with vapor. We find that vapor deposition results in growth of stacking disordered ice only under conditions of extreme supersaturation, for which a nonequilibrium liquid layer completely wets the surface of ice. Such extreme conditions have been used to produce stacking disordered frost ice in experiments and may be plausible in the summer polar mesosphere. Growth of ice from vapor at moderate supersaturations in the temperature range relevant to cirrus and mixed-phase clouds, from 200 to 260 K, produces exclusively the stable hexagonal ice polymorph. Cubic ice is disfavored with respect to hexagonal ice not only by a small penalty in the bulk free energy (3.6 ± 1.5 J mol(-1) at 260 K) but also by a large free energy penalty at the ice-vapor interface (89.7 ± 12.8 J mol(-1) at 260 K). The latter originates in higher vibrational entropy of the hexagonal-terminated ice-vapor interface. We predict that the free energy penalty against the cubic ice interface should decrease strongly with temperature, resulting in some degree of stacking disorder in ice grown from vapor in the tropical tropopause layer, and in polar stratospheric and noctilucent clouds. Our findings support and explain the evolution of the morphology of ice crystals from hexagonal to trigonal symmetry with decreasing temperature, as reported by experiments and in situ measurements in clouds. We conclude that selective growth of the elusive cubic ice polymorph by manipulation of the interfacial properties can likely be achieved at the ice-liquid interface but not at the ice-vapor interface.

Surface freezing of water

Freezing, melting, evaporation and condensation of water are essential ingredients for climate and eventually life on Earth. In the present work, we show how surface freezing of supercooled water in an open container is conditioned and triggered-exclusively-by humidity in air. Additionally, a change of phase is demonstrated to be triggered on the water surface forming surface ice crystals prior to freezing of bulk. The symmetry of the surface crystal, as well as the freezing point, depend on humidity, presenting at least three different types of surface crystals. Humidity triggers surface freezing as soon as it overpasses a defined value for a given temperature, generating a plurality of nucleation nodes. An evidence of simultaneous nucleation of surface ice crystals is also provided.

Free energy contributions and structural characterization of stacking disordered ices

Crystallization of ice from deeply supercooled water and amorphous ices - a process of fundamental importance in the atmosphere, interstellar space, and cryobiology - results in stacking disordered ices with a wide range of metastabilities with respect to hexagonal ice. The structural origin of this high variability, however, has not yet been elucidated. Here we use molecular dynamics simulations with the mW water model to characterize the structure of ice freshly grown from supercooled water at temperatures from 210 to 270 K, the thermodynamics of stacking faults, line defects, and interfaces, and to elucidate the interplay between kinetics and thermodynamics in determining the structure of ice. In agreement with experiments, the ice grown in the simulations is stacking disordered with a random distribution of cubic and hexagonal layers, and a cubicity that decreases with growth temperature. The former implies that the cubicity of ice is determined by processes at the ice/liquid interface, without memory of the structure of buried ice layers. The latter indicates that the probability of building a cubic layer at the interface decreases upon approaching the melting point of ice, which we attribute to a more efficient structural equilibration of ice at the liquid interface as the driving force for growth wanes. The free energy cost for creating a pair of cubic layers in ice is 8.0 J mol(-1) in experiments, and 9.7 ± 1.9 J mol(-1) for the mW water model. This not only validates the simulations, but also indicates that dispersion in cubicity is not sufficient to explain the large energetic variability of stacking disordered ices. We compute the free energy cost of stacking disorder, line defects, and interfaces in ice and conclude that a characterization of the density of these defects is required to predict the degree of metastability and vapor pressure of atmospheric ices.

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.

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.

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.

Low-density liquid water is the mother of ice: on the relation between mesostructure, thermodynamics and ice crystallization in solutions

Predicting the temperature and extent of ice freezing in aqueous solutions is crucial for areas as diverse as cryobiology and materials design. It has long been recognized that the thermodynamics of liquid water controls the temperature and kinetics of ice crystallization. Parameterizations of the freezing temperatures in terms of the water activity of the solution have been successfully established, but the fundamental origin of the thermodynamic control of the non-equilibrium crystallization of ice has remained elusive. Here we use large-scale molecular simulations to elucidate the relationship between the structure, thermodynamics, and ice crystallization temperatures for solutions of mW water and a strongly hydrophilic solute that mimics LiCI ions. Fast cooling of solutions with up to 20 mol% ions results in the formation of nanosegregated glasses with domains of low-density amorphous ice and an ion-rich vitrified solution. Slow cooling of the mixtures results in nucleation and growth of ice within the domains of four-coordinated liquid water. The temperature of crystallization Tx coincides with the temperature of appearance of nanoscopic domains of four-coordinated liquid water in the mixture, T(L). We use the insight provided by the simulations to derive a thermodynamic expression for the crystallization temperature as a function of the water activity, T(X)(a(W)), analogous to the dependence of the melting temperature, T(m)(a(W)). The simple expression derived in this work provides a good account of the experimental freezing temperatures of water and the well-known steepest dependence of Tx on solute concentration compared to that of T(m).

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.

Coarse-Graining of TIP4P/2005, TIP4P-Ew, SPC/E, and TIP3P to Monatomic Anisotropic Water Models Using Relative Entropy Minimization

Coarse-grained models are becoming a competitive alternative for modeling processes that occur over time and length scales beyond the reach of fully atomistic molecular simulations. Ideally, coarse-grained models should not only achieve high computational efficiency but also provide accurate predictions and fundamental insight into the role of molecular interactions, the characteristic behavior, and properties of the system they model. In this work we derive a series of monatomic coarse-grained water models mX(REM) from the most popular atomistic water models X = TIP3P, SPC/E, TIP4P-Ew, and TIP4P/2005, using the relative entropy minimization (REM) method. Each coarse-grained water molecule is represented by a single particle that interacts through short-ranged anisotropic interactions that encourage the formation of "hydrogen-bonded" structures. We systematically investigate the features of the coarse-grained models in reproducing over 20 structural, dynamic, and thermodynamic properties of the reference atomistic water models-including the existence and locus of the characteristic density anomaly. The mX(REM) coarse-grained models reproduce quite faithfully the radial and angular distribution function of water, produce a temperature of maximum density (TMD), and stabilize the ice I crystal. Moreover, the ratio between the TMD and the melting temperature of the crystal in the mX(REM) models and liquid-ice equilibrium properties show reasonable agreement with the results of the corresponding atomistic models. The mX(REM) models, however, severely underestimate the cohesive energy of the condensed water phases. We investigate which specific limitations of the coarse-grained models arise from the REM methodology, from the monatomic nature of the models, and from the Stillinger-Weber interaction potential form. Our analysis indicates that a small compromise in the accuracy of structural properties can result in a significant increase of the overall accuracy and representability of the coarse-grained water models. We evaluate the accuracy of the atomistic and the monatomic anisotropic coarse-grained water models, including the mW water model, in reproducing experimental water properties. We find that mW and mTIP4P/2005(REM) score closer to experiment than widely used atomistic water models. We conclude that monatomic models of water with short-range, anisotropic "hydrogen-bonding" three-body interactions can be competitive in accuracy with fully atomistic models for the study of a wide range of properties and phenomena at less than 1/100th of the computational cost.

Heterogeneous nucleation of ice on carbon surfaces

Atmospheric aerosols can promote the heterogeneous nucleation of ice, impacting the radiative properties of clouds and Earth's climate. The experimental investigation of heterogeneous freezing of water droplets by carbonaceous particles reveals widespread ice freezing temperatures. It is not known which structural and chemical characteristics of soot account for the variability in ice nucleation efficiency. Here we use molecular dynamics simulations to investigate the nucleation of ice from liquid water in contact with graphitic surfaces. We find that atomically flat carbon surfaces promote heterogeneous nucleation of ice, while molecularly rough surfaces with the same hydrophobicity do not. Graphitic surfaces and other surfaces that promote ice nucleation induce layering in the interfacial water, suggesting that the order imposed by the surface on liquid water may play an important role in the heterogeneous nucleation mechanism. We investigate a large set of graphitic surfaces of various dimensions and radii of curvature and find that variations in nanostructures alone could account for the spread in the freezing temperatures of ice on soot in experiments. We conclude that a characterization of the nanostructure of soot is needed to predict its ice nucleation efficiency.