Electromelting of confined monolayer ice

Electric-field-induced structure and dynamics of an ethanol-water mixture in hydrophobic-hydrophilic nanochannels

Understanding the behavior of confined fluid mixtures under external electric fields is essential for advancing scientific knowledge and improving a wide range of technological applications, from energy systems to biological processes. An electric field has been widely used to investigate the phase transition of water and modification of interfacial water at the nanoscale. However, the structure and dynamics of the interfacial layer in complex confined fluid mixtures, such as ethanol-water mixtures, remain unexplored under the influence of an electric field. In this study, we explore the structural and dynamic behavior of binary ethanol-water mixtures confined within slit-like hydrophilic (mica) and hydrophobic (graphene) nanochannels under an external electric field using classical molecular dynamics (MD) simulations. We find two distinct interfacial water layers near the hydrophilic mica surface, and a more pronounced sharp peak appears near the hydrophobic graphene sheet with increasing electric field. The density maxima of the -OH and -CH3 groups of ethanol shift towards and away from the graphene surface with an increasing electric field. Our simulations reveal that the electric field strongly impacts the inter and intralayer hydrogen bonding among water and ethanol molecules. The diffusion coefficient of water slightly increases with the electric field and then reduces with an electric field for a lower concentration of ethanol. This finding reveals that the electric field influences the desorption of interfacial water near the hydrophilic mica surface, which can be an implication for diverse technological applications like modifying surface wettability.

Characterizing the Orderliness of Interfacial Water through Stretching Vibrations

Spatial orderliness, which is the orderly structure of molecules, differs significantly between interfacial water and bulk water. Understanding this property is essential for various applications in both natural and engineered environments. However, the subnanometer thickness of interfacial water presents challenges for direct and rapid characterization of its structural orderliness. Herein, through molecular dynamics simulations and infrared spectral analysis of interfacial water in a graphene slit pore, we reveal a hyperbolic tangent relationship between the water ordering and its O-H stretching information in the infrared spectrum. Specifically, O-H symmetric stretching dominated in the highly ordered water structure, while a transition to the asymmetric stretching corresponded to an increase in the degree of disorder. Thus, the O-H stretching behavior could serve as a useful and quick assessment of the orderliness of interfacial water. This work provided insights into interfacial water's unique molecular network and structural dynamics and identified the stretching vibrations' key role in its degree of order, providing insight for fields such as nanotechnology, biology, and material science.

Nanofluidic sensing inspired by the anomalous water dynamics in electrical angstrom-scale channels

Manipulation of confined water dynamics by voltage keeps great importance for diverse applications. However, limitations on the membrane functions, voltage-control range, and unclear dynamics need to be addressed. Herein, we report an anomalous electrically controlled gating phenomenon on cation-intercalated multi-layer Ti3C2 membranes and reveal the confined water dynamics. The water permeation rate was improved rapidly following the application and rise of voltage and finally reached a maximum rate at 0.9 V. The permeation rate starts to decrease from 0.9 V. Below 0.9 V, the electric field affects the charge and polarity of water molecules and then leads to ordered and denser rearrangement in the two-dimensional (2D) channel to accelerate the permeation rate. Above 0.9 V, with the assistance of metal cations, the surge in current induced aggregation of water molecules into clusters, thereby limiting the water mobility. Based on these findings, a high-performance humidity sensor was developed by simultaneously optimizing the response and recovery speeds through electric manipulation. This work provides flexible strategies in intelligent membrane design and nanofluidic sensing.

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.

Ferroelectricity of ice nanotube forests grown in three-dimensional graphene: the electric field effect

Generating diverse ferroelectric ice nanotubes (NTs) efficiently has always been challenging, but matters in nanomaterial synthesis and processing technology. In the present work, we propose a method of growing ice NT forests in a single cooling process. A three-dimensional (3D) graphene structure was selected to behave as a representative container in which a batch of (5, 0) ice NTs was formed simultaneously under the cooling process from molecular dynamics simulation. Other similar 3D graphene structures but with different hole configurations, like uniform triangle or both triangle and pentagon, were also tested, revealing that ice NTs with different tube indices, i.e. both (3, 0) and (5, 0), could also be formed at the same time. Intriguingly, the orientations of the dipole moments of the water molecules of an ice NT formed were independent of each other, making the net ferroelectricity of the whole system weakened or even cancelled. An electric field could help change the orientation of the water molecules of the already obtained ice NTs and even twist the tube to be a spiral (5, 1) one if it was applied during the cooling process, such that the net ferroelectricity was greatly improved. The underlying physical mechanism of all phase transition phenomena, including the improvement of the ferroelectricity under an electric field, were explored in depth from the phase transition curves and structural point of view. The obtained results are of significant application value for improving the preparation efficiency of nano-ferroelectric materials, which are prosperous in nano-devices.

Modeling the influence of the external electric fields on water viscosity inside carbon nanotubes

Equilibrium molecular dynamics simulations were performed to explore the effects of external electric fields and confinement on water properties inside various carbon nanotubes (CNTs). Using different GHz electric field frequencies as well as various constant electric field strengths, the radial distribution function and density profile were investigated, by which the impact of the electric fields and confinement on the water structure are revealed. The results indicated water molecules inside the CNT form layered structures due to topological confinement applying external electric fields can disturb this ordered water molecules structure and increase the viscosity of confined water, particularly in the case of CNTs with a radius less than 13.5 Å. Conversely, for CNTs with a radius greater than13.5 Å, the viscosity decreases under the influence of external oscillating or constant electric fields. How dose the synergism of confinement and external electric fields affect the water properties inside the CNTs?

Enhanced Ion Rejection in Carbon Nanotubes by a Lateral Electric Field

Reverse osmosis membranes hold great promise for dealing with global water scarcity. However, the trade-off between ion selectivity and water permeability is a serious obstacle to desalination. Herein, we introduce an effective strategy to enhance the desalination performance of the membrane. A series of molecular dynamics simulations manifest that an additional lateral electric field significantly promotes ion rejection in carbon nanotubes (CNTs) under the drive of longitudinal pressure. Specifically, with the increase in the electric field, the ion flux shows a deep linear decay, while the water flux decreases only slightly, resulting in a linear increase in ion rejection. The energy barriers of ions around the CNT inlet are obtained by calculating the potentials of mean force to explain enhanced ion rejection. The lateral electric field uniformly raises the energy barriers of ions by pushing them away from the CNT inlet, corresponding to the enhanced ion velocity in the field direction. Furthermore, with the increase in CNT diameter, there is a significant increase in the flux of both ions and water; however, the lateral electric field can also obviously enhance the ion rejection in wider CNTs. Consequently, the enhancement of ion rejection by lateral electric fields should be universal for different CNT diameters, which opens a new avenue for selective permeation and may have broad implications for desalination devices with large pore sizes.

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.

Effects of channel size, wall wettability, and electric field strength on ion removal from water in nanochannels

Molecular dynamics simulations are employed to estimate the effect of nanopore size, wall wettability, and the external field strength on successful ion removal from water solutions. It is demonstrated that the presence of ions, along with the additive effect of an external electric field, constitute a multivariate environment that affect fluidic interactions and facilitate, or block, ion drift to the walls. The potential energy is calculated across every channel case investigated, indicating possible ion localization, while electric field lines are presented, to reveal ion routing throughout the channel. The electric field strength is the dominant ion separation factor, while wall wettability strength, which characterizes if the walls are hydrophobic or hydrophilic has not been found to affect ion movement significantly at the scale studied here. Moreover, the diffusion coefficient values along the three dimensions are reported. Diffusion coefficients have shown a decreasing tendency as the external electric field increases, and do not seem to be affected by the degree of wall wettability at the scale investigated here.

Ferroelectric 2D ice under graphene confinement

We here report on the direct observation of ferroelectric properties of water ice in its 2D phase. Upon nanoelectromechanical confinement between two graphene layers, water forms a 2D ice phase at room temperature that exhibits a strong and permanent dipole which depends on the previously applied field, representing clear evidence for ferroelectric ordering. Characterization of this permanent polarization with respect to varying water partial pressure and temperature reveals the importance of forming a monolayer of 2D ice for ferroelectric ordering which agrees with ab-initio and molecular dynamics simulations conducted. The observed robust ferroelectric properties of 2D ice enable novel nanoelectromechanical devices that exhibit memristive properties. A unique bipolar mechanical switching behavior is observed where previous charging history controls the transition voltage between low-resistance and high-resistance state. This advance enables the realization of rugged, non-volatile, mechanical memory exhibiting switching ratios of 10⁶, 4 bit storage capabilities and no degradation after 10,000 switching cycles.

Ferroelectric ordering of water has been at the heart of intense debates due to its importance in enhancing our understanding of the condensed matter. Here, the authors observe ferroelectric properties of water ice in a two dimensional phase under confinement between two graphene layers.

Two-dimensional monolayer salt nanostructures can spontaneously aggregate rather than dissolve in dilute aqueous solutions

It is well known that NaCl salt crystals can easily dissolve in dilute aqueous solutions at room temperature. Herein, we reported the first computational evidence of a novel salt nucleation behavior at room temperature, i.e., the spontaneous formation of two-dimensional (2D) alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement. Microsecond-scale classical molecular dynamics (MD) simulations showed that NaCl or LiCl, initially fully dissolved in confined water, can spontaneously nucleate into 2D monolayer nanostructures with either ordered or disordered morphologies. Notably, the NaCl nanostructures exhibited a 2D crystalline square-unit pattern, whereas the LiCl nanostructures adopted non-crystalline 2D hexagonal ring and/or zigzag chain patterns. These structural patterns appeared to be quite generic, regardless of the water and ion models used in the MD simulations. The generic patterns formed by 2D monolayer NaCl and LiCl nanostructures were also confirmed by ab initio MD simulations. The formation of 2D salt structures in dilute aqueous solution at room temperature is counterintuitive. Free energy calculations indicated that the unexpected spontaneous salt nucleation behavior can be attributed to the nanoscale confinement and strongly compressed hydration shells of ions.

Aqueous solutions under nanoscale confinement exhibit interesting physicochemical properties. This work reports evidence on the spontaneous formation of two-dimensional alkali chloride crystalline/non-crystalline nanostructures in dilute aqueous solution under nanoscale confinement by computer simulations.

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.

Structure, dynamics, and morphology of nanostructured water confined between parallel graphene surfaces and in carbon nanotubes by applying magnetic and electric fields

Water molecules experience certain changes in their properties when they feel an external magnetic or electric field. These changes are significant in different applications, such as biological and biotechnological processes, nano-pumping, and water treatment. In this work, we have performed molecular dynamics (MD) simulations to investigate the different thermodynamics, structure, and dynamics of water molecules confined between two parallel surfaces and also confined in carbon nanotubes (CNTs). We have also applied different electric and magnetic fields in different directions to the confined molecules. In the graphene system, no polygonal shape was formed in either low or high electric fields, whereas rhombic and pentagonal shapes were formed in low and high magnetic fields. In the CNT system, applying electric fields in all three dimensions made the pentagonal shape disappear and the confined water molecules formed a ring shape when the electric field was applied in the axial direction. Applying the electric field perpendicular to the graphene surfaces increases the self-diffusion of the confined molecules, whereas applying the electric and magnetic fields along the CNT axis decreases the self-diffusion of the confined water molecules. In the graphene system, applying the electric field perpendicular to the graphene surfaces decreases the average number of hydrogen bonds (〈HB〉) whereas the magnetic field has little effect on the 〈HB〉. In the CNT system, applying E_x also leads to a smaller number of HBs. Also, applying the magnetic field along the x-direction (along the CNT direction) leads to a greater number of HBs than the other fields.

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.

Phase-dependent friction of nanoconfined water meniscus

A water meniscus naturally forms under ambient conditions at the point of contact between a nanoscale tip and an atomically flat substrate. Here, we study the effect of the phase state of this nanoscale meniscus-consisting of coexisting monolayer, bilayer and trilayer phase domains-on the frictional behavior during tip sliding by means of molecular dynamics simulations. While the meniscus experiences a domain-by-domain liquid-to-solid phase transition induced by lateral compression, we observe an evident transition in measured friction curves from continuous sliding to stick-slip and meanwhile a gradual increase in friction forces. Moreover, the stick-slip friction can be modulated by varying lattice orientation of the monolayer ice domain in the meniscus, choosing the sliding direction or applying in-plane strains to the substrate. Our results shed light on the rational design of high-performance micro- and nano-electromechanical systems relying on hydration lubrication.

Electric field-induced gas dissolving in aqueous solutions

Gas dissolution or accumulation regulating in an aqueous environment is important but difficult in various fields. Here, we performed all-atom molecular dynamics simulations to study the dissolution/accumulation of gas molecules in aqueous solutions. It was found that the distribution of gas molecules at the solid-water interface is regulated by the direction of the external electric field. Gas molecules attach and accumulate to the interface with an electric field parallel to the interface, while the gas molecules depart and dissolve into the aqueous solutions with a vertical electric field. The above phenomena can be attributed to the redistribution of water molecules as a result of the change of hydrogen bonds of water molecules at the interface as affected by the electric field. This finding reveals a new mechanism of regulating gas accumulation and dissolution in aqueous solutions and can have tremendous applications in the synthesis of drugs, the design of microfluidic device, and the extraction of natural gas.

Water electrification: Principles and applications

Deep engineering of liquid water by charge and impurity injection, charged support, current flow, hydrophobic confinement, or applying a directional field has becoming increasingly important to the mankind toward overcoming energy and environment crisis. One can mediate the processes or temperatures of molecular evaporation for clean water harvesting, HO bond dissociation for H₂ fuel generation, solidification for living-organism cryopreservation, structure stiffening for bioengineering, etc., with mechanisms being still puzzling. We show that the framework of "hydrogen bonding and electronic dynamics" has substantiated the progress in the fundamental issues and the aimed engineering. The segmental disparity of the coupled hydrogen bond (O:HO or HB with ":" being lone pair of oxygen) resolves their specific-heat curves and turns out a quasisolid phase (QS, bound at -15 and 4 °C). Electrification shows dual functionality that not only aligns, orders, polarizes water molecules but also stretches the O:HO bond. The O:HO segmental cooperative relaxation and polarization shift the QS boundary through Einstein's relation, ΔΘ_Dx ∝ Δω_x, resulting in a gel-like, viscoelastic, and stable supersolid phase with raised melting point T_m and lowered temperatures for vaporization T_V and ice nucleation T_N. The supersolidity and electro structure ordering provide additional forces to reinforce Armstrong's water bridge. QS dispersion and the secondary effect of electrification such as compression define the T_N for Dufour's electro-freezing. The T_V depression, surface stress disruption, and electrostatic attraction raise Asakawa's molecular evaporability. Composition of opposite, compatible fields eases the HO dissociation and soil wetting. Progress evidences not only the essentiality of the coupled O:HO bond theory but also the feasibility of engineering water and solutions by programmed electrification.

Electric field triggered release of gas from a quasi-one-dimensional hydrate in the carbon nanotube

We systematically investigate the effects of an axial electric field on the formation and decomposition of quasi-one-dimensional nitrogen gas hydrates within a single-walled carbon nanotube (SWNT) by using molecular dynamics (MD) simulations. We find that the nitrogen hydrate in the SWNT undergoes a series of structure phase transitions with increasing electric field. Corresponding to the structure transition, the nitrogen gas releases from the carbon nanotube in the electric field range of 1 V nm^-1 to 2 V nm^-1. However, nitrogen molecules are trapped as guest molecules, forming a molecule wire, in the ice nanotube when the electric field is less than 1 V nm^-1 or larger than 2 V nm^-1. Our simulations indicate that the nanotube is an excellent tiny gas tank that can be used to trap gas molecules and control their release triggered sensitively by electric signals. The key to this phenomenon is the change in orientations of water dipoles induced by the electric field, which leads to the structural change in the hydrogen-bonding network and the change in the diffusion coefficient of the water molecules. Our findings here may help understanding the mechanism of the electrorelease of gas from hydrates confined in the nanoscale space.

Bidirectional regulation of configuration of the carbon nanotube containing a water droplet

Carbon nanotube complexes are known for their miraculous mechanical and electronic properties that are crucial for nano-electromechanical systems (MEMS). In this study, through molecular dynamics simulations we found for the first time that the electric field and temperature can be used to co-regulate a reversible change of cross-sectional configuration of single-wall carbon nanotubes (SWCNTs). We showed that the electric field can help induce the collapse of an SWCNT when it contains a water droplet, while the increase of temperature can quickly recover its configuration. This controllable bistability of SWCNTs is promising for the design of nanodevices such as electromechanical switches in NEMS.

Electric Field Induced Dewetting of Hydrophobic Nanocavities at Ambient Temperature

The understanding of water dewetting in nanoporous materials is of great importance in various fields of science and technology. Herein, we report molecular dynamics simulation results of dewetting of water droplet in hydrophobic nanocavities between graphene walls under the influence of electric field. At ambient temperature, the rate of dewetting induced by electric field is significantly large. Whereas, it is a very low rate of dewetting induced by high temperature (423 K) due to the strong interaction of the hydrogen-bonding networks of water droplets in nanocavities. In addition, the electric filed induced formation of a water column has been found in a vacuum chamber. When the electric field is turned off, the water column will transform into a water droplet. Importantly, the results demonstrate that the rate of electric field-induced dewetting increases with growth of the electric field. Overall, our results suggest that electric field may have a great potential application for nanomaterial dewetting.

Field-enhanced selectivity in nanoconfined ionic transport

Fluid transport confined in nanochannels shows ultrafast permeation and highly efficient separation performance. However, the size-controlled selectivity of hydrated ions with a similar valence and size, such as alkali ions, is well below 5. We propose in this work to boost ion selectivity through the interaction with the wall of flow channels, which can be enhanced by applying an external electric field across the channel. Molecular simulations show that for ions diffusing near the walls of a graphene nanochannel, the hydration shells are perturbed, endowing the contrast in ion-wall interactions to modify the ion-specific free energy landscape. The trapping/hopping nature of ion diffusion near the wall leads to the conclusion that the diffusivity depends on the free energy barriers rather than the hydration size. This effect can be magnified by elevating the field strength, yielding more than ∼10-fold enhancement in the diffusivity-specific selectivity. With recent experimental advances in external electric field control and local electric field modulation near the surface, this work demonstrates a possible route to achieve high selectivity of alkali ions in nanofluidics, and explore the molecular structures and dynamics of hydrated ions near a surface.

Graphynes for Water Desalination and Gas Separation

Selective transport of mass through membranes, so-called separation, is fundamental to many industrial applications, e.g., water desalination and gas separation. Graphynes, graphene analogs yet containing intrinsic uniformly distributed pores, are excellent candidates for highly permeable and selective membranes owing to their extreme thinness and high porosity. Graphynes exhibit computationally determined separation performance far beyond experimentally measured values of commercial state-of-the-art polyamide membranes; they also offer advantages over other atomically thin membranes like porous graphene in terms of controllability in pore geometry. Here, recent progress in proof-of-concept computational research into various graphynes for water desalination and gas separation is discussed, and their theoretically predicted outstanding permeability and selectivity are highlighted. Challenges associated with the future development of graphyne-based membranes are further analyzed, concentrating on controlled synthesis of graphyne, maintenance of high structural stability to withstand loading pressures, as well asthe demand for accurate computational characterization of separation performance. Finally, possible directions are discussed to align future efforts in order to push graphynes and other 2D material membranes toward practical separation applications.

Electrolytes under Inhomogeneous Nanoconfinement: Water Structuring-Mediated Local Ion Accumulation

The behaviors of aqueous electrolytes confined in nanoscale spaces impact a broad range of biological processes and industrial applications. Our current microscopic understanding of confined electrolytes relies primarily on ideal model systems featuring homogeneous nanoconfinement. Here, we investigate the structure and dynamics of various electrolytes subject to inhomogeneous nanoconfinement, i.e., confined in two-dimensional nanochannels with gradually varying local channel heights, by means of molecular dynamics simulations. Our results reveal unexpected local ion accumulation in the inhomogeneous space occurring at boundaries between coexisting structured water phases, including trilayer, four-layer, and bulk-like waters. This contrasts markedly with the intuition that hydrated ions are more favorable to weakly confined regimes due to a steric exclusion effect. We further show that the location and intensity of the water structuring-mediated ion accumulation are sensitive to the nanochannel's geometry and surface wettability. The revealed anomalous ion behaviors in inhomogeneous nanoconfinement should help to improve our understanding of the microscopic mechanism underlying the operation of biological ion channels and to develop functional nanofluidic devices.

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.

Separation of water-alcohol mixtures using carbon nanotubes under an electric field

Carbon nanotubes (CNTs) are a promising candidate for separation membranes because of their ability to transport substances at very high flow rates. However, there is a tradeoff between achieving a high selectivity using small pore sizes and the reduction of water flux. Here, using molecular dynamics simulations, we report that CNTs can effectively separate water-methanol mixtures under an electric field. Without an electric field and under piston pressure, both water and methanol flow through a CNT, resulting in no separation effect. In contrast, under an electric field and high piston pressure, CNTs allow selective water permeation while rejecting the permeation of methanol molecules. This separation effect is caused by the ordered structures of water molecules in the CNT. A high filtering effect is observed under the conditions of high methanol concentration in the solution or even with large-diameter CNTs up to 3.39 nm. As long as the ordered structure of water in the CNTs can be maintained, the strong filtering effect can be maintained.

Supersolidity of undercoordinated and hydrating water

Supersolidity of ice, which was proposed in 2013 and intensively verified since then [C. Q. Sun et al., Density, Elasticity, and Stability Anomalies of Water Molecules with Fewer than Four Neighbors, J. Phys. Chem. Lett., 2013, 4, 2565-2570; C. Q. Sun et al., Density and phonon-stiffness anomalies of water and ice in the full temperature range, J. Phys. Chem. Lett., 2013, 4, 3238-3244], refers to the water molecules being polarized by molecular undercoordination, which is associated with the skin of bulk ice, nanobubbles, and nanodroplets (often called confinement), or by the electrostatic field of ions in salt solutions [X. Zhang et al., Mediating relaxation and polarization of hydrogen-bonds in water by NaCl salting and heating, Phys. Chem. Chem. Phys., 2014, 16(45), 24666-24671; C. Q. Sun et al., (H, Li)Br and LiOH solvation bonding dynamics: molecular nonbond interactions and solute extraordinary capabilities, J. Phys. Chem. B, 2018, 122(3), 1228-1238]. From the perspective of hydrogen bond (O:H-O or HB with ":" representing the lone pairs on O2-) cooperative relaxation and polarization, this review features the recent progress and recommends future trends in understanding the bond-electron-phonon correlation in the supersolid phase. Supersolidity is characterized by a shorter and stiffer H-O bond, longer and softer O:H nonbond, deeper O 1s energy band, and longer photoelectron and phonon lifetimes. The supersolid phase is less dense, viscoelastic, and mechanically and thermally more stable. Furthermore, O:H-O bond cooperative relaxation offsets the boundaries of structural phases and increases the melting point while lowering the freezing temperature of ice, which is known as supercooling and superheating.

Why different water models predict different structures under 2D confinement

Experiments of nanoconfined water between graphene sheets at high pressure suggest that it forms a square ice structure (Algara-Siller et al., Nature, 2015, 519, 443). Molecular dynamics (MD) simulations have been used to attempt to recreate this structure, but there have been discrepancies in the structure formed by the confined water depending on the simulation set-up that was employed and particularly on the choice of water model. Here, using classical molecular dynamics simulations, we have systematically investigated the effect that three different water models (SPC/E, TIP4P/2005 and TIP5P) have on the structure of water confined between two rigid graphene sheets with a 0.9 nm separation. We show that the TIP4P/2005 and the TIP5P water models form a hexagonal AA-stacked structure, whereas the SPC/E model forms a rhombic AB-stacked structure. Our work demonstrates that the formation of these structures is driven by differences in the strength of hydrogen bonds predicted by the three water models, and that the nature of the graphene/water interaction only mildly affects the phase diagram. Considering the available experimental data and first-principle simulations we conclude that, among the models tested, the TIP4P/2005 and TIP5P force fields are for now the most reliable when simulating water under confinement. © 2018 Wiley Periodicals, Inc.

Phase transition-like behavior of the water monolayer close to the polarized surface of a nanotube

By molecular dynamics simulations, we have investigated effects of temperature on the dynamical behavior of water layers at the charged surface of a nanotube. The behavior of the first water monolayer at the charged surface is very different from that of bulk water. There are three different temperature regions for the axial diffusion coefficient and they increase in different ways (linearly or exponentially) with temperature. The dipole distribution of water molecules was chosen as the order parameter to analyze the phase transition-like behavior. The simulation results indicate that the transition from ordered water to disordered water is continuous, which has not been found in the bulk counterpart. The mechanism behind the unexpected phenomenon was also investigated.

Electrically controlled water permeation through graphene oxide membranes

Controlled transport of water molecules through membranes and capillaries is important in areas as diverse as water purification and healthcare technologies^1-7. Previous attempts to control water permeation through membranes (mainly polymeric ones) have concentrated on modulating the structure of the membrane and the physicochemical properties of its surface by varying the pH, temperature or ionic strength^3,8. Electrical control over water transport is an attractive alternative; however, theory and simulations^9-14 have often yielded conflicting results, from freezing of water molecules to melting of ice^14-16 under an applied electric field. Here we report electrically controlled water permeation through micrometre-thick graphene oxide membranes^17-21. Such membranes have previously been shown to exhibit ultrafast permeation of water^17,22 and molecular sieving properties^18,21, with the potential for industrial-scale production. To achieve electrical control over water permeation, we create conductive filaments in the graphene oxide membranes via controllable electrical breakdown. The electric field that concentrates around these current-carrying filaments ionizes water molecules inside graphene capillaries within the graphene oxide membranes, which impedes water transport. We thus demonstrate precise control of water permeation, from ultrafast permeation to complete blocking. Our work opens up an avenue for developing smart membrane technologies for artificial biological systems, tissue engineering and filtration.

Anomalous cation diffusion in salt-doped confined bilayer ice

The diffusive dynamics of aqueous electrolyte solutions in nanoconfined spaces has attracted considerable attention due to their potential applications in desalination, biosensors and supercapacitors. Here we show by molecular dynamics simulations that lithium and sodium ions diffuse at a rate at least an order of magnitude higher than that of water molecules when the ions are trapped in an ice bilayer confined between two parallel plates. This novel picture is in sharp contrast to the prevailing view that the diffusion rate of ions is comparable to or even lower than that of water in both bulk and confined solutions. The predicted high ion mobility stems from frequent lateral hopping of ions along the coordination sites inside the hydrogen-bonding network connecting the two water layers of the ice bilayer. This anomalous diffusion should provide new insights into the physics of confined aqueous electrolytes.

Separation of water-ethanol solutions with carbon nanotubes and electric fields

Bioethanol has been used as an alternative energy source for transportation vehicles to reduce the use of fossil fuels. The separation of water-ethanol solutions from fermentation processes is still an important issue in the production of anhydrous ethanol. Using molecular dynamics simulations, we investigate the effect of axial electric fields on the separation of water-ethanol solutions with carbon nanotubes (CNTs). In the absence of an electric field, CNT-ethanol van der Waals interactions allow ethanol to fill the CNTs in preference to water, i.e., a separation effect for ethanol. However, as the CNT diameter increases, this ethanol separation effect significantly decreases owing to a decrease in the strength of the van der Waals interactions. In contrast, under an electric field, the energy of the electrostatic interactions within the water molecule structure induces water molecules to fill the CNTs in preference to ethanol, i.e., a separation effect for water. More importantly, the electrostatic interactions are dependent on the water molecule structure in the CNT instead of the CNT diameter. As a result, the separation effect observed under an electric field does not diminish over a wide CNT diameter range. Moreover, CNTs and electric fields can be used to separate methanol-ethanol solutions too. Under an electric field, methanol preferentially fills CNTs over ethanol in a wide CNT diameter range.

Fundamental interfacial mechanisms underlying electrofreezing

This article reviews the fundamental interfacial mechanisms underlying electrofreezing (promotion of ice nucleation via the application of an electric field). Electrofreezing has been an active research topic for many decades, with applications in food preservation, cryopreservation, cryogenics and ice formation. There is substantial literature detailing experimental and simulations-based studies, which aim to understand the complex mechanisms underlying accelerated ice nucleation in the presence of electric fields and electrical charge. This work provides a critical review of all such studies. It is noted that application-focused studies of electrofreezing are excluded from this review; such studies have been previously reviewed in literature. This review focuses only on fundamental studies, which analyze the physical mechanisms underlying electrofreezing. Topics reviewed include experimental studies on electrofreezing (DC and AC electric fields), pyroelectricity-based control of freezing, molecular dynamics simulations of electrofreezing, and thermodynamics-based explanations of electrofreezing. Overall, it is seen that electrofreezing can enable disruptive advancements in the control of liquid-to-solid phase change, and that our current understanding of the underlying mechanisms can be significantly improved through further studies of various interfacial effects coming into play.

A controllable water signal transistor

We performed molecular dynamics simulations to study the regulating ability of water chains confined in a Y-shaped nanochannel. It was shown that a signal at the molecular level could be controlled by two other charge-induced signals when the water chains were confined in a Y-shaped nanochannel, demonstrating promising applications as water signal transistors in nanosignal systems. The mechanism of a water signal transistor is similar to a signal logic device. This remarkable ability to control the water signal is attributed to the strong dipole-ordering of the water chains in the nanochannel. The controllable water signal process of the Y-shaped nanochannel provides opportunities for future application in the design of molecular-scale signal devices.

Effects of Confinement on the Dielectric Response of Water Extends up to Mesoscale Dimensions

The extent of confinement effects on water is not clear in the literature. While some properties are affected only within a few nanometers from the wall surface, others are affected over long length scales, but the range is not clear. In this work, we have examined the dielectric response of confined water under the influence of external electric fields along with the dipolar fluctuations at equilibrium. The confinement induces a strong anisotropic effect which is evident up to 100 nm channel width, and may extend to macroscopic dimensions. The root-mean-square fluctuations of the total orientational dipole moment in the direction perpendicular to the surfaces is 1 order of magnitude smaller than the value attained in the parallel direction and is independent of the channel width. Consequently, the isotropic condition is unlikely to be recovered until the channel width reaches macroscopic dimensions. Consistent with dipole moment fluctuations, the effect of confinement on the dielectric response also persists up to channel widths considerably beyond 100 nm. When an electric field is applied in the perpendicular direction, the orientational relaxation is 3 orders of magnitude faster than the dipolar relaxation in the parallel direction and independent of temperature.

Structural and configurational properties of nanoconfined monolayer ice from first principles

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

Water in Inhomogeneous Nanoconfinement: Coexistence of Multilayered Liquid and Transition to Ice Nanoribbons

Phase behavior and the associated phase transition of water within inhomogeneous nanoconfinement are investigated using molecular dynamics simulations. The nanoconfinement is constructed by a flat bottom plate and a convex top plate. At 300 K, the confined water can be viewed as a coexistence of monolayer, bilayer, and trilayer liquid domains to accommodate the inhomogeneous confinement. With increasing liquid density, the confined water with uneven layers transforms separately into two-dimensional ice crystals with unchanged layer number and rhombic in-plane symmetry for oxygen atoms. The monolayer water undergoes the transition first into a puckered ice nanoribbon, and the bilayer water transforms into a rhombic ice nanoribbon next, followed by the transition of trilayer water into a trilayer ice nanoribbon. The sequential localized liquid-to-solid transition within the inhomogeneous confinement can also be achieved by gradually decreasing the temperature at low liquid densities. These findings of phase behaviors of water under the inhomogeneous nanoconfinement not only extend the phase diagram of confined water but also have implications for realistic nanofluidic systems and microporous materials.

Water–methanol separation with carbon nanotubes and electric fields

Methanol is used in various applications, such as fuel for transportation vehicles, fuel cells, and in chemical industrial processes. Conventionally, separation of methanol from aqueous solution is by distillation. However, this method consumes a large amount of energy; hence development of a new method is needed. In this work, molecular dynamics simulations are performed to investigate the effect of an electric field on water–methanol separation by carbon nanotubes (CNTs) with diameters of 0.81 to 4.07 nm. Without an electric field, methanol molecules fill the CNTs in preference to water molecules. The preference of methanol to occupy the CNTs over water results in a separation effect. This separation effect is strong for small CNT diameters and significantly decreases with increasing diameter. In contrast, under an electric field, water molecules strongly prefer to occupy the CNTs over methanol molecules, resulting in a separation effect for water. More interestingly, the separation effect for water does not decrease with increasing CNT diameter. Formation of water structures in CNTs induced by an electric field has an important role in the separation of water from methanol.