Water adsorption in disordered mesoporous silica (Vycor) at 300 K and 650 K: a Grand Canonical Monte Carlo simulation study of hysteresis.

This numerical simulation paper focuses on the adsorption/desorption of water in disordered mesoporous silica glasses (Vycor-like). The numerical adsorbent was previously obtained by off lattice method, and was shown to reproduce quite well the micro- and mesotextural properties of real Vycor, as well as morphological (pore size distribution) and topological (pore interconnections) disorder. The water-water interactions are described by the SPC model while water-silica interactions are calculated in the framework of the PN-TrAZ model. The water adsorption/desorption isotherms and the configurational energies are calculated by the Grand Canonical Monte Carlo simulation method. The low pressure results compare well with experiments, showing the good transferability of the intermolecular potential. It is shown that if the hysteresis loop observed in the adsorption/desorption isotherm is considered as a true phase transition (which is actually still an open question in the case of disordered porous materials), then it is possible to calculate the grand potential by applying the thermodynamic integration scheme. The grand potential is shown to be multivalued for low (subcritical) temperature, and continuous for high (supercritical) temperature. A coexistence point is found within the hysteresis loop, actually close to the vertical desorption line. Below the equilibrium chemical potential, the gaslike branch is stable whereas the liquidlike branch is metastable. The situation is reversed above the coexistence point.