MD simulation studies for effect of membrane structures and dynamics on gas permeation properties through microporous amorphous silica membranes

Imaginary amorphous silica membranes were prepared on a computer and gas permeation simulations were conducted using a dual control plane non-equilibrium molecular dynamics (DCP-NEMD) method. The Melt-Quench technique was employed to prepare various types of imaginary amorphous silica membranes which had different densities (from 1.3 to 2.2 g/cm3) and different mean pore sizes. Helium was adopted as a permeating gas species and its permeability was calculated at temperatures from 300 to 800 K. The Knudsen diffusion-like temperature dependencies of permeability could be observed for densities below 1.7 g/cm3, while the activated diffusion for the higher density models. We have also examined the effect of 3-body membrane potential parameters on membrane dynamics and gas permeation properties. The larger thermal vibration of oxygen atoms both in siloxane bonds and silanol groups on membranes could be observed for greater γ1 parameter is the SW potential function, which might result in the change of activation energy for gas permeation.