A molecular dynamics simulation study of surface effects on gas permeation in free-standing polyimide membranes

Abstract

A 141100-atom model of a glassy ODPA–ODA polyimide free-standing membrane, corresponding to a thickness of two average radii of gyration for the 40-mers chains, has been studied using molecular dynamics (MD) simulations. Due to the large-scale of the fully atomistic model, a parallelized particle-mesh technique using an iterative solution of the Poisson equation had to be implemented for the efficient evaluation of the electrostatic interactions. With flattened-chain configurations, the density was found to adjust itself naturally in the middle of the membrane to ∼95% of the ODPA–ODA experimental value. At the free-standing surfaces, the density profile became sigmoïdal, indicating surface roughness. For comparison, two isotropic bulk models, one at the “normal” density as obtained for ODPA–ODA under ambient conditions and the other one at 95% of the normal-density, were built. Small gas probes were inserted into all three models in order to investigate whether the interfacial structure of the glassy free-standing membrane can influence penetrant transport. Gas diffusion in the low-density part of the interface was found to be very fast. The limiting value for the gas diffusion coefficient D membrane is only attained when the probes enter more dense regions in the membrane. Indeed, D membrane compares well with D bulk obtained for the 95%-density bulk system, i.e. about twice that in the normal-density bulk. Solubility is larger in the membrane than in both bulk models, thus suggesting an effect of chain flattening in addition to the density. Adsorption is particularly high at the free-standing interfaces.