Development and applications of molecular modeling techniques for the characterization of porous materials.

Identifier:  TDX:295

Handle: https://hdl.handle.net/20.500.11797/TDX295

Authors:  Figueroa Gerstenmaier, Susana

Abstract:
Porous materials are widely used in many branches of modern science and technology, such as catalysis, separation of mixtures, purification of fluids and fabrication of membranes. A successful application of porous solids requires a precise characterization of their surface and structural properties, as well as a good understanding of the physical and chemical behavior of fluids inside the pores. Some materials, such as zeolites, have well defined porous structures, but others, such as porous oxides, carbons and controlled-porous glasses, are quite amorphous. Therefore, a proper characterization of this kind of materials is an important topic, and more often than not, a complicated one. For many years, gas adsorption has been used to study properties of porous solids, since it is fast, simple and informative. Many methods were developed to extract information about porosity and surface properties of materials from adsorption isotherm data. In the last two decades, with the aid of the increasingly faster computers, the use of molecular modeling techniques has been gaining relevance. In this context, the general objective of this thesis is to develop tools at the molecular level using statistical mechanics for the characterization of adsorbent materials.<br/>After a brief introduction on the topic (chapter 1), chapter 2 is devoted to a review of the basic methodology employed in this work. In chapter 3 we have implemented the Fundamental-Measure density functional theory (FMT) due to Kierlik and Rosinberg to describe the adsorption of Lennard-Jones molecules in cylindrical pores. To our best knowledge, this is the first time that this theory is applied to a cylindrical geometry. The accuracy of the theory in predicting adsorption isotherms and density profiles is checked by comparison with Grand Canonical Monte Carlo simulations for a wide range of pore sizes, showing very good agreement in all cases. In addition, the theory has been applied to the adsorption in slit-like pores to study the influence of the pore geometry on this property. The results indicate that the confinement of the cylindrical geometry introduces significant differences in the shape of the adsorption isotherms and density profiles. These differences are relevant for the characterization of porous materials. Our results indicate that a layering behavior takes place in the smallest cylindrical pore considered, while the adsorption in a planar pore of the same size needs a much higher chemical potential to achieve a significant adsorption. As the pore size increases, the influence of the geometry becomes less important, although a certain shift in the capillary condensation transition can still be observed. Additionally, for wider pores, we obtain multilayer adsorption with capillary condensation at high chemical potentials, with the same qualitative behavior observed for both geometries. When the diameter size reaches the limit where the curvature effects are not of further relevance, the cylindrical pores reduce to the same quantitative behavior of the slit-like pores. The formation of a thin adsorbed layer at intermediate and large pore sizes seems to correspond to a second order thermodynamic phase transition, for the range of parameters used and the thermodynamic conditions studied. However, the results found seem to indicate some relationship between this behavior and the prewetting transition observed in semi-infinite geometries, especially in the vicinity of the critical end point of the prewetting line. The effect of the confinement is very important in this crossover behavior. From the comparison of Fundamental-Measure density functional theory calculations versus non-local density functional theory results, we conclude that the FMT is an excellent tool for the study of the behavior of fluids in confined <br/>cylindrical geometries.<br/>In chapter 4 we have applied the FMT in conjunction with a regularization method to estimate the pore-size distribution (PSD) of model porous glasses. We have chosen this particular material because it was developed with molecular modeling techniques, and a direct comparison can be made with the theory used here. An additional advantage of these model materials, versus experimental ones, is that in this case the size and shape of the pores is well known, as well as the position of the atoms in the surface, making it a perfect material to check the accuracy of the theoretical characterization methods available. Since there are several solutions of the adsorption integral equation compatible with the experimental adsorption isotherm, and several factors can hide defects of the molecular model, we have done the characterization in a systematic manner: we have first checked the accuracy of the FMT and the independent pore model for predicting the 'experimental' adsorption isotherms using the geometrical PSD already known for the materials. This has been done with individual cylindrical and slit-like pores. Secondly, once the adsorption isotherm was successfully reconstructed, we inverted the integral adsorption isotherm with a regularization procedure. The accuracy of the inversion method has also been checked before estimating the PSD of the different materials. Finally, once the method has been proved to be correct, we used it to estimate the PSD of four materials. We have also studied the influence of choosing different values of molecular parameters for the fluid-fluid and the solid-fluid interaction on the adsorption behavior of these systems. We have obtained that the independent pore model is adequate for the four materials investigated here. The slit-like geometry seems to represent the overall adsorption behavior better than the cylindrical geometry. As far as the PSD obtained with our procedure is concerned, the distributions obtained by inversion of the integral are in better agreement with the geometrical distributions than the ones calculated with the Barrett-Joyner-Halenda (BJH) method. The locus of the peak is at the same pore size, and all of them are unimodal, while the BJH distributions show a maximum systematically located at smaller pores, underestimating the PSD of the material, and they are not unimodal. Regarding the geometry of the individual pores that form the material, we can say that, although the PSD is broader than the geometrical ones, the adsorption predicted by a collection of individual slit-like pores is in almost quantitative agreement with the 'experimental' adsorption isotherm.<br/>Finally, in chapter 5 we have characterized three different samples of g­alumina, one of them without treatment and the others two calcined in a furnace during several hours at 823 and 1,023K. For this we have measured adsorption isotherms of nitrogen at 77.35K in a Micromeritics ASAP 2000 apparatus. Additionally, we have used the PSD's provided by the software of the experimental equipment using the BJH method. We have calculated theoretical isotherms by the FMT approach. We have inverted the adsorption integral equations with the regularization method and, finally, we have obtained the PSD's for our three samples of alumina, and the corresponding adsorption isotherms. In this way we have observed the influence of the calcination of alumina on its PSD. Moreover, we have tested the accuracy of the FMT/Regularization method in a systematic way. When we compared the PSD's obtained with the corresponding BJH distributions, we verified that in the two first cases (untreated alumina and alumina calcined at 823K) the BJH method underestimated the size of the pores, giving PSD's shifted to smaller sizes. In the case of alumina calcined at 1,023K, in which the sintering process has produced the disappearance of the smallest pores, favoring the wider ones, the BJH PSD's and the FMT/regularization PSD's perform very similar. With this, we corroborated the known fact that the BJH method is quite accurate in the macroporous region. Finally, we have predicted an adsorption isotherm of a different fluid (ethane) at a different temperature (333K) in one of the characterized materials (untreated alumina) with the aim of establishing the robustness of the PSD obtained. The agreement obtained shows that it is possible to use this characterization method and extrapolate the results at other conditions, provided that a enough number of different pore sizes are used to calculate the desirable isotherm, and the solid-fluid interaction parameters are well chosen.