Osmotic Transport at the Aqueous Graphene and hBN Interfaces: Scaling Laws from a Unified, First-Principles Description

Interfacial Assembly of 2D Graphene-Derived Ion Channels for Water-Based Green Energy Conversion

The utilization of sustained and green energy is believed to alleviate increasing menace of global environmental concerns and energy dilemma. Interfacial assembly of 2D graphene-derived ion channels (2D-GDICs) with tunable ion/fluid transport behavior enables efficient harvesting of renewable green energy from ubiquitous water, especially for osmotic energy harvesting. In this review, various interfacial assembly strategies for fabricating diverse 2D-GDICs are summarized and their ion transport properties are discussed. This review analyzes how particular structure and charge density/distribution of 2D-GDIC can be modulated to minimize internal resistance of ion/fluid transport and enhance energy conversion efficiency, and highlights stimuli-responsive functions and stability of 2D-GDIC and further examines the possibility of integrating 2D-GDIC with other energy conversion systems. Notably, the presented preparation and applications of 2D-GDIC also inspire and guide other 2D materials to fabricate sophisticated ion channels for targeted applications. Finally, potential challenges in this field is analyzed and a prospect to future developments toward high-performance or large-scale real-word applications is offered.

Mono-component bacterial cellulose heterogeneous membrane mediated by ionic liquids for osmotic energy harvesting

The massive reserves of osmotic energy existing in estuary will be highly desired as promising energy source that avails to solve the problem of energy shortage and environment deterioration. The ion transport membrane is core component optimized through composite membrane heterostructure to maximize the osmotic energy harvesting but suffer from gaps and resistance increase, which limit their practical applications. Here we demonstrate mono-component heterogeneous regenerated bacterial cellulose (RBC) membranes fabricated by subtle regenerated technique through Ionic Liquids (ILs). Such membranes obtain heterogeneous nature by the difference in fiber intertwining states due to the different treatment conditions on both sides. It achieves osmotic energy conversion with maximum power density of 0.70 W·m-2at 100-fold, which provides ingenious strategy for excellent performance and low-cost osmotic energy harvesting. By minimizing pores and maximizing the surface charges, energy barriers can be lowered, ion permeable and selective transport channels for energy harvesting device can be increased, as supported by the numerical simulation. This is the first time the construction strategy for mono-component heterogeneous membrane mediated by ILs for osmotic energy harvesting is proposed, which averts gaps between the layers of different materials effectively and provides theoretical guidance for subsequent in-depth research on mono-component ion-selective heterogeneous membrane.

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.

First-Principles Prediction of Amorphous Silica Nanoparticle Surface Charge: Effect of Size, pH, and Ionic Strength

The quantification of surface charge properties of silica nanoparticles is essential for several applications. To determine these properties, many experimental and theoretical methods have been introduced, which are time-consuming and/or challenging to use. In this study, a first-principles approach is developed to determine the surface charge properties of amorphous silica nanoparticles against the nanoparticle size, pH, and ionic strength without relying on experimental data. An amorphous silica nanoparticle of 1.34 nm diameter is simulated by using integrated molecular dynamics and Monte Carlo methods. A detailed analysis of the nanoparticle structure is provided by analyzing the types of silanol groups on the surface. Moreover, a model is developed to estimate the probability distribution of the surface silanol groups based on the nearest neighbor distances and the diameter of the nanoparticle to determine the number of surface silanols on larger nanoparticles. Thereafter, a computational chemistry approach is used to calculate the acid dissociation constants of the corresponding deprotonation reactions. The calculated constants and the point of zero charge value are in excellent agreement with experiments. The surface charge properties of the nanoparticle with various diameters are then estimated by using a mean-field model at different pH and ionic strength values. The results of the developed model are compared to the Poisson-Boltzmann equation as a reference model. The developed model predictions agree well with the reference model for low and mid-electrolyte concentrations (1 and 10 mM) and small nanoparticles (smaller than 100 nm). However, the developed model seems to qualitatively predict the surface charge properties more accurately than the Poisson-Boltzmann model for high electrolyte concentrations.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.

Complex coupling between surface charge and thermo-osmotic phenomena

Thermo-osmotic flows, generated at liquid-solid interfaces by thermal gradients, can be used to produce electric currents from waste heat on charged surfaces. The two key parameters controlling the thermo-osmotic current are the surface charge and the interfacial enthalpy excess due to liquid-solid interactions. While it has been shown that the contribution from water to the enthalpy excess can be crucial, how this contribution is affected by surface charge remained to be understood. Here, we start by discussing how thermo-osmotic flows and induced electric currents are related to the interfacial enthalpy excess. We then use molecular dynamics simulations to investigate the impact of surface charge on the interfacial enthalpy excess, for different distributions of the surface charge, and two different wetting conditions. We observe that surface charge has a strong impact on enthalpy excess, and that the dependence of enthalpy excess on surface charge depends largely on its spatial distribution. In contrast, wetting has a very small impact on the charge-enthalpy coupling. We rationalize the results with simple analytical models, and explore their consequences for thermo-osmotic phenomena. Overall, this work provides guidelines to search for systems providing optimal waste heat recovery performance.

Coupled Interactions at the Ionic Graphene-Water Interface

We compute ionic free energy adsorption profiles at an aqueous graphene interface by developing a self-consistent approach. To do so, we design a microscopic model for water and put the liquid on an equal footing with the graphene described by its electronic band structure. By evaluating progressively the electronic and dipolar coupled electrostatic interactions, we show that the coupling level including mutual graphene and water screening permits one to recover remarkably the precision of extensive quantum simulations. We further derive the potential of mean force evolution of several alkali cations.

The dielectric function profile across the water interface through surface-specific vibrational spectroscopy and simulations

The dielectric properties of interfacial water on subnanometer length scales govern chemical reactions, carrier transfer, and ion transport at interfaces. Yet, the nature of the interfacial dielectric function has remained under debate as it is challenging to access the interfacial dielectric with subnanometer resolution. Here we use the vibrational response of interfacial water molecules probed using surface-specific sum-frequency generation (SFG) spectra to obtain exquisite depth resolution. Different responses originate from water molecules at different depths and report back on the local interfacial dielectric environment via their spectral amplitudes. From experimental and simulated SFG spectra at the air/water interface, we find that the interfacial dielectric constant changes drastically across an ∼1 Å thin interfacial water region. The strong gradient of the interfacial dielectric constant leads, at charged planar interfaces, to the formation of an electric triple layer that goes beyond the standard double-layer model.

Effects of surface rigidity and metallicity on dielectric properties and ion interactions at aqueous hydrophobic interfaces

Using classical molecular dynamics simulations we investigate the dielectric properties at interfaces of water with graphene, graphite, hexane and water vapor. For graphite we compare metallic and non-metallic versions. At the vapor-liquid water and hexane-water interfaces the laterally averaged dielectric profiles are significantly broadened due to interfacial roughness and only slightly anisotropic. In contrast, at the rigid graphene surface the dielectric profiles are strongly anisotropic and the perpendicular dielectric profile exhibits pronounced oscillations and sign changes. The interfacial dielectric excess, characterized by the shift of the dielectric-dividing-surface with respect to the Gibbs-dividing-surface, is positive for all surfaces, showing that water has an enhanced dielectric response at hydrophobic surfaces. The dielectric-dividing-surface positions vary significantly among the different surfaces, which points to pronounced surface-specific dielectric behavior. The interfacial repulsion of a chloride ion is shown to be dominated by electrostatic interactions for the soft fluid-fluid interfaces and by non-electrostatic Lennard-Jones interactions for the rigid graphene-water interface. A linear tensorial dielectric model for the ion-interface interaction with sharp dielectric interfaces located on the dielectric-dividing-surface positions works well for graphene but fails for vapor and hexane, because these interfaces are smeared out. The repulsion of chloride from the metallic and non-metallic graphite versions differs very little, which reflects the almost identical interfacial water structure and can be understood based on linear continuum dielectric theory. Interface flexibility shows up mostly in the non-linear Coulomb part of the ion-interface interaction, which changes significantly close to the interfaces and signals the breakdown of linear dielectric continuum theory.

Unraveling the Molecular Interface and Boundary Problems in an Electrical Double Layer and Electroosmotic Flow

In a nanofluidic system, the electroosmotic flow (EOF) is a complex fluid transport mechanism, where the formation of an electrical double layer (EDL) occurs ubiquitously at the dissimilar atomic interface. Several studies have suggested various interface boundaries to calculate the EDL thickness. However, the physical origin of the interface boundary and its effects on the flow properties is not yet clearly understood. Combining the theoretical framework and molecular dynamics (MD) simulations, we show the effects of different interfacial boundaries on the EDL thickness and EOF characteristics. Implemented interface boundaries exhibit the EDL thickness-boundary relation, i.e., the EDL thickness from MD simulations shows the tendency of converging toward the continuum approximation. Furthermore, inserting these values of EDL thicknesses into the continuum equation shows the convergence of flow transition of the molecular state to a neutral from an electrical violation phase, which takes a parabolic to plug-like shape in the velocity profile. Different interface boundaries also affect the hydrodynamic properties (viscosity and electroviscosity) of EOF, which varies from the bulk to interface region, as well as the fluid flow. Therefore, we can infer that, at the molecular level, the dissimilar atomic boundary and hydrodynamic properties dominate the electrokinetic flow. Our simulation results and theoretical model provide fundamental insightful information and guidelines for the EOF study based on the atomic interface and dynamic structure-based hydrodynamic property.