Phase Diagram of Water Confined by Graphene

Water molecules in boron nitride interlayer space: ice and hydrolysis in super confinement

Development of nano-sized channels and filters in the recent years has made the role of water immensely important as water molecules affect their performance and durability. Here, we take advantage of molecular dynamics and density functional theory methods to demonstrate the shift in behavior of water molecules confined between hexagonal boron nitride (HBN) sheets spaced at 3.0 to 6.5 Å. Our results demonstrate that lower interlayer spaces cause higher amounts of charge transferred between the species, while at extreme degrees of confinement, these interactions cause the disintegration of trapped water molecules. Consequently, the inner face of the HBN sheets is functionalized with hydroxyl groups, releasing hydrogens in the form of protons that travel the interlayer space by Grotthuss mechanism. This is the first-hand evidence of a mechanical form of hydrolysis that corresponds with a nucleophilic attack (on boron atoms) to relieve water from extreme confined conditions. This process unveils a previously unknown behavior of water within extremely confined spaces and reveals new considerations concerning nanofilters and nanochannels.

Structure and dynamics of amphiphilic patchy cubes in a nanoslit under shear

Patchy nanocubes are intriguing materials with simple shapes and space-filling and multidirectional bonding properties. Previous studies have revealed various mesoscopic structures such as colloidal crystals in the solid regime and rod-like or fractal-like aggregates in the liquid regime of the phase diagram. Recent studies have also shown that mesoscopic structural properties, such as an average cluster size M and orientational order, in amphiphilic nanocube suspensions are associated with macroscopic viscosity changes, mainly owing to differences in cluster shape among patch arrangements. Although many studies have been conducted on the self-assembled structures of nanocubes in bulk, little is known about their self-assembly in nanoscale spaces or structural changes under shear. In this study, we investigated mixtures of one- and two-patch amphiphilic nanocubes confined in two flat parallel plates at rest and under shear using molecular dynamics simulations coupled with multiparticle collision dynamics. We considered two different patch arrangements for the two-patch particles and two different slit widths H to determine the degree of confinement in constant volume fractions in the liquid regime of the phase diagram. We revealed two unique cluster morphologies that have not been previously observed under bulk conditions. At rest, the size of the rod-like aggregates increased with decreasing H, whereas that of the fractal-like aggregates remained constant. Under weak shear with strong confinement, the rod-like aggregates maintained a larger M than the fractal-like aggregates, which were more rigid and maintained a larger M than the rod-like aggregates under bulk conditions.

Observation of the early stages of environmental contamination in graphene by friction force

Exposure to ambient air contaminates the surface of graphene sheets. Contamination may arise from different sources, and its nature alters the frictional behavior of the material. These changes in friction enable the observation of the early stages of contaminants' adsorption in graphene. Using a friction force microscope, we show that molecular adsorption initiates at the edges and mechanical defects in the monolayer. Once the monolayer is covered, the contaminants spread over the additional graphene layers. With this method, we estimate the contamination kinetics. In monolayer graphene, the surface area covered with adsorbed molecules increases with time of air exposure at a rate of 10-14 m2/s, while in bilayer graphene, it is one order of magnitude smaller. Finally, as the contaminants cover the additional graphene layers, friction no longer has a difference concerning the number of graphene layers.

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

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

Characterization of Water Structure and Phase Behavior within Metal-Organic Nanotubes

Water behavior under nanoconfinement varies significantly from that in the bulk but also depends on the nature of the pore walls. Hybrid compound offers the ideal system to explore water behavior in complex materials, so a model metal-organic nanotube (UMONT) material was utilized to explore the behavior of water between 100 and 293 K. Single-crystal X-ray and neutron diffraction revealed the formation of a filled Ice-I arrangement that was previously predicted to only occur under high pressures. 17O NMR spectra suggest that the onset of melting for the water in the UMONT channels occurs at 98 K and the presence of ice-like water up to 293 K, indicating that the complete ice-water transition does not occur before dehydration of the material. Overall, the water behavior differs significantly from hydrophobic single-walled carbon nanotubes indicating precise control over water can be achieved through rational design of hybrid materials.

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.

Hydrophobic Nanoconfinement Enhances CO2 Conversion to H2CO3

Understanding the formation of H₂CO₃ in water from CO₂ is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO₂ to H₂CO₃ in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO₂ conversion to H₂CO₃ in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO₂ concentration.

Impacts of Surface Adsorption on Water Uptake within a Metal Organic Nanotube Material

The confinement-dependent properties of solvents, particularly water, within nanoporous spaces impart unique physical and chemical behavior compared to those of the bulk. This has previously been demonstrated for a U(VI)-based metal organic nanotube that displays ice-like arrays of water molecules within the 1-D pore space and complete selectivity to H₂O over all other solvents and isotopologues. Based upon our previous work on D₂O and HTO adsorption processes, we suggested that the water uptake was controlled by a two-step process: (1) surface adsorption via hydrogen bonding to hydrophilic amine and carboxylic groups and (2) diffusion of the water into the hydrophobic 1-D nanochannels. The current study seeks to evaluate this hypothesis and expand our existing kinetic model for the water diffusion step to account for the initial surface adsorption process. Vapor sorption experiments, paired with thermogravimetric and Fourier-transform infrared analyses, yielded uptake data that were fit using a Langmuir model for the surface-adsorption step of the mechanism. The water adsorption curve was designated a type IV Brunauer-Emmett-Teller isotherm, which indicated that our original hypothesis was correct. Additional work with binary solvent systems enabled us to evaluate the uptake in a range of conditions and determine that the uptake is not controlled by the vapor pressure but is instead completely dependent on the relative humidity of the system.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

Nanoconfinement facilitates reactions of carbon dioxide in supercritical water

The reactions of CO2 in water under extreme pressure-temperature conditions are of great importance to the carbon storage and transport below Earth's surface, which substantially affect the carbon budget in the atmosphere. Previous studies focus on the CO2(aq) solutions in the bulk phase, but underground aqueous solutions are often confined to the nanoscale, and nanoconfinement and solid-liquid interfaces may substantially affect chemical speciation and reaction mechanisms, which are poorly known on the molecular scale. Here, we apply extensive ab initio molecular dynamics simulations to study aqueous carbon solutions nanoconfined by graphene and stishovite (SiO2) at 10 GPa and 1000 ~ 1400 K. We find that CO2(aq) reacts more in nanoconfinement than in bulk. The stishovite-water interface makes the solutions more acidic, which shifts the chemical equilibria, and the interface chemistry also significantly affects the reaction mechanisms. Our findings suggest that CO2(aq) in deep Earth is more active than previously thought, and confining CO2 and water in nanopores may enhance the efficiency of mineral carbonation.

Three-Dimensional Confinement of Water: H2O Exhibits Long-Range (>50 nm) Structure while D2O Does Not

Water is the liquid of life thanks to its three-dimensional adaptive hydrogen (H)-bond network. Confinement of this network may lead to dramatic structural changes influencing chemical and physical transformations. Although confinement effects occur on a <1 nm length scale, the upper length scale limit is unknown. Here, we investigate 3D-confinement over lengths scales ranging from 58-140 nm. By confining water in zwitterionic liposomes of different sizes and measuring the change in H-bond network conformation using second harmonic scattering (SHS), we determined long-range confinement effects in light and heavy water. D₂O displays no detectable 3D-confinement effects <58 nm (<3 × 10⁶ D₂O molecules). H₂O is distinctly different. The vesicle enclosed inner H-bond network has a different conformation compared to the outside network and the SHS response scales with the volume of the confining space. H₂O displays confinement effects over distances >100 nm (>2 × 10⁷ H₂O molecules).

Correlation between ordering and shear thinning in confined OMCTS liquids

Despite the extensive research that has been conducted for decades on the behavior of confined liquids, detailed knowledge of this phenomenon, particularly in the mixed/boundary lubrication regime, remains limited. This can be attributed to several factors including the difficulty of direct experimental observations of the behavior of lubricant molecules under non-equilibrium conditions, the high computational cost of molecular simulations to reach steady state, and the low signal-to-noise ratio at extremely low shear rates corresponding to actual operating conditions. To this end, we studied the correlation between the structure formation and shear viscosity of octamethylcyclotetrasiloxane confined between two mica surfaces in a mixed/boundary lubrication regime. Three different surface separations corresponding to two-, three-, and five-layered structures were considered to analyze the effect of confinement. The orientational distributions with one specific peak for n=2 and two distributions, including a parallel orientation with the surface normal for n&gt;2, were observed at rest. The confined liquids exhibited a distinct shear-thinning behavior independent of surface separations for a relatively low shear rate, \dot{\gamma}&lt;10^{8} s^{-1}. However, the shear viscosities at \dot{\gamma}&lt;10^{8} s^{-1} depended on the number of layered structures. Newtonian behavior was observed with a further increase in the shear rate. Furthermore, we found a strong correlation between the degree of molecular orientation and the shear viscosity of the confined liquids. The magnitude of the shear viscosity of the confined liquids can primarily be determined by the degree of molecular orientation, and shear-thinning originates from the vanishing of specific orientational distributions with increasing shear rate.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

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

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

Study on Flow Distribution Pattern and Conductivity of Porous Media in Bioretention Cells

To evaluate the long-term performance of bioretention cell (BRC), a study was undertaken to assess the flow distribution and conductivity. Despite initial conductivity of the original medium being the common predictor of hydraulic performance, most of the BRCs are affected by conductivity variations during actual operation. This happen due to the fact that microbial behavior plays an important role in the conductivity variations. This linkage may occur when bacteria as inert colloids transports between particles and biodegrades dissolved pollutants, either promoting or retarding flow distribution and conductivity in BRC. Flow distribution was determined by numerical simulation and tracer test, and the correlation between conductivity and flow distribution was revealed by conductivity experiment coupled with flow distribution analysis. Results revealed a non-uniform flow distribution in BRC, and seepage flow in submerged zone was virtually impossible push flow. Conductivity had an inversely proportional relationship with hydraulic efficiency where hydraulic efficiency reached the highest value (0.297) under a low hydraulic conductivity (0.000107 m/s, approximately K/K_ini = 0.79). Primary cause of hydraulic capacity reduction was the initial permeability decrease due to medium structure changes. Results revealed a sharp upward trend followed by a slight decrease, and then, stabilized to a stable infiltration stage. Permeation process of sewage influent was similar to the one of potable water where the permeability reduced to 0.000102 m/s after 450 h and declined continuously. Thus, it is clear that flow distribution and conductivity in bioretention must be estimated more accurately on a microscopic scale.

Anomalously low friction of confined monolayer water with a quadrilateral structure

In this work, we explored how the structure of monolayer water confined between two graphene sheets is coupled to its dynamic behavior. Our molecular dynamics simulations show that there is a remarkable interrelation between the friction of confined water with two walls and its structure under extreme confinement. When the water molecules formed a regular quadrilateral structure, the friction coefficient is dramatically reduced. Such a low-friction coefficient can be attributed to the formation of long-range ordered hydrogen bond network, which not only decreases the structure corrugation in the direction perpendicular to the walls but also promotes the collective motion of the confined water. The regular quadrilateral structure can be formed only if the number density of confined water falls within a certain range. Higher number density results in larger structure corrugations, which increases the friction, while smaller number density leads to an irregular hydrogen bond network in which the collective motion cannot play the role. We demonstrated that there are four distinct stages in the diagram of the friction coefficient vs the number density of confined water. This research clearly established the connection between the dynamic characteristics of confined monolayer water and its structure, which is beneficial to further understand the mechanism of the high-speed water flow through graphene nanocapillaries observed in recent experiments.

Confinement-Controlled Aqueous Chemistry within Nanometric Slit Pores

In this Focus Review, we put the spotlight on very recent insights into the fascinating world of wet chemistry in the realm offered by nanoconfinement of water in mechanically rather rigid and chemically inert planar slit pores wherein only monolayer and bilayer water lamellae can be hosted. We review the effect of confinement on different aspects such as hydrogen bonding, ion diffusion, and charge defect migration of H⁺(aq) and OH^-(aq) in nanoconfined water depending on slit pore width. A particular focus is put on the strongly modulated local dielectric properties as quantified in terms of anisotropic polarization fluctuations across such extremely confined water films and their putative effects on chemical reactions therein. The stunning findings disclosed only recently extend wet chemistry in particular and solvation science in general toward extreme molecular confinement conditions.

Mechanical hydrolysis imparts self-destruction of water molecules under steric confinement

Decoding behavioral aspects associated with the water molecules in confined spaces such as an interlayer space of two-dimensional nanosheets is key for the fundamental understanding of water-matter interactions and identifying unexpected phenomena of water molecules in chemistry and physics. Although numerous studies have been conducted on the behavior of water molecules in confined spaces, their reach stops at the properties of the planar ice-like formation, where van der Waals interactions are the predominant interactions and many questions on the confined space such as the possibility of electron exchange and excitation state remain unsettled. We used density functional theory and reactive molecular dynamics to reveal orbital overlap and induction bonding between water molecules and graphene sheets under much less pressure than graphene fractures. Our study demonstrates high amounts of charge being transferred between water and the graphene sheets, as the interlayer space becomes smaller. As a result, the inner face of the graphene nanosheets is functionalized with hydroxyl and epoxy functional groups while released hydrogen in the form of protons either stays still or traverses a short distance inside the confined space via the Grotthuss mechanism. We found signatures of a new hydrolysis mechanism in the water molecules, i.e. mechanical hydrolysis, presumably responsible for relieving water from extremely confined conditions. This phenomenon where water reacts under extreme confinement by disintegration rather than forming ice-like structures is observed for the first time, illustrating the prospect of treating ultrafine porous nanostructures as a driver for water splitting and material functionalization, potentially impacting the modern design of nanofilters, nanochannels, nano-capacitators, sensors, and so on.

Binding Free Energy Calculations of Bilayer Graphenes Using Molecular Dynamics

Bilayer graphenes are dimeric assemblies of single graphene layers bound together by π-complexation interactions. Controlling these assemblies can be complicated, as the layered compounds disperse in solvents or aggregate into higher columnar configurations and clusters. One way to assess the interactions that contribute to the stability of the layered compounds is to use molecular simulation. We perform pulling molecular dynamics on bilayer graphenes with different sizes and obtain the normal and shear force profiles of dissociation. We generate pathways of dissociation along the two directions and calculate the binding free energies of the structures with umbrella sampling simulations. We show that the dissociation process is direction-dependent. Along the shear direction, we compute the same free energy for the different samples, which validates the consistency of our simulations. We notice that the dissociation is less adiabatic on the normal than the shear direction, having an entropic contribution to the Gibbs energy. This contribution is more enhanced for the larger bilayer graphenes.

Role of Nanoscale Interfacial Proximity in Contact Freezing in Water

Contact freezing is a mode of atmospheric ice nucleation in which a collision between a dry ice nucleating particle (INP) and a water droplet results in considerably faster heterogeneous nucleation. The molecular mechanism of such an enhancement is, however, still a mystery. While earlier studies had attributed it to collision-induced transient perturbations, recent experiments point to the pivotal role of nanoscale proximity of the INP and the free interface. By simulating the heterogeneous nucleation of ice within INP-supported nanofilms of two model water-like tetrahedral liquids, we demonstrate that such nanoscale proximity is sufficient for inducing rate increases commensurate with those observed in contact freezing experiments, but only if the free interface has a tendency to enhance homogeneous nucleation. Water is suspected of possessing this latter property, known as surface freezing propensity. Our findings therefore establish a connection between the surface freezing propensity and kinetic enhancement during contact nucleation. We also observe that faster nucleation proceeds through a mechanism markedly distinct from classical heterogeneous nucleation, involving the formation of hourglass-shaped crystalline nuclei that conceive at either interface and that have a lower free energy of formation due to the nanoscale proximity of the interfaces and the modulation of the free interfacial structure by the INP. In addition to providing valuable insights into the physics of contact nucleation, our findings can assist in improving the accuracy of heterogeneous nucleation rate measurements in experiments and in advancing our understanding of ice nucleation on nonuniform surfaces such as organic, polymeric, and biological materials.

Controlling water structure and behavior: design principles from metal organic nanotubular materials

Water exhibits unique and unexpected behavioral and structural changes when confined to the nanoscale, notably within the pores of metal–organic nanotubes.

Water as a tuneable solvent: a perspective

Water is the most sustainable solvent, but its polarity limits the solubility of non-polar solutes. Confining water in hydrophobic nanopores could be a way to modulate water solvent properties and enable using water as tuneable solvent (WaTuSo).

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

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

Ice Nucleation of Confined Monolayer Water Conforms to Classical Nucleation Theory

We confirmed that monolayer water confined by parallel graphene sheets spontaneously crystallizes from a structurally and dynamically heterogeneous liquid phase under moderate supercooling via direct molecular dynamics simulation. Square-lattice-like geometric order is observed at the early stage of nucleation and is preserved during the entire nucleus growth process. The diffusion coefficient and free energy profile in the cluster space extracted from a Bayesian trajectory analysis agree well with the classical nucleation theory (CNT) prediction and yield thermodynamic quantities exhibiting linear temperature dependence. The effectiveness of maximum cluster size as the descriptor of ice nucleation dynamics in the CNT framework can be attributed to the dynamical time scale decoupling and strong structural pattern dependence of density fluctuation in the liquid phase.

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.

Novel pressure-induced topological phase transitions of supercooled liquid and amorphous silicene

This molecular dynamics (MD) simulation carries a detailed analysis of a pressure-induced structural transition supercooled liquid and amorphous silicene (a-silicene). Low-density models of supercooled liquid and a-silicene containing 10 000 atoms are obtained by rapid cooling processes from the melts. Then, an a-silicene model at T  =  1000 K, a supercooled liquid model at T  =  1500 K and a liquid silicon model at T  =  2000 K have been isothermally compressed step by step up to a high density in order to observe the pressure-induced structural changes. Specifically 'Cairo tiling' pentagonal and square lattices of silicene are discovered in our calculations. Structural properties of those penta-silicene and tetra-silicene models have been carefully analyzed through the radial distribution functions, interatomic distances, bond-angle distributions under high-pressure condition. The dependence of pressure on formation behaviors is calculated via pressure-volume and energy-density relationships. The first order transition from low-density supercooled liquid/amorphous silicene to high-density penta-silicene and continuous transition from low-density liquid to high-density tetra-silicene are discussed. Atomic mechanism and sp³/sp² hybridization evolution are inspected whereas the role of low-membered ring defects/boundary promises remarkable application and advanced research in future.

Chemistry in confined space through the eyes of surface science-2D porous materials

There are a rapidly growing number of studies showing exciting new opportunities in the way confinement effects on surfaces affect the properties of materials and their chemistry. These effects have been observed recently under two-dimensional (2D) van der Waals materials such as a graphene and boron nitride and for the case of supported 2D-porous oxides, including silicates, aluminosilicates and zeolite nanosheets. This review summarizes the current state of the art in this area of research and how confinement effects in 2D systems relate to those found in 3D porous and layered materials. The focus of this review is put in 2D-materials with inherent porosity, such as 2D-porous oxides. An outlook is also given for the future of this exciting emerging area.

Water flow modeling through a graphene-based nanochannel: theory and simulation

Understanding the behavior of water molecule transport through artificial nano-channels is essential in designing novel nanofluidic devices that could be used especially in nanofiltration processes. In this study, using nonequilibrium molecular dynamics (MD) simulations, we simulated the water flow through different graphene-based channels to investigate the influences of some key factors such as the channel thickness and applied pressure on the water flow. It was demonstrated that the water flow was enhanced by increasing the applied pressure and channel thickness. Our results indicated that a third order polynomial curve could describe the variation of the water flow as a function of the channel thickness and the applied pressure. In addition, we improved the hydrodynamics equation used to consider the water flow through nano-channels, by adding two terms to describe the slip effect and the entrance/exit effect, in which the first term increased the water flow rate, while the second term reduced it. This study may be helpful in designing high-performance graphene-based membranes with some practical applications such as desalination.

Nanoconfined Water within Graphene Slit Pores Adopts Distinct Confinement-Dependent Regimes

In view of the increasing importance of nanoconfined aqueous solutions for various technological applications, it has become necessary to understand how strong confinement affects the properties of water at the level of molecular and even electronic structure. By performing extensive ab initio simulations of two-dimensionally nanoconfined water lamellae between graphene sheets subject to different interlayer spacings, we find new regimes at interlayer distances of 10 Å and less where water can be described neither to behave like interfacial water nor to be bulklike at the level of its H-bonding characteristics and electronic structure properties. It is expected that this finding will offer new opportunities to tune both diffusive and reactive processes taking place in aqueous environments that are strongly confined by chemically inert hard walls.

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.

Adsorption of water, methanol, and their mixtures in slit graphite pores

The behavior of water, methanol, and water-methanol mixtures confined in narrow slit graphite pores as a function of pore size was investigated by Monte Carlo, hybrid Monte Carlo, and Molecular Dynamics simulations. Interactions were described using TIP4P/2005 for water, OPLS/2016 for methanol, and cross interactions fitted to excess water/methanol properties over the whole range of concentrations, which provide a rather accurate description of water-methanol mixtures. As expected for hydrophobic pores, whereas pure methanol is adsorbed already from the gas phase, pure water only enters the pore at pressures well beyond bulk saturation for all pore sizes considered. When adsorbed from a mixture, however, water adsorbs at much lower pressures due to the formation of hydrogen bonds with previously adsorbed methanol molecules. For all studied compositions and pore sizes, methanol adsorbs preferentially over water at liquid-vapor equilibrium conditions. In pure components, both water and methanol are microscopically structured in layers, the number of layers increasing with pore size. This is also the case in adsorbed mixtures, in which methanol has a higher affinity for the walls. This becomes more evident as the pore widens. Diffusion of pure water is higher than that of pure methanol for all pore sizes due to the larger size of the methyl group. In mixtures, both components present similar diffusivities at all pore sizes, which is explained in terms of the coupling of molecular movements due to strong hydrogen bonding between methanol and water molecules. This is particularly evident in very narrow pores, in which pure methanol diffusion is completely impeded on the time scale of our simulations, but the presence of a small amount of water molecules facilitates alcohol diffusion following a single-file mechanism. Additionally, our results indicate that pure water diffusivities display a non-monotonous dependence of pore size, due to effects of confinement (proximity to a fluid-solid-fluid transition induced by confinement as reported in previous work) and the dynamic anomalies of water.

Temperature Effects on Water-Mediated Interactions at the Nanoscale

We perform molecular dynamics simulations to study the effects of temperature on the water-mediated interactions between nanoscale apolar solutes. Specifically, we calculate the potential of mean force (PMF) between two graphene plates immersed in water at 240 ≤ T ≤ 400 K and P = 0.1 MPa. These are thermodynamic conditions relevant to cooling- and heating-induced protein denaturation. It is found that both cooling and heating tend to suppress the attraction, and ultimate collapse, of the graphene plates. However, the underlying role played by water upon heating and cooling is different. Isobaric heating reduces the strength and range of the interactions between the plates. Instead, isobaric cooling stabilizes the plates separations that can accommodate an integer number of water layers between the graphene plates. In particular, the energy barriers separating these plate separations increase linearly with 1/ T. We also explore the sensitivity of the plates PMF to the water model employed. In the case of the TIP4P/2005 model, water confined between the plates crystallizes into a defective bilayer ice at low temperatures, whereas in the case of the SPC/E model, water remains in the liquid state at same thermodynamic conditions. The effects of varying water-graphene interactions on the plates PMF are also studied.