Phase Transition in Monolayer Water Confined in Janus Nanopore

Does twist angle affect the properties of water confined inside twisted bilayer graphene?

Graphene nanoslit pores are used for nanofluidic devices, such as, in water desalination, ion-selective channels, ionic transistors, sensing, molecular sieving, blue energy harvesting, and protein sequencing. It is a strenuous task to prepare nanofluidic devices, because a small misalignment leads to a significant alteration in various properties of the devices. Here, we focus on the rotational misalignment between two parallel graphene sheets. Using molecular dynamics simulation, we probe the structure and dynamics of monolayer water confined inside graphene nanochannels for a range of commensurate twist angles. With SPC/E and TIP4P/2005 water models, our simulations reveal the independence of the equilibrium number density- n ∼ 13 nm^-2 for SPC/E and n ∼ 11.5 nm^-2 for TIP4P/2005- across twists. Based on the respective densities of the water models, the structure and dielectric constant are invariant of twist angles. The confined water structure at this density shows square ice ordering for SPC/E water only. TIP4P/2005 shows ordering at the vicinity of a critical density (n ∼ 12.5 nm^-2). The average perpendicular dielectric constant of the confined water remains anomalously low (∼2 for SPC/E and ∼6 for TIP4P/2005) for the studied twist angles. We find that the friction coefficient of confined water molecules varies for small twist angles, while becoming independent for twists greater than 5.1°. Our results indicate that a small, angular misalignment will not impair the dielectric properties of monolayer water within a graphene slit-pore, but can significantly influence its dynamics.

Anisotropy of the Proton Kinetic Energy as a Tool for Capturing Structural Transition in Water Confined in a Graphene Nanoslit Pore

The proton dynamics of a 2D water monolayer confined inside a graphene slit pore is studied in Cartesian and molecular frames of reference using molecular dynamics simulations. The vibrational density of states of the proton was calculated versus temperature and was further used to deduce the mean kinetic energy of the hydrogen atoms, Ke(H), in both frames of reference. The directional components of Ke(H) are in good agreement with experimental observations for bulk as well as nanoconfined water. Nonetheless, while in the molecular frame of reference the effect of temperature on the anisotropy ratios of Ke(H) (the ratio between its directional components) are practically invariant between the 2D and 3D cases, those in the Cartesian frame of reference reveal a rather notable reduction across 200 K, indicating the occurrence of an order-disorder transition. This result is further supported by the calculated entropy and enthalpy of the confined water molecules. Overall, it is shown that Ke(H) anisotropy ratios may serve as a valuable order parameter for detecting structural transformations in hydrogen bonds containing molecular systems.

Influence of the extent of hydrophobicity on water organization and dynamics on 2D graphene oxide surfaces

Graphene oxide (GO) nanomaterials are being extensively explored for a wide spectrum of applications, ranging from water desalination to fuel cell applications due to their tunable mechanical, thermal, and electrical...

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

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

Dimensionality dependence of the Kauzmann temperature: A case study using bulk and confined water

The Kauzmann temperature (T_K) of a supercooled liquid is defined as the temperature at which the liquid entropy becomes equal to that of the crystal. The excess entropy, the difference between liquid and crystal entropies, is routinely used as a measure of the configurational entropy, whose vanishing signals the thermodynamic glass transition. The existence of the thermodynamic glass transition is a widely studied subject, and of particular recent interest is the role of dimensionality in determining the presence of a glass transition at a finite temperature. The glass transition in water has been investigated intensely and is challenging as the experimental glass transition appears to occur at a temperature where the metastable liquid is strongly prone to crystallization and is not stable. To understand the dimensionality dependence of the Kauzmann temperature in water, we study computationally bulk water (three-dimensions), water confined in the slit pore of the graphene sheet (two-dimensions), and water confined in the pore of the carbon nanotube of chirality (11,11) having a diameter of 14.9 Å (one-dimension), which is the lowest diameter where amorphous water does not always crystallize into nanotube ice in the supercooled region. Using molecular dynamics simulations, we compute the entropy of water in bulk and under reduced dimensional nanoscale confinement to investigate the variation of the Kauzmann temperature with dimension. We obtain a value of T_K (133 K) for bulk water in good agreement with experiments [136 K (C. A. Angell, Science 319, 582-587 (2008) and K. Amann-Winkel et al., Proc. Natl. Acad. Sci. U. S. A. 110, 17720-17725 (2013)]. However, for confined water, in two-dimensions and one-dimension, we find that there is no finite temperature Kauzmann point (in other words, the Kauzmann temperature is 0 K). Analysis of the fluidicity factor, a measure of anharmonicity in the oscillation of normal modes, reveals that the Kauzmann temperature can also be computed from the difference in the fluidicity factor between amorphous and ice phases.

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.

DNA Translocation through Vertically Stacked 2D Layers of Graphene and Hexagonal Boron Nitride Heterostructure Nanopore

Cost-effective, fast, and reliable DNA sequencing can be enabled by advances in nanopore-based methods, such as the use of atomically thin graphene membranes. However, strong interaction of DNA bases with graphene leads to undesirable effects such as sticking of DNA strands to the membrane surface. While surface functionalization is one way to counter this problem, here, we present another solution based on a heterostructure nanopore system, consisting of a monolayer of graphene and hexagonal boron nitride (hBN) each. Molecular dynamics studies of DNA translocation through this heterostructure nanopore revealed a surprising and crucial influence of the heterostructure layer order in controlling the base specific signal variability. Specifically, the heterostructure with graphene on top of hBN had nearly 3-10× lower signal variability than the one with hBN on top of graphene. Simulations point to the role of differential underside sticking of DNA bases as a possible reason for the observed influence of the layer order. Our studies can guide the development of experimental systems to study and exploit DNA translocation through two-dimensional heterostructure nanopores for single molecule sequencing and sensing applications.

Disentangling the liquid phase exfoliation of two-dimensional materials: an "in silico" perspective

Liquid Phase Exfoliation (LPE) is one of the most successful synthetic roots for the preparation of two-dimensional (2D) materials from their bulk counterparts. In recent years, significant progress has been accomplished for the development and modification of LPE techniques. However, precise identification of the hierarchical steps of the molecular mechanism of LPE remains to some extent elusive. Additionally, the a priori choice of suitable solvents for successful exfoliation and dispersion of various layered materials poses a challenge for both academia and industry. Computational methods, particularly Molecular Dynamics (MD) simulations with classical force-fields have contributed a great deal towards the understanding of the underlying mechanism of LPE, providing remarkable insights into the molecular-level details of the solvent-material interactions at the nanoscale and predicting "good" and "bad" solvents for exfoliation as well as stabilization of the dispersed state. With an intention to build up a unified understanding, in this perspective article, we summarize the recent advancements of molecular simulation techniques employed to decipher the mechanism of LPE, pointing out the key features of molecular interactions and identifying several thermodynamic parameters governing the phenomena. In addition, we outline the necessary characteristics of solvent molecules, essential for their use as "good" solvents towards LPE. Also, we highlight the limitations of simulation methods for the modelling of LPE. We believe that this article will be beneficial for the selection of solvents for the synthesis of novel 2D materials via LPE and will also provide a comprehensive view to computational material scientists towards the development of novel simulation protocols for investigating and analysing such complex molecular events.

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.

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.

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.

Thick Two-Dimensional Water Film Confined between the Atomically Thin Mica Nanosheet and Hydrophilic Substrate

The interesting properties of water molecules confined in a two-dimensional (2D) environment have aroused great attention. However, the study of 2D-confined water at the hydrophilic-hydrophilic interface is largely unexplored due to the lack of appropriate system. In this work, the behavior of water molecules confined between an atomically thin mica nanosheet and a hydrophilic SiO₂/Si substrate was investigated using an atomic force microscope in detail at ambient conditions. The confined water molecules aggregated as droplets when the relative humidity (RH) of the environment was 11%. A large-area 2D water film with a uniform thickness of ∼2 nm was observed when the mica flake was incubated at 33% RH for 1 h before being mechanically exfoliated on a SiO₂/Si substrate. Interestingly, the water film showed ordered edges with a predominant angle of 120°, which was the same with the lattice orientation of the mica nanosheet on top of it. The water film showed a fluidic behavior at the early stage and reached a stable state after 48 h under ambient conditions. The surface properties of the upper mica nanosheet and the underlying substrate played a crucial role in manipulating the behavior of confined water molecules. When the surface of the upper mica nanosheet was modified by Na⁺, Ni²⁺_, and aminopropyltriethoxysilane (APS), only some small water droplets were observed instead of a water film. The surface of the underlying SiO₂/Si substrate was functionalized by hydrophilic APS and hydrophobic octadecyltrimethoxysiliane (OTS). The small water droplets were imaged on a hydrophobic OTS-SiO₂/Si substrate, while the water film with regular edges was maintained on a hydrophilic APS-SiO₂/Si substrate. Our results might provide an alternative molecular view for investigating structures and properties of confined water molecules in 2D environments.