Dominance of Dispersion Interactions and Entropy over Electrostatics in Determining the Wettability and Friction of Two-Dimensional MoS2 Surfaces

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.

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.

High-Performance Broadband Image Sensing Photodetector Based on MnTe/WS2 van der Waals Epitaxial Heterostructures

Two-dimensional transition metal dichalcogenide (TMDC) heterostructure is receiving considerable attention due to its novel electronic, optoelectronic, and spintronic devices with design-oriented and functional features. However, direct design and synthesis of high-quality TMDC/MnTe heterostructures remain difficult, which severely impede further investigations of semiconductor/magnetic semiconductor devices. Herein, the synthesis of high-quality vertically stacked WS2/MnTe heterostructures is realized via a two-step chemical vapor deposition method. Raman, photoluminescence, and scanning transmission electron microscopy characterizations reveal the high-quality and atomically sharp interfaces of the WS2/MnTe heterostructure. WS2/MnTe-based van der Waals field effect transistors demonstrate high rectification behavior with rectification ratio up to 106, as well as a typical p-n electrical transport characteristic. Notably, the fabricated WS2/MnTe photodetector exhibits sensitive and broadband photoresponse ranging from UV to NIR with a maximum responsivity of 1.2 × 103 A/W, a high external quantum efficiency of 2.7 × 105%, and fast photoresponse time of ∼50 ms. Moreover, WS2/MnTe heterostructure photodetectors possess a broadband image sensing capability at room temperature, suggesting potential applications in next-generation high-performance and broadband image sensing photodetectors.

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.

Biomolecular Adsorprion at ZnS Nanomaterials: A Molecular Dynamics Simulation Study of the Adsorption Preferences, Effects of the Surface Curvature and Coating

The understanding of interactions between nanomaterials and biological molecules is of primary importance for biomedical applications of nanomaterials, as well as for the evaluation of their possible toxic effects. Here, we carried out extensive molecular dynamics simulations of the adsorption properties of about 30 small molecules representing biomolecular fragments at ZnS surfaces in aqueous media. We computed adsorption free energies and potentials of mean force of amino acid side chain analogs, lipids, and sugar fragments to ZnS (110) crystal surface and to a spherical ZnS nanoparticle. Furthermore, we investigated the effect of poly-methylmethacrylate (PMMA) coating on the adsorption preferences of biomolecules to ZnS. We found that only a few anionic molecules: aspartic and glutamic acids side chains, as well as the anionic form of cysteine show significant binding to pristine ZnS surface, while other molecules show weak or no binding. Spherical ZnS nanoparticles show stronger binding of these molecules due to binding at the edges between different surface facets. Coating of ZnS by PMMA changes binding preferences drastically: the molecules that adsorb to a pristine ZnS surface do not adsorb on PMMA-coated surfaces, while some others, particularly hydrophobic or aromatic amino-acids, show high binding affinity due to binding to the coating. We investigate further the hydration properties of the ZnS surface and relate them to the binding preferences of biomolecules.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

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

Recent Advancements in the Recovery and Reuse of Organic Solvents Using Novel Nanomaterial-Based Membranes for Renewable Energy Applications

The energy crisis in the world is increasing rapidly owing to the shortage of fossil fuel reserves. Climate change and an increase in global warming necessitates a change in focus from petroleum-based fuels to renewable fuels such as biofuels. The remodeling of existing separation processes using various nanomaterials is of a growing interest to industrial separation methods. Recently, the design of membrane technologies has been the most focused research area concerning fermentation broth to enhance performance efficiency, while recovering those byproducts to be used as value added fuels. Specifically, the use of novel nano material membranes, which brings about a selective permeation of the byproducts, such as organic solvent, from the fermentation broth, positively affects the fermentation kinetics by eliminating the issue of product inhibition. In this review, which and how membrane-based technologies using novel materials can improve the separation performance of organic solvents is considered. In particular, technical approaches suggested in previous studies are discussed with the goal of emphasizing benefits and problems faced in order to direct research towards an optimized membrane separation performance for renewable fuel production on a commercial scale.

Wettability and Its Influencing Factors of Tight Sandstones in Coal Measures in Ordos Basin, China

With the recent development in exploration technology, extensive attention has been given to the tight sandstone gas reservoirs in coal measures. As a surface property of tight sandstone, wettability is a very important factor that controls the distribution and transport of gas and water inside the sandstone; thus, it plays a vital role in the sandstone gas-recovery process. In this study, a series of tests, namely, the Amott-Harvey wettability test, contact angle test, ζ-potential test, X-ray diffraction, thin-section analysis, scanning electron microscopy, and mercury intrusion porosimetry, were conducted to investigate the basic features of tight sandstones in coal measures and the effects of petrological characteristics, fluid properties, the underground environment, sedimentation, and diagenesis on the wetting behavior of sandstone. The impacts of three major forces, namely, the structural, electrostatic, and dispersion forces, on sandstone wettability were also discussed. The results showed that high surface roughness, high textural and compositional maturity, low organic matter content, a strong hydrodynamic sedimentary environment, weak compaction and cementation, high temperature, and high solution pH, as well as a low ionic concentration, enhance the hydrophilicity of the sandstone.

Polyaromatic cores for the exfoliation of popular 2D materials

Two-dimensional (2D) nanomaterials have attracted interest from the scientific community due to their unique properties. The production of these materials has been carried out by diverse methodologies, the liquid phase exfoliation being the most promising one due to its simplicity and potential scalability. The use of several stabilizers allows to obtain dispersions of these 2D nanomaterials in solvents with low boiling points. Herein we describe a general exfoliation method for different 2D materials employing a biphasic water/dichloromethane system and two different (poly)aromatic hydrocarbons (PAHs). This method allows us to obtain dispersions of the exfoliated 2D materials with high concentrations in the organic solvent. Due to the low boiling point of dichloromethane, and therefore its easy removal, the obtained dispersions can be employed as additives for different composites. We corroborate that the exfoliation efficiency is improved due to the π-π and van der Waals interactions between the PAHs and the layers of the 2D materials.

Patchwork-Structured Heterointerface of 1T-WS2/a-WO3 with Sustained Hydrogen Spillover as a Highly Efficient Hydrogen Evolution Reaction Electrocatalyst

Using tungsten disulfide (WS2) as a hydrogen evolution reaction (HER) electrocatalyst brought on several ways to surpass its intrinsic catalytic activity. This study introduces a nanodomain tungsten oxide (WO3) interface to 1T-WS2, opening a new route for facilitating the transfer of a proton to active sites, thereby enhancing the HER performance. After H2S plasma sulfurization on the W layer to realize nanocrystalline 1T-WS2, subsequent O2 plasma treatment led to the formation of amorphous WO3 (a-WO3), resulting in a patchwork-structured heterointerface of 1T-WS2/a-WO3 (WSO). Addition of a hydrophilic interface (WO3) facilitates the hydrogen spillover effect, which represents the transfer of absorbed protons from a-WO3 to 1T-WS2. Moreover, the faster response of the cathodic current peak (proton insertion) in cyclic voltammetry is confirmed by the higher degree of oxidation. The rationale behind the faster proton insertion is that the introduced a-WO3 works as a proton channel. As a result, WSO-1.2 (the ratio of 1T-WS2 to a-WO3) exhibits a remarkable HER activity in that 1T-WS2 consumes more protons provided by the channel, showing an overpotential of 212 mV at 10 mA/cm2. Density functional theory calculations also show that the WO3 phase gives higher binding energies for initial proton adsorption, while the 1T-WS2 phase shows reduced HER overpotential. This improved catalytic performance demonstrates a novel strategy for water splitting to actively elicit the related reaction via efficient proton transport.

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.

Membrane Insertion of MoS2 Nanosheets: Fresh vs. Aged

Fresh two-dimensional molybdenum disulfide (MoS₂) absorbs the hydrocarbon contaminations in the ambient air and makes surface aging. To understand how the surface aging influences the interactions between MoS₂ and biomolecules is important in the biomedical applications. Here, employing all-atom molecular dynamics simulations, we investigated the interactions of the fresh and aged MoS₂ nanosheets with the lipid membranes of different components. Our results demonstrate that both the fresh and aged MoS₂ nanosheets can spontaneously insert into the bilayer membranes. However, the fresh MoS₂ nanosheet displays significantly stronger interaction and then has a larger penetration depth than the aged counterpart, regardless of the lipid components. The calculations of potential mean forces through the umbrella sampling further confirm that the insertion of fresh MoS₂ into the lipid membranes is more energetically favorable. Moreover, we found that the fresh MoS₂ nanosheet can cause a larger damage to the integrity of lipid membranes than the aged one. This work provides insightful understandings of the surface-aging-dependent interactions of the MoS₂ nanosheets with biomembranes, which could facilitate the design of novel MoS₂-based nanodevices with advanced surface properties.

Water friction in nanofluidic channels made from two-dimensional crystals

Membrane-based applications such as osmotic power generation, desalination and molecular separation would benefit from decreasing water friction in nanoscale channels. However, mechanisms that allow fast water flows are not fully understood yet. Here we report angstrom-scale capillaries made from atomically flat crystals and study the effect of confining walls' material on water friction. A massive difference is observed between channels made from isostructural graphite and hexagonal boron nitride, which is attributed to different electrostatic and chemical interactions at the solid-liquid interface. Using precision microgravimetry and ion streaming measurements, we evaluate the slip length, a measure of water friction, and investigate its possible links with electrical conductivity, wettability, surface charge and polarity of the confining walls. We also show that water friction can be controlled using hybrid capillaries with different slip lengths at opposing walls. The reported advances extend nanofluidics' toolkit for designing smart membranes and mimicking manifold machinery of biological channels.

Does Expanding or Contracting MgO Lattice Really Help with Corrosion Resistance of Mg Surface: Insights from Molecular Dynamics Simulations

In a humid environment, water droplets on the solid surface can act as a medium to accelerate corrosion. If the solid material has hydrophobic properties, the surface of the material will remain "clean" and corrosion may be retarded to a certain extent. In theory, MgO itself is a hydrophilic material, and we can apply additional stress or strain to change its lattice constant and adjust the wetting behavior of water on the MgO surface, resulting in changes of corrosion resistance. In order to study the effects of MgO lattice expansion or contraction on the wetting behavior of nano-water, molecular dynamics simulations have been performed in this work. It is found that the changes of the lattice constants on the MgO surface can significantly change the wetting tendency. It will alter the interaction forces between water molecules and MgO surfaces, which in turn changes the atomic density profiles, the orientation of OH bonds, and hydrogen bond networks. The contraction of MgO can actually result in the increase of wetting angles of nano-water droplets on the MgO surface and gradually exhibits hydrophobic properties.

Electrochemical Determination of Hydroxyurea in a Complex Biological Matrix Using MoS2-Modified Electrodes and Chemometrics

Hydroxyurea, an oral medication with important clinical benefits in the treatment of sickle cell anemia, can be accurately determined in plasma with a transition metal dichalcogenide-based electrochemical sensor. We used a two-dimensional molybdenum sulfide material (MoS₂) selectively electrodeposited on a polycrystalline gold electrode via tailored waveform polarization in the gold electrical double layer formation region. The electro-activity of the modified electrode depends on the electrical waveform parameters used to electro-deposit MoS₂. The concomitant oxidation of the MoS₂ material during its electrodeposition allows for the tuning of the sensor’s specificity. Chemometrics, utilizing mathematical procedures such as principal component analysis and multivariable partial least square regression, were used to process the electrochemical data generated at the bare and the modified electrodes, thus allowing the hydroxyurea concentrations to be predicted in human plasma. A limit-of-detection of 22 nM and a sensitivity of 37 nA cm⁻² µM⁻¹ were found to be suitable for pharmaceutical and clinical applications.

Relating Selectivity and Separation Performance of Lamellar Two-Dimensional Molybdenum Disulfide (MoS2) Membranes to Nanosheet Stacking Behavior

Increased demand for highly selective and energy-efficient separations processes has stimulated substantial interest in emerging two-dimensional (2D) nanomaterials as a potential platform for next-generation membranes. However, persistently poor separation performance continues to hinder the viability of many novel 2D-nanosheet membranes in desalination applications. In this study, we examine the role of the lamellar structure of 2D membranes on their performance. Using self-fabricated molybdenum disulfide (MoS₂) membranes as a platform, we show that the separation layer of 2D nanosheet frameworks not only fails to demonstrate water-salt selectivity but also exhibits low rejection toward dye molecules. Moreover, the MoS₂ membranes possess a molecular weight cutoff comparable to its underlying porous support, implying negligible selectivity of the MoS₂ layer. By tuning the nanochannel size through intercalation with amphiphilic molecules and analyzing mass transport in the lamellar structure using Monte Carlo simulations, we reveal that small imperfections in the stacking of MoS₂ nanosheets result in the formation of catastrophic microporous defects. These defects lead to a precipitous reduction in the selectivity of the lamellar structure by negating the interlayer sieving mechanism that prevents the passage of large penetrants. Notably, the imperfect stacking of nanosheets in the MoS₂ membrane was further verified using 2D X-ray diffraction measurements. We conclude that developing a well-controlled fabrication process, in which the lamellar structure can be carefully tuned, is critical to achieving defect-free and highly selective 2D desalination membranes.

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

Despite relevance to water purification and renewable energy conversion membranes, the molecular mechanisms underlying water slip are poorly understood. We disentangle the static and dynamical origin of water slippage on graphene, hBN and MoS2 by means of large-scale ab initio molecular dynamics. Accounting for the role of the electronic structure of the interface is essential to determine that water slips five and eleven times faster on graphene compared to hBN and to MoS2, respectively. Intricate changes in the water energy landscape as well as in the density correlations of the fluid provide, respectively, the main static and dynamical origin of water slippage. Surprisingly, the timescales of the density correlations are the same on graphene and hBN, whereas they are longer on MoS2 and yield a 100% slowdown in the flow of water on this material. Our results pave the way for an in silico first principles design of materials with enhanced water slip, through the modification of properties connected not only to the structure, but also to the dynamics of the interface.

Parameterization of Molybdenum Disulfide Interacting with Water Using the Free Energy Perturbation Method

Water contact angles (WCA) are often used to parametrize force field parameters of novel 2D nanomaterials, such as molybdenum disulfide (MoS₂), which has emerged as a promising nanomaterial in many biomedical applications due to its unique and impressive properties. However, there is a wide range of water-MoS₂ contact angles in the literature depending on the aging process on the surface of a MoS₂ nanosheet and/or substrate material. In this study, we revisit and optimize existing parameters for the basal plane of MoS₂ with two popular water models, TIP3P and SPC/E, using the wide range of WCAs from various experiments. We develop and deploy the free energy perturbation method for parametrizing MoS₂ with experimentally determined WCAs for both fresh and aged surfaces. Energy decomposition analysis on the simulation trajectories reveals that MoS₂-water interaction is dominated by van der Waals interaction, which mainly comes from the top layer of MoS₂. We conclude that to describe both fresh and aged MoS₂ surfaces it is convenient to only adjust the Lennard-Jones parameter ε_S (the depth of the potential well of a sulfur atom), which displays a surprisingly linear correlation with WCAs.

Determining the Surface Tension of Two-Dimensional Nanosheets by a Low-Rate Advancing Contact Angle Measurement

Because of their atomic thinness, two-dimensional (2D) nanosheets need be bound to a substrate or be dispersed in material in various applications. The surface tension (ST) of a 2D nanosheet is critical for analyzing the physicochemical interactions between 2D nanosheets and other materials. To date, the determination of the ST of 2D nanosheets has relied mainly on the contact angle (CA) method. However, because of the difficulty in measuring the thermodynamically significant Young?s CA, which is the only meaningful CA that can be used to determine the ST, significant differences exist in reported STs of 2D nanosheets. In this study, we obtained such unique Young?s CAs on graphene, boron nitride, molybdenum disulfide, and tungsten disulfide nanosheets by a low-rate advancing contact angle measurement using a rigorously designed experimental setup. By interpreting the CA with Neumann?s equation of state, we determined the STs of these four nanosheets to be 29.7 ? 0.6, 30.9 ? 0.7, 27.8 ? 0.7, and 29.1 ? 0.8 mJ/m², respectively. The surface energies of these 2D nanosheets were estimated to be in the range 95?120 mJ/m² by considering the contribution of ST and surface entropy. The accuracy of these determined STs was validated by the exfoliation and dispersion of 2D nanosheets in liquids with a series of STs. The study may have important implications for understanding the physicochemical interactions between 2D nanosheets and other materials and the development of 2D nanosheet-based devices.

Temperature regulation of the contact angle of water droplets on the solid surfaces

We investigate theoretically the stability of the wetting property, i.e., the contact angle values, as a function of the temperature. We find that the estimated temperature coefficient of the contact angle for the water droplets on an ordered water monolayer on a 100 surface of face-center cubic (FCC) is about one order of magnitude larger than that on a hydrophobic hexagonal surface in the temperature range between 290 K and 350 K, using molecular dynamics simulations. As temperature rises, the number of hydrogen bonds between the ordered water monolayer and the water droplet will increase, which therefore enhances the hydrophilicity of the ordered water monolayer at the FCC model surface. Our work thus provides an easily controllable and reversible way to control the degree of hydrophobicity of various solid surfaces exhibiting a similar wetting property of water droplets on the ordered water monolayer as such particular FCC (100) surfaces.

Liquids with Lower Wettability Can Exhibit Higher Friction on Hexagonal Boron Nitride: The Intriguing Role of Solid-Liquid Electrostatic Interactions

We investigate the wetting and frictional behavior of polar (water and ethylene glycol) and nonpolar (diiodomethane) liquids on the basal plane of hexagonal boron nitride (hBN) using molecular dynamics simulations. Our results for the wettability of water on the hBN basal plane (contact angle 81°) are in qualitative agreement with the experimentally deduced mild hydrophilicity of the hBN basal plane (contact angle 66°). We find that water exhibits the lowest wettability, as quantified by the highest contact angle, but the highest friction coefficient of (1.9 ± 0.4) × 10⁵ N-s/m³ on the hBN basal plane among the three liquids considered. This intriguing finding is explained in terms of the competition between dispersion and electrostatic interactions operating between the hBN basal plane and the three liquids. We find that electrostatic interactions do not affect the wetting behavior appreciably, as quantified by a less than 3° change in the respective contact angles of the three liquids considered. On the other hand, electrostatic interactions are found to increase the friction coefficients of the three liquids in contact with hBN to different extents, indicating that despite the increased friction of water on hBN, relative to that on graphene, nonpolar liquids may exhibit similar friction coefficients on hBN and graphene. Our findings reveal that the increase in the friction coefficient, upon incorporation of solid-liquid electrostatic interactions, is brought about by a greater increase in the solid-liquid mean-squared total lateral force, as compared to a smaller reduction in the decorrelation time of the solid-liquid force.

Exfoliation of Two-Dimensional Materials: The Role of Entropy

Liquid-phase exfoliation (LPE) is the best-known method for the synthesis of two-dimensional (2D) nanosheets. Compared to enthalpy, entropy is hardly considered to be a factor in choosing energy-efficient solvents and has not even been verified to be negligible. In this Letter, we explore the entropy contribution in LPE by performing molecular dynamics (MD) simulation of the structural flexibility effect in graphene, hexagonal boron nitride (hBN), and molybdenum disulfide (MoS₂). Our results show that surface vibration favors the exfoliation of graphene and hBN and destabilizes the reaggregation of nanosheets in water at 300 K, whereas the opposite is found for MoS₂. The entropy change is found to be 41%, 48%, and 4% of the enthalpy gain for graphene, hBN, and MoS₂ in LPE, respectively, and 64%, 32%, and 56% in reaggregation, which amounts to a step advancement for solvent screening in LPE of 2D materials.

Formation and stability of water clusters at the molybdenum disulfide interface: a molecular dynamics simulation investigation

In this work structural properties and dynamic behavior of a water nano droplet on the molybdenum disulfide were considered. The simulation results show that water molecules form polygon clusters on the interface, and most of which are hexagonal. Structures of water clusters at the interface are seen in two forms of curved and flattened polygons, which result in the formation of hydrogen bonds between and in the adjacent layers, respectively. Most of the clusters have circular flattened structures. Calculations of the lifetime of hydrogen bonds of water molecules at the interface also show that hydrogen bonds between water molecules at the interface have a low stability. This leads to the permanent formation and breaking down of hydrogen bonds of water molecules which can cause movement of water molecules and, consequently, the displacement of the center of mass and droplet motion. Considering the changes in the center of mass of a water droplet at the MoS₂ interface display, the water droplet has a significant spontaneous motion.

Black Phosphorus: Degradation Favors Lubrication

Due to its innate instability, the degradation of black phosphorus (BP) with oxygen and moisture was considered the obstacle for its application in ambient conditions. Here, a friction force reduced by about 50% at the degraded area of the BP nanosheets was expressly observed using atomic force microscopy due to the produced phosphorus oxides during degradation. Energy-dispersive spectrometer mapping analyses corroborated the localized concentration of oxygen on the degraded BP flake surface where friction reduction was observed. Water absorption was discovered to be essential for the degraded characteristic as well as the friction reduction behavior of BP sheets. The combination of water molecules as well as the resulting chemical groups (P-OH bonds) that are formed on the oxidized surface may account for the friction reduction of degraded BP flakes. It is indicated that, besides its layered structure, the ambient degradation of BP significantly favors its lubrication behavior.

Adsorption of common solvent molecules on graphene and MoS2 from first-principles

Solvents are an essential element in the production and processing of two-dimensional (2D) materials. For example, the liquid-phase exfoliation of layered materials requires a solvent to prevent the resulting monolayers from re-aggregating, while solutions of functional atoms and molecules are routinely used to modify the properties of the layers. It is generally assumed that these solvents do not interact strongly with the layer and so their effects can be neglected. Yet experimental evidence has suggested that explicit atomic-scale interactions between the solvent and layered material may play a crucial role in exfoliation and cause unintended electronic changes in the layer. Little is known about the precise nature of the interaction between the solvent molecules and the 2D layer. Here, we use density functional theory calculations to determine the adsorption configuration and binding energy of a variety of common solvent molecules, both polar and non-polar, on two of the most popular 2D materials, namely, graphene and MoS₂. We show that these molecules are physisorbed on the surface with negligible charge transferred between them. We find that the adsorption strength of the different molecules is independent of the polar nature of the solvent. However, we show that the molecules induce a significant charge rearrangement at the interface after adsorption as a result of polar bonds in the molecule.

Anisotropic Polymer Adsorption on Molybdenite Basal and Edge Surfaces and Interaction Mechanism With Air Bubbles

The anisotropic surface characteristics and interaction mechanisms of molybdenite (MoS₂) basal and edge planes have attracted much research interest in many interfacial processes such as froth flotation. In this work, the adsorption of a polymer depressant [i.e., carboxymethyl cellulose (CMC)] on both MoS₂ basal and edge surfaces as well as their interaction mechanisms with air bubbles have been characterized by atomic force microscope (AFM) imaging and quantitative force measurements. AFM imaging showed that the polymer coverage on the basal plane increased with elevating polymer concentration, with the formation of a compact polymer layer at 100 ppm CMC; however, the polymer adsorption was much weaker on the edge plane. The anisotropy in polymer adsorption on MoS₂ basal and edge surfaces coincided with water contact angle results. Direct force measurements using CMC functionalized AFM tips revealed that the adhesion on the basal plane was about an order of magnitude higher than that on the edge plane, supporting the anisotropic CMC adsorption behaviors. Such adhesion difference could be attributed to their difference in surface hydrophobicity and surface charge, with weakened hydrophobic attraction and strengthened electrostatic repulsion between the polymers and edge plane. Force measurements using a bubble probe AFM showed that air bubble could attach to the basal plane during approach, which could be effectively inhibited after polymer adsorption. The edge surface, due to the negligible polymer adsorption, showed similar interaction behaviors with air bubbles before and after polymer treatment. This work provides useful information on the adsorption of polymers on MoS₂ basal/edge surfaces as well as their interaction mechanism with air bubbles at the nanoscale, with implications for the design and development of effective polymer additives to mediate the bubble attachment on solid particles with anisotropic surface properties in mineral flotation and other engineering processes.

Water wettability of graphene: interplay between the interfacial water structure and the electronic structure

Wetting phenomena are ubiquitous and impact a wide range of applications. Simulations so far have largely relied on classical potentials. Here, we report the development of an approach that combines density-functional theory (DFT)-based calculations with classical wetting theory that allows practical but sufficiently accurate determination of the water contact angle (WCA). As a benchmark, we apply the approach to the graphene and graphite surfaces that recently received considerable attention. The results agree with and elucidate the experimental data. For metal-supported graphene where electronic interactions play a major role, we demonstrate that doping of graphene by the metal substrate significantly alters the wettability. In addition to theory, we report new experimental measurements of the WCA and the force of adhesion that corroborate the theoretical results. We demonstrate a correlation between the force of adhesion and WCA, and the use of the atomic force microscope (AFM) technique as an alternative measure for wettability at the nanoscale. The present work not only provides a detailed understanding of the wettability of graphene, including the role of electrons, but also sets the stage for studying the wettability alteration mechanism when sufficiently accurate force fields may not be available.

Ab Initio Molecular Dynamics and Lattice Dynamics-Based Force Field for Modeling Hexagonal Boron Nitride in Mechanical and Interfacial Applications

Hexagonal boron nitride (hBN) is an up-and-coming two-dimensional material, with applications in electronic devices, tribology, and separation membranes. Herein, we utilize density-functional-theory-based ab initio molecular dynamics (MD) simulations and lattice dynamics calculations to develop a classical force field (FF) for modeling hBN. The FF predicts the crystal structure, elastic constants, and phonon dispersion relation of hBN with good accuracy and exhibits remarkable agreement with the interlayer binding energy predicted by random phase approximation calculations. We demonstrate the importance of including Coulombic interactions but excluding 1-4 intrasheet interactions to obtain the correct phonon dispersion relation. We find that improper dihedrals do not modify the bulk mechanical properties and the extent of thermal vibrations in hBN, although they impact its flexural rigidity. Combining the FF with the accurate TIP4P/Ice water model yields excellent agreement with interaction energies predicted by quantum Monte Carlo calculations. Our FF should enable an accurate description of hBN interfaces in classical MD simulations.

Under-Oil Switchable Superhydrophobicity to Superhydrophilicity Transition on TiO2 Nanotube Arrays

Recently, smart interfacial materials that can reversibly transit between the superhydrophobicity and superhydrophilicity have aroused much attention. However, all present performances happen in air, and to realize such a smart transition in complex environments, such as oil, is still a challenge. Herein, TiO₂ nanotube arrays with switchable transition between the superhydrophobicity and superhydrophilicity in oil are reported. The switching can be observed by alternation of UV irradiation and heating process, and the smart controllability can be ascribed to the cooperative effect between the surface nanostructures and the chemical composition variation. By using the controllable wetting performances, some applications such as under-oil droplet-based microreaction and water-removal from oil were demonstrated on our surface. This paper reports a surface with smart water wettability in oil, which could start some fresh ideas for wetting control on interfacial materials.

Exfoliation of graphene and fluorographene in molecular and ionic liquids

We use molecular dynamics simulations to study the exfoliation of graphene and fluorographene in molecular and ionic liquids using computer experiments in which one layer of the 2D nanomaterial is peeled, in vacuum or with solvent present.

Doping-Driven Wettability of Two-Dimensional Materials: A Multiscale Theory

Engineering molecular interactions at two-dimensional (2D) materials interfaces enables new technological opportunities in functional surfaces and molecular epitaxy. Understanding the wettability of 2D materials represents the crucial first step toward quantifying the interplay between the interfacial forces and electric potential of 2D materials interfaces. Here we develop the first theoretical framework to model the wettability of the doped 2D materials by properly bridging the multiscale physical phenomena at the 2D interfaces, including (i) the change of 2D materials surface energy (atomistic scale, several angstroms), (ii) the molecular reorientation of liquid molecules adjacent to the interface (molecular scale, 10⁰-10¹ nm), and (iii) the electrical double layer (EDL) formed in the liquid phase (mesoscopic scales, 10⁰-10⁴ nm). The latter two effects are found to be the major mechanisms responsible for the contact angle change upon doping, while the surface energy change of a pure 2D material has no net effect on the wetting property. When the doping level is electrostatically tuned, we demonstrate that 2D materials with high quantum capacitances (e.g., transition metal dichalcogenides, TMDCs) possess a wider range of tunability in the interfacial tension, under the same applied gate voltage. Furthermore, practical considerations such as defects and airborne contamination are also quantitatively discussed. Our analysis implies that the doping level can be another variable to modulate the wettability at 2D materials interfaces, as well as the molecular packing behavior on a 2D material-coated surface, essentially facilitating the interfacial engineering of 2D materials.

Imidazolium Ionic Liquid Mediates Black Phosphorus Exfoliation while Preventing Phosphorene Decomposition

Forthcoming applications in electronics and optoelectronics make phosphorene a subject of vigorous research efforts. Solvent-assisted exfoliation of phosphorene promises affordable delivery in industrial quantities for future applications. We demonstrate, using equilibrium, steered and umbrella sampling molecular dynamics, that the 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIM][BF₄] ionic liquid is an excellent solvent for phosphorene exfoliation. The presence of both hydrophobic and hydrophilic moieties, as well as substantial shear viscosity, allows [EMIM][BF₄] simultaneously to facilitate separation of phosphorene sheets and to protect them from getting in direct contact with moisture and oxygen. The exfoliation thermodynamics is moderately unfavorable, which indicates that an external stimulus is necessary. Unexpectedly, [EMIM][BF₄] does not coordinates phosphorene by π-electron stacking with the imidazole ring. Instead, the solvation proceeds via hydrophobic side chains, while polar imidazole rings form an electrostatically stabilized protective layer. The simulations suggest that further efforts in solvent engineering for phosphorene exfoliation should concentrate on use of weakly coordinating ions and grafting groups that promote stronger dispersion interactions and on elongation of nonpolar chains.

Interaction Mechanisms between Air Bubble and Molybdenite Surface: Impact of Solution Salinity and Polymer Adsorption

The surface characteristics of molybdenite (MoS₂) such as wettability and surface interactions have attracted much research interest in a wide range of engineering applications, such as froth flotation. In this work, a bubble probe atomic force microscope (AFM) technique was employed to directly measure the interaction forces between an air bubble and molybdenite mineral surface before/after polymer (i.e., guar gum) adsorption treatment. The AFM imaging showed that the polymer coverage on the surface of molybdenite could achieve ∼5.6, ∼44.5, and ∼100% after conditioning in 1, 5, and 10 ppm polymer solution, respectively, which coincided with the polymer coverage results based on contact angle measurements. The electrolyte concentration and surface treatment by polymer adsorption were found to significantly affect bubble-mineral interaction and attachment. The experimental force results on bubble-molybdenite (without polymer treatment) agreed well with the calculations using a theoretical model based on the Reynolds lubrication theory and augmented Young-Laplace equation including the effect of disjoining pressure. The overall surface repulsion was enhanced when the NaCl concentration decreased from 100 to 1 mM, which inhibited the bubble-molybdenite attachment. After conditioning the molybdenite surface in 1 ppm polymer solution, it was more difficult for air bubbles to attach to the molybdenite surface due to the weakened hydrophobic interaction with a shorter decay length. Increasing the polymer concentration to 5 ppm effectively inhibited bubble attachment on mineral surface, which was mainly due to the much reduced hydrophobic interaction as well as the additional steric repulsion between the extended polymer chains and bubble surface. The results provide quantitative information on the interaction mechanism between air bubbles and molybdenite mineral surfaces on the nanoscale, with useful implications for the development of effective polymer depressants and fundamental understanding of bubble-solid interactions in mineral flotation. The methodologies used in this work can be readily extended to studying similar interfacial interactions in many other engineering applications such as froth flotation deinking and bitumen extraction in oil sands industry.

Revisiting the droplet simulation approach to derive force-field parameters for water on molybdenum disulfide from wetting angle measurements

Owing to its peculiar electronic properties, molybdenum disulfide (MoS₂) has been the subject of a growing number of studies in the recent years. In applications, this material and other transition metal dichalcogenides (TMDs) may have to interact with a liquid or polymer phase as well as solutions of biomolecules. It is therefore of primary importance to understand the wetting and adhesion properties of TMDs. Starting from existing models, we derive Lennard-Jones parameters for the interaction between water and the basal plane of MoS₂ that are consistent with recent wetting experiments. Molecular dynamics simulations indicate that a stack of only two MoS₂ monolayers is necessary to capture the wetting behavior of bulk MoS₂. It is found that the Coulomb interaction between water and monolayer and bilayer MoS₂ plays no role in the related interfacial thermodynamics. Calculations with the optimized parameters show that the depth of the well of the interaction potential between water and bulk MoS₂ is of the order of 8.2 kJ/mol. Such a value is comparable with what was found for graphite and consistent with the fact that the wetting angles of water on graphite and MoS₂ are almost equal. The derivation of the force-field parameters is performed using a methodology which, contrary to previous studies, makes a consistent use of droplet calculations. The results of our work should find application in further simulation studies on the wetting behavior of TMDs and other dispersive materials.