Comparative features of ar adsorption on smooth and amorphous surfaces examined by density functional theory

We analyze argon adsorption isotherms and isosteric heat of adsorption on graphitized and nongraphitized carbon black and silica surfaces by means of nonlocal density functional theory (NLDFT). It is shown that in the case of graphitized carbon black the behavior of the adsorbed phase is nearly identical to that in the bulk phase at a distance larger than about 3-4 molecular diameters from the surface. At a smaller distance argon forms solid-like molecular layers at a temperature at least 3.5 K above the triple point, with the interlayer distance being markedly smaller than the argon collision diameter. In the case of defected or amorphous surfaces adsorbed argon is liquid-like below its triple point. Our extension of the Tarazona NLDFT to amorphous solids (NLDFT-AS) and the Kierlik and Rosinberg version of NLDFT excellently fit argon adsorption isotherms and properly predict the isosteric heat of adsorption. We showed that the surface roughness affects the calculated heat of adsorption, which allowed us to adjust the width of the diffuse zone of the nongraphitized carbon black and the silica surface.     

Probing the mechanism of water adsorption in carbon micropores with multitemperature isotherms and water preadsorption experiments

The phenomenon of water adsorption in carbon micropores is examined through the study of water adsorption equilibrium in molecular sieving carbon. Adsorption and desorption isotherms are obtained over a wide range of concentrations from less than 0.1% to beyond 80% of the vapor pressure. Evidence is provided in support of a proposed bimodal water adsorption mechanism that involves the interaction of water molecules with functional groups at low relative pressures and the adsorption of water molecules between graphene layers at higher pressures. Decomposition of the equilibrium isotherm data through application of the extended cooperative multimolecular sorption theory, together with favorable quantitative comparison, provides support for the proposed adsorption mechanism. Additional support is obtained from a multitemperature study of water equilibrium. Temperatures of 20, 50, and 60 degrees C were probed in this investigation in order to provide isosteric heat of adsorption data for water interaction with the carbon molecular sieve. At low loading, the derived isosteric heat of adsorption is estimated to be 69 kJ/mol. This value is indicative of the adsorption of water to functional groups. At higher loading, the isosteric heat of adsorption decreases with increasing loading and approaches the heat of condensation, indicative of adsorption between graphene layers. Further support for the proposed adsorption mechanism is derived from carbon dioxide adsorption experiments on carbon molecular sieve that is preadsorbed with various amounts of water. Significant exclusion of carbon dioxide occurs, and a quantitative analysis that is based on the proposed bimodal water adsorption mechanism is employed in this investigation.     

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.     

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.     

Adsorption isotherms of 2,2,4-trimethylpentane and toluene vapors on hydrocarbon adsorber and light-off catalyst

Two monolithic hydrocarbon adsorbers and a monolithic light-off catalyst were selected as adsorbents, and the adsorptive capacity of a hydrocarbon for the adsorbents was measured by using a precise volumetric adsorption apparatus. 2,2,4-Trimethylpentane and toluene vapors were chosen as adsorbates. Equilibrium experiments were carried out at three different temperatures of 303.15, 323.15, and 343.15 K. Adsorption data of each hydrocarbon was fitted to the well-known isotherms such as the Langmuir equation and the Freundlich equation. The Freundlich isotherm predicted equilibrium data better than the Langmuir isotherm. Furthermore, the surface energetic heterogeneity of the adsorbents was evaluated using the isosteric heat of adsorption based on Clausius-Clapeyron equation. The surface energetic heterogeneity of the adsorbents depended on the precious metal (PM) loading and H-ZSM5 loading.     

Monte carlo simulation and pore-size distribution analysis of the isosteric heat of adsorption of methane in activated carbon

The isosteric heat of adsorption of methane in an activated carbon adsorbent has been modeled by Monte Carlo simulation, using a pore-size distribution (PSD) to relate simulation results for pores of different sizes to the experimental adsorbent. Excellent fits were obtained between experimental and simulated isosteric heats of adsorption of methane in BPL activated carbon. The PSD was then used to predict the adsorption of methane and ethane in the same carbon adsorbent, with good results. The PSD derived from isosteric heat data was shown to be richer in information than PSDs obtained by the more conventional method of fitting to isotherm data.     

Explanation of the unusual peak of calorimetric heat in the adsorption of nitrogen, argon and methane on graphitized thermal carbon black

Heats of adsorption and adsorption isotherms of argon, nitrogen and methane on a perfect graphitic surface and a defective graphitic surface are studied with a Grand Canonical Monte Carlo Simulation (GCMC). For the perfect surface, the isosteric heat versus loading shows a typical pattern of adsorption of simple fluids on graphite. Depending on adsorbate, degree of graphitization and temperature, a spike in the heat curve versus loading is observed when the first layer is mostly covered with adsorbate molecules. The heat spike is observed for argon and nitrogen at 77 K while for argon at 87.3 K it is no longer present. These simulation results are consistent with the experimental data of J. Rouquerol, S. Partyka and F. Rouquerol, J. Chem. Soc., Faraday Trans. 1, 1977, 73, 306. In the case of methane we observe heat spikes at low temperatures, 84.5, 92.5 and 104 K. The heat spike shifts to higher loading with temperature and it then disappears at high temperatures. These observations are in qualitative agreement with the experimental data of A. Inaba, Y. Koga and J. A. Morrison, J. Chem. Soc., Faraday Trans. 2, 1986, 82, 1635. In all cases where heat spikes are observed, the GCMC simulation results indicate that the heat spike is associated with the squeezing of molecules into the already dense first layer, and the rearrangement of molecules to form a highly structured fluid of this layer. While this squeezing into the first layer is happening, molecules continue to adsorb onto the relatively sparse second layer.     

Hierarchical modeling N2 adsorption on the surface of and within a C60 crystal: from quantum mechanics to molecular simulation

The adsorption of N(2) on the surface of, and within, a C(60) face-centered cubic crystal has been studied using a hierarchical approach. First, an ab initio potential between N(2) and C(60) is obtained from a recently developed quantum mechanical hybrid method, and then the adsorption behavior is predicted using Monte Carlo simulation. On the crystal surface, N(2) adsorption isotherm at 77.3 K is of type II. The adsorption simulated with the ab initio potential is slightly greater than that with the empirical Steele potential derived from experimental N(2) adsorption on planar graphite, and both are in fairly good agreement with measured results. With increasing pressure, N(2) molecules are found to sequentially occupy three favorable sites: the octahedral sites, the tetrahedral sites, and the top of C(60) molecules. Finally multiple layers form and wetting occurs as the bulk N(2) saturation pressure is reached. The isosteric heat of adsorption exhibits two maxima and finally approaches the enthalpy of vaporization of bulk N(2). Within the crystal, the N(2) adsorption isotherm at 77.3 K is of type I, and the use of ab initio potential leads to significantly greater adsorption than the Steele potential. N(2) molecules are observed to intercalate only the octahedral sites, and the isosteric heat of adsorption is nearly a constant. As in our previous work of N(2) and O(2) adsorption in the C(168) schwarzite (Jiang et al. J. Phys. Chem. B 2004, 108, 9852), this work demonstrates the importance of an accurate adsorbate-adsorbent interaction potential in the determination of gas adsorption behavior.     

Grand canonical monte carlo simulation study of water adsorption in silicalite at 300 K

This molecular simulation work focuses on the adsorption of water in a priori hydrophobic silicalite-1, a microporous ordered silica. The water-water interactions are described with the SPC model, while water-silica interactions are calculated in the framework of the PN-TrAZ model. The water adsorption isotherm at 300 K, the configurational energies, and the isosteric heat of adsorption are calculated by the grand canonical Monte Carlo (GCMC) simulation method. The thermodynamic integration scheme allows one to calculate the grand potential along the adsorption isotherm. The adsorption results are compared with experiments, showing only qualitative agreement. Indeed, the simulations do not reproduce the expected hydrophobicity of silicalite (Eroshenko, V.; Regis, R.-C.; Soulard, M.; Patarin, J. C. R. Phys. 2002, 3, 111). This indicates that common models used to describe confined polar molecules are far from being operative. In this work, it is shown, on the basis of periodic ab initio calculations, that confined water molecules in silicalite have a dipole value roughly 10% smaller than that in the bulk liquid phase, indicating that the environment felt by a confined water molecule in silicalite pores is not equivalent to that in the bulk liquid. This suggests that effective intermolecular potentials parametrized for bulk water are inefficient to describe ultraconfined water molecules. Reducing the SPC water dipole moment by 5% (i.e., decreasing water partial charges in magnitude) in GCMC calculations does allow reproducing the experimental water/silicalite isotherm at 300 K.     

Heterogeneity of multiwalled carbon nanotubes based on adsorption of simple aromatic compounds from aqueous solutions

The surface heterogeneity of multiwalled carbon nanotubes (MWCNTs) is studied on the basis of adsorption isotherms from dilute aqueous phenol and dopamine solutions at various pH values. The generalized Langmuir–Freundlich isotherm equation was applied to investigate the cooperative effect of the surface heterogeneity and the lateral interactions between the adsorbates. The theoretical isosteric heats of adsorption were obtained assuming that the heat of adsorption profile reveals both the energetic heterogeneity of the adsorption system and the strength of the interactions between the neighboring molecules. The adsorption energy distribution functions were calculated by using algorithm based on a regularization method. The great advantage of this method is that the regularization makes no assumption about the shape of the obtained energy distribution functions. Analysis of the isosteric heats of adsorption for MWCNTs showed that the influence of the surface heterogeneity is much stronger than the role of the lateral interactions. The most typical adsorption heat is 20–22 kJ/mol for both phenol and dopamine. After purification of nanotubes, heat value for phenol dropped to 16–17 kJ/mol. The range of the energy distribution is only slightly influenced by the surface chemistry of the nanotubes in the aqueous conditions.     

On the Henry constant and isosteric heat at zero loading in gas phase adsorption

The Henry constant and the isosteric heat of adsorption at zero loading are commonly used as indicators of the strength of the affinity of an adsorbate for a solid adsorbent. It is assumed that (i) they are observable in practice, (ii) the Van Hoff's plot of the logarithm of the Henry constant versus the inverse of temperature is always linear and the slope is equal to the heat of adsorption, and (iii) the isosteric heat of adsorption at zero loading is either constant or weakly dependent on temperature. We show in this paper that none of these three points is necessarily correct, first because these variables might not be observable since they are outside the range of measurability; second that the linearity of the Van Hoff plot breaks down at very high temperature, and third that the isosteric heat versus loading is a strong function of temperature. We demonstrate these points using Monte Carlo integration and Monte Carlo simulation of adsorption of various gases on a graphite surface. Another issue concerning the Henry constant is related to the way the adsorption excess is defined. The most commonly used equation is the one that assumes that the void volume is the volume extended all the way to a boundary passing through the centres of the outermost solid atoms. With this definition the Henry constant can become negative at high temperatures. Although adsorption at these temperatures may not be practical because of the very low value of the Henry constant, it is more useful to define the Henry constant in such a way that it is always positive at all temperatures. Here we propose the use of the accessible volume; the volume probed by the adsorbate when it is in nonpositive regions of the potential, to calculate the Henry constant.     

水在NanZSM-5型分子筛中吸附的研究: 分子模拟

利用分子动力学(MD)模拟退火的方法和巨正则系综Monte Carlo模拟方法(GCMC)研究了水在NanZSM-5型分子筛中的吸附行为, 计算结果与文献中报道的实验结果吻合较好. 在此基础上, 进一步预测了水在不同硅铝比的NanZSM-5型分子筛中的吸附性质, 计算结果显示: 分子筛骨架上的硅铝比会显著影响水分子的吸附量和吸附等温线, 随着硅铝比的降低, 水的吸附量增加; 水分子的吸附位置是在钠离子和铝原子的周围, 平均每个钠离子周围吸附4个水分子, 而当水的吸附量增大时, 水分子与分子筛骨架上的氧原子之间发生了氢键作用; 在吸附量相同的条件下, 水的吸附热随着硅铝比的降低而升高.     

Semiempirical model for adsorption of polyatomics

Multisite-occupancy adsorption is described by using a new formalism based on the occupation balance approximation (Roma, F.; Ramirez-Pastor, A. J.; Riccardo, J. L. J. Chem. Phys. 2001, 114, 10932). In this framework, the adsorption isotherm is characterized by a correction function C, which relates to the conditional probability of finding the ith empty site to a lattice with i - 1 already vacant sites. A simple semiempirical adsorption isotherm is proposed by approximating C as a combination of the correction functions corresponding to exact 1-D calculations and the Guggenheim-DiMarzio approximation, with adequate weights. Results are compared with corresponding ones from Monte Carlo simulations.     

The effects of energy sites on adsorption of Lennard-Jones fluids and phase transition in carbon slit pore of finite length a computer simulation study

A Monte Carlo simulation method is used to study the effects of adsorption strength and topology of sites on adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length. Argon is used as a model adsorbate, while the adsorbent is modeled as a finite carbon slit pore whose two walls composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. Impurities having well depth of interaction greater than that of carbon atom are assumed to be grafted onto the surface. Different topologies of the impurities; corner, centre, shell and random topologies are studied. Adsorption isotherms of argon at 87.3 K are obtained for pore having widths of 1, 1.5 and 3 nm using a Grand Canonical Monte Carlo simulation (GCMC). These results are compared with isotherms obtained for infinite pores. It is shown that the surface heterogeneity affects significantly the overall adsorption isotherm, particularly the phase transition. Basically it shifts the onset of adsorption to lower pressure and the adsorption isotherms for these four impurity models are generally greater than that for finite pore. The positions of impurities on solid surface also affect the shape of the adsorption isotherm and the phase transition. We have found that the impurities allocated at the centre of pore walls provide the greatest isotherm at low pressures. However when the pressure increases the impurities allocated along the edges of the graphene layers show the most significant effect on the adsorption isotherm. We have investigated the effect of surface heterogeneity on adsorption hysteresis loops of three models of impurity topology, it shows that the adsorption branches of these isotherms are different, while the desorption branches are quite close to each other. This suggests that the desorption branch is either the thermodynamic equilibrium branch or closer to it than the adsorption branch.     

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

为研究影响碳基吸附剂吸附超临界温度气体的主要因素,选择石墨化热解碳黑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,体现了两种吸附剂不同的表面能量分布;甲烷分子间作用能随吸附量的变化特点反映了超临界温度甲烷以类似于压缩气体状态聚集的特点和吸附剂结构上的差异。碳基吸附剂的比表面积和微孔容积是影响其储存甲烷容量的重要因素。     

Comparative simulation study of nitrogen and ammonia adsorption on graphitized and nongraphitized carbon blacks

Grand canonical Monte Carlo simulation is used to study the adsorption of nitrogen at 77 K and ammonia at 240 K to represent weakly polar and polar molecules, respectively, on infinite and finite graphite surfaces. These graphite surfaces were modeled with different percentages of carbons removed (defects) from the top graphite layer. Increasing the number of defects increases the adsorption and the isosteric heat of nitrogen at low pressure. At moderate pressures the amount adsorbed is less due to the disruption in the packing of the nitrogen in the first layer. In contrast, the adsorption of ammonia at all pressures is reduced as the percentage of defects is increased. This is due to the disruption in ammonia bonding caused by the defects. The condensation-like step change in the ammonia isotherm on the perfect graphite surface is not observed for any of these surfaces with defects even for the case of only 10% defects. At high percentage of defects the adsorption isotherm is close to Henry law behavior for much of the pressure range. The adsorption on finite surfaces shows that the amount adsorbed for both molecules decreases compared with that of the infinite surfaces, resulting from interaction potentials with the surface and other fluid molecules at the edge. The decrease is much greater for the ammonia adsorption because the bonding between ammonia molecules is disrupted, meaning that the adsorption cannot follow the mechanism of condensation seen for the infinite surface.     

Room temperature water splitting at the surface of magnetite

An array of surface science measurements has revealed novel water adsorption behavior at the Fe(3)O(4)(001) surface. Following room temperature exposure to water, a low coverage of hydrogen atoms is observed, with no associated water hydroxyl group. Mild annealing of the hydrogenated surface leads to desorption of water via abstraction of surface oxygen atoms, leading to a reduction of the surface. These results point to an irreversible splitting of the water molecule. The observed phenomena are discussed in the context of recent DFT calculations (Mulakaluri, N.; Pentcheva, R.; Scheffler, M. J. Phys. Chem. C 2010, 114, 11148), which show that the Jahn-Teller distorted surface isolates adsorbed H in a geometry that could kinetically hinder recombinative desorption. In contrast, the adsorption geometry facilitates interaction between water hydroxyl species, which are concluded to leave the surface following a reactive desorption process, possibly via the creation of O(2).