A grand canonical Monte-Carlo simulation study of xenon adsorption in a Vycor-like porous matrix

We have performed atomistic Grand Canonical Monte-Carlo (GCMC) simulations of adsorption of xenon in a Vycor-like matrix at 195 K. The disordered mesoporous network is obtained by applying a numerical 3D off-lattice reconstruction procedure to a simulation box originally containing silicon and oxygen atoms of a non-porous silica solid. In order to reduce the computational cost, we have applied a homothetic decrease of the simulation box dimensions which preserves the morphology and the topology of the pore network (the average pore dimension is then around 30 Angstrom). The surface chemistry is obtained in a realistic fashion by saturating all dangling bonds with hydrogen atoms. Small angle scattering spectra calculated on different numerical samples have evidenced a departure from Porod's law due to surface roughness. The simulated isotherms calculated on such disordered connected porous networks, show the capillary condensation phenomenon. The shape of the adsorption curves differs from that obtained for simple pore geometries. The analysis of the adsorbed quantity distribution indicates partial molecular-film formation depending on the local surface curvature and roughness.