Wetting and interfacial properties of water nanodroplets in contact with graphene and monolayer boron-nitride sheets

Modulating the Water Contact Angle Using Surface Roughness: Interfacial Properties of Hexagonal Boron Nitride Surfaces

Hexagonal boron nitride (hBN) exhibits immense potential in H2O-related technologies, but its interaction with H2O, especially on rough surfaces, remains unclear. This study unravels the influence of surface roughness and force field selection on hBN wettability using molecular dynamics (MD) simulations. We leverage quantum mechanical calculations to accurately capture the hBN surface charge distribution and combine it with free energy calculations via MD simulations for the hBN-H2O interfaces. Incorporating surface roughness into the model yields results in close agreement with the experimental contact angle of 66° for H2O using FF-3 force fields, validating the simulation approach. However, this approach can yield an unrealistic water contact angle (WCA) of 0° for FF-2 force fields, highlighting the crucial role of force field selection and realistic surface representations. We further dissect the impact of roughness on the WCA, identifying the individual contributions of electrostatic and Lennard-Jones interactions to the work of adhesion. This research investigates the combined impact of surface roughness and force fields on interfacial properties, providing new possibilities for the advancement and optimization of desalination.

1D and 2D Boron Nitride Nano Structures: A Critical Analysis for Emerging Applications in the Field of Nanocomposites

Boron nitride (BN) with its 1D and 2D nano derivatives have gained immense popularity in both the field of research and applications. These nano derivatives have proved to be one of the most promising fillers which can be incorporated in polymers to form nanocomposites with excellent properties. These materials have been around for 25 years whereas significant research has been done in this field for only the past decade. There are many interesting properties which are imparted to the nanocomposites wherein thermal stability, large energy band gap, resistance to oxidation, excellent thermal conductivity, chemical inertness, and exceptional mechanical properties are just a few worthy of mention. Hexagonal boron nitride (h-BN) was selected as the parent material by most researchers reviewed in this paper through which 2D derivative Boron nitride nanosheets (BNNS) and 1D derivative Boron nitride nanotubes (BNNTs) are synthesized. This review will focus on the in-depth properties of h-BN and further will concisely focus on BNNS and BNNTs for their various properties. A detailed discussion of the addition of BNNS and BNNTs into polymers to form nanocomposites, their synthesis, properties, and applications is followed by a summary determining the most suitable synthesizing processes and the materials, keeping in mind the current challenges.

Probing structural superlubricity of two-dimensional water transport with atomic resolution

Low-dimensional water transport can be drastically enhanced under atomic-scale confinement. However, its microscopic origin is still under debate. In this work, we directly imaged the atomic structure and transport of two-dimensional water islands on graphene and hexagonal boron nitride surfaces using qPlus-based atomic force microscopy. The lattice of the water island was incommensurate with the graphene surface but commensurate with the boron nitride surface owing to different surface electrostatics. The area-normalized static friction on the graphene diminished as the island area was increased by a power of ~-0.58, suggesting superlubricity behavior. By contrast, the friction on the boron nitride appeared insensitive to the area. Molecular dynamic simulations further showed that the friction coefficient of the water islands on the graphene could reduce to <0.01.

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.

Understanding the Intrinsic Water Wettability of Hexagonal Boron Nitride

The water wettability of hexagonal boron nitride (hBN) has attracted a lot of research interest in the past 15 years. Experimentally, the static water contact angle (WCA) has been widely utilized to characterize the intrinsic water wettability of hBN. In the current study, we have investigated the effect of airborne hydrocarbons and defects on both static and dynamic WCAs of hBN. Our results showed that the static WCA is impacted by defects, which suggests that previously reported static WCAs do not characterize the intrinsic water wettability of hBN since the state-of-the-art hBN samples always have relatively high defect density. Instead, we found that the advancing WCA of freshly exfoliated hBN is not affected by the defects and airborne hydrocarbons. As a result, the advancing WCA on freshly exfoliated hBN, determined to be 79 ± 3°, best represents the intrinsic water wettability of hBN. A qualitative model has been proposed to describe the effect of airborne hydrocarbons and defects on the static and dynamic WCA of hBN, which is well supported by the experimental results. The finding here has important implications for the water wettability of 2D materials.

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.

The effect of water on gold supported chiral graphene nanoribbons: rupture of conjugation by an alternating hydrogenation pattern

In the last few years we have observed a breakpoint in the development of graphene-derived technologies, such as liquid phase filtering and their application to electronics. In most of these cases, they imply exposure of the material to solvents and ambient moisture, either in the fabrication of the material or the final device. The present study demonstrates the sensitivity of graphene nanoribbon (GNR) zigzag edges to water, even in extremely low concentrations. We have addressed the unique reactivity of (3,1)-chiral GNR with moisture on Au(111). Water shows a reductive behaviour, hydrogenating the central carbon of the zigzag segments. By combining scanning tunnelling microscopy (STM) with simulations, we demonstrate how their reactivity reaches a thermodynamic limit when half of the unit cells are reduced, resulting in an alternating pattern of hydrogenated and pristine unit cells starting from the terminal segments. Once a quasi-perfect alternation is reached, the reaction stops regardless of the water concentration. The hydrogenated segments limit the electronic conjugation of the GNR, but the reduction can be reversed both by tip manipulation and annealing. Selective tip-induced dehydrogenation allowed the stabilization of radical states at the edges of the ribbons, while the annealing of the sample completely recovered the original, pristine GNR.

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.

Diameter-Optimum Spreading for the Impinging of Water Nanodroplets on Solid Surfaces

The impinging of water nanodroplets on solid surfaces is crucial to many nanotechnologies. Through large-scale molecular dynamics simulations, the size effect on the spreading of water nanodroplets after impinging on hydrophilic, graphite, and hydrophobic surfaces under low impinging velocities has been systematically studied. The spreading rates of nanodroplets first increase and then decrease and gradually become constant with the increase of nanodroplet diameter. The nanodroplets with the diameters of 17-19 nm possess the highest spreading rates because of the combined effect of the strongest interfacial interaction and the strongest surface interaction within water molecules. The highest water molecule densities, hydrogen bond numbers, and dielectric constants of interface and surface layers mainly contribute to the lowest interface work of adhesion and surface tension values at optimal diameters. These results unveil the nonmonotonic characteristics of spreading velocity, interface work of adhesion and surface tension with nanodroplet diameter for nanodroplets on solid surfaces.

Structural Dynamics in the Presence of Water of Graphene Bilayers with Defects

The dynamics of a bilayer of graphene containing one mono-vacancy in the top layer has been investigated in the framework of DFTB in the absence and in the presence of water. Due to the speed of the code, we can describe details of the behavior, which are not directly accessible experimentally and cannot be treated by DFT or classical molecular dynamics. The presence of water enhances the displacement of carbon atoms in the perpendicular direction to the surface. Our results explain very well a variety of experimental findings. In particular, the stabilization of the Jahn-Teller distortion by hydrogenation of one of the carbon atoms at the edge of a mono-vacancy has been elucidated. This work is the first analysis of the behavior of a graphene vacancy at room temperature in contact with water based on a quantum mechanical molecular dynamics method, where both graphene and solvent are treated at the same level.

Potential and Progress of 2D Materials in Photomedicine for Cancer Treatment

Over the last decades, photomedicine has made a significant impact and progress in treating superficial cancer. With tremendous efforts many of the technologies have entered clinical trials. Photothermal agents (PTAs) have been considered as emerging candidates for accelerating the outcome from photomedicine based cancer treatment. Besides various inorganic and organic candidates, 2D materials such as graphene, boron nitride, and molybdenum disulfide have shown significant potential for photothermal therapy (PTT). The properties such as high surface area to volume, biocompatibility, stability in physiological media, ease of synthesis and functionalization, and high photothermal conversion efficiency have made 2D nanomaterials wonderful candidates for PTT to treat cancer. The targeting or localized activation could be achieved when PTT is combined with chemotherapies, immunotherapies, or photodynamic therapy (PDT) to provide better outcomes with fewer side effects. Though significant development has been made in the field of phototherapeutic drugs, several challenges have restricted the use of PTT in clinical use and hence they have not yet been tested in large clinical trials. In this review, we attempted to discuss the progress, properties, applications, and challenges of 2D materials in the field of PTT and their application in photomedicine.

Interfacial Liquid Water on Graphite, Graphene, and 2D Materials

The optical, electronic, and mechanical properties of graphite, few-layer, and two-dimensional (2D) materials have prompted a considerable number of applications. Biosensing, energy storage, and water desalination illustrate applications that require a molecular-scale understanding of the interfacial water structure on 2D materials. This review introduces the most recent experimental and theoretical advances on the structure of interfacial liquid water on graphite-like and 2D materials surfaces. On pristine conditions, atomic-scale resolution experiments revealed the existence of 1-3 hydration layers. Those layers were separated by ∼0.3 nm. The experimental data were supported by molecular dynamics simulations. However, under standard working conditions, atomic-scale resolution experiments revealed the presence of 2-3 hydrocarbon layers. Those layers were separated by ∼0.5 nm. Linear alkanes were the dominant molecular specie within the hydrocarbon layers. Paradoxically, the interface of an aged 2D material surface immersed in water does not have water molecules on its vicinity. Free-energy considerations favored the replacement of water by alkanes.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

A multiscale approach for electronic transport simulation of carbon nanostructures in aqueous solvent

Theoretical works addressing electronic nano-devices operating in an aqueous environment often neglect solvent effects. In order to assess the role played by the polarization effects on the electronic transport properties of solvated graphene, for example in possible bio-sensing applications, we have used here a combination of polarizable force-field molecular dynamics, hybrid quantum mechanics/molecular mechanics (QM/MM) approach, density functional theory, and non-equilibrium Green's function method. We considered different solvation conditions, the presence of defects in graphene, as well as various choices for the partitions between the quantum and classical regions in QM/MM, in which we explicitly account for polarization effects. Our results show that the polarization effects on graphene lead to changes in the structure of interfacial water molecules which are more pronounced in the vicinity of defects. The presence of water leads to increased scattering due to the long-range charge interactions with graphene. At the same time, changes in the conductance due to polarization or salt concentration are found to be small, paving the way for robust electronic nano-devices operating in aqueous environments.

Ultrahigh resistance of hexagonal boron nitride to mineral scale formation

Formation of mineral scale on a material surface has profound impact on a wide range of natural processes as well as industrial applications. However, how specific material surface characteristics affect the mineral-surface interactions and subsequent mineral scale formation is not well understood. Here we report the superior resistance of hexagonal boron nitride (hBN) to mineral scale formation compared to not only common metal and polymer surfaces but also the highly scaling-resistant graphene, making hBN possibly the most scaling resistant material reported to date. Experimental and simulation results reveal that this ultrahigh scaling-resistance is attributed to the combination of hBN’s atomically-smooth surface, in-plane atomic energy corrugation due to the polar boron-nitrogen bond, and the close match between its interatomic spacing and the size of water molecules. The latter two properties lead to strong polar interactions with water and hence the formation of a dense hydration layer, which strongly hinders the approach of mineral ions and crystals, decreasing both surface heterogeneous nucleation and crystal attachment.

Scale formation may have detrimental effects on the properties and functions of materials’ surfaces. Here the authors report the high scaling resistance of hexagonal boron nitride and relate it to the atomic level structure and interaction with water molecules.

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.

Advancements and Applications in the Composites of Silk Fibroin and Graphene-Based Materials

Silk fibroin and three kinds of graphene-based materials (graphene, graphene oxide, and reduced graphene oxide) have been widely investigated in biomedical fields. Recently, the hybrid composites of silk fibroin and graphene-based materials have attracted much attention owing to their combined advantages, i.e., presenting outstanding biocompatibility, mechanical properties, and excellent electrical conductivity. However, maintaining bio-toxicity and biodegradability at a proper level remains a challenge for other applications. This report describes the first attempt to summarize the hybrid composites’ preparation methods, properties, and applications to the best of our knowledge. We strongly believe that this review will open new doors for coming researchers.

Wettability of graphene oxide functionalized with N-alkylamines: a molecular dynamics study

The wettability of graphene oxide functionalized with N-alkylamines was studied by molecular dynamics simulations. Six different N-alkylamines and two functionalization degrees were reviewed. The nucleophilic ring-opening reaction mechanism between the N-alkylamines and epoxy functional groups of graphene oxide was considered to generate the atomistic models. Water contact angles increased with both the alkyl chain length and substitution degree. The Wenzel model was used to access the effect of both the surface roughness and alkyl chain length on wettability. The results indicated that functionalization introduces an important increase of surface roughness but its effect on wettability is countered by the alkyl chain length.

Role of Interfacial Bonding in Tribochemical Wear

Tribochemical wear of contact materials is an important issue in science and engineering. Understanding the mechanisms of tribochemical wear at an atomic scale is favorable to avoid device failure, improve the durability of materials, and even achieve ultra-precision manufacturing. Hence, this article reviews some of the latest developments of tribochemical wear of typical materials at micro/nano-scale that are commonly used as solid lubricants, tribo-elements, or structural materials of the micro-electromechanical devices, focusing on their universal mechanisms based on the studies from experiments and numerical simulations. Particular focus is given to the fact that the friction-induced formation of interfacial bonding plays a critical role in the wear of frictional systems at the atomic scale.

On the effects of induced polarizability at the water-graphene interface via classical charge-on-spring models

Molecular models of the water-graphene interaction are essential to describe graphene in condensed phases. Different challenges are associated with the generation of these models, in particular π-π and dispersion interactions; thus quantum and classical models have been developed and due to the numerical efficiency of the latter, they have been extensively employed. In this work, we have systematically studied, via molecular dynamics, two polarizable graphene models, denominated CCCP and CCCPD, employing the charge-on-spring model of the GROMOS forcefield, both being compatible with the polarizable water models COS/G2 and COS/D2, respectively. These models were compared with non-polarizable graphene and SPC water models. We focused the study on the water-graphene interface in two distinct systems and under the influence of an electric field: one composed of graphene immersed in water and the other composed of graphene with a water droplet above it. In the former, the orientation of water close to the graphene layer is affected by polarizable graphene in comparison to non-polarizable graphene. This effect is emphasised when an electric field is applied. In the latter, carbon polarizability reduced water contact angles, but graphene retained its hydrophobicity and the computed angles are within the experimental data. Given the significant extra computational cost, the use of polarizable models instead of the traditional fixed-charged approach for the graphene-water interaction may be justified when polarizability effects are relevant, for example, when applying relatively strong fields or in very anisotropic systems, such as the vacuum-bulk interface, as these models are more responsive to such conditions.

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.

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

Atomic-scale defects are ubiquitous in nanomaterials, yet their role in modulating fluid flow is inadequately understood. Hexagonal boron nitride (hBN) is an important two-dimensional material with applications in desalination and osmotic power. Although pristine hBN offers higher friction to the flow of water than graphene, we show here that certain defects can enhance water slippage on hBN. Using classical molecular dynamics simulations assisted by quantum-mechanical density functional theory, we compute the friction coefficient of water on hBN containing various vacancies (B, N, BN, B₂N, and B₃N) and the Stone-Wales defect. By investigating two defect concentrations, we obtain friction coefficients ranging from 0.4 to 2.6 times that of pristine hBN, leading to a maximum water slip length of 18.1 nm on hBN with a N vacancy or a Stone-Wales defect. Our work informs the use of defects to tune water flow and reveals defective hBN as an alternative high-slip surface to graphene.

Competition between Hydrogen Bonding and Dispersion Force in Water Adsorption and Epoxy Adhesion to Boron Nitride: From the Flat to the Curved

Hexagonal boron nitride (h-BN) is a material with excellent thermal conductivity and electrical insulation, used as an additive to various matrices. To increase the affinity of h-BN to them, hydrogen bonds should be formed at the interface. In reality, however, they are not formed; the N atoms are not capable of accepting hydrogen bonds due to the delocalization of their lone pair electrons over the B-N π bonds. To make it form hydrogen bonds, one may need to break the planarity of h-BN so that the orbital overlap in the B-N π bonds can be reduced. This idea is verified with first-principles calculations on the adsorption of a water molecule on hypothetical h-BN surfaces, the planarity of which is broken. One can do it in silico but not in vitro. BN nanotubes (BNNTs) are considered as a more realistic BN surface with nonplanarity. The hydrogen bond is shown to become stronger as the curvature of the tube increases. On the contrary, the strength of the dispersion force acting at the interface becomes weaker. In water adsorption, these two interactions are in competition with each other. However, in epoxy adhesion, the interaction due to dispersion forces is overwhelmingly stronger than that due to hydrogen bonding. The smaller the curvature of the surface, the smaller the distance between more atoms at the interface; thus, the interaction due to dispersion forces maximized.

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.

Effect of Graphene on Ice Polymorph

Recently, ice with stacking disorder structure, consisting of random sequences of cubic ice (Ic) and hexagonal ice (Ih) layers, was reported to be more stable than pure Ih/Ic. Due to a much lower free energy barrier of heterogeneous nucleation, in practice, the freezing process of water is controlled by heterogeneous nucleation triggered by an external medium. Therefore, we carry out molecular dynamic simulations to explore how ice polymorphism depends on the lattice structure of the crystalline substrates on which the ice is grown, focusing on the primary source of atmospheric aerosols, carbon materials. It turns out that, during the nucleation stage, the polymorph of ice nuclei is strongly affected by graphene substrates. For ice nucleation on graphene, we find Ih is the dominant polymorph. This can be attributed to structural similarities between graphene and basal face of Ih. Our results also suggest that the substrate only affects the polymorph of ice close to the graphene surface, with the preference for Ih diminishing as the ice layer grows.

Aqueous Affinity and Interfacial Dynamics of Anisotropic Buckled Black Phosphorous

The structure of black phosphorous (BP) is similar to the honeycomb arrangement of graphene, but the layered BP is found to be buckled and highly anisotropic. The buckled surface structure affects interfacial molecule mobility and plays a vital role in various nanomaterial applications. The BP is also known for wettability, droplet formation, stability, and hydrophobicity in the aqueous environment. However, there is a gap concerning the structural and dynamical behavior of water molecules, which is available in abundance for other monoatomic and polyatomic two-dimensional (2D) materials. Motivated by the technological importance, we try to bridge the gap by explaining the surface anisotropy-facilitated behavior of water molecules on bilayer BP using classical and first principles molecular dynamics (MD) simulations. From our classical MD study, we find three distinct layers of water molecules. The water layer closest to the interface is L1, followed by L2 and L3/bulk perpendicular to the BP surface. Water molecules in the L1 layer experience some structural disintegration in hydrogen bond (HB) phenomena compared to the bulk. There is a loss of HB donor-acceptor count per water molecule. The average HB count decreases because of an elevated rate of HB formation and deformation; this would affect the dynamic properties in terms of HB lifetime. Therefore, we observe the reduced lifetime of HB in the layer in close contact with BP, which again complements our finding on the diffusion coefficient of water molecules in distinct layers. Water diffuses relatively faster with diffusion coefficient 3.25 × 10^-9 m² s^-1 in L1, followed by L2 and L3. The BP layer shows moderate hydrophobic nature. Our results also indicate the anisotropic behavior as the diffusion along the x-direction is faster than that along the y-direction. The gap in the slope of the x and y components of mean-squared displacement (MSD) complements the pinning effect in an aqueous environment. We observe blue-shifted and red-shifted libration and O-H stretching modes from the calculated power spectra for the L1 water molecules compared to the L2 and L3 molecules from first principles MD simulations. Our analysis may help understand the physical phenomena that occur during the surface wetting of the predroplet formation process observed experimentally.

Surface thermodynamics at the nanoscale

Fluid interfaces with nanoscale radii of curvature are generating great interest, both for their applications and as tools to probe our fundamental understanding. One important question is what is the smallest radius of curvature at which the three main thermodynamic combined equilibrium equations are valid: the Kelvin equation for the effect of curvature on vapor pressure, the Gibbs-Thomson equation for the curvature-induced freezing point depression, and the Ostwald-Freundlich equation for the curvature-induced increase in solubility. The objective of this Perspective is to provide conceptual, molecular modeling, and experimental support for the validity of these thermodynamic combined equilibrium equations down to the smallest interfacial radii of curvature. Important concepts underpinning thermodynamics, including ensemble averaging and Gibbs's treatment of bulk phase heterogeneities in the region of an interface, give reason to believe that these equations might be valid to smaller scales than was previously thought. There is significant molecular modeling and experimental support for all three of the Kelvin equation, the Gibbs-Thomson equation, and the Ostwald-Freundlich equation for interfacial radii of curvature from 1 to 4 nm. There is even evidence of sub-nanometer quantitative accuracy for the Kelvin equation and the Gibbs-Thomson equation.

ReaxFF reactive molecular dynamics simulations to study the interfacial dynamics between defective h-BN nanosheets and water nanodroplets

In this work, the authors have developed a reactive force field (ReaxFF) to investigate the effect of water molecules on the interfacial interactions with vacancy defective hexagonal boron nitride (h-BN) nanosheets by introducing parameters suitable for the B/N/O/H chemistry. Initially, molecular dynamics simulations were performed to validate the structural stability and hydrophobic nature of h-BN nanosheets. The water molecule dissociation mechanism in the vicinity of vacancy defective h-BN nanosheets was investigated, and it was shown that the terminal nitrogen and boron atoms bond with a hydrogen atom and hydroxyl group, respectively. Furthermore, it is predicted that the water molecules arrange themselves in layers when compressed in between two h-BN nanosheets, and the h-BN nanosheet fracture nucleates from the vacancy defect site. Simulations at elevated temperatures were carried out to explore the water molecule trajectory near the functionalized h-BN pores, and it was observed that the intermolecular hydrogen bonds lead to agglomeration of water molecules near these pores when the temperature was lowered to room temperature. The study was extended to observe the effect of pore sizes and temperatures on the contact angle made by a water nanodroplet on h-BN nanosheets, and it was concluded that the contact angle would be less at higher temperatures and larger pore sizes. This study provides important information for the use of h-BN nanosheets in nanodevices for water desalination and underwater applications, as these h-BN nanosheets possess the desired adsorption capability and structural stability.

Associated molecular liquids at the graphene monolayer interface

We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.

Graphene-facilitated synthesized vertically aligned hexagonal boron nitride nanowalls and their gas adsorption properties

The capability of hexagonal boron nitride (h-BN) to adsorb gas atoms may stimulate various promising applications in environment remediation and energy storage, while the interactivity with gas molecules yet remains challenging due to its inherent chemical inertness. In this article, we report a feasible and effective route for the scalable synthesis of vertically aligned h-BN nanowalls assisted by reduced graphene oxide (rGO) without metallic catalysts. The average thickness of the fine h-BN nanowalls is few-atomic layers about 3.7 nm, that grow on the large substrate-like flakes transformed from the pristine rGO. The hierarchical h-BN nanowalls exhibit an enhanced gas adsorption performance, not only through physisorption owing to the synergistic combination of different porous geometries, but also through chemisorption via the open edge groups. Moreover, it demonstrates a significantly enhanced adsorption of CO₂ over CH₄ as compared to the h-BN nanosheets with similar sizes. Density functional theory calculations reveal that the -OH edge groups can effectively increase the adsorption capability towards CO₂, accompanied by a shortened adsorption distance when the gas molecule is energetically stabilized. The wetting characteristics of h-BN nanowalls was further examined by contact angle goniometry.

Screening effect of monolayer van der Waals crystals on surface deicing: a molecular simulation study

Our extensive molecular dynamics simulations reveal a significant screening effect of monolayer graphene and hexagonal boron nitride (h-BN) on surface deicing of substrates with different degrees of hydrophilicity, including superhydrophilic (SHP) and superhydrophobic (SHB) substrates. Compared with bare surfaces, graphene and h-BN reduce the interfacial shear strength and the normal detaching strength of ice on an SHP substrate but increase the shear and detaching strengths on hydrophobic and SHB substrates. However, the shear and detaching strengths of ice become approximately unified on all of the surfaces, when interface ice layers melt into liquid water, demonstrating the screening capability from graphene and h-BN that weakens the influence of substrates on ice adhesion. Graphene and h-BN coatings suppress ice premelting on the SHP surface and change the dielectric constant of interface ice or water. This work could deepen our understanding of the role of van der Waals crystals in deicing coating.

Atomic Steps Induce the Aligned Growth of Ice Crystals on Graphite Surfaces

Heterogeneous ice nucleation on atmospheric aerosols strongly affects the earth's climate, and at the microscopic level, surface-irregularity-induced ice crystallization behaviors are common but crucial. Because of the lack of visual evidence and effective experimental methods, the mechanism of atomic-structure-dependent ice formation on aerosol surfaces is poorly understood. Here we chose highly oriented pyrolytic graphite (HOPG) to represent soot (a primary aerosol), and environmental scanning electron microscopy (ESEM) was performed for in situ observations of ice formation. We found that hexagonal ice crystals show an aligned growth pattern via a two-stage pathway with one a axis coinciding with the direction of atomic step edges on the HOPG surface. Additionally, the ice crystals grow at a noticeably higher speed along this direction. This study reveals the role of atomic surface defects in heterogeneous ice nucleation and may pave the way to control icing-related processes in practical applications.

Enhanced nanoparticle rejection in aligned boron nitride nanotube membranes

The rejection of particles with different charges and sizes, ranging from a few Ångstroms to tens of nanometers, is key to a wide range of industrial applications, from wastewater treatment to product purification in biotech processes. Carbon nanotubes (CNTs) have long held the promise to revolutionize filtration, with orders of magnitude higher fluxes compared to commercial membranes. CNTs, however, can only reject particles and ions wider than their internal diameter. In this work, the fabrication of aligned boron nitride nanotube (BNNT) membranes capable of rejecting nanoparticles smaller than their internal diameter is reported for the first time. This is due to a mechanism of charge-based rejection in addition to the size-based one, enabled by the BNNTs surface structure and chemistry and elucidated here using high fidelity molecular dynamics and Brownian dynamics simulations. This results in ∼40% higher rejection of the same particles by BNNT membranes than CNT ones with comparable nanotube diameter. Furthermore, since permeance is proportional to the square of the nanotubes' diameter, using BNNT membranes with ∼30% larger nanotube diameter than a CNT membrane with comparable rejection would result in up to 70% higher permeance. These results open the way to the design of more effective nanotube membranes, capable of high particle rejection and, at the same time, high water permeance.

Mechanical strength of a nanoporous bicrystalline h-BN nanomembrane in a water submerged state

Due to superior water permeability, structural stability, and adsorption capability, h-BN nanosheets are emerging as an efficient membrane for water desalination. In order to cater to the demand for potable water, large size membranes are required to maintain a high desalination rate from water purification systems. These large size membranes usually contain polycrystals with an offset in their mechanical properties from pristine h-BN nanosheets. In this article, molecular dynamics based simulations were performed in conjunction with a hybrid interatomic potential (reactive force field, TIP3P, and Lennard Jones) to simulate the mechanical strength of nanoporous single and bicrystalline h-BN nanosheets under water submerged conditions. The interaction between the atomic configuration of grain boundary atoms and nanopores in the presence of water molecules helps in investigating the viability of defective h-BN nanomembranes for underwater applications. Higher dislocation density enhances the mechanical strength of nanoporous bicrystalline h-BN nanosheets containing twin nanopores, which makes them a better substitute for water submerged applications as compared to the pristine nanosheets. The mechanical strength of nanoporous single crystalline h-BN nanosheets deteriorates with an increase in the number of nanopores, whereas a contrasting trend was observed with bicrystalline h-BN nanosheets.

Solvent-Solvent Correlations across Graphene: The Effect of Image Charges

Wetting experiments show pure graphene to be weakly hydrophilic, but its contact angle (CA) also reflects the character of the supporting material. Measurements and molecular dynamics simulations on suspended and supported graphene often reveal a CA reduction due to the presence of the supporting substrate. A similar reduction is consistently observed when graphene is wetted from both sides. The effect has been attributed to transparency to molecular interactions across the graphene sheet; however, the possibility of substrate-induced graphene polarization has also been considered. Computer simulations of CA on graphene have so far been determined by ignoring the material's conducting properties. We improve the graphene model by incorporating its conductivity according to the constant applied potential molecular dynamics. Using this method, we compare the wettabilities of suspended graphene and graphene supported by water by measuring the CA of cylindrical water drops on the sheets. The inclusion of graphene conductivity and concomitant polarization effects leads to a lower CA on suspended graphene, but the CA reduction is significantly bigger when the sheets are also wetted from the opposite side. The stronger adhesion is accompanied by a profound change in the correlations among water molecules across the sheet. While partial charges on water molecules interacting across an insulator sheet attract charges of the opposite sign, apparent attraction among like charges is manifested across the conducting graphene. The change is associated with graphene polarization, as the image charges inside the conductor attract equally signed partial charges of water molecules on both sides of the sheet. Additionally, using a nonpolar liquid (diiodomethane), we affirm a detectable wetting translucency when liquid-liquid forces are dominated by dispersive interactions. Our findings are important for predictive modeling toward a variety of applications including sensors, fuel cell membranes, water filtration, and graphene-based electrode materials in high-performance supercapacitors.

Exploring the limits of sensitivity for strain gauges of graphene and hexagonal boron nitride decorated with metallic nanoislands

The purpose of this work is to clarify the mechanism of piezoresistance in a class of ultra-sensitive strain gauges based on metallic films on 2D substrates ("2D/M" films). The metals used are gold or palladium deposited as ultrathin films (≤16 nm). These films transition from a regime of subcontiguous growth to a percolated morphology with increasing nominal thickness. The 2D substrates are either single-layer graphene or hexagonal boron nitride (hBN). By using either a conductor (graphene) or an insulator (hBN), it is possible to de-couple the relative contributions of the metal and the 2D substrate from the overall piezoresistance of the composite structure. Here, we use a combination of measurements including electron microscopy, automated image analysis, temperature-dependent conductivity, and measurements of gauge factor of the films as they are bent over a 1 μm step edge (0.0001% or 1 ppm). Our observations are enumerated as follows: (1) of the four permutations of metal and 2D substrate, all combinations except hBN/Au are able to resolve 1 ppm strain (considered extraordinary for strain gauges) at some threshold thickness of metal; (2) for non-contiguous (i.e., unpercolated) films of metal on hBN, changes in resistance for these small step strains cannot be detected; (3) for percolated films on hBN, changes in resistance upon strain can be resolved only for palladium and not for gold; (4) graphene does not exhibit detectable changes in resistance when subjected to step strains of either 1 or 10 ppm, but does so upon the deposition of any amount of gold or palladium, even for nominal thicknesses below the threshold for percolation. Our observations reveal unexpected complexity in the properties of these simple composite materials, and ways in which these materials might be combined to exhibit even greater sensitivity.

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.

Closest-Packing Water Monolayer Stably Intercalated in Phyllosilicate Minerals under High Pressure

The directional hydrogen-bond (HB) network and nondirectional van der Waals (vdW) interactions make up the specificity of water. Directional HBs could construct an ice-like monolayer in hydrophobic confinement even in the ambient regime. Here, we report a water monolayer dominated by vdW interactions confined in a phyllosilicate interlayer under high pressure. Surprisingly, it was in a thermodynamically stable state coupled with bulk water at the same pressure (P) and temperature (T), as revealed by the thermodynamic integration approach on the basis of molecular dynamics (MD) simulations. Both classical and ab initio MD simulations showed water O atoms were stably trapped and exhibited an ordered hexagonal closest-packing arrangement, but OH bonds of water reoriented frequently and exhibited a specific two-stage reorientation relaxation. Strikingly, hydration in the interlayer under high pressure had no relevance with surface hydrophilicity rationalized by the HB forming ability, which, however, determines wetting in the ambient regime. Intercalated water molecules were trapped by vdW interactions, which shaped the closest-packing arrangement and made hydration energetically available. The high pressure-volume term largely drives hydration, as it compensates the entropy penalty which is restricted by a relatively lower temperature. This vdW water monolayer should be ubiquitous in the high pressure but low-temperature regime.

On the wetting translucency of hexagonal boron nitride

When a drop sits on an atomically thin coating supported by a hydrophilic material, it is possible that the underlying substrate influences the equilibrium contact angle. Such behavior is known as the wetting translucency effect.

Insights into water permeability and Hg2+ removal using two-dimensional nanoporous boron nitride

Industrial wastewater containing Hg²⁺, when discharged into nature, will pose a serious threat to ecological security.

A Bespoke Force Field To Describe Biomolecule Adsorption at the Aqueous Boron Nitride Interface

Reliable manipulation of the interface between 2D nanomaterials and biomolecules represents a current frontier in nanoscience. The ability to resolve the molecular-level structures of these biointerfaces would provide a fundamental data set that is needed to enable systematic and knowledge-based progress in this area. These structures are challenging to obtain via experiment alone, and molecular simulations offer a complementary approach to address this problem. Compared with graphene, the interface between hexagonal boron nitride (h-BN) and biomolecules is relatively understudied at present. While several force fields are currently available for modeling the h-BN/water interface, there is a lack of a suitable force field that can describe the interactions between h-BN, liquid water, and biomolecules. Here, we use density functional theory calculations to create a force field, BoNi-CHARMM, to describe biomolecular interactions at the aqueous h-BN interface. Verifying our force field presents an additional challenge, given the scarcity of available experimental data for these interfaces. We test our force field against experimental evidence regarding the water/surface contact angle and confirm that the force field provides experimentally consistent values. We also present preliminary data regarding predictions of the free energy of adsorption of a selection of amino acids at the aqueous h-BN interface, revealing arginine and tryptophan to be among the strongest binders. This force field provides an opportunity to initiate a systematic progression in our current understanding of how to capture the intermolecular interactions at the h-BN biointerface.

Highly Thermally Conductive Polymer Composite Originated from Assembly of Boron Nitride at an Oil-Water Interface

Thermally conductive polymer packaging material is of great significance for the thermal management of electronics. Inorganic thermally conductive fillers have been demonstrated as a convenient approach to achieve this goal by sacrificing the lightweight and processability of the polymer. To address this problem, an effective 3D boron nitride (BN) network was constructed as a heat conduction pathway in a polystyrene (PS) matrix based on an oil-water interface assembly in this work. Styrene oil droplets were stabilized by BN sheets in the water phase to form Pickering emulsions, and then in situ polymerization was trigged to synthesize PS microspheres with ultrathin BN layer-covered surfaces (PS@BN microspheres). Composite substrates were fabricated through hot-compressing the PS@BN microspheres to form BN networks based on the original microsphere template. Benefited from the network structure, the maximum thermal conductivity of the composite substrate reached 0.94 W/mK at 33.3 wt % BN, which is 626% folds of that of pure PS. It was also demonstrated that the storage modulus and thermal stability of the composite substrate were dramatically improved by the BN network. The reported composite substrate and its fabrication strategy are promising in the development of thermal management of electronics.

Passivation of Layered Gallium Telluride by Double Encapsulation with Graphene

Layered semiconductor gallium telluride (GaTe) undergoes a rapid structural transition to a degraded phase in ambient conditions, limiting its utility in devices such as optical switches. In this work, we demonstrate that the degradation process in GaTe flakes can be slowed down dramatically via encapsulation with graphene. Through examining Raman signatures of degradation, we show that the choice of substrate significantly impacts the degradation rate and that the process is accelerated by the transfer of GaTe to hydrophilic substrates such as SiO₂/Si. We find that double encapsulation with both top and bottom graphene layers can extend the lifetime of the material for several weeks. The photoresponse of flakes encapsulated in this way is only reduced by 17.6 ± 0.4% after 2 weeks, whereas unencapsulated flakes display no response after this time. Our results demonstrate the potential for alternative, van der Waals material-based passivation strategies in unstable layered materials and highlight the need for careful selection of substrates for 2D electronic devices.

Liquefaction of water on the surface of anisotropic two-dimensional atomic layered black phosphorus

The growth and wetting of water on two-dimensional(2D) materials are important to understand the development of 2D material based electronic, optoelectronic, and nanomechanical devices. Here, we visualize the liquefaction processes of water on the surface of graphene, MoS2 and black phosphorus (BP) via optical microscopy. We show that the shape of the water droplets forming on the surface of BP, which is anisotropic, is elliptical. In contrast, droplets are rounded when they form on the surface of graphene or MoS2, which do not possess orthometric anisotropy. Molecular simulations show that the anisotropic liquefaction process of water on the surface of BP is attributed to the different binding energies of H2O molecules on BP along the armchair and zigzag directions. The results not only reveal the anisotropic nature of water liquefaction on the BP surface but also provide a way for fast and nondestructive determination of the crystalline orientation of BP.