Structure of water confined between two parallel graphene plates

Collective Vibration Decoupling of Confined Water in Membrane Channels

We previously reported that the asymmetric IR absorption of monolayer water confined within two-dimensional nanochannels is capable of nonthermally inducing a unidirectional flow [Zhang, Q. L. Phys. Rev. Lett. 2024, 132, 184003], while the reason for the difference in the collective vibration IR spectrum between the confined water (CW) and bulk water is still not fully understood. Here, using molecular dynamics simulations, we systematically demonstrated that the CW in narrow graphene membrane channels will appear as a predominant fingerprint-peak and a subpeak in the collective vibration spectrum band. A comparison with the calculated IR spectrum for the CW in the channels with different interlayer spacings revealed that the double-peaked pattern originates from the decoupling of the CW's collective vibration. The highlight spectral intensity of the fingerprint-peak is attributed to the low-cost out-of-plane vibration (wag mode) of the CW molecules. These findings help us understand the physical origins of the unique IR spectra of CW in nanochannels, thereby providing a robust theoretical support for the regulation of the CW's structure and dynamics properties by a remote terahertz stimulation.

Microscopic insights into the effects of interfacial dynamics and nanoconfinement on characteristics of calcium carbonate clusters within two-dimensional nanochannels

Herein, the interfacial effects on calcium carbonate clustering within two-dimensional (2D) graphene nanochannels were systematically investigated using molecular dynamics simulations. The distribution characteristics of the ions at the interface can be attributed to the ordered water layers within the 2D nanochannels. The orientation of CO32- is approximately perpendicular to the interface, which can be attributed to hydrogen bonding and its association with Ca2+ at the interface region. The results show that characteristics of CaCO3 clusters can be affected by ion dynamics at the interface and nanoconfinement, although they prefer to locate in the bulk-like region. Due to nanoconfinement, ion dynamics are slowed down, especially in the direction perpendicular to the graphene surface. Due to the distribution and orientation characteristics of CO32- in the interface region, particularly considering the hydration dynamics of Ca2+ and CO32-, the association between Ca2+ and CO32- ions in CaCO3 clusters at the interface can be promoted as Ca2+ moves from the interface region to the bulk-like region. The ion dynamics and coordination characteristics of CaCO3 near the interface region within 2D nanochannels facilitate the formation of CaCO3 clusters with highly coordinated Ca2+-CO32- structures, which might favor the nucleation of aragonite. The results provide insight into the effects of nanoconfinement and interfacial water layers on biomineral nucleation and offer theoretical insights into the new preparation methods of novel inorganic functional materials.

Ultrahigh-Flux Water Nanopumps Generated by Asymmetric Terahertz Absorption

Controlling active transport of water through membrane channels is essential for advanced nanofluidic devices. Despite advancements in water nanopump design using techniques like short-range invasion and subnanometer-level control, challenges remain facilely and remotely realizing massive waters active transport. Herein, using molecular dynamic simulations, we propose an ultrahigh-flux nanopump, powered by frequency-specific terahertz stimulation, capable of unidirectionally transporting massive water through asymmetric-wettability membrane channels at room temperature without any external pressure. The key physics behind this terahertz-powered water nanopump is revealed to be the energy flow resulting from the asymmetric optical absorption of water.

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.

A new mechanism of the interfacial water film dominating low ice friction

It is generally accepted that ice is slippery due to an interfacial water film wetting the ice surface. Despite the current progress in research, the mechanism of low ice friction is not clear, and especially little is known about the behavior of this surface water film under shear and how the sheared interfacial water film influences ice friction. In our work, we investigated the ordering and diffusion coefficient of the interfacial water film and the friction of ice sliding on an atomically smooth solid substrate at the atomic level using molecular dynamics simulations. There are two layers of water molecules at the ice-solid interface that exhibit properties very different from bulk ice. The ice-adjacent water layer is ice-like, and the solid-adjacent water layer is liquid-like. This liquid-like layer behaves in the manner of "confined water," with high viscosity while maintaining fluidity, leading to the slipperiness of the ice. Furthermore, we found that the interfacial water exhibits shear thinning behavior, which connects the structure of the interfacial water film to the coefficient of friction of the ice surface. We propose a new ice friction mechanism based on shear thinning that is applicable to this interfacial water film structure.

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.

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.

Homogenization of Continuum-Scale Transport Properties from Molecular Dynamics Simulations: An Application to Aqueous-Phase Methane Diffusion in Silicate Channels

Silica-based materials including zeolites are commonly used for wide-ranging applications including separations and catalysis. Substrate transport rates in these materials can significantly influence the efficiency of such applications. Two factors that contribute to transport rates include (1) the porosity of the silicate matrix and (2) nonbonding interactions between the diffusing species and the silicate surface. These contributions generally emerge from disparate length scales, namely, "microscopic" (roughly nanometer-scale) and "macroscopic" (roughly micron-scale), respectively. Here, we develop a simulation framework to estimate the simultaneous impact of these factors on methane mass transport in silicate channels. Specifically, we develop a model of methane transport using homogenization theory to obtain transport parameters valid at length scales of hundreds to thousands of nanometers. These parameters implicitly reflect interactions taking place at fractions of a nanometer. The inputs to the homogenization analysis are data from extensive molecular dynamics simulations that incorporate atomistic-scale interactions, processed to yield local diffusion coefficients and mean force potentials. With this model, we demonstrate how nuances in silicate hydration and silica/methane interactions impact methane diffusion rates in silicate materials, including the effects of silicate surface chemistry such as the presence of silanol groups. The molecular dynamics simulations indicate that methane diffusivity at the silica surface is lower than the bulk-like rates observed at the center of channels of sufficient width. However, potentials of mean force generally evidence attractive methane/silica interactions that enhance diffusion overall. By simultaneously accounting for both of these effects, we show that the effective diffusion coefficient for the nanoporous silicate can be approximately double the value of estimates assuming fully bulk-like behavior in the channel. This study therefore demonstrates the importance of determining diffusion coefficients and potentials of mean force at an atomistic level when estimating transport properties in bulk materials. Importantly, we provide a simple homogenization framework to incorporate these molecular-scale attributes into bulk material transport estimates. This hybrid homogenization/molecular dynamics approach will be of general use for describing small-molecule transport in materials with detailed molecular interactions.

Comparative first principles-based molecular dynamics study of catalytic mechanism and reaction energetics of water oxidation reaction on 2D-surface

The study of the water-splitting process, which can proceed in 2e^- as well as 4e^- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e^- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol^-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.

Temperature induced dynamics of water confined between graphene and MoS2

Water trapped between MoS₂ and graphene assumes a form of ice composed of two planar hexagonal layers with a non-tetrahedral geometry. Additional water does not wet these ice layers but forms three-dimensional droplets. Here, we have investigated the temperature induced dewetting dynamics of the confined ice and water droplets. The ice crystals gradually decrease in size with increasing substrate temperature and completely vanish at about 80 °C. Further heating to 100 °C induces changes in water droplet density, size, and shape through droplet coalescence and dissolution. However, even prolonged annealing at 100 °C does not completely dry the interface. The dewetting dynamics are controlled by the graphene cover. Thicker graphene flakes allow faster water diffusion as a consequence of the reduction of graphene's conformity along the ice crystal's edges, which leaves enough space for water molecules to diffuse along the ice edges and evaporate to the environment through defects in the graphene cover.

Hierarchical thermal transport in nanoconfined water

The structure of nanoconfined fluids is particularly non-uniform owing to the wall interaction, resulting in the distinctive characteristic of thermal transport compared to bulk fluids. We present the molecular simulations on the thermal transport of water confined in nanochannels with a major investigation of its spatial distribution under the effects of wall interaction. The results show that the thermal conductivity of nanoconfined water is inhomogeneous and its layered distribution is very similar to the density profile. The layered thermal conductivity is the coupling result of inhomogeneous density and energy distributions that are generally diametrical, and their contributions to the thermal conductivity compensate with each other. However, the accumulative effect of water molecules is really dominating, resulting in a high thermal conductivity in the high-density layers with the low-energy molecules, and vice versa. Moreover, it is found that the adsorptive and repulsive interactions from solid walls have different roles in the hierarchical thermal transport in nanoconfined water. The adsorptive interaction is only responsible for the layered distribution of thermal conductivity, while the repulsive interaction is responsible for the overall thermal conductivity; accordingly, the thermal conductivity is independent of the strength of water-solid interactions. The identified hierarchical thermal transport in nanoconfined water and its underlying mechanisms have a great significance for the understanding of nanoscale thermal transport and even the mass and energy transport of nanoconfined fluids.

Fast increase of nanofluidic slip in supercooled water: the key role of dynamics

Nanofluidics is an emerging field offering innovative solutions for energy harvesting and desalination. The efficiency of these applications depends strongly on liquid-solid slip, arising from a favorable ratio between viscosity and interfacial friction. Using molecular dynamics simulations, we show that wall slip increases strongly when water is cooled below its melting point. For water on graphene, the slip length is multiplied by up to a factor of five and reaches 230 nm at the lowest simulated temperature, T ∼ 225 K; experiments in nanopores can reach much lower temperatures and could reveal even more drastic changes. The predicted fast increase in water slip can also be detected at supercoolings reached experimentally in bulk water, as well as in droplets flowing on anti-icing surfaces. We explain the anomalous slip behavior in the supercooled regime by a decoupling between viscosity and bulk density relaxation dynamics, and we rationalize the wall-type dependence of the enhancement in terms of interfacial density relaxation dynamics. While providing fundamental insights on the molecular mechanisms of hydrodynamic transport in both interfacial and bulk water in the supercooled regime, this study is relevant to the design of anti-icing surfaces, could help explain the subtle phase and dynamical behaviors of supercooled confined water, and paves the way to explore new behaviors in supercooled nanofluidic systems.

Nanoconfinement-Mediated Water Treatment: From Fundamental to Application

Safe and clean water is of pivotal importance to all living species and the ecosystem on earth. However, the accelerating economy and industrialization of mankind generate water pollutants with much larger quantity and higher complexity than ever before, challenging the efficacy of traditional water treatment technologies. The flourishing researches on nanomaterials and nanotechnologies in the past decade have generated new understandings on many fundamental processes and brought revolutionary upgrades to various traditional technologies in almost all areas, including water treatment. An indispensable step toward the real application of nanomaterials in water treatment is to confine them in large processable substrate to address various inherent issues, such as spontaneous aggregation, difficult operation and potential environmental risks. Strikingly, when the size of the spatial restriction provided by the substrate is on the order of only one or several nanometers, referred to as nanoconfinement, the phase behavior of matter and the energy diagram of a chemical reaction could be utterly changed. Nevertheless, the relationship between such changes under nanoconfinement and their implications for water treatment is rarely elucidated systematically. In this Critical Review, we will briefly summarize the current state-of-the-art of the nanomaterials, as well as the nanoconfined analogues (i.e., nanocomposites) developed for water treatment. Afterward, we will put emphasis on the effects of nanoconfinement from three aspects, that is, on the structure and behavior of water molecules, on the formation (e.g., crystallization) of confined nanomaterials, and on the nanoenabled chemical reactions. For each aspect, we will build the correlation between the nanoconfinement effects and the current studies for water treatment. More importantly, we will make proposals for future studies based on the missing links between some of the nanoconfinement effects and the water treatment technologies. Through this Critical Review, we aim to raise the research attention on using nanoconfinement as a fundamental guide or even tool to advance water treatment technologies.

Nanoconfined Fluids: What Can We Expect from Them?

Nanoconfined fluids (NCFs), which are confined in nanospaces, exhibit distinctive nanoscale effects, including surface effects, small-size effects, quantum effects, and others. The continuous medium hypothesis in fluid mechanics is not valid in this context because of the comparable characteristic length of spaces and molecular mean free path, and accordingly, the classical continuum theories developed for the bulk fluids usually cannot describe the mass and energy transport of NCFs. In this Perspective, we summarize the nanoscale effects on the thermodynamics, mass transport, flow dynamics, heat transfer, phase change, and energy transport of NCFs and highlight the related representative works. The applications of NCFs in the fields of membrane separation, oil and gas production, energy harvesting and storage, and biological engineering are especially indicated. Currently, the theoretical description framework of NCFs is still missing, and it is expected that this framework can be established by adopting the classical continuum theories with the consideration of nanoscale effects.

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).

Frenkel line crossover of confined supercritical fluids

We investigate the temperature evolution of dynamics and structure of partially confined Lennard Jones (LJ) fluids in supercritical phase along an isobaric line in the P-T phase diagram using molecular dynamics simulations. We compare the Frenkel line (FL) crossover features of partially confined LJ fluids to that of the bulk LJ fluids in supercritical phase. Five different spacings have been chosen in this study and the FL crossover characteristics have been monitored for each of these spacings for temperatures ranging from 240 K to 1500 K keeping the pressure fixed at 5000 bar. We characterize the FL crossover using density of states (DoS) function and find that partially confined supercritical fluids (SCF) exhibit a progressive shift of FL crossover point to higher temperatures for smaller spacings. While the DoS perpendicular to the walls shows persistent oscillatory modes, the parallel component exhibits a smooth crossover from an oscillatory to non-oscillatory characteristics representative of FL crossover. We find that the vanishing of peaks in DoS parallel to the walls indicates that the SCF no longer supports shear mode excitations and could serve as an identifier of the FL crossover for confined systems just as is done for the bulk. Layer heights of density profiles, self-diffusivity and the peak heights of radial distribution function parallel to the walls also feature the FL crossover consistent with the DoS criteria. Surprisingly, self-diffusivity undergoes an Arrhenius to super-Arrhenius crossover at low temperatures for smaller spacings as a result of enhanced structural order evidenced via pair-excess entropy. This feature, typical of glass-forming liquids and binary supercooled liquids, is found to develop from the glass-like characteristic slowdown and strong caging in confined supercritical fluid, evidenced via mean squared displacement and velocity autocorrelation function respectively, over intermediate timescales.

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.

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.

Understanding the effect of nanoconfinement on the structure of water hydrogen bond networks

Dynamic hydrogen bond trails in water confined between two phospholipid membranes traced by the information flow model.