The special features of equilibrium adsorption of argon on homogeneous and inhomogeneous surfaces

Comparative patterns of equilibrium adsorption of argon on the surface of graphitized thermal carbon black (GCB) and the inhomogeneous surfaces of nongraphitized carbon black and silica at 77 and 87.3 K were considered. It was shown that argon acquires the properties of a special phase with a layered structure and exhibits two-dimensional phase transitions with the formation of crystal-like layers near the homogeneous surface of GCB even at a temperature exceeding the triple point. However, already at a distance of three-four molecular diameters from the surface, adsorbed argon behaves as a bulk phase in a weak external field. The defect surface of nongraphitized carbon black and the amorphous surface structure of silica destroy the longrange order of adsorbed argon and lower its solidification temperature. Therefore, argon adsorbed at a temperature of 77 K, i.e., below the triple point, exhibits the properties of a supercooled liquid. The applicability of density functional theory to describe argon isotherms and heat of adsorption on inhomogeneous surfaces was demonstrated.     

Modeling of adsorption on nongraphitized carbon surface: GCMC simulation studies and comparison with experimental data

We model nongraphitized carbon black surfaces and investigate adsorption of argon on these surfaces by using the grand canonical Monte Carlo simulation. In this model, the nongraphitized surface is modeled as a stack of graphene layers with some carbon atoms of the top graphene layer being randomly removed. The percentage of the surface carbon atoms being removed and the effective size of the defect (created by the removal) are the key parameters to characterize the nongraphitized surface. The patterns of adsorption isotherm and isosteric heat are particularly studied, as a function of these surface parameters as well as pressure and temperature. It is shown that the adsorption isotherm shows a steplike behavior on a perfect graphite surface and becomes smoother on nongraphitized surfaces. Regarding the isosteric heat versus loading, we observe for the case of graphitized thermal carbon black the increase of heat in the submonolayer coverage and then a sharp decline in the heat when the second layer is starting to form, beyond which it increases slightly. On the other hand, the isosteric heat versus loading for a highly nongraphitized surface shows a general decline with respect to loading, which is due to the energetic heterogeneity of the surface. It is only when the fluid-fluid interaction is greater than the surface energetic factor that we see a minimum-maximum in the isosteric heat versus loading. These simulation results of isosteric heat agree well with the experimental results of graphitization of Spheron 6 (Polley, M. H.; Schaeffer, W. D.; Smith, W. R. J. Phys. Chem. 1953, 57, 469; Beebe, R. A.; Young, D. M. J. Phys. Chem. 1954, 58, 93). Adsorption isotherms and isosteric heat in pores whose walls have defects are also studied from the simulation, and the pattern of isotherm and isosteric heat could be used to identify the fingerprint of the surface.     

甲烷在石墨化热解碳黑和活性炭上的吸附

为研究影响碳基吸附剂吸附超临界温度气体的主要因素,选择石墨化热解碳黑BP280和Ajax活性炭,分析超临界温度高压甲烷在其上的吸附平衡。应用容积法,在压力0~20.5 MPa、温度253 K~313 K测定甲烷的吸附平衡数据,并由等量吸附线标绘和亨利定律常数确定等量吸附热。引入通用吸附等温方程,再由方程的Langmuir标绘确定最大吸附容量,进而通过方程的线性化计算吸附平衡态中甲烷分子的作用能。结果表明,甲烷在两种吸附剂上的最大吸附容量均随温度而变化,并都小于液态甲烷的密度;甲烷在碳黑和活性炭上的等量吸附热分别为11.9 kJ/mol~12.5 kJ/mol和17.5 kJ/mol~22.5 kJ/mol,体现了两种吸附剂不同的表面能量分布;甲烷分子间作用能随吸附量的变化特点反映了超临界温度甲烷以类似于压缩气体状态聚集的特点和吸附剂结构上的差异。碳基吸附剂的比表面积和微孔容积是影响其储存甲烷容量的重要因素。     

Modeling of adsorption of gases on graphite surfaces accounting for the solid-fluid nonadditivity correction

A series of graphitized carbon blacks have been studied using argon and nitrogen adsorption at their boiling points. Analysis of adsorption isotherms was performed with nonlocal density functional theory (NLDFT) accounting for the Axilrod-Teller equation to describe the effect of nonadditivity of the gas-solid interaction. In our previous study [Ustinov, E. A. J. Chem. Phys. 2010, 132, 194703] we have shown that the nonadditivity effect decreases the attractive component of Ar-Ar interaction in the first molecular layer adjacent to the graphite surface by about 23%. This is a source of a large error (up to 40%) when a standard NLDFT is applied to fitting the low-temperature Ar adsorption isotherm on a graphitized carbon black. A new approach that incorporates the Axilrod-Teller equation into the standard NLDFT diminishes the relative error from 40 to 4%, which suggests that the nonadditivity correction should not be ignored in most adsorption systems including crystalline and amorphous solids. The present study is an extension of our approach to N(2) adsorption isotherms at 77.3 K on graphitized carbon blacks. We show that the approach allows to reliably determine the gas-solid molecular parameters, the gas-solid nonadditivity coefficient, the Henry coefficient, and the specific surface area. The surface areas of different carbon blacks determined with the N(2) at 77.35 K and Ar at 87.29 K are very close to each other, though in the former case the values proved to be slightly smaller presumably due to nonspherical shape of the nitrogen molecule. A comparison with the Brunauer, Emmett, and Teller method is provided.     

Features of nitrogen adsorption on nonporous carbon and silica surfaces in the framework of classical density functional theory

Equilibrium adsorption data of nitrogen on a series of nongraphitized carbon blacks and nonporous silica at 77 K were analyzed by means of classical density functional theory to determine the solid-fluid potential. The behavior of this potential profile at large distance is particularly considered. The analysis of nitrogen adsorption isotherms seems to indicate that the adsorption in the first molecular layer is localized and controlled mainly by short-range forces due to the surface roughness, crystalline defects, and functional groups. At distances larger than approximately 1.3-1.5 molecular diameters, the adsorption is nonlocalized and appears as a thickening of the adsorbed film with increasing bulk pressure in a relatively weak adsorption potential field. It has been found that the asymptotic decay of the potential obeys the power law with the exponent being -3 for carbon blacks and -4 for silica surface, which signifies that in the latter case the adsorption potential is mainly exerted by surface oxygen atoms. In all cases, the absolute value of the solid-fluid potential is much smaller than that predicted by the Lennard-Jones pair potential with commonly used solid-fluid molecular parameters. The effect of surface heterogeneity on the heat of adsorption is also discussed.     

Discussion of thermodynamic data obtained by adsorption gas chromatography of hydrocarbons on carbograph 4

Summary A new graphitized carbon black, Carbograph 4, with a specific surface area of 240 m² g⁻¹ was studied in terms of enthalpy, entropy and free energy of adsorption, determined by gas chromatographic method on a series of alkanes (C₂−C₆) and on benzene. The changes in the isosteric heat and entropy of adsorption when a non-polar stationary phase (squalane) was added to the adsorbent were also investigated. The data obtained are discussed and compared to those present in the literature for other graphitized carbon blacks. Some examples of separations obtained with Carbograph 4 are reported.     

Discussion of thermodynamic data obtained by adsorption gas chromatography of hydrocarbons on a new graphitized carbon black with a high specific surface area

Summary A new graphitized carbon black (Carbograph 5) with a specific surface area (560 m² g⁻¹) greater than those of commerically available graphitized carbons was studied by gas chromatography to determine the enthalpy, entropy, and free energy of adsorption of a series of alkanes (C₂−C₆). The adsorption properties were also investigated by considering changes in the isosteric heats and entropies of adsorption when a nonpolar stationary phase (squalane) was added to the adsorbent. The data obtained are discussed and compared with literature values for other graphitized carbon blacks.     

Adsorption of n-alkanes on a graphitized carbon black and its activation products in dry air

Summary Adsorption of n-alkanes (n-heptane to n-decane) on a graphitized carbon black and its activation products in dry air was carried out by a gas chromatographic technique. Adsorption runs were performed at finite surface coverage, to obtain the surface area and the London component of the surface free energy of the solids from the adsorption isotherms. On the other hand, the adsorption of n-alkanes carried out at zero surface coverage gave the differential heats of adsorption. The results show that, after the activation in dry air of the graphitized carbon black used, the surface area has a linear relationship with the degree of activation, whereas the surface heterogeneity changes in a parabolic fashion.     

Adsorption of benzene on graphitized thermal carbon black: reduction of the quadrupole moment in the adsorbed phase

The performance of intermolecular potential models on the adsorption of benzene on graphitized thermal carbon black at various temperatures is investigated. Two models contain only dispersive sites, whereas the other two models account explicitly for the dispersive and electrostatic sites. Using numerous data in the literature on benzene adsorption on graphitized thermal carbon black at various temperatures, we have found that the effect of surface mediation on interaction between adsorbed benzene molecules must be accounted for to describe correctly the adsorption isotherm as well as the isosteric heat. Among the two models with partial charges tested, the WSKS model of Wick et al. that has only six dispersive sites and three discrete partial charges is better than the very expensive all-atom model of Jorgensen and Severance. Adsorbed benzene molecules on graphitized thermal carbon black have a complex orientation with respect to distance from the surface and also with respect to loading. At low loadings, they adopt the parallel configuration relative to the graphene surface, whereas at higher loadings (still less than monolayer coverage) some molecules adopt a slant orientation to maximize the fluid-fluid interaction. For loadings in the multilayer region, the orientation of molecules in the first layer is influenced by the presence of molecules in the second layer. The data that are used in this article come from the work of Isirikyan and Kiselev, Pierotti and Smallwood, Pierce and Ewing, Belyakova, Kiselev, and Kovaleva, and Carrott et al.     

Modeling of adsorption and nucleation in infinite cylindrical pores by two-dimensional density functional theory

In this paper, we present an analysis of argon adsorption in cylindrical pores having amorphous silica structure by means of a nonlocal density functional theory (NLDFT). In the modeling, we account for the radial and longitudinal density distributions, which allow us to consider the interface between the liquidlike and vaporlike fluids separated by a hemispherical meniscus in the canonical ensemble. The Helmholtz free energy of the meniscus was determined as a function of pore diameter. The canonical NLDFT simulations show the details of density rearrangement at the vaporlike and liquidlike spinodal points. The limits of stability of the smallest bridge and the smallest bubble were also determined with the canonical NLDFT. The energy of nucleation as a function of the bulk pressure and the pore diameter was determined with the grand canonical NLDFT using an additional external potential field. It was shown that the experimentally observed reversibility of argon adsorption isotherms at its boiling point up to the pore diameter of 4 nm is possible if the potential barrier of 22kT is overcome due to density fluctuations.     

The critical size of clusters of water molecules on a carbon surface

The experimental adsorption isotherms of water and nitrogen vapors on graphitized carbonaceous adsorbents with large pore size prepared from ultradispersive technical carbon black have been compared with those on the surface of non-porous graphitized carbon black. The saturation value of water vapor adsorption has been shown to be proportional to the concentration of primary adsorption centers. At low concentrations of these centers the saturation value corresponds to the formation of fractions of a dense monolayer on the surface. The maximum size of clusters of water molecules on a carbonaceous adsorbent surface has been estimated.Translated fromIzvestiya Akademii Nauk, Seriya Khimicheskaya, No. 1, pp. 54–56, January, 1993.     

Modeling of Adsorption in Finite Cylindrical Pores by Means of Density Functional Theory

Adsorption of argon at its boiling point in finite cylindrical pores is considered by means of the non-local density functional theory (NLDFT) with a reference to MCM-41 silica. The NLDFT was adjusted to amorphous solids, which allowed us to quantitatively describe argon adsorption isotherm on nonporous reference silica in the entire bulk pressure range. In contrast to the conventional NLDFT technique, application of the model to cylindrical pores does not show any layering before the phase transition in conformity with experimental data. The finite pore is modeled as a cylindrical cavity bounded from its mouth by an infinite flat surface perpendicular to the pore axis. The adsorption of argon in pores of 4 and 5 nm diameters is analyzed in canonical and grand canonical ensembles using a two-dimensional version of NLDFT, which accounts for the radial and longitudinal fluid density distributions. The simulation results did not show any unusual features associated with accounting for the outer surface and support the conclusions obtained from the classical analysis of capillary condensation and evaporation. That is, the spontaneous condensation occurs at the vapor-like spinodal point, which is the upper limit of mechanical stability of the liquid-like film wetting the pore wall, while the evaporation occurs via a mechanism of receding of the semispherical meniscus from the pore mouth and the complete evaporation of the core occurs at the equilibrium transition pressure. Visualization of the pore filling and empting in the form of contour lines is presented.     

On the physical adsorption of gases on carbon materials from molecular simulation

In this paper we present a series of work covering a range of aspects relating molecular simulation to experiment. The importance of surface mediation type effects to the adsorption of simple and complex gases is demonstrated. Coupled with the adsorption of simple gases is their projection area when used for surface area determination. The pressure dependence of a projection area is demonstrated for argon at 77 and 87.3 K. A simple model is used to account for the degree of graphitisation of a surface is demonstrated and used to account for the isosteric heat behaviour of non-graphitised carbon blacks. Turning from surfaces to porous solids, an alternative treatment of experiment data (either sub or super critical) is presented that avoids the ambiguity of excess amounts adsorbed. Using this method one is able to obtain pore size distributions and amounts adsorbed without relying on such things as helium expansion volumes. Since this type of method is usually applied to composite solids we also demonstrate the correct method for calculating the heat of adsorption using independent sets of simulations. The final topic covered in this paper is an example of the information that can be gained from the heat capacity of an adsorbed phase.     

Thermodynamic analysis of ordered and disordered monolayer of argon adsorption on graphite

We presented a detailed thermodynamic analysis of argon adsorption on a graphitized carbon black with a kinetic Monte Carlo scheme. In this study, we particularly paid attention to the formation of a hexagonal two-dimensional molecular layer on a graphite surface and discuss conditions of its stability and thermodynamic properties of the adsorbed phase as a function of loading. It is found that the simulation results are substantially affected by the dimensions of the simulation box when the monolayer forms a hexagonal ordered structure. This is due to the fact that the lattice constant is constrained by the dimensions of the surface. To circumvent this, we presented a thermodynamic technique, which allows for the variation of the box size as a function of loading, to determine the "intrinsic" lattice constant (rather than apparent average value because of the fixed dimensions of the simulation box) and the thermodynamic functions for the adsorbed phase: the Helmholtz free energy, the chemical potential, and the surface tension. The tangential and normal pressures as a function of the distance from the surface are also discussed.     

Two-dimensional order-disorder transition of argon monolayer adsorbed on graphitized carbon black: kinetic Monte Carlo method

We present results of application of the kinetic Monte Carlo technique to simulate argon adsorption on a graphite surface at temperatures below and above the triple point. We show that below the triple point the densification of the adsorbed layer with loading results in the rearrangement of molecules to form a hexagonal structure, which is accompanied by the release of an additional heat, associated with this disorder-order transition. This appears as a spike in the plot of the heat of adsorption versus loading at the completion of a monolayer on the surface. To describe the details of the adsorbed phase, we analyzed thermodynamic properties and the effects of temperature on the order-disorder transition of the first layer.     

A model of the two-stage condensation mechanism of water adsorption on nonporous carbon adsorbents

The mechanism of adsorption of water molecules on nonporous carbon adsorbents has been suggested in terms of two different states of adsorbed water; stretched liquid water and water that occupies an intermediate state between the liquid and vapor. Two stages of adsorption were distinguished: condensation and pre-condensation that assumes the formation of molecular associates. The BET model was used to describe the pre-condensation stage. The equations of the adsorption isotherm for water vapor in the region of condensation process and the expression for the determination of the specific hydrophilic surface of adsorbents were found. Examination of the experimental data on adsorption of water vapor on nongraphitized samples of carbon adsorbents shows that in the region of polymolecular adsorption, all isotherms fall into a common curve determined by the equation of the stretched liquid film and can be calculated regardless of the properties of individual liquid water. The equation for adsorption of water vapor on the hydrophobic surface was obtained. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 10, pp. 1933–1939, October, 1998.     

Cooperative and Competitive Adsorption of Ethylene,Ethane, Nitrogen and Argon on Graphitized Carbon Black and in Slit Pores

In this paper we investigate the mixture adsorption of ethylene, ethane, nitrogen and argon on graphitized thermal carbon black and in slit pores by means of the Grand Canonical Monte Carlo simulations. Pure component adsorption isotherms on graphitized thermal carbon black are first characterized with the GCMC method, and then mixture simulations are carried out over a wide range of pore width, temperature, pressure and composition to investigate the cooperative and competitive adsorption of all species in the mixture. Results of mixture simulations are compared with the experimental data of ethylene and ethane (Friederich and Mullins, 1972) on Sterling FTG-D5 (homogeneous carbon black having a BET surface area of 13 m²/g) at 298 K and a pressure range of 1.3–93 kPa. Because of the co-operative effect, the Henry constant determined by the traditional chromatography method is always greater than that obtained from the volumetric method.     

Adsorption of Carbon Dioxide on Activated Carbon

The adsorption of CO2 on a raw activated carbon A and three modified activated carbon samples B, C, and D at temperatures ranging from 303 to 333 K and the thermodynamics of adsorption have been investigated using a vacuum adsorption apparatus in order to obtain more information about the effect of CO2 on removal of organic sulfur-containing compounds in industrial gases. The active ingredients impregnated in the carbon samples show significant influence on the adsorption for CO2 and its volumes adsorbed on modified carbon samples B, C, and D are all larger than that on the raw carbon sample A. On the other hand, the physical parameters such as surface area, pore volume, and micropore volume of carbon samples show no influence on the adsorbed amount of CO2. The Dubinin-Radushkevich (D-R) equation was the best model for fitting the adsorption data on carbon samples A and B, while the Preundlich equation was the best fit for the adsorption on carbon samples C and D. The isosteric heats of adsorption on carbon samples A, B, C, and D derived from the adsorption isotherms using the Clapeyron equation decreased slightly increasing surface loading. The heat of adsorption lay between 10.5 and 28.4 kJ/mol, with the carbon sample D having the highest value at all surface coverages that were studied. The observed entropy change associated with the adsorption for the carbon samples A, B, and C (above the surface coverage of 7 ml/g) was lower than the theoretical value for mobile adsorption. However, it was higher than the theoretical value for mobile adsorption but lower than the theoretical value for localized adsorption for carbon sample D.