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.

Testing multiple mapping conditioning mixing for Monte Carlo probability density function simulations

Mitarai [Phys. Fluids 17, 047101 (2005)] compared turbulent combustion models against homogeneous direct numerical simulations with extinction/recognition phenomena. The recently suggested multiple mapping conditioning (MMC) was not considered and is simulated here for the same case with favorable results. Implementation issues crucial for successful MMC simulations are also discussed.

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.

Pore size distribution analysis of activated carbons: Application of density functional theory using nongraphitized carbon black as a reference system

The application of nonlocal density functional theory (NLDFT) to determine pore size distribution (PSD) of activated carbons using a nongraphitized carbon black, instead of graphitized thermal carbon black, as a reference system is explored. We show that in this case nitrogen and argon adsorption isotherms in activated carbons are precisely correlated by the theory, and such an excellent correlation would never be possible if the pore wall surface was assumed to be identical to that of graphitized carbon black. It suggests that pore wall surfaces of activated carbon are closer to that of amorphous solids because of defects of crystalline lattice, finite pore length, and the presence of active centers.. etc. Application of the NLDFT adapted to amorphous solids resulted in quantitative description of N-2 and Ar adsorption isotherms on nongraphitized carbon black BP280 at their respective boiling points. In the present paper we determined solid-fluid potentials from experimental adsorption isotherms on nongraphitized carbon black and subsequently used those potentials to model adsorption in slit pores and generate a corresponding set of local isotherms, which we used to determine the PSD functions of different activated carbons. (c) 2005 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 Comparison of Bayesian and Sampling Theory Inferences in a Probit Model

HE PROBIT MODEL IS A POPULAR DEVICE for explaining binary choice decisions in econometrics. It has been used to describe choices such as labor force participation, travel mode, home ownership, and type of education. These and many more examples can be found in papers by Amemiya (1981) and Maddala (1983). Given the contribution of economics towards explaining such choices, and given the nature of data that are collected, prior information on the relationship between a choice probability and several explanatory variables frequently exists. Bayesian inference is a convenient vehicle for including such prior information. Given the increasing popularity of Bayesian inference it is useful to ask whether inferences from a probit model are sensitive to a choice between Bayesian and sampling theory techniques. Of interest is the sensitivity of inference on coefficients, probabilities, and elasticities. We consider these issues in a model designed to explain choice between fixed and variable interest rate mortgages. Two Bayesian priors are employed: a uniform prior on the coefficients, designed to be noninformative for the coefficients, and an inequality restricted prior on the signs of the coefficients. We often know, a priori, whether increasing the value of a particular explanatory variable will have a positive or negative effect on a choice probability. This knowledge can be captured by using a prior probability density function (pdf) that is truncated to be positive or negative. Thus, three sets of results are compared:those from maximum likelihood (ML) estimation, those from Bayesian estimation with an unrestricted uniform prior on the coefficients, and those from Bayesian estimation with a uniform prior truncated to accommodate inequality restrictions on the coefficients.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Analysis of adsorption data of graphitized thermal carbon black with a DFT-lattice gas theory

In this paper we analyzed the adsorption of gases and vapors on graphitised thermal carbon black by using a modified DFT-lattice theory, in which we assume that the behavior of the first layer in the adsorption film is different from those of second and higher layers. The effects of various parameters on the topology of the adsorption isotherm were first investigated, and the model was then applied in the analysis of adsorption data of numerous substances on carbon black. We have found that the first layer in the adsorption film behaves differently from the second and higher layers in such a way that the adsorbate-adsorbate interaction energy in the first layer is less than that of second and higher layers, and the same is observed for the partition function. Furthermore, the adsorbate-adsorbate and adsorbate-adsorbent interaction energies obtained from the fitting are consistently lower than the corresponding values obtained from the viscosity data and calculated from the Lorentz-Berthelot rule, respectively.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Semiquantum versus semiclassical mechanics for simple nonlinear systems

Quantum mechanics has been formulated in phase space, with the Wigner function as the representative of the quantum density operator, and classical mechanics has been formulated in Hilbert space, with the Groenewold operator as the representative of the classical Liouville density function. Semiclassical approximations to the quantum evolution of the Wigner function have been defined, enabling the quantum evolution to be approached from a classical starting point. Now analogous semiquantum approximations to the classical evolution of the Groenewold operator are defined, enabling the classical evolution to be approached from a quantum starting point. Simple nonlinear systems with one degree of freedom are considered, whose Hamiltonians are polynomials in the Hamiltonian of the simple harmonic oscillator. The behavior of expectation values of simple observables and of eigenvalues of the Groenewold operator are calculated numerically and compared for the various semiclassical and semiquantum approximations.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

Facile ketene-ketene and ketene-ketenimine rearrangements: A study of the 1,3-migration of alpha-substituents interconverting alpha-imidoylketenes and alpha-oxoketenimines, a pseudopericyclic reaction

Theoretical calculations (B3LYP/6-311+G(3df,2p)//B3LYP/6-31G*) of the 1,3 migration of NR2 transforming alpha-oxoketenimines 1 to alpha-imidoylketenes 3 and vice versa indicate that this process is a pseudo-pericyclic reaction with a low activation energy (NH2 97 kJ mol(-1), N(CH3)(2) 62 kJ mol(-1)). The oxoketenimines were found to be more stable (by 18-35 kJ mol(-1)) which is in line with experimental observations. The hindered amine rotation in the amide and amidine moieties adjacent to the cumulenes are important in the migration of the NR2 group, as one of the rotation transition states is close to the 1,3 migration pathway. This gives an interesting potential energy surface with a valley-ridge inflection (VRI) between the orthogonal hindered amine rotation and 1,3 migration transition states. The imidoylketene may also undergo ring closure to an azetinone 5; however, this is metastable, and under the conditions that allow the 1,3-migration, the oxoketenimine 1 will be favored. The imine NH E/Z-interconversion of the ketenimine group takes place by inversion and has a low activation barrier (similar to40 kJ mol(-1)). In all the amidines examined the E/Z-interconversion of the imine function was predicted to be by rotation with a high barrier (>80 kJ mol(-1)), in contrast to all other reported imine E/Z-interconversions which are by inversion.

High-pressure adsorption capacity and structure of CO2 in carbon slit pores: Theory and simulation

We present new simulation results for the packing of single-center and three-center models of carbon dioxide at high pressure in carbon slit pores. The former shows a series of packing transitions that are well described by our density functional theory model developed earlier. In contrast, these transitions are absent for the three-center model. Analysis of the simulation results shows that alternations of flat-lying molecules and rotated molecules can occur as the pore width is increased. The presence or absence of quadrupoles has negligible effect on these high-density structures.