Use of Boundary-Driven Nonequilibrium Molecular Dynamics for Determining Transport Diffusivities of Multicomponent Mixtures in Nanoporous Materials

J Phys Chem B. 2022 Feb 10;126(5):1085-1100. doi: 10.1021/acs.jpcb.1c09159. Epub 2022 Feb 1.

Abstract

The boundary-driven molecular modeling strategy to evaluate mass transport coefficients of fluids in nanoconfined media is revisited and expanded to multicomponent mixtures. The method requires setting up a simulation with bulk fluid reservoirs upstream and downstream of a porous media. A fluid flow is induced by applying an external force at the periodic boundary between the upstream and downstream reservoirs. The relationship between the resulting flow and the density gradient of the adsorbed fluid at the entrance/exit of the porous media provides for a direct path for the calculation of the transport diffusivities. It is shown how the transport diffusivities found this way relate to the collective, Onsager, and self-diffusion coefficients, typically used in other contexts to describe fluid transport in porous media. Examples are provided by calculating the diffusion coefficients of a Lennard-Jones (LJ) fluid and mixtures of differently sized LJ particles in slit pores, a realistic model of methane in carbon-based slit pores, and binary mixtures of methane with hypothetical counterparts having different attractions to the solid. The method is seen to be robust and particularly suited for the study of study of transport of dense fluids and liquids in nanoconfined media.