Artifacts and misinterpretations in gas physisorption measurements and characterization of porous solids

This contribution provides a critical review of gas physisorption in the textural characterization of porous solids, with the focus on the artifacts in experimental data that lead to serious misinterpretation of the results derived from the analysis of adsorption isotherms. Apart from the problems related to the determination and interpretation of the BET area, we paid particular attention to the issues associated with the determination of pore size distribution; for example, the choice of the correct branch of the hysteresis loop and the network effects. Pitfalls in the analyses using either the classical macroscopic or the advanced microscopic (DFT, GCMC) methodology are addressed. The ultimate aim is to provide guidance for proper calculations and correct interpretation of physisorption data.

Effects of Molecular Cross-Sectional Areas of Adsorbed Nitrogen on the Brunauer-Emmett-Teller Analysis for Carbon-Based Slit Pores

Brunauer-Emmett-Teller (BET) surface area obtained by nitrogen adsorption is a commonly adopted value to characterize the specific surface area for porous materials. In the BET method, the widely applied cross-sectional area of nitrogen is 16.2 Ų, which has been found to be an oversimplified assumption. The adsorption isotherms of nitrogen simulated with the 2CLJ (Lennard-Jones) + 3q molecular model at 77 K was utilized to determine the cross-sectional area and its behavior as a function of pressure, pore size, and solid affinity. The cross-sectional area shows a negative relation with the pressure and varies with pore size and solid affinity. The overestimation of the BET surface area might arise from the inaccurately determined monolayer adsorption capacity and the cross-sectional area.

Monolayer Gas Adsorption on Graphene-Based Materials: Surface Density of Adsorption Sites and Adsorption Capacity

Surface density of adsorption sites on an adsorbent (including affinity-based sensors) is one of the basic input parameters in modeling of process kinetics in adsorption based devices. Yet, there is no simple expression suitable for fast calculations in current multiscale models. The published experimental data are often application-specific and related to the equilibrium surface density of adsorbate molecules. Based on the known density of adsorbed gas molecules and the surface coverage, both of these in equilibrium, we obtained an equation for the surface density of adsorption sites. We applied our analysis to the case of pristine graphene and thus estimated molecular dynamics of adsorption on it. The monolayer coverage was determined for various pressures and temperatures. The results are verified by comparison with literature data. The results may be applicable to modeling of the surface density of adsorption sites for gas adsorption on other homogeneous crystallographic surfaces. In addition to it, the obtained analytical expressions are suitable for training artificial neural networks determining the surface density of adsorption sites on a graphene surface based on the known binding energy, temperature, mass of adsorbate molecules and their affinity towards graphene. The latter is of interest for multiscale modelling.

The physisorption mechanism of SO2 on graphitized carbon

Sulfur dioxide (SO2) in flue gases emitted from fossil fuel power plants dramatically reduces the CO2 capture efficiency via adsorption, which is due to the potential reaction of SO2 with basic functional groups on the adsorbent. Physisorption rather than chemisorption is preferred, because adsorbents can be more easily regenerated by either reducing the pressure or increasing the temperature. Carbon is a suitable adsorbent for SO2 capture and widely used, and therefore it is important to study SO2 adsorption onto carbon with the Monte Carlo simulation to provide microscopic details to demarcate the roles of the basal plane of the graphene layer and the functional groups in adsorption. SO2 is a polar molecule like water, as they both carry partial charges, but they interact differently with functional groups. Instead of 3D-clusters in the case of water, SO2 is localized around the functional groups and spreads over the basal plane to form 2D-molecular layers because of the strong dispersive interactions with graphite. The results indicate that the functional group has a negligible effect on the enhancement of adsorption and its role is to localize 2D-clusters of SO2 molecules. For non-graphitized carbon, we have found that the greater loadings at low pressure compared to the highly graphitized carbon is due to the presence of defects (crevices) on the basal plane surface. Finally, to describe better the experimental data, we have found that the reduction in the interactions between adsorbed molecules in the first layer is because of the repulsion of their dipoles pointing normal to the surface, a phenomenon called surface mediation and is widely used in the description of gas adsorption on surfaces.

Monte Carlo optimization scheme to determine the physical properties of porous and nonporous solids

A new method, based on a Monte Carlo scheme, is developed to determine physical properties of nonporous and porous solids. In the case of nonporous solids, we calculate the surface area. This surface area is found as the sum of areas of patches of different surface energy on the solid, which is assumed to take a patchwise topology (i.e., adsorption sites of the same energy are grouped together in one patch). As a result of this assumption, we derive not only the surface area, but also the accessible volume and the surface energy distribution. In the case of porous solids, the optimization method is used to derive the surface area and the pore size distribution simultaneously. The derivation of these physical properties is based on adsorption data from a volumetric apparatus. We test this novel idea with the inversion problem of deriving surface areas of patches of different energies for a number of nonporous solids. The method is also tested with the derivation of the pore size distribution of some porous solid models. The results are very encouraging and demonstrate the great potential of this method as an alternative to the usual deterministic optimization algorithms which are known to be sensitive to the choice of the initial guess of the parameters. Since the geometrical parameters are physical quantities (i.e., only positive values are accepted), we also propose a scheme to enforce the positivity constraint of the solution.

New method to determine surface area and its energy distribution for nonporous solids: a computer simulation and experimental study

We present in this article a new method to determine the "geometrical" surface area of nonporous solids. This method is based on the total number of molecules dosed into the adsorption cell and a knowledge of the distribution of molecules between the gas phase and the surface phase. By matching this experimental amount with the corresponding theoretical equation, we can derive not only the surface area but also its energy distribution and the void volume of the adsorption cell. The method avoids the limitations of other methods presented in the literature. The BET method, for example, involves unrealistic assumptions and necessitates the choice of a molecular projection area. Our method does not suffer from these assumptions or limitations and is self-consistent, from the measurement of adsorption data to the final analysis of the surface area. The novelties of the method are the following: (i) it is valid over the complete range of reduced pressure, (ii) it does not require a molecular projection area, (iii) beside the total surface area, we also derive its energy distribution, and (iv) the helium expansion method (or any equivalent method) is not required to determine the void volume.

Ar, CCl(4) and C(6)H(6) adsorption outside and inside of the bundles of multi-walled carbon nanotubes-simulation study

This is the first paper reporting the results of systematic study of the adsorption of Ar, C(6)H(6) and CCl(4) on the bundles of closed and opened multi-walled carbon nanotubes. Using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations, we also study the effect of the introducing defects in the external and internal walls of osculating and separated nanotubes on Ar diffusion and on adsorption of all three adsorbates. The Ar diffusion coefficients obtained are very sensitive to the presence of defects. Simulated isotherms are discussed to show the relation between the shapes of the high resolution alpha(s)-plots and the mechanisms of adsorption. From obtained data, as well as from geometric considerations, from the VEGA ZZ package, and from simulations (ASA), the values of surface areas of all nanotubes are calculated and compared with those obtained using the most popular adsorption methods (BET, alpha(s) and the A,B,C-points). We show that the adsorption value for the C-point of the isotherm should be taken for the calculation of the specific surface area of carbon nanotubes to obtain a value which approaches the absolute geometric surface area. A fully packed monolayer is not created at the A-, B- or C-points of the isotherm; however, the number of molecules adsorbed at the latter point is closest to the number of molecules in the monolayer as calculated via the ASA method, the VEGA ZZ package or from geometric considerations.