Confined fluid dynamics in a viscoelastic, amorphous, and microporous medium: Study of a kerogen by molecular simulations and the generalized Langevin equation

We combine the use of molecular dynamics simulations and the generalized Langevin equation to study the diffusion of a fluid adsorbed within kerogen, the main organic phase of shales. As a class of microporous and amorphous materials that can exhibit significant adsorption-induced swelling, the dynamics of the kerogen's microstructure is expected to play an important role in the confined fluid dynamics. This role is investigated by conducting all-atom simulations with or without solid dynamics. Whenever the dynamics coupling between the fluid and solid is accounted for, we show that the fluid dynamics displays some qualitative differences compared to bulk fluids, which can be modulated by the amount of adsorbed fluid owing to adsorption-induced swelling. We highlight that working with the memory kernel, the central time correlation function of the generalized Langevin equation, allows the fingerprint of the dynamics of the solid to appear on that of the fluid. Interestingly, we observe that the memory kernels of fluid diffusion in kerogen qualitatively behave as those of tagged particles in supercooled liquids. We emphasize the importance of reproducing the velocity-force correlation function to validate the memory kernel numerically obtained as confinement enhances the numerical instabilities. This route is interesting as it opens the way for modeling the impact of fluid concentration on the diffusion coefficient in such ultra-confining cases.