Ab initio nanofluidics: disentangling the role of the energy landscape and of density correlations on liquid/solid friction

Elucidating the transport of water and ions in the nanochannel of covalent organic frameworks by molecular dynamics

Inspired by biological channels, achieving precise separation of ion/water and ion/ion requires finely tuned pore sizes at molecular dimensions and deliberate exposure of charged groups. Covalent organic frameworks (COFs), a class of porous crystalline materials, offer well-defined nanoscale pores and diverse structures, making them excellent candidates for nanofluidic channels that facilitate ion and water transport. In this study, we perform molecular simulations to investigate the structure and kinetics of water and ions confined within the typical COFs with varied exposure of charged groups. The COFs exhibit vertically arrayed nanochannels, enabling diffusion coefficients of water molecules within COFs to remain within the same order of magnitude as in the bulk. The motion of water molecules manifests in two distinct modes, creating a mobile hydration layer around acid groups. The ion diffusion within COFs displays a notable disparity between monovalent (M+) and divalent (M2+) cations. As a result, the selectivity of M+/M2+ can exceed 100, while differentiation among M+ is less pronounced. In addition, our simulations indicate a high rejection (R > 98%) in COFs, indicating their potential as ideal materials for desalination. The chemical flexibility of COFs indicates that would hold significant promise as candidates for advanced artificial ion channels and separation membranes.

Impart of Heterogeneous Charge Polarization and Distribution on Friction at Water-Graphene Interfaces: a Density-Functional-Theory based Machine Learning Study

Accurately characterizing friction behaviors at water-solid interfaces remains a challenge because of the dynamic nature of water molecules and temporal variations in solid surface charges. By using a density-functional-theory (DFT) based machine learning (ML) technique and long-time ML-parametrized molecular dynamics simulations, we have systematically investigated water-induced charge polarization and redistribution on graphene, as well as its impact on friction at water-graphene interfaces. Heterogeneous charge polarization and distribution are observed for water-covered graphene accompanied by the formation of electric double layers (EDLs). The introduction of defects into graphene significantly enhances the heterogeneity in charge polarization and distribution. Compared to pristine graphene, defected graphene exhibits reduced friction at water-graphene interfaces due to stronger charge heterogeneity, resulting in lower surface charge density and the inverse relationship between slip length and surface charge density for EDLs. Our results highlight the pivotal roles of defects and charge heterogeneity in reducing friction at water-graphene interfaces.

Passive fractionating mechanism for oil spill using shear-wettability modulation

Oil spillage and organic solvent leakage have been a frequent occurrence in recent years, which pose a significant threat not only to the aquatic ecosystems but also result in substantial economic burdens. This has necessitated the search for materials capable of separating oil from water at enhanced efficiency with superior mechanical and thermal properties. In this study, we conduct a set of systematic molecular dynamics simulations to investigate the potential of two-dimensional graphene-like channels under extreme confinement to achieve efficient oil-water separation. Effective modulation of the wetting characteristics of graphene-like surfaces juxtaposed with unconventional flow behavior at the nanoscale unveils differential interaction of water and oil molecules towards the wall, thereby resulting in distinct separation zones for varying compositions of the oil-water mixture. Such separation zones have been observed to be highly correlated with mixture temperature, which provides effective separation pathways across diverse environmental conditions. Our study offers a paradigm shift in oil-water separation strategies, which not only provides deeper insights into the equilibrium and dynamic behavior of a two-phase mixture but also holds immense implications for the development of smart, wettability-based oil separation devices.

Near-room-temperature water-mediated densification of bulk van der Waals materials from their nanosheets

The conventional fabrication of bulk van der Waals (vdW) materials requires a temperature above 1,000 °C to sinter from the corresponding particulates. Here we report the near-room-temperature densification (for example, ∼45 °C for 10 min) of two-dimensional nanosheets to form strong bulk materials with a porosity of <0.1%, which are mechanically stronger than the conventionally made ones. The mechanistic study shows that the water-mediated activation of van der Waals interactions accounts for the strong and dense bulk materials. Initially, water adsorbed on two-dimensional nanosheets lubricates and promotes alignment. The subsequent extrusion closes the gaps between the aligned nanosheets and densifies them into strong bulk materials. Water extrusion also generates stresses that increase with moulding temperature, and too high a temperature causes intersheet misalignment; therefore, a near-room-temperature moulding process is favoured. This technique provides an energy-efficient alternative to design a wide range of dense bulk van der Waals materials with tailored compositions and properties.

Predicting Quantum-Mechanical Partial Charges in Arbitrarily Long Boron Nitride Nanotubes to Accurately Simulate Nanoscale Water Transport

Single-walled boron nitride nanotubes (BNNTs) have been explored for various applications, ranging from water desalination to osmotic power harvesting. However, no simulation work so far has modeled the changes in the partial charge distribution when a flat sheet is rolled into a tube, hindering the ability to perform accurate molecular dynamics (MD) simulations of water flow through BNNTs. To address this knowledge gap, we employ electronic density functional theory (DFT) calculations to precisely estimate quantum-mechanically derived partial charges on boron (B) and nitrogen (N) atoms in BNNTs of varying lengths and diameters. We observe a spatially varying charge distribution inside both armchair and zigzag nanotubes of finite lengths. Performing DFT calculations for longer BNNTs is computationally intractable, even with state-of-the-art computing resources. To solve this issue, we devise a charge assignment scheme to predict partial charges for longer BNNTs using DFT data for shorter nanotubes, thus overcoming the need to perform more expensive DFT calculations. We show that these charges reproduce the electrostatic potential predicted from first-principles simulations. Subsequently, we carried out MD simulations to predict the effect of the charge distribution inside BNNTs on water flow enhancement via them. We find that using uniform charges leads to an underprediction in flow enhancement, as compared to using quantum-mechanical charges for both armchair and zigzag BNNTs. We also incorporate atomic vibrations into our simulations and show that these vibrations lead to a reduction in the water flow through aperiodic BNNTs. Our work demonstrates the requirement of a quantum-mechanical charge assignment scheme for BNNTs and evolves a framework to assign charges to nanotubes of arbitrary length, thus allowing realistic MD simulations of long BNNTs using accurate partial charges.

Accuracy of TIP4P/2005 and SPC/Fw Water Models

Water is used as the main solvent in model systems containing bioorganic molecules. Choosing the right water model is an important step in the study of the biophysical and biochemical processes that occur in cells. In the present work, we perform molecular dynamics simulations using two distinct force fields for water: the rigid model TIP4P/2005, where only intermolecular interactions are considered, and the flexible model SPC/Fw, where intramolecular interactions are also taken into account. The simulations aim to determine the effect of the inclusion of intramolecular interactions on the accuracy of calculated properties of bulk water (density and thermal expansion coefficient, self-diffusion coefficients, shear viscosity, radial distribution functions, and dielectric constant), as compared to experimental results, over a temperature range between 250 and 370 K. We find that the results of the rigid model present the smallest deviations relative to experiments for most of the calculated quantities, except for the shear viscosity of supercooled water and the water dielectric constant, where the flexible model presents better agreement with experiments.

Development of Heteroatomic Constant Potential Method with Application to MXene-Based Supercapacitors

We describe a method for modeling constant-potential charges in heteroatomic electrodes, keeping pace with the increasing complexity of electrode composition and nanostructure in electrochemical research. The proposed "heteroatomic constant potential method" (HCPM) uses minimal added parameters to handle differing electronegativities and chemical hardnesses of different elements, which we fit to density functional theory (DFT) partial charge predictions in this paper by using derivative-free optimization. To demonstrate the model, we performed molecular dynamics simulations using both HCPM and conventional constant potential method (CPM) for MXene electrodes with Li-TFSI/AN (lithium bis(trifluoromethane sulfonyl)imide/acetonitrile)-based solvent-in-salt electrolytes. Although the two methods show similar accumulated charge storage on the electrodes, the results indicated that HCPM provides a more reliable depiction of electrode atom charge distribution and charge response compared with CPM, accompanied by increased cationic attraction to the MXene surface. These results highlight the influence of elemental composition on electrode performance, and the flexibility of our HCPM opens up new avenues for studying the performance of diverse heteroatomic electrodes including other types of MXenes, two-dimensional materials, metal-organic frameworks (MOFs), and doped carbonaceous electrodes.

Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models

In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.

Electron cooling in graphene enhanced by plasmon-hydron resonance

Evidence is accumulating for the crucial role of a solid's free electrons in the dynamics of solid-liquid interfaces. Liquids induce electronic polarization and drive electric currents as they flow; electronic excitations, in turn, participate in hydrodynamic friction. Yet, the underlying solid-liquid interactions have been lacking a direct experimental probe. Here we study the energy transfer across liquid-graphene interfaces using ultrafast spectroscopy. The graphene electrons are heated up quasi-instantaneously by a visible excitation pulse, and the time evolution of the electronic temperature is then monitored with a terahertz pulse. We observe that water accelerates the cooling of the graphene electrons, whereas other polar liquids leave the cooling dynamics largely unaffected. A quantum theory of solid-liquid heat transfer accounts for the water-specific cooling enhancement through a resonance between the graphene surface plasmon mode and the so-called hydrons-water charge fluctuations-particularly the water libration modes, which allows for efficient energy transfer. Our results provide direct experimental evidence of a solid-liquid interaction mediated by collective modes and support the theoretically proposed mechanism for quantum friction. They further reveal a particularly large thermal boundary conductance for the water-graphene interface and suggest strategies for enhancing the thermal conductivity in graphene-based nanostructures.

A numerical model for water hydration on nanosurfaces: from friction to hydrophilicity and hydrophobicity

Fluidic transport down to the nanometer scale is of great importance for a wide range of applications such as energy harvesting, seawater desalination, and water treatment and may help to understand many biological processes. In this work, we studied the interfacial friction of liquid water on a series of nanostructures through molecular dynamics (MD) simulations. Our results reveal that the friction coefficient of the water-solid interface cannot be described using a previously reported simple function of the free energy corrugation. Considering that the water-solid friction is firmly correlated with the microscopic water motion, we proposed a probability parameter P(d, t) to classify water motion modes on a surface. We demonstrate that this parameter can be used to accurately predict the water-solid friction by simply monitoring the water binding time on a nanosurface. More importantly, according to the relationship between P(d, t) and friction, we found that the friction coefficient can be used as an indicative criterion for quantitatively assessing hydrophobic or hydrophilic materials, where the borderline is roughly 2 × 10⁵ N s m^-3. That is if the water-solid friction is less than 2 × 10⁵ N s m^-3, the surface is considered hydrophobic. But if the friction is larger than this value, the surface is hydrophilic. The present findings could help to better understand fluidic transport at the nanoscale and guide the future design of functional materials, such as super-hydrophobic and super-hydrophilic surfaces by structure engineering.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Classical Quantum Friction at Water-Carbon Interfaces

Friction at water-carbon interfaces remains a major puzzle with theories and simulations unable to explain experimental trends in nanoscale waterflow. A recent theoretical framework─quantum friction (QF)─proposes to resolve these experimental observations by considering nonadiabatic coupling between dielectric fluctuations in water and graphitic surfaces. Here, using a classical model that enables fine-tuning of the solid's dielectric spectrum, we provide evidence from simulations in general support of QF. In particular, as features in the solid's dielectric spectrum begin to overlap with water's librational and Debye modes, we find an increase in friction in line with that proposed by QF. At the microscopic level, we find that this contribution to friction manifests more distinctly in the dynamics of the solid's charge density than that of water. Our findings suggest that experimental signatures of QF may be more pronounced in the solid's response rather than liquid water's.

Co-Ion Desorption as the Main Charging Mechanism in Metallic 1T-MoS2 Supercapacitors

Metallic 1T-MoS₂ is a promising electrode material for supercapacitor applications. Its layered structure allows the efficient intercalation of ions, leading to experimental volumetric capacitance as high as 140 F/cm³. Molecular dynamics could in principle be used to characterize its charging mechanism; however, unlike conventional nanoporous carbon, 1T-MoS₂ is a multicomponent electrode. The Mo and S atoms have very different electronegativities so that 1T-MoS₂ cannot be simulated accurately using the conventional constant potential method. In this work, we show that controlling the electrochemical potential of the atoms allows one to recover average partial charges for the elements in agreement with electronic structure calculations for the material at rest, without compromising the ability to simulate systems under an applied voltage. The simulations yield volumetric capacitances in agreement with experiments. We show that due to the large electronegativity of S, the co-ion desorption is the main charging mechanism at play during the charging process. This contrasts drastically with carbon materials for which ion exchange and counterion adsorption usually dominate. In the future, our method can be extended to the study of a wide range of families of 2D layered materials such as MXenes.

Water Flow in Single-Wall Nanotubes: Oxygen Makes It Slip, Hydrogen Makes It Stick

Experimental measurements have reported ultrafast and radius-dependent water transport in carbon nanotubes which are absent in boron nitride nanotubes. Despite considerable effort, the origin of this contrasting (and fascinating) behavior is not understood. Here, with the aid of machine learning-based molecular dynamics simulations that deliver first-principles accuracy, we investigate water transport in single-wall carbon and boron nitride nanotubes. Our simulations reveal a large, radius-dependent hydrodynamic slippage on both materials, with water experiencing indeed a ≈5 times lower friction on carbon surfaces compared to boron nitride. Analysis of the diffusion mechanisms across the two materials reveals that the fast water transport on carbon is governed by facile oxygen motion, whereas the higher friction on boron nitride arises from specific hydrogen-nitrogen interactions. This work not only delivers a clear reference of quantum mechanical accuracy for water flow in single-wall nanotubes but also provides detailed mechanistic insight into its radius and material dependence for future technological application.

Connection between water's dynamical and structural properties: Insights from ab initio simulations

SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.

Dielectric response of thin water films: a thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent's contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid's boundary, even for very thin (∼1 nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

Free energy calculations from molecular simulations reveal that water's interfacial dielectric response is well-described by bulk properties.

Osmotic Transport at the Aqueous Graphene and hBN Interfaces: Scaling Laws from a Unified, First-Principles Description

Osmotic transport in nanoconfined aqueous electrolytes provides alternative venues for water desalination and "blue energy" harvesting. The osmotic response of nanofluidic systems is controlled by the interfacial structure of water and electrolyte solutions in the so-called electrical double layer (EDL), but a molecular-level picture of the EDL is to a large extent still lacking. Particularly, the role of the electronic structure has not been considered in the description of electrolyte/surface interactions. Here, we report enhanced sampling simulations based on ab initio molecular dynamics, aiming at unravelling the free energy of prototypical ions adsorbed at the aqueous graphene and hBN interfaces, and its consequences on nanofluidic osmotic transport. Specifically, we predicted the zeta potential, the diffusio-osmotic mobility, and the diffusio-osmotic conductivity for a wide range of salt concentrations from the ab initio water and ion spatial distributions through an analytical framework based on Stokes equation and a modified Poisson-Boltzmann equation. We observed concentration-dependent scaling laws, together with dramatic differences in osmotic transport between the two interfaces, including diffusio-osmotic flow and current reversal on hBN but not on graphene. We could rationalize the results for the three osmotic responses with a simple model based on characteristic length scales for ion and water adsorption at the surface, which are quite different on graphene and on hBN. Our work provides fundamental insights into the structure and osmotic transport of aqueous electrolytes on 2D materials and explores alternative pathways for efficient water desalination and osmotic energy conversion.

Distinct Chemistries Explain Decoupling of Slip and Wettability in Atomically Smooth Aqueous Interfaces

Despite essentially identical crystallography and equilibrium structuring of water, nanoscopic channels composed of hexagonal boron nitride and graphite exhibit an order-of-magnitude difference in fluid slip. We investigate this difference using molecular dynamics simulations, demonstrating that its origin is in the distinct chemistries of the two materials. In particular, the presence of polar bonds in hexagonal boron nitride, absent in graphite, leads to Coulombic interactions between the polar water molecules and the wall. We demonstrate that this interaction is manifested in a large typical lateral force experienced by a layer of oriented hydrogen atoms in the vicinity of the wall, leading to the enhanced friction in hexagonal boron nitride. The fluid adhesion to the wall is dominated by dispersive forces in both materials, leading to similar wettabilities. Our results rationalize recent observations that the difference in frictional characteristics of graphite and hexagonal boron nitride cannot be explained on the basis of the minor differences in their wettabilities.

Anomalously low friction of confined monolayer water with a quadrilateral structure

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

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.

Liquid-Solid Slip on Charged Walls: The Dramatic Impact of Charge Distribution

Nanofluidic systems show great promise for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counterions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight into the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g., for blue energy or waste heat harvesting.