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.     

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.     

Lithium adsorption on graphite from density functional theory calculations

The structural, energetic, and electronic properties of the Li/graphite system are studied through density functional theory (DFT) calculations using both the local spin density approximation (LSDA), and the gradient-corrected Perdew-Burke-Ernzerhof (PBE) approximation to the exchange-correlation energy. The calculations were performed using plane waves basis, and the electron-core interactions are described using pseudopotentials. We consider a disperse phase of the adsorbate comprising one Li atom for each 16 graphite surface cells, in a slab geometry. The close contact between the Li nucleus and the graphene plane results in a relatively large binding energy (larger than 1.1 eV). A detailed analysis of the electronic charge distribution, density difference distribution, and band structures indicates that one valence electron is entirely transferred from the atom to the surface, which gives rise to a strong interaction between the resulting lithium ion and the cloud of pi electrons in the substrate. We show that it is possible to explain the differences in the binding of Li, Na, and K adatoms on graphite considering the properties of the corresponding cation/aromatic complexes.     

NO adsorption on MoS(x) clusters: a density functional theory study

The density functional theory (DFT) method has been used to investigate NO probe molecule adsorption on the stoichiometric (Mo(16)S(32)) and nonstoichiometric (Mo(16)S(34) and Mo(16)S(29)) clusters. The calculated adsorption energies indicate that the stoichiometric cluster has stronger NO affinity than the nonstoichiometric surfaces. It is also found that mononitrosyl adsorption is favored at low NO coverage, while dinitrosyl (germinal) and (NO)(2) dimer adsorption at high NO coverage are possible. Strong repulsive interaction has been found for the adsorbed dinitrosyl and (NO)(2) dimer species. In addition, the computed NO stretching frequencies for the mononitrosyl and dinitrosyl species agree well with the experimental data, while those of the dimer species are much lower than the suggested experimental data.     

A density functional theory for Lennard-Jones fluids in cylindrical pores and its applications to adsorption of nitrogen on MCM-41 materials

A density functional theory (DFT) constructed from the modified fundamental-measure theory and the modified Benedict-Webb-Rubin equation of state is presented. The Helmholtz free energy functional due to attractive interaction is expressed as a functional of attractive weighted-density in which the weight function is a mean-field-like type. An obvious advantage of the present theory is that it reproduces accurate bulk properties such as chemical potential, bulk pressure, vapor-liquid interfacial tension, and so forth when compared with molecular simulations and experiments with the same set of molecular parameters. Capabilities of the present DFT are demonstrated by its applicability to adsorption of argon and nitrogen on, respectively, a model cylindrical pore and mesoporous MCM-41 materials. Comparison of the theoretical results of argon in the model cylindrical pore with those from the newly published molecular simulations indicates that the present DFT predicts accurate average densities in the pore, slightly overestimates the pore pressure, and correctly describes the effect of the fluid-pore wall interaction on average densities and pressures in the pore. Application to adsorption of nitrogen on MCM-41 at 77.4 K shows that the present DFT predicts density profiles and adsorption isotherms in good agreement with those from molecular simulations and experiments. In contrast, the hysteresis loop of adsorption calculated from the mean-field theory shifts toward the low pressure region because a low bulk saturated pressure is produced from the mean-field equation of state. The present DFT offers a good way to describe the adsorption isotherms of porous materials as a function of temperature and pressure.     

A density functional theory study of hydrogen adsorption in MOF-5

Ab initio molecular dynamics in the generalized gradient approximation to density functional theory and ground-state relaxations are used to study the interaction between molecular hydrogen and the metal-organic framework with formula unit Zn4O(O2C-C6H4-CO2)3. Five symmetrically unique adsorption sites are identified, and calculations indicate that the sites with the strongest interaction with hydrogen are located near the Zn4O clusters. Twenty total adsorption sites are found around each Zn4O cluster, but after 16 of these are populated, the interaction energy at the remaining four sites falls off significantly. The adsorption of hydrogen on the pore walls creates an attractive potential well for hydrogen in the center of the pore. The effect of the framework on the physical structure and electronic structure of the organic linker is calculated, suggesting ways by which the interaction between the framework and hydrogen could be modified.     

Filling and emptying transitions in cylindrical channels: a density functional approach

The authors use density functional theory in a square gradient approximation to investigate capillary condensation and evaporation in cylindrical channels of finite lengths. The model allows them to systematically explore the effect of the pore's length, width, and surface fields on the location of the transition between "empty" (vapor-filled) and "full" (liquid-filled) states. In general, their results indicate that decreasing the length of the channel drastically reduces the range of pore widths where a transition between liquidlike and vaporlike configurations may be observed. For the wide pores, the transition occurs at very low pressures where the liquid is no longer stable, while for the narrow pores, the transition is hindered by the solid-fluid interactions that favor the vapor phase in lyophobic pores. For the limited range of sizes where the transition can occur, the authors' results confirm the existence of two competing minima that may explain the density oscillations observed in a computer simulation of nanochannels. Comparisons between these results with those generated using a phenomenological model based on the capillary approximation indicate that this simplified approach yields fairly good predictions for medium size pores. However, the capillary approach fails to properly describe the properties of the very small and very large pores.     

Electrical potential in a cylindrical double layer: a functional theory approach

Employing an iterative method in functional theory, the electrical potential distribution for the case of a cylindrical surface is solved. Although the analytical result derived is of an iterative nature, the second-order solution is found to be sufficiently accurate under conditions of practical significance. For the case of constant surface potential, the radius and the surface potential of a cylindrical surface can be estimated based on the extreme of the electrical potential distribution. The effects of the key parameters, including the number and the valence of the ions on a surface, the length of a particle, the relative permittivity of the liquid phase, the temperature, and the concentration of electrolyte on the surface potential, are examined. The general behavior of these effects is similar to that for a spherical surface, except that the surface potential of a cylindrical surface is independent of the electrolyte concentration. The present approach is also applicable to the case where a cylindrical surface remains at a constant charge density.     

Characterization of methoxy adsorption on some transition metals: a first principles density functional theory study

Based on the gradient-density functional theory, calculation results of methoxy adsorption on Au(111), Ag(111), Cu(111), Pt(111), Pd(111), Ni(111), Rh(111), and Fe(100) surfaces are presented, and a consistent picture for some key physical properties determining the reactivity of metals appears. These eight metals belong to two groups: either with filled d electrons (group IB) or with unfilled but more than half filled d electrons (group VIII). The calculated adsorption energies are quite in agreement with the experimental data as well as the previous theoretical calculation results. Importantly, using the analysis of B. Hammer and J. K. Norskov, Nature (London) 376, 232 (1995) and in Chemisorption and Reactivity on Supported Clusters and Thin Films, edited by R. M. Lambert and G. Pacchioni (Kluwer Academic, Dordrecht, 1997), pp. 285-351, the binding energies have selectively been linearly correlated to the d-band center and to the size of the metal d-band orbital overlapping with the adsorbate (coupling matrix element) for these two groups of metals. And by analyzing the nature of the adsorption bonding, the possible reason of this difference is suggested.     

Acetone adsorption on ice investigated by X-ray spectroscopy and density functional theory

Using a combination of X-ray photoemission and near-edge X-ray absorption spectroscopy (NEXAFS) as well as density-functional theory (DFT), we have investigated the adsorption of acetone on ice in the temperature range from 218 to 245 K. The adsorption enthalpy determined from experiment (45 kJ mol(-1)) agrees well with the adsorption energy predicted by theory (41 to 44 kJ mol(-1)). Oxygen K-edge NEXAFS spectra indicate that the presence of acetone at the ice surface does not induce the formation of a pre-melted layer at temperatures up to 243 K. DFT calculations show that the energetically most favored adsorption geometry for acetone on ice is with the molecular plane almost parallel to the surface.     

Hydrogen bonding versus coordination of adsorbate molecules on Ti-silicalites: a density functional theory study

The present study discusses the results of theoretical calculations obtained at the B3LYP/ 6-31G level on the structural, electronic, and energetic properties of Ti-silicalites. Particularly, the relevance of 5T cluster models, either H- or OH-terminated, in large-scale calculations has been critically considered. It was shown that an open surface structure with one OH group and a closed-bulk structure with no bonded OH group at the Ti site are responsible for the observed UV-vis properties of Ti-silicalite materials. Both water and methanol can preferably interact with Ti-silicalites through the H-bonding mechanism, while ammonia can form either H-bonded or coordination complexes. The calculations support the existence of highly dispersed Ti sites in a tetrahedral environment only in Ti-silicalites because an increase in the coordination number of the Ti site by next-neighbor lattice oxygens is the energetically less favorable process.     

Halide adsorption on single-crystal silver substrates: dynamic simulations and ab initio density functional theory

We investigate the static and dynamic behaviors of a Br adlayer electrochemically deposited onto single-crystal Ag(100) using an off-lattice model of the adlayer. Unlike previous studies using a lattice-gas model, the off-lattice model allows adparticles to be located at any position within a two-dimensional approximation to the substrate. Interactions with the substrate are approximated by a corrugation potential. Using density functional theory (DFT) to calculate surface binding energies, a sinusoidal approximation to the corrugation potential is constructed. A variety of techniques, including Monte Carlo and Langevin simulations, are used to study the behavior of the adlayer. The lateral root-mean-square (rms) deviation of the adparticles from the binding sites is presented along with equilibrium coverage isotherms, and the thermally activated Arrhenius barrier-hopping model used in previous dynamic Monte Carlo simulations is tested.     

Equilibrium structure and Ti-catalyzed H2 desorption in NaAlH4 nanoparticles from density functional theory

Improving the hydrogen ab- and desorption kinetics in complex hydrides is essential if these materials are to be used as reversible hydrogen storage media in the transport sector. Although reductions in particle size and the addition of titanium based compounds have been found to improve the kinetics significantly, the physical understanding remains elusive. Density functional theory is used to calculate the energy of the potential low energy surfaces of NaAlH(4) to establish the equilibrium particle shape, and furthermore to determine the deposition energy of Ti/TiH(2) and the substitutional energy for Ti@Al and Ti@Na-sites on the exposed facets. The substitutional processes are energetically preferred and the Na-vacancy formation energy is found to be strongly reduced in the presence of Ti. The barrier for H(2) desorption is found to depend significantly on surface morphology and in particular on the presence of Ti, where the activation energy for H(2) desorption on NaAlH(4){001} surfaces can drop to 0.98 eV--in good agreement with the experimentally observed activation energy for dehydrogenation.     

Lifting the Pt[100] surface reconstruction through oxygen adsorption: a density functional theory analysis

The adsorption of atomic oxygen on unreconstructed Pt[100]-(1 x 1) and reconstructed Pt[100]-(5 x 1) was modeled using density-functional theory in an attempt to understand the relative stability of the unreconstructed phase as a function of oxygen coverage. Our calculations showed that at zero temperature the (5 x 1) is more stable than the unreconstructed (1 x 1) phase at zero oxygen coverage. However, oxygen absorption on the Pt[100]-(5 x 1) phase removed the reconstruction, reversing the phase stability. Using thermochemical analysis, we show desorption of oxygen corresponding to a temperature near 730 K, consistent with experimentally observed desorption peaks for oxygen covered (1 x 1) surfaces. These results have ramifications for understanding the full Pt[100](1 x 1)-->Pt[100]hex-R0.7 degrees surface phase transition.     

Pt adsorption on the PbTiO3(110) polar surface: a density functional theory study

The monolayer (ML) and submonolayer Pt on both terminations of PbTiO₃(110) polar surface have been studied by using density functional theory (DFT) with projector‐augmented wave(PAW) potential and a supercell approach. The most favored ML Pt arrangements on PbTiO and O₂ terminations are the hollow site and the short‐bridge site, respectively. By examining the geometries of different ML arrangements, we know that the dominant impetus for stability of the favored adsorption site for PbTiO termination is the Pt–Ti interaction (mainly from covalent bonding), while that for O₂ termination is the Pt–O interaction (mainly from ionic bonding). In addition, the appearance of the gap electronic states in the outermost layers of each termination indicates that a channel for charge transfer between adsorbed layer and substrate is formed. Moreover, the interface hybridization between Pt 5d and O 2p orbitals is also observed, especially for ML Pt on O₂ termination. The stability sequences for various arrangements of 1/2 ML Pt adsorption conform well with those of ML Pt adsorption, and the most stable arrangement is energetically more favorable than the corresponding ML coverage in the view of adsorption energy maximization. The behavior, i.e. the increase in adsorption energy with decrease in coverage, indicates that Pt? Pt interactions weaken those between Pt and the substrate. Copyright © 2009 John Wiley & Sons, Ltd.     

CO2 adsorption on TiO2(101) anatase: a dispersion-corrected density functional theory study

Adsorption, diffusion, and dissociation of CO(2) on the anatase (101) surface were investigated using dispersion-corrected density functional theory. On the oxidized surface several different local minima were identified of which the most stable corresponds to a CO(2) molecule adsorbed at a five-fold coordinated Ti site in a tilted configuration. Surface diffusion is characterized by relatively small activation barriers. Preferential diffusion takes place along Ti rows and involves a cartwheel type of motion. The presence of a bridging oxygen defect or a surface interstitial Ti atom allows creation of several new strong binding configurations the most stable of which have bent CO(2) structures with simultaneous bonding to two surface Ti atoms. Subsurface oxygen vacancy or interstitial Ti defects are found to enhance the bonding of CO(2) molecules to the surface. CO(2) dissociation from these defect sites is calculated to be exothermic with barriers less than 21 kcal/mol. The use of such defects for catalytic activation of CO(2) on anatase (101) surface would require a mechanism for their regeneration.     

Collision-induced dissociation and density functional theory studies of CO adsorption over zirconium oxide cluster ions: oxidative and nonoxidative adsorption

Zirconium oxide cluster cations and anions are produced by laser ablation and reacted with CO in a fast flow reactor. The CO adsorption products Zr(x)O(y)CO(+) are observed for most of the generated cationic clusters (Zr(x)O(y)(+) = Zr(2)O(5,6)(+), Zr(3)O(7,8)(+), Zr(4)O(9,10)(+)...) while only specific anionic systems (Zr(x)O(y)(-) = Zr(3)O(7)(-), Zr(4)O(9)(-)...) absorb CO to produce Zr(x)O(y)CO(-). To study how the CO molecule is adsorbed on the clusters, the Zr(x)O(y)CO(±) products are mass-selected by a time-of-flight mass spectrometer (TOF-MS) and collided with a crossed helium beam. The fragment ions from collision-induced dissociation (CID) are detected by a secondary TOF-MS. Loss of CO and CO(2) is observed upon the collision of the helium beam with Zr(x)O(y)CO(+) and Zr(x)O(2x+1)CO(-), respectively. Density functional theory calculations indicate that oxidative and nonoxidative adsorption of CO takes place over Zr(3)O(7)(-) and Zr(3)O(7)(+), respectively. This is consistent with the CID experiments.