Adsorption in aerogel materials has been studied for a number of years both by simulation and experiment. The extremely low density of these materials has made them particularly attractive for studies of the effects of weak confinement on fluid critical properties, but adsorption under conditions nearer to the fluid triple point has also been of considerable interest. Unlike porous materials such as vycor glasses and xerogels, low density aerogels have very low bulk moduli (generally in the range 1-100 bar). As a result they therefore display nontrivial mechanical deformations during adsorption and desorption processes, especially near to the capillary rise and desorption shoulder [1].
Due to the large length scales in aerogels, computational studies of adsorption have generally used lattice-based density functional theory approaches (a la Kierlik, et al. [2,3]) rather than molecular simulation. All such studies to date have kept the gel structure itself (generally obtained by DLCA or similar simulation [4]) frozen, which is a poor approximation for these materials; the purpose of this study is to address this approximation.
We have previously developed a simple coarse-grained model of silica sol-gel materials which was applied to investigate the formation and structural and mechanical properties of low density aerogels [5]. This model is fully flexible and can account for bond breaking and creation upon deformation of the gel structure. This model is structurally reasonable, and for appropriate choice of bond-formation control parameters, quantitative results for bulk modulus and the behavior under mechanical loading can be obtained.
In this work, the commonly used lattice-based fluid adsorption model is coupled to our coarse-grained off-lattice aerogel model using thermodynamic perturbation theory. This provides a complete free energy surface describing the behavior of the gel/fluid system at varying compression and fluid chemical potential. The resulting isotherms and dependence of gel compression on gas pressure are in excellent qualitative agreement with experimental findings, which are intuitively explained as motion on the underlying free energy surface. Perhaps most interestingly, in rigid aerogel models the adsorption isotherms are continuous in the capillary rise region, but when the gel is allowed to relax a sharp jump, corresponding to a first-order transition, is observed.
[1] Reichenauer and Scherer, J. Non-Cryst. Solids 277 (2000) 162.
[2] Kierlik et al., J. Phys.: Cond. Mat. 14 (2002) 9295
[3] Kierlik et al., PRL 87 (2001) 055701
[4] Salazar and Gelb, Mol. Phys. 102 (2004) 1015
[5] Gelb, J. Phys. Chem. C, 111 (2007) 15792