Confined Water Determines Transport Properties of Guest Molecules in Narrow Pores

High-Pressure CO2 and CH4 Transport in Smooth Crystalline Silica Mesopores: A Molecular Dynamics Study

In this study, we use molecular dynamics (MD) simulation to study pressure-driven CO2 and CH4 flows and their slippage behaviors in β-cristobalite mesopores. The result illustrates that both CO2 and CH4 have an apparent adsorption layer on pore surface. However, significant differences in gas slippage are observed: CH4 flow shows considerable slippage, while it is negligible for CO2 flow. This disparity is attributed to the collective effect of gas molecular configurations and surface structure. The linear molecular structure of CO2 allows it to align perpendicular to the surface, even penetrating into the surface. Notably, the perpendicular orientation of CO2 molecules is energetically favored near the center of the equilateral triangle formed by adjacent oxygen atoms on β-cristobalite surface. Conversely, the symmetric molecular structure of CH4, coupled with its larger size, prevents its penetration into pore surfaces. Therefore, despite smooth crystalline surfaces, CO2 topological accessible plane is much more curved than that of CH4. Consequently, CO2 displays hesitating motions undergoing rotational movements, which significantly hinders its slippage. This study highlights the collective influences of gas molecular characteristics and surface structure on gas slippage, affording important insights into gas sequestration and the development of functional materials for gas separation.

Molecular Insight into CO2 Improving Oil Mobility in Shale Inorganic Nanopores Containing Water Films

CO2 injection into shale reservoirs has been recognized as one of the most promising techniques for enhanced oil recovery and carbon capture, utilization, and storage. However, the omnipresent nanopores and the water films formed near the pore walls affect the understanding of mechanisms of CO2 regulating crude oil mobility in shale nanopores. In this work, we employ molecular dynamics simulations to study the occurrence and flow of CO2 and octane (nC8) mixtures in quartz nanopores containing water films, to illustrate the impact mechanisms of CO2 on nC8 mobility. The results indicate that nC8 exists between water films, and CO2 is mainly miscible with nC8 in the pore center, and a small portion of it accumulates at the interface between nC8 and the water film. CO2 significantly decreases the apparent viscosity of nC8 in both the bulk nC8 region and the nC8-water interface region, improving nC8 fluidity. As the percentage of CO2 in the CO2 and nC8 mixtures increases from 0 to 50%, the mean flow velocities of nC8 in the bulk phase region and the nC8-water interface region increase by 92.85 and 60.64%, respectively. Three major microscopic mechanisms of CO2 improving nC8 fluidity in quartz nanopores with water films are summarized: (i) CO2 reduces friction between nC8 and the water film by increasing the angle between nC8 molecules and the plane of the water film; (ii) CO2 widens the distance between nC8 molecules, causing the volume expansion of nC8 and its viscosity reduction; (iii) CO2 significantly increases the most probable and average velocities of nC8 molecules, thus improving their mobility. Our results enhance the comprehension of how CO2 facilitates oil flow in water-bearing shale reservoirs and the exploitation of unconventional oil resources.

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. In this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

Control of the Structural Charge Distribution and Hydration State upon Intercalation of CO2 into Expansive Clay Interlayers

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO₂ preferentially into their interlayers and therefore potentially act as a material for CO₂ separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO₂ into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO₂ intercalation. Our free energy calculations indicate that CO₂ favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO₂/H₂O ratio in clay interlayers as compared to that in bulk water. CO₂ molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO₂ intercalation. The structural charge distribution significantly affects the amount of CO₂ intercalated in the interlayers. These results provide a mechanistic understanding of CO₂ intercalation in clays for CO₂ separation, capture, and storage.

Breakthrough pressure of oil displacement by water through the ultra-narrow kerogen pore throat from the Young-Laplace equation and molecular dynamic simulations

As one common unconventional reservoir, shale plays a pivotal role to compensate the depletion of conventional oil resources. There are numerous nanoscale pores and ultra-narrow pore throats (sub 2 nm) in shale media. To displace oil through ultra-narrow pore throats by water, one needs to overcome excessively-high capillary pressure. Understanding the water-oil two-phase displacement process through pore throats is critical to numerical simulation on tight/shale oil exploitation and ultimate oil recovery estimation. In this work, we use molecular dynamics simulations to investigate oil (represented by n-octane) displacement by water through a ~2 nm kerogen (represented by Type II-C kerogen) pore throat. Besides, the applicability of the Young-Laplace equation to the ultra-narrow kerogen pore throat has been assessed. We find that although the Type II-C kerogen is generally oil-wet, water has an excellent displacement efficiency without the oil film on the substrate, thanks to the hydrogen bonding formed between water and heteroatoms (such as O, N, and S) on the kerogen surface. Unlike previous studies, the capillary pressure obtained from the widely used Young-Laplace equation shows a good agreement with the breakthrough pressure obtained from MD simulations for the ∼2 nm kerogen pore throat. Our work indicates that explicitly considering intermolecular interactions as well as atomistic and molecular level characteristics is imperative to study the two-phase displacement process through ultra-narrow pore throats.

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (∼20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. The dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule's local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

Homogenization of Continuum-Scale Transport Properties from Molecular Dynamics Simulations: An Application to Aqueous-Phase Methane Diffusion in Silicate Channels

Silica-based materials including zeolites are commonly used for wide-ranging applications including separations and catalysis. Substrate transport rates in these materials can significantly influence the efficiency of such applications. Two factors that contribute to transport rates include (1) the porosity of the silicate matrix and (2) nonbonding interactions between the diffusing species and the silicate surface. These contributions generally emerge from disparate length scales, namely, "microscopic" (roughly nanometer-scale) and "macroscopic" (roughly micron-scale), respectively. Here, we develop a simulation framework to estimate the simultaneous impact of these factors on methane mass transport in silicate channels. Specifically, we develop a model of methane transport using homogenization theory to obtain transport parameters valid at length scales of hundreds to thousands of nanometers. These parameters implicitly reflect interactions taking place at fractions of a nanometer. The inputs to the homogenization analysis are data from extensive molecular dynamics simulations that incorporate atomistic-scale interactions, processed to yield local diffusion coefficients and mean force potentials. With this model, we demonstrate how nuances in silicate hydration and silica/methane interactions impact methane diffusion rates in silicate materials, including the effects of silicate surface chemistry such as the presence of silanol groups. The molecular dynamics simulations indicate that methane diffusivity at the silica surface is lower than the bulk-like rates observed at the center of channels of sufficient width. However, potentials of mean force generally evidence attractive methane/silica interactions that enhance diffusion overall. By simultaneously accounting for both of these effects, we show that the effective diffusion coefficient for the nanoporous silicate can be approximately double the value of estimates assuming fully bulk-like behavior in the channel. This study therefore demonstrates the importance of determining diffusion coefficients and potentials of mean force at an atomistic level when estimating transport properties in bulk materials. Importantly, we provide a simple homogenization framework to incorporate these molecular-scale attributes into bulk material transport estimates. This hybrid homogenization/molecular dynamics approach will be of general use for describing small-molecule transport in materials with detailed molecular interactions.

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

Osmotic transport in nanoconfined aqueous electrolytes provides alternative venues for water desalination and "blue energy" harvesting. The osmotic response of nanofluidic systems is controlled by the interfacial structure of water and electrolyte solutions in the so-called electrical double layer (EDL), but a molecular-level picture of the EDL is to a large extent still lacking. Particularly, the role of the electronic structure has not been considered in the description of electrolyte/surface interactions. Here, we report enhanced sampling simulations based on ab initio molecular dynamics, aiming at unravelling the free energy of prototypical ions adsorbed at the aqueous graphene and hBN interfaces, and its consequences on nanofluidic osmotic transport. Specifically, we predicted the zeta potential, the diffusio-osmotic mobility, and the diffusio-osmotic conductivity for a wide range of salt concentrations from the ab initio water and ion spatial distributions through an analytical framework based on Stokes equation and a modified Poisson-Boltzmann equation. We observed concentration-dependent scaling laws, together with dramatic differences in osmotic transport between the two interfaces, including diffusio-osmotic flow and current reversal on hBN but not on graphene. We could rationalize the results for the three osmotic responses with a simple model based on characteristic length scales for ion and water adsorption at the surface, which are quite different on graphene and on hBN. Our work provides fundamental insights into the structure and osmotic transport of aqueous electrolytes on 2D materials and explores alternative pathways for efficient water desalination and osmotic energy conversion.

Application of the dynamic mean field theory to fluid transport in slit pores

We explore the applicability of the lattice model and dynamic mean field theory as a computationally efficient tool to study transport across heterogeneous porous media, such as mixed matrix membranes. As a starting point and to establish some basic definitions of properties analogous to those in the off-lattice systems, we consider transport across simple models of porous materials represented by a slit pore in a chemical potential gradient. Using this simple model, we investigate the distribution of density and flux under steady state conditions, define the permeability across the system, and explore how this property depends on the length of the pore and the solid-fluid interactions. Among other effects, we observe that the flux in the system goes through a maximum as the solid-fluid interaction is varied from weak to strong. This effect is dominated by the behavior of the fluid near the walls and is also confirmed by off-lattice molecular dynamics simulations. We further extend this study to explore transport across heterogeneous slit pore channels composed of two solids with different values of solid-fluid interaction strengths. We demonstrate that the lattice models and dynamic mean field theory provide a useful framework to pose questions on the accuracy and applicability of the classical theories of transport across heterogeneous porous systems.

Aqueous films on pore surfaces mediate adsorption and transport of gases through crowded nanopores

Interactions of trapped reservoir gases within organic-rich and brine-bearing sedimentary rocks have direct relevance to many geoenergy applications. Extracting generalizable information from experimental campaigns is hindered by the fact that geological systems are extremely complex. However, modern computational tools offer the opportunity of studying systems with controlled complexity, in an effort to better understand the mechanisms at play. Employing molecular dynamics, we examine here adsorption and transport of gases containing CH₄ and either CO₂ or H₂S within amorphous silica nanopores filled with benzene. We explicitly quantify the effect of small amounts of water/brines at geological temperature and pressure conditions. Because of wetting, the presence of brines lessens the adsorption capacity of the aromatic-filled pore. The simulation results show salt-specific effects on the transport properties of the gases when either KCl or CaCl₂ brines are considered, although adsorption was not affected. The acid gases considered either facilitate or hinder CH₄ transport depending on whether they are more or less preferentially adsorbed within the pore as compared to benzene, and this effect is mediated by the presence of water/brines. Our simulation results could be used to extract thermodynamic quantities that in the future will help to optimize transport of various gases through organic-rich and brine-bearing sedimentary rocks, which is likely to have a positive impact on both hydrocarbon production and carbon sequestration applications. As a first step, a phenomenological model is presented here, which allows one to predict permeability based on interatomic energies.

A Novel Modeling Approach to Stochastically Evaluate the Impact of Pore Network Geometry, Chemistry and Topology on Fluid Transport

Fine-grained sandstones, siltstones, and shales have become increasingly important to satisfy the ever-growing global energy demands. Of particular current interest are shale rocks, which are mudstones made up of organic and inorganic constituents of varying pore sizes. These materials exhibit high heterogeneity, low porosity, varying chemical composition and low pore connectivity. Due to the complexity and the importance of such materials, many experimental, theoretical and computational efforts have attempted to quantify the impact of rock features on fluids diffusivity and ultimately on permeability. In this study, we introduce a stochastic kinetic Monte Carlo approach developed to simulate fluid transport. The features of this approach allow us to discuss the applicability of 2D vs 3D models for the calculation of transport properties. It is found that a successful model should consider realistic 3D pore networks consisting of pore bodies that communicate via pore throats, which however requires a prohibitive amount of computational resources. To overcome current limitations, we present a rigorous protocol to stochastically generate synthetic 3D pore networks in which pore features can be isolated and varied systematically and individually. These synthetic networks do not correspond to real sample scenarios but are crucial to achieve a systematic evaluation of the pore features on the transport properties. Using this protocol, we quantify the contribution of the pore network’s connectivity, porosity, mineralogy, and pore throat width distribution on the diffusivity of supercritical methane. A sensitivity analysis is conducted to rank the significance of the various network features on methane diffusivity. Connectivity is found to be the most important descriptor, followed by pore throat width distribution and porosity. Based on such insights, recommendations are provided on possible technological approaches to enhance fluid transport through shale rocks and equally complex pore networks. The purpose of this work is to identify the significance of various pore network characteristics using a stochastic KMC algorithm to simulate the transport of fluids. Our findings could be relevant for applications that make use of porous media, ranging from catalysis to radioactive waste management, and from environmental remediation to shale gas production.

Evidence of Facilitated Transport in Crowded Nanopores

Fluid transport in nature often occurs through crowded nanopores, where a number of phenomena can affect it, because of fluid-fluid and fluid-solid interactions, as well as the presence of organic compounds filling the pores and their structural fluctuations. Employing molecular dynamics, we probe here the transport of fluid mixtures (CO₂-CH₄ and H₂S-CH₄) through silica nanopores filled with benzene. Both CO₂ and H₂S are strongly adsorbed within the organic-filled pore, partially displacing benzene. Unexpectedly, CO₂/H₂S adsorption facilitates CH₄ transport. Analysis of the trajectories suggests that both CO₂ and H₂S act as vehicle-like carriers and might swell benzene, generating preferential transport pathways within the crowded pore. The results are useful for identifying unexpected transport mechanisms and for developing engineering approaches that could lead to storage of CO₂ in caprocks.

Structure and dynamics of ethane confined in silica nanopores in the presence of CO2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C2H6) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C2H6 and CO2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C2H6, C2H6:CO2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C2H6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C2H6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO2 molecules are more strongly attracted to the pore surface compared to C2H6. The C2H6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO2 molecules play in enhancing C2H6 diffusion ("molecular lubrication") at high CO2:C2H6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels

We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.

Effects of water on the stochastic motions of propane confined in MCM-41-S pores

Quasi-elastic neutron scattering (QENS) and molecular dynamics simulations (MDS) reveal the effects of water on the structure and dynamics of propane confined in 1.5 nm wide cylindrical pores of MCM-41-S.

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling

CRC-MD simulations show that nanopores in shales bounded by clay minerals have a strong preference for CO₂ relative to CH₄.

Nanomodel visualization of fluid injections in tight formations

The transport and phase change of a complex fluid mixture under nanoconfinement is of fundamental importance in nanoscience, and limits the recovery efficiency from tight oil reservoirs (<10%). Herein, through experiments and supporting theory we characterize the transport and phase change of a nanoconfined complex fluid mixture. Our nanofluidic platform, nanomodel, replicates shale reservoirs in terms of mean pore size (∼100 nm), permeability (∼μD) and porosity (∼10%). We screen conditions for the most promising shale EOR strategies, directly quantifying their pore-scale efficiency and underlying mechanisms. We find that immiscible gas (N2) flooding presents a prohibitively large capillary pressure threshold (∼2 MPa). Miscible (CO2) gas flooding eliminates this threshold leading to film-wise stable oil displacement with high recovery efficiency. Strong capillary forces present barriers as well as opportunities for recovery strategies unique to nanoporous reservoirs by transitioning from a miscible to an immiscible condition locally within the reservoir. These results quantify the fundamental transport and phase change mechanisms applicable to nanoconfined complex fluids, with direct implications in unconventional oil as well as nanoporous media more broadly.

Confinement Effect on Molecular Conformation of Alkanes in Water-Filled Cavitands: A Combined Quantum/Classical Density Functional Theory Study

The depletion force exerted on an alkane molecule from surrounding solvent may greatly alter its conformation. Such a behavior is closely related to the selective molecular recognition, molecular sensors, self-assembly, and so on. Herein, we report a multiscale theoretical study on the conformational change of a single alkane molecule confined in water-filled cavitands, in which the quantum and classical density functional theories (DFTs) are combined to determine the grand potential of alkane-water system. Specifically, the intrinsic free energy of the alkane molecule is tackled by quantum DFT, while the solvent effect arising from the solvent density inhomogeneity in confined space is addressed by classical DFT. By varying the alkane chain length, pore size, and wettability of inner pore surface, we find that pore confinement and hydrophilic inner surface facilitate the alkane conformational change from extended state to helical state, which becomes more significant as the alkane chain length increases. Our findings, which are in line with previous experimental observations, provide not only the microscopic mechanism but also theoretical guidance for elaborately manipulating molecular conformation at the nanoscale.

Nanoscale Phase Measurement for the Shale Challenge: Multicomponent Fluids in Multiscale Volumes

Hydrocarbon recovery from shale reservoirs provides an increasing share of world energy. These resources are multicomponent fluid mixtures within multiscale geometries, and understanding their associated phase-change thermodynamics presents an array of challenges for experimentalists, theorists, operators, and policy makers. Here, we quantify hydrocarbon mixture phase behavior via direct imaging of connected channels spanning 4 orders of magnitude (10 nm to 10 μm) with supporting density functional theory. The methane/propane mixture dew point shifts, with early condensation of heavy components in nanopores because of a combination of capillarity and competitive surface adsorption. The bubble point in nanoconfinement is found to be deeply suppressed (∼3-fold), to below the bulk dew point of the original mixture, because of the exchange of mixture components with larger connected volumes. The trapping of the heaviest components of hydrocarbon mixtures within the smallest connected pores has implications for shale operations, reserve estimation, and ultimately energy security.

Solvent Transport Behavior of Shear Aligned Graphene Oxide Membranes and Implications in Organic Solvent Nanofiltration

Solvent transport in membranes composed of stacked sheets of graphene oxide (GO) with molecular scale channels and a complex arrangement of hydrophobic and hydrophilic domains is not well understood. Here, we observe that the interlayer space between GO sheets expands in different solvents without disturbing the membrane integrity and is typically larger in aqueous media compared to nonaqueous media. However, the membranes have a tighter molecule sieving feature in aqueous media as demonstrated by lower permeance and higher solute rejection arising from interfacial water layers "sticking" to charged polar groups. As a result of this polar interaction, the permeance of polar solvents in GO membrane scales inversely to the polarity of the solvent, which is contrary to other polymeric and ceramic hydrophilic membranes and also scales inversely to the viscosity of solvents as per continuum expectations. We highlight the extended solvent-handling space of GO membranes, such as in polar protic, polar aprotic, and nonpolar solvents, demonstrating versatility over a commercial nanofiltration membrane, and we predict exciting new applications in advanced separation engineering.

Density Fluctuation in Aqueous Solutions and Molecular Origin of Salting-Out Effect for CO2

Using molecular dynamics simulation, we studied the density fluctuations and cavity formation probabilities in aqueous solutions and their effect on the hydration of CO₂. With increasing salt concentration, we report an increased probability of observing a larger than the average number of species in the probe volume. Our energetic analyses indicate that the van der Waals and electrostatic interactions between CO₂ and aqueous solutions become more favorable with increasing salt concentration, favoring the solubility of CO₂ (salting in). However, due to the decreasing number of cavities forming when salt concentration is increased, the solubility of CO₂ decreases. The formation of cavities was found to be the primary control on the dissolution of gas, and is responsible for the observed CO₂ salting-out effect. Our results provide the fundamental understanding of the density fluctuation in aqueous solutions and the molecular origin of the salting-out effect for real gas.

A kinetic Monte Carlo approach to study fluid transport in pore networks

The mechanism of fluid migration in porous networks continues to attract great interest. Darcy's law (phenomenological continuum theory), which is often used to describe macroscopically fluid flow through a porous material, is thought to fail in nano-channels. Transport through heterogeneous and anisotropic systems, characterized by a broad distribution of pores, occurs via a contribution of different transport mechanisms, all of which need to be accounted for. The situation is likely more complicated when immiscible fluid mixtures are present. To generalize the study of fluid transport through a porous network, we developed a stochastic kinetic Monte Carlo (KMC) model. In our lattice model, the pore network is represented as a set of connected finite volumes (voxels), and transport is simulated as a random walk of molecules, which "hop" from voxel to voxel. We simulated fluid transport along an effectively 1D pore and we compared the results to those expected by solving analytically the diffusion equation. The KMC model was then implemented to quantify the transport of methane through hydrated micropores, in which case atomistic molecular dynamic simulation results were reproduced. The model was then used to study flow through pore networks, where it was able to quantify the effect of the pore length and the effect of the network's connectivity. The results are consistent with experiments but also provide additional physical insights. Extension of the model will be useful to better understand fluid transport in shale rocks.

Facile Anhydrous Proton Transport on Hydroxyl Functionalized Graphane

Graphane functionalized with hydroxyl groups is shown to rapidly conduct protons under anhydrous conditions through a contiguous network of hydrogen bonds. Density functional theory calculations predict remarkably low barriers to diffusion of protons along a 1D chain of surface hydroxyls. Diffusion is controlled by the local rotation of hydroxyl groups, a mechanism that is very different from that found in 1D water wires in confined nanopores or in bulk water. The proton mean square displacement in the 1D chain was observed to follow Fickian diffusion rather than the expected single-file mobility. A charge analysis reveals that the charge on the proton is essentially equally shared by all hydrogens bound to oxygens, effectively delocalizing the proton.

Molecular dynamics simulations of propane in slit shaped silica nano-pores: direct comparison with quasielastic neutron scattering experiments

MD simulations reveal the origin of anomalous pressure dependence of propane diffusion in silica mesopores.

Impact of solvent granularity and layering on tracer hydrodynamics in confinement

Classic hydrodynamic arguments establish that when a spherical tracer particle is suspended between parallel walls, tracer-wall coupling mediated by the solvent will cause the tracer to exhibit position-dependent diffusivity. We investigate how the diffusivity profiles of confined tracers are impacted by the diameter size-ratio of the tracer to solvent: starting from the classic limit of infinite size-ratio (i.e., continuum solvent), we consider size-ratios of four or less to examine how hydrodynamic predictions are disrupted for systems where the tracer and solvent are of similar scale. We use computer simulations and techniques based on the Fokker-Planck formalism to calculate the diffusivity profiles of hard-sphere tracer particles in hard-sphere solvents, focusing on the dynamics perpendicular to the walls. Given wall separations of several tracer diameters, we first consider confinement between hard walls, where anisotropic structuring at the solvent lengthscale generates inhomogeneity in the tracer free-energy landscape and undermines hydrodynamic predictions locally. We then introduce confining planes that we term transparent walls, which restrict tracer and solvent center-accessibilities while completely eliminating static anisotropy, and reveal position-dependent signatures in tracer diffusivity solely attributable to confinement. With or without suppressing static heterogeneity, we find that tracer diffusivity increasingly deviates on a local basis from hydrodynamic predictions at smaller size-ratios. However, hydrodynamic theory still approximately captures spatially-averaged dynamics across the pores even for very small tracer-solvent size-ratios over a wide range of solvent densities and wall separations.