Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface

We present a new version of non-local density functional theory (NL-DFT) adapted to description of vapor adsorption isotherms on amorphous materials like non-porous silica. The novel feature of this approach is that it accounts for the roughness of adsorbent surface. The solid–fluid interaction is described in the same framework as in the case of fluid–fluid interactions, using the Weeks–Chandler–Andersen (WCA) scheme and the Carnahan–Starling (CS) equation for attractive and repulsive parts of the Helmholtz free energy, respectively. Application to nitrogen and argon adsorption isotherms on non-porous silica LiChrospher Si-1000 at their boiling points, recently published by Jaroniec and co-workers, has shown an excellent correlative ability of our approach over the complete range of pressures, which suggests that the surface roughness is mostly the reason for the observed behavior of adsorption isotherms. From the analysis of these data, we found that in the case of nitrogen adsorption short-range interactions between oxygen atoms on the silica surface and quadrupole of nitrogen molecules play an important role. The approach presented in this paper may be further used in quantitative analysis of adsorption and desorption isotherms in cylindrical pores such as MCM-41 and carbon nanotubes.

Improvement of the Derjaguin-Broekhoff-de Boer theory for capillary condensation/evaporation of nitrogen in mesoporous systems and its implications for pore size analysis of MCM-41 silicas and related materials

In this work, we propose an improvement of the classical Derjaguin-Broekhoff-de Boer (DBdB) theory for capillary condensation/evaporation in mesoporous systems. The primary idea of this improvement is to employ the Gibbs-Tolman-Koenig-Buff equation to predict the surface tension changes in mesopores. In addition, the statistical film thickness (so-called t-curve) evaluated accurately on the basis of the adsorption isotherms measured for the MCM-41 materials is used instead of the originally proposed t-curve (to take into account the excess of the chemical potential due to the surface forces). It is shown that the aforementioned modifications of the original DBdB theory have significant implications for the pore size analysis of mesoporous solids. To verify our improvement of the DBdB pore size analysis method (IDBdB), a series of the calcined MCM-41 samples, which are well-defined materials with hexagonally ordered cylindrical mesopores, were used for the evaluation of the pore size distributions. The correlation of the IDBdB method with the empirically calibrated Kruk-Jaroniec-Sayari (KJS) relationship is very good in the range of small mesopores. So, a major advantage of the IDBdB method is its applicability for small mesopores as well as for the mesopore range beyond that established by the KJS calibration, i.e., for mesopore radii greater than similar to4.5 nm. The comparison of the IDBdB results with experimental data reported by Kruk and Jaroniec for capillary condensation/evaporation as well as with the results from nonlocal density functional theory developed by Neimark et al. clearly justifies our approach. Note that the proposed improvement of the classical DBdB method preserves its original simplicity and simultaneously ensures a significant improvement of the pore size analysis, which is confirmed by the independent estimation of the mean pore size by the powder X-ray diffraction method.

Equilibrium adsorption in cylindrical mesopores: A modified Broekhoff and de boer theory versus density functional theory

This paper presents a thermodynamic analysis of capillary condensation phenomena in cylindrical pores. Here, we modified the Broekhoff and de Boer (BdB) model for cylindrical pores accounting for the effect of the pore radius on the potential exerted by the pore walls. The new approach incorporates the recently published standard nitrogen and argon adsorption isotherm on nonporous silica LiChrospher Si-1000. The developed model is tested against the nonlocal density functional theory (NLDFT), and the criterion for this comparison is the condensation/evaporation pressure versus the pore diameter. The quantitative agreement between the NLDFT and the refined version of the BdB theory is ascertained for pores larger than 2 nm. The modified BdB theory was applied to the experimental adsorption branch of adsorption isotherms of a number of MCM-41 samples to determine their pore size distributions (PSDs). It was found that the PSDs determined with the new BdB approach coincide with those determined with the NLDFT (also using the experimental adsorption branch). As opposed to the NLDFT, the modified BdB theory is very simple in its utilization and therefore can be used as a convenient tool to obtain PSDs of all mesoporous solids from the analysis of the adsorption branch of adsorption isotherms of any subcritical fluids.

Modeling of adsorption in finite cylindrical pores by means of density functional theory

Adsorption of argon at its boiling point infinite 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.

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.

Comparison of nitrogen adsorption at 77 K on non-porous silica and pore wall of MCM-41 materials by means of density functional theory

Nitrogen adsorption on a surface of a non-porous reference material is widely used in the characterization. Traditionally, the enhancement of solid-fluid potential in a porous solid is accounted for by incorporating the surface curvature into the solid-fluid Potential of the flat reference surface. However, this calculation procedure has not been justified experimentally. In this paper, we derive the solid-fluid potential of mesoporous MCM-41 solid by using solely the adsorption isotherm of that solid. This solid-fluid potential is then compared with that of the non-porous reference surface. In derivation of the solid-fluid potential for both reference surface and mesoporous MCM-41 silica (diameter ranging front 3 to 6.5 nm) we employ the nonlocal density functional theory developed for amorphous solids. It is found that, to out, surprise, the solid-fluid potential of a porous solid is practically the same as that for the reference surface, indicating that there is no enhancement due to Surface curvature. This requires further investigations to explain this unusual departure from our conventional wisdom of curvature-induced enhancement. Accepting the curvature-independent solid-fluid potential derived from the non-porous reference surface, we analyze the hysteresis features of a series of MCM-41 samples. (c) 2005 Elsevier Inc. All rights reserved.

Modeling hydrogen separation in high temperature silica membrane systems

In this work, a working model is proposed of molecular sieve silica (MSS) multistage membrane systems for CO cleanup at high temperatures (up to 500 degrees C) in a simulated fuel cell fuel processing system. Gases are described as having little interactions with each other relative to the pore walls due to low isosteric heat of adsorption on silica surfaces and high temperatures. The Arrhenius function for activated transport of pure gases was used to predict mixture concentration in the permeate and retentate streams. Simulation predicted CO could be reduced to levels below the required 50 ppmv for polymer electrolyte membrane fuel cell anodes at a stage H-2/CO selectivity of higher than 40 in 4 series membrane units. Experimental validation showed predicting mixture concentrations required only pure gas permeation data. This model has significant application for setting industrial stretch targets and as a robust basis for complex membrane model configurations. (c) 2006 American Institute of Chemical Engineers.

NMR microscopy: A tool for measuring monomer diffusion in supercritical CO2

The diffusion of styrene into linear low density polyethylene in a solution of supercritical CO2 was investigated using NMR microimaging. For both pure styrene and styrene dissolved in supercritical CO2, the diffusion was found to follow Fickian kinetics. Supercritical CO2 was found to enhance the rate and extent of diffusion of styrene into the substrate by up to three times under the conditions of this investigation, compared to pure styrene. NMR imaging was used to measure the concentration profiles of the styrene penetrants in real time, and the results were fitted to a Fickian model for diffusion. At a CO2 pressure of 150 bar and temperature of 40 degrees C, the diffusion coefficient of a 30 wt-% solution of styrene into LLDPE was calculated to be 1 X 10(-11) m(2). s(-1). This is significantly faster than the diffusion coefficient measured for pure styrene diffusion at 40 degrees C (3 x 10(-12) m(2). s(-1)). The diffusion coefficients determined by gravimetric analysis were slightly higher than those determined by the imaging method. This was probably due to residual styrene and/or polystyrene adhering to the surface of the films in the gravimetric technique.

Thermal Cycling Stability of Silica Membranes for Gas Separation

Hydrogen is being seen as an alternative energy carrier to conventional hydrocarbons to reduce greenhouse gas emissions. High efficiency separation technologies to remove hydrogen from the greenhouse gas, carbon dioxide, are therefore in growing demand. Traditional thermodynamic separation systems utilise distillation, absorption and adsorption, but are limited in efficiency at compact scales. Molecular sieve silica (MSS) membranes can perform this separation as they have high permselectivity of hydrogen to carbon dioxide, but their stability under thermal cycling is not well reported. In this work we exposed a standard MSS membrane and a carbonised template MSS (CTMSS) membrane to thermal cycling from 100 to 450°C. The standard MSS and carbonised template CTMSS membranes both showed permselectivity of helium to nitrogen dropping from around 10 to 6 in the first set of cycles, remaining stable until the last test. The permselectivity drop was due to small micropore collapse, which occurred via structure movement during cycling. Simulating single stage membrane separation with a 50:50 molar feed of H2:CO2, H2 exiting the permeate stream would start at 79% and stabilise at 67%. Higher selectivity membranes showed less of a purity drop, indicating the margin at which to design a stable membrane separation unit for CO2 capture.

Functionalisation of nanoparticles as electrolytes in fuel cells

Inorganic metal oxide materials are generally poor proton conductors as conductivities are lower than 10-5-10-6 S.cm-1. However, by functionalising Silica, Zirconia or Titania, proton conduction increases by up to 5 orders of magnitude. Hence, functionalised nanomaterials are becoming very competitive against conventional electrolyte materials such as Nafion. In this work, sol-gel processes are employed to produce silica phosphate, zirconia phosphate and titania phosphate functionalised nanoparticles. Furthermore, conductivities at hydrate conditions are investigated, and nanoparticle formation and functionalisation effects on proton conductivity are discussed. Results show conductivities up to 10-1 S.cm-1 (95% RH). Proton conduction increases with the functionalisation content, however heat treatment of nanoparticles locks the functionality in the crystal phase, thus inhibiting proton conduction. Controlling the mesopore phase allows for high proton conduction at hydrated conditions, clearly indicating facilitated ion transport through the pore channels.

Hydrothermal stable templated molecular sieve silica (TMSS) membranes for gas separation

In this work we compare the hydrothermal stability performance of a Templated Molecular Sieve Silica (TMSS) membrane against a standard, non-templated Molecular Sieve Silica (MSS) membrane. The tests were carried under dry and wet (steam) conditions for single gas (He, H2, CO and CO2) at 1-2 atm membrane pressure drop at 200oC. Single gas TMSS membrane H2, permeance and H2/CO permselectivity was found to be 2.05 x 10-8 mols.m-2.s-1.Pa-1 and 15, respectively. The MSS membrane showed similar selectivity, but increased overall flux. He permeance through membranes decayed at a rate of 4-5 x 10-10 mols.m-2.s-1.Pa-1 per day regardless of membrane ambience (dry or wet). Although H2/CO permselectivity of the TMSS membrane slightly improved from 15 to 18 after steam testing, the MSS membrane resulted in significant reduction from 16 to 8.3. In addition, membrane regeneration after more than 50 days resulted in the TMSS membrane reverting to its original permeation levels while no significant improvements were observed for the MSS membra ne. Results showed that the TMSS membrane had enhanced hydrothermal stability and regeneration ability.

Hydrothermal Stability of Modified Silica Membranes for Gas Separation

A new class of hybrid molecular sieve silica (MSS) membranes is developed and tested against standard and organic templated membranes. The hybrid membrane is synthesized by the standard sol-gel process, integrating a template (methyltriethoxysilane - MTES) and a C6 surfactant (triethylhexylammonium bromide) into the silica film matrix. After hydro treatment under a relative humidity of 96% for 50h, the hybrid membrane shows no changes in its gas separation capabilities or energy of mobility. The structural characteristics and integrity of the hybrid membrane are retained due to a high concentration of organophilic functional groups and alkoxides observed using 29 Si NMR. In contrast, the structural integrity of the membranes prepared with non-templated films deteriorated during the hydro treatment due to a large percentage of silanol groups (Si-OH) which react with water. The hybrid membranes underwent a decrease in the H2/CO2 selectivity of only 1% whereas for the non-templated membrane a 21% decrease was observed. The transport mechanism of the hybrid membranes is activated as permeation increased with temperature. The activation energy for the permeation of H2 is positive while negative for CO2. The H2 permeation obtained was 3x 10 -8 mol.m -2 .s -1 .Pa -1 and permselectivities for H2/CO2 and H2/N2 varied between 1-7 and 31-34, respectively.

Hydrostability and Scaling Up Molecular Sieve Silica (MSS) Membranes for H2/CO Separation in Fuel Cell Systems

MSS membranes are a good candidate for CO cleanup in fuel cell fuel processing systems due to their ability to selectively permeate H2 over CO via molecular sieving. Successfully scaled up tubular membranes were stable under dry conditions to 400°C with H2 permeance as high as 2 x 10-6 mol.m-2.s^-1.Pa^-1 at 200 degrees C and H2/CO selectivity up to 6.4, indicating molecular sieving was the dominant mechanism. A novel carbonised template molecular sieve silica (CTMSS) technology gave the scaled up membranes resilience in hydrothermal conditions up to 400 degrees C in 34% steam and synthetic reformate, which is required for use in fuel cell CO cleanup systems.

Oxidation mechanism of Si3N4-bonded SiC ceramics by CO, CO2 and steam

This paper presents a theoretical and experimental investigation into the oxidation reactions of Si3N4-bonded SiC ceramics. Such ceramics which contain a small amount of silicon offer increased oxidation and wear resistance and are widely used as lining refractories in blast furnaces. The thermodynamics of oxidation reactions were studied using the JANAF tables. The weight gain was measured using a thermogravimetric analysis technique to study the kinetics. The temperature range of oxidation measurements is from 1073 to 1573 K and the oxidation atmosphere is water vapour, pure CO and CO-CO2 gas mixtures with various CO-to-CO2 ratios. Thermodynamic simulations showed that the oxidation mechanism of Si3N4-bonded SiC ceramics is passive oxidation and all components contribute to the formation of a silica film. The activated energies of the reactions follow the sequence Si3N4>SiC>Si. The kinetic study revealed that the oxidation of Si3N4-bonded SiC ceramics occurred in a mixed regime controlled by both interface reaction and diffusion through the silica film. Under the atmosphere conditions prevailing in the blast furnace, this ceramic is predicted to be passively oxidized with the chemical reaction rate becoming more dominant as the CO concentration increases. (C) 1998 Chapman & Hall.

Capillary network models for transport in packed beds: Considerations of pore aspect ratio

The conventional analysis for the estimation of the tortuosity factor for transport in porous media is modified here to account for the effect of pore aspect ratio. Structural models of the porous medium are also constructed for calculating the aspect ratio as a function of porosity. Comparison of the model predictions with the extensive data of Currie (1960) for the effective diffusivity of hydrogen in packed beds shows good agreement with a network model of randomly oriented intersecting pores for porosities upto about 50 percent, which is the region of practical interest. The predictions based on this network model are also found to be in better agreement with the data of Currie than earlier expressions developed for unconsolidated and grainy media.