Phase transitions induced by nanoconfinement in liquid water

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.

Complex dynamics of partially freezable confined water revealed by combined experimental and computational studies

We investigate water dynamics in mesoporous silica across partial crystallization by combining broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR), and molecular dynamics simulations (MDS). Exploiting the fact that not only BDS but also NMR field-cycling relaxometry and stimulated-echo experiments provide access to dynamical susceptibilities in broad frequency and temperature ranges, we study both the fully liquid state above the melting point Tm and the dynamics of coexisting water and ice phases below this temperature. It is found that partial crystallization leads to a change in the temperature dependence of rotational correlation times τ, which occurs in addition to previously reported dynamical crossovers of confined water and depends on the pore diameter. Furthermore, we observe that dynamical susceptibilities of water are strongly asymmetric in the fully liquid state, whereas they are much broader and nearly symmetric in the partially frozen state. Finally, water in the nonfreezable interfacial layer below Tm does not exhibit a much debated dynamical crossover at ∼220 K. We argue that its dynamics is governed by a static energy landscape, which results from the interaction with the bordering silica and ice surfaces and features a Gaussian-like barrier distribution. Consistently, our MDS analysis of the motional mechanism reveals a hopping motion of water in thin interfacial layers. The rotational correlation times of the confined ice phases follow Arrhenius laws. While the values of τ depend on the pore diameter, freezable water in various types of confinements and mixtures shows similar activation energies of Ea ≈ 0.43 eV.

Probing structural superlubricity of two-dimensional water transport with atomic resolution

Low-dimensional water transport can be drastically enhanced under atomic-scale confinement. However, its microscopic origin is still under debate. In this work, we directly imaged the atomic structure and transport of two-dimensional water islands on graphene and hexagonal boron nitride surfaces using qPlus-based atomic force microscopy. The lattice of the water island was incommensurate with the graphene surface but commensurate with the boron nitride surface owing to different surface electrostatics. The area-normalized static friction on the graphene diminished as the island area was increased by a power of ~-0.58, suggesting superlubricity behavior. By contrast, the friction on the boron nitride appeared insensitive to the area. Molecular dynamic simulations further showed that the friction coefficient of the water islands on the graphene could reduce to <0.01.

Estimating Position-Dependent and Anisotropic Diffusivity Tensors from Molecular Dynamics Trajectories: Existing Methods and Future Outlook

Confinement can substantially alter the physicochemical properties of materials by breaking translational isotropy and rendering all physical properties position-dependent. Molecular dynamics (MD) simulations have proven instrumental in characterizing such spatial heterogeneities and probing the impact of confinement on materials' properties. For static properties, this is a straightforward task and can be achieved via simple spatial binning. Such an approach, however, cannot be readily applied to transport coefficients due to lack of natural extensions of autocorrelations used for their calculation in the bulk. The prime example of this challenge is diffusivity, which, in the bulk, can be readily estimated from the particles' mobility statistics, which satisfy the Fokker-Planck equation. Under confinement, however, such statistics will follow the Smoluchowski equation, which lacks a closed-form analytical solution. This brief review explores the rich history of estimating profiles of the diffusivity tensor from MD simulations and discusses various approximate methods and algorithms developed for this purpose. Besides discussing heuristic extensions of bulk methods, we overview more rigorous algorithms, including kernel-based methods, Bayesian approaches, and operator discretization techniques. Additionally, we outline methods based on applying biasing potentials or imposing constraints on tracer particles. Finally, we discuss approaches that estimate diffusivity from mean first passage time or committor probability profiles, a conceptual framework originally developed in the context of collective variable spaces describing rare events in computational chemistry and biology. In summary, this paper offers a concise survey of diverse approaches for estimating diffusivity from MD trajectories, highlighting challenges and opportunities in this area.

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.

Accurate prediction of solvent flux in sub-1-nm slit-pore nanosheet membranes

Nanosheet-based membranes have shown enormous potential for energy-efficient molecular transport and separation applications, but designing these membranes for specific separations remains a great challenge due to the lack of good understanding of fluid transport mechanisms in complex nanochannels. We synthesized reduced MXene/graphene hetero-channel membranes with sub-1-nm pores for experimental measurements and theoretical modeling of their structures and fluid transport rates. Our experiments showed that upon complete rejection of salt and organic dyes, these membranes with subnanometer channels exhibit remarkably high solvent fluxes, and their solvent transport behavior is very different from their homo-structured counterparts. We proposed a subcontinuum flow model that enables accurate prediction of solvent flux in sub-1-nm slit-pore membranes by building a direct relationship between the solvent molecule-channel wall interaction and flux from the confined physical properties of a liquid and the structural parameters of the membranes. This work provides a basis for the rational design of nanosheet-based membranes for advanced separation and emerging nanofluidics.

HBCalculator: A Tool for Hydrogen Bond Distribution Calculations in Molecular Dynamics Simulations

Hydrogen bonds, crucial noncovalent interactions in molecular systems, significantly impact biological, chemical, and energy-related processes; therefore, characterizing hydrogen bond information is of importance to fundamental studies. This work introduces the HBCalculator, a Tcl-based tool integrated with VMD for calculating 1D and 2D distributions of hydrogen bond density and strength. The tool facilitates spatial analysis, overcoming limitations in existing packages that lack direct spatial distribution output. By employing HBCalculator in MD simulations, three systems of cellulose/water and graphene/water interfaces, were tested to showcase its functionality. The 1D and 2D hydrogen bond distributions reveal insights into interfacial properties, reflecting the impact of material hydrophobicity. The simplicity of usage, along with its potential for diverse molecular systems, positions HBCalculator as a valuable tool for researchers exploring hydrogen bond networks.

Ice-Water Equilibrium in Nanoscale Confinement

We show that 2D ^{2}H NMR spectra enable valuable insights into the nature of an ice-water equilibrium in nanoscale confinement, which extends over a broad temperature range. In particular, 2D ^{2}H NMR line-shape analysis allows us to determine the timescale on which the coexisting ice and water phases exchange molecules. For D_{2}O in a silica nanopore with a diameter of 5.4 nm, we find that the residence time of a water molecule in either phase is characterized by an NMR exchange time of τ_{X}=5.7  ms at 220 K. Thus, the ice-water equilibrium is highly dynamic, which is an important aspect for an understanding of deeply cooled confined and, possibly, bulk waters.

Surface polarization enhances ionic transport and correlations in electrolyte solutions nanoconfined by conductors

Layered materials that perform mixed electron and ion transport are promising for energy harvesting, water desalination, and bioinspired functionalities. These functionalities depend on the interaction between ionic and electronic charges on the surface of materials. Here we investigate ion transport by an external electric field in an electrolyte solution confined in slit-like channels formed by two surfaces separated by distances that fit only a few water layers. We study different electrolyte solutions containing monovalent, divalent, and trivalent cations, and we consider walls made of non-polarizable surfaces and conductors. We show that considering the surface polarization of the confining surfaces can result in a significant increase in ionic conduction. The ionic conductivity is increased because the conductors' screening of electrostatic interactions enhances ionic correlations, leading to faster collective transport within the slit. While important, the change in water's dielectric constant in confinement is not enough to explain the enhancement of ion transport in polarizable slit-like channels.

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.

Neural network predicts ion concentration profiles under nanoconfinement

Modeling the ion concentration profile in nanochannel plays an important role in understanding the electrical double layer and electro-osmotic flow. Due to the non-negligible surface interaction and the effect of discrete solvent molecules, molecular dynamics (MD) simulation is often used as an essential tool to study the behavior of ions under nanoconfinement. Despite the accuracy of MD simulation in modeling nanoconfinement systems, it is computationally expensive. In this work, we propose neural network to predict ion concentration profiles in nanochannels with different configurations, including channel widths, ion molarity, and ion types. By modeling the ion concentration profile as a probability distribution, our neural network can serve as a much faster surrogate model for MD simulation with high accuracy. We further demonstrate the superior prediction accuracy of neural network over XGBoost. Finally, we demonstrated that neural network is flexible in predicting ion concentration profiles with different bin sizes. Overall, our deep learning model is a fast, flexible, and accurate surrogate model to predict ion concentration profiles in nanoconfinement.

Water Self-Dissociation is Insensitive to Nanoscale Environments

Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.

Temperature-pressure phase diagram of confined monolayer water/ice at first-principles accuracy with a machine-learning force field

Understanding the phase behaviour of nanoconfined water films is of fundamental importance in broad fields of science and engineering. However, the phase behaviour of the thinnest water film - monolayer water - is still incompletely known. Here, we developed a machine-learning force field (MLFF) at first-principles accuracy to determine the phase diagram of monolayer water/ice in nanoconfinement with hydrophobic walls. We observed the spontaneous formation of two previously unreported high-density ices, namely, zigzag quasi-bilayer ice (ZZ-qBI) and branched-zigzag quasi-bilayer ice (bZZ-qBI). Unlike conventional bilayer ices, few inter-layer hydrogen bonds were observed in both quasi-bilayer ices. Notably, the bZZ-qBI entails a unique hydrogen-bonding network that consists of two distinctive types of hydrogen bonds. Moreover, we identified, for the first time, the stable region for the lowest-density [Formula: see text] monolayer ice (LD-48MI) at negative pressures (<-0.3 GPa). Overall, the MLFF enables large-scale first-principle-level molecular dynamics (MD) simulations of the spontaneous transition from the liquid water to a plethora of monolayer ices, including hexagonal, pentagonal, square, zigzag (ZZMI), and hexatic monolayer ices. These findings will enrich our understanding of the phase behaviour of the nanoconfined water/ices and provide a guide for future experimental realization of the 2D ices.

Unexpected Behavior in Thermal Conductivity of Confined Monolayer Water

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

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

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

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.

Electrically Suppressed Outflow of Confined Liquid in Hydrophobic Nanopores

Confining liquid in a hydrophobic nanoenvironment has enabled a broad spectrum of applications in biomedical sensors, mechanical actuators, and energy storage and converters, where the outflow of confined liquid is spontaneous and fast due to the intrinsic hydrophobic nature of nanopores with extremely low interfacial friction, challenging design capacity and control tolerance of structures and devices. Here, we present a facile approach of suppressing the outflow of water confined in hydrophobic nanopores with an electric field. Extensive molecular dynamics simulations show that the presence of an electric field could significantly strengthen hydrogen bonds and retard degradations of the associated networks during the outflow. The outflow deformation and strength are extracted to quantitatively characterize the electrical suppression to outflow and agree well with simulations. This study proposes a practical means of impeding the fast liquid outflow in hydrophobic nanopores, potentially useful for devising nanofluidics-based functional structures and devices with controllable performance.

2D Material Nanofiltration Membranes: From Fundamental Understandings to Rational Design

Since the discovery of 2D materials, 2D material nanofiltration (NF) membranes have attracted great attention and are being developed with a tremendously fast pace, due to their energy efficiency and cost effectiveness for water purification. The most attractive aspect for 2D material NF membranes is that, anomalous water and ion permeation phenomena have been constantly observed because of the presence of the severely confined nanocapillaries (<2 nm) in the membrane, leading to its great potential in achieving superior overall performance, e.g., high water flux, high rejection rates of ions, and high resistance to swelling. Hence, fundamental understandings of such water and ion transport behaviors are of great significance for the continuous development of 2D material NF membranes. In this work, the microscopic understandings developed up to date on 2D material NF membranes regarding the abnormal transport phenomena are reviewed, including ultrafast water and ion permeation rates with the magnitude several orders higher than that predicted by conventional diffusion behavior, ion dehydration, ionic Coulomb blockade, ion-ion correlations, etc. The state-of-the-art structural designs for 2D material NF membranes are also reviewed. Discussion and future perspectives are provided highlighting the rational design of 2D material membrane structures in the future.

Continuous Water Filling in a Graphene Nanochannel: A Molecular Dynamics Study

Low dimensional materials especially carbon materials hold high promise in the fields of water purification, mineral separation, energy harvesting/conversion, and so on. The fluidic devices fabricated by direct synthesis, lithography, or self-assembly of low dimensional materials provide opportunities for exploring the novel properties and applications of nanoconfined transport. Here, continuous filling of water and acetone molecules into a graphene nanochannel is investigated. A stairlike nonlinear dependence of the number of filling water molecules on interlayer distance d is found when d < 1 nm due to the existence of out-plane layered and in-plane ordered monolayer structure, while near-linear dependence is found for acetone because of the freely rotating configurations along with varying d during the filling process. The entropy, potential energy, and free energy of the confined system during the continuous filling are analyzed to understand the structural evolution of water. The energy-costs are discussed depending on the structure evolution of water during the filling, which is crucial to understanding the swelling and capillary condensation widely existing in the angstrom/nanometer-scale separation membranes.

The physics of cement cohesion

Cement is the most produced material in the world. A major player in greenhouse gas emissions, it is the main binding agent in concrete, providing a cohesive strength that rapidly increases during setting. Understanding how such cohesion emerges is a major obstacle to advances in cement science and technology. Here, we combine computational statistical mechanics and theory to demonstrate how cement cohesion arises from the organization of interlocked ions and water, progressively confined in nanoslits between charged surfaces of calcium-silicate-hydrates. Because of the water/ions interlocking, dielectric screening is drastically reduced and ionic correlations are proven notably stronger than previously thought, dictating the evolution of nanoscale interactions during cement hydration. By developing a quantitative analytical prediction of cement cohesion based on Coulombic forces, we reconcile a fundamental understanding of cement hydration with the fully atomistic description of the solid cement paste and open new paths for scientific design of construction materials.

First-Principles Molecular Dynamics Simulations of the Spontaneous Freezing Transition of 2D Water in a Nanoslit

As with bulk ices, two-dimensional (2D) ices exhibit diverse crystalline structures, and the majority of these 2D structures have been predicted based on classical molecular dynamics (MD) simulations. Here, the spontaneous freezing transition of 2D liquid water within hydrophobic nanoslits is demonstrated for the first time using first-principles MD simulations. Various 2D ices are observed under different lateral pressure and temperature conditions. Notably, the liquid water confined to a 6.0 Å-wide nanoslit can spontaneously freeze into a monolayer ice consisting of an array of zigzag water chains at 2.5 GPa and 250 K. Moreover, within an 8.0 Å-wide nanoslit and at 4.0 GPa and 300 K, a previously unreported bilayer ice forms spontaneously that has a structure resembling that of the double surface layers of bulk ice-VII. Both 2D crystalline ices do not obey the ice rule, suggesting first-principles simulation can access a certain phase space that is not easily approached using classical simulations.

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.

Machine learning for reparameterization of four-site water models: TIP4P-BG and TIP4P-BGT

Parameterizing an effective water model is a challenging issue because of the difficulty in maintaining a comprehensive balance among the diverse physical properties of water with a limited number of parameters. The advancement in machine learning provides a promising path to search for a reliable set of parameters. Based on the TIP4P water model, hence, about 6000 molecular dynamics (MD) simulations for pure water at 1 atm and in the range of 273-373 K are conducted here as the training data. The back-propagation (BP) neural network is then utilized to construct an efficient mapping between the model parameters and four crucial physical properties of water, including the density, vaporization enthalpy, self-diffusion coefficient and viscosity. Without additional time-consuming MD simulations, this mapping operation could result in sufficient and accurate data for high-population genetic algorithm (GA) to optimize the model parameters as much as possible. Based on the proposed parameterizing strategy, TIP4P-BG (a conventional four-site water model) and TIP4P-BGT (an advanced model with temperature-dependent parameters) are established. Both the water models exhibit excellent performance with a reasonable balance among the four crucial physical properties. The relevant mean absolute percentage errors are 3.53% and 3.08%, respectively. Further calculations on the temperature of maximum density, isothermal compressibility, thermal expansion coefficient, radial distribution function and surface tension are also performed and the resulting values are in good agreement with the experimental values. Through this water modeling example, the potential of the proposed data-driven machine learning procedure has been demonstrated for parameterizing a MD-based material model.

Snap-through in Graphene Nanochannels: With Application to Fluidic Control

Recent studies on the structure and transport behaviors of water confined within lamellar graphene have attracted intense interest in filtration technology, but the mechanism of water transport in complex membrane nanostructures remains an open question. For example, similar systems but at much larger scales have indicated that the instabilities of an elastic structure, such as snap-through, play an essential role in controlling the fluid flow. Graphene sheets, which have an atomic thickness, often appear highly wrinkled in nanofluidic devices and so are vulnerable to elastic instabilities. However, it remains unclear how does the flexible wrinkled structure affect the transport of water and filtration efficiency or whether such an effect can be exploited in devices. In this work, we explore the flow-induced snap-through in graphene nanochannels by combining molecular simulations with the theoretical analysis. We further demonstrate its applications to passive control of fluid flow and to ion/molecule selection. By introducing a flexible arch embedded within a graphene nanochannel, we observe the "snap" of the arched graphene wall from one stable state to another by varying the fluid flux (i.e., velocity); the critical velocity of this snap transition is found to depend nonmonotonically on the geometric size of the channel and the arch. We also demonstrate reversible snap-through by fixing the end parts of the flexible arch. These results suggest the potential of flow-induced snap-through in graphene-based nanochannels for ion/molecule selection applications in, for example, the design of a foul-resistant, easy-to-clean, reusable filter membrane.

Molecular Interactions between Asphaltene and Surfactants in a Hydrocarbon Solvent: Application to Asphaltene Dispersion

Heavy oil and bitumen supply the vast majority of energy resources in Canada. Different methods can be implemented to produce oil from such unconventional resources. Surfactants are employed as an additive to water/steam to improve an injected fluid’s effectiveness and enhance oil recovery. One of the main fractions in bitumen is asphaltene, which is a non-symmetrical molecule. Studies of interactions between surfactants, anionic, and non-anionic, and asphaltene have been very limited in the literature. In this paper, we employed molecular dynamics (MD) simulation to theoretically focus on the interactions between surfactant molecules and different types of asphaltene molecules observed in real oil sands. Both non-anionic and anionic surfactants showed promising results in terms of dispersant efficiency; however, their performance depends on the asphaltene architecture. Moreover, a hydrogen/carbon (H/C) ratio of asphaltenes plays an inevitable role in asphaltene aggregation behavior. A higher H/C ratio resulted in decreasing asphaltene aggregation tendency. The results of these studies will give a deep understanding of the interactions between asphaltene and surfactant molecules.

Wettability and confinement size effects on stability of water conveying nanotubes

This study investigates the wettability and confinement size effects on vibration and stability of water conveying nanotubes. We present an accurate assessment of nanotube stability by considering the exact mechanics of the fluid that is confined in the nanotube. Information on the stability of nanotubes in relation to the fluid viscosity, the driving force of the fluid flow, the surface wettability of the nanotube, and the nanotube size is missing in the literature. For the first time, we explore the surface wettability dependence of the nanotube natural frequencies and stability. By means of hybrid continuum-molecular mechanics (HCMM), we determined water viscosity variations inside the nanotube. Nanotubes with different surface wettability varying from super-hydrophobic to super-hydrophilic nanotubes were studied. We demonstrated a multiphase structure of nanoconfined water in nanotubes. Water was seen as vapor at the interface with the nanotube, ice shell in the middle, and liquid water in the nanotube core. The average velocity of water flow in the nanotube was obtained strongly depend on the surface wettability and the confinement size. In addition, we report the natural frequencies of the nanotube as functions of the applied pressure and the nanotube size. Mode divergence and flutter instabilities were observed, and the activation of these instabilities strongly depended on the nanotube surface wettability and size. This work gives important insights into understanding the stability of nanotubes conveying fluids depending on the operating pressures and the wettability and size of confinement. We revealed that hydrophilic nanotubes are generally more stable than hydrophobic nanotubes when conveying fluids.

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.

Microscopic Pressure Tensor in Cylindrical Geometry: Pressure of Water in a Carbon Nanotube

The microscopic pressure tensor plays an important role in understanding the mechanical stability, transport, and high-pressure phenomena of confined phases. The lack of an exact formulation to account for the long-range Coulombic contribution to the local pressure tensor in cylindrical geometries prevents the characterization of molecular fluids confined in cylindrical pores. To address this problem, we first derive the local cylindrical pressure tensor for Lennard-Jones fluids based on the Harasima (H) definition, which is expected to be compatible with the Ewald summation method. The test of the H-definition pressure equations in a homogeneous system shows that the radial and azimuthal pressure have unphysical radial dependence near the origin, while the axial pressure gives physically meaningful values. We propose an alternative contour definition that is more appropriate for cylindrical geometry and show that it leads to physically realistic results for all three pressure tensor components. With this definition, the radial and azimuthal pressures are of Irving-Kirkwood (IK) type, and the axial pressure is of Harasima type. Because of the practical interest in the axial pressure, we develop a Harasima/Ewald (H/E) method for calculating the long-range Coulombic contribution to the local axial pressure for rigid molecules. As an application, the axial pressure profile of water inside and outside a (20, 20) single-wall carbon nanotube is determined. The H/E method is compared to the IK method, which assumes a spherically truncated Coulombic potential. Detailed analysis of the pressure profile by both methods shows that the water confined in the nanotube is in a stretched state overall in the axial direction.

Refractive index of nanoconfined water reveals its anomalous physical properties

Despite extensive studies on the distinctive properties of water confined in a nanospace, the underlying mechanism and significance of the lengthscale involved in the confinement effects are still subjects of controversy. The dielectric constant and the refractive index in particular are key parameters in modeling and understanding nanoconfined water, yet experimental evidence is lacking. We report the measurement of the refractive indices of water in 10-100 nm spaces by exploiting the confinement of water and localized surface plasmons in a physicochemically well-defined nanocavity. The results revealed significantly low values and the scaling behavior of the out-of-plane refractive index n⊥ of confined water. They are attributed to the polarization suppression at the interfaces and the long-range correlation in electronic polarization facilitated by the strengthened H-bonding network. Using the refractive index as a sensing probe, we also observed anomalous stability of water structures over a wide range of temperature. Our measurement results provide essential feedback information for benchmarking water models and molecular interactions under nanoconfinement. This study also opens up a new methodology of using plasmon resonance in characterizing nanoconfined molecules and chemical reactions, and thus gives us fundamental insight into confinement effects.

Two-Phase Equilibrium Conditions in Nanopores

It is known that thermodynamic properties of a system change upon confinement. To know how, is important for modelling of porous media. We propose to use Hill's systematic thermodynamic analysis of confined systems to describe two-phase equilibrium in a nanopore. The integral pressure, as defined by the compression energy of a small volume, is then central. We show that the integral pressure is constant along a slit pore with a liquid and vapor in equilibrium, when Young and Young-Laplace's laws apply. The integral pressure of a bulk fluid in a slit pore at mechanical equilibrium can be understood as the average tangential pressure inside the pore. The pressure at mechanical equilibrium, now named differential pressure, is the average of the trace of the mechanical pressure tensor divided by three as before. Using molecular dynamics simulations, we computed the integral and differential pressures, p ^ and p, respectively, analysing the data with a growing-core methodology. The value of the bulk pressure was confirmed by Gibbs ensemble Monte Carlo simulations. The pressure difference times the volume, V, is the subdivision potential of Hill, ( p - p ^ ) V = ϵ . The combined simulation results confirm that the integral pressure is constant along the pore, and that ϵ / V scales with the inverse pore width. This scaling law will be useful for prediction of thermodynamic properties of confined systems in more complicated geometries.

Closest-Packing Water Monolayer Stably Intercalated in Phyllosilicate Minerals under High Pressure

The directional hydrogen-bond (HB) network and nondirectional van der Waals (vdW) interactions make up the specificity of water. Directional HBs could construct an ice-like monolayer in hydrophobic confinement even in the ambient regime. Here, we report a water monolayer dominated by vdW interactions confined in a phyllosilicate interlayer under high pressure. Surprisingly, it was in a thermodynamically stable state coupled with bulk water at the same pressure (P) and temperature (T), as revealed by the thermodynamic integration approach on the basis of molecular dynamics (MD) simulations. Both classical and ab initio MD simulations showed water O atoms were stably trapped and exhibited an ordered hexagonal closest-packing arrangement, but OH bonds of water reoriented frequently and exhibited a specific two-stage reorientation relaxation. Strikingly, hydration in the interlayer under high pressure had no relevance with surface hydrophilicity rationalized by the HB forming ability, which, however, determines wetting in the ambient regime. Intercalated water molecules were trapped by vdW interactions, which shaped the closest-packing arrangement and made hydration energetically available. The high pressure-volume term largely drives hydration, as it compensates the entropy penalty which is restricted by a relatively lower temperature. This vdW water monolayer should be ubiquitous in the high pressure but low-temperature regime.

Carbon Nanocomposite Membrane Electrolytes for Direct Methanol Fuel Cells-A Concise Review

A membrane electrolyte that restricts the methanol cross-over while retaining proton conductivity is essential for better electrochemical selectivity in direct methanol fuel cells (DMFCs). Extensive research carried out to explore numerous blends and composites for application as polymer electrolyte membranes (PEMs) revealed promising electrochemical selectivity in DMFCs of carbon nanomaterial-based polymer composites. The present review covers important literature on different carbon nanomaterial-based PEMs reported during the last decade. The review emphasises the proton conductivity and methanol permeability of nanocomposite membranes with carbon nanotubes, graphene oxide and fullerene as additives, assessing critically the impact of each type of filler on those properties.

Replica exchange MD simulations of two-dimensional water in graphene nanocapillaries: rhombic versus square structures, proton ordering, and phase transitions

The hydrogen bond patterns, proton ordering, and phase transitions of monolayer ice in two-dimensional hydrophobic confinement are fundamentally different from those found for bulk ice. To investigate the behavior of quasi-2D ice, we perform molecular dynamics simulations of water confined between fixed graphene plates at a distance of 0.65 nm. While experimental results are still limited and theoretical investigations are often based on a single, often empirically based force field model, this work presents a systematic study modeling the water-graphene interaction by effective Lennard-Jones potentials previously derived from high-level ab initio CCSD(T) calculations of water adsorbed on graphene [Phys. Chem. Chem. Phys., 2013, 15, 4995]. For the water-water interaction different water force fields, i.e. SPCE, TIP3P, TIP4P, TIP4P/ICE, and TIP5P, are used. The water occupancy of the graphene capillary at a pressure of 1000 MPa is determined to be between 13.5 and 13.9 water molecules per square nanometer, depending on the choice of the water force field. Based on these densities, we explore the structure and dynamics of quasi-2D water for temperatures ranging from 200 K to about 600 K for each of the five force fields. To ensure complete sampling of the configurational space and to overcome the barriers separating metastable structures, these simulations are based on the replica exchange molecular dynamics technique. We report different tetragonal hydrogen bond patterns, which are classified as nearly square or as rhombic. While many of these arrangements are essentially flat, in some cases puckered arrangements are found, too. Also the proton ordering of the quasi-2D water structures is considered, allowing us to identify them as ferroelectric, ferrielectric or antiferroelectric. For temperatures between 200 K and 400 K we find several second-order phase transitions from one ice structure to another, changing in many cases both the arrangements of the oxygen atoms and the proton ordering. For temperatures between 400 K and 600 K there are melting-like transitions from a monolayer of ice to a monolayer of liquid water. These first-order phase transitions have a latent heat between 3.4 and 4.0 kJ mol^-1. Both the values of the transition temperatures and of the latent heats display considerable model dependence for the five different water models investigated here.

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

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

Crystallization behavior of a confined CuZr metallic liquid film with a sandwich-like structure

Despite the fact that its crystal state is thermodynamically stable, Cu₆₄Zr₃₆ alloy is prone to form metastable glass at a high cooling rate. However, the confinement can induce nano-crystallization with a novel sandwich-like hierarchical structure consisting of pure Cu layers, pure Zr layers and mixed layers by conducting molecular dynamics simulations. The liquid-to-crystal transition temperature and interatomic repulsion softness display abnormal oscillations, instead of monotonous variation, as the wall-wall separation increases. When the confinement size is 10 Å and 12 Å, the transition temperature reaches a maximum, resulting from the pending new sandwich layer. The atomic movement and dynamical heterogeneity are demonstrated to play a vital role in the abnormal oscillation behavior of physical properties of the nano confined metallic glass. The sandwich-like structure can alter the Cu-Zr bond fraction, which eventually influences the liquid-to-crystal transition temperature and interatomic repulsion softness. Our findings provide a deep insight into the hierarchical nanostructures and its liquid-to-crystal transition characteristics under confinement at the atomic level.

Molecular dynamics study of nanoconfined TIP4P/2005 water: how confinement and temperature affect diffusion and viscosity

In the past few decades great effort has been devoted to the study of water confined in hydrophobic geometries at the nanoscale (tubes and slit pores) due to the multiple technological applications of such systems, ranging from drug delivery to water desalination devices. To our knowledge, neither numerical/theoretical nor experimental approaches have so far reached a consensual understanding of structural and transport properties of water under these conditions. In this work, we present molecular dynamics simulations of TIP4P/2005 water under different nanoconfinements (slit pores or nanotubes, with two degrees of hydrophobicity) within a wide temperature range. It has been found that water is more structured near the less hydrophobic walls, independently of the confining geometries. Meanwhile, we observe an enhanced diffusion coefficient of water in both hydrophobic nanotubes. Finally, we propose a confined Stokes-Einstein relation to obtain the viscosity from diffusivity, whose result strongly differs from the Green-Kubo expression that has been used in previous works. While viscosity computed with the Green-Kubo formula (applied for anisotropic and confined systems) strongly differs from that of the bulk, viscosity computed with the confined Stokes-Einstein relation is not so much affected by the confinement, independently of its geometry. We discuss the shortcomings of both approaches, which could explain this discrepancy.

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.

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.

Molecular dynamics simulations of nano-confined methanol and methanol-water mixtures between infinite graphite plates: Structure and dynamics

Molecular dynamics simulations are used to investigate microscopic structures and dynamics of methanol and methanol-water binary mixture films confined between hydrophobic infinite parallel graphite plate slits with widths, H, in the range of 7-20 Å at 300 K. The initial geometric densities of the liquids were chosen to be the same as bulk methanol at the same temperature. For the two narrowest slit widths, two smaller initial densities were also considered. For the nano-confined system with H = 7 Å and high pressure, a solid-like hexagonal arrangement of methanol molecules arranged perpendicular to the plates is observed which reflects the closest packing of the molecules and partially mirrors the structure of the underlying graphite structure. At lower pressures and for larger slit widths, in the contact layer, the methanol molecules prefer having the C-O bond oriented parallel to the walls. Layered structures of methanol parallel to the wall were observed, with contact layers and additional numbers of central layers depending on the particular slit width. For methanol-water mixtures, simulations of solutions with different composition were performed between infinite graphite slits with H = 10 and 20 Å at 300 K. For the nanoslit with H = 10 Å, in the solution mixtures, three layers of molecules form, but for all mole fractions of methanol, methanol molecules are excluded from the central fluid layer. In the nanopore with H = 20 Å, more than three fluid layers are formed and methanol concentrations are enhanced near the confining plates walls compared to the average solution stoichiometry. The self-diffusion coefficients of methanol and water molecules in the solution show strong dependence on the solution concentration. The solution mole fractions with minimal diffusivity are the same in confined and non-confined bulk methanol-water mixtures.

Fluidity and phase transitions of water in hydrophobic and hydrophilic nanotubes

We put water flow under scrutiny to report radial distributions of water viscosity within hydrophobic and hydrophilic nanotubes as functions of the water-nanotube interactions ([Formula: see text]), surface wettability (θ), and nanotube size (R) using a proposed hybrid continuum-molecular mechanics. Based on the computed viscosity data, [Formula: see text] phase diagram of the phase transitions of confined water in nanotubes is developed. It is revealed that water exhibits different multiphase structures, and the formation of one of these structures depends on [Formula: see text] R parameters. A drag of water flow at the first water layer is revealed, which is conjugate to sharp increase in the viscosity and formation of an ice phase under severe confinement (R ≤ 3.5 nm) and strong water-nanotube interaction conditions. A vapor/vapor-liquid phase is observed at hydrophobic and hydrophilic interfaces. A state of confinement is revealed at which water exhibits different multiphase structures under the same flow rate. The derived viscosity functions are used to accurately determine factors of flow enhancement/inhibition of confined water.

Cosolvent Impurities in SWCNT Nanochannel Confinement: Length Dependence of Water Dynamics Investigated with Atomistic Simulations

The advent of nanotechnology has seen a growing interest in the nature of fluid flow and transport under nanoconfinement. The present study leverages fully atomistic molecular dynamics (MD) simulations to study the effect of nanochannel length and intrusion of molecules of the organic solvent, hexafluoro-2-propanol (HFIP), on the dynamical characteristics of water within it. Favorable interactions of HFIP with the nanochannels comprised of single-walled carbon nanotubes traps them over time scales greater than 100 ns, and confinement confers small but distinguishable spatial redistribution between neighboring HFIP pairs. Water molecules within the nanochannels show clear signatures of dynamical slowdown relative to bulk water even for pure systems. The presence of HFIP causes further rotational and translational slowdown in waters when the nanochannel dimension falls below a critical length of 30 Å. The enhanced slowdown in the presence of HFIP is quantified from characteristic relaxation parameters and diffusion coefficients in the absence and presence of HFIP. It is finally seen that the net flow of water between the ends of the nanochannel shows a decreasing dependence with nanochannel length only when the number of HFIP molecules is small. These results lend insights into devising ways of modulating solvent properties within nanochannels with cosolvent impurities.

Comparative Study of the Effects of Temperature and Pressure on the Water-Mediated Interactions between Apolar Nanoscale Solutes

We perform molecular dynamics simulations to study the effects of temperature and pressure on the water-mediated interaction (WMI) between two nanoscale (apolar) graphene plates at 240 ≤ T ≤ 400 K and -100 ≤ P ≤ 1200 MPa. These are thermodynamic conditions relevant to, for example, cooling-, heating-, compression-, and decompression-induced protein denaturation. We find that at all ( T, P) studied, the potential of mean force between the graphene plates, as a function of plate separation r, exhibits local minima at specific plate separations r = r _n that can accommodate n water layers ( n = 0,1,2,3). In particular, our results show that isobaric cooling and isothermal compression have a similar effect on WMI between the plates; both processes tend to suppress the attraction and ultimate collapse of the graphene plates by kinetically trapping the plates at the metastable states with r = r _n ( n > 0). In addition, isobaric heating and isothermal decompression also have a similar effect; both processes tend to reduce the range and strength of the interactions between the graphene plates. Interestingly, at low temperatures, the WMI between the plates is affected by crystallization. However, crystallization depends deeply on the water model considered, SPC/E and TIP4P/2005 water models, with the crystallization occurring at different ( T, P) conditions, into different forms of ice.