Prediction of high-pressure adsorption equilibrium of supercritical gases using density functional theory

In this paper, we present the results of the prediction of the high-pressure adsorption equilibrium of supercritical. gases (Ar, N-2, CH4, and CO2) on various activated carbons (BPL, PCB, and Norit R1 extra) at various temperatures using a density-functional-theory-based finite wall thickness (FWT) model. Pore size distribution results of the carbons are taken from our recent previous work 1,2 using this approach for characterization. To validate the model, isotherms calculated from the density functional theory (DFT) approach are comprehensively verified against those determined by grand canonical Monte Carlo (GCMC) simulation, before the theoretical adsorption isotherms of these investigated carbons calculated by the model are compared with the experimental adsorption measurements of the carbons. We illustrate the accuracy and consistency of the FWT model for the prediction of adsorption isotherms of the all investigated gases. The pore network connectivity problem occurring in the examined carbons is also discussed, and on the basis of the success of the predictions assuming a similar pore size distribution for accessible and inaccessible regions, it is suggested that this is largely related to the disordered nature of the carbon.

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.

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.

Capillary phenomena in the framework of the two-dimensional density functional theory

We present results of application of the density functional theory (DFT) to adsorption and desorption in finite and infinite cylindrical pores accounting for the density distribution in radial and axial directions. Capillary condensation via formation of bridges is considered using canonical and grand canonical versions of the 2D DFT. The potential barrier of nucleation is determined as a function of the bulk pressure and the pore diameter. In the framework of the conventional assumptions on intermolecular interactions both 1D and 2D DFT versions lead to the same results and confirm the classical scenario of condensation and evaporation: the condensation occurs at the vapor-like spinodal point, and the evaporation corresponds to the equilibrium transition pressure. The analysis of experimental data on argon and nitrogen adsorption on MCM-41 samples seems to not completely corroborate this scenario, with adsorption branch being better described by the equilibrium pressure - diameter dependence. This points to the necessity of the further development of basic representations on the hysteresis phenomena.

Van der Waals-corrected density functional theory: benchmarking for hydrogen-nanotube and nanotube-nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H-2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H-2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

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.

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.

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.

Prediction of the equilibrium structures and photomagnetic properties of the Prussian blue analogue RbMn[Fe(CN)6] by density functional theory

A periodic density functional theory method using the B3LYP hybrid exchange-correlation potential is applied to the Prussian blue analogue RbMn[Fe(CN)6] to evaluate the suitability of the method for studying, and predicting, the photomagnetic behavior of Prussian blue analogues and related materials. The method allows correct description of the equilibrium structures of the different electronic configurations with regard to the cell parameters and bond distances. In agreement with the experimental data, the calculations have shown that the low-temperature phase (LT; Fe(2+)(t(6)2g, S = 0)-CN-Mn(3+)(t(3)2g e(1)g, S = 2)) is the stable phase at low temperature instead of the high-temperature phase (HT; Fe(3+)(t(5)2g, S = 1/2)-CN-Mn(2+)(t(3)2g e(2)g, S = 5/2)). Additionally, the method gives an estimation for the enthalpy difference (HT LT) with a value of 143 J mol(-1) K(-1). The comparison of our calculations with experimental data from the literature and from our calorimetric and X-ray photoelectron spectroscopy measurements on the Rb0.97Mn[Fe(CN)6]0.98 x 1.03 H2O compound is analyzed, and in general, a satisfactory agreement is obtained. The method also predicts the metastable nature of the electronic configuration of the high-temperature phase, a necessary condition to photoinduce that phase at low temperatures. It gives a photoactivation energy of 2.36 eV, which is in agreement with photoinduced demagnetization produced by a green laser.

Metastable de-excitation spectroscopy and density functional theory study of the selective oxidation of crotyl alcohol over Pd(111)

The extremely surface sensitive technique of metastable de-excitation spectroscopy (MDS) has been utilized to probe the bonding and reactivity of crotyl alcohol over Pd(111) and provide insight into the selective oxidation pathway to crotonaldehyde. Auger de-excitation (AD) of metastable He (23S) atoms reveals distinct features associated with the molecular orbitals of the adsorbed alcohol, corresponding to emission from the hydrocarbon skeleton, the O n nonbonding, and C═C π states. The O n and C═C π states of the alcohol are reversed when compared to those of the aldehyde. Density functional theory (DFT) calculations of the alcohol show that an adsorption mode with both C═C and O bonds aligned somewhat parallel to the surface is energetically favored at a substrate temperature below 200 K. Density of states calculations for such configurations are in excellent agreement with experimental MDS measurements. MDS revealed oxidative dehydrogenation of crotyl alcohol to crotonaldehyde between 200 and 250 K, resulting in small peak shifts to higher binding energy. Intramolecular changes lead to the opposite assignment of the first two MOs in the alcohol versus the aldehyde, in accordance with DFT and UPS studies of the free molecules. Subsequent crotonaldehyde decarbonylation and associated propylidyne formation above 260 K could also be identified by MDS and complementary theoretical calculations as the origin of deactivation and selectivity loss. Combining MDS and DFT in this way represents a novel approach to elucidating surface catalyzed reaction pathways associated with a “real-world” practical chemical transformation, namely the selective oxidation of alcohols to aldehydes.

Hybrid density functional theory description of N- and C-doping of NiO

The large intrinsic bandgap of NiO hinders its potential application as a photocatalyst under visible-light irradiation. In this study, we have performed first-principles screened exchange hybrid density functional theory with the HSE06 functional calculations of N- and C-doped NiO to investigate the effect of doping on the electronic structure of NiO. C-doping at an oxygen site induces gap states due to the dopant, the positions of which suggest that the top of the valence band is made up primarily of C 2p-derived states with some Ni 3d contributions, and the lowest-energy empty state is in the middle of the gap. This leads to an effective bandgap of 1.7 eV, which is of potential interest for photocatalytic applications. N-doping induces comparatively little dopant-Ni 3d interactions, but results in similar positions of dopant-induced states, i.e., the top of the valence band is made up of dopant 2p states and the lowest unoccupied state is the empty gap state derived from the dopant, leading to bandgap narrowing. With the hybrid density functional theory (DFT) results available, we discuss issues with the DFT corrected for on-site Coulomb description of these systems.