Atomistic molecular dynamics simulations of gas diffusion and solubility in rubbery amorphous hydrocarbon polymers

Diffusion coefficients D of H2, He, O2, N2, and CO2 in different rubbery amorphous polymeric matrices were estimated by atomistic molecular dynamics simulations at 298 K using the Einstein relationship, and compared with the relevant experimental values, where available. The simulated diffusion coefficients D of all the gases in all polymers considered almost regularly decreased with increasing molecular gas volumes and increasing polymer glass transition temperature. Further, solubility coefficients and heats of solution were obtained for all gases from Grand Canonical Monte Carlo simulations, which were also used to calculate sorption isotherms. In general, there is a good agreement between experimental and simulated values of diffusion and solubility coefficients for all gases considered.