Influence of system size in molecular dynamics simulations of gas permeation in glassy polymers

The dependence of gas permeation on system size in glassy polymers has been tested by creating several models containing different numbers of molecules of the same chain length for the ODPA-ODA homopolyimide and by subsequently determining the permeation characteristics of helium. Eight "standard size" 4150-atom, a 6225-atom, and a much more expensive computationally 56025-atom systems were generated using hybrid pivot Monte Carlo (PMC)-molecular dynamics (MD) single-chain sampling. Following the careful relaxation of the polyimide matrices, helium atoms were then inserted into these systems. and MD simulations were carried out at the same applied external conditions of constant temperature and pressure tensor. Average densities for the pure matrices all fell within 0.7% of the experimental value. Energetic and structural properties as well as solubilities and characterization of the void space were also found to be number independent, thus showing that the preparation procedure gives reproducible and reliable results. Helium diffusion was analyzed over periods up to 20 ns using different approaches, such as observation of the individual trajectories, mean-square displacements, distributions of penetrant displacements components, and van Hove correlation functions. No number dependence could be detected, whether the gas molecules were in the anomalous or in the Fickian regime.