Modulating Water Slip Using Atomic-Scale Defects: Friction on Realistic Hexagonal Boron Nitride Surfaces

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.

Experimental and Theoretical Insights into Interfacial Properties of 2D Materials for Selective Water Transport Membranes: A Critical Review

Interfacial properties, such as wettability and friction, play critical roles in nanofluidics and desalination. Understanding the interfacial properties of two-dimensional (2D) materials is crucial in these applications due to the close interaction between liquids and the solid surface. The most important interfacial properties of a solid surface include the water contact angle, which quantifies the extent of interactions between the surface and water, and the water slip length, which determines how much faster water can flow on the surface beyond the predictions of continuum fluid mechanics. This Review seeks to elucidate the mechanism that governs the interfacial properties of diverse 2D materials, including transition metal dichalcogenides (e.g., MoS2), graphene, and hexagonal boron nitride (hBN). Our work consolidates existing experimental and computational insights into 2D material synthesis and modeling and explores their interfacial properties for desalination. We investigated the capabilities of density functional theory and molecular dynamics simulations in analyzing the interfacial properties of 2D materials. Specifically, we highlight how MD simulations have revolutionized our understanding of these properties, paving the way for their effective application in desalination. This Review of the synthesis and interfacial properties of 2D materials unlocks opportunities for further advancement and optimization in desalination.

Investigation of Tribological Behavior and Lubrication Mechanisms of Zinc Oxide under Poly α-olefin Lubrication Enhanced by the Electric Field

The electric field induces complex effects on the tribological properties of zinc oxide (ZnO) under lubricated conditions, particularly at the nanoscale, where the friction process and mechanism remain unclear. In this paper, the tribological behaviors of ZnO under the lubrication of poly α-olefins (PAO) were investigated by molecular dynamics (MD) simulations with reactive force field (ReaxFF). The results reveal a significant enhancement in the tribological performances of ZnO with the application of the electric field, resulting in a 58.6% reduction in the coefficient of friction (COF) from 0.193 at 0 V/Å to 0.080 at 0.1 V/Å. This improvement can be attributed to the weakening of interfacial interaction, evidenced by a reduction in the number of C-O covalent bonds under the influence of the electric field, along with the formation of an adsorption film due to applied load and shear effects. Notably, the effect of the electric field and applied load extends the impact of interface slip on the tribological performance of ZnO. Overall, this study provides a comprehensive understanding of the impact of the electric field on reducing the friction of ZnO-based structured models, shedding light on explaining their tribological properties and lubrication mechanisms.

Water Electric Field Induced Modulation of the Wetting of Hexagonal Boron Nitride: Insights from Multiscale Modeling of Many-Body Polarization

Understanding the behavior of water contacting two-dimensional materials, such as hexagonal boron nitride (hBN), is important in practical applications, including seawater desalination and energy harvesting. Water, being a polar solvent, can strongly polarize the hBN surface via the electric fields that it generates. However, there is a lack of molecular-level understanding about the role of polarization effects at the hBN/water interface, including its effect on the wetting properties of water. In this study, we develop a theoretical framework that introduces an all-atomistic polarizable force field to accurately model the interactions of water molecules with hBN surfaces. The force field is then utilized to self-consistently describe the water-induced polarization of hBN using the classical Drude oscillator model, including predicting the hBN-water binding energies which are found to be in excellent agreement with diffusion Monte Carlo (DMC) predictions. By carrying out molecular dynamics (MD) simulations, we demonstrate that the polarizable force field yields a water contact angle on multilayered hBN which is in close agreement with the recent experimentally reported values. Conversely, an implicit modeling of the hBN-water polarization energy utilizing a Lennard-Jones (LJ) potential, a commonly utilized approximation in previous MD simulation studies, leads to a considerably lower water contact angle. This difference in the predicted contact angles is attributed to the significant energy-entropy compensation resulting from the incorporation of polarization effects at the hBN-water interface. Our work highlights the importance of self-consistently modeling the hBN-water polarization energy and offers insights into the wetting-related interfacial phenomena of water on polarizable materials.

Anisotropic Interfacial Force Field for Interfaces of Water with Hexagonal Boron Nitride

This study introduces an anisotropic interfacial potential that provides an accurate description of the van der Waals (vdW) interactions between water and hexagonal boron nitride (h-BN) at their interface. Benchmarked against the strongly constrained and appropriately normed functional, the developed force field demonstrates remarkable consistency with reference data sets, including binding energy curves and sliding potential energy surfaces for various configurations involving a water molecule adsorbed atop the h-BN surface. These findings highlight the significant improvement achieved by the developed force field in empirically describing the anisotropic vdW interactions of the water/h-BN heterointerfaces. Utilizing this anisotropic force field, molecular dynamics simulations demonstrate that atomically flat, pristine h-BN exhibits inherent hydrophobicity. However, when atomic-step surface roughness is introduced, the wettability of h-BN undergoes a significant change, leading to a hydrophilic nature. The calculated water contact angle (WCA) for the roughened h-BN surface is approximately 64°, which closely aligns with experimental WCA values ranging from 52° to 67°. These findings indicate the high probability of the presence of atomic steps on the surfaces of the experimental h-BN samples, emphasizing the need for further experimental verification. The development of the anisotropic interfacial force field for accurately describing interactions at the water/h-BN heterointerfaces is a significant advancement in accurately simulating the wettability of two-dimensional (2D) materials, offering a reliable tool for studying the dynamic and transport properties of water at these interfaces, with implications for materials science and nanotechnology.

Liquid-activated quantum emission from pristine hexagonal boron nitride for nanofluidic sensing

Liquids confined down to the atomic scale can show radically new properties. However, only indirect and ensemble measurements operate in such extreme confinement, calling for novel optical approaches that enable direct imaging at the molecular level. Here we harness fluorescence originating from single-photon emitters at the surface of hexagonal boron nitride for molecular imaging and sensing in nanometrically confined liquids. The emission originates from the chemisorption of organic solvent molecules onto native surface defects, revealing single-molecule dynamics at the interface through the spatially correlated activation of neighbouring defects. Emitter spectra further offer a direct readout of the local dielectric properties, unveiling increasing dielectric order under nanometre-scale confinement. Liquid-activated native hexagonal boron nitride defects bridge the gap between solid-state nanophotonics and nanofluidics, opening new avenues for nanoscale sensing and optofluidics.

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.

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

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

Surface Roughness Explains the Observed Water Contact Angle and Slip Length on 2D Hexagonal Boron Nitride

Hexagonal boron nitride (hBN) is a two-dimensional (2D) material that is currently being explored in a number of applications, such as atomically thin coatings, water desalination, and biological sensors. In many of these applications, the hBN surface comes into intimate contact with water. In this work, we investigate the wetting and frictional behavior of realistic 2D hBN surfaces with atomic-scale defects and roughness. We combine density functional theory calculations of the charge distribution inside hBN with free energy calculations using molecular dynamics simulations of the hBN-water interface. We find that the presence of surface roughness, but not that of vacancy defects, leads to remarkable agreement with the experimentally observed water contact angle of 66° on freshly synthesized, uncontaminated hBN. Not only that, the inclusion of surface roughness predicts with exceptional accuracy the experimental water slip length of ∼1 nm on hBN. Our results underscore the importance of considering realistic models of 2D materials with surface roughness while modeling nanomaterial-water interfaces in molecular simulations.

How Grain Boundaries and Interfacial Electrostatic Interactions Modulate Water Desalination via Nanoporous Hexagonal Boron Nitride

To fulfill the increasing demand for drinking water, researchers are currently exploring nanoporous two-dimensional materials, such as hexagonal boron nitride (hBN), as potential desalination membranes. A prominent, yet unsolved challenge is to understand how such membranes will perform in the presence of defects or surface charge in the membrane material. In this work, we study the effect of grain boundaries (GBs) and interfacial electrostatic interactions on the desalination performance of bicrystalline nanoporous hBN using classical molecular dynamics simulations supported by quantum-mechanical density functional theory (DFT) calculations. We investigate three different nanoporous bicrystalline hBN configurations, with symmetric tilt GBs having misorientation angles of 13.2, 21.8, and 32.2°. Using lattice dynamics calculations, we find that grain boundaries alter the areas and shapes of nanopores in bicrystalline hBN, as compared to the nanopores in monocrystalline hBN. We observe that, although bicrystalline nanoporous hBN with a misorientation angle of 13.2° shows an improved water flow rate by ∼30%, it demonstrates reduced Na⁺ ion rejection by ∼6%, as compared to monocrystalline hBN. We also uncover the role of the nanopore shape in water desalination, finding that more elongated pores with smaller sizes (in 21.8- and 32.2°-misoriented bicrystalline hBN) can match water permeation through less elongated pores of slightly larger sizes, with a concomitant ∼3-4% decrease in Na⁺ rejection. Simulations also predict that the water flow rate is significantly affected by interfacial electrostatic interactions. Indeed, the water flow rate is the highest when altered partial charges on B and N atoms were determined using DFT calculations, as compared to when no partial charges or bulk partial charges (i.e., charged hBN) were considered. Overall, our work on water/ion transport through nanopores in bicrystalline hBN indicates that the presence of GBs and surface charge can lead, respectively, to a decrease in the ion rejection and water permeation performance of hBN membranes.