Pore size characterization of mesoporous materials by a thermodynamic approach: A curvature-dependent solid-fluid potential

A thermodynamic analysis of nitrogen adsorption in cylindrical pores of MCM-41 and SBA-15 samples at 77 K is presented within the framework of the Broekhoff and de Boer (BdB) theory. We accounted for the effect of the solid surface curvature on the potential exerted by the pore walls. The developed model is in quantitative agreement with the non-local density functional theory (NLDFT) for pores larger than 2 tun. This modified BdB theory accounting for the Curvature Dependent Potential (CDP-BdB) was applied to determine the pore size distribution (PSD) of a number of MCM-41 and SBA-15 samples on the basis of matching the equilibrium theoretical isotherm against the adsorption branch of the experimental isotherm. In all cases investigated the PSDs determined with the new approach are very similar to those determined with the non-local density functional theory also using the same basis of matching of theoretical isotherm against the experimental adsorption branch. The developed continuum theory is very simple in its utilization, suggesting that CDP-BdB could be used as an alternative tool to obtain PSD for mesoporous solids from the analysis of adsorption branch of adsorption isotherms of any sub-critical fluids.

A SPH solver for simulating paramagnetic solid fluid interaction in the presence of an external magnetic field

Acknowledgment The first two authors wish to express their sincerest thanks to Iran National Science Foundation (INSF) for supporting this work under Contract Number 92021291.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Multicomponent adsorption on activated carbons under supercritical conditions

Adsorption of binary mixtures onto activated carbon Norit R1 for the system nitrogen-methane-carbon dioxide was investigated over the pressure range up to 15 MPa. A new model is proposed to describe the experimental data. It is based on the assumption that an activated carbon can be characterized by the distribution function of elements of adsorption volume (EAV) over the solid-fluid potential. This function may be evaluated from pure component isotherms using the equality of the chemical potentials in the adsorbed phase and in the bulk phase for each EAV. In the case of mixture adsorption a simple combining rule is proposed, which allows determining the adsorbed phase density and its composition in the EAV at given pressure and compositions of the bulk phase. The adsorbed concentration of each adsorbate is the integral of its density over the set of EAV. The comparison with experimental data on binary mixtures has shown that the approach works reasonably well. In the case of high-pressure binary mixture adsorption, when only total amount adsorbed was measured, the proposed model allows reliably determining partial amounts of the adsorbed components. (C) 2004 Elsevier Inc. All rights reserved.

Application of a generalized thermodynamic approach to characterize mesoporous materials

Equilibrium adsorption and desorption in mesoporous adsorbents is considered on the basis of rigorous thermodynamic analysis, in which the curvature-dependent solid-fluid potential and the compressibility of the adsorbed phase are accounted for. The compressibility of the adsorbed phase is considered for the first time in the literature in the framework of a rigorous thermodynamic approach. Our model is a further development of continuum thermodynamic approaches proposed by Derjaguin and Broekhoff and de Boer, and it is based on a reference isotherm of a non-porous material having the same chemical structure as that of the pore wall. In this improved thermodynamic model, we incorporated a prescription for transforming the solid-fluid potential exerted by the flat reference surface to the potential inside cylindrical and spherical pores. We relax the assumption that the adsorbed film density is constant and equal to that of the saturated liquid. Instead, the density of the adsorbed fluid is allowed to vary over the adsorbed film thickness and is calculated by an equation of state. As a result, the model is capable to describe the adsorption-desorption reversibility in cylindrical pores having diameter less than 2 nm. The generalized thermodynamic model may be applied to the pore size characterization of mesoporous materials instead of much more time-consuming molecular approaches. (c) 2005 Elsevier B.V. All rights reserved.

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.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.

Effects of potential models on the adsorption of ethane and ethylene on graphitized thermal carbon black. Study of two-dimensional critical temperature and isosteric heat versus loading

Adsorption of ethylene and ethane on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers is studied in detail to investigate the packing efficiency, the two-dimensional critical temperature, and the variation of the isosteric heat of adsorption with loading and temperature. Here we used a Monte Carlo simulation method with a grand canonical Monte Carlo ensemble. A number of two-center Lennard-Jones (LJ) potential models are investigated to study the impact of the choice of potential models in the description of adsorption behavior. We chose two 2C-LJ potential models in our investigation of the (i) UA-TraPPE-LJ model of Martin and Siepmann (J. Phys. Chem. B 1998,102, 25692577) for ethane and Wick et al. (J. Phys. Chem. B 2000,104, 8008-8016) for ethylene and (ii) AUA4-LJ model of Ungerer et al. (J. Chem. Phys. 2000,112, 5499-5510) for ethane and Bourasseau et al. (J. Chem. Phys. 2003, 118, 3020-3034) for ethylene. These models are used to study the adsorption of ethane and ethylene on graphitized thermal carbon black. It is found that the solid-fluid binary interaction parameter is a function of adsorbate and temperature, and the adsorption isotherms and heat of adsorption are well described by both the UA-TraPPE and AUA models, although the UA-TraPPE model performs slightly better. However, the local distributions predicted by these two models are slightly different. These two models are used to explore the two-dimensional condensation for the graphitized thermal carbon black, and these values are 110 K for ethylene and 120 K for ethane.

Impact of potential models on adsorption of linear molecules on carbon black

In this paper, we investigate the effects of potential models on the description of equilibria of linear molecules (ethylene and ethane) adsorption on graphitized thermal carbon black. GCMC simulation is used as a tool to give adsorption isotherms, isosteric heat of adsorption and the microscopic configurations of these molecules. At the heart of the GCMC are the potential models, describing fluid-fluid interaction and solid-fluid interaction. Here we studied the two potential models recently proposed in the literature, the UA-TraPPE and AUA4. Their impact in the description of adsorption behavior of pure components will be discussed. Mixtures of these components with nitrogen and argon are also studied. Nitrogen is modeled a two-site plus discrete charges while argon as a spherical particle. GCMC simulation is also used for generating simulation mixture isotherms. It is found that co-operation between species occurs when the surface is fractionally covered while competition is important when surface is fully loaded.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

Capillary coexistence and criticality in mesopores: Modification of the Kelvin theory

The classical model of capillary equilibrium in cylindrical pores is modified here by the introduction of molecular concepts and the solid fluid interaction potential. The new approach accurately predicts capillary coexistence and criticality, with results quantitatively matching those from density functional theory for nitrogen adsorption, while also predicting condensation pressures in agreement with reported experimental findings for MCM-41. The larger critical pore size for nitrogen adsorption in these materials, however, suggests a modification of the potential function parameters, evaluated here from data for hydroxylated silica.

A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.