Application of density functional theory to capillary phenomena in cylindrical mesopores with radial and longitudinal density distributions

In this paper, we applied a version of the nonlocal density functional theory (NLDFT) accounting radial and longitudinal density distributions to study the adsorption and desorption of argon in finite as well as infinite cylindrical nanopores at 87.3 K. Features that have not been observed before with one-dimensional NLDFT are observed in the analysis of an inhomogeneous fluid along the axis of a finite cylindrical pore using the two-dimensional version of the NLDFT. The phase transition in pore is not strictly vapor-liquid transition as assumed and observed in the conventional version, but rather it exhibits a much elaborated feature with phase transition being complicated by the formation of solid phase. Depending on the pore size, there are more than one phase transition in the adsorption-desorption isotherm. The solid formation in finite pore has been found to be initiated by the presence of the meniscus. Details of the analysis of the extended version of NLDFT will be discussed in the paper. (C) 2004 American Institute of Physics.

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.

Study of hexane adsorption in nanoporous MCM-41 silica

We study here the adsorption of hexane on nanoporous MCM-41 silica at 303, 313, and 323 K, for various pore diameters between 2.40 and 4.24 nm. Adsorption equilibria, measured thermogravimetrically, show that all the isotherms, that are somewhat akin to those of type V, exhibit remarkably sharp capillary adsorption phase transition steps and are reversible. The position of the phase transition step gradually shifts from low to high relative pressure with an increase in the temperature as well as the pore sizes. The isosteric heats of adsorption derived from the equilibrium information using the Clapeyron equation reveal a gradual decrease with increasing adsorbed amount because of the surface heterogeneity but approach a constant value near the phase transition. A decrease in the pore size results in an increase in the isosteric heat of adsorption because of the increased dispersion forces. A simple strategy, based on the Broekhoff and De Boer adsorption theory, successfully interprets the hexane adsorption isotherms for the different pore size MCM-41 samples. The parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting the monolayer region prior to capillary condensation and the experimental phase transition simultaneously, for some pore sizes. Subsequently, the parameters are used to predict the adsorption isotherm on other pore size samples, which showed good agreement with experimental data.

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.

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.

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.

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory

Adsorption of argon and nitrogen at their respective boiling points in cylindrical pores of MCM-41 type silica-like adsorbents is studied by means of a non-local density functional theory (NLDFT), which is modified to deal with amorphous solids. By matching the theoretical results of the pore filling pressure versus pore diameter against the experimental data, we arrive at a conclusion that the adsorption branch (rather than desorption) corresponds to the true thermodynamic equilibrium. If this is accepted, we derive the optimal values for the solid–fluid molecular parameters for the system amorphous silica–Ar and amorphous silica–N2, and at the same time we could derive reliably the specific surface area of non-porous and mesoporous silica-like adsorbents, without a recourse to the BET method. This method is then logically extended to describe the local adsorption isotherms of argon and nitrogen in silica-like pores, which are then used as the bases (kernel) to determine the pore size distribution. We test this with a number of adsorption isotherms on the MCM-41 samples, and the results are quite realistic and in excellent agreement with the XRD results, justifying the approach adopted in this paper.

Batchelor's theory extended to elongated cylindrical or ellipsoidal particles

Extensions to Batchelor's theory have been derived to take into account different shaped particles while relating extensional viscosity enhancement to three parameters - shape, volume fraction and particle aspect ratio. The extended theory now allows calculation of the extensional viscosity enhancement, at a given volume fraction of particles, for either ellipsoidal or cylindrical particles. The formula improves the predictive capability of Batchelor's theory when compared with measurements found in the literature for different rod-like polymer solutions. (c) 2005 Elsevier B.V. All rights reserved.

Modeling nitrogen adsorption in spherical pores of siliceous materials by density functional theory

Adsorption of nitrogen in spherical pores of FDU-1 silica at 77 K is considered by means of a nonlocal density functional theory (NLDFT) accounting for a disordered structure of pore walls. Pore size distribution analysis of various FDU-1 samples subject to different temperatures of calcination revealed three distinct groups of pores. The principal group of pores is identified as ordered spherical mesopores connected with each other by smaller interconnecting pores and irregular micropores present in the mesopore walls. To account for the entrances (connecting pores) into spherical mesopores, a concept of solid mass distribution with respect to the apparent density was introduced. It is shown that the introduction of the aforementioned distribution was sufficient to quantitatively describe experimental adsorption isotherms over the entire range of relative pressures spanning six decades.

Modeling of adsorption and nucleation in infinite cylindrical pores by two-dimensional density functional theory

In this paper, we present an analysis of argon adsorption in cylindrical pores having amorphous silica structure by means of a nonlocal density functional theory (NLDFT). In the modeling, we account for the radial and longitudinal density distributions, which allow us to consider the interface between the liquidlike and vaporlike fluids separated by a hemispherical meniscus in the canonical ensemble. The Helmholtz free energy of the meniscus was determined as a function of pore diameter. The canonical NLDFT simulations show the details of density rearrangement at the vaporlike and liquidlike spinodal points. The limits of stability of the smallest bridge and the smallest bubble were also determined with the canonical NLDFT. The energy of nucleation as a function of the bulk pressure and the pore diameter was determined with the grand canonical NLDFT using an additional external potential field. It was shown that the experimentally observed reversibility of argon adsorption isotherms at its boiling point up to the pore diameter of 4 nm is possible if the potential barrier of 22kT is overcome due to density fluctuations.

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.

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.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.