Coarse-Grained Molecular Simulation of Penetrant Diffusion in a Glassy Polymer Using Reverse and Kinetic Monte Carlo

A coarse-grained method for simulating diffusion of a small molecule within a glassy polymer was developed and implemented. The method builds on our previous work in which molecular-level jump rates between likely sorption states were calculated with multidimensional transition-state theory incorporating explicit chain motions that accompany each jump. In this work we first use a reverse Monte Carlo approach to generate large microstructures of sorption states and jump paths whose size, connectivity, and rate constant distributions match those found in detailed molecular simulations of methane in glassy atactic polypropylene. Next we simulate diffusion of isolated penetrant molecules in these microstructures using kinetic Monte Carlo. Over small to moderate times, mean-squared displacement increases sublinearly (anomalous diffusion) in structures of either low or moderate connectivity and with either uniform rate constants or a distribution of rate constants. At long times, regular diffusion is observed i...