A New Method for the Generation of Realistic Atomistic Models of Siliceous MCM-41

A new method is outlined for constructing realistic models of the mesoporous amorphous silica adsorbent, MCM-41. The procedure uses the melt-quench molecular dynamics technique. Previous methods are either computationally expensive or overly simplified, missing key details necessary for agreement with experimental data. Our approach enables a whole family of models spanning a range of pore widths and wall thicknesses to be efficiently developed and yet sophisticated enough to allow functionalisation of the surface – necessary for modelling systems such as self-assembled monolayers on mesoporous supports (SAMMS), used in nuclear effluent clean-up.

The models were validated in two ways. The first method involved the construction of adsorption isotherms from grand canonical Monte Carlo simulations, which were in line with experimental data. The second method involved computing isosteric heats at zero coverage and Henry law coefficients for small adsorbate molecules. The values obtained for carbon dioxide gave good agreement with experimental values.

We use the new method to explore the effect of increasing the preparation quench rate, pore diameter and wall thickness on low pressure adsorption. Our results show that tailoring a material to have a narrow pore diameter can enhance the physisorption of gas species to MCM-41 at low pressure.