Modeling Self-Diffusion of Simple Fluids in Nanopores

Position-Dependent Dynamics Explain Pore-Averaged Diffusion in Strongly Attractive Adsorptive Systems

Using molecular simulations, we investigate the relationship between the pore-averaged and position-dependent self-diffusivity of a fluid adsorbed in a strongly attractive pore as a function of loading. Previous work (Krekelberg, W. P.; Siderius, D. W.; Shen, V. K.; Truskett, T. M.; Errington, J. R. Connection between thermodynamics and dynamics of simple fluids in highly attractive pores. Langmuir 2013, 29, 14527-14535, doi: 10.1021/la4037327) established that pore-averaged self-diffusivity in the multilayer adsorption regime, where the fluid exhibits a dense film at the pore surface and a lower density interior pore region, is nearly constant as a function of loading. Here we show that this puzzling behavior can be understood in terms of how loading affects the fraction of particles that reside in the film and interior pore regions as well as their distinct dynamics. Specifically, the insensitivity of pore-averaged diffusivity to loading arises from the approximate cancellation of two factors: an increase in the fraction of particles in the higher diffusivity interior pore region with loading and a corresponding decrease in the particle diffusivity in that region. We also find that the position-dependent self-diffusivities scale with the position-dependent density. We present a model for predicting the pore-average self-diffusivity based on the position-dependent self-diffusivity, which captures the unusual characteristics of pore-averaged self-diffusivity in strongly attractive pores over several orders of magnitude.

Transport Diffusion of Light Gases in Polyethylene Using Atomistic Simulations

We explore the temperature dependence of the self-, corrected-, and transport-diffusivities of CO₂, CH₄, and N₂ in a polyethylene (PE) polymer membrane through equilibrium molecular dynamics simulations. We also investigate the morphology of the polymer membrane based on the intermolecular radial distribution function, free volume, and pore size distribution analysis. The results indicate the existence of 1.5-3 Å diameter pores in the PE membrane, and with the increase in the temperature, the polymer swells linearly with changing slope at 450 K in the absence of gas and exponentially in the presence of gas. The gas adsorption isotherms extracted via a two-step methodology, considering the dynamics and structural transitions in the polymer matrix upon gas adsorption, were fitted using a "two-mode sorption" model. Our results suggest that CO₂ adsorbs strongly, whereas N₂ shows weak adsorption in PE. The results demonstrate that CO₂ is more soluble, whereas N₂ is least soluble. Further, it is found that an increase in the temperature negatively impacts the solubility of CO₂ and CH₄ but positively for N₂; this reverse solubility behavior is due to increased availability of pores accessible to N₂, which are kinetically closed at the lowest temperatures. The reported self-diffusivities of the gases from our simulations are on the order of 10^-6 cm²/s, consistent with the experimental evidence, whereas transport-diffusivities are 2 orders of magnitude higher than self-diffusivities. Furthermore, the temperature dependence of the self-diffusivity follows Arrhenius behavior, whereas the transport-diffusivity follows non-Arrhenius behavior having different activation energies in low and high temperature regions. Also, it is seen that loading has little effect on the self- and corrected-diffusion coefficients of all gases in the PE membrane.

Influence of Structural Heterogeneity on Diffusion of CH4 and CO2 in Silicon Carbide-Derived Nanoporous Carbon

We investigate the influence of structural heterogeneity on the transport properties of simple gases in a Hybrid Reverse Monte Carlo (HRMC) constructed model of silicon carbide-derived carbon (SiC-DC). The energy landscape of the system is determined based on free energy analysis of the atomistic model. The overall energy barriers of the system for different gases are computed along with important properties, such as Henry constant and differential enthalpy of adsorption at infinite dilution, and indicate hydrophobicity of the SiC-DC structure and its affinity for CO₂ and CH₄ adsorption. We also study the effect of molecular geometry, pore structure and energy heterogeneity considering different hopping scenarios for diffusion of CO₂ and CH₄ through ultramicropores using the Nudged Elastic Band (NEB) method. It is shown that the energy barrier of a hopping molecule is very sensitive to the shape of the pore entry. We provide evidence for the influence of structural heterogeneity on self-diffusivity of methane and carbon dioxide using molecular dynamics simulation, based on a maximum in the variation of self-diffusivity with loading. A comparison of the MD simulation results with self-diffusivities from quasi-elastic neutron scattering (QENS) measurements and, with macroscopic uptake-based low-density transport coefficients, reveals the existence of internal barriers not captured in MD simulation and QENS experiments. Nevertheless, the simulation and macroscopic uptake-based diffusion coefficients agree within a factor of 2-3, indicating that our HRMC model structure captures most of the important energy barriers affecting the transport of CH₄ in the nanostructure of SiC-DC.

Friction between solids and adsorbed fluids is spatially distributed at the nanoscale

The widespread developments in the use of nanomaterials in catalysis, adsorption, and nanofluidics present significant new challenges in achieving optimal adsorbed fluid flow characteristics. Here we demonstrate, using molecular dynamics simulations of nanoconfined fluids, that at nanoscales, fluid-solid friction is not restricted to a sharp interface as is commonly assumed; instead it is distributed over the whole adsorbed fluid phase, and is strongest in an interfacial region that is not negligible in comparison to the system size. Our simulations yield position-dependent dynamical fluid-solid friction coefficients, and lead to a modification of conventional hydrodynamics, incorporating distributed momentum loss in the fluid due to fluid-solid interaction. The results demonstrate that the usual concepts of slip length or interfacial friction coefficient are meaningful only for uniform fluids, and lose their significance for adsorbates in nanospaces, which are intrinsically inhomogeneous. We show that static friction coefficients, based on equilibrium density distributions, follow the same spatial dependence as the dynamical coefficients. These results open up possibilities for tailoring nanomaterials and surfaces to engineer low friction pathways for adsorbed fluid flow by tuning the potential energy landscape.

Connection between thermodynamics and dynamics of simple fluids in highly attractive pores

Using molecular simulations, we investigate the structural and diffusive dynamics properties of a model fluid in highly absorptive cylindrical pores. At subcritical temperatures, self-diffusivity displays three distinct regimes as a function of average pore density ρ: (1) a decrease in self-diffusivity with increasing ρ at low ρ, (2) constant self-diffusivity with respect to varying ρ at moderate density, and (3) a decrease in self-diffusivity with increasing ρ at high density. These regimes are closely linked to the thermodynamic properties of the fluid in the pore, specifically, the adsorption isotherm, isosteric heat of adsorption, and the density profile. We show that these three diffusivity regimes qualitatively correspond to three distinct adsorption regimes: monolayer formation, multilayer adsorption, and pore filling, respectively. In addition, we find that the self-diffusivity is a universal function of the local film density in the monolayer formation regime at subcritical temperatures. The results of this work suggest a potential means to estimate the self-diffusivity over a broad pressure range using a limited number of experiments.