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

Contributions of Edge and Internal Atoms to the Friction of Two-Dimensional Heterojunctions

The weak interlayer interaction and strong intralayer bonding in Van der Waals (vdW) layered materials give rise to ultralow friction at incommensurate contact interfaces, a phenomenon known as structural superlubricity. This phenomenon is complicated by the interplay between atomic degrees of freedom, twist angle, and normal force. In this Letter, we exploit naturally occurring cracks in vdW crystals and microfabrication techniques to qualitatively separate the contributions of edge, complete moiré tile, and incomplete moiré rim regions. It is observed that friction from the incomplete moiré rim region, scaling linearly with the rim area, plays a dominant role in determining the twist angle dependence of friction. Interestingly, despite lower friction contributions from complete moirés, friction from a complete moiré tile is independent of moiré size; consequently, friction from the moiré tile scales linearly with the total number of moirés. By integrating the contributions from edge, moiré tile, and incomplete moiré rim, a friction law for vdW heterojunctions is proposed. Finally, friction changes from positively to negatively correlating with normal force, contingent on the suppression of moiré ridge and the dissipation from edge atoms.

Exploring hydrophobicity or hydrophilicity of borophene surface via reactive molecular dynamics simulation

Borophene, a novel two-dimensional material unveiled in 1998, has garnered significant interest among researchers due to its distinct mechanical and electrical characteristics. Efforts to experimentally synthesize borophene continue to captivate researchers' interest in recent years. Given the current lack of experimental studies on the interaction between water and the borophene surface, molecular dynamics simulation offers a valuable approach for predicting the substance's reactivity with water. Additionally, such simulations can assess the hydrophilicity and hydrophobicity of borophene, providing valuable insights into its properties. In our current research, we utilized reactive molecular dynamics simulation to investigate the wetting behavior of borophene. Our findings reveal that the borophene surface exhibits hydrophobic characteristics, demonstrating anisotropic wettability. Specifically, the water contact angle was calculated to be 149.11° along the zigzag direction and 148.4° along the armchair direction. The contour map of the interaction energy between a water molecule and the borophene surface revealed a notable energy barrier in the zigzag direction. This barrier contributes to the asymmetric spreading of the water droplet on the surface. Density profiles and radial pair distribution function (RDF) diagrams of the water droplet on the borophene surface further corroborated the hydrophobic nature of borophene by indicating a significant distance between the water droplet and the surface. Moreover, analysis of the number of hydrogen bonds demonstrated that borophene efficiently utilizes nearly all its capacity to form hydrogen bonds. Additionally, we compared the wettability of borophene with that of other two-dimensional materials, such as various graphene allotropes and phosphorene, which have been subjects of recent investigation.

Ab Initio Molecular Dynamics Investigation on the Permeation of Sodium and Chloride Ions Through Nanopores in Graphene and Hexagonal Boron Nitride Membranes

Nanoporous membranes promise energy-efficient water desalination. Hexagonal boron nitride (h-BN), like graphene, exhibits outstanding physical and chemical properties, making it a promising candidate for water treatment. We employed Car-Parrinello molecular dynamics simulations to establish an accurate modeling of Na+ and Cl- permeation through hydrogen passivated nanopores in graphene and h-BN membranes. We demonstrate that ion separation works well for the h-BN system by imposing a barrier of 0.13 eV and 0.24 eV for Na+ and Cl- permeation, respectively. In contrast, for permeation of the graphene nanopore, the Cl- ion faces a minimum of energy of 0.68 eV in the nanopore plane and is prone toward blockade of the nanopore, while the Na+ ion experiences a slight minimum of 0.03 eV. Overall, the desalination performance of h-BN nanopores surpasses that of their graphene counterparts.

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.

Strength, number, and kinetics of hydrogen bonds for water confined inside boron nitride nanotubes

Water has shown a myriad of highly interesting properties and behaviors, such as very low friction, phase transition under unexpected conditions, massive property alterations, etc. inside strong nanoconfinements of few-nanometer to sub-nanometer diameters. Water-water hydrogen bonding is one of the most important factors dictating such water behavior and properties inside such strong nanoconfinements. In this paper, we employ Reactive Force Field (ReaxFF) molecular dynamics (MD) simulations for studying multiple facets of such water-water hydrogen bonds (HBs) inside boron-nitride nanotubes (BNNTs) having diameters ranging from a few nanometers to sub-nanometers. First, the strength of the water-water HB interactions, as a function of the HB configuration, is quantified by studying the corresponding PMF (potential of mean force). For water present in extreme confinements (BNNTs with sub-nanometric diameters), we see completely isolated HB basins. On the other hand, for bulk water the HB basin is connected (via a saddle point) to a nearby second PMF well. Therefore, our analysis successfully distinguishes the HB characteristics between the cases of water in extreme confinement and bulk water. Second, we study the kinetics of such water-water HBs: HBs formed by a given pair of water molecules in extreme confinements show a much larger probability of remaining intact once formed or re-forming after they have been broken. Both these results, which shed new light on water-water hydrogen bonding inside strong nanoconfinements, can be explained by the single-file structure formed by the water molecules in extreme BNNT confinements.

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.

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.

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.

Covalently bonded interface in polymer/boron nitride nanosheet composite toward enhanced mechanical and thermal behaviour

This experimental study aimed to enhance the mechanical and thermal properties of BN (hexagonal boron nitride) nanosheet-reinforced high-density polyethylene by functionalizing its interface. The challenges associated with this nanocomposites are its poor dispersion and weak interface. Accordingly, to improve the load transfer at the interface, BN nanosheets were chemically modified with silane functional groups ((3-aminopropyl)tri-ethoxy silane), making it possible to form covalent bonds between the maleic anhydride-grafted polyethylene and nanosheet. Consequently, three different types of nanocomposite samples were fabricated based on the covalently bonded or non-bonded interface. Two nanocomposite configurations featured a non-bonded interface between the nanofiller and PE matrix (p-BN/PE and (silane functionalized) s-BN/PE). In contrast, the third configuration had a covalently bonded interface (silane-functionalized h-BN + maleic anhydride-grafted PE, i.e., PE-g-BN). According to the zeta potential analysis, the silane-functionalized BN nanosheets were stable suspensions and uniformly dispersed in the polymer matrix. The tensile and flexure strength of the nanocomposites showed over 100% improvement due to the covalently bonded interface. The lamellae structure of PE in the bonded interface samples was responsible for achieving higher mechanical strength in the nanocomposites. Furthermore, the thermal conductivity of the nanocomposites was significantly affected by the type of interfacial bonding, BN wt%, and operating temperature.

Introducing Structures from Hexagonal Borophene to Nitrophene and Their Thermal Conductivity Investigation Using a Reactive Molecular Dynamics Simulation

In the present work, the thermal conductivity (TC) of hexagonal structures of boron nitride and borophene was investigated by a reactive molecular dynamics (MD) simulation. Also, to figure out the effect of the boron and nitrogen in the hexagonal structure, five other hypothetical structures were created (in addition to the structure of boron nitride and borophene) and their structures were represented by the symbol BxNy, where x refers to the number of boron atoms and y refers to the number of nitrogen atoms. In this regard, B6N0 refers to borophene, B3N3 is boron nitride, and B0N6 is called nitrophene. The TC of B6N0 and B3N3 structures was calculated and compared with the literature values. Besides these two compounds, the five other structures have not been experimentally synthesized yet, so the TC of the five other hypothetical structures were predicted in the present work. The lowest TC belonged to B3N3, and the highest one was for B0N6. Based on the inherent potential of reactive MD simulation, during TC calculation, atoms' coordination and partial charges are changed and new bonds, rings, or even defects were automatically created on the surfaces. The coordination contour map showed that in B3N3, the atoms have collective movements like a large and single wave, while B0N6 and B6N0 have small group movements as vibrations. So, it became clear that the higher stability of structures caused more curved movements. In addition, the contour map of partial charges is calculated, and the results showed that the high differences in partial charge between atoms in the structure cause high TC, while small charge differences in the structure inhibit heat transfer and cause lower TC.

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.

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.