Molecular Simulation of Excess Isotherm and Excess Enthalpy Change in Gas-Phase Adsorption

Heat of adsorption, adsorption stress, and optimal storage of methane in slit and cylindrical carbon pores predicted by classical density functional theory

Temperature, pressure and pore-size dependences of the heat of adsorption, adsorption stress, and adsorption capacity of methane in simple models of slit and cylindrical carbon pores are studied using classical density functional theory (CDFT) and grand-canonical Monte-Carlo (MC) simulation. Studied properties depend nontrivially on the bulk pressure and the size of the pores. Heat of adsorption increases with loading, but only for sufficiently narrow pores. While the increase is advantageous for gas storage applications, it is less significant for cylindrical pores than for slits. Adsorption stress and the average adsorbed fluid density show oscillatory dependence on the pore size and increase with bulk pressure. Slit pores exhibit larger amplitude of oscillations of the normal adsorption stress with pore size increase than cylindrical pores. However, the increase of the magnitude of the adsorption stress with bulk pressure increase is more significant for cylindrical than for slit pores. Adsorption stress appears to be negative for a wide range of pore sizes and external conditions. The pore size dependence of the average delivered density of the gas is analyzed and the optimal pore sizes for storage applications are estimated. The optimal width of slit pore appears to be almost independent of storage pressure at room temperature and pressures above 10 bar. Similarly to the case of slit pores, the optimal radius of cylindrical pores does not exhibit much dependence on the storage pressure above 15 bar. Both optimal width and optimal radii of slit and cylindrical pores increase as the temperature decreases. A comparison of the results of CDFT theory and MC simulations reveals subtle but important differences in the underlying fluid models employed by the approaches. The differences in the high-pressure behaviour between the hard-sphere 2-Yukawa and Lennard-Jones models of methane, employed by the CDFT and MC approaches, respectively, result in an overestimation of the heat of adsorption by the CDFT theory at higher loadings. However, both adsorption stress and adsorption capacity appear to be much less sensitive to the differences between the models and demonstrate excellent agreement between the theory and the computer experiment.

Determination of Isosteric Heat of Adsorption by Quenched Solid Density Functional Theory

The heat of adsorption is one of the most important parameters characterizing energetic heterogeneity of the adsorbent surface. Heats of adsorption are either determined directly by calorimetry or calculated from adsorption isotherms measured at different temperatures using the thermodynamic Clausius-Clapeyron equation. Here, we present a method for calculating the isosteric heat of adsorption that requires as input only a single adsorption isotherm measured at one temperature. The proposed method is implemented with either nonlocal (NLDFT) or quenched solid (QSDFT) density functional theory models of adsorption that are currently widely used for calculating pore size distributions in various micro- and mesoporous solids. The pore size distribution determined from the same experimental isotherm is used for predicting the isosteric heat. The QSDFT method has advantages of taking into account two factors contributing to the structural heterogeneity of adsorbents: the molecular level roughness of the surface and the pore size distribution. The method is illustrated with examples of low temperature nitrogen and argon adsorption on selected samples of carbons of different degree of graphitization and MCM-41 mesoporous silicas of different pore size. The isosteric heat predictions from the NLDFT and QSDFT methods are compared against relevant experiments and the results of Monte Carlo (MC) simulations, with good agreement found in the cases where the surface model adequately reflects the pore surface roughness. Analyses with the QSDFT method show that the isosteric heat of adsorption significantly depends of the molecular level roughness of the adsorbent surface, which is ignored in NLDFT and MC models. The proposed QSDFT method with further verification can be used for calculating the isosteric heat as an additional parameter characterizing the adsorbent surface in parallel with routine calculations of the pore size distribution from a single adsorption isotherm.

Development of equations for differential and integral enthalpy change of adsorption for simulation studies

We present equations to calculate the differential and integral enthalpy changes of adsorption for their use in Monte Carlo simulation. Adsorption of a system of N molecules, subject to an external potential energy, is viewed as one of transferring these molecules from a reference gas phase (state 1) to the adsorption system (state 2) at the same temperature and equilibrium pressure (same chemical potential). The excess amount adsorbed is the difference between N and the hypothetical amount of gas occupying the accessible volume of the system at the same density as the reference gas. The enthalpy change is a state function, which is defined as the difference between the enthalpies of state 2 and state 1, and the isosteric heat is defined as the negative of the derivative of this enthalpy change with respect to the excess amount of adsorption. It is suitable to determine how the system behaves for a differential increment in the excess phase adsorbed under subcritical conditions. For supercritical conditions, use of the integral enthalpy of adsorption per particle is recommended since the isosteric heat becomes infinite at the maximum excess concentration. With these unambiguous definitions we derive equations which are applicable for a general case of adsorption and demonstrate how they can be used in a Monte Carlo simulation. We apply the new equations to argon adsorption at various temperatures on a graphite surface to illustrate the need to use the correct equation to describe isosteric heat of adsorption.

Net adsorption: a thermodynamic framework for supercritical gas adsorption and storage in porous solids

The thermodynamic treatment of adsorption phenomena is based on the Gibbs dividing surface, which is conceptually clear for a flat surface. On a flat surface, the primary extensive property is the area of the solid. As applications became more significant, necessitating microporous solids, early researchers such as McBain and Coolidge implemented the Gibbs definition by invoking a reference state for microporous solids. The mass of solid is used as a primary extensive property because surface area loses its physical meaning for microporous solids. A reference state is used to fix the hypothetical hyperdividing surface typically using helium as a probe molecule, resulting in the commonly used excess adsorption; experimentalists measure this reference state for each new sample. Molecular simulations, however, provide absolute adsorption. Theoreticians perform helium simulations to convert absolute to excess adsorption, mimicking experiments for comparison. This current structure of adsorption thermodynamics is rigorous (if the conditions for reference state helium measurements are completely disclosed) but laborious. In addition, many studies show that helium, or any other probe molecule for that matter, does adsorb, albeit to a small extent. We propose a novel thermodynamic framework, net adsorption, which completely circumvents the use of probe molecules to fix the reference state for each microporous sample. Using net adsorption, experimentalists calibrate their apparatus only once without any sample in the system. Theoreticians can directly calculate net adsorption; no additional simulations with a probe gas are necessary. Net adsorption also provides a direct indication of the density enhancement achieved (by using an adsorbent) over simple compression for gas (e.g., hydrogen) storage applications.

Computer simulation of argon adsorption on graphite surface from subcritical to supercritical conditions: the behavior of differential and integral molar enthalpies of adsorption

We investigate in detail the computer simulation of argon adsorption on a graphite surface over a very wide range of temperature, from below the triple point to well above the critical point. Adsorption over such a wide temperature range has not been reported previously in the form of adsorption isotherms and enthalpy change during adsorption. The adsorption isotherms can be classified broadly into four categories: below the triple point, the isotherms show stepwise character (a strict layering mechanism) with 2D condensation; type II (according to the IUPAC classification) is followed by isotherms at temperatures above the triple point and below the critical point and a sharp spike is seen for isotherms in the neighborhood of the critical point; and finally the typical behavior of a maximum is observed for isotherms above the critical point. For the isosteric heat, the heat curve (plotted against loading) remains finite for subcritical conditions but is infinite (singularity) at the maximum in excess loading for supercritical adsorption. For the latter case, a better representation of the energy change is the use of the integral molecular enthalpy because this does not exhibit a singularity as in the case of isosteric heat. We compare the differential and integral molecular enthalpies for the subcritical and supercritical adsorptions.