Gas adsorption in active carbons and the slit-pore model 2: Mixture adsorption prediction with DFT and IAST

We use a fast density functional theory (a "slab-DFT") and the polydisperse independent ideal slit-pore model to predict gas mixture adsorption in active carbons. The DFT is parametrized by fitting to pure gas isotherms generated by Monte Carlo simulation of adsorption in model graphitic slit-pores. Accurate gas molecular models are used in our Monte Carlo simulations with gas-surface interactions calibrated to a high surface area carbon, rather than a low surface area carbon as in all previous work of this type, as described in part 1 of this work. We predict the adsorption of binary mixtures of carbon dioxide, methane, and nitrogen on two active carbons up to about 30 bar at near-ambient temperatures. We compare two sets of results; one set obtained using only the pure carbon dioxide adsorption isotherm as input to our pore characterization process, and the other obtained using both pure gas isotherms as input. We also compare these results with ideal adsorbed solution theory (IAST). We find that our methods are at least as accurate as IAST for these relatively simple gas mixtures and have the advantage of much greater versatility. We expect similar results for other active carbons and further performance gains for less ideal mixtures.     

Molecular simulation investigation into the performance of Cu-BTC metal-organic frameworks for carbon dioxide-methane separations

We report a molecular simulation study for Cu-BTC metal-organic frameworks as carbon dioxide-methane separation devices. For this study we have computed adsorption and diffusion of methane and carbon dioxide in the structure, both as pure components and mixtures over the full range of bulk gas compositions. From the single component isotherms, mixture adsorption is predicted using the ideal adsorbed solution theory. These predictions are in very good agreement with our computed mixture isotherms and with previously reported data. Adsorption and diffusion selectivities and preferential sitings are also discussed with the aim to provide new molecular level information for all studied systems.     

Monte Carlo Simulations of Adsorbed Solutions in Heterogeneous Porous Materials

We present results of a Monte Carlo simulation study of binary mixtures of ethane and methane in silica gel. The molecular model treats the adsorbent as a matrix of silica microspheres. The adsorption isotherms, adsorption selectivities and isosteric heats of adsorption have been determined for these systems. The results are compared with predictions from the ideal adsorbed solution (IAS) theory and with experiment. The heats of adsorption are accurately described by the IAS theory. The adsorption isotherms are accurately described by the IAS theory at low bulk pressure but the IAS theory overpredicts the density at high bulk pressure. This latter effect is opposite to that observed in bulk mixtures of this type where nonidealities generally lead to a density increase on mixing. The pressure dependence of the selectivity does not exhibit a maximum at low pressure. We discuss this effect in terms of the adsorbent microstructure.     

Effects of molecular siting and adsorbent heterogeneity on the ideality of adsorption equilibria

The ideal adsorbed solution (IAS) theory is the benchmark for the prediction of mixed-gas adsorption equilibria from pure-component isotherms. In this work, we use atomistic grand canonical Monte Carlo simulations to test the effects of molecular siting and adsorbent energetic heterogeneity on the applicability of the IAS theory. Pure-component isotherms generated by atomistic simulation are used to predict binary isobaric isotherms using the IAS theory. These predicted isotherms are compared with those obtained by a full atomistic simulation of the binary mixture. Binary mixtures of argon, methane, and CF4 in silicalite are found to obey IAS theory, while benzene/methane and cyclohexane/methane in silicalite are nonideal. The mixture of argon and CF4 is ideal despite the large difference in the sizes of the two species. This contradicts previous hypotheses in the literature, which state that mixtures of species of unequal size do not adsorb ideally. The nonideal behavior of the benzene/methane and cyclohexane/methane systems occurs because of adsorbent heterogeneity in these systems, which depends on both sorbent and sorbate. In addition, we use a lattice gas model with parameters derived from atomistic simulation to demonstrate analytically that a sufficiently energetically heterogeneous adsorbent will result in the breakdown of IAS theory even in the absence of interactions between sorbates.     

Examining the accuracy of ideal adsorbed solution theory without curve-fitting using transition matrix Monte Carlo simulations

Ideal adsorbed solution theory (IAST) is a well-known approach to predicting multicomponent adsorption isotherms in microporous materials from experimental or simulation data for single-component adsorption. A limitation in practical applications of IAST is that useful calculations often require extrapolation of fitted single-component isotherms beyond the range for which data are available. We introduce a molecular simulation approach in which the intrinsic accuracy of IAST can be examined in a context that avoids any need to perform curve fitting with single-component data. Our approach is based on using transition matrix Monte Carlo to define single-component adsorption isotherms for arbitrary bulk-phase pressures from a single simulation. We apply our approach to several light gas mixtures in silica zeolites and a carbon nanotube to examine the intrinsic accuracy of IAST for these model systems.     

Molecular simulation of binary mixture adsorption in buckytubes and MCM-41

MCM-41 and buckytubes are novel porous materials with controllable pore sizes and narrow pore size distributions. Buckytubes are carbon tubes with internal diameters in the range 1–5 urn. The structure of each tube is thought to be similar to one or more graphite sheets rolled up in a helical manner. MCM-41 is one member of a new family of highly uniform mesoporous silicate materials produced by Mobil, whose pore size can be accurately controlled in the range 1.5–10 nm. We present grand canonical Monte Carlo (GCMC) simulations of single fluid and binary mixture adsorption in a model buckytube, and nonlocal density functional theory (DFT) calculations of trace pollutant separation in a range of buckytubes and MCM-41 pores. Three adsorbed fluids are considered; methane, nitrogen and propane. The GCMC studies show that the more strongly adsorbed pure fluid is adsorbed preferentially from an equimolar binary mixture. Ideal adsorbed solution theory (IAST) is shown to give good qualitative agreement with GCMC when predicting binary mixture separations. The DFT results demonstrate the very large increases in trace pollutant separation that can be achieved by tuning the pore size, structure, temperature and pressure of the MCM-41 and buckytube adsorbent systems to their optimal values.     

高比表面活性碳微球分离H2中少量CO2

以实验数据为依据, 结合双Langmuir模型研究了用高比表面活性碳微球材料分离H2中少量CO2的行为. 在实验中, 用高精度的IGA-003重力吸附仪测定了温度为298、273 和268 K, 压力在0-1.8 MPa范围内CO2、H2及n(CO2):n(H2)=1:9混合物在活性碳微球中的吸附等温线. 比较不同吸附模型的计算结果与实验数据, 结果表明, 双Langmuir模型与实验结果拟合得较好; 而且通过结合理想吸附溶液理论, 该模型可以准确地计算不同的混合物体系(包括H2-CO2体系)的吸附量和吸附选择性. 利用该模型求解了不同温度下各组分的分吸附量, 得到了CO2的吸附选择性;在268 K和1.7 MPa下, CO2的吸附选择性可达到73.4, 表明活性碳微球是一种优秀的吸附H2中少量CO2的材料.     

Adsorption equilibria of binary gas mixtures on graphitized carbon black

Adsorption equilibria for binary gas mixtures (methane-carbon dioxide, methane-ethane, and carbon dioxide-ethane) on the graphitized carbon black STH-2 were measured by the open flow method at 293.2 K. The experimental pressure range was (0 to 1.6) MPa. The extended Langmuir (EL) model and the ideal adsorption solution theory (IAST) have been adopted to predict the equilibria of binary gas mixtures. The results indicate that gas mixtures adsorbed on the homogeneous surface of STH-2 exhibit the nonideal behavior, which is mainly induced by adsorbate-adsorbate interactions. The real adsorption solution theory (RAST) has been used to analyze the property of the adsorbed mixtures. The activity coefficients have been correlated with the Wilson equation. The investigation demonstrates that the nonideality of adsorbed phase is completely dissimilar with the bulk liquid phase. The adsorption of the heavier component would benefit the adsorption of the lighter component.     

Surface phase separation in complex mixed adsorbing systems: an interface-bulk coupling effect

The interfacial thermodynamics and structure of ternary mixtures of the type A+B+solvent are investigated. According to the Gibbs phase rule, the coupling between the bulk phase and the interfacial region-which is related to the reversibility of the adsorption of the corresponding species-is a determinant as to whether phase separation can be observed at the interface. For an n-component adsorbing solution, at least one of the species has to adsorb irreversibly over the experimental time scales in order not to fix more intensive variables than those required to observe surface phase separation. We present results for a lattice model planar interface consisting of the ternary mixture A+B+solvent. The solvent molecules and the type A molecules have fixed chemical potentials at the interface since they are equilibrated with a bulk solution. In contrast, the type B molecules are irreversibly adsorbed at the interface and do not equilibrate with the bulk. Mean-field theory is compared with Monte Carlo simulation. Interestingly, the spinodal line in the interaction-composition plane shows a reentrant on the B-rich phase side. We discuss the implications of these results for surface phase separation of adsorbing mixtures of proteins and low-molecular-weight surfactants.     

New thermodynamic/electrostatic models of adsorption and tension equilibria of aqueous ionic surfactant mixtures: application to sodium dodecyl sulfate/sodium dodecyl sulfonate systems

The nonideal adsorbed solution (NAS) theory has been formally extended to adsorption at the air/water interface from aqueous mixtures of ionic surfactants, explicitly accounting for the surface potential of the adsorbed monolayer with the Gouy-Chapman theory. This new ionic NAS (iNAS) theory is thermodynamically consistent and, when coupled to a micellization model, is valid for concentrations below and above the mixed cmc. Counterion binding is incorporated into the model using two fractional binding parameters, beta(sigma) for the adsorbed monolayer and beta(m) for the micelles. The regular solution theory is used to model the nonideal interactions within the adsorbed monolayer and within the mixed micelles. New tension data for an equimolar mixture of sodium dodecyl sulfate (SDS) and sodium dodecyl sulfonate (SDSn) at two salinities fit this model well when mixing is ideal. The total surface densities, the surface compositions, and the surface potentials for the mixed monolayers are calculated. When there is no added salt, at total surfactant concentrations below the mixed cmc, the adsorbed monolayer is enriched in SDSn, but at total concentrations at and above the mixed cmc, the adsorbed monolayer is nearly an equimolar mixture. In the presence of 100 mM NaCl, the adsorbed monolayer is nearly an equimolar mixture, independent of the total surfactant concentration.     

Testing the accuracy of correlations for multicomponent mass transport of adsorbed gases in metal-organic frameworks: diffusion of H2/CH4 mixtures in CuBTC

Mass transport of chemical mixtures in nanoporous materials is important in applications such as membrane separations, but measuring diffusion of mixtures experimentally is challenging. Methods that can predict multicomponent diffusion coefficients from single-component data can be extremely useful if these methods are known to be accurate. We present the first test of a method of this kind for molecules adsorbed in a metal-organic framework (MOF). Specifically, we examine the method proposed by Skoulidas, Sholl, and Krishna (SSK) ( Langmuir, 2003, 19, 7977) by comparing predictions made with this method to molecular simulations of mixture transport of H 2/CH 4 mixtures in CuBTC. These calculations provide the first direct information on mixture transport of any species in a MOF. The predictions of the SSK approach are in good agreement with our direct simulations of binary diffusion, suggesting that this approach may be a powerful one for examining multicomponent diffusion in MOFs. We also use our molecular simulation data to test the ideal adsorbed solution theory method for predicting binary adsorption isotherms and a method for predicting mixture self-diffusion coefficients.     

Modeling adsorption in metal-organic frameworks with open metal sites: propane/propylene separations

We present a new approach for modeling adsorption in metal-organic frameworks (MOFs) with unsaturated metal centers and apply it to the challenging propane/propylene separation in copper(II) benzene-1,3,5-tricarboxylate (CuBTC). We obtain information about the specific interactions between olefins and the open metal sites of the MOF using quantum mechanical density functional theory. A proper consideration of all the relevant contributions to the adsorption energy enables us to extract the component that is due to specific attractive interactions between the π-orbitals of the alkene and the coordinatively unsaturated metal. This component is fitted using a combination of a Morse potential and a power law function and is then included into classical grand canonical Monte Carlo simulations of adsorption. Using this modified potential model, together with a standard Lennard-Jones model, we are able to predict the adsorption of not only propane (where no specific interactions are present), but also of propylene (where specific interactions are dominant). Binary adsorption isotherms for this mixture are in reasonable agreement with ideal adsorbed solution theory predictions. We compare our approach with previous attempts to predict adsorption in MOFs with open metal sites and suggest possible future routes for improving our model.     

Adsorption and selectivity of carbon dioxide with methane and nitrogen in slit-shaped carbonaceous micropores: Simulation and experiment

We have used the grand canonical Monte Carlo method to study the adsorption and selectivity of mixtures of carbon dioxide with methane and nitrogen at high (i.e., ambient) temperatures in model slit pores with graphitic surfaces. Experimental data, including new high pressure measurements for carbon dioxide and methane on a non-porous graphitic standard, were used to test the potential models. The mixture simulations predict that carbon dioxide is preferentially adsorbed in both systems. The results are discussed in terms of competing energetic and entropic effects and the underlying molecular mechanisms.     

Fluids in nanospaces: molecular simulation studies to find out key mechanisms for engineering

We have analyzed various phenomena that occur in nanopores, focusing on elucidating their key mechanisms, to advance the effective engineering use of nanoporous materials. As ideal experimental systems, molecular simulations can effectively provide information at the molecular level that leads to mechanistic insight. In this short review, several of our recent results are presented. The first topic is the critical point depression of Lennard-Jones fluid in silica slit pores due to finite size effects, studied by our original Monte Carlo (MC) technique. We demonstrate that the first layers of adsorbed molecules in contact with the pore walls act as a “fluid wall” and impose extra finite size effects on the fluid confined in the central portion of the pore. We next present a new kernel for pore size distribution (PSD) analysis, based entirely on molecular simulation, which consists of local isotherms for nitrogen adsorption in carbon slit pores at 77 K. The kernel is obtained by combining grand canonical Monte Carlo (GCMC) method and open pore cell MC method that was developed in the previous study. We show that overall trends of the PSDs of activated carbons calculated with our new kernel and with conventional kernel from non-local density functional theory are nearly the same; however, apparent difference can be seen between them. As the third topic, we apply a free energy analysis method with the aid of GCMC simulations to investigate the gating behavior observed in a porous coordination polymer, and propose a mechanism for the adsorption-induced structural transition based on both the theory of equilibrium and kinetics. Finally, we construct an atomistic silica pore model that mimics MCM-41, which has atomic-level surface roughness, and perform molecular simulations to understand the mechanism of capillary condensation with hysteresis. We calculate the work required for the gas–liquid transition from the simulation data, and show that the adsorption branch with hysteresis for MCM-41 arise from spontaneous capillary condensation from a metastable state.     

Optimization of slitlike carbon nanopores for storage of hythane fuel at ambient temperatures

Carbons with slitlike pores can serve as effective host materials for storage of hythane fuel, a bridge between the petrol combustion and hydrogen fuel cells. We have used grand canonical Monte Carlo simulation for the modeling of the hydrogen and methane mixture storage at 293 K and pressure of methane and hydrogen mixture up to 2 MPa. We have found that these pores serve as efficient vessels for the storage of hythane fuel near ambient temperatures and low pressures. We find that, for carbons having optimized slitlike pores of size H congruent with 7 A (pore width that can accommodate one adsorbed methane layer), and bulk hydrogen mole fraction >or=0.9, the volumetric stored energy exceeds the 2010 target of 5.4 MJ dm(-3) established by the U.S. FreedomCAR Partnership. At the same condition, the content of hydrogen in slitlike carbon pores is approximately = 7% by energy. Thus, we have obtained the composition corresponding to hythane fuel in carbon nanospaces with greatly enhanced volumetric energy in comparison to the traditional compression method. We proposed the simple system with added extra container filled with pure free/adsorbed methane for adjusting the composition of the desorbed mixture as needed during delivery. Our simulation results indicate that light slit pore carbon nanomaterials with optimized parameters are suitable filling vessels for storage of hythane fuel. The proposed simple system consisting of main vessel with physisorbed hythane fuel, and an extra container filled with pure free/adsorbed methane will be particularly suitable for combustion of hythane fuel in buses and passenger cars near ambient temperatures and low pressures.     

Gas adsorption in active carbons and the slit-pore model 1: Pure gas adsorption

We describe procedures based on the polydisperse independent ideal slit-pore model, Monte Carlo simulation and density functional theory (a 'slab-DFT') for predicting gas adsorption and adsorption heats in active carbons. A novel feature of this work is the calibration of gas-surface interactions to a high surface area carbon, rather than to a low surface area carbon as in all previous work. Our models are used to predict the adsorption of carbon dioxide, methane, nitrogen, and hydrogen up to 50 bar in several active carbons at a range of near-ambient temperatures based on an analysis of a single 293 K carbon dioxide adsorption isotherm. The results demonstrate that these models are useful for relatively simple gases at near-critical or supercritical temperatures.     

Gibbs ensemble Monte Carlo simulation of supercritical CO2 adsorption on NaA and NaX zeolites

Adsorption of supercritical carbon dioxide on two kinds of zeolites with identical chemical composition but different pore structure (NaA and NaX) was studied using the Gibbs ensemble Monte Carlo simulation. The model frameworks for the two zeolites with SiAl ratio being unity have been chosen as the solid structures in the simulation. The adsorption behaviors of supercritical CO2 on the NaA and NaX zeolites, based on the adsorption isotherms and isosteric heats of adsorption, were discussed in detail and were compared with the available experimental results. A good agreement between the simulated and experimental results is obtained for both the adsorbed amount and the bulk phase density. The intermediate configurational snapshots and the radial distribution functions between zeolite and adsorbed CO2 molecules were collected in order to investigate the preferable adsorption locations and the confined structure behavior of CO2. The structure behaviors of the adsorbed CO2 molecules show various performances, as compared with the bulk phase, due to the confined effect in the zeolite pores.     

Adsorption equilibrium and kinetics of copper ions and phenol onto modified adsorbents

The adsorption equilibrium and kinetics of single and binary component copper ions and phenol onto powdered activated carbon (PAC), alginate beads and alginate-activated carbon beads (AAC) were studied. Adsorption equilibrium data for single component copper ions and phenol onto the adsorbents could be represented by the Langmuir equation. Multicomponent equilibrium data were correlated by the extended Langmuir and ideal adsorbed solution theory (IAST). The IAST gave the best fit to our data. The amount of copper ions adsorbed onto the AAC beads in the binary component was greater than that of phenol. The internal diffusion coefficients were determined by comparing the experimental concentration curves with those predicted from surface diffusion and pore diffusion model.     

Adsorption of benzene and propene in zeolite MCM-22: a grand canonical Monte Carlo study

The GCMC (grand canonical Monte Carlo) simulation technique was used to predict the competition adsorption characteristics of benzene and propene in different pore systems of MCM-22. The nine-site model of benzene was used, which proved to be effective and efficient. The zeolite was divided into three adsorption sites following a simulated annealing method. It is found that benzene and propene have the same preferential adsorption site and a similar adsorption order in different sites. Moreover, the pure and mixture isotherms of the three sites are drawn. From the isotherms, we obtained a selectivity reversal of the mixture isotherms of benzene and propene in different sites. It is also noted that the competition adsorption in the three adsorption sites for the two adsorbates can fall into three successive steps and the adsorption order of propene in mixture in these three sites is S3→S1→S2. A new model is presented to predict the benzene and propene adsorption equilibrium in MCM-22. This approach yields better multicomponent equilibrium predictions than ideal adsorbed solution theory (IAST). Isotherms at different mole fraction of benzene in gas phase indicate an advantage to increase the feed radio of benzene and propene. Thus, this work is helpful for a better understanding of the adsorption mechanism of benzene and propene in MCM-22 and hence the relation of the catalytic properties of the zeolite to its structure.     

Multicomponent Equilibrium Adsorption on Heterogeneous Adsorbents

A generalized method for prediction of multicomponent adsorption is suggested based on representing that adsorbent volume as energetically inhomogeneous. The method depends on extending the Polanyi potential theory to mixture adsorption. The main feature of the method is that, at constant partial pressure and temperature the composition of an adsorbed phase is not uniform over its volume. Results of applying this theory to non-porous adsorbents have been considered. The prediction ability of the theory is confirmed for the strongly non-ideal system acetone–chloroform–graphitized carbon black. It was shown that the departure from ideal behavior of adsorbed phase is quite close to that for the liquid mixture. Another system considered was oxygen–nitrogen–anatase at 78 K. Although this mixture is ideal, it has been found that there is significant variation in composition over the adsorbed layer due to the difference in the interactions of the quadrupolar N₂ molecule and nonpolar O₂ molecule with the anatase surface.