Molecular Simulation of Argon Adsorption and Diffusion in a Microporous Carbon with Poroelastic Couplings

Neglected for a long time in molecular simulations of fluid adsorption and transport in microporous carbons, adsorption-induced deformations of the matrix have recently been shown to have important effects on both sorption isotherms and diffusion coefficients. Here we investigate in detail the behavior of a recently proposed 3D-connected mature kerogen model, as a generic model of aromatic microporous carbon with atomic H/C ∼ 0.5, in both chemical and mechanical equilibrium with argon at 243 K over an extended pressure range. We show that under these conditions the material exhibits some viscoelasticity, and simulations of hundreds of nanoseconds are required to accurately determine the equilibrium volumes and sorption loadings. We also show that neglecting matrix internal deformations and swelling can lead to underestimations of the loading by up to 19% (swelling only) and 28% (swelling and internal deformations). The volume of the matrix is shown to increase up to about 8% at the largest pressure considered (210 MPa), which induces an increase of about 33% of both pore volume and specific surface area via the creation of additional pores, yet does not significantly change the normalized pore size distribution. Volume swelling is also rationalized by using a well-known linearized microporomechanical model. Finally, we show that self-diffusivity decreases with applied pressure, following an almost perfectly linear evolution with the free volume. Quantitatively, neglecting swelling and internal deformations tends to reduce the computed self-diffusivities.