Assessment of the pairwise additive approximation and evaluation of many-body terms for water clusters

We have tested the ability of four commonly used density functionals (three of which are semilocal and one of which is nonlocal) to outperform accurate pairwise additive approximations in the prediction of binding energies for a series of water clusters ranging in size from dimer to pentamer. Comparison to results obtained with the Weizmann-1 (W1) level of wave function theory shows that while all density functionals are capable of outperforming the accurate pairwise data, the choice of basis set used is crucial to the performance of the method, and if a poor choice of basis set is made the errors obtained with density functional theory (DFT) can be larger than those obtained with the simple pairwise approximation. We have also compared the binding energies and many-body terms determined with DFT to those obtained with W1, and have found that all semilocal functionals have significant errors in the many-body components of the full interactions energy. Despite this limitation, however, we have found that, of the four functionals tested, PBE1W/MG3S is the most accurate for predicting the binding energies of the clusters.     

Pairwise additive model for the He-MgO(100) interaction

We develop a model, based on pairwise additive He-Mg and He-O interactions, for the potential energy of He adsorbates above a rigid MgO(100) surface. The attractive long-range He-Mg and He-O interactions are assumed to have the form C(6)/r(6), with the C(6) coefficients determined from atomic data within the context of the Slater-Kirkwood approximation. The repulsive short-range He-Mg and He-O interactions are assumed to have the form C(p)/r(p), with the exponent p and the C(p) coefficients taken as adjustable parameters. We find that for p = 9, the C(p) coefficients can be chosen so that the laterally averaged He-MgO(100) pairwise additive interaction supports low-lying selective adsorption states, some of whose energies agree very well with the states' apparent energies inferred from experimental measurements. However, for realistic values of the adjustable parameters that define our model, the lateral corrugation of the model pairwise additive He-MgO(100) potential energy surface far exceeds the corrugation that has been inferred both from experimental measurements and from density functional calculations of the short-range He-MgO(100) interaction.     

Evaluation of Desorption Energy Distributions from TPD Spectra on a Heterogeneous Solid Surface

The condensation approximation (CA) and numerical regularization procedure (RP) methods used to solve a Fredholm integral equation of the first kind describing the adsorption equilibria on a heterogeneous solid surface under isothermal conditions have been adopted in the present study to evaluate desorption energy distributions from temperature-programmed desorption (TPD) spectra. From comparisons of the computational results obtained by means of these methods on the basis of simulated TPD spectra, it follows that the CA gives stable solutions for wide desorption energy distributions and it can be used successfully for calculations from wide and clear resolved peaks in the TPD spectra. The use of the RP is more advantageous for acquisition of the distributions from closely related narrow peaks in the TPD spectra. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the Role of Charge Transfer in Many-Body Non-Covalent Interactions**

Charge transfer is one of the mechanisms involved in non-covalent interactions. In molecular dimers, its contribution to pairwise interaction energies has been studied extensively using a variety of interaction energy decomposition schemes. In polar interactions such as hydrogen bonds, it can contribute ten or several tens of percent of the interaction energy. Less is known about its importance in higher-order interactions in many-body systems, mainly because of the lack of methods applicable to this problem. In this work, we extend our method for the quantification of the charge-transfer energy based on constrained DFT to many-body cases and apply it to model trimers extracted from molecular crystals. Our calculations show that charge transfer can account for a large fraction of the total three-body interaction energy. This also has implications for DFT calculations of many-body interactions in general as it is known that many DFT functionals struggle to describe charge-transfer effects correctly.     

Pattern of adsorption isotherms in Ono-Kondo coordinates

The Ono-Kondo lattice density functional theory is used to analyze adsorbate-adsorbate interactions for supercritical systems. In prior work, this approach has been used to study intermolecular interactions in subcritical adsorbed phases, and this has included the study of adsorbate-adsorbate repulsions in the regime of adsorption compression. In this paper, we present the general pattern of adsorption isotherms in Ono-Kondo coordinates; this has not been done in the past. For this purpose, experimental isotherms for adsorption of supercritical fluids (including nitrogen, methane, and carbon dioxide) are plotted in Ono-Kondo coordinates. In addition, we performed Grand Canonical Monte Carlo simulations of adsorption for Lennard-Jones molecules and plotted isotherms in Ono-Kondo coordinates. Our results indicate a pattern of isotherms with regimes of adsorbate-adsorbate attractions at low surface coverage and adsorbate-adsorbate repulsions at high surface coverage. When the generalized Ono-Kondo model is used over a wide range of pressures - from low pressures of the Henry's law regime to supercritical pressures - the slope of the isotherm varies from positive at low pressures to negative at high pressures. The linear sections of these graphs show when the adsorbate-adsorbate interaction energies are approximately constant. When these linear sections have negative slopes, it indicates that the system is in a state of adsorption compression.     

Benchmarking the accuracy of coverage-dependent models: adsorption and desorption of benzene on Pt (1 1 1) and Pt3Sn (1 1 1) from first principles

Bimetallic catalysts have demonstrated properties favorable for upgrading biofuel through catalytic hydrodeoxygenation. However, the design and optimization of such bimetallic catalysts requires the ability to construct accurate, predictive models of these systems. To generate a model that predicts the kinetic behavior of benzene adsorbed on Pt (1 1 1) and a Pt₃Sn (1 1 1) surface alloy (Pt₃Sn (1 1 1)), the adsorption of benzene was studied for a wide range of benzene coverages on both surfaces using density functional theory (DFT) calculations. The adsorption energy of benzene was found to correlate linearly with benzene coverage on Pt (1 1 1) and Pt₃Sn (1 1 1); both surfaces exhibited net repulsive lateral interactions. Through an analysis of the d-band properties of the metal surface, it was determined that the coverage dependence is a consequence of the electronic interactions between benzene and the surface. The linear coverage dependence of the adsorption energy allowed us to quantify the influence of the lateral interactions on the heat of adsorption and temperature programmed desorption (TPD) spectra using a mean-field model. A comparison of our simulated TPD to experiment showed that this mean-field model adequately reproduces the desorption behavior of benzene on Pt (1 1 1) and Pt₃Sn (1 1 1). In particular, the TPD correctly exhibits a broadening desorption peak as the initial coverage of benzene increases on Pt (1 1 1) and a low temperature desorption peak on Pt₃Sn (1 1 1). However, due to the sensitivity of the TPD peak temperature to the desorption energy, precise alignment of experimental and theoretical TPD spectra demands an accurate calculation of the adsorption energy. Therefore, an analysis of the effect of the exchange-correlation functional on TPD modeling is presented. Through this work, we show the necessity of incorporating lateral interactions into theoretical models in order to correctly predict experimental behavior.     

Reflection absorption infrared spectroscopy and temperature-programmed desorption studies of the adsorption and desorption of amorphous and crystalline water on a graphite surface

Reflection absorption infrared spectroscopy (RAIRS) and temperature-programmed desorption (TPD) have been used to perform a detailed investigation of the adsorption of water on highly oriented pyrolytic graphite (HOPG) at 90 K. RAIRS shows that water is physisorbed on HOPG at all coverages, as expected. Experiments at higher surface temperatures show marked changes in the O-H stretching region of the spectrum which can be assigned to the observation of the amorphous to crystalline ice phase transition. The infrared signature of both phases of solid water has been determined on HOPG and can be used to identify the phase of the ice. TPD spectra show the desorption of multilayers of crystalline ice. At high exposures a small bump appears in the TPD spectrum, on the low temperature side of the main peak, which is attributed to the amorphous to crystalline phase transition. At very low exposures of water, it is possible to distinguish the desorption of water from two- and three-dimensional islands and hence to determine the growth mode of water on the HOPG surface. Isothermal TPD studies have also been performed and show that the desorption of water does not obey perfect zero-order kinetics. Desorption orders, derived directly from the TPD spectra, confirm this observation. Desorption energies and preexponential factors have also been determined for this adsorption system.     

Physisorption of N2, O2, and CO on fully oxidized TiO2(110)

Physisorption of N(2), O(2), and CO was studied on fully oxidized TiO(2)(110) using beam reflection and temperature-programmed desorption (TPD) techniques. Sticking coefficients for all three molecules are nearly equal (0.75 +/- 0.05) and approximately independent of coverage suggesting that adsorption occurs via a precursor-mediated mechanism. Excluding multilayer coverages, the TPD spectra for all three adsorbates exhibit three distinct coverage regimes that can be interpreted in accord with previous theoretical studies of N(2) adsorption. At low coverages (0-0.5 N(2)/Ti(4+)), N(2) molecules bind head-on to five-coordinated Ti(4+) ions. The adsorption occurs preferentially on the Ti(4+) sites that do not have neighboring adsorbates. This arrangement minimizes the repulsive interactions between the adsorbed molecules along the Ti(4+) rows resulting in a relatively small shift of the TPD peak (105 --> 90 K) with increasing coverage. At higher N(2) coverages (0.5-1.0 N(2)/Ti(4+)) the nearest-neighbor Ti(4+) sites become occupied. The close proximity of the adsorbates results in strong repulsion thus giving rise to a significant shift of the TPD leading edges (90 --> 45 K) with increasing coverage. For N(2)/Ti(4+) > 1, an additional low-temperature peak (approximately 43 K) is present and is ascribed to N(2) adsorption on bridge-bonded oxygen rows. The results for O(2) and CO are qualitatively similar. The repulsive adsorbate-adsorbate interactions are largest for CO, most likely due to alignment of CO dipole moments. The coverage-dependent binding energies of O(2), N(2), and CO are determined by inverting TPD profiles.     

Water-enhanced low-temperature CO oxidation and isotope effects on atomic oxygen-covered Au(111)

Water-oxygen interactions and CO oxidation by water on the oxygen-precovered Au(111) surface were studied by using molecular beam scattering techniques, temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. Water thermally desorbs from the clean Au(111) surface with a peak temperature of approximately 155 K; however, on a surface with preadsorbed atomic oxygen, a second water desorption peak appears at approximately 175 K. DFT calculations suggest that hydroxyl formation and recombination are responsible for this higher temperature desorption feature. TPD spectra support this interpretation by showing oxygen scrambling between water and adsorbed oxygen adatoms upon heating the surface. In further support of these experimental findings, DFT calculations indicate rapid diffusion of surface hydroxyl groups at temperatures as low as 75 K. Regarding the oxidation of carbon monoxide, if a C (16)O beam impinges on a Au(111) surface covered with both atomic oxygen ( (16)O) and isotopically labeled water (H 2 (18)O), both C (16)O (16)O and C (16)O (18)O are produced, even at surface temperatures as low as 77 K. Similar experiments performed by impinging a C (16)O beam on a Au(111) surface covered with isotopic oxygen ( (18)O) and deuterated water (D 2 (16)O) also produce both C (16)O (16)O and C (16)O (18)O but less than that produced by using (16)O and H 2 (18)O. These results unambiguously show the direct involvement and promoting role of water in CO oxidation on oxygen-covered Au(111) at low temperatures. On the basis of our experimental results and DFT calculations, we propose that water dissociates to form hydroxyls (OH and OD), and these hydroxyls react with CO to produce CO 2. Differences in water-oxygen interactions and oxygen scrambling were observed between (18)O/H 2 (16)O and (18)O/D 2 (16)O, the latter producing less scrambling. Similar differences were also observed in water reactivity toward CO oxidation, in which less CO 2 was produced with (16)O/D 2 (16)O than with (16)O/H 2 (16)O. These differences are likely due to primary kinetic isotope effects due to the differences in O-H and O-D bond energies.     

Implementation of new TPD analysis techniques in the evaluation of second order desorption kinetics of cyanogen from Cu(001)

The interactions of cyanide species with a copper (001) surface were studied with temperature programmed desorption (TPD) and X-ray photoelectron spectroscopy (XPS). Adsorbed cyanide species (CN(a)) undergo recombinative desorption evolving molecular cyanogen (C(2)N(2)). As the adsorbed CN species charge upon adsorption, mutually repulsive dipolar interactions lead to a marked desorption energy reduction with increasing CN(a) coverages. Two new TPD analysis approaches were developed, which used only accurately discernible observables and which do not assume constant desorption energies, E(d), and pre-exponential values, ν. These two approaches demonstrated a linear variation of E(d) with instantaneous coverage. The first approach involved an analysis of the variations of desorption peak asymmetry with initial CN coverages. The second quantitative approach utilized only temperatures and intensities of TPD peaks, together with deduced surface coverages at the peak maxima, also as a function of initial surface coverages. Parameters derived from the latter approach were utilized as initial inputs for a comprehensive curve fit analysis technique. Excellent fits for all experimental desorption curves were produced in simulations. The curve fit analysis confirms that the activation energy of desorption of 170-180 kJ/mol at low coverage decreases by up to 14-15 kJ/mol at CN saturation.     

Six-dimensional potential energy surface for H2 at Ru(0001)

The six-dimensional (6D) potential energy surface (PES) for the H(2) molecule interacting with a clean Ru(0001) surface has been computed accurately for the first time. Density functional theory (DFT) and a pseudopotential based periodic plane-wave approach have been used to calculate the electronic interactions between the molecule and the surface. Two different generalized gradient approximation (GGA) exchange-correlation functionals, PW91 and RPBE, have been adopted. Based on the DFT/GGA calculated potential energies, an analytical 6D PES has been constructed using the corrugation reducing procedure. A very accurate representation of the DFT/GGA data has been achieved, with an average error in the interpolation of about 3 meV and a maximum error not larger than about 30 meV. The top site is found to be the most reactive site for both functionals used, but PW91 predicts a higher reactivity than RPBE, with lower-energy and earlier-located dissociation barriers. The energetic corrugation displayed by the RPBE PES is larger than the PW91 PES while the geometric corrugation is smaller. The differences between the two PESs increase as the distance of the molecular center of mass to the surface decreases. A direct comparison with experimental investigations on H(2)/Ru(0001) could shed light on the suitability of these XC potentials often used in DFT calculations.     

Experimental and theoretical studies of surface nitrate species on Ag/Al2O3 using DRIFTS and DFT

Surface nitrate (NO3(-)) species on the Ag/Al2O3 play an important role in the selective catalytic reduction (SCR) of NOx. In this study, the formation and configuration of surface nitrate NO3(-)(ads) species on Ag/Al2O3 and Al2O3 in the oxidation of NO have been studied using in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and density functional theory (DFT) calculations. Different nitrates species (bridging, bidentate and monodentate) were observed by in situ DRIFTS and validated by DFT calculations results. Attention was especially focused on the proposal of two different bidentate nitrates species (a normal bidentate and an isolated bidentate). In addition, the thermal stability of different surface nitrate species was discussed based on the adsorption energies calculations, DRIFTS, and temperature-programmed desorption (TPD) results. It was suggested that the decomposition and desorption of the surface nitrate species could be controlled by kinetics.     

New Method for the Temperature‐ Programmed Desorption (TPD) of Ammonia Experiment for Characterization of Zeolite Acidity: A Review

In this review, a method for the temperature‐programmed desorption (TPD) of ammonia experiment for the characterization of zeolite acidity and its improvement by simultaneous IR measurement and DFT calculation are described. First, various methods of ammonia TPD are explained, since the measurements have been conducted under the concepts of kinetics, equilibrium, or diffusion control. It is however emphasized that the ubiquitous TPD experiment is governed by the equilibrium between ammonia molecules in the gas phase and on the surface. Therefore, a method to measure quantitatively the strength of the acid site (?H upon ammonia desorption) under equilibrium‐controlled conditions is elucidated. Then, a quantitative relationship between ?H and H₀ function is proposed, based on which the acid strength ?H can be converted into the H₀ function. The identification of the desorption peaks and the quantitative measurement of the number of acid sites are then explained. In order to overcome a serious disadvantage of the method (i.e., no information is provided about the structure of acid sites), the simultaneous measurement of IR spectroscopy with ammonia TPD, named IRMS‐TPD (infrared spectroscopy/mass spectrometry–temperature‐programmed desorption), is proposed. Based on this improved measurement, Brønsted and Lewis acid sites were differentiated and the distribution of Brønsted OH was revealed. The acidity characterized by IRMS‐TPD was further supported by the theoretical DFT calculation. Thus, the advanced study of zeolite acidity at the molecular level was made possible. Advantages and disadvantages of the ammonia TPD experiment are discussed, and understanding of the catalytic cracking activity based on the derived acidic profile is explained.     

Enantiospecific desorption of chiral compounds from chiral Cu(643) and achiral Cu(111) surfaces

Temperature-programmed desorption (TPD) experiments have been conducted to investigate enantiospecific desorption from chiral single-crystal surfaces. The (643) and (six four three) planes of face-centered cubic metals such as Cu have kinked and stepped structures which are nonsuperimposable mirror images of one another and therefore are chiral. These chiral surfaces are denoted Cu(643)(R) and Cu(643)(S). We have observed that the desorption energies of (R)-3-methylcyclohexanone and (R)- and (S)-propylene oxides from the Cu(643)(R) and Cu(643)(S) surfaces depend on the relative handedness of the adsorbate/substrate combination. Since the (643) surface is comprised of terraces with local (111) orientation which are separated by kinked monatomic steps, it is instructive to perform TPD experiments with these chiral compounds on the achiral Cu(111) surface. These experiments have given some insight into the adsorption sites for the chiral molecules on the Cu(643) surfaces. There are several high-temperature features in the TPD spectra of the chiral compounds that only appear in the spectra from the (643) surfaces and thus are attributed to molecules adsorbed at or near the kinked steps. In addition there are lower temperature desorption features observed on the Cu(643) surfaces which occur in the same temperature range as desorption features observed on the Cu(111) surface. These features observed on the (643) surfaces are attributed to desorption from the flat (111) terraces.     

甲醛在CeO2(111)表面吸附的密度泛函理论研究

采用基于第一性原理的密度泛函理论和周期平板模型, 研究了甲醛在以桥氧为端面的CeO2(111)稳定表面上的吸附行为. 通过对不同覆盖度, 不同吸附位的甲醛吸附构型、吸附能及电子态密度的分析发现, 甲醛在CeO2(111)表面存在化学吸附与物理吸附两种情况. 化学吸附结构中甲醛的碳、氧原子分别与表面的氧、铈原子发生相互作用, 形成CH2O2物种; 吸附能随着覆盖度的增加而减小. 与自由甲醛分子相比, 物理吸附的甲醛构型变化不大, 其吸附能较小. 利用CNEB(climbing nudged elastic band)方法计算了甲醛在CeO2(111)表面的初步解离反应活化能(约1.71 eV), 远高于甲醛脱附能垒, 这与甲醛在清洁CeO2(111)表面程序升温脱附实验中产物主要为甲醛的结果相一致.     

Enantioselective separation on a naturally chiral surface

Kinked-stepped, high Miller index surfaces of metal crystals are chiral and, therefore, exhibit enantiospecific properties. Previous temperature-programmed desorption (TPD) spectra have shown that the desorption energies of R-3-methylcyclohexanone (R-3-MCHO) on the chiral Cu(643)(R) and Cu(643)(S) surfaces are enantiospecific (J. Am. Chem. Soc. 2002, 124, 2384). Here, a comparison of the TPD spectra from Cu(111), Cu(221), Cu(533), Cu(653)(R&S), and Cu(643)(R&S) surfaces reveals that the enantiospecific desorption occurs from the chiral kink sites on the Cu(643) surfaces. Titration of the chiral kink sites with I atoms confirms this assignment of desorption features in the TPD spectra. Finally, the enantiospecific difference in the desorption energies of R- and S-3-MCHO has been used as the basis for demonstration of an enantioselective, kinetic separation of racemic 3-MCHO into its purified components during adsorption and desorption on the Cu(643)(R&S) surfaces.     

Adsorption energies, inter-adsorbate interactions, and the two binding sites within monolayer benzene on Ag(111)

The adsorption of monolayer and multilayer benzene on the Ag(111) surface was characterized using temperature programmed desorption (TPD). TPD spectra revealed two broad peaks at approximately 205 and approximately 150 K at submonolayer coverage and a sharper, multilayer peak at 140 K. Analysis of the coverage-dependent shape and shift of the two submonolayer peaks has resulted in their assignment to desorption from two different binding geometries on threefold-hollow sites with symmetries C(3v)(sigma d) and C(3v)(sigma v). The TPD peak analysis incorporated inter-adsorbate repulsive interaction that resulted from the local dipole moment at the adsorption site induced by the adsorbate-surface charge transfer bonding. The analysis has yielded desorption energies of 54.9 +/- 0.8 and 50.4 +/- 0.4 kJ/mol for the C(3v)(sigma d) and C(3v)(sigma v) configurations, respectively. The interface dipole and polarizability of the benzene-silver complex have been determined to be 5.4 +/- 1.8 D and 14 +/- 10 A3, respectively. Repulsive interactions in the monolayer were found to lower the desorption energy from the zero-coverage value by 14.8 kJ/mol. Leading edge analysis of the multilayer peak yielded a desorption energy of 40.9 +/- 0.7 kJ/mol.     

最小二乘法计算苯、噻吩和正辛烷在NaY上程序升温脱附活化能

采用程序升温脱附(TPD)技术测定了苯、噻吩和正辛烷在NaY上以不同升温速率升温时的TPD谱图. 利用TPD谱图的峰形和其微分曲线判断了程序升温脱附过程中的脱附级数. 提出了一种利用最小二乘法计算吸附剂/催化剂的脱附活化能及其动力学参数的方法. 以这些TPD谱图为基础, 分别采用传统TPD计算模型、最小二乘法以及一阶微分曲线法计算了苯、噻吩和正辛烷在NaY上的脱附活化能和动力学参数. 结果表明, 最小二乘法对在不同线性升温速率时的程序升温脱附活化能的计算结果是一致的.     

Adsorption and reaction of 1-epoxy-3-butene on Pt(111): implications for heterogeneous catalysis of unsaturated oxygenates

High-resolution electron energy loss spectroscopy (HREELS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations were used to study the adsorption and reaction of 1-epoxy-3-butene (EpB) on Pt(111). These investigations were conducted to help elucidate mechanisms for improving olefin hydrogenation selectivity in reactions of unsaturated oxygenates. EpB dosed to Pt(111) at 91 K adsorbs molecularly on the surface through the vinyl group with apparent rehybridization to a di-sigma-bound state. By 233 K, however, EpB undergoes epoxide ring opening to form an aldehyde intermediate, which further decomposes upon heating to yield gas phase products CO, H2, and propylene. Comparison of the HREELS and TPD data to experiments performed with 2-butenal (crotonaldehyde) shows that EpB and 2-butenal decompose through related pathways. However, the EpB-derived aldehyde intermediate clearly has a unique structure, features of which have been elucidated by DFT calculations. In conjunction with previous surface science studies of EpB chemistry, these results can help explain selectivity trends for reactions of EpB on Pt catalysts and bimetallic PtAg catalysts, with indications that the enhanced olefin hydrogenation selectivity of PtAg catalysts likely originates from a bifunctional effect.