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.

Surfactant-Aided Stabilization of Individual Carbon Nanotubes in Water around the Critical Micelle Concentration.

Surfactants are widely used to disperse single-walled carbon nanotubes (SWCNTs) and other nanomaterials for liquid-phase processing and characterization. Traditional techniques, however, demand high surfactant concentrations, often in the range of 1-2 wt/v% of the solution. Here, we show that optimal dispersion efficiency can be attained at substantially lower surfactant concentrations of approximately 0.08 wt/v%, near the critical micelle concentration. This unexpected observation is achieved by introducing "bare" nanotubes into water containing the anionic surfactant sodium deoxycholate (DOC) through a superacid-surfactant exchange process that eliminates the need for ultrasonication. Among the diverse ionic surfactants and charged biopolymers explored, DOC exhibits the highest dispersion efficiency, outperforming sodium cholate, a structurally similar bile salt surfactant containing just one additional oxygen atom compared to DOC. Employing all-atomistic molecular dynamics simulations, we unravel that the greater stabilization by DOC arises from its higher binding affinity to nanotubes and a substantially larger free energy barrier that resists nanotube rebundling. Further, we find that this barrier is nonelectrostatic in nature and does not obey the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloidal stability, underscoring the important role of nonelectrostatic dispersion and hydration interactions at the nanoscale, even in the case of ionic surfactants like DOC. These molecular insights advance our understanding of surfactant chemistry at the bare nanotube limit and suggest low-energy, surfactant-efficient solution processing of SWCNTs and potentially other nanomaterials.

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.

Discretized hexagonal boron nitride quantum emitters and their chemical interconversion.

Quantum emitters in two-dimensional hexagonal boron nitride (hBN) are of significant interest because of their unique photophysical properties, such as single-photon emission at room temperature, and promising applications in quantum computing and communications. The photoemission from hBN defects covers a wide range of emission energies but identifying and modulating the properties of specific emitters remain challenging due to uncontrolled formation of hBN defects. In this study, more than 2000 spectra are collected consisting of single, isolated zero-phonon lines (ZPLs) between 1.59 and 2.25 eV from diverse sample types. Most of ZPLs are organized into seven discretized emission energies. All emitters exhibit a range of lifetimes from 1 to 6 ns, and phonon sidebands offset by the dominant lattice phonon in hBN near 1370 cm-1. Two chemical processing schemes are developed based on water and boric acid etching that generate or preferentially interconvert specific emitters, respectively. The identification and chemical interconversion of these discretized emitters should significantly advance the understanding of solid-state chemistry and photophysics of hBN quantum emission.

Breakdown of the Nernst-Einstein relation in carbon nanotube porins.

For over 100 years, the Nernst-Einstein relation has linked a charged particle's electrophoretic mobility and diffusion coefficient. Here we report experimental measurements of diffusion and electromigration of K+ ions in narrow 0.8-nm-diameter single-walled carbon nanotube porins (CNTPs) and demonstrate that the Nernst-Einstein relation in these channels breaks down by more than three orders of magnitude. Molecular dynamics simulations using polarizable force fields show that K+ ion diffusion in CNTPs in the presence of a single-file water chain is three orders of magnitude slower than bulk diffusion. Intriguingly, the simulations also reveal a disintegration of the water chain upon application of electric fields, resulting in the formation of distinct K+-water clusters, which then traverse the CNTP at high velocity. Finally, we show that although individual ion-water clusters still obey the Nernst-Einstein relation, the overall relation breaks down because of two distinct mechanisms for ion diffusion and electromigration.

Gas Separations using Nanoporous Atomically Thin Membranes: Recent Theoretical, Simulation, and Experimental Advances.

Porous graphene and other atomically thin 2D materials are regarded as highly promising membrane materials for high-performance gas separations due to their atomic thickness, large-scale synthesizability, excellent mechanical strength, and chemical stability. When these atomically thin materials contain a high areal density of gas-sieving nanoscale pores, they can exhibit both high gas permeances and high selectivities, which is beneficial for reducing the cost of gas-separation processes. Here, recent modeling and experimental advances in nanoporous atomically thin membranes for gas separations is discussed. The major challenges involved, including controlling pore size distributions, scaling up the membrane area, and matching theory with experimental results, are also highlighted. Finally, important future directions are proposed for real gas-separation applications of nanoporous atomically thin membranes.

Direct Chemical Vapor Deposition Synthesis of Porous Single-Layer Graphene Membranes with High Gas Permeances and Selectivities.

Single-layer graphene containing molecular-sized in-plane pores is regarded as a promising membrane material for high-performance gas separations due to its atomic thickness and low gas transport resistance. However, typical etching-based pore generation methods cannot decouple pore nucleation and pore growth, resulting in a trade-off between high areal pore density and high selectivity. In contrast, intrinsic pores in graphene formed during chemical vapor deposition are not created by etching. Therefore, intrinsically porous graphene can exhibit high pore density while maintaining its gas selectivity. In this work, the density of intrinsic graphene pores is systematically controlled for the first time, while appropriate pore sizes for gas sieving are precisely maintained. As a result, single-layer graphene membranes with the highest H2 /CH4 separation performances recorded to date (H2 permeance > 4000 GPU and H2 /CH4 selectivity > 2000) are fabricated by manipulating growth temperature, precursor concentration, and non-covalent decoration of the graphene surface. Moreover, it is identified that nanoscale molecular fouling of the graphene surface during gas separation where graphene pores are partially blocked by hydrocarbon contaminants under experimental conditions, controls both selectivity and temperature dependent permeance. Overall, the direct synthesis of porous single-layer graphene exploits its tremendous potential as high-performance gas-sieving membranes.

Predicting Gas Separation through Graphene Nanopore Ensembles with Realistic Pore Size Distributions.

The development of nanoporous single-layer graphene membranes for gas separation has prompted increasing theoretical investigations of gas transport through graphene nanopores. However, computer simulations and theories that predict gas permeances through individual graphene nanopores are not suitable to describe experimental results, because a realistic graphene membrane contains a large number of nanopores of diverse sizes and shapes. With this need in mind, here, we generate nanopore ensembles in silico by etching carbon atoms away from pristine graphene with different etching times, using a kinetic Monte Carlo algorithm developed by our group for the isomer cataloging problem of graphene nanopores. The permeances of H2, CO2, and CH4 through each nanopore in the ensembles are predicted using transition state theory based on classical all-atomistic force fields. Our findings show that the total gas permeance through a nanopore ensemble is dominated by a small fraction of large nanopores with low energy barriers of pore crossing. We also quantitatively predict the increase of the gas permeances and the decrease of the selectivities between the gases as functions of the etching time of graphene. Furthermore, by fitting the theoretically predicted selectivities to the experimental ones reported in the literature, we show that nanopores in graphene effectively expand as the temperature of permeation measurement increases. We propose that this nanopore "expansion" is due to the desorption of contaminants that partially clog the graphene nanopores. In general, our study highlights the effects of the pore size and shape distributions of a graphene nanopore ensemble on its gas separation properties and calls into attention the potential effect of pore-clogging contamination in experiments.

Uncovering a Universal Molecular Mechanism of Salt Ion Adsorption at Solid/Water Interfaces.

Solid/water interfaces, in which salt ions come in close proximity to solids, are ubiquitous in nature. Because water is a polar solvent and salt ions are charged, a long-standing puzzle involving solid/water interfaces is how do the electric fields exerted by the salt ions and the interfacial water molecules polarize the charge distribution in the solid and how does this polarization, in turn, influence ion adsorption at any solid/water interface. Here, using state-of-the-art polarizable force fields derived from quantum chemical simulations, we perform all-atomistic molecular dynamics simulations to investigate the adsorption of various ions comprising the well-known Hofmeister series at the graphene/water interface, including comparing with available experimental data. Our findings reveal that, in vacuum, the ionic electric field-induced polarization of graphene results in a significantly large graphene-ion polarization energy, which drives all salt ions to adsorb to graphene. On the contrary, in the presence of water molecules, we show that the ions and the water molecules exert waves of molecular electric fields on graphene which destructively interfere with each other. This remarkable phenomenon is shown to cause a water-mediated screening of more than 85% of the graphene-ion polarization energy. Finally, by investigating superhydrophilic and superhydrophobic model surfaces, we demonstrate that this phenomenon occurs universally at all solid/water interfaces and results in a significant weakening of the ion-solid interactions, such that ion specific effects are governed primarily by a competition between the ion-water and water-water interactions, irrespective of the nature of the solid/water interface.

Diameter Dependence of Water Filling in Lithographically Segmented Isolated Carbon Nanotubes.

Although the structure and properties of water under conditions of extreme confinement are fundamentally important for a variety of applications, they remain poorly understood, especially for dimensions less than 2 nm. This problem is confounded by the difficulty in controlling surface roughness and dimensionality in fabricated nanochannels, contributing to a dearth of experimental platforms capable of carrying out the necessary precision measurements. In this work, we utilize an experimental platform based on the interior of lithographically segmented, isolated single-walled carbon nanotubes to study water under extreme nanoscale confinement. This platform generates multiple copies of nanotubes with identical chirality, of diameters from 0.8 to 2.5 nm and lengths spanning 6 to 160 µm, that can be studied individually in real time before and after opening, exposure to water, and subsequent water filling. We demonstrate that, under controlled conditions, the diameter-dependent blue shift of the Raman radial breathing mode (RBM) between 1 and 8 cm-1 measures an increase in the interior mechanical modulus associated with liquid water filling, with no response from exterior water exposure. The observed RBM shift with filling demonstrates a non-monotonic trend with diameter, supporting the assignment of a minimum of 1.81 ± 0.09 cm-1 at 0.93 ± 0.08 nm with a nearly linear increase at larger diameters. We find that a simple hard-sphere model of water in the confined nanotube interior describes key features of the diameter-dependent modulus change of the carbon nanotube and supports previous observations in the literature. Longer segments of 160 µm show partial filling from their ends, consistent with pore clogging. These devices provide an opportunity to study fluid behavior under extreme confinement with high precision and repeatability.

Analytical Prediction of Gas Permeation through Graphene Nanopores of Varying Sizes: Understanding Transitions across Multiple Transport Regimes.

Nanoporous graphene is a promising candidate material for gas separation membranes, due to its atomic thickness and low cross-membrane transport resistance. The mechanisms of gas permeation through graphene nanopores, in both the large and small pore size limits, have been reported in the literature. However, mechanistic insights into the crossover from the small pore size limit to the large pore size limit are still lacking. In this study, we develop a comprehensive theoretical framework to predict gas permeance through graphene nanopores having a wide range of diameters using analytical equations. We formulate the transport kinetics associated with the direct impingement from the bulk and with the surface diffusion from the adsorption layer on graphene and then combine them to predict the overall gas permeation rate using a reaction network model. We also utilize molecular dynamics simulations to validate and calibrate our theoretical model. We show that the rates of both the direct impingement and the surface diffusion pathways need to be corrected using different multiplicative factors, which are functions of temperature, gas kinetic diameter, and pore diameter. Further, we find a minor spillover pathway that originates from the surface adsorption layer, but is not included in our theoretical model. Finally, we utilize the corrected model to predict the permeances of CO2, CH4, and Ar through graphene nanopores. We show that as the pore diameter increases, gas transport through graphene nanopores can transition from being translocation dominated (pore diameter < 0.7 nm), to surface pathway dominated (pore diameter 1-2 nm), and finally to direct pathway dominated (pore diameter > 4 nm). The various gas permeation mechanisms outlined in this study will be particularly useful for the rational design of membranes made out of two-dimensional materials such as graphene for gas separation applications.

Theory of Surface Forces in Multivalent Electrolytes.

Aqueous electrolyte solutions containing multivalent ions exhibit various intriguing properties, including attraction between like-charged colloidal particles, which results from strong ion-ion correlations. In contrast, the classical Derjaguin-Landau-Verwey-Overbeek theory of colloidal stability, based on the Poisson-Boltzmann mean-field theory, always predicts a repulsive electrostatic contribution to the disjoining pressure. Here, we formulate a general theory of surface forces, which predicts that the contribution to the disjoining pressure resulting from ion-ion correlations is always attractive and can readily dominate over entropic-induced repulsions for solutions containing multivalent ions, leading to the phenomenon of like-charge attraction. Ion-specific short-range hydration interactions, as well as surface charge regulation, are shown to play an important role at smaller separation distances but do not fundamentally change these trends. The theory is able to predict the experimentally observed strong cohesive forces reported in cement pastes, which result from strong ion-ion correlations involving the divalent calcium ion.

Multi-scale approach for modeling stability, aggregation, and network formation of nanoparticles suspended in aqueous solutions.

Suspensions of nanoparticles (NPs) in aqueous solutions hold promise in many research fields, including energy applications, water desalination, and nanomedicine. The ability to tune NP interactions, and thereby to modulate the NP self-assembly process, holds the key to rationally synthesize NP suspensions. However, traditional models obtained by coupling the DLVO (Derjaguin, Landau, Verwey, and Overbeek) theory of NP interactions, or suitable modifications of it, with the kinetic theory of colloidal aggregation are inadequate to precisely model NP self-assembly because they neglect hydration forces and discrete-size effects predominant at the nanoscale. By synergistically blending molecular dynamics and stochastic dynamics simulations with continuum theories, we develop a multi-scale (MS) model, which is able to accurately predict suspension stability, timescales for NP aggregation, and macroscopic properties (e.g., the thermal conductivity) of bare and surfactant-coated NP suspensions, in good agreement with the experimental data. Our results enable the formulation of design rules for engineering NP aqueous suspensions in a wide range of applications.

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) × 105 N-s/m3 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.

Addressing the isomer cataloguing problem for nanopores in two-dimensional materials.

The presence of extended defects or nanopores in two-dimensional (2D) materials can change the electronic, magnetic and barrier membrane properties of the materials. However, the large number of possible lattice isomers of nanopores makes their quantitative study a seemingly intractable problem, confounding the interpretation of experimental and simulated data. Here we formulate a solution to this isomer cataloguing problem (ICP), combining electronic-structure calculations, kinetic Monte Carlo simulations, and chemical graph theory, to generate a catalogue of unique, most-probable isomers of 2D lattice nanopores. The results demonstrate remarkable agreement with precise nanopore shapes observed experimentally in graphene and show that the thermodynamic stability of a nanopore is distinct from its kinetic stability. Triangular nanopores prevalent in hexagonal boron nitride are also predicted, extending this approach to other 2D lattices. The proposed method should accelerate the application of nanoporous 2D materials by establishing specific links between experiment and theory/simulations, and by providing a much-needed connection between molecular design and fabrication.

Stable, Temperature-Dependent Gas Mixture Permeation and Separation through Suspended Nanoporous Single-Layer Graphene Membranes.

Graphene membranes with nanometer-scale pores could exhibit an extremely high permeance and selectivity for the separation of gas mixtures. However, to date, no experimental measurements of gas mixture separation through nanoporous single-layer graphene (SLG) membranes have been reported. Herein, we report the first measurements of the temperature-dependent permeance of gas mixtures in an equimolar mixture feed containing H2, He, CH4, CO2, and SF6 from 22 to 208 °C through SLG membranes containing nanopores formed spontaneously during graphene synthesis. Five membranes were fabricated by transfer of CVD graphene from catalytic Cu film onto channels framed in impermeable Ni. Two membranes exhibited gas permeances on the order of 10-6 to 10-5 mol m-2 s-1 Pa-1 as well as gas mixture selectivities higher than the Knudsen effusion selectivities predicted by the gas effusion mechanism. We show that a new steric selectivity mechanism explains the permeance data and selectivities. This mechanism predicts a mean pore diameter of 2.5 nm and an areal pore density of 7.3 × 1013 m-2, which is validated by experimental observations. A third membrane exhibited selectivities lower than the Knudsen effusion selectivities, suggesting a combination of effusion and viscous flow. A fourth membrane exhibited increasing permeance values as functions of temperature from 27 to 200 °C, and a CO2/SF6 selectivity > 20 at 200 °C, suggestive of activated translocation through molecular-sized nanopores. A fifth membrane exhibited no measurable permeance of any gas above the detection limit of our technique, 2 × 10-7 mol m-2 s-1 Pa-1, indicating essentially a molecularly impermeable barrier. Overall, these data demonstrate that SLG membranes can potentially provide a high mixture separation selectivity for gases, with CVD synthesis alone resulting in nanometer-scale pores useful for gas separation. This work also shows that temperature-dependent permeance measurements on SLG can be used to reveal underlying permeation mechanisms.

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.

Molecular Rotors for Universal Quantitation of Nanoscale Hydrophobic Interfaces in Microplate Format.

Hydrophobic self-assembly pairs diverse chemical precursors and simple formulation processes to access a vast array of functional colloids. Exploration of this design space, however, is stymied by lack of broadly general, high-throughput colloid characterization tools. Here, we show that a narrow structural subset of fluorescent, zwitterionic molecular rotors, dialkylaminostilbazolium sulfonates [DASS] with intermediate-length alkyl tails, fills this major analytical void by quantitatively sensing hydrophobic interfaces in microplate format. DASS dyes supersede existing interfacial probes by avoiding off-target fluorogenic interactions and dye aggregation while preserving hydrophobic partitioning strength. To illustrate the generality of this approach, we demonstrate (i) a microplate-based technique for measuring mass concentration of small (20-200 nm), dilute (submicrogram sensitivity) drug delivery nanoparticles; (ii) elimination of particle size, surfactant chemistry, and throughput constraints on quantifying the complex surfactant/metal oxide adsorption isotherms critical for environmental remediation and enhanced oil recovery; and (iii) more reliable self-assembly onset quantitation for chemically and structurally distinct amphiphiles. These methods could streamline the development of nanotechnologies for a broad range of applications.

Schizophrenic Diblock-Copolymer-Functionalized Nanoparticles as Temperature-Responsive Pickering Emulsifiers.

Stimuli-responsive pickering emulsions have received considerable attention in recent years, and the utilization of temperature as a stimulus has been of particular interest. Previous efforts have led to responsive systems that enable the formation of stable emulsions at room temperature, which can subsequently be triggered to destabilize with an increase in temperature. The development of a thermoresponsive system that exhibits the opposite response, however, i.e., one that can be triggered to form stable emulsions at elevated temperatures and subsequently be induced to phase separate at lower temperatures, has so far been lacking. Here, we describe a system that accomplishes this goal by leveraging a schizophrenic diblock copolymer that exhibits both an upper and a lower critical solution temperature. The diblock copolymer was conjugated to 20 nm silica nanoparticles, which were subsequently demonstrated to stabilize O/W emulsions at 65 °C and trigger phase separation upon cooling to 25 °C. The effects of particle concentration, electrolyte concentration, and polymer architecture were investigated, and facile control of emulsion stability was demonstrated for multiple oil types. Our approach is likely to be broadly adaptable to other schizophrenic diblock copolymers and find significant utility in applications such as enhanced oil recovery and liquid-phase heterogeneous catalysis, where stable emulsions are desired only at elevated temperatures.