Enhanced fluidity of water in superhydrophobic nanotubes: estimating viscosity using jump-corrected confined Stokes-Einstein approach

Accurately predicting the viscosity of water confined within nanotubes is vital for various technological applications. Traditional methods have failed in this regard, necessitating a novel approach. We introduced the jump-corrected confined Stokes-Einstein (JCSE) method and now employ the same to estimate the viscosity and diffusion in superhydrophobic nanotubes. Our study covers a temperature range of 230-300 K and considers three nanotube diameters. Results show that water inside superhydrophobic nanotubes exhibits a significantly lower viscosity and higher diffusion than those inside hydrophobic nanotubes. Narrower nanotubes and lower temperatures accentuate these effects. Furthermore, water inside superhydrophobic nanotubes display a lower viscosity than bulk water, with the difference increasing at lower temperatures. This reduction is attributed to weaker water-water interactions caused by a lower water density in the interfacial region. These findings highlight the importance of interfacial water density and its influence on nanotube viscosity, shedding light on nanoscale fluid dynamics and opening avenues for diverse applications.

Modeling water transport properties in carbon nanotubes: Interplay between force-field flexibility and geometrical parameters

Modeling water and other liquids in computational simulations requires a large set of parameters. Many works have been devoted to finding new, improved water models, with almost all of them designed for bulk systems. Here, we use carbon nanotubes as a play model to investigate the effects of introducing flexibility in water force fields during molecular dynamics simulations of nanoconfined water. We explore six different models to show that viscosity, diffusion, and dipole orientation are vastly influenced by the flexibility and the family of force fields used. Particularly, we found the level of confinement (decreasing the nanotube's diameter) to increase discrepancies in the description of the dipole alignment. In smaller (10,10) nanotubes, the flexible version of the transferable intermolecular potential with three points (TIP3P/Fs) features a high directionality, while its rigid counterpart shows a more distributed dipole orientation. Both viscosity and diffusion are also extremely dependent on the force-field family, with the flexible version of the simple point charge (SPC/Fw) featuring the lower confidence interval.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Molecular-Level Insights into Unique Behavior of Water Molecules Confined in the Heterojunction between One- and Two-Dimensional Nanochannels

With the increasing importance of nanoconfined water in various heterostructures, it is quite essential to clarify the influence of nanoconfinement on the unique properties of water molecules in the pivotal heterojunction. In this work, we reported a series of classical molecular dynamics (MD) simulations to explore nanoconfined water in the subnanometer-sized and nanometer-sized heterostructures by adjusting one-dimensional (1-D) carbon nanotubes with different diameters and two-dimensional (2-D) graphene sheets with different interlayer distances. Our simulation results demonstrated that water molecules in the 1-D/2-D heterojunction show an obvious structural rearrangement associated with the remarkable breaking and formation of hydrogen bonds (HBs), and such rearrangements in the subnanometer-sized systems are much more pronounced than those in the nanometer-sized ones. When water molecules in the 1-D/2-D heterojunctions migrate from 2-D to 1-D confinements, the ordered multi-layer structure in the 2-D confinement are completely destroyed and then transform into different circular HB networks near the nanotube orifice for better connecting to the single-file or helical HB network in the 1-D nanotubes. Furthermore, water molecules in the 1-D/2-D heterojunctions can form stronger HBs with those water molecules further away from the 1-D confinement, leading to an asymmetrical orientational distribution near the orifice. More importantly, our comparison results revealed that the 1-D confinement plays a more important role than the 2-D confinement in determining both the structures and dynamics of water molecules in the 1-D/2-D heterojunction.

Translational Jump-Diffusion of Hydroxide Ion in Anion Exchange Membrane: Deciphering the Nature of Vehicular Diffusion

Earlier, ab initio and reactive force-field-based molecular dynamics (MD) simulation studies suggested an overwhelming contribution of the vehicular diffusion in the total diffusion of hydroxide ions rather than structural diffusion. But does the vehicular diffusion occur via small-step displacement? This question is important to have an understanding of the real characteristics of vehicular diffusion. To answer this question, we perform a classical molecular dynamics simulation of a system containing a hydroxide ion exchange membrane polymer and hydroxide ion at different hydration levels and temperatures using the same molecular force field (Dubey, V. Chem. Phys. Lett. 2020, 755, 137802), which successfully captured the microscopic structure and dynamics of the system. We use the translational jump-diffusion approach, used previously in supercooled water for understanding the origin of breakdown of the Stokes-Einstein relation, to calculate the jump-diffusion coefficient of hydroxide ion and water in the anion exchange membrane. We have seen a significant role of hydration level and temperature in the mechanism of vehicular diffusion. In overhydrated membrane, both hydroxide ions and water molecules diffuse via both small- and large-step displacement. With decreasing hydration level and temperature, the diffusion is increasingly governed by the jump-diffusion mechanism. The larger contribution of jump-diffusion comes from the stronger caging of the diffusing species by the solvent at lower hydration levels and temperature. These results, therefore, suggest that the hydration level and temperature of the hydroxide ion exchange membrane determine the detailed mechanism of the vehicular diffusion of hydroxide ion, especially whether the diffusion follows hydrodynamics or not.

Rheological Characteristics of Ionic Liquids under Nanoconfinement

While the dynamic properties of ionic liquids (ILs) in nanoconfinement play a crucial role in the performance of IL-based electrochemical and mechanical devices, experimental work mostly falls short at reporting "solid-like" versus "liquid-like" behavior of confined ILs. The present work is the first to conduct frequency-sweep oscillatory-shear rheology on IL nanofilms, reconciling the solid-versus-liquid debate and revealing the importance of shear rate in the behavior. We disentangle and analyze the viscoelasticity of nanoconfined ILs and shed light on their relaxation mechanisms. Furthermore, a master curve describes the scaling of the dynamic behavior of four (non-hydrogen-bonding) ILs under nanoconfinement and reveals the role of the compressibility of the flow units.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations

Nanofluids-dispersions of nanometer-sized particles in a liquid medium-have been proposed for a wide variety of thermal management applications. It is known that a solid-like nanolayer of liquid of typical thicknesses of 0.5-1 nm surrounding the colloidal nanoparticles can act as a thermal bridge between the nanoparticle and the bulk liquid. Yet, its effect on the nanofluid viscosity has not been elucidated so far. In this article, we compute the local viscosity of the nanolayer using equilibrium molecular dynamics based on the Green-Kubo formula. We first assess the validity of the method to predict the viscosity locally. We apply this methodology to the calculation of the local viscosity in the immediate vicinity of a metallic nanoparticle for a wide range of solid-liquid interaction strength, where a nanolayer of thickness 1 nm is observed as a result of the interaction with the nanoparticle. The viscosity of the nanolayer, which is found to be higher than its corresponding bulk value, is directly dependent on the solid-liquid interaction strength. We discuss the origin of this viscosity enhancement and show that the liquid density increment alone cannot explain the values of the viscosity observed. Rather, we suggest that the solid-like structure of the distribution of the liquid atoms in the vicinity of the nanoparticle contributes to the nanolayer viscosity enhancement. Finally, we observe a failure of the Stokes-Einstein relation between viscosity and diffusion close to the wall, depending on the liquid-solid interaction strength, which we rationalize in terms of the hydrodynamic slip.

Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach

A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.

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.

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS₂, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.

Neutron spin echo monitoring of segmental-like diffusion of water confined in the cores of carbon nanotubes

Following the stream of increasing scientific interest in condensed-matter systems under ultra-hydrophobic confinement, the present work reports the first incoherent neutron spin echo assessment of the dynamics of water axially confined inside single-wall carbon nanotubes of diameter d∼ 1.4 nm. At the time scale of nanoseconds, two water populations are retrieved, whose relative proportion matches the one expected for a concentric shell + chain arrangement with cylindrical symmetry. The time dependence of the mean square displacement related to the external component is found to be subdiffusive, with peculiar resemblance to segmental diffusion typical of entangled polymeric systems.

Characterization of the non-Arrhenius behavior of supercooled liquids by modeling nonadditive stochastic systems

The characterization of the formation mechanisms of amorphous solids is a large avenue for research, since understanding its non-Arrhenius behavior is challenging to overcome. In this context, we present one path toward modeling the diffusive processes in supercooled liquids near glass transition through a class of nonhomogeneous continuity equations, providing a consistent theoretical basis for the physical interpretation of its non-Arrhenius behavior. More precisely, we obtain the generalized drag and diffusion coefficients that allow us to model a wide range of non-Arrhenius processes. This provides a reliable measurement of the degree of fragility of the system and an estimation of the fragile-to-strong transition in glass-forming liquids, as well as a generalized Stokes-Einstein equation, leading to a better understanding of the classical and quantum effects on the dynamics of nonadditive stochastic systems.

Dynamical aspects of supercooled TIP3P-water in the grooves of DNA

We investigate by molecular dynamics simulations the mobility of the water located at the DNA minor and major grooves. We employ the TIP3P water model, and our system is analyzed for a range of temperatures 190-300 K. For high temperatures, the water at the grooves shows an Arrhenius behavior similar to that observed in the bulk water. At lower temperatures, a departure from the bulk behavior is observed. This slowing down in the dynamics is compared with the dynamics of the hydrogen of the DNA at the grooves and with the autocorrelation functions of the water hydrogen bonds. Our results indicate that the hydrogen bonds of the water at the minor grooves are highly correlated, which suggests that this is the mechanism for the slow dynamics at this high confinement.

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.

Fluidity and phase transitions of water in hydrophobic and hydrophilic nanotubes

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

Breakdown of continuum model for water transport and desalination through ultrathin graphene nanopores: insights from molecular dynamics simulations

In the quest for identifying a graphene membrane for efficient water desalination, molecular dynamics simulations were performed for the pressure-driven flow of salty water across a multilayer graphene membrane.

Single-file mobility of water-like fluid in a generalized Frenkel-Kontorova model

In this work, we used a generalized Frenkel-Kontorova model to study the mobility of water molecules inside carbon nanotubes with small radius at low temperatures. Our simulations show that the mobility of confined water decreases monotonically increasing the amplitude of the substrate potential at fixed commensurations. On the other hand, the mobility of the water molecules shows a non-monotonic behavior when varying the commensuration. This result indicates that the mobility of the confined fluid presents different behavior regimes depending on the amplitude of the water-nanotube interaction. In order to qualitatively understand these results, we study analytically the driven Frenkel-Kontorova model at finite temperatures. This analysis allows us to obtain the curves of the mobility versus commensurations, at fixed substrate potentials. Such curves show the existence of three regimes of mobility behavior as a function of the commensuration ratio. Additionally, our study indicates a nontrivial and strong dependence of the mobility with a quantity that can be interpreted as an effective amplitude of the substrate potential, depending on the bare amplitude of the substrate potential, the commensuration ratio, and temperature.